Atomistic modeling of structure II gas hydrate mechanics: Compressibility and equations of state

This work uses density functional theory (DFT) to investigate the poorly characterized structure II gas hydrates, for various guests (empty, propane, butane, ethane-methane, propane-methane), at the atomistic scale to determine key structure and mechanical properties such as equilibrium lattice volume and bulk modulus. Several equations of state (EOS) for solids (Murnaghan, Birch-Murnaghan, Vinet, Liu) were fitted to energy-volume curves resulting from structure optimization simulations. These EOS, which can be used to characterize the compressional behaviour of gas hydrates, were evaluated in terms of their robustness. The three-parameter Vinet EOS was found to perform just as well if not better than the four-parameter Liu EOS, over the pressure range in this study. As expected, the Murnaghan EOS proved to be the least robust. Furthermore, the equilibrium lattice volumes were found to increase with guest size, with double-guest hydrates showing a larger increase than single-guest hydrates, which has significant implications for the widely used van der Waals and Platteeuw thermodynamic model for gas hydrates. Also, hydrogen bonds prove to be the most likely factor contributing to the resistance of gas hydrates to compression; bulk modulus was found to increase linearly with hydrogen bond density, resulting in a relationship that could be used predictively to determine the bulk modulus of various structure II gas hydrates. Taken together, these results fill a long existing gap in the material chemical physics of these important clathrates.