Force matching and iterative Boltzmann inversion coarse grained force fields for ZIF-8.

Despite the intense activity at electronic and atomistic resolutions, coarse grained (CG) modeling of metal-organic frameworks remains largely unexplored. One of the main reasons for this is the lack of adequate CG force fields. In this work, we present iterative Boltzmann inversion and force matching (FM) force fields for modeling ZIF-8 at three different coarse grained resolutions. Their ability to reproduce structure, elastic tensor, and thermal expansion is evaluated and compared with that of MARTINI force fields considered in previous work [Alvares et al., J. Chem. Phys. 158, 194107 (2023)]. Moreover, MARTINI and FM are evaluated for their ability to depict the swing effect, a subtle phase transition ZIF-8 undergoes when loaded with guest molecules. Overall, we found that all our force fields reproduce structure reasonably well. Elastic constants and volume expansion results are analyzed, and the technical and conceptual challenges of reproducing them are explained. Force matching exhibits promising results for capturing the swing effect. This is the first time these CG methods, widely applied in polymer and biomolecule communities, are deployed to model porous solids. We highlight the challenges of fitting CG force fields for these materials.

Pore-Networked Gels: Permanently Porous Ionic Liquid Gels with Linked Metal-Organic Polyhedra Networks.

Porous liquids (PLs) are attractive materials because of their capability to combine the intrinsic porosity of microporous solids and the processability of liquids. Most of the studies focus on the synthesis of PLs with not only high porosity but also low viscosity by considering their transportation in industrial plants. However, a gap exists between PLs and solid adsorbents for some practical cases, where the liquid characteristics and mechanical stability without leakage are simultaneously required. Here, we fill in this gap by demonstrating a new concept of pore-networked gels, in which the solvent phase is trapped by molecular networks with accessible porosity. To achieve this, we fabricate a linked metal-organic polyhedra (MOPs) gel, followed by exchanging the solvent phase with a bulky liquid such as ionic liquids (ILs); the dimethylformamide solvent trapped inside the as-synthesized gel is replaced by the target IL, 1-butyl-3-methylimidazolium tetrafluoroborate, which in turn cannot enter MOP pores due to their larger molecular size. The remaining volatile solvents in the MOP cavities can then be removed by thermal activation, endowing the obtained IL gel (Gel_IL) with accessible microporosity. The CO2 capacities of the gels are greatly enhanced compared to the neat IL. The exchange with the IL also exerts a positive influence on the final gel performances such as mechanical properties and low volatility. Besides ILs, various functional liquids are shown to be amenable to this strategy to fabricate pore-networked gels with accessible porosity, demonstrating their potential use in the field of gas adsorption or separation.

Coarse-grained modeling of zeolitic imidazolate framework-8 using MARTINI force fields.

In this contribution, the well-known MARTINI particle-based coarse graining approach is tested for its ability to model the ZIF-8 metal-organic framework. Its capability to describe structure, lattice parameters, thermal expansion, elastic constants and amorphization is evaluated. Additionally, the less coarsened models were evaluated for reproducing the swing effect and the host-guest interaction energies were analyzed. We find that MARTINI force fields successfully capture the structure of the Metal-Organic Framework (MOF) for different degrees of coarsening, with the exception of the MARTINI 2.0 models for the less coarse mapping. MARTINI 2.0 models predict more accurate values of C11 and C12, while MARTINI 3.0 has a tendency to underestimate them. Among the possibilities tested, the choice of bead flavors within a particular MARTINI version appears to have a less critical impact in the simulated properties of the empty framework. None of the coarse-grained (CG) models investigated were able to capture the amorphization nor the swing effect within the scope of MD simulations. A perspective on the importance of having a proper Lennard-Jones (LJ) parametrization for modeling guest-MOF and MOF-MOF interactions is highlighted.

Gas-sieving zeolitic membranes fabricated by condensation of precursor nanosheets.

The synthesis of molecular-sieving zeolitic membranes by the assembly of building blocks, avoiding the hydrothermal treatment, is highly desired to improve reproducibility and scalability. Here we report exfoliation of the sodalite precursor RUB-15 into crystalline 0.8-nm-thick nanosheets, that host hydrogen-sieving six-membered rings (6-MRs) of SiO4 tetrahedra. Thin films, fabricated by the filtration of a suspension of exfoliated nanosheets, possess two transport pathways: 6-MR apertures and intersheet gaps. The latter were found to dominate the gas transport and yielded a molecular cutoff of 3.6 Å with a H2/N2 selectivity above 20. The gaps were successfully removed by the condensation of the terminal silanol groups of RUB-15 to yield H2/CO2 selectivities up to 100. The high selectivity was exclusively from the transport across 6-MR, which was confirmed by a good agreement between the experimentally determined apparent activation energy of H2 and that computed by ab initio calculations. The scalable fabrication and the attractive sieving performance at 250-300 °C make these membranes promising for precombustion carbon capture.

Structure of the Polymer Backbones in polyMOF Materials.

The molecular connectivity of polymer-metal-organic framework (polyMOF) hybrid materials was investigated using density functional theory calculations and solid-state NMR spectroscopy. The architectural constraints that dictate the formation of polyMOFs were assessed by examining poly(1,4-benzenedicarboxylic acid) (pbdc) polymers in two archetypical MOF lattices (UiO-66 and IRMOF-1). Modeling of the polyMOFs showed that in the IRMOF-1-type lattice, six, seven, and eight methylene (-CH2-) groups between 1,4-benzenedicarboxylate (terephthalate, bdc2-) units can be accommodated without significant distortions, while in the UiO-66-type lattice, an optimal spacing of seven methylene groups between bdc2- units is needed to minimize strain. Solid-state NMR supports these predictions and reveals pronounced spectral differences for the same polymer in the two polyMOF lattices. With seven methylene groups, polyUiO-66-7a shows 7 ± 3% of uncoordinated terephthalate linkers, while these are undetectable (<4%) in the corresponding polyIRMOF-1-7a. In addition, NMR-detected backbone mobility is significantly higher in the polyIRMOF-1-7a than in the corresponding polyUiO-66-7a, again indicative of taut chains in the latter.

Polymer Infiltration into Metal-Organic Frameworks in Mixed-Matrix Membranes Detected in Situ by NMR.

Solid-state NMR has been used to study mixed-matrix membranes (MMMs) prepared with a metal-organic framework (MOF, UiO-66) and two different high molecular weight polymers (PEO and PVDF). 13C and 1H NMR data provide overwhelming evidence that most UiO-66 organic linkers are within 1 nm of PEO, which indicates that PEO is homogeneously distributed throughout the MOF. Systematic changes in MOF 13C NMR peak positions and 1H NMR line widths, as well as dramatic reductions in the MOF 1H T1ρ relaxation times, are observed as the PEO content increases, and when the pores have been filled, a further increase in PEO results in the formation of semicrystalline PEO outside the UiO-66 particles. In contrast, similar studies on PVDF MMMs show that the polymer contacts only a small fraction (<20%) of the MOF linkers. Simulations confirm that PEO penetrates into UiO-66 more easily than does PVDF. These studies are among the first to provide experimental insights into MOF-polymer interactions in an MMM.

A new kind of atlas of zeolite building blocks.

We have analyzed structural motifs in the Deem database of hypothetical zeolites to investigate whether the structural diversity found in this database can be well-represented by classical descriptors, such as distances, angles, and ring sizes, or whether a more general representation of the atomic structure, furnished by the smooth overlap of atomic position (SOAP) method, is required to capture accurately structure-property relations. We assessed the quality of each descriptor by machine-learning the molar energy and volume for each hypothetical framework in the dataset. We have found that a SOAP representation with a cutoff length of 6 Å, which goes beyond near-neighbor tetrahedra, best describes the structural diversity in the Deem database by capturing relevant interatomic correlations. Kernel principal component analysis shows that SOAP maintains its superior performance even when reducing its dimensionality to those of the classical descriptors and that the first three kernel principal components capture the main variability in the dataset, allowing a 3D point cloud visualization of local environments in the Deem database. This "cloud atlas" of local environments was found to show good correlations with the contribution of a given motif to the density and stability of its parent framework. Local volume and energy maps constructed from the SOAP/machine learning analyses provide new images of zeolites that reveal smooth variations of local volumes and energies across a given framework and correlations between the contributions to volume and energy associated with each atom-centered environment.

An In-Situ Neutron Diffraction and DFT Study of Hydrogen Adsorption in a Sodalite-Type Metal-Organic Framework, Cu-BTTri.

Herein we present a detailed study of the hydrogen adsorption properties of Cu-BTTri, a robust crystalline metal-organic framework containing open metal-coordination sites. Diffraction techniques, carried out on the activated framework, reveal a structure that is different from what was previously reported. Further, combining standard hydrogen adsorption measurements with in-situ neutron diffraction techniques provides molecular level insight into the hydrogen adsorption process. The diffraction experiments unveil the location of four D2 adsorption sites in Cu-BTTri and shed light on the structural features that promote hydrogen adsorption in this material. Density functional theory (DFT), used to predict the location and strength of binding sites, corroborate the experimental findings. By decomposing binding energies in different sites in various energetic contributions, we show that van der Waals interactions play a crucial role, suggesting a possible route to enhancing the binding energy around open metal coordination sites.

Revealing the Transient Concentration of CO2 in a Mixed-Matrix Membrane by IR Microimaging and Molecular Modeling.

Through IR microimaging the spatially and temporally resolved development of the CO2 concentration in a ZIF-8@6FDA-DAM mixed matrix membrane (MMM) was visualized during transient adsorption. By recording the evolution of the CO2 concentration, it is observed that the CO2 molecules propagate from the ZIF-8 filler, which acts as a transport "highway", towards the surrounding polymer. A high-CO2 -concentration layer is formed at the MOF/polymer interface, which becomes more pronounced at higher CO2 gas pressures. A microscopic explanation of the origins of this phenomenon is suggested by means of molecular modeling. By applying a computational methodology combining quantum and force-field based calculations, the formation of microvoids at the MOF/polymer interface is predicted. Grand canonical Monte Carlo simulations further demonstrate that CO2 tends to preferentially reside in these microvoids, which is expected to facilitate CO2 accumulation at the interface.

Proton Transport in a Highly Conductive Porous Zirconium-Based Metal-Organic Framework: Molecular Insight.

The water stable UiO-66(Zr)-(CO2H)2 MOF exhibits a superprotonic conductivity of 2.3×10(-3)  S cm(-1) at 90 °C and 95 % relative humidity. Quasi-elastic neutron scattering measurements combined with aMS-EVB3 molecular dynamics simulations were able to probe individually the dynamics of both confined protons and water molecules and to further reveal that the proton transport is assisted by the formation of a hydrogen-bonded water network that spans from the tetrahedral to the octahedral cages of this MOF. This is the first joint experimental/modeling study that unambiguously elucidates the proton-conduction mechanism at the molecular level in a highly conductive MOF.

Lithium solvation in dimethyl sulfoxide-acetonitrile mixtures.

We present molecular dynamics simulation results pertaining to the solvation of Li(+) in dimethyl sulfoxide-acetonitrile binary mixtures. The results are potentially relevant in the design of Li-air batteries that rely on aprotic mixtures as solvent media. To analyze effects derived from differences in ionic size and charge sign, the solvation of Li(+) is compared to the ones observed for infinitely diluted K(+) and Cl(-) species, in similar solutions. At all compositions, the cations are preferentially solvated by dimethyl sulfoxide. Contrasting, the first solvation shell of Cl(-) shows a gradual modification in its composition, which varies linearly with the global concentrations of the two solvents in the mixtures. Moreover, the energetics of the solvation, described in terms of the corresponding solute-solvent coupling, presents a clear non-ideal concentration dependence. Similar nonlinear trends were found for the stabilization of different ionic species in solution, compared to the ones exhibited by their electrically neutral counterparts. These tendencies account for the characteristics of the free energy associated to the stabilization of Li(+)Cl(-), contact-ion-pairs in these solutions. Ionic transport is also analyzed. Dynamical results show concentration trends similar to those recently obtained from direct experimental measurements.

Excess protons in water-acetone mixtures. II. A conductivity study.

In the present work we complement a previous simulation study [R. Semino and D. Laria, J. Chem. Phys. 136, 194503 (2012)] on the disruption of the proton transfer mechanism in water by the addition of an aprotic solvent, such as acetone. We provide experimental measurements of the mobility of protons in aqueous-acetone mixtures in a wide composition range, for water molar fractions, xw, between 0.05 and 1.00. Furthermore, new molecular dynamics simulation results are presented for rich acetone mixtures, which provide further insight into the proton transport mechanism in water-non-protic solvent mixtures. The proton mobility was analyzed between xw 0.05 and 1.00 and compared to molecular dynamics simulation data. Results show two qualitative changes in the proton transport composition dependence at xw ∼ 0.25 and 0.8. At xw < 0.25 the ratio of the infinite dilution molar conductivities of HCl and LiCl, Λ(0)(HCl).Λ(0)(LiCl)(-1), is approximately constant and equal to one, since the proton diffusion is vehicular and equal to that of Li(+). At xw ∼ 0.25, proton mobility starts to differ from that of Li(+) indicating that above this concentration the Grotthuss transport mechanism starts to be possible. Molecular dynamics simulation results showed that at this threshold concentration the probability of interconversion between two Eigen structures starts to be non-negligible. At xw ∼ 0.8, the infinite molar conductivity of HCl concentration dependence qualitatively changes. This result is in excellent agreement with the analysis presented in the previous simulation work and it has been ascribed to the interchange of water and acetone molecules in the second solvation shell of the hydronium ion.

Nucleation of zeolitic imidazolate frameworks: from molecules to nanoparticles.

We have studied the clusters involved in the initial stages of nucleation of Zeolitic Imidazolate Frameworks, employing a wide range of computational techniques. In the pre-nucleating solution, the prevalent cluster is the ZnIm4 cluster (formed by a zinc cation, Zn2+, and four imidazolate anions, Im-), although clusters such as ZnIm3, Zn2Im7, Zn2Im7, Zn3Im9, Zn3Im10, or Zn4Im12 have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm3, have a tetrahedrally coordinated Zn2+ cation. Small ZnxImy clusters are less stable than the ZnIm4 cluster. The first cluster that is found to be more stable than ZnIm4 is the Zn41Im88 cluster, which is a disordered cluster with glassy structure. Bulk-like clusters do not begin to be more stable than glassy clusters until much larger sizes, since the larger cluster we have studied (Zn144Im288) is still less stable than the glassy Zn41Im88 cluster, suggesting that Ostwald's rule (the less stable polymorph crystallizes first) could be fulfilled, not for kinetic, but for thermodynamic reasons. Our results suggest that the first clusters formed in the nucleation process would be glassy clusters, which then undergo transformation to any of the various crystal structures possible, depending on the kinetic routes provided by the synthesis conditions. Our study helps elucidate the way in which the various species present in solution interact, leading to nucleation and crystal growth.

Excess protons in water-acetone mixtures.

Using molecular dynamics experiments, we analyze equilibrium and dynamical characteristics related to the solvation of excess protons in water-acetone mixtures. Our approach is based on the implementation of an extended valence-bond Hamiltonian, which incorporates translocation of the excess charge between neighboring water molecules. Different mixtures have been analyzed, starting from the pure water case down to solutions with a water molar fraction x(w) = 0.25. In all cases, we have verified that the structure of the first solvation shell of the H(3)O(+) moiety remains practically unchanged, compared to the one observed in pure water. This shell is composed by three water molecules acting as hydrogen bond acceptors, with no evidence of hydrogen bond donor-like connectivity. Moreover, the increment in the acetone concentration leads to a gradual stabilization of Eigen-like [H(3)O[middle dot](H(2)O)(3)](+) configurations, in detriment of Zundel-like [H[middle dot](H(2)O)(2)](+) ones. Rates of proton transfer and proton diffusion coefficients have been recorded at various water-acetone relative concentrations. In both cases, we have found a transition region, in the vicinity of x(w) ∼ 0.8, where the concentration dependences of the two magnitudes change at a quantitative level. A crude estimate shows that, at this tagged concentration, the volumes "occupied" by the two solvents become comparable. The origins of this transition separating water-rich from acetone-rich realms is rationalized in terms of modifications operated in the nearby, second solvation shell, which in the latter solutions, normally includes at least, one acetone molecule. Our results would suggest that one possible mechanism controlling the proton transfer in acetone-rich solutions is the exchange of one of these tagged acetone molecules, by nearby water ones. This exchange would give rise to Zundel-like structures, exhibiting a symmetric, first solvation shell composed exclusively by water molecules, and would facilitate the transfer between neighboring water molecules along the resonant complex.

Excess protons in mesoscopic water-acetone nanoclusters.

We carried out molecular dynamics simulation experiments to examine equilibrium and dynamical characteristics of the solvation of excess protons in mesoscopic, [m:n] binary polar clusters comprising m = 50 water molecules and n = 6, 25, and 100 acetone molecules. Contrasting from what is found in conventional macroscopic phases, the characteristics of the proton solvation are dictated, to a large extent, by the nature of the concentration fluctuations prevailing within the clusters. At low acetone contents, the overall cluster morphology corresponds to a segregated aqueous nucleus coated by an external aprotic phase. Under these circumstances, the proton remains localized at the surface of the water core, in a region locally deprived from acetone molecules. At higher acetone concentrations, we found clear evidence of the onset of the mixing process. The cluster structures present aqueous domains with irregular shape, fully embedded within the acetone phase. Still, the proton remains coordinated to the aqueous phase, with its closest solvation shell composed exclusively by three water molecules. As the relative concentration of acetone increases, the time scales characterizing proton transfer events between neighboring water molecules show considerable retardations, stretching into the nanosecond time domain already for n ~ 25. In water-rich aggregates, and similarly to what is found in the bulk, proton transfers are controlled by acetone/water exchange processes taking place at the second solvation shell of the proton. As a distinctive feature of the transfer mechanism, translocation pathways also include diffusive motions of the proton from the surface down into inner regions of the underlying water domain.

Putting Forward NUS-8-CO2H/PIM-1 as a Mixed Matrix Membrane for CO2 Capture.

Mixed matrix membranes (MMMs) composed of NUS-8 metal-organic framework (MOF) nanosheets dispersed into a polymer of intrinsic microporosity 1 (PIM-1) polymer matrix are known to be promising candidates for CO2/N2 separation because of a solubility-driven separation mechanism. In this work, we predict that a chemical functionalization of the organic linker of NUS-8 by a CO2-philic function confers an even better separation performance to the resulting MMM. Our simulations revealed that the NUS-8-CO2H/PIM-1 composite exhibits a 3-fold increase in CO2/N2 selectivity versus the NUS-8/PIM-1 analogue while achieving a high CO2 permeability (6700 barrer). We demonstrated that this improved level of performance is due to an increase both in the total MOF/polymer interfacial pore volume and in the CO2-affinity due to the chemical functionalization. These results suggest that an appropriate choice of chemical functionalization of a MOF is a promising strategy to improve gas separation performances for MMM composites that exhibit a solubility-driven separation mechanism.

Ranking the synthesizability of hypothetical zeolites with the sorting hat.

Zeolites are nanoporous alumino-silicate frameworks widely used as catalysts and adsorbents. Even though millions of siliceous networks can be generated by computer-aided searches, no new hypothetical framework has yet been synthesized. The needle-in-a-haystack problem of finding promising candidates among large databases of predicted structures has intrigued materials scientists for decades; yet, most work to date on the zeolite problem has been limited to intuitive structural descriptors. Here, we tackle this problem through a rigorous data science scheme-the "Zeolite Sorting Hat"-that exploits interatomic correlations to discriminate between real and hypothetical zeolites and to partition real zeolites into compositional classes that guide synthetic strategies for a given hypothetical framework. We find that, regardless of the structural descriptor used by the Zeolite Sorting Hat, there remain hypothetical frameworks that are incorrectly classified as real ones, suggesting that they might be good candidates for synthesis. We seek to minimize the number of such misclassified frameworks by using as complete a structural descriptor as possible, thus focusing on truly viable synthetic targets, while discovering structural features that distinguish real and hypothetical frameworks as an output of the Zeolite Sorting Hat. Further ranking of the candidates can be achieved based on thermodynamic stability and/or their suitability for the desired applications. Based on this workflow, we propose three hypothetical frameworks differing in their molar volume range as the top targets for synthesis, each with a composition suggested by the Zeolite Sorting Hat. Finally, we analyze the behavior of the Zeolite Sorting Hat with a hierarchy of structural descriptors including intuitive descriptors reported in previous studies, finding that intuitive descriptors produce significantly more misclassified hypothetical frameworks, and that more rigorous interatomic correlations point to second-neighbor Si-O distances around 3.2-3.4 Å as the key discriminatory factor.

Rational design of mixed-matrix metal-organic framework membranes for molecular separations.

Conventional separation technologies to separate valuable commodities are energy intensive, consuming 15% of the worldwide energy. Mixed-matrix membranes, combining processable polymers and selective adsorbents, offer the potential to deploy adsorbent distinct separation properties into processable matrix. We report the rational design and construction of a highly efficient, mixed-matrix metal-organic framework membrane based on three interlocked criteria: (i) a fluorinated metal-organic framework, AlFFIVE-1-Ni, as a molecular sieve adsorbent that selectively enhances hydrogen sulfide and carbon dioxide diffusion while excluding methane; (ii) tailoring crystal morphology into nanosheets with maximally exposed (001) facets; and (iii) in-plane alignment of (001) nanosheets in polymer matrix and attainment of [001]-oriented membrane. The membrane demonstrated exceptionally high hydrogen sulfide and carbon dioxide separation from natural gas under practical working conditions. This approach offers great potential to translate other key adsorbents into processable matrix.

Insights into the Enhancement of MOF/Polymer Adhesion in Mixed-Matrix Membranes via Polymer Functionalization.

MOF-based mixed-matrix membranes (MMMs) prepared using standard routes often exhibit poor adhesion between polymers and MOFs. Herein, we report an unprecedented systematic exploration on polymer functionalization as the key to achieving defect-free MMMs. As a case study, we explored computationally MMMs based on the combination of the prototypical UiO-66(Zr) MOF with polymer of intrinsic porosity-1 (PIM-1) functionalized with various groups including amidoxime, tetrazole, and N-((2-ethanolamino)ethyl)carboxamide. Distinctly, the amidoxime-derivative PIM-1/UiO-66(Zr) MMM was predicted to express the desired enhanced MOF/polymer interfacial interactions and thus subsequently prepared and evaluated experimentally. Prominently, high-resolution transmission electron microscopy confirmed optimal adhesion between the two components in contrast to the nanometer-sized voids/defects shown by the pristine PIM-1/UiO-66(Zr) MMM. Notably, single-gas permeation measurements further corroborated the need of optimal MOF/polymer adhesion in order to effectively enable the MOF to play a role in the gas transport of the resulting MMM.