Molecular dynamics study of thermal-driven methane hydrate dissociation

Comparing brute force to transition path sampling for gas hydrate nucleation with a flat interface: comments on time reversal symmetry

Fluid to solid nucleation is often investigated with the rare event method transition path sampling (TPS). I claim that the inherent irreversibility of solid nucleation, even at stationary conditions, calls into question TPS's applicability for determining solid nucleation mechanisms, especially for pre-critical behavior. Even when applied to a phenomenon which displays time reversal asymmetry like solid nucleation, TPS is a good means of exploring phase space and giving trends in post-critical structure, and its ability to facilitate nucleation rate and free energy calculations remains outstanding. Forward-only splitting and ratcheting methods such as forward flux sampling are more attractive for understanding nucleation mechanisms as they do not require time reversal symmetry, but at low driving forces may suffer from the same limitations as brute force: they may never make it to the first ratchet. Here I briefly summarize the TPS method and gas hydrate nucleation simulation literature, focusing on topics within both to facilitate a comparison of brute force hydrate nucleation to transition path sampling of hydrate nucleation. Perhaps anecdotally, the brute force technique results in more crystalline trajectories despite having higher driving forces than TPS. I maintain this difference is because of the inherent irreversibility of hydrate nucleation, meaning its pre-critical behavior cannot accurately be determined by the melting trajectories that comprise approximately half of the configurations in TPS's path ensemble.

Methane hydrate phase equilibrium considering dissolved methane concentrations and interfacial geometries from molecular simulations

Natural gas hydrates, mainly existing in permafrost and on the seabed, are expected to be a new energy source with great potential. The exploitation technology of natural gas hydrates is one of the main focuses of hydrate-related studies. In this study, a large-size liquid aqueous solution wrapping a methane hydrate system was established and molecular dynamics simulations were used to investigate the phase equilibrium conditions of methane hydrate at different methane concentrations and interfacial geometries. It is found that the methane concentration of a solution significantly affects the phase equilibrium of methane hydrates. Different methane concentrations at the same temperature and pressure can lead to hydrate formation or decomposition. At the same temperature and pressure, in a system reaching equilibrium, the size of spherical hydrate clusters is coupled to the solution concentration, which is proportional to the Laplace pressure at the solid-liquid interface. Lower solution concentrations reduce the phase equilibrium temperature of methane hydrates at the same pressure; as the concentration increases, the phase equilibrium temperature gradually approaches the actual phase equilibrium temperature. In addition, the interfacial geometry of hydrates affects the thermodynamic stability of hydrates. The spherical hydrate particles have the highest stability for the same volume. Through this study, we provide a stronger foundation to understand the principles driving hydrate formation/dissociation relevant to the exploitation of methane hydrates.

Facilitation of Hydrate Dissociation and Structural Evolution by Major Marine Anions under Static Electric Fields

Electric fields have been proven to be capable of significantly affecting the equilibrium state of hydrates. In this study, the thermodynamic properties and structural changes of methane hydrate (MH) in various anion solutions in an electric field at 0.7 V/nm were investigated by molecular dynamics simulations. The presence of anions significantly enhances the instability of methane hydrates under electric fields, leading to a staged dissociation process. First, the anions coexist with MH to form a temporary metastable structure under the action of an electric field. Then, the migration of anions causes the dissociation of nearby hydrates and the formation of flow channels in the hydrate layer, which leads to the complete dissociation of MH after a period. The promotive effects of F-, Br-, I-, and Cl- ions were close, while SO42- was relatively weak. The anions are still in hydration shells in the MH phase, but the structure of the hydration shells differs slightly from that in solution (the coordination numbers of I- and SO42- ions increased). The migration resistances of multiple anions to cross the surface of the hydrate layer are similar. However, inside the hydrate phase, the anions with a larger radius have a higher migration resistance. It is difficult for SO42- ions to migrate inside the hydrate phase, and they tend to form a metastable structure on the hydrate surface. Combining our previous studies, SrCl2 solution has the best hydrate promotion under an electric field environment.

Molecular dynamics of the spontaneous generation mechanism of natural gas hydrates during methane nanobubble rupture

Natural gas hydrates have garnered significant attention as a potential new source of alternative energy, and understanding their formation mechanism is of paramount importance for efficient utilization and pipeline transportation. However, there is no consensus among academics on the formation mechanism of natural gas hydrates. In this paper, we propose a method for promoting the rapid formation of natural gas hydrates based on the rupture of methane nanobubbles, which creates local high temperature and pressure to facilitate the mixing of methane and water. The rapid decrease in system temperature and pressure during the process further enhances the formation of gas hydrates. Using molecular dynamics simulations, we theoretically verify the formation of natural gas hydrates. Our results indicate that the instantaneous rupture of methane nanobubbles induced by shock waves leads to a dramatic increase in the local molecular motion velocity around the bubbles. This results in extreme local high temperature and high pressure, leading to complete mixing of methane and water and rapid formation of gas hydrates during the cooling and pressure drop of the mixture. We confirm our findings by analyzing F3-order parameters, F4-order parameters, and water cage statistics.

Large-Scale Nonequilibrium Molecular Studies of Thermal Hydrate Dissociation

The energy content of methane hydrate reservoirs (MHRs) is at least twice that of conventional fossil fuels. So, there is considerable interest in their commercial development by heating, among other dissociation mechanisms. However, a few researchers have highlighted the potentially uncontrollable release of methane from MHRs, which could occur because of global warming. Therefore, it is crucial to understand the kinetics of thermal hydrate dissociation to safely develop these resources and prevent the release of this greenhouse gas into the environment. Although there have been several molecular studies of thermal dissociation, most of these use small simulation domains that cannot capture the transient nature of the process. To address this limitation, we performed coarse-grained molecular dynamics (CGMD) simulations on a significantly larger domain with a hundred times more hydrate unit cells than those used in previous studies. We monitored the kinetics of dissociation using an image-processing algorithm and observed the dynamics of the process while maintaining a thermal gradient at the dissociation front. For the first time, we report the formation of an unstable secondary dissociation path that triggers gas bubbles within the solid hydrate. The kinetics of thermal dissociation appears to occur in three stages. In the first stage, the energy of the system increases until it exceeds the activation energy, and dissociation is initiated. Consistent dissociation occurs in the second stage, whereas the third stage involves the dissociation of the remaining hydrates across a nonplanar and heterogeneous interface.

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.

Molecular insight into carbon dioxide hydrate formation from saline solution

Carbon dioxide hydrate has been intensively investigated in recent years because of its potential use as gas and heat storage materials. To understand the hydrate formation mechanisms, the crystallization of CO₂ hydrate from NaCl solutions was simulated at a molecular level. The influence of temperature, pressure, salt concentration and CO₂ concentration on CO₂ hydrate formation was evaluated. Results showed that the amount of the newly formed hydrate cages pressure went through a fast linear growth period followed by a relatively stable period. Pressure had little effect on CO₂ hydrate formation and temperature had a significant influence. The linear growth rate was greatly reduced as the temperature dropped from 255 to 235 K. The salt ion pairs could inhibit CO₂ hydrate formation, suggesting that we should choose the lower salinity areas if we want to storage CO₂ as gas hydrates in the seabed sediments. The observations in this study can provide theoretical support for the micro mechanism of hydrate formation, and provide a theoretical reference for the technology of hydrate based CO₂ storage.

In the process of the carbon dioxide hydrate formation in NaCl solution, it could form 5¹², 5¹²6² and 5¹²6³ cages, and the 5¹²6² cage and 5¹² cage number ratio was slightly above 3 : 1.

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.

Fluctuation-dissipation analysis of nonequilibrium thermal transport at the hydrate dissociation interface

Herein, a nonequilibrium molecular-dynamics simulation in an NVE ensemble was performed to investigate the spontaneous dissociation of methane-hydrate when it came in contact with liquid water. The nonequilibrium in the interface region is linked to the dissociation process of the hydrate near the interface according to the Onsager's hypothesis. The simulated thickness of the interface was found to be close to the acoustic phonon mean path of methane hydrate and agreed with the reference value. The normalized heat flow autocorrelation function was introduced to study fluctuation-dissipation in terms of the thickness and moving velocity of the interface and the Stefan number. This helped to clearly identify three distinct hydrate-decomposition regimes dominated by sensible heat, latent heat and an intrinsically unstable lattice framework. It was found that the fluctuation-dissipation theory could express the nonequilibrium nature in the front two stages before the threshold was reached, and the dissociation rate increased in the latter stage; this was different from the case of thermal-driven dissociation. The Stefan number decreased rapidly with dissociation in the initial stage and then fluctuated in the intermediate stage; this was analogous to the fluctuation characteristics of the heat flow autocorrelation function. The Stefan number effect shows that thermal dissipation drives the hydrate dissociation and correlates fluctuation to the nonequilibrium nature. It was also found that a small Stefan number was enough to break up the residual hydrate soon after the threshold was achieved. The transient interfacial thermal resistance of the interfacial region was obtained as a typical value in the range of 10^-7-10^-9 m² K W^-1. This justifies that fluctuation-dissipation exists in the nonequilibrium process of hydrate dissociation either in terms of heat flux, as observed in this study, or the diffusion of guest molecules, as reported in other studies.

Decomposition Characterizations of Methane Hydrate Confined inside Nanoscale Pores of Silica Gel below 273.15 K

The formation and decomposition of gas hydrates in nanoscale sediments can simulate the accumulation and mining process of hydrates. This paper investigates the Raman spectra of water confined inside the nanoscale pores of silica gel, the decomposition characterizations of methane hydrate that formed from the pore water, and the intrinsic relationship between them. The results show that pore water has stronger hydrogen bonds between the pore water molecules at both 293 K and 223 K. The structure of pore water is conducive to the nucleation of gas hydrate. Below 273.15 K, the decomposition of methane hydrate formed from pore water was investigated at atmospheric pressure and at a constant volume vessel. We show that the decomposition of methane hydrate is accompanied by a reformation of the hydrate phase: The lower the decomposition temperature, the more times the reformation behavior occurs. The higher pre-decomposition pressure that the silica gel is under before decomposition is more favorable to reformation. Thus, reformation is the main factor in methane hydrate decomposition in nanoscale pores below 273.15 K and is attributed to the structure of pore water. Our results provide experimental data for exploring the control mechanism of hydrate accumulation and mining.

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.

A modeling study of methane hydrate decomposition in contact with the external surface of zeolites

Methane hydrate dissociates on the external surface of siliceous zeolites with methane absorbed by the solid and water forming a liquid-like phase.

Effects of ensembles on methane hydrate nucleation kinetics

By performing molecular dynamics simulations to form a hydrate with a methane nano-bubble in liquid water at 250 K and 50 MPa, we report how different ensembles, such as the NPT, NVT, and NVE ensembles, affect the nucleation kinetics of the methane hydrate. The nucleation trajectories are monitored using the face-saturated incomplete cage analysis (FSICA) and the mutually coordinated guest (MCG) order parameter (OP). The nucleation rate and the critical nucleus are obtained using the mean first-passage time (MFPT) method based on the FS cages and the MCG-1 OPs, respectively. The fitting results of MFPT show that hydrate nucleation and growth are coupled together, consistent with the cage adsorption hypothesis which emphasizes that the cage adsorption of methane is a mechanism for both hydrate nucleation and growth. For the three different ensembles, the hydrate nucleation rate is quantitatively ordered as follows: NPT > NVT > NVE, while the sequence of hydrate crystallinity is exactly reversed. However, the largest size of the critical nucleus appears in the NVT ensemble, rather than in the NVE ensemble. These results are helpful for choosing a suitable ensemble when to study hydrate formation via computer simulations, and emphasize the importance of the order degree of the critical nucleus.

A Theoretical Study of the Hydration of Methane, from the Aqueous Solution to the sI Hydrate-Liquid Water-Gas Coexistence

Monte Carlo and molecular dynamics simulations were done with three recent water models TIP4P/2005 (Transferable Intermolecular Potential with 4 Points/2005), TIP4P/Ice (Transferable Intermolecular Potential with 4 Points/ Ice) and TIP4Q (Transferable Intermolecular Potential with 4 charges) combined with two models for methane: an all-atom one OPLS-AA (Optimal Parametrization for the Liquid State) and a united-atom one (UA); a correction for the C–O interaction was applied to the latter and used in a third set of simulations. The models were validated by comparison to experimental values of the free energy of hydration at 280, 300, 330 and 370 K, all under a pressure of 1 bar, and to the experimental radial distribution functions at 277, 283 and 291 K, under a pressure of 145 bar. Regardless of the combination rules used for σ_C,O, good agreement was found, except when the correction to the UA model was applied. Thus, further simulations of the sI hydrate were performed with the united-atom model to compare the thermal expansivity to the experiment. A final set of simulations was done with the UA methane model and the three water models, to study the sI hydrate-liquid water-gas coexistence at 80, 230 and 400 bar. The melting temperatures were compared to the experimental values. The results show the need to perform simulations with various different models to attain a reliable and robust molecular image of the systems of interest.

Effects of thermodynamic inhibitors on the dissociation of methane hydrate: a molecular dynamics study

We investigate the effects of methanol and NaCl, which are known as thermodynamic hydrate inhibitors, on the dissociation kinetics of methane hydrate in aqueous solutions by using molecular dynamics simulations. It is shown that the dissociation rate is not constant but changes with time. The dissociation rate in the initial stage is increased by methanol whereas it is decreased by NaCl. This difference arises from the opposite effects of the two thermodynamic inhibitors on the hydration free energy of methane. The dissociation rate of methane hydrate is increased by the formation of methane bubbles in the aqueous phase because the bubbles absorb surrounding methane molecules. It is found that both methanol and NaCl facilitate the bubble formation. However, their mechanisms are completely different from each other. The presence of ions enhances the hydrophobic interactions between methane molecules. In addition, the ions in the solution cause a highly non-uniform distribution of dissolved methane molecules. These two effects result in the easy formation of bubbles in the NaCl solution. In contrast, methanol assists the bubble formation because of its amphiphilic character.

Nucleation rate analysis of methane hydrate from molecular dynamics simulations

Clathrate hydrates are solid crystalline structures most commonly formed from solutions that have nucleated to form a mixed solid composed of water and gas. Understanding the mechanism of clathrate hydrate nucleation is essential to grasp the fundamental chemistry of these complex structures and their applications. Molecular dynamics (MD) simulation is an ideal method to study nucleation at the molecular level because the size of the critical nucleus and formation rate occur on the nano scale. Various analysis methods for nucleation have been developed through MD to analyze nucleation. In particular, the mean first-passage time (MFPT) and survival probability (SP) methods have proven to be effective in procuring the nucleation rate and critical nucleus size for monatomic systems. This study assesses the MFPT and SP methods, previously used for monatomic systems, when applied to analyzing clathrate hydrate nucleation. Because clathrate hydrate nucleation is relatively difficult to observe in MD simulations (due to its high free energy barrier), these methods have yet to be applied to clathrate hydrate systems. In this study, we have analyzed the nucleation rate and critical nucleus size of methane hydrate using MFPT and SP methods from data generated by MD simulations at 255 K and 50 MPa. MFPT was modified for clathrate hydrate from the original version by adding the maximum likelihood estimate and growth effect term. The nucleation rates calculated by MFPT and SP methods are within 5%, and the critical nucleus size estimated by the MFPT method was 50% higher, than values obtained through other more rigorous but computationally expensive estimates. These methods can also be extended to the analysis of other clathrate hydrates.

Reaction coordinate of incipient methane clathrate hydrate nucleation

Nucleation from solution is a ubiquitous phenomenon with relevance to myriad scientific disciplines, including pharmaceuticals, biomineralization, and disease. One prominent example is the nucleation of clathrate hydrates, multicomponent crystalline inclusion compounds relevant to the energy industry where they block pipelines and also constitute a potential vast energy resource. Despite their importance, the molecular mechanism of incipient hydrate formation remains unknown. Herein, we employ advanced molecular simulation tools (pB histogram, equilibrium path sampling) to provide a statistical-mechanical basis for extracting physical insight into the molecular steps by which clathrates form. Through testing the Mutually Coordinated Guest (MCG) order parameter, we demonstrate that both guest (methane) and host (water) structuring are crucial to accurately describe the nucleation of hydrates and determine a critical nucleus size of MCG-1 = 16 at 255 K and 500 bar. Equipped with a validated (and novel) reaction coordinate, subsequent equilibrium path sampling simulations yield the free energy barrier and nucleation rate. The resulting quantitative nucleation process is described by the MCG clustering mechanism. This constitutes a significant advance in the field of hydrates research, as the fitness of a molecular descriptor has never been statistically verified. More broadly, this work has significance to a wide range of multicomponent nucleation contexts wherein the formation mechanism depends on contributions from both solute and solvent.

Two-component order parameter for quantifying clathrate hydrate nucleation and growth

Methane clathrate hydrate nucleation and growth is investigated via analysis of molecular dynamics simulations using a new order parameter. This order parameter (OP), named the Mutually Coordinated Guest (MCG) OP, quantifies the appearance and connectivity of molecular clusters composed of guests separated by water clusters. It is the first two-component OP used for quantifying hydrate nucleation and growth. The algorithm for calculating the MCG OP is described in detail. Its physical motivation and advantages compared to existing methods are discussed.

Nucleation of gas hydrates within constant energy systems

The early stage of formation of gas hydrates has recently attracted attention as amorphous intermediate gas hydrate structures have been observed, apparently contrary to a classical model of nucleation and some experimental observations. To date, essentially all reported molecular simulations of the nucleation of gas hydrates have been under constant temperature conditions, which does not consider the possible impacts of heat transfer on the nucleation processes. Here we show, using constant energy molecular simulations, that the nuclei at an early stage of the hydrate formation have relatively more crystalline order in comparison with those observed in previous isothermal (NPT or NVT) work. The current work suggests a more transient role for intermediate amorphous structures during hydrate nucleation, thereby providing a stronger link between molecular simulation and experimental observations. Our NVE results nevertheless support the two-step nucleation mechanism proposed in previous simulation studies under constant temperature conditions which features the initial formation of amorphous hydrate-like structures.