All-Atom Molecular Dynamics of Pure Water–Methane Gas Hydrate Systems under Pre-Nucleation Conditions: A Direct Comparison between Experiments and Simulations of Transport Properties for the Tip4p/Ice Water Model

• Published in: Molecules
• Main Authors: Guerra, André, Mathews, Samuel, Marić, Milan, Servio, Phillip, Rey, Alejandro D.
• Format: Text
• Language: English
• Published: MDPI 2022
• Subjects: Article; Methane hydrate
• Online Access: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC9370622/; https://doi.org/10.3390/molecules27155019

(1) Background: New technologies involving gas hydrates under pre-nucleation conditions such as gas separations and storage have become more prominent. This has necessitated the characterization and modeling of the transport properties of such systems. (2) Methodology: This work explored methane hydrate systems under pre-nucleation conditions. All-atom molecular dynamics simulations were used to quantify the performance of the TIP4P/2005 and TIP4P/Ice water models to predict the viscosity, diffusivity, and thermal conductivity using various formulations. (3) Results: Molecular simulation equilibrium was robustly demonstrated using various measures. The Green–Kubo estimation of viscosity outperformed other formulations when combined with TIP4P/Ice, and the same combination outperformed all TIP4P/2005 formulations. The Green–Kubo TIP4P/Ice estimation of viscosity overestimates (by 84% on average) the viscosity of methane hydrate systems under pre-nucleation conditions across all pressures considered (0–5 MPag). The presence of methane was found to increase the average number of hydrogen bonds over time (6.7–7.8%). TIP4P/Ice methane systems were also found to have 16–19% longer hydrogen bond lifetimes over pure water systems. (4) Conclusion: An inherent limitation in the current water force field for its application in the context of transport properties estimations for methane gas hydrate systems. A re-parametrization of the current force field is suggested as a starting point. Until then, this work may serve as a characterization of the deviance in viscosity prediction.