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.

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.

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.

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.

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.

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.

Multifaceted Study of the Interactions between CPO-27-Ni and Polyurethane and Their Impact on Nitric Oxide Release Performance.

A multifaceted study involving focused ion beam scanning electron microscopy techniques, mechanical analysis, water adsorption measurements, and molecular simulations is employed to rationalize the nitric oxide release performance of polyurethane films containing 5, 10, 20, and 40 wt % of the metal-organic framework (MOF) CPO-27-Ni. The polymer and the MOF are first demonstrated to exhibit excellent compatibility. This is reflected in the even distribution and encapsulation of large wt % MOF loadings throughout the full thickness of the films and by the rather minimal influence of the MOF on the mechanical properties of the polymer at low wt %. The NO release efficiency of the MOF is attenuated by the polymer and found to depend on wt % of MOF loading. The formation of a fully connected network of MOF agglomerates within the films at higher wt % is proposed to contribute to a more complex guest transport in these formulations, resulting in a reduction of NO release efficiency and film ductility. An optimum MOF loading of 10 wt % is identified for maximizing NO release without adversely impacting the polymer properties. Bactericidal efficacy of released NO from the films is demonstrated against Pseudomonas aeruginosa, with a >8 log10 reduction in cell density observed after a contact period of 24 h.

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.

Modeling of Gas Transport through Polymer/MOF Interfaces: A Microsecond-Scale Concentration Gradient-Driven Molecular Dynamics Study.

Membrane-based separation technologies offer a cost-effective alternative to many energy-intensive gas separation processes, such as distillation. Mixed matrix membranes (MMMs) composed of polymers and metal-organic frameworks (MOFs) have attracted a great deal of attention for being promising systems to manufacture durable and highly selective membranes with high gas fluxes and high selectivities. Therefore, understanding gas transport through these MMMs is of significant importance. There has been longstanding speculation that the gas diffusion behavior at the interface formed between the polymer matrix and MOF particles would strongly affect the global performance of the MMMs due to the potential presence of nonselective voids or other defects. To shed more light on this paradigm, we have performed microsecond long concentration gradient-driven molecular dynamics (CGD-MD) simulations that deliver an unprecedented microscopic picture of the transport of H2 and CH4 as single components and as a mixture in all regions of the PIM-1/ZIF-8 membrane, including the polymer/MOF interface. The fluxes of the permeating gases are computed and the impact of the polymer/MOF interface on the H2/CH4 permselectivity of the composite membrane is clearly revealed. Specifically, we show that the poor compatibility between PIM-1 and ZIF-8, which manifests itself by the presence of nonselective void spaces at their interface, results in a decrease of the H2/CH4 permselectivity for the corresponding composite membrane as compared to the performances simulated for PIM-1 and ZIF-8 individually. We demonstrate that CGD-MD simulations based on an accurate atomistic description of the polymer/MOF composite is a powerful tool for characterization and understanding of gas transport and separation mechanisms in MMMs.

Understanding How Ligand Functionalization Influences CO2 and N2 Adsorption in a Sodalite Metal-Organic Framework.

In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal-organic frameworks (MOFs) known as Cu-BTT (BTT-3 = 1,3,5-benzenetristetrazolate) and Cu-BTTri (BTTri-3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu-BTTri compared to its tetrazole counterpart, Cu-BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu-BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu-BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn-Teller effect in Cu-BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials' potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu-BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.

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.

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.

Understanding of the Graphene Oxide/Metal-Organic Framework Interface at the Atomistic Scale.

An atomistic model of the metal-organic framework (MOF) ZIF-8/graphene oxide (GO) interface has been constructed using a combination of density functional theory calculations and force-field-based molecular dynamics simulations. Two microscopic models of GO were constructed integrating basal plane and both basal and edge plane functional groups, called GO-OH and GO-CO2H, respectively. Analysis of the MOF/GO site-to-site interactions, surface coverage, and GO conformation/stiffness and a full characterization of the interfacial region is provided with a special emphasis on the influence of the chemical composition of GO. It was evidenced that the structure of the GO/ZIF-8 composite at the interface is stabilized by a relatively homogeneous set of interactions between the hydrogen atoms of the -NH and -OH terminal functions of ZIF-8 and the oxygen atoms of the epoxy, hydroxyl, and carboxylic groups of GO, leading to an optimal coverage of the MOF surface by GO. Such a scenario implies a significant distortion of the first GO layer brought into contact with the MOF surface, leading to an interfacial region with a relatively small width. This computational exploration strongly suggests that a very good compatibility between these two components would lead, in turn, to the preparation of defect-free ZIF-8/GO films. These predictions are correlated with an experimental effort that consists of successfully prepared homogeneous MOF/GO films that were further characterized by transmission electron microscopy and mechanical testing.

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.

Toward an Understanding of the Microstructure and Interfacial Properties of PIMs/ZIF-8 Mixed Matrix Membranes.

A study integrating advanced experimental and modeling tools was undertaken to characterize the microstructural and interfacial properties of mixed matrix membranes (MMMs) composed of the zeolitic imidazolate framework ZIF-8 nanoparticles (NPs) and two polymers of intrinsic microporosity (PIM-1 and PIM-EA-TB). Analysis probed both the initial ZIF-8/PIM-1 colloidal suspensions and the final hybrid membranes. By combination of dynamic light scattering (DLS) and transmission electron microscopy (TEM) analytical and imaging techniques with small-angle X-ray scattering (SAXS), the colloidal suspensions were shown to consist mainly of two distinct kinds of particles, namely, polymer aggregates of about 200 nm in diameter and densely packed ZIF-8-NP aggregates of a few 100 nm in diameter with a 3 nm thick polymer top-layer. Such aggregates are likely to impart the granular texture of ZIF-8/PIMs MMMs as shown by SEM-XEDS analysis. At the molecular scale, modeling studies showed that the surface coverage of ZIF-8 NPs by both polymers appears not to be optimal with the presence of microvoids at the interfaces that indicates only a moderate compatibility between the polymer and ZIF-8. This study shows that the microstructure of MMMs results from a complex interplay between the ZIF-8/PIM compatibility, solvent, surface chemistry of the ZIF-8 NPs, and the physicochemical properties of the polymers such as molecular structure and rigidity.

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.