Molecular dynamics simulation of CO2 hydrates: Prediction of three phase coexistence line

J.M. Míguez (1) , M.M. Conde (1) , Jean-Philippe Torre (1) , F.J. Blas (2) , M.M. Piñeiro (3) , C. Vega (4)

1 LFCR - Laboratoire des Fluides Complexes et leurs Réservoirs
2 Universidad de Huelva
3 Universidade de Vigo
4 UCM - Universidad Complutense de Madrid = Complutense University of Madrid [Madrid]

The three phase equilibrium line (hydrate-liquid water-liquid carbon dioxide) has been estimated for the water + carbon dioxide binary mixture using molecular dynamics simulation and the direct coexistence technique. Both molecules have been represented using rigid nonpolarizable models. TIP4P/2005 and TIP4P/Ice were used for the case of water, while carbon dioxide was considered as a three center linear molecule with the parameterizations of MSM, EPM2, TraPPE, and ZD. The influence of the initial guest occupancy fraction on the hydrate stability has been analyzed first in order to determine the optimal starting configuration for the simulations, paying attention to the influence of the two different cells existing in the sI hydrate structure. The three phase coexistence temperature was then determined for a pressure range from 2 to 500 MPa. The qualitative shape of the equilibrium curve estimated is correct, including the high pressure temperature maximum that determines the hydrate re-entrant behaviour. However, in order to obtain quantitative agreement with experimental results, a positive deviation from the classical Lorentz-Berthelot combining rules must be considered.

Mots clés: CO2 -; Gas hydrates

Domaines: Chimie théorique et/ou physique

J.M. Míguez, M.M. Conde, Jean-Philippe Torre, F.J. Blas, M.M. Piñeiro, et al.. Molecular dynamics simulation of CO2 hydrates: Prediction of three phase coexistence line. Journal of Chemical Physics, 2015, 142 (12), pp.124505. ⟨10.1063/1.4916119⟩. ⟨hal-01815919⟩