Simulation of the THF hydrate–water interfacial free energy from computer simulation

Abstract

In this work, the tetrahydrofuran (THF) hydrate–water interfacial free energy is determined at 500 bar, at one point of the univariant two-phase coexistence line of the THF hydrate, by molecular dynamics simulation. The mold integration–host methodology, an extension of the original mold integration technique to deal with hydrate–fluid interfaces, is used to calculate the interfacial energy. Water is described using the well-known TIP4P/Ice model, and THF is described using a rigid version of the TraPPE model. We have recently used the combination of these two models to accurately describe the univariant two-phase dissociation line of the THF hydrate in a wide range of pressures from computer simulation [Algaba et al., J. Chem. Phys. 160, 164718 (2024)]. The THF hydrate–water interfacial free energy predicted in this work is compared with the only experimental data available in the literature. The value obtained, 27(2) mJ/m2, is in excellent agreement with the experimental data taken from the literature, 24(8) mJ/m2. To the best of our knowledge, this is the first time that the THF hydrate–water interfacial free energy is predicted from computer simulation. This work confirms that the mold integration technique can be used with confidence to predict the solid–fluid interfaces of complex structures, including hydrates that exhibit sI and sII crystallographic structures.