ABSTRACT
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.
ABSTRACT
Molecular dynamics simulations in the equilibrium isobaric-isothermal (NPT) ensemble were used to examine the strength of an external electric field required to dissociate the methane hydrate sI structure. The water molecules were modeled using the four-site TIP4P/Ice analytical potential and methane was described as a simple Lennard-Jones interaction site. A series of simulations were performed at T = 260 K with P = 80 bars and at T = 285 K with P = 400 bars with an applied electric field ranging from 1.0 V nm(-1) to 5.0 V nm(-1). For both (T,P) conditions, applying a field greater than 1.5 V nm(-1) resulted in the orientation of the water molecules such that an ice Ih-type structure was formed, from which the methane was segregated. When the simulations were continued without the external field, the ice-like structures became disordered, resulting in two separate phases: gas methane and liquid water.
