We report a joint experimental and molecular simulation study of water intrusion in silicalite-1 and ferrerite zeolites. The main conclusion of this study is that water condensation takes place through a genuine first-order phase transition, provided that the interconnected pores structure is 3-dimensional. In the extreme confinement situation (ferrierite zeolite), condensation takes place through a continuous transition, which is explained by a shift of both the first-order transition line and the critical point with increasing confinement. The present findings are at odds with the common belief that conventional phase transitions cannot take place in microporous solids such as zeolites. The most important features of the intrusion/extrusion process can be understood in terms of equilibrium thermodynamics considerations. We believe that these findings are very general for hydrophobic solids, i.e. for both nonwetting as well as wetting water–solid interface systems.

Calibration of forcefields for molecular simulation should account for the measurement uncertainty of the reference dataset and for the model inadequacy, i.e., the inability of the force-field/simulation pair to reproduce experimental data within their uncertainty range. In all rigour, the resulting uncertainty of calibrated force-field parameters is a source of uncertainty for simulation predictions. Various calibration strategies and calibration models within the Bayesian calibration/prediction framework are explored in the present article. In the case of Lennard-Jones potential for Argon, we show that prediction uncertainty for thermodynamical and transport properties, albeit very small, is larger than statistical simulation uncertainty.

We provide a comprehensive depiction of the behaviour of a nanodroplet of ≃20 water molecules confined in the pores of a series of 3D-connected isostructural zeolites with varying acidity, by means of molecular simulations. Both grand canonical Monte Carlo simulations using classical interatomic forcefields and first-principles Car–Parrinello molecular dynamics were used in order to characterise the behaviour of confined water by computing a range of properties, from thermodynamic quantities to electronic properties such as dipole moment, including structural and dynamical information. From the thermodynamic point of view, we have identified the all-silica zeolite as hydrophobic, and the cationic zeolites as hydrophilic; the condensation transition in the first case was demonstrated to be of first order. Furthermore, in-depth analysis of the dynamical and electronic properties of water showed that water in the hydrophobic zeolite behaves as a nanodroplet trying to close its hydrogen-bond network onto itself, with a few short-lived dangling OH groups, while water in hydrophilic zeolites “opens up” to form weak hydrogen bonds with the zeolite oxygen atoms. Finally, the dipole moment of confined water is studied and the contributions of water self-polarisation and the zeolite electric field are discussed.