Rotationally invariant local bond order parameters for accurate determination of hydrate structures

TL;DR

This paper introduces a new combination of local bond order parameters,  , that more accurately distinguish hydrate structures from liquid water in molecular simulations, improving upon previous methods.

Contribution

The study identifies and validates a novel   combination of bond order parameters for better hydrate-liquid phase differentiation in simulations.

Findings

The   combination outperforms previous options in hydrate detection.

Applicable to multiple hydrate types including sI and sII structures.

Enhances accuracy in hydrate structure analysis in molecular simulations.

Abstract

Averaged local bond order parameters based on spherical harmonics, also known as Lechner and Dellago order parameters, are routinely used to determine crystal structures in molecular simulations. Among different options, the combination of the $\overline{q}_{4}$ and $\overline{q}_{6}$ parameters is one of the best choices in the literature since allows one to distinguish, not only between solid- and liquid-like particles but also between different crystallographic phases, including cubic and hexagonal phases. Recently, Algaba et al. [J. Colloid Interface Sci. 623, 354, (2022)] have used the Lechner and Dellago order parameters to distinguish hydrate- and liquid-like water molecules in the context of determining the carbon dioxide hydrate-water interfacial free energy. According to the results, the preferred combination previously mentioned is not the best option to differentiate between hydrate- and liquid-like water molecules. In this work, we revisit and extend the use of these parameters to deal with systems in which clathrate hydrates phases coexist with liquid phases of water. We consider carbon dioxide, methane, tetrahydrofuran, nitrogen, and hydrogen hydrates that exhibit sI and sII crystallographic structures. We find that the $\overline{q}_{3}$ and $\overline{q}_{12}$ combination is the best option possible between a large number of possible different pairs to distinguish between hydrate- and liquid-like water molecules in all cases.