Molecular Simulation of the Impact of Defects on Electrolyte Intrusion in Zeolites.

We have investigated through molecular simulation the intrusion of electrolytes in two representative pure-silica zeolites, silicalite-1 and chabazite, in which point defects were introduced in varying amounts. We distinguish between two types of defects, considering either "weak" or "strong" silanol nest defects, resulting in different hydration behaviors. In the presence of weak defects, the hydration process occurs through a homogeneous nucleation process, while with strong defects, we observe an initial adsorption followed by a filling of the nanoporous volume at a higher pressure. However, we show that electrolytes do not penetrate the zeolites, and these defects appear to have only marginal influence on the thermodynamics of electrolyte intrusion. While replacing pure water by the electrolyte solution shifts the intrusion pressure toward higher values because of the drop of water saturation vapor pressure, an increase in hydrophilicity of the framework due to point defects has the opposite effect, showing that controlling the amount of defects in zeolites is crucial for storage energy applications.

Alchemical Osmostat for Monte Carlo Simulation: Sampling Aqueous Electrolyte Solution in Open Systems.

Molecular simulations involving electrolytes are usually performed at a fixed amount of salt ions in the simulation box, reproducing macroscopic concentration. Although this statement is valid in the bulk, the concentration of an electrolyte confined in nanoporous materials such as MOFs or zeolites is greatly affected and remains a priori unknown. The nanoporous material in equilibrium with the bulk electrolyte exchange water and ions at a given chemical potential Δµ in the semi-grand-canonical ensemble, that must be calibrated in order to determine the concentration in the nanoporous material. In this work, we propose an algorithm based on nonequilibrium candidate Monte Carlo (NCMC) moves to ultimately perform MC simulations in contact with a saline reservoir. First, we adapt the Widom insertion technique to calibrate the chemical potential by alchemically transmuting water molecules into ions by using NCMC moves. The chemical potential defines a Monte Carlo osmostat in the semi-grand-constant volume and temperature ensemble (Δµ, N, V, T) to be added in a Monte Carlo simulation where the number of ions fluctuates. In order to validate the method, we adapted the NCMC move to determine the free energy of water solvation and subsequently explore thermodynamics of electrolyte solvation at infinite dilution in water. Finally, we implemented the osmostat in MC simulations initialized with bulk water that are driven toward electrolyte solutions of similar concentration as the saline reservoir. Our results demonstrate that alchemical osmostat for MC simulation is a promising tool for use to sample electrolyte insertion in nanoporous materials.

Forced intrusion of water and aqueous solutions in microporous materials: from fundamental thermodynamics to energy storage devices.

We review the high pressure forced intrusion studies of water in hydrophobic microporous materials such as zeolites and MOFs, a field of research that has emerged some 15 years ago and is now very active. Many of these studies are aimed at investigating the possibility of using these systems as energy storage devices. A series of all-silica zeolites (zeosil) frameworks were found suitable for reversible energy storage because of their stability with respect to hydrolysis after several water intrusion-extrusion cycles. Several microporous hydrophobic zeolite imidazolate frameworks (ZIFs) also happen to be quite stable and resistant towards hydrolysis and thus seem very promising for energy storage applications. Replacing pure water by electrolyte aqueous solutions enables to increase the stored energy by a factor close to 3, on account of the high pressure shift of the intrusion transition. In addition to the fact that aqueous solutions and microporous silica materials are environmental friendly, these systems are thus becoming increasingly interesting for the design of new energy storage devices. This review also addresses the theoretical approaches and molecular simulations performed in order to better understand the experimental behavior of nano-confined water. Molecular simulation studies showed that water condensation takes place through a genuine first-order phase transition, provided that the interconnected pores structure is 3-dimensional and sufficiently open. In an extreme confinement situations such as in ferrierite zeosil, condensation seem to take place through a continuous supercritical crossing from a diluted to a dense fluid, on account of the fact that the first-order transition line is shifted to higher pressure, and the confined water critical point is correlatively shifted to lower temperature. These molecular simulation studies suggest that the most important features of the intrusion/extrusion process can be understood in terms of equilibrium thermodynamics considerations.

Interplay between defects, disorder and flexibility in metal-organic frameworks.

Metal-organic frameworks are a novel family of chemically diverse materials, which are of interest across engineering, physics, chemistry, biology and medicine-based disciplines. Since the development of the field in its current form more than two decades ago, priority has been placed on the synthesis of new structures. However, more recently, a clear trend has emerged in shifting the emphasis from material design to exploring the chemical and physical properties of structures already known. In particular, although such nanoporous materials were traditionally seen as rigid crystalline structures, there is growing evidence that large-scale flexibility, the presence of defects and long-range disorder are not the exception in metal-organic frameworks, but the rule. Here we offer some perspective into how these concepts are perhaps inescapably intertwined, highlight recent advances in our understanding and discuss how a consideration of the interfaces between them may lead to enhancements of the materials' functionalities.

Adsorption deformation of microporous composites.

We study here the behavior of flexible adsorbent materials, or soft porous crystals, when used in practical applications as nanostructured composites such as core-shell particles or mixed matrix membranes. Based on simple models and the well-established laws of elasticity, we demonstrate how the presence of a binder results in an attenuation of the adsorption-induced stress and deformation. In the case where the adsorbent undergoes adsorption-induced structural transitions, such as the gate opening phenomenon occurring in some metal-organic frameworks, we show that the presence of the binder will result in shifts of the adsorption-induced transition pressures.

Softening upon Adsorption in Microporous Materials: A Counterintuitive Mechanical Response.

We demonstrate here that microporous materials can exhibit softening upon adsorption of guest molecules, at low to intermediate pore loading, in parallel to the pore shrinking that is well-known in this regime. This novel and counterintuitive mechanical response was observed through molecular simulations of both model pore systems (such as slit pore) and real metal-organic frameworks. It is contrary to common belief that adsorption of guest molecules necessarily leads to stiffening due to increased density, a fact that we show is the high-loading limit of a more complex behavior: a nonmonotonic softening-then-stiffening.

Hydrothermal Breakdown of Flexible Metal-Organic Frameworks: A Study by First-Principles Molecular Dynamics.

Flexible metal-organic frameworks, also known as soft porous crystals, have been proposed for a vast number of technological applications, because they respond by large changes in structure and properties to small external stimuli, such as adsorption of guest molecules and changes in temperature or pressure. While this behavior is highly desirable in applications such as sensing and actuation, their extreme flexibility can also be synonymous with decreased stability. In particular, their performance in industrial environments is limited by a lack of stability at elevated temperatures and in the presence of water. Here, we use first-principles molecular dynamics to study the hydrothermal breakdown of soft porous crystals. Focusing on the material MIL-53(Ga), we show that the weak point of the structure is the bond between the metal center and the organic linker and elucidate the mechanism by which water lowers the activation free energy for the breakdown. This allows us to propose strategies for the synthesis of MOFs with increased heat and water stability.

What makes zeolitic imidazolate frameworks hydrophobic or hydrophilic? The impact of geometry and functionalization on water adsorption.

We demonstrate, by means of Grand Canonical Monte Carlo simulation on different members of the ZIF family, how topology, geometry, and linker functionalization drastically affect the water adsorption properties of these materials, tweaking the ZIF materials from hydrophobic to hydrophilic. We show that adequate functionalization of the linkers allows one to tune the host-guest interactions, even featuring dual amphiphilic materials whose pore space features both hydrophobic and hydrophilic regions. Starting from an initially hydrophobic material (ZIF-8), various degrees of hydrophilicity could be obtained, with a gradual evolution from a type V adsorption isotherm in the liquid phase to a type I isotherm in the gas phase. This behavior is similar to what was described earlier in families of hydrophobic all-silica zeolites, with hydrophilic "defects" of various strength, such as silanol nests or the presence of extra-framework cations.

Metal-organic frameworks with wine-rack motif: what determines their flexibility and elastic properties?

We present here a framework for the analysis of the full tensors of second-order elastic constants of metal-organic frameworks, which can be obtained by ab initio calculations. We describe the various mechanical properties one can derive from such tensors: directional Young's modulus, shear modulus, Poisson ratio, and linear compressibility. We then apply this methodology to four different metal-organic frameworks displaying a wine-rack structure: MIL-53(Al), MIL-47, MIL-122(In), and MIL-140A. From these results, we shed some light into the link between mechanical properties, geometric shape, and compliance of the framework of these porous solids. We conclude by proposing a simple criterion to assess the framework compliance, based on the lowest eigenvalue of its second-order elastic tensor.

Adsorption induced transitions in soft porous crystals: an osmotic potential approach to multistability and intermediate structures.

Soft porous crystals are flexible metal-organic frameworks that respond to physical stimuli (temperature, pressure, and gas adsorption) by large changes in their structure and unit cell volume. We propose here a thermodynamic treatment, based on the osmotic ensemble, of the interplay between guest adsorption and host deformation, where the bare host material can undergo elastic deformation, as well as structural transitions between metastable phases in the case of a multistable material. We show that in addition to structural transitions between metastable phases of bistable or multistable host frameworks, a new guest-stabilized host phase can be created when the size of the adsorbate is larger than the empty material's pore size. We then confront the findings of our approach with experimental data for systems exhibiting phenomena such as gate opening and breathing.

Investigating the Pressure-Induced Amorphization of Zeolitic Imidazolate Framework ZIF-8: Mechanical Instability Due to Shear Mode Softening.

We provide the first molecular dynamics study of the mechanical instability that is the cause of pressure-induced amorphization of zeolitic imidazolate framework ZIF-8. By measuring the elastic constants of ZIF-8 up to the amorphization pressure, we show that the crystal-to-amorphous transition is triggered by the mechanical instability of ZIF-8 under compression, due to shear mode softening of the material. No similar softening was observed under temperature increase, explaining the absence of temperature-induced amorphization in ZIF-8. We also demonstrate the large impact of the presence of adsorbate in the pores on the mechanical stability and compressibility of the framework, increasing its shear stability. This first molecular dynamics study of ZIF mechanical properties under variations of pressure, temperature, and pore filling opens the way to a more comprehensive understanding of their mechanical stability, structural transitions, and amorphization.

Anisotropic elastic properties of flexible metal-organic frameworks: how soft are soft porous crystals?

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.

Understanding adsorption-induced structural transitions in metal-organic frameworks: from the unit cell to the crystal.

Breathing transitions represent recently discovered adsorption-induced structural transformations between large-pore and narrow-pore conformations in bi-stable metal-organic frameworks such as MIL-53. We present a multiscale physical mechanism of the dynamics of breathing transitions. We show that due to interplay between host framework elasticity and guest molecule adsorption, these transformations on the crystal level occur via layer-by-layer shear. We construct a simple Hamiltonian that describes the physics of host-host and host-guest interactions on the level of unit cells and reduces to one effective dimension due to the long-range elastic cell-cell interactions. We then use this Hamiltonian in Monte Carlo simulations of adsorption-desorption cycles to study how the behavior of unit cells is linked to the transition mechanism at the crystal level through three key physical parameters: the transition energy barrier, the cell-cell elastic coupling, and the system size.

How can a hydrophobic MOF be water-unstable? Insight into the hydration mechanism of IRMOFs.

We report an ab initio molecular dynamics study of the hydration process in a model IRMOF material. At low water content (one molecule per unit cell), water physisorption is observed on the zinc cation but the free⇄bound equilibrium strongly favors the free state. This is consistent with the hydrophobic nature of the host matrix and its type-V isotherm observed in a classical Monte Carlo simulation. At higher loading, a water cluster can be formed at the Zn(4)O site and this is shown to stabilize the water-bound state. This structure rapidly transforms into a linker-displaced state, where water has fully displaced one arm of a linker and which corresponds to the loss of the material's fully ordered structure. Thus an overall hydrophobic MOF material can also become water unstable, a feature that has not been fully understood until now.

Experiment and theory of low-pressure nitrogen adsorption in organic layers supported or grafted on inorganic adsorbents: toward a tool to characterize surfaces of hybrid organic/inorganic systems.

We report experimental nitrogen adsorption isotherms of organics-coated silicas, which exhibit a low-pressure desorption branch that does not meet the adsorption branch upon emptying of the pores. To address the physical origin of such a hysteresis loop, we propose an equilibrium thermodynamic model that enables one to explain this phenomenon. The present model assumes that, upon adsorption, a small amount of nitrogen molecules penetrate within the organic layer and reach adsorption sites that are located on the inorganic surface, between the grafted or adsorbed organic molecules. The number of accessible adsorption sites thus varies with the increasing gas pressure, and then we assume that it stays constant upon desorption. Comparison with experimental data shows that our model captures the features of nitrogen adsorption on such hybrid organic/inorganic materials. In particular, in addition to predicting the shape of the adsorption isotherm, the model is able to estimate, with a reasonable number of adjustable parameters, the height of the low-pressure hysteresis loop and to assess in a qualitative fashion the local density of the organic chains at the surface of the material.

Predicting mixture coadsorption in soft porous crystals: experimental and theoretical Study of CO2/CH4 in MIL-53(Al).

We present a synergistic experimental and theoretical study of CO(2)/CH(4) mixture coadsorption in breathing metal-organic framework MIL-53(Al). Mixture adsorption experiments were performed and their results were analyzed by comparing them to predictions made from pure-component adsorption data using the Osmotic Framework Adsorption Solution Theory (OFAST). This analytical model, fully validated for the first time, was then used to predict coadsorption properties as a function of temperature, pressure, and mixture composition. The phase diagrams obtained show a surprising non-monotonic behavior.