We provide a general thermodynamic framework for the understanding of guest-induced structural transitions in hybrid organic-inorganic materials. The method is based on the analysis of experimental adsorption isotherms. It allows the determination of the free energy differences between host structures involved in guest-induced transitions, especially hard to obtain experimentally. We discuss the general case of adsorption in flexible materials and show how a few key quantities, such as pore volumes and adsorption affinities, entirely determine the phenomenology of adsorption, including the occurrence of structural transitions. On the basis of adsorption thermodynamics, we then propose a taxonomy of guest-induced structural phase transitions and the corresponding isotherms. In particular, we derive generic conditions for observing a double structural transition upon adsorption, often resulting in a two-step isotherm. Finally, we show the wide applicability and the robustness of the model through three case studies of topical hybrid organic-inorganic frameworks: the hysteretic hydrogen adsorption in Co(1,4-benzenedipyrazolate), the guest-dependent gate-opening in Cu(4,4'-bipyridine)(2,5-dihydroxybenzoate)2 and the CO2-induced "breathing" of hybrid material MIL-53.

The parameters of the anisotropic united atoms potential for linear alkanes proposed by Toxvaerd [S. Toxvaerd, J. Chem. Phys. 107, 5197 (1997)] have been optimized on the basis of selected equilibrium properties (vapor pressures, vaporization enthalpies, and liquid densities) of ethane, n-pentane, and n-dodecane. The optimized parameters for the CH2 and CH3 groups form a regular sequence with those of methane and the force centers are found between the carbon and hydrogen atoms, as expected. The resulting potential, called AUA4, has been compared with Toxvaerd’s potential (AUA3) by using several molecular simulation methods (Gibbs ensemble Monte Carlo, thermodynamic integration, and molecular dynamics). An investigation performed at temperatures ranging from 140 to 700 K and with various chain lengths up to 20 carbon atoms has shown AUA4 to provide systematic improvements of vapor pressures, vaporization enthalpies, and liquid densities for pure n-alkanes. Significant improvements have been also noticed on the critical temperatures of n-alkanes, estimated from coexistence density curves, and on the equilibrium properties of CO2–n-alkane binary mixtures. Self-diffusion coefficients of n-hexane, however, are slightly improved by the new potential, but still exceed experimental measurements at low temperature. As we have only optimized the intermolecular potential in the present study, it is suggested that further optimization of the intramolecular potentials of the anisotropic united atoms model could allow simultaneous prediction of thermodynamic properties and of transport coefficients, particularly in very dense liquids.

We performed ab initio calculations of the elastic constants of five flexible metal-organic frameworks (MOFs): MIL-53(Al), MIL-53(Ga), MIL-47, and the square and lozenge structures of DMOF-1. Tensorial analysis of the elastic constants reveals a highly anisotropic elastic behavior, some deformation directions exhibiting very low Young's modulus and shear modulus. This anisotropy can reach a 400:1 ratio between the most rigid and weakest directions, in stark contrast to the case of nonflexible MOFs such as MOF-5 and ZIF-8. In addition, we show that flexible MOFs can display extremely large negative linear compressibility. These results uncover the microscopic roots of stimuli-induced structural transitions in flexible MOFs, by linking the local elastic behavior of the material and its multistability.

The use of the osmotic thermodynamic model, combined with a series of methane and carbon dioxide gas adsorption experiments at various temperatures, has allowed to shed some new light on the fascinating phase behavior of flexible MIL-53(Al) metal-organic frameworks. A generic temperature-loading phase diagram has been derived; it is shown that the breathing effect in MIL-53 is a very general phenomenon, which should be observed in a limited temperature range regardless of the guest molecule. In addition, the previously proposed stress model for the structural transitions of MIL-53 is shown to be transferable from xenon to methane adsorption. The stress model also provides a theoretical framework for understanding the existence of lp/np phase mixtures at pressures close to the breathing transition pressure, without having to invoke an inhomogeneous distribution of the adsorbate in the porous sample.

We report a joint experimental and molecular simulation study of water condensation in silicalite-1 zeolite. A sample was synthesized using the fluoride route and was found to contain essentially no defects. A second sample synthesized using the hydroxide route was found to contain a small amount of silanol groups. The thermodynamics of water condensation was studied in these two samples, as well as in a commercial sample, in order to understand the effect of local defects on water adsorption. The molecular simulation study enabled us to qualitatively reproduce the experimentally observed condensation thermodynamics features. A shift and a rounding of the condensation transition was observed with an increasing hydrophilicity of the local defect, but the condensation transition was still observed above the water saturation vapor pressure P0. Both experiments and simulations agree on the fact that a small water uptake can be observed at very low pressure, but that the bulk liquid does not form from the gas phase below P0. The picture that emerges from the observed water condensation mechanism is the existence of a heterogeneous internal surface that is overall hydrophobic, despite the existence of hydrophilic "patches". This heterogeneous surface configuration is thermodynamically stable in a wide range of reduced pressures (from P/P0 = 0.2 to a few thousands), until the condensation transition takes place.

Gas adsorption in pores of flexible metal-organic frameworks (MOF) induces elastic deformation and structural transitions associated with stepwise expansion and contraction of the material, known as breathing transitions between large pore (lp) and narrow pore (np) phases. We present here a simple yet instructive model for the physical mechanism of this enigmatic phenomenon considering the adsorption-induced stress exerted on the material as a stimulus that triggers breathing transitions. The proposed model implies that the structural transitions in MOFs occur when the stress reaches a certain critical threshold. We showcase this model by drawing on the example of Xe adsorption in MIL-53 (Al) at 220 K, which exhibits two consecutive hysteretic breathing transitions between lp and np phases. We also propose an explanation for the experimentally observed coexistence of np and lp phases in MIL-53 materials.

We report grand canonical Monte Carlo simulations of the gas and liquid phase adsorption of water in silicalite-1 zeolite. Simple but effective models and simulation methods, found useful for studying gas adsorption in nanoporous materials, have been extended to describe the intrusion/extrusion cycle of water in this hydrophobic solid. The picture of water confined to hydrophobic spaces of nanoscopic dimensions that emerges from this study is one of a strongly depleted and highly inhomogeneous fluid.

Porous metal-organic frameworks (MOFs) are a topical class of materials that display an extremely large range of crystal structures and host-guest properties, potentially giving them a major impact in many areas of science and technology.