Lattice constants of pure methane and carbon dioxide hydrates at low temperatures. Implementing quantum corrections to classical molecular dynamics studies

Dynamic Dissociation Behaviors of sII Hydrates in Liquid Water by Heating: A Molecular Dynamics Simulation Approach

An understanding of the dynamic behavior of subtle hydrate dissociation in the liquid water phase is fundamental for gas production from marine hydrate reservoirs. Molecular dynamics simulations are performed in this study to investigate the dissociation kinetics of pure propane and binary propane + methane sII hydrates in a liquid water environment. The results show that faster hydrate dissociation rates are observed at higher initial temperatures. The hydrate phase dissociates from the cluster surface to the inside in a layer-by-layer manner under the simulation temperature conditions, which is similar to the behavior of sI hydrates and is independent of the hydrate crystal type. Compared to the binary sII hydrate, the pure sII hydrate dissociates more easily under the same initial temperature conditions, which can be attributed to the stabilizing effect of guest molecules in the hydrate cages. The empty cages collapse in one step, in contrast to the two-step pathway induced by the guest-host interaction. In addition, a hydrocarbon phase forms in the binary hydrate dissociation system instead of nanobubbles. These results can provide molecular-level insights into the dynamic mechanism of hydrate dissociation and theoretical guidance for gas recovery by thermal injection from marine hydrate reservoirs.

A novel hybrid method for the calculation of methane hydrate-water interfacial tension along the three-phase (hydrate-liquid water-vapor) equilibrium line

We use a novel hybrid method to explore the temperature dependence of the solid-liquid interfacial tension of a system that consists of solid methane hydrate and liquid water. The calculated values along the three-phase (hydrate-liquid water-vapor) equilibrium line are obtained through the combination of available experimental measurements and computational results that are based on approaches at the atomistic scale, including molecular dynamics and Monte Carlo. An extensive comparison with available experimental and computational studies is performed, and a critical assessment and re-evaluation of previously reported data is presented.

Local structure and distortions of mixed methane-carbon dioxide hydrates

A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.

Structure and thermodynamics of empty clathrate hydrates below the freezing point of water

Recently prepared as a new H2O phase, ice XVI was obtained by degassing a Ne-sII clathrate hydrate under vacuum, however very little is known of that crystalline solid under temperatures (T ≤ 220 K) and pressures (p ≤ 5000 bar) relevant for the Earth's environment and geochemistry. In this work, atomically detailed calculations using long time-scale molecular simulations, seldom paralelled before, are employed to probe empty sII clathrate hydrates. It is found that the volumetric response to an applied pressure-temperature gradient is accurately described by the Parsafar and Mason equation of state with an accuracy of at least 99.7%. Structural deformation induced upon the crystals is interpreted by monitoring the unit cell length and isobaric thermal expansivity, whilst benchmarked against previous neutron diffraction measurements of ice XVI and hexagonal ice under room pressure conditions; a critical comparison is established with other sII guest occupied lattices (CH4, CO2 and CnH2n+2 with n = 2, 3, 4), often found in permafrost regions and in the margins of continental shelves. Such an analysis reveals that empty sII frameworks are slightly more stable to thermal deformation than their sI analogues and that hexagonal ice is the structurally most stable of the condensed H2O phases addressed here. Of paramount importance for the oil and natural gas industries, heat capacities obtained from enthalpy profiles are identical for the sI and sII empty clathrates up to 2000 bar and diverge by only ∼7.3% at 5000 bar. The canonical tetrahedral symmetry of water-bonded networks is analysed in terms of an angular and a distance order parameters, which are observed to decrease (increase) as pressure (temperature) increases (decreases). The results now being reported constitute a landmark for future studies dealing with high-pressure and very low-temperature conditions, characteristic of the Earth's permafrost environment and other planetary interiors, whilst contributing to expand our knowledge regarding the recently discovered ice XVI phase.