Acoustic-propagation properties of methane clathrate hydrates from non-equilibrium molecular dynamics

Given methane hydrates' importance in marine sediments, as well as the widespread use of seabed acoustic-signaling methods in oil and gas exploration, the elastic characterization of these materials is particularly relevant. A greater understanding of the properties governing phonon, sound, and acoustic propagation would help to better classify methane-hydrate deposits, aiding in their discovery. Recently, we have published a new nonequilibrium molecular-dynamics (NEMD) methodology to recreate longitudinal and transverse perturbations, observing their propagation through a crystalline lattice by various metrics, to study the underlying S- and P-wave velocities (achieving excellent agreement with experiment) [Melgar et al., J. Phys. Chem. 122(5), 3006-3013 (2018); ibid. 150, 084101 (2019)]. Here, we apply these NEMD methods to methane-clathrate systems to study acoustic-propagation characteristics, as well as the lattice elastic behavior. In so doing, we determine S- and P-wave velocities in excellent accord with experiment; we also ascertain the allowable magnitude range of acoustic perturbation and establish a threshold for lattice breakup and hydrate decomposition. Interestingly, upon dissociation, we observe the formation of methane nanobubbles, which agrees with previous studies on the microscopic fundamentals of hydrate dissociation by various means.