RESUMEN
With the specter of accelerating climate change, securing access to potable water has become a critical global challenge. Atmospheric water harvesting (AWH) through metal-organic frameworks (MOFs) emerges as one of the promising solutions. The standard numerical methods applied for rapid and efficient screening for optimal sorbents face significant limitations in the case of water adsorption (slow convergence and inability to overcome high energy barriers). To address these challenges, we employed grand canonical transition matrix Monte Carlo (GC-TMMC) methodology and proposed an efficient interpolation scheme that significantly reduces the number of required simulations while maintaining accuracy of the results. Through the example of water adsorption in three MOFs: MOF-303, MOF-LA2-1, and NU-1000, we show that the extrapolation of the free energy landscape allows for prediction of the adsorption properties over a continuous range of pressure and temperature. This innovative and versatile method provides rich thermodynamic information, enabling rapid, large-scale computational screening of sorbents for adsorption, applicable for a variety of sorbents and gases. As the presented methodology holds strong applicative potential, we provide alongside this paper a modified version of the RASPA2 code with a ghost swap move implementation and a Python library designed to minimize the user's input for analyzing data derived from the TMMC simulations.
RESUMEN
The hase behavior of confined fluids adsorbed in nanopores differs significantly from their bulk counterparts and depends on the chemical and structural properties of the confining structures. In general, phase transitions in nanoconfined fluids are reflected in stepwise adsorption isotherms with a pronounced hysteresis. Here, we show experimental evidence and an in silico interpretation of the reversible stepwise adsorption isotherm which is observed when methane is adsorbed in the rigid, crystalline metal-organic framework IRMOF-1 (MOF-5). In a very narrow range of pressures, the adsorbed fluid undergoes a structural and highly cooperative reconstruction and transition between low-density and high-density nanophases, as a result of the competition between the fluid-framework and fluid-fluid interactions. This mechanism evolves with temperature: below 110 K, a reversible stepwise isotherm is observed, which is a result of the bimodal distribution of the coexisting nanophases. This temperature may be considered as a critical temperature of methane confined to nanopores of IRMOF-1. Above 110 K, as the entropy contribution increases, the isotherm shape transforms to a common continuous S-shaped form that is characteristic to a gradual densification of the adsorbed phase as the pressure increases.