Simulation of the carbon dioxide hydrate-water interfacial energy

Abstract

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. 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. 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. 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.

20 pages, 7 figures

References in corpus (5)

• Accurate determination of crystal structures based on averaged local bond order parameters
• On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice
• Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water
• The mold integration method for the calculation of the crystal-fluid interfacial free energy from simulations
• Computer simulation study of surface wave dynamics at the crystal--melt interface

Cited by in corpus (8)

• Solubility of carbon dioxide in water: some useful results for hydrate nucleation
• Homogeneous nucleation rate of carbon dioxide hydrate formation under experimental condition from Seeding simulations
• Direct calculation of the planar NaCl-aqueous solution interfacial free energy at the solubility limit
• Effect of pressure on the carbon dioxide hydrate-water interfacial free energy along its dissociation line
• Simulation of the THF hydrate-water interfacial free energy from computer simulation
• Prediction of the univariant two-phase coexistence line of the tetrahydrofuran hydrate from computer simulation
• Rotationally invariant local bond order parameters for accurate determination of hydrate structures
• Prediction of the three-phase coexistence line of the ethane hydrate from molecular simulation