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: André Guerra, Samuel Mathews, Milan Marić, Phillip Servio, Alejandro D. Rey
• Format: Article in Journal/Newspaper
• Language: English
• Published: MDPI AG 2022
• Subjects: water; methane; hydrates; viscosity; diffusivity; molecular dynamics; Organic chemistry; QD241-441; Methane hydrate
• Online Access: https://doi.org/10.3390/molecules27155019; https://doaj.org/article/57803440f74c4eb8ac3c7819692ceb01

(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.