Some porous crystalline solids change their structure upon guest inclusion. Unlocking the potential of these solids for a wide variety of applications requires full characterisation of the response to adsorption and the underlying framework–guest interactions. Here, we introduce an approach to understanding gas uptake in porous metal-organic frameworks (MOFs) by loading liquefied gases at GPa pressures inside the Zn-based framework ZIF-8. An integrated experimental and computational study using high-pressure crystallography, grand canonical Monte Carlo (GCMC) and periodic DFT simulations has revealed six symmetry-independent adsorption sites within the framework and a transition to a high-pressure phase. The cryogenic high-pressure loading method offers a different approach to obtaining atomistic detail on guest molecules. The GCMC simulations provide information on interaction energies of the adsorption sites allowing to classify the sites by energy. DFT calculations reveal the energy barrier of the transition to the high-pressure phase. This combination of techniques provides a holistic approach to understanding both structural and energetic changes upon adsorption in MOFs.

The incorporation of coordinatively unsaturated metal sites (cus's), also known as open metal sites, into metal− organic frameworks (MOFs), significantly enhances the uptake of certain gases, such as CO 2 and CH 4 , especially at low loadings when fluid−framework interactions play the predominant role. However, due to the considerably enhanced, localized guest interactions with the cus's, it remains a challenge to predict correctly adsorption isotherms and mechanisms in MOFs with cus's using grand-canonical Monte Carlo (GCMC) simulations based on generic classical force fields. To address this problem, we carefully investigated several well-established semiempirical model potentials and used a multiobjective genetic algorithm to parametrize them using accurate ab initio data as reference. The Carra− Konowalow potential, a modified Buckingham potential, in combination with the MMSV potential for the cus's gives not only adsorption isotherms in very good agreement with experiments but also correctly captures the adsorption mechanisms, including adsorption on the cus's, for CO 2 in CPO-27-Mg and CH 4 in CuBTC. Moreover, the parameters obtained also give quantitative predictions of CH 4 adsorption in PCN-14, another MOF with Cu cus's, which is an important step for developing transferable force fields that reliably predict adsorption in MOFs with cus's.

Grace under pressure: In the first high‐pressure crystallographic study on the metal–organic framework MOF‐5, increasing pressure initially results in the pressure‐transmitting fluid being squeezed into the pores. Further pressure increase causes a large reduction in pore content as solvent is evacuated from the pores, until a complete loss of crystallinity is observed at pressures above 3.2 GPa.

A new method to obtain improved structural parameters by supplementing gas-phase electron diffraction (GED) data with restraints based on the results of ab initio calculations is proposed. The procedure involves the use of ab initio parameters with estimated uncertainties as additional observations; this allows previously fixed parameters to refine, with all geometrical parameters included in the final refinement. The refinement of the molecular structure of 2,5-dichloropyrimidine is used as an example to illustrate the principle of this technique. In this simple case, the effects are not very great, but this new approach allowed refinement of all structural parameters. The nine independent structural parameters (r R structure) were found to be: 6)-H( 10)] ) 109.9(12) pm, ∠[N(1)C( 2)N(3)] ) 127.9(4)°, ∠[C(2)N(3)C( 4)] ) 116.3(7)°, and ∠[N(3)C(4)H( 8)] ) 117.2(5)°. All structural parameters were found to be in good agreement with both ab initio and crystallographic values, which are presented for comparison.

Ab initio molecular dynamics (AIMD) simulations have been used to predict structural transitions of the breathing metal-organic framework (MOF) MIL-53(Sc) in response to changes in temperature over the range 100-623 K and adsorption of CO2 at 0-0.9 bar at 196 K. The method has for the first time been shown to predict successfully both temperature-dependent structural changes and the structural response to variable sorbate uptake of a flexible MOF. AIMD employing dispersion-corrected density functional theory accurately simulated the experimentally observed closure of MIL-53(Sc) upon solvent removal and the transition of the empty MOF from the closed-pore phase to the very-narrow-pore phase (symmetry change from P2(1)/c to C2/c) with increasing temperature, indicating that it can directly take into account entropic as well as enthalpic effects. We also used AIMD simulations to mimic the CO2 adsorption of MIL-53(Sc) in silico by allowing the MIL-53(Sc) framework to evolve freely in response to CO2 loadings corresponding to the two steps in the experimental adsorption isotherm. The resulting structures enabled the structure determination of the two CO2-containing intermediate and large-pore phases observed by experimental synchrotron X-ray diffraction studies with increasing CO2 pressure; this would not have been possible for the intermediate structure via conventional methods because of diffraction peak broadening. Furthermore, the strong and anisotropic peak broadening observed for the intermediate structure could be explained in terms of fluctuations of the framework predicted by the AIMD simulations. Fundamental insights from the molecular-level interactions further revealed the origin of the breathing of MIL-53(Sc) upon temperature variation and CO2 adsorption. These simulations illustrate the power of the AIMD method for the prediction and understanding of the behavior of flexible microporous solids.