Unveiling the methane hydrate–water interfacial free energy through direct molecular simulation at coexistence conditions

Abstract

Hypothesis: The interfacial free energy between hydrate and water phases is a key thermodynamic parameter that governs both nucleation kinetics and crystal growth of gas hydrates. In these strategic materials—crystalline inclusion compounds where hydrogen-bonded water cages encapsulate small guest molecules such as methane (CH4) and carbon dioxide (CO2) -this interfacial energy plays a crucial role in determining phase stability and formation pathways. Given the significance of gas hydrates in energy storage, CO2 sequestration, and climate-related processes, accurately determining their interfacial energies is essential for advancing both fundamental understanding and technological applications. Despite its importance, the hydrate–water interfacial energy remains poorly constrained due to substantial experimental uncertainties and the limitations of indirect estimation methods. For example, reported experimental values for CH4 hydrate span a wide range from 28 to 40 mJ/m2. Interestingly, some studies suggest these values are comparable to the interfacial free energy of the hexagonal ice (ice Ih)--water interface, approximately 32 mJ/m2, hinting at potential analogies between clathrate hydrate and ice interfaces. Calculations: In this work, we present a direct molecular simulation of the CH4 hydrate–water interfacial free energy using two novel and independent extensions of the mold integration method. These extensions are specifically designed to induce the formation of a thin, planar hydrate–water interface and to compute the reversible work required to create it. For this purpose, we employ the TIP4P/Ice force field—one of the most reliable water models available—known for accurately reproducing the melting temperature of ice Ih under ambient conditions.