Thermodynamic properties of methane/water interface predicted by molecular dynamics simulations

Interfacial Tensions, Solubilities, and Transport Properties of the H2/H2O/NaCl System: A Molecular Simulation Study

Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H2), such as underground H2 storage and H2O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H2O, (b) self-diffusivities of infinitely diluted H2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O, and (c) solubilities of H2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O. The force fields used are the TIP4P/2005 for H2O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H2. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.

Hierarchical triphase diffusion photoelectrodes for photoelectrochemical gas/liquid flow conversion

Photoelectrochemical device is a versatile platform for achieving various chemical transformations with solar energy. However, a grand challenge, originating from mass and electron transfer of triphase-reagents/products in gas phase, water/electrolyte/products in liquid phase and catalyst/photoelectrode in solid phase, largely limits its practical application. Here, we report the simulation-guided development of hierarchical triphase diffusion photoelectrodes, to improve mass transfer and ensure electron transfer for photoelectrochemical gas/liquid flow conversion. Semiconductor nanocrystals are controllably integrated within electrospun nanofiber-derived mat, overcoming inherent brittleness of semiconductors. The mechanically strong skeleton of free-standing mat, together with satisfactory photon absorption, electrical conductivity and hierarchical pores, enables the design of triphase diffusion photoelectrodes. Such a design allows photoelectrochemical gas/liquid conversion to be performed continuously in a flow cell. As a proof of concept, 16.6- and 4.0-fold enhancements are achieved for the production rate and product selectivity of methane conversion, respectively, with remarkable durability.

Resonance-Induced Reduction of Interfacial Tension of Water-Methane and Improvement of Methane Solubility in Water

Molecular dynamics simulations were performed to study the effect of the periodic oscillating electric field on the interface between water and methane. We propose a new strategy that utilizes oscillating electric fields to reduce the interfacial tension (IFT) between water and methane and increase the solubility of methane in water simultaneously. These are attributed to the hydrogen bond resonance induced by an electric field with a frequency close to the natural frequency of the hydrogen bond. The resonance breaks the hydrogen bond network among water molecules to the maximum, which destroys the hydration shell and reduces the cohesive action of water, thus resulting in the decrease of IFT and the increase of methane solubility. As the frequency of the electric field is close to the optimum resonant frequency of hydrogen bonds, IFT decreases from 56.43 to 5.66 mN/m; water and methane are miscible because the solubility parameter of water reduces from 47.63 to 2.85 MPa^1/2, which is close to that of methane (3.43 MPa^1/2). Our results provide a new idea for reducing the water-gas IFT and improving the solubility of insoluble gas in water and theoretical guidance in the fields of natural gas exploitation, hydrate generation, and nanobubble nucleation.

A novel hybrid method for the calculation of methane hydrate-water interfacial tension along the three-phase (hydrate-liquid water-vapor) equilibrium line

We use a novel hybrid method to explore the temperature dependence of the solid-liquid interfacial tension of a system that consists of solid methane hydrate and liquid water. The calculated values along the three-phase (hydrate-liquid water-vapor) equilibrium line are obtained through the combination of available experimental measurements and computational results that are based on approaches at the atomistic scale, including molecular dynamics and Monte Carlo. An extensive comparison with available experimental and computational studies is performed, and a critical assessment and re-evaluation of previously reported data is presented.

Comparative Analysis of Four Neural Network Models on the Estimation of CO2-Brine Interfacial Tension

During the CO₂ injection of geological carbon sequestration and CO₂-enhanced oil recovery, the contact of CO₂ with underground salt water is inevitable, where the interfacial tension (IFT) between gas and liquid determines whether the projects can proceed smoothly. In this paper, three traditional neural network models, the wavelet neural network (WNN) model, the back propagation (BP) model, and the radical basis function model, were applied to predict the IFT between CO₂ and brine with temperature, pressure, monovalent cation molality, divalent cation molality, and molar fraction of methane and nitrogen impurities. A total of 974 sets of experimental data were divided into two data groups, the training group and the testing group. By optimizing the WNN model (I_WNN), a most stable and precise model is established, and it is found that temperature and pressure are the main parameters affecting the IFT. Through the comparison of models, it is found that I_WNN and BP models are more suitable for the IFT evaluation between CO₂ and brine.

Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations

Accurate prediction of gas solubility in a liquid is crucial in many areas of chemistry, and a detailed understanding of the molecular mechanism of the gas solvation continues to be an active area of research. Here, we extend the idea of the constant chemical potential molecular dynamics (CμMD) approach to the calculation of the gas solubility in the liquid under constant gas chemical potential conditions. As a representative example, we utilize this method to calculate the isothermal solubility of carbon dioxide in water. Additionally, we provide microscopic insight into the mechanism of solvation that preferentially occurs in areas of the surface where the hydrogen network is broken.

Molecular dynamics simulations of interfacial properties of the CO2-water and CO2-CH4-water systems

Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO₂ and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO₂-water interface has a lower interfacial tension than the CH₄-water interface. A greater number of CO₂ molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO₂ and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO₂. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO₂-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO₂ capture and sequestration processes are discussed.

Squeezing Oil into Water under Pressure: Inverting the Hydrophobic Effect

The molecular structure of dense homogeneous fluid water-methane mixtures has been determined for the first time using high-pressure neutron-scattering techniques at 1.7 and 2.2 GPa. A mixed state with a fully H-bonded water network is revealed. The hydration shell of the methane molecules is, however, revealed to be pressure-dependent with an increase in the water coordination between 1.7 and 2.2 GPa. In parallel, ab initio molecular dynamics simulations have been performed to provide insight into the microscopic mechanisms associated with the phenomenon of mixing. These calculations reproduce the observed phase change from phase separation to mixing with increasing pressure. The calculations also reproduce the experimentally observed structural properties. Unexpectedly, the simulations show mixing is accompanied by a subtle enhancement of the polarization of methane. Our results highlight the key role played by fine electronic effects on miscibility and the need to readjust our fundamental understanding of hydrophobicity to account for these.

Confined methane-water interfacial layers and thickness measurements using in situ Raman spectroscopy

Gas-liquid interfaces broadly impact our planet, yet confined interfaces behave differently than unconfined ones. We report the role of tangential fluid motion in confined methane-water interfaces. The interfaces are created using microfluidics and investigated by in situ 1D, 2D and 3D Raman spectroscopy. The apparent CH₄ and H₂O concentrations are reported for Reynolds numbers (Re), ranging from 0.17 to 8.55. Remarkably, the interfaces are comprised of distinct layers of thicknesses varying from 23 to 57 μm. We found that rarefaction, mixture, thin film, and shockwave layers together form the interfaces. The results indicate that the mixture layer thickness (δ) increases with Re (δ ∝ Re), and traditional transport theory for unconfined interfaces does not explain the confined interfaces. A comparison of our results with thin film theory of air-water interfaces (from mass transfer experiments in capillary microfluidics) supports that the hydrophobicity of CH₄ could decrease the strength of water-water interactions, resulting in larger interfacial thicknesses. Our findings help explain molecular transport in confined gas-liquid interfaces, which are common in a broad range of societal applications.

Factors governing the behaviour of aqueous methane in narrow pores

All-atom equilibrium molecular dynamics simulations were employed to investigate the behaviour of aqueous methane confined in 1-nm-wide pores obtained from different materials. Models for silica, alumina and magnesium oxide were used to construct the slit-shaped pores. The results show that methane solubility in confined water strongly depends on the confining material, with silica yielding the highest solubility in the systems considered here. The molecular structure of confined water differs within the three pores, and density fluctuations reveal that the silica pore is effectively less 'hydrophilic' than the other two pores considered. Comparing the water fluctuation autocorrelation function with local diffusion coefficients of methane across the hydrated pores we observed a direct proportional coupling between methane and water dynamics. These simulation results help to understand the behaviour of gas in water confined within narrow subsurface formations, with possible implications for fluid transport.

Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations

The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.

Vlugt TJ. Compressibility, thermal expansion coefficient and heat capacity of CH4 and CO2 hydrate mixtures using molecular dynamics simulations

Understanding the thermal and mechanical properties of CH₄ and CO₂ hydrates is essential for the replacement of CH₄ with CO₂ in natural hydrate deposits as well as for CO₂ sequestration and storage.

Two-component order parameter for quantifying clathrate hydrate nucleation and growth

Methane clathrate hydrate nucleation and growth is investigated via analysis of molecular dynamics simulations using a new order parameter. This order parameter (OP), named the Mutually Coordinated Guest (MCG) OP, quantifies the appearance and connectivity of molecular clusters composed of guests separated by water clusters. It is the first two-component OP used for quantifying hydrate nucleation and growth. The algorithm for calculating the MCG OP is described in detail. Its physical motivation and advantages compared to existing methods are discussed.

Insight into the molecular mechanism of water evaporation via the finite temperature string method

The process of water's evaporation at its liquid/air interface has proven challenging to study experimentally and, because it constitutes a rare event on molecular time scales, presents a challenge for computer simulations as well. In this work, we simulated water's evaporation using the classical extended simple point charge model water model, and identified a minimum free energy path for this process in terms of 10 descriptive order parameters. The measured free energy change was 7.4 kcal/mol at 298 K, in reasonable agreement with the experimental value of 6.3 kcal/mol, and the mean first-passage time was 1375 ns for a single molecule, corresponding to an evaporation coefficient of 0.25. In the observed minimum free energy process, the water molecule diffuses to the surface, and tends to rotate so that its dipole and one O-H bond are oriented outward as it crosses the Gibbs dividing surface. As the water molecule moves further outward through the interfacial region, its local density is higher than the time-averaged density, indicating a local solvation shell that protrudes from the interface. The water molecule loses donor and acceptor hydrogen bonds, and then, with its dipole nearly normal to the interface, stops donating its remaining hydrogen bond. At that point, when the final, accepted hydrogen bond is broken, the water molecule is free. We also analyzed which order parameters are most important in the process and in reactive trajectories, and found that the relative orientation of water molecules near the evaporating molecule, and the number of accepted hydrogen bonds, were important variables in reactive trajectories and in kinetic descriptions of the process.