Prediction of the three-phase coexistence line of the ethane hydrate from molecular simulation

Abstract

We investigate the three-phase coexistence line of ethane (C${}_2$H${}_6$) hydrate through molecular dynamics simulations using the direct coexistence approach. In this framework, C${}_2$H${}_6$ sI hydrate, aqueous, and pure guest phases are constructed within a single simulation box, allowing us to monitor their mutual stability. From the temporal evolution of the potential energy, we identify the equilibrium temperature (T${}_3$) at which all three phases coexist, across pressures ranging from 1000 to 4000 bar, in accordance with available experimental data. Simulations are performed with the GROMACS package (version 2016, double precision) in the $NPT$ ensemble. Water and C${}_2$H${}_6$ molecules are represented using the TIP4P/Ice and TraPPE-UA models, respectively, while unlike non-bonded interactions are computed with the Lorentz-Berthelot combining rule. Dispersive Lennard-Jones and Coulomb interactions are truncated at 1.6 nm, with long-range Coulombic contributions treated via Particle-Mesh Ewald summation. The predicted three-phase coexistence line shows excellent agreement with experimental measurements within the investigated pressure range. These results demonstrate the suitability of the direct coexistence methodology, combined with established molecular models, for reproducing hydrate dissociation behavior in systems that have received little prior computational attention.

11 pages, 1 figures

References in corpus (16)

• Canonical sampling through velocity-rescaling
• Solubility of carbon dioxide in water: some useful results for hydrate nucleation
• Simulation of the carbon dioxide hydrate-water interfacial energy
• Three phase equilibria of the methane hydrate in NaCl solutions: A simulation study
• Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate
• Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate
• Homogeneous nucleation rate of carbon dioxide hydrate formation under experimental condition from Seeding simulations
• Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates
• Effect of pressure on the carbon dioxide hydrate-water interfacial free energy along its dissociation line
• Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation
• Simulation of the THF hydrate-water interfacial free energy from computer simulation
• Three-phase equilibria of CO${}_2$ hydrate from computer simulation in presence of NaCl
• Prediction of the univariant two-phase coexistence line of the tetrahydrofuran hydrate from computer simulation
• Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation. II. Effect of multiple occupancy
• Rotationally invariant local bond order parameters for accurate determination of hydrate structures
• Dissociation Line and Driving Force for Nucleation of the Multiple Occupied Hydrogen Hydrate from Computer Simulation