Mechanisms for thermal conduction in methane hydrate

Geologic Constraints on the Formation and Evolution of Saturn's Mid-Sized Moons

Saturn's mid-sized icy moons have complex relationships with Saturn's interior, the rings, and with each other, which can be expressed in their shapes, interiors, and geology. Observations of their physical states can, thus, provide important constraints on the ages and formation mechanism(s) of the moons, which in turn informs our understanding of the formation and evolution of Saturn and its rings. Here, we describe the cratering records of the mid-sized moons and the value and limitations of their use for constraining the histories of the moons. We also discuss observational constraints on the interior structures of the moons and geologically-derived inferences on their thermal budgets through time. Overall, the geologic records of the moons (with the exception of Mimas) include evidence of epochs of high heat flows, short- and long-lived subsurface oceans, extensional tectonics, and considerable cratering. Curiously, Mimas presents no clear evidence of an ocean within its surface geology, but its rotation and orbit indicate a present-day ocean. While the moons need not be primordial to produce the observed levels of interior evolution and geologic activity, there is likely a minimum age associated with their development that has yet to be determined. Uncertainties in the populations impacting the moons makes it challenging to further constrain their formation timeframes using craters, whereas the characteristics of their cores and other geologic inferences of their thermal evolutions may help narrow down their potential histories. Disruptive collisions may have also played an important role in the formation and evolution of Saturn's mid-sized moons, and even the rings of Saturn, although more sophisticated modeling is needed to determine the collision conditions that produce rings and moons that fit the observational constraints. Overall, the existence and physical characteristics of Saturn's mid-sized moons provide critical benchmarks for the development of formation theories.

Thermal Conductivity of a Fluid-Filled Nanoporous Material: Underlying Molecular Mechanisms and the Rattle Effect

Nanoporous materials are central to the energy and environmental crisis, with key applications in adsorption, separation, and catalysis. While confinement and surface effects on fluids severely confined in their porosity are well documented, the thermal behavior of nanoporous solids subjected to fluid adsorption remains puzzling in many aspects. With striking phenomena such as the so-called rattle effect, through which fluid/solid collisions decrease the overall thermal conductivity, the solid thermal conductivity and, more generally, heat transfer and dispersion in these complex systems challenge classical approaches (e.g., mixing rules including effective medium approaches fail to capture such effects as shown here). In particular, a robust molecular framework to describe the crossover between the decrease in thermal conductivity through the rattle effect in very narrow pores and the increase in thermal conductivity when replacing vacuum with a fluid phase in larger pores is still missing. Here, using a prototypical model of fluid-filled nanoporous materials (a Lennard-Jones phase confined in an all-silica zeolite), we perform a molecular simulation study to shed light on the parameters that govern the rattle effect in nanoporous solids. First, by varying the fluid/fluid, fluid/solid, and solid/solid interaction strengths as well as the fluid number density and mass density, we unravel the ingredients that lead to the essential coupling between fluid adsorption and phonon transport. Second, despite this complex interplay, inspired by pioneering molecular approaches on the rattle effect, we show that all data obey a simple statistical physics model that relies on the change in the speed of sound due to the fluid adsorbed density and the decrease in phonon lifetime due to scattering by fluid molecules. This framework, which provides a simple formalism to rationalize the thermal behavior of this class of solid/fluid composites, points to a decrease in thermal conductivity upon fluid confinement (up to 30% in some cases). Such an effect paves the way for the design of novel applications involving fluids in interaction with nanoporous materials.

Effect of pressure on the carbon dioxide hydrate-water interfacial free energy along its dissociation line

We investigate the effect of pressure on the carbon dioxide (CO2) hydrate-water interfacial free energy along its dissociation line using advanced computer simulation techniques. In previous works, we have determined the interfacial energy of the hydrate at 400 bars using the TIP4P/Ice and TraPPE molecular models for water and CO2, respectively, in combination with two different extensions of the Mold Integration technique [J. Colloid Interface Sci. 623, 354 (2022) and J. Chem. Phys. 157, 134709 (2022)]. Results obtained from computer simulation, 29(2) and 30(2) mJ/m2, are found to be in excellent agreement with the only two measurements that exist in the literature, 28(6) mJ/m2 determined by Uchida et al. [J. Phys. Chem. B 106, 8202 (2002)] and 30(3) mJ/m2 determined by Anderson et al. [J. Phys. Chem. B 107, 3507 (2002)]. Since the experiments do not allow to obtain the variation of the interfacial energy along the dissociation line of the hydrate, we extend our previous studies to quantify the effect of pressure on the interfacial energy at different pressures. Our results suggest that there exists a correlation between the interfacial free energy values and the pressure, i.e., it decreases with the pressure between 100 and 1000 bars. We expect that the combination of reliable molecular models and advanced simulation techniques could help to improve our knowledge of the thermodynamic parameters that control the interfacial free energy of hydrates from a molecular perspective.

Simulation of the carbon dioxide hydrate-water interfacial energy

HYPOTHESIS

Carbon dioxide hydrates are ice-like nonstoichiometric inclusion solid compounds with importance to global climate change, and gas transportation and storage. The thermodynamic and kinetic mechanisms that control carbon dioxide nucleation critically depend on hydrate-water interfacial free energy. Only two independent indirect experiments are available in the literature. Interfacial energies show large uncertainties due to the conditions at which experiments are performed. Under these circumstances, we hypothesize that accurate molecular models for water and carbon dioxide combined with computer simulation tools can offer an alternative but complementary way to estimate interfacial energies at coexistence conditions from a molecular perspective.

CALCULATIONS

We have evaluated the interfacial free energy of carbon dioxide hydrates at coexistence conditions (three-phase equilibrium or dissociation line) implementing advanced computational methodologies, including the novel Mold Integration methodology. Our calculations are based on the definition of the interfacial free energy, standard statistical thermodynamic techniques, and the use of the most reliable and used molecular models for water (TIP4P/Ice) and carbon dioxide (TraPPE) available in the literature.

FINDINGS

We find that simulations provide an interfacial energy value, at coexistence conditions, consistent with the experiments from its thermodynamic definition. Our calculations are reliable since are based on the use of two molecular models that accurately predict: (1) The ice-water interfacial free energy; and (2) the dissociation line of carbon dioxide hydrates. Computer simulation predictions provide alternative but reliable estimates of the carbon dioxide interfacial energy. Our pioneering work demonstrates that is possible to predict interfacial energies of hydrates from a truly computational molecular perspective and opens a new door to the determination of free energies of hydrates.

Role of mechanical deformation in the thermal transport of sI-type methane hydrate

Natural gas hydrates (NGHs) are rising as an unconventional energy resource. The fundamental thermal characteristics of NGHs are of importance for natural gas exploitation from permafrost and oceanic sediments that are geomechanically deformed. Here, utilizing classic molecular dynamics simulations with all-atom (AA) and coarse-grained (CG) models of the methane guest molecule, the effects of mechanical strain on the thermal conductivity of sI-type methane hydrate are for the first time examined. Upon triaxial tension and compression, methane hydrate exhibits strong asymmetry in the stress responses. As the triaxial loads go from compression to tension, a reduction trend in the thermal conductivity is revealed for methane hydrate with both AA and CG models of methane, within a maximum reduction of over 44%. This reduction is because triaxial strain from compression to tension softens the phonon modes. Interestingly, there is a sudden rise in thermal conductivity at critical triaxial strain of 0.06, originating from that, at which, the phonon modes are hardened and the peaks of radial distribution functions are shifted back. This study provides important information on the thermal conductivity of methane hydrate, which is helpful for the practical production of natural gas from geo-deformed NGH-bearing sediments via a heating technique as well as evaluating their stability.

Reduced phase stability and faster formation/dissociation kinetics in confined methane hydrate

The mechanisms involved in the formation/dissociation of methane hydrate confined at the nanometer scale are unraveled using advanced molecular modeling techniques combined with a mesoscale thermodynamic approach. Using atom-scale simulations probing coexistence upon confinement and free energy calculations, phase stability of confined methane hydrate is shown to be restricted to a narrower temperature and pressure domain than its bulk counterpart. The melting point depression at a given pressure, which is consistent with available experimental data, is shown to be quantitatively described using the Gibbs-Thomson formalism if used with accurate estimates for the pore/liquid and pore/hydrate interfacial tensions. The metastability barrier upon hydrate formation and dissociation is found to decrease upon confinement, therefore providing a molecular-scale picture for the faster kinetics observed in experiments on confined gas hydrates. By considering different formation mechanisms-bulk homogeneous nucleation, external surface nucleation, and confined nucleation within the porosity-we identify a cross-over in the nucleation process; the critical nucleus formed in the pore corresponds either to a hemispherical cap or to a bridge nucleus depending on temperature, contact angle, and pore size. Using the classical nucleation theory, for both mechanisms, the typical induction time is shown to scale with the pore volume to surface ratio and hence the pore size. These findings for the critical nucleus and nucleation rate associated with such complex transitions provide a means to rationalize and predict methane hydrate formation in any porous media from simple thermodynamic data.

Gas hydrates in sustainable chemistry

This review includes the current state of the art understanding and advances in technical developments about various fields of gas hydrates, which are combined with expert perspectives and analyses.

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.

Energy and Environmental Analysis of Membrane-Based CH4-CO2 Replacement Processes in Natural Gas Hydrates

Natural gas hydrates are the largest reservoir of natural gas worldwide. This paper proposes and analyzes the CH4-CO2 replacement in the hydrate phase and pure methane collection through the use of membrane-based separation. The investigation uses a 1 L lab reactor, in which the CH4 hydrates are formed in a quartz sand matrix partially saturated with water. CH4 is subsequently dissociated with a CO2 stream supplied within the sediment inside the reactor. An energy and environmental analysis was carried out to prove the sustainability of the process. Results show that the process energy consumption constitutes 4.75% of the energy stored in the recovered methane. The carbon footprint of the CH4-CO2 exchange process is calculated as a balance of the CO2 produced in the process and the CO2 stored in system. Results provide an estimated negative value, equal to 0.004 moles sequestrated, thus proving the environmental benefit of the exchange process.

Non-equilibrium molecular-dynamics study of electromagnetic-field-induced propane-hydrate dissociation

Non-equilibrium molecular-dynamics simulations have been performed for dissolution of planar propane-hydrate/water interfaces in externally-applied electromagnetic (e/m) fields in the microwave to far infrared range (∼2.45-200 GHz) at electric-field intensities up to 2.0 V/nm and at roughly 20 K over/under temperatures vis-à-vis the zero-field propane-hydrate melting point. Upon e/m-field application, there is a field-frequency threshold above which the dissociation rate drops significantly, with a plateau therein for larger-frequencies. It was found that higher intensity and lower frequency facilitates dissociation. Except in the presence of a thermal driving-force, the 10 GHz frequency shows more substantial rate-enhancement effect vis-à-vis static electric fields or, indeed, lower-frequency e/m fields.

Molecular-dynamics study of propane-hydrate dissociation: Fluctuation-dissipation and non-equilibrium analysis

Equilibrium and non-equilibrium molecular-dynamics (MD) simulations have been performed to investigate thermal-driven break-up of planar propane-hydrate interfaces in contact with liquid water over the 260-320 K range. Two types of hydrate-surface water-lattice molecular termination were adopted, at the hydrate edge with water, for comparison: a 001-direct surface cleavage and one with completed cages. Statistically significant differences in melting temperatures and initial break-up rates were observed between both interface types. Dissociation rates were observed to be strongly dependent on temperature, with higher rates at larger over-temperatures vis-à-vis melting. A simple coupled mass and heat transfer model, developed previously, was applied to fit the observed dissociation profiles, and this helps us to identify clearly two distinct hydrate-decomposition régimes; following a highly temperature-dependent break-up phase, a second well-defined stage is essentially independent of temperature, in which the remaining nanoscale, de facto two-dimensional system's lattice framework is intrinsically unstable. Further equilibrium MD-analysis of the two-phase systems at their melting point, with consideration of the relaxation times gleaned from the auto-correlation functions of fluctuations in a number of enclathrated guest molecules, led to statistically significant differences between the two surface-termination cases; a consistent correlation emerged in both cases between the underlying, non-equilibrium, thermal-driven dissociation rates sampled directly from melting with that from an equilibrium-MD fluctuation-dissipation approach.

Direct measurement of individual phonon lifetimes in the clathrate compound Ba7.81Ge40.67Au5.33

Engineering lattice thermal conductivity requires to control the heat carried by atomic vibration waves, the phonons. The key parameter for quantifying it is the phonon lifetime, limiting the travelling distance, whose determination is however at the limits of instrumental capabilities. Here, we show the achievement of a direct quantitative measurement of phonon lifetimes in a single crystal of the clathrate Ba_7.81Ge_40.67Au_5.33, renowned for its puzzling ‘glass-like’ thermal conductivity. Surprisingly, thermal transport is dominated by acoustic phonons with long lifetimes, travelling over distances of 10 to 100 nm as their wave-vector goes from 0.3 to 0.1 Å⁻¹. Considering only low-energy acoustic phonons, and their observed lifetime, leads to a calculated thermal conductivity very close to the experimental one. Our results challenge the current picture of thermal transport in clathrates, underlining the inability of state-of-the-art simulations to reproduce the experimental data, thus representing a crucial experimental input for theoretical developments.

Phonon lifetime is a fundamental parameter of thermal transport however its determination is challenging. Using inelastic neutron scattering and the neutron resonant spin-echo technique, Lory et al. determine the acoustic phonon lifetime in a single crystal of clathrate Ba7.81Ge40.67Au5.33.

Elasticity and Stability of Clathrate Hydrate: Role of Guest Molecule Motions

Molecular dynamic simulations were performed to determine the elastic constants of carbon dioxide (CO₂) and methane (CH₄) hydrates at one hundred pressure–temperature data points, respectively. The conditions represent marine sediments and permafrost zones where gas hydrates occur. The shear modulus and Young’s modulus of the CO₂ hydrate increase anomalously with increasing temperature, whereas those of the CH₄ hydrate decrease regularly with increase in temperature. We ascribe this anomaly to the kinetic behavior of the linear CO₂ molecule, especially those in the small cages. The cavity space of the cage limits free rotational motion of the CO₂ molecule at low temperature. With increase in temperature, the CO₂ molecule can rotate easily, and enhance the stability and rigidity of the CO₂ hydrate. Our work provides a key database for the elastic properties of gas hydrates, and molecular insights into stability changes of CO₂ hydrate from high temperature of ~5 °C to low decomposition temperature of ~−150 °C.

Microscopic Origin of Strain Hardening in Methane Hydrate

It has been reported for a long time that methane hydrate presents strain hardening, whereas the strength of normal ice weakens with increasing strain after an ultimate strength. However, the microscopic origin of these differences is not known. Here, we investigated the mechanical characteristics of methane hydrate and normal ice by compressive deformation test using molecular dynamics simulations. It is shown that methane hydrate exhibits strain hardening only if the hydrate is confined to a certain finite cross-sectional area that is normal to the compression direction. For normal ice, it does not present strain hardening under the same conditions. We show that hydrate guest methane molecules exhibit no long-distance diffusion when confined to a finite-size area. They appear to serve as non-deformable units that prevent hydrate structure failure, and thus are responsible for the strain-hardening phenomenon.

Massively parallel molecular dynamics simulation of formation of ice-crystallite precursors in supercooled water: incipient-nucleation behavior and role of system size

Ice-crystallite precursor formation in supercooled water was studied via molecular dynamics for systems ranging from ∼10^{6} to 8.6×10^{6} molecules, using a tetrahedrally biased single-site "mW" model. This has established system-size effects in the early onset of nucleation, so as to study often-transient precursors' beguiling propensity to "flicker" into instantaneous locally ordered molecular arrangements redolent of ice. In addition, the adoption of solidlike and liquidlike bimodal local configurational-energy distributions was observed, characteristic of early nucleation. Larger systems favored a higher probability of precursor formation, although such ones were not usually longer lived relative to those in smaller systems (which themselves are rather transient). It was concluded tentatively that subtle effects of differences in systemwide density fluctuations and accessible lower-frequency modes tend to favor precursor formation in larger systems, although not necessarily the precursor's kinetic stability.

Thermal Conductivity of Supercooled Water: An Equilibrium Molecular Dynamics Exploration

The thermal conductivity of both supercooled and ambient-temperature water at atmospheric pressure has been computed over the 140-270 K temperature range for three popular water models via equilibrium molecular dynamics in the Green-Kubo setting. No strong temperature dependence of thermal conductivity was observed. The underlying phonon modes contributing to thermal conduction processes have been examined in the present work, and it has been established that (translational) acoustic modes dominate in supercooled water.

Molecular Dynamics Study of Thermal Conduction in Carbon Dioxide Hydrates

Equilibrium molecular dynamics simulations that use the Green–Kubo method for sI CO2-hydrate systems from medium to full occupancy were performed to estimate the corresponding thermal conductivities at temperatures that range from 233.15K to 278.15K and pressures that range from 3MPa to 100MPa. Specific potential models for water and CO2were adopted. The effects of guest occupancy ratios and outside thermobaric conditions on CO2hydrate thermal conductivity were studied. The thermal mechanism was also analyzed. The thermal conductivities of hydrates of CH4, C2H6, N2, and O2were estimated. The size ratio of guest diameter to cavity diameter provided an adequate basis for understanding the thermal conductivities of gas hydrates.

Stability, adsorption, and diffusion of CH₄, CO₂, and H₂ in clathrate hydrates

We present a study of the adsorption and diffusion of CH₄, CO₂, and H₂ molecules in clathrate hydrates using ab initio van der Waals density functional formalism [M. Dion, Phys. Rev. Lett. 92, 246401 (2004)10.1103/PhysRevLett.92.246401]. We find that the adsorption energy is dominated by van der Waals interactions and that, without them, gas hydrates would not be stable. We calculate the maximum adsorption capacity as well as the maximum hydrocarbon size that can be adsorbed. The relaxation of the host lattice is essential for a good description of the diffusion activation energies, which are estimated to be of the order of 0.2, 0.4, and 1.0 eV for H₂, CO₂, and CH₄, respectively.