Simulation of the carbon dioxide hydrate-water interfacial energy

Hypothesis: Carbon dioxide hydrates are ice-like nonstoichiometric inclusion solid compounds with importance to global climate change, and gas transportation and storage. The thermodynamic and kinetic mechanisms that control carbon dioxide nucleation critically depend on hydrate-water interfacial free energy. Only two independent indirect experiments are available in the literature. Interfacial energies show large uncertainties due to the conditions at which experiments are performed. Under these circumstances, we hypothesize that accurate molecular models for water and carbon dioxide combined with computer simulation tools can offer an alternative but complementary way to estimate interfacial energies at coexistence conditions from a molecular perspective. Calculations: We have evaluated the interfacial free energy of carbon dioxide hydrates at coexistence conditions (three-phase equilibrium or dissociation line) implementing advanced computational methodologies, including the novel Mold Integration methodology. Our calculations are based on the definition of the interfacial free energy, standard statistical thermodynamic techniques, and the use of the most reliable and used molecular models for water (TIP4P/Ice) and carbon dioxide (TraPPE) available in the literature. Findings: We find that simulations provide an interfacial energy value, at coexistence conditions, consistent with the experiments from its thermodynamic definition. Our calculations are reliable since are based on the use of two molecular models that accurately predict: (1) The ice-water interfacial free energy; and (2) the dissociation line of carbon dioxide hydrates. Computer simulation predictions provide alternative but reliable estimates of the carbon dioxide interfacial energy. Our pioneering work demonstrates that is possible to predict interfacial energies of hydrates from a truly computational molecular perspective and opens a new door to the determination of free energies of hydrates.We thank Pedro J. Pérez for the critical reading of the manuscript. We also acknowledge Centro de Supercomputación de Galicia (CESGA, Santiago de Compostela, Spain) and MCIA (Mésocentre de Calcul Intensif Aquitain) of the Universités de Bordeaux and Pau et Pays de l’Adour (France) for providing access to computing facilities. We thank financial support from the Ministerio de Economía, Industria y Competitividad (FIS2017- 89361-C3-1-P), Junta de Andalucía (P20-00363), and Universidad de Huelva (P.O. FEDER UHU-1255522), all three cofinanced by EU FEDER funds. J.A. acknowledges Contrato Predoctoral de Investigación from XIX Plan Propio de Investigación de la Universidad de Huelva and a FPU Grant (Ref. FPU15/03754) from Ministerio de Educación, Cultura y Deporte. J. A., J. M. M., and F. J. B. thankfully acknowledge the computer resources at Magerit and the technical support provided by the Spanish Supercomputing Network (RES) (Project QCM- 2018–2- 0042). Funding for open access charge: Universidad de Huelva / CBU

Solubility of Methane in Water: Some Useful Results for Hydrate Nucleation

In this paper, the solubility of methane in water along the 400 bar isobar is determined by computer simulations using the TIP4P/Ice force field for water and a simple LJ model for methane. In particular, the solubility of methane in water when in contact with the gas phase and the solubility of methane in water when in contact with the hydrate has been determined. The solubility of methane in a gas–liquid system decreases as temperature increases. The solubility of methane in a hydrate–liquid system increases with temperature. The two curves intersect at a certain temperature that determines the triple point T3 at a certain pressure. We also determined T3 by the three-phase direct coexistence method. The results of both methods agree, and we suggest 295(2) K as the value of T3 for this system. We also analyzed the impact of curvature on the solubility of methane in water. We found that the presence of curvature increases the solubility in both the gas–liquid and hydrate–liquid systems. The change in chemical potential for the formation of hydrate is evaluated along the isobar using two different thermodynamic routes, obtaining good agreement between them. It is shown that the driving force for hydrate nucleation under experimental conditions is higher than that for the formation of pure ice when compared at the same supercooling. We also show that supersaturation (i.e., concentrations above those of the planar interface) increases the driving force for nucleation dramatically. The effect of bubbles can be equivalent to that of an additional supercooling of about 20 K. Having highly supersaturated homogeneous solutions makes possible the spontaneous formation of the hydrate at temperatures as high as 285 K (i.e., 10K below T3). The crucial role of the concentration of methane for hydrate formation is clearly revealed. Nucleation of the hydrate can be either impossible or easy and fast depending on the concentration of methane which seems to play the leading role in the understanding of the kinetics of hydrate formation

Transport properties of the square-well fluid from molecular dynamics simulation

In this work, we have calculated self-diffusion and shear viscosity, two of the most important transport properties, of the spherical square-well (SW) fluid interacting with potential range λ=1.5σ. To this end, we have used a combination of molecular dynamics simulation and the continuous version of the square-well (CSW) intermolecular potential recently proposed by Zerón et al. [Mol. Phys. 116, 3355 (2018)]. In addition to that, we have also determined a number of equilibrium properties, including internal energy, compressibility factor, radial distribution function and coordination number. All properties are evaluated in a wide range of temperatures and densities, including subcritical and supercritical thermodynamic conditions. Results obtained in this work show an excellent agreement with available data reported in the literature and demonstrate that the CSW intermolecular potential can be used in molecular dynamics simulations to emulate SW transport properties with confidence.</p