Atomistic-geometric simulations to investigate the mechanical stability of monocrystalline sI methane hydrates under pressure

Abstract

Gas hydrate mechanical stability under pressure is critically important in energy supply, global warming, and carbon-neutral technologies. The stability of these polyhedral guest–host crystals under increasing pressure is affected by host cage type and face connectivity as well as guest gas occupancy. The geometry-imposed cage connectivity generates crystal lattices that include inclusion-matrix material composite structures. In this paper, we integrate Density Functional Theory simulations with a polyhedral-inspired composite material model that quantifies stability limits, failure modes, and the impact of the type of cage occupancy. DFT reveals the existence of two failure mechanisms under increasing pressure: (i) a multistep lattice breakdown under total occupancy and under only large cage occupancy and (ii) a single-step breakdown under zero occupancy as well as with only small cage occupancy. The DFT-composite model predicts optimal occupancy pathways to generate strength and critical occupancy pathways to promote decomposition.