Calculation of Liquid Water-Hydrate-Methane Vapor Phase Equilibria from Molecular Simulations

Lars Jensen, Kaj Thomsen, Nicolas von Solms, S. Wierzchowski, M.R. Walsh, C.A. Koh, E.D. Sloan, D.T. Wu, A.K. Sum

Research output: Contribution to journalJournal articleResearchpeer-review


Monte Carlo simulation methods for determining fluid- and crystal-phase chemical potentials are used for the first time to calculate liquid water-methane hydrate-methane vapor phase equilibria from knowledge of atomistic interaction potentials alone. The water and methane molecules are modeled using the TIP4P/ice potential and a united-atom Lennard-Jones potential. respectively. The equilibrium calculation method for this system has three components, (i) thermodynamic integration from a supercritical ideal gas to obtain the fluid-phase chemical potentials. (ii) calculation of the chemical potential of the zero-occupancy hydrate system using thermodynamic integration from an Einstein crystal reference state, and (iii) thermodynamic integration to obtain the water and guest molecules' chemical potentials as a function of the hydrate occupancy. The three-phase equilibrium curve is calculated for pressures ranging from 20 to 500 bar and is shown to follow the Clapeyron behavior, in agreement with experiment; coexistence temperatures differ from the latter by 4-16 K in the pressure range studied. The enthalpy of dissociation extracted from the calculated P-T curve is within 2% of the experimental value at corresponding conditions. While computationally intensive, simulations such as these are essential to map the thermodynamically stable conditions for hydrate systems.
Original languageEnglish
JournalJournal of Physical Chemistry Part B: Condensed Matter, Materials, Surfaces, Interfaces & Biophysical
Issue number17
Pages (from-to)5775-5782
Publication statusPublished - 2010


Dive into the research topics of 'Calculation of Liquid Water-Hydrate-Methane Vapor Phase Equilibria from Molecular Simulations'. Together they form a unique fingerprint.

Cite this