Search

The microscopic interfacial structures for a series of metal-organic framework/polymer composites consisting of the Zr-based UiO-66 coupled with different polymers are systematically explored by applying a computational methodology that integrates density functional theory calculations and force field-based molecular dynamics simulations. These predictions are correlated with experimental findings to unravel the structure-compatibility relationship of the MOF/polymer pairs. The relative contributions of the intermolecular MOF/polymer interactions and the flexibility/rigidity of the polymer with respect to the microscopic structure of the interface are rationalized, and their impact on the compatibility of the two components in the resulting composite is discussed. The most compatible pairs among those investigated involve more flexible polymers, i.e. polyvinylidene fluoride (PVDF) and polyethylene glycol (PEG). These polymers exhibit an enhanced contact surface, due to a better adaptation of their configuration to the MOF surface. In these cases, the irregularities at the MOF surface are filled by the polymer, and even some penetration of the terminal groups of the polymer into the pores of the MOF can be observed. As a result, the affinity between the MOF and the polymer is very high; however, the pores of the MOF may be sterically blocked due to the strong MOF/polymer interactions, as evidenced by UiO-66/PEG composites. In contrast, composites involving polymers that exhibit higher rigidity, such as the polymer of intrinsic microporosity-1 (PIM-1) or polystyrene (PS), present interfacial microvoids that contribute to a decrease in the contact surface between the two components, thus reducing the MOF/polymer affinity.

We employ all-atom well-tempered metadynamics simulations to study the mechanistic details of both the early stages of nucleation and crystal decomposition for the benchmark metal–organic framework (MOF) ZIF-8. To do so, we developed and validated a force field that reliably models the modes of coordination bonds via a Morse potential functional form and employs cationic and anionic dummy atoms to capture coordination symmetry. We also explored a set of physically relevant collective variables and carefully selected an appropriate subset for our problem at hand. After a rapid increase of the Zn–N connectivity, we observe the evaporation of small clusters in favor of a few large clusters, which leads to the formation of an amorphous highly connected aggregate. Zn (MIm) 4 2 − and Zn (MIm) 3 − complexes are observed with lifetimes in the order of a few picoseconds, while larger structures, such as four-, five-, and six-membered rings, have substantially longer lifetimes of a few nanoseconds. The free ligands act as "templating agents" for the formation of sodalite cages. ZIF-8 crystal decomposition results in the formation of a vitreous phase. Our findings contribute to a fundamental understanding of MOF's synthesis that paves the way to controlling synthesis products. Furthermore, our developed force field and methodology can be applied to model solution processes that require coordination bond reactivity for other ZIFs besides ZIF-8. [ABSTRACT FROM AUTHOR]

We have performed a data science study of Monte Carlo (MC) simulation trajectories to understand factors that can accelerate the formation of zeolite nanoporous crystals, a process that can take days or even weeks. In previous work, MC simulations predicted and experiments confirmed that using a secondary organic structure-directing agent (OSDA) accelerates the crystallization of all-silica LTA zeolite, with experiments finding a three-fold speedup [Bores et al., Phys. Chem. Chem. Phys. 24, 142-148 (2022)]. However, it remains unclear what physical factors cause the speed-up. Here, we apply data science to analyze the simulation trajectories to discover what drives accelerated zeolite crystallization in MC simulations going from a one-OSDA synthesis (1OSDA) to a two-OSDA version (2OSDA). We encoded simulation snapshots using the smooth overlap of atomic positions approach, which represents all two- and three-body correlations within a given cutoff distance. Principal component analyses failed to discriminate datasets of structures from 1OSDA and 2OSDA simulations, while the Support Vector Machine (SVM) approach succeeded at classifying such structures with an area-under-curve (AUC) score of 0.99 (where AUC = 1 is a perfect classification) with all three-body correlations and as high as 0.94 with only two-body correlations. SVM decision functions reveal relatively broad/narrow histograms for 1OSDA/2OSDA datasets, suggesting that the two simulations differ strongly in information heterogeneity. Informed by these results, we performed pair (2-body) entropy calculations during crystallization, resulting in entropy differences that semi-quantitatively account for the speedup observed in the previous MC simulations. We conclude that altering synthesis conditions in ways that substantially change the entropy of labile silica networks may accelerate zeolite crystallization, and we discuss possible approaches for achieving such acceleration., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

The current paradigm considers the control of the MOF/polymer interface mostly for achieving a good compatibility between the two components to ensure the fabrication of continuous mixed-matrix metal-organic framework (MMMOF) membranes. Here, we unravel that the interfacial pore shape nanostructure plays a key role for an optimum molecular transport. The prototypical ultrasmall pore AlFFIVE-1-Ni MOF was assembled with the polymer PIM-1 to design a composite with gradually expanding pore from the MOF entrance to the MOF/polymer interfacial region. Concentration gradient-driven molecular dynamics simulations demonstrated that this pore nanostructuring enables an optimum guided path for the gas molecules at the MOF/polymer interface that decisively leads to an acceleration of the molecular transport all along the MMMOF membrane. This numerical prediction resulted in the successful fabrication of a [001]-oriented nanosheets AlFFIVE-1-Ni/PIM-1 MMMOF membrane exhibiting an excellent CO 2 permeability, better than many MMMs, and ideally associated with a sufficiently high CO 2 /CH 4 selectivity that makes this membrane very promising for natural gas/biogas purification.

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., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

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 CO 2 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.

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., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

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 ZnIm 4 cluster (formed by a zinc cation, Zn 2+ , and four imidazolate anions, Im - ), although clusters such as ZnIm 3 , Zn 2 Im 7 , Zn 2 Im 7 , Zn 3 Im 9 , Zn 3 Im 10 , or Zn 4 Im 12 have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm 3 , have a tetrahedrally coordinated Zn 2+ cation. Small Zn x Im y clusters are less stable than the ZnIm 4 cluster. The first cluster that is found to be more stable than ZnIm 4 is the Zn 41 Im 88 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 (Zn 144 Im 288 ) is still less stable than the glassy Zn 41 Im 88 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.

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., Competing Interests: The authors declare no competing interests., (This journal is © The Royal Society of Chemistry.)

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.

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 CO 2 /N 2 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 CO 2 -philic function confers an even better separation performance to the resulting MMM. Our simulations revealed that the NUS-8-CO 2 H/PIM-1 composite exhibits a 3-fold increase in CO 2 /N 2 selectivity versus the NUS-8/PIM-1 analogue while achieving a high CO 2 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 CO 2 -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.

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.

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 SiO 4 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 H 2 /N 2 selectivity above 20. The gaps were successfully removed by the condensation of the terminal silanol groups of RUB-15 to yield H 2 /CO 2 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 H 2 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.

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 log 10 reduction in cell density observed after a contact period of 24 h.

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 (-CH 2 -) groups between 1,4-benzenedicarboxylate (terephthalate, bdc 2- ) units can be accommodated without significant distortions, while in the UiO-66-type lattice, an optimal spacing of seven methylene groups between bdc 2- 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.

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 H 2 and CH 4 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 H 2 /CH 4 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 H 2 /CH 4 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., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

In this work, a detailed study is conducted to understand how ligand substitution influences the CO 2 and N 2 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 CO 2 and N 2 binding sites on the internal surface of Cu-BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu 2+ sites and surface-bound CO 2 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 CO 2 and N 2 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 CO 2 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 CO 2 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., Competing Interests: The authors declare no competing financial interest.

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.

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). 13 C and 1 H 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 13 C NMR peak positions and 1 H NMR line widths, as well as dramatic reductions in the MOF 1 H T 1ρ 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.

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 D 2 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.

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-CO 2 H, 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.

Through IR microimaging the spatially and temporally resolved development of the CO 2 concentration in a ZIF-8@6FDA-DAM mixed matrix membrane (MMM) was visualized during transient adsorption. By recording the evolution of the CO 2 concentration, it is observed that the CO 2 molecules propagate from the ZIF-8 filler, which acts as a transport "highway", towards the surrounding polymer. A high-CO 2 -concentration layer is formed at the MOF/polymer interface, which becomes more pronounced at higher CO 2 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 CO 2 tends to preferentially reside in these microvoids, which is expected to facilitate CO 2 accumulation at the interface., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The microscopic interfacial structures for a series of metal-organic framework/polymer composites consisting of the Zr-based UiO-66 coupled with different polymers are systematically explored by applying a computational methodology that integrates density functional theory calculations and force field-based molecular dynamics simulations. These predictions are correlated with experimental findings to unravel the structure-compatibility relationship of the MOF/polymer pairs. The relative contributions of the intermolecular MOF/polymer interactions and the flexibility/rigidity of the polymer with respect to the microscopic structure of the interface are rationalized, and their impact on the compatibility of the two components in the resulting composite is discussed. The most compatible pairs among those investigated involve more flexible polymers, i.e. polyvinylidene fluoride (PVDF) and polyethylene glycol (PEG). These polymers exhibit an enhanced contact surface, due to a better adaptation of their configuration to the MOF surface. In these cases, the irregularities at the MOF surface are filled by the polymer, and even some penetration of the terminal groups of the polymer into the pores of the MOF can be observed. As a result, the affinity between the MOF and the polymer is very high; however, the pores of the MOF may be sterically blocked due to the strong MOF/polymer interactions, as evidenced by UiO-66/PEG composites. In contrast, composites involving polymers that exhibit higher rigidity, such as the polymer of intrinsic microporosity-1 (PIM-1) or polystyrene (PS), present interfacial microvoids that contribute to a decrease in the contact surface between the two components, thus reducing the MOF/polymer affinity.

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.

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., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

An innovative computational methodology integrating density functional theory calculations and force field-based molecular dynamics simulations was developed to provide a first microscopic model of the interactions at the metal-organic framework (MOF) surface/polymer interface. This was applied to the case of the composite formed by the polymer of intrinsic microporosity, PIM-1, and the zeolitic imidazolate framework, ZIF-8, as a model system. We found that the structure of the composite at the interface is the result of both the chemical affinity between PIM-1 and ZIF-8 and the rigidity of the polymer. Specifically, there is a preferential interaction between the -CN groups of PIM-1 and the NH terminal functions of the organic linker at the ZIF-8 surface. Additionally, the resulting conformation of the polymer gives rise to interfacial microvoids at the vicinity of the MOF surface. The porosity, rigidity, and density of the interfacial polymer were analyzed and compared to those for the bulk polymer. It was shown that the polymer still feels the impact of the MOF surface even at long distances above 15-20 Å. Further, both the polydispersity of the polymer and the flexibility of the MOF surface were revealed to only slightly affect the properties of the MOF/interface. This work, which delivers a microscopic picture of the MOF surface/polymer interactions at the interface, would lead, in turn, to the understanding of the compatibility in MOF-based mixed-matrix membranes.

We have studied 1:1 inclusion complexes of two imidazole-based ionic liquids within β-cyclodextrin: 1-dodecyl-3-methylimidazolium and 1-butyl-3-methylimidazolium. By means of an adaptive biasing force scheme, we obtained the free energy profile along two different pathways, differing in the orientations of the head-to-tail vector with respect to the primary-secondary rim axis. Regarding 1-dodecyl-3-methylimidazolium, we found one minimum energy structure for each pathway, in which the hydrophobic tail remains embedded within the cyclodextrin, while the headgroup lies ∼11-12 Å from one of the rims; the structure where the polar head lies near the primary rim is the most stable. The analysis of the free energy of encapsulation of 1-butyl-3-methylimidazolium shows two minima for each insertion pathway, each of them associated with configurations where the imidazolium head lies close to one of the polar rims. As such, the most stable structure corresponds to one where the hydrophobic tail lies embedded within the cyclodextrin, while its head is localized near the secondary rim. The results are interpreted in terms of a simple model which captures the essential features that control the encapsulation process. A comparison with available experimental data is presented.

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.

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.

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.

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.

Using molecular dynamics experiments, we analyze the association of "Janus" 6-amino-6-deoxy-2O-carboxymethyl-beta-cyclodextrins (JCD) in aqueous solutions. In JCD dimers, the free energy associated with the primary-rim-secondary-rim docking shows a stable minimum of approximately -45 kcal mol(-1). Trimers in solution are also remarkably stable, exhibiting minimal distortions in their spatial and orientational distributions. The resulting geometrical docking shows the incipient characteristics of flexible nanotubes in solution, with eventual water interchange between the central channel and the bulk at the junctions between monomers. Structural and dynamical properties of the trapped water filling the nanotube are dictated to a large extent by the charge density at the rims.