Despite the multiple length and time scales over which fluid-mineral interactions occur, interfacial phenomena control the exchange of matter and impact the nature of multiphase flow, as well as the reactivity of C–O–H fluids in geologic systems. In general, the properties of confined fluids, and their influence on porous geologic phenomena are much less well understood compared to those of bulk fluids. We used equilibrium molecular dynamics simulations to study fluid systems composed of propane and water, at different compositions, confined within cylindrical pores of diameter ∼16 Å carved out of amorphous silica. The simulations are conducted within a single cylindrical pore. In the simulated system all the dangling silicon and oxygen atoms were saturated with hydroxyl groups and hydrogen atoms, respectively, yielding a total surface density of 3.8 −OH/nm². Simulations were performed at 300 K, at different bulk propane pressures, and varying the composition of the system. The structure of the confined fluids was quantified in terms of the molecular distribution of the various molecules within the pore as well as their orientation. This allowed us to quantify the hydrogen bond network and to observe the segregation of propane near the pore center. Transport properties were quantified in terms of the mean square displacement in the direction parallel to the pore axis, which allows us to extract self-diffusion coefficients. The diffusivity of propane in the cylindrical pore was found to depend on pressure, as well as on the amount of water present. It was found that the propane self-diffusion coefficient decreases with increasing water loading because of the formation of water bridges across the silica pores, at sufficiently high water content, which hinder propane transport. The rotational diffusion, the lifespan of hydrogen bonds, and the residence time of water molecules at contact with the silica substrate were quantified from the simulated trajectories using the appropriate autocorrelation functions. The simulations contribute to a better understanding of the molecular phenomena relevant to the behavior of fluids in the subsurface.

Metal-Organic Frameworks (MOF) have been identified as highly efficient nanoporous adsorbents for CO₂ storage. In particular, Mg-MOF-74 has been shown to promise exceptionally high CO₂ sorption. Although several studies have reported adsorption isotherms of CO₂ in Mg-MOF-74, the effect of inter-crystalline spacing in Mg-MOF-74 on the sorption of CO₂ has not been addressed. These effects have been shown to be profound for a quadrupolar molecule like CO₂ in the case of silicalite (Phys. Chem. Chem. Phys. 22 (2020) 13951). Here, we report the effects of inter-crystalline spacing on the adsorption of CO₂ in Mg-MOF-74, studied using grand canonical Monte Carlo (GCMC) simulations. The inter-crystalline spacing is found to enhance adsorption at the crystallite surfaces. Larger inter-crystalline spacing up to twice the kinetic diameter of CO₂ results in higher adsorption and larger crystallite sizes suppress adsorption. Magnitudes of the inter-crystalline space relative to the kinetic diameter of the adsorbed fluid and the surface to volume ratio of the adsorbent crystallites are found to be important factors determining the adsorption amounts. The results of this study suggest that the ideal Mg-MOF-74 sample for CO₂ storage applications should have smaller crystallites separated from each other with an inter-crystalline space of approximately twice the kinetic diameter of CO₂.

Geo-fluids consisting of C-O-H volatiles are the main mode of transport of mass and energy throughout the lithosphere and are commonly found confined in pores, grain boundaries and fractures. The confinement of these fluids by porous media at the length scales of a few nanometers gives rise to numerous physical and chemical properties that deviate from the bulk behavior. Studying the structural and dynamical properties of these confined fluids at the length and time scales of nanometers and picoseconds respectively forms an important component of understanding their behavior. To study confined fluids, non-destructive penetrative probes are needed. Nuclear magnetic resonance (NMR) by virtue of its ability to monitor longitudinal and transverse magnetization relaxations of spins, and chemical shifts brought about by the chemical environment of a nucleus, and measuring diffusion coefficient provides a good opportunity to study dynamics and chemical structure at the molecular length and time scales. Another technique that gives insights into the dynamics and structure at these length and time scales is neutron scattering (NS). This is because the wavelength and energies of cold and thermal neutrons used in scattering experiments are in the same range as the spatial features and energies involved in the dynamical processes occurring at the molecular level. Molecular Dynamics (MD) simulations on the other hand help with the interpretation of the NMR and NS data. Simulations can also supplement the experiments by calculating quantities not easily accessible to experiments. Thus using NMR, NS and MD simulations in conjunction, a complete description of the molecular structure and dynamics of confined geo-fluids can be obtained. In the current review, our aim is to show how a synergistic use of these three techniques has helped shed light on the complex behavior of water, CO2, and low molecular weight hydrocarbons. After summarizing the theoretical backgrounds of the techniques, we will discuss some recent examples of the use of NMR, NS, and MD simulations to the study of confined fluids.