Water Capture Mechanisms at Zeolitic Imidazolate Framework Interfaces

Current Status of the MB-pol Data-Driven Many-Body Potential for Predictive Simulations of Water Across Different Phases

Developing a molecular-level understanding of the properties of water is central to numerous scientific and technological applications. However, accurately modeling water through computer simulations has been a significant challenge due to the complex nature of the hydrogen-bonding network that water molecules form under different thermodynamic conditions. This complexity has led to over five decades of research and many modeling attempts. The introduction of the MB-pol data-driven many-body potential energy function marked a significant advancement toward a universal molecular model capable of predicting the structural, thermodynamic, dynamical, and spectroscopic properties of water across all phases. By integrating physics-based and data-driven (i.e., machine-learned) components, which correctly capture the delicate balance among different many-body interactions, MB-pol achieves chemical and spectroscopic accuracy, enabling realistic molecular simulations of water, from gas-phase clusters to liquid water and ice. In this review, we present a comprehensive overview of the data-driven many-body formalism adopted by MB-pol, highlight the main results and predictions made from computer simulations with MB-pol to date, and discuss the prospects for future extensions to data-driven many-body potentials of generic and reactive molecular systems.

Entropy of Liquid Water as Predicted by the Two-Phase Thermodynamic Model and Data-Driven Many-Body Potentials

We investigate the entropy of liquid water at ambient conditions using the two-phase thermodynamic (2PT) model, applied to both common pairwise-additive water models and the MB-pol and MB-pol(2023) data-driven many-body potentials. Our simulations demonstrate that the 2PT model yields entropy values in semiquantitative agreement with experimental data when using MB-pol and MB-pol(2023). Additionally, our analyses indicate that the entropy values predicted by pairwise-additive water models may benefit from error compensation between the inherent approximations of the 2PT model and the known limitations of these water models in describing many-body interactions. Despite its approximate nature, the simplicity of the 2PT model makes it a valuable tool for estimating relative entropy changes of liquid water across various environments, especially when combined with water models that provide a consistently robust representation of the structural, thermodynamic, and dynamical properties of liquid water.

The Surface of Colloidal Metal-Organic Framework Nanoparticles Revealed by Vibrational Sum Frequency Scattering Spectroscopy

Solvation shells strongly influence the interfacial chemistry of colloidal systems, from the activity of proteins to the colloidal stability and catalysis of nanoparticles. Despite their fundamental and practical importance, solvation shells have remained largely undetected by spectroscopy. Furthermore, their ability to assemble at complex but realistic interfaces with heterogeneous and rough surfaces remains an open question. Here, we apply vibrational sum frequency scattering spectroscopy (VSFSS), an interface-specific technique, to colloidal nanocrystals with porous metal-organic frameworks (MOFs) as a case study. Due to the porous nature of the solvent-particle boundary, MOF particles challenge conventional models of colloidal and interfacial chemistry. Their multiweek colloidal stability in the absence of conventional surface ligands suggests that stability may arise in part from solvation forces. Spectra of colloidally stable Zn(2-methylimidazolate)2 (ZIF-8) in polar solvents indicate the presence of ordered solvation shells, solvent-metal binding, and spontaneous ordering of organic bridging linkers within the MOF. These findings help explain the unexpected colloidal stability of MOF colloids, while providing a roadmap for applying VSFSS to wide-ranging colloidal nanocrystals in general.

Flat-Histogram Monte Carlo Simulation of Water Adsorption in Metal-Organic Frameworks

Molecular simulations of water adsorption in porous materials often converge slowly due to sampling bottlenecks that follow from hydrogen bonding and, in many cases, the formation of water clusters. These effects may be exacerbated in metal-organic framework (MOF) adsorbents, due to the presence of pore spaces (cages) that promote the formation of discrete-size clusters and hydrophobic effects (if present), among other reasons. In Grand Canonical Monte Carlo (MC) simulations, these sampling challenges are typically manifested by low MC acceptance ratios, a tendency for the simulation to become stuck in a particular loading state (i.e., macrostates), and the persistence of specific clusters for long periods of the simulation. We present simulation strategies to address these sampling challenges, by applying flat-histogram MC (FHMC) methods and specialized MC move types to simulations of water adsorption. FHMC, in both Transition-matrix and Wang-Landau forms, drives the simulation to sample relevant macrostates by incorporating weights that are self-consistently adjusted throughout the simulation and generate the macrostate probability distribution (MPD). Specialized MC moves, based on aggregation-volume bias and configurational bias methods, separately address low acceptance ratios for basic MC trial moves and specifically target water molecules in clusters; in turn, the specialized MC moves improve the efficiency of generating new configurations which is ultimately reflected in improved statistics collected by FHMC. The combined strategies are applied to study the adsorption of water in CuBTC and ZIF-8 at 300 K, through examination of the MPD and the adsorption isotherm generated by histogram reweighting. A key result is the appearance of nontrivial oscillations in the MPD, which we show to be associated with water clusters in the adsorption system. Additionally, we show that the probabilities of certain clusters become similar in value near the boundaries of the isotherm hysteresis loop, indicating a strong connection between cluster formation/destruction and the thermodynamic limits of stability. For a hydrophobic MOF, the FHMC results show that the phase transition from low density to high density is suppressed to water pressure far above the bulk-fluid saturation pressure; this is consistent with results presented elsewhere. We also compare our FHMC simulation isotherm to one measured by a different technique but with ostensibly the same molecular interactions and comment on observed differences and the need for follow-up work. The simulation strategies presented here can be applied to the simulation of water in other MOFs using heuristic guidelines laid out in our text, which should facilitate the more consistent and efficient simulation of water adsorption in porous materials in future applications.

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal-organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal-organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

Molecular driving forces for water adsorption in MOF-808: A comparative analysis with UiO-66

Metal-organic frameworks (MOFs), with their unique porous structures and versatile functionality, have emerged as promising materials for the adsorption, separation, and storage of diverse molecular species. In this study, we investigate water adsorption in MOF-808, a prototypical MOF that shares the same secondary building unit (SBU) as UiO-66, and elucidate how differences in topology and connectivity between the two MOFs influence the adsorption mechanism. To this end, molecular dynamics simulations were performed to calculate several thermodynamic and dynamical properties of water in MOF-808 as a function of relative humidity (RH), from the initial adsorption step to full pore filling. At low RH, the μ3-OH groups of the SBUs form hydrogen bonds with the initial water molecules entering the pores, which triggers the filling of these pores before the μ3-OH groups in other pores become engaged in hydrogen bonding with water molecules. Our analyses indicate that the pores of MOF-808 become filled by water sequentially as the RH increases. A similar mechanism has been reported for water adsorption in UiO-66. Despite this similarity, our study highlights distinct thermodynamic properties and framework characteristics that influence the adsorption process differently in MOF-808 and UiO-66.

Co(II)-doped hybrid Zn(II) tetraborate complexes, [ZnxCo(1-x)(1,3-dap)B4O7] (1,3-dap = 1,3-diaminopropane): BET analysis and N2/H2O/D2O adsorption studies

A series of mono/bimetallic isostructural hybrid tetraborates of the general formula [ZnxCo(1-x)(1,3-dap)B4O7] has been prepared using a solvothermal method. Their adsorption/desorption curves for H2O and D2O demonstrate that these materials have a stronger affinity for H2O than for D2O and enrich the D2O content of D2O/H2O mixtures.

Tannic acid and carboxymethyl chitosan-based multi-functional double-layered hydrogel with pH-stimulated response behavior for smart real-time infection monitoring and wound treatment

The dramatic increase of drug-resistant pathogenic bacteria has seriously effect on human health, appealing the needs of developing theranostic platforms with stimuli-responsive materials to realize the accurate bacterial diagnostics and therapeutics. Herein, a tannic acid and carboxymethyl chitosan-based multifunctional ZIF-90@i-PPOPs-phenol red double-layered hydrogel with stimuli-responsiveness and antibacterial activity was fabricated. The inner layer hydrogel (ZIF-90@i-PPOPs-based TFC hydrogels) was fabricated based on ZIF-90@i-PPOPs, integrate tannic acid and carboxymethyl chitosan linked by formylphenylboronic acid (FPBA), which exhibited outstanding injectable, biodegradability and antibacterial activity. The outer layer hydrogel (PR@PAM hydrogels) were constructed from polyacrylamide (PAM) and pH indicator phenol red, owning porous structure and excellent tissue adhesion. Due to the weakly acidic microenvironment within wound, the inner-layer hydrogel was stimulus-responsively decomposed, resulting in the accurate delivery of the positively charged ZIF-90@i-PPOPs to the lesion site to capture and kill bacteria by enhanced Zn2+ and ROS release. Meantime, the outer-layer hydrogel could real-timely monitor the pH changes to evaluate the wound recovery status. These double-layered hydrogels possessed precisely pH monitoring capacity, excellent antibacterial ability and negligible side effect to normal tissue in vivo, implying the high potential of the suggested hydrogels as theranostic platform for antibacterial treatment.

Molecular Insights into the Influence of Ions on the Water Structure. I. Alkali Metal Ions in Solution

In this study, we explore the impact of alkali metal ions (Li+, Na+, K+, Rb+, and Cs+) on the hydration structure of water using molecular dynamics simulations carried out with MB-nrg potential energy functions (PEFs). Our analyses include radial distribution functions, coordination numbers, dipole moments, and infrared spectra of water molecules, calculated as a function of solvation shells. The results collectively indicate a highly local influence of all of the alkali metal ions on the hydrogen-bond network established by the surrounding water molecules, with the smallest and most densely charged Li+ ion exerting the most pronounced effect. Remarkably, the MB-nrg PEFs demonstrate excellent agreement with available experimental data for the position and size of the first solvation shells, underscoring their potential as predictive models for realistic simulations of ionic aqueous solutions across various thermodynamic conditions and environments.

Accuracy of TIP4P/2005 and SPC/Fw Water Models

Water is used as the main solvent in model systems containing bioorganic molecules. Choosing the right water model is an important step in the study of the biophysical and biochemical processes that occur in cells. In the present work, we perform molecular dynamics simulations using two distinct force fields for water: the rigid model TIP4P/2005, where only intermolecular interactions are considered, and the flexible model SPC/Fw, where intramolecular interactions are also taken into account. The simulations aim to determine the effect of the inclusion of intramolecular interactions on the accuracy of calculated properties of bulk water (density and thermal expansion coefficient, self-diffusion coefficients, shear viscosity, radial distribution functions, and dielectric constant), as compared to experimental results, over a temperature range between 250 and 370 K. We find that the results of the rigid model present the smallest deviations relative to experiments for most of the calculated quantities, except for the shear viscosity of supercooled water and the water dielectric constant, where the flexible model presents better agreement with experiments.

Visualizing Water Monomers and Chiral OH-(H2O) Complexes Infiltrated in a Macroscopic Hydrophobic Teflon Matrix

Insights into the interaction of fluoroalkyl groups with water are crucial to understanding the polar hydrophobicity of fluorinated compounds, such as Teflon. While an ordered hydrophobic-like 2D water layer has been demonstrated to be present on the surface of macroscopically hydrophobic fluorinated polymers, little is known about how the water infiltrates into the Teflon and what is the molecular structure of the water infiltrated into the Teflon. Using highly sensitive femtosecond sum frequency generation vibrational spectroscopy (SFG-VS), we observe for the first time that monomeric H2O and chiral OH-(H2O) complexes are present in macroscopically hydrophobic Teflon. The species are inhomogeneously distributed inside the Teflon matrix and at the Teflon surface. No water clusters or single-file water "wires" are observed in the matrix. SFG free induction decay (SFG-FID) experiments demonstrate that the OH oscillators of physically absorbed molecular water at the surface dephase on the time scale of <230 fs, whereas the water monomers and hydrated hydroxide ions infiltrated in the Teflon matrix dephase much more slowly (680-830 fs), indicating that the embedded monomeric H2O and OH-(H2O) complexes are decoupled from the outer environment. Our findings can well interpret ultrafast water permeation through fluorous nanochannels and the charging mechanism of Teflon, which may tailor the desired applications of organofluorines.

Effects of Defects and Presence of Open-Metal Sites on the Structure and Dynamics of Water in Hydrophobic Zeolitic Imidazolate Frameworks

Most of the chemistry in nanoporous materials with small pore sizes and windows takes place on the outer surface, which is in direct contact with the substrate/solvent, rather than within the pores and channels. Here, we report the results of our comprehensive atomistic molecular dynamics (MD) simulations to decipher the interaction of water with a realistic finite ∼5.1 nm nanoparticle (NP) model of ZIF-8, with edges containing undercoordinated Zn metal sites, vs a conventionally employed pristine crystalline bulk (CB) model. The hydrophobic interior surface of the CB model imparts significant dynamical behavior on water molecules with (i) increasing diffusivity from the surface toward the center of the pores and (ii) confined water, at low concentration, showing similar diffusivity to that of the bulk water. On the other hand, water molecules adsorbed on the surface of the NP model exhibit a range of characteristics, including "coordinated", "confined", and "bulk-like" behavior. Some of the water molecules form coordinative bonds with the undercoordinated Zn metal centers and act as nucleation sites for the water droplets to form, facilitating diffusion into the pores. However, diffusion of water molecules is limited to the areas near the surface and not all the way to the core of the NP model. Our atomistic MD simulations provide insights into the stability of ZIFs in aqueous solutions despite hydrolysis of their outer surface. Such insights are helpful in designing more robust nanoporous materials for applications in humid environments.

Elucidating the Competitive Adsorption of H2O and CO2 in CALF-20: New Insights for Enhanced Carbon Capture Metal-Organic Frameworks

In light of the pressing need for efficient carbon capture solutions, our study investigates the simultaneous adsorption of water (H2O) and carbon dioxide (CO2) as a function of relative humidity in CALF-20, a highly scalable and stable metal-organic framework (MOF). Advanced computer simulations reveal that due to their similar interactions with the framework, H2O and CO2 molecules compete for the same binding sites, occupying similar void regions within the CALF-20 pores. This competition results in distinct thermodynamic and dynamical behaviors of H2O and CO2 molecules, depending on whether one or both guest species are present. Notably, the presence of CO2 molecules forces the H2O molecules to form more connected hydrogen-bond networks within smaller regions, slowing water reorientation dynamics and decreasing water entropy. Conversely, the presence of water speeds up the reorientation of CO2 molecules, decreases the CO2 entropy, and increases the propensity for CO2 to be adsorbed within the framework due to stronger water-mediated interactions. Due to the competition for the same void spaces, both H2O and CO2 molecules exhibit slower diffusion when molecules of the other guest species are present. These findings offer valuable strategies and insights into enhancing the differential affinity of H2O and CO2 for MOFs specifically designed for carbon capture applications.

Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models

In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.

Ultrafast Water H-Bond Rearrangement in a Metal-Organic Framework Probed by Femtosecond Time-Resolved Infrared Spectroscopy

We investigated the water H-bond network and its dynamics in Ni₂Cl₂BTDD, a prototypical MOF for atmospheric water harvesting, using linear and ultrafast IR spectroscopy. Utilizing isotopic labeling and infrared spectroscopy, we found that water forms an extensive H-bonding network in Ni₂Cl₂BTDD. Further investigation with ultrafast spectroscopy revealed that water can reorient in a confined cone up to ∼50° within 1.3 ps. This large angle reorientation indicates H-bond rearrangement, similar to bulk water. Thus, although the water H-bond network is confined in Ni₂Cl₂BTDD, different from other confined systems, H-bond rearrangement is not hindered. The picosecond H-bond rearrangement in Ni₂Cl₂BTDD corroborates its reversibility with minimal hysteresis in water sorption.

Structure and thermodynamics of water adsorption in NU-1500-Cr

Metal-organic frameworks (MOFs) are a class of materials with diverse chemical and structural properties, and have been shown to effectively adsorb various types of guest molecules. The mechanism of water adsorption in NU-1500-Cr, a high-performance atmospheric water harvesting MOF, is investigated using a combination of molecular dynamics simulations and infrared spectroscopy. Calculations of thermodynamic and dynamical properties of water as a function of relative humidity allow for following the adsorption process from the initial hydration stage to complete filling of the MOF pores. Initial hydration begins at the water molecules that saturate the open Cr³⁺ sites of the framework, which is then followed by the formation of water chains that extend along the channels connecting the hexagonal pores of the framework. Water present in these channels gradually coalesces and fills the hexagonal pores sequentially after the channels are completely hydrated. The development of hydrogen-bond networks inside the MOF pores as a function of relative humidity is characterized at the molecular level using experimental and computational infrared spectroscopy. A detailed analysis of the OH-stretch vibrational band indicates that the low-frequency tail stems from strongly polarized hydrogen-bonded water molecules, suggesting the presence of some structural disorder in the experimental samples. Strategies for designing efficient water harvesting MOFs are also proposed based on the mechanism of water adsorption in NU-1500-Cr.

Imaging Orientation of a Single Molecular Hierarchical Self-Assembled Sheet: The Combined Power of a Vibrational Sum Frequency Generation Microscopy and Neural Network

In this work, we determined the tilt angles of molecular units in hierarchical self-assembled materials on a single-sheet level, which were not available previously. This was achieved by developing a fast line-scanning vibrational sum frequency generation (VSFG) hyperspectral imaging technique in combination with neural network analysis. Rapid VSFG imaging enabled polarization resolved images on a single sheet level to be measured quickly, circumventing technical challenges due to long-term optical instability. The polarization resolved hyperspectral images were then used to extract the supramolecular tilt angle of a self-assembly through a set of spectra-tilt angle relationships which were solved through neural network analysis. This unique combination of both novel techniques offers a new pathway to resolve molecular level structural information on self-assembled materials. Understanding these properties can further drive self-assembly design from a bottom-up approach for applications in biomimetic and drug delivery research.

Prospects for Simultaneously Capturing Carbon Dioxide and Harvesting Water from Air

Mitigation of anthropogenic climate change is expected to require large-scale deployment of carbon dioxide removal strategies. Prominent among these strategies is direct air capture with sequestration (DACS), which encompasses the removal and long-term storage of atmospheric CO₂  by purely engineered means. Because it does not require arable land or copious amounts of freshwater, DACS is already attractive in the context of sustainable development, but opportunities to improve its sustainability still exist. Leveraging differences in the chemistry of CO₂  and water adsorption within porous solids, here, the prospect of simultaneously removing water alongside CO₂  in direct air capture operations is investigated. In many cases, the co-adsorbed water can be desorbed separately from chemisorbed CO₂  molecules, enabling efficient harvesting of water from air. Depending upon the material employed and process conditions, the desorbed water can be of sufficiently high purity for industrial, agricultural, or potable use and can thus improve regional water security. Additionally, the recovered water can offset a portion of the costs associated with DACS. In this Perspective, molecular- and process-level insights are combined to identify routes toward realizing this nascent yet enticing concept.

Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks

Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.

Revealing the Molecular Physics of Lattice Self-Assembly by Vibrational Hyperspectral Imaging

Lattice self-assemblies (LSAs), which mimic protein assemblies, were studied using a new nonlinear vibrational imaging technique called vibrational sum-frequency generation (VSFG) microscopy. This technique successfully mapped out the mesoscopic morphology, microscopic geometry, symmetry, and ultrafast dynamics of an LSA formed by β-cyclodextrin (β-CD) and sodium dodecyl sulfate (SDS). The spatial imaging also revealed correlations between these different physical properties. Such knowledge shed light on the functions and mechanical properties of LSAs. In this Feature Article, we briefly introduce the fundamental principles of the VSFG microscope and then discuss the in-depth molecular physics of the LSAs revealed by this imaging technique. The application of the VSFG microscope to the artificial LSAs also paved the way for an alternative approach to studying the structure-dynamic-function relationships of protein assemblies, which were essential for life and difficult to study because of their various and complicated interactions. We expect that the hyperspectral VSFG microscope could be broadly applied to many noncentrosymmetric soft materials.