Dissociation temperature of gas hydrates through isenthalpic-isobaric molecular dynamics simulations

Molecular simulations are a powerful tool to understand phenomena and obtain properties of gas hydrate systems. The direct coexistence method (DCM) in the NVT or NPT ensembles, the most commonly used method to determine hydrate dissociation temperatures, can be computationally expensive due to the need for several long simulations. Through an extensive set of simulations, we report here the details of the DCM within the NPH (isobaric-isenthalpic) ensemble, which require fewer and shorter trajectories. The dissociation pressure of methane hydrates is obtained for pressures of 4, 8, 15, 30, and 50 MPa. The values are in agreement with other literature simulations and experimental data. The results are further validated with the calculation of the enthalpy of dissociation, with a value of 50 kJ/mol of methane, also in agreement with the literature. The complexity of a multiphase and multicomponent system presents challenges lacking in simpler water/ice systems. These are found to be dependent on energy conservation. The optimal set of parameters to achieve it is also reported, including a smaller time step and the use of double precision, along with an analysis of some factors that could affect the convergence of the method. Although these parameters require more computational cost, the NPH ensemble is successful in providing the dissociation temperature of gas hydrates in fewer simulations than other ensembles and with productions lasting only 500 ns.