The coexistence temperature of hydrogen clathrates: A molecular dynamics study.

Extensive molecular dynamics simulations in the equilibrium isobaric-isothermal (NPT) ensemble were developed to determine the coexistence temperatures of the water hydrogen mixture using the direct coexistence method. The water molecules were modeled using the four-site TIP4P/Ice analytical potential, and the hydrogen molecules were described using a three-site potential. The simulations were performed at different pressures (p = 900, 1500, 3000, and 4000 bars). At each pressure, a series of simulations were developed at different temperatures (from 230 to 270 K). Our results followed a line parallel to the experimental coexistence temperatures and underestimated these temperatures by approximately 25 K in the investigated range. The final configurations could or could not contain a fluid phase depending on the pressure, in accordance with the phase diagram. In addition, we explored the dynamics of the H2 molecules through clathrate hydrate cages and observed different behaviors of the H2 molecules in the small cages and the large cages of the sII structure.