Melting and superheating of sI methane hydrate: Molecular dynamics study

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.

Molecular Dynamics Simulation of the Three-Phase Equilibrium Line of CO2 Hydrate with OPC Water Model

The three-phase coexistence line of the CO2 hydrate was determined using molecular dynamics (MD) simulations. By using the classical and modified Lorentz-Berthelot (LB) parameters, the simulations were carried out at 10 different pressures from 3 to 500 MPa. For the OPC water model, simulations with the classic and the modified LB parameters both showed negative deviations from the experimental values. For the TIP4P/Ice water model, good agreement with experimental equilibrium data can be achieved when the LB parameter is adjusted based on the solubility of CO2 in water. Our results also show that the influence of the water model on the equilibrium prediction is much larger than the CO2 model. Current simulations indicated that the H2O-H2O and H2O-CO2 cross-interactions' parameters might contribute equally to the accurate prediction of T3. According to our simulations, the prediction of T3 values showed relatively higher accuracy while using the combination of TIP4P/Ice water and EPM2 CO2 with modified LB parameter. Furthermore, varied χ values are recommended for accurate T3 estimation over a wide pressure range. The knowledge obtained in this study will be helpful for further accurate MD simulation of the process of CO2/CH4 replacement.

Solubility of Methane in Water: Some Useful Results for Hydrate Nucleation

In this paper, the solubility of methane in water along the 400 bar isobar is determined by computer simulations using the TIP4P/Ice force field for water and a simple LJ model for methane. In particular, the solubility of methane in water when in contact with the gas phase and the solubility of methane in water when in contact with the hydrate has been determined. The solubility of methane in a gas-liquid system decreases as temperature increases. The solubility of methane in a hydrate-liquid system increases with temperature. The two curves intersect at a certain temperature that determines the triple point T₃ at a certain pressure. We also determined T₃ by the three-phase direct coexistence method. The results of both methods agree, and we suggest 295(2) K as the value of T₃ for this system. We also analyzed the impact of curvature on the solubility of methane in water. We found that the presence of curvature increases the solubility in both the gas-liquid and hydrate-liquid systems. The change in chemical potential for the formation of hydrate is evaluated along the isobar using two different thermodynamic routes, obtaining good agreement between them. It is shown that the driving force for hydrate nucleation under experimental conditions is higher than that for the formation of pure ice when compared at the same supercooling. We also show that supersaturation (i.e., concentrations above those of the planar interface) increases the driving force for nucleation dramatically. The effect of bubbles can be equivalent to that of an additional supercooling of about 20 K. Having highly supersaturated homogeneous solutions makes possible the spontaneous formation of the hydrate at temperatures as high as 285 K (i.e., 10K below T₃). The crucial role of the concentration of methane for hydrate formation is clearly revealed. Nucleation of the hydrate can be either impossible or easy and fast depending on the concentration of methane which seems to play the leading role in the understanding of the kinetics of hydrate formation.

All-Atom Molecular Dynamics of Pure Water-Methane Gas Hydrate Systems under Pre-Nucleation Conditions: A Direct Comparison between Experiments and Simulations of Transport Properties for the Tip4p/Ice Water Model

(1) Background: New technologies involving gas hydrates under pre-nucleation conditions such as gas separations and storage have become more prominent. This has necessitated the characterization and modeling of the transport properties of such systems. (2) Methodology: This work explored methane hydrate systems under pre-nucleation conditions. All-atom molecular dynamics simulations were used to quantify the performance of the TIP4P/2005 and TIP4P/Ice water models to predict the viscosity, diffusivity, and thermal conductivity using various formulations. (3) Results: Molecular simulation equilibrium was robustly demonstrated using various measures. The Green–Kubo estimation of viscosity outperformed other formulations when combined with TIP4P/Ice, and the same combination outperformed all TIP4P/2005 formulations. The Green–Kubo TIP4P/Ice estimation of viscosity overestimates (by 84% on average) the viscosity of methane hydrate systems under pre-nucleation conditions across all pressures considered (0–5 MPag). The presence of methane was found to increase the average number of hydrogen bonds over time (6.7–7.8%). TIP4P/Ice methane systems were also found to have 16–19% longer hydrogen bond lifetimes over pure water systems. (4) Conclusion: An inherent limitation in the current water force field for its application in the context of transport properties estimations for methane gas hydrate systems. A re-parametrization of the current force field is suggested as a starting point. Until then, this work may serve as a characterization of the deviance in viscosity prediction.

The performance of OPC water model in prediction of the phase equilibria of methane hydrate

Molecular dynamics (MD) simulations were performed to determine the three-phase coexistence line of sI methane hydrates. The MD simulations were carried out at four different pressures (4, 10, 40 and 100 MPa) by using direct phase coexistence method. In current simulations, water was described by either TIP4P/Ice or OPC models and methane was described as a simple Lennard-Jones (LJ) interaction site. Lorentz-Berthelot combining rules were used to calculate the parameters of the cross interactions. For OPC model, positive deviations from the energetic Lorentz-Berthelot rule were also considered based on the solubility of methane in water. For TIP4P/Ice water model, the obtained three phase coexistence temperatures showed good agreement with experiment data at higher pressures, which is consistent with previous predictions. For OPC water model, simulations using the classic and the modified LB parameters both showed negative deviations to the experimental values. Our results also indicated that the deviation of the T3 prediction by OPC model not much correlated with the predicted melting point of ice. At 4 MPa, the modified OPC model showed outstanding prediction of hydrate equilibrium temperature, even better than the prediction by TIP4P/Ice. The relative higher accuracy in biomolecular MD of OPC model suggests that this model may have a better performance in hydrate MD simulations of biomolecule-based additives.

Thermodynamic Features of the Intensive Formation of Hydrocarbon Hydrates

This paper presents the results of a study on the influence of pressure and temperature of the gas–water medium on the process of hydrocarbon gas hydrate formation occurring at the phase interface. Herein, a mathematical model is proposed to determine the optimum ratios of pressure, gas temperatures, water, and gas bubble sizes in the bubbling, gas ejection, or mixing processes. As a result of our work, we determined that gas hydrate in these processes is formed at the gas–water interface, that is, on the boundary surface of gas bubbles. Moreover, there is a gas temperature range where the hydrate formation rate reaches its maximum. These study findings can be used to optimize various technological processes associated with the production of gas hydrates in the industry.

Does Confinement Enable Methane Hydrate Growth at Low Pressures? Insights from Molecular Dynamics Simulations

Natural methane hydrates are estimated to be the largest source of unexploited hydrocarbon fuel. The ideal conditions for methane hydrate formation are low temperatures and high pressures. On the other hand, recent experimental studies suggest that porous materials, thanks to their confinement effects, can enable methane hydrate formation at milder conditions, although there has not been a consensus on this. A number of studies have investigated methane hydrate growth in confinement by employing molecular simulations; however, these were carried out at either very high pressures or very low temperatures. Therefore, the effects of confinement on methane hydrate growth at milder conditions have not yet been elaborated by molecular simulations. In order to address this, we carried out a systematic study by performing molecular dynamics (MD) simulations of methane water systems. Using a direct phase coexistence approach, microsecond-scale MD simulations in the isobaric-isothermal (NPT) ensemble were performed in order to study the behavior of methane hydrates in the bulk and in confined nanospaces of hydroxylated silica pores at external pressures ranging from 1 to 100 bar and a simulation temperature corresponding to a 2 °C experimental temperature. We validated the combination of the TIP4P/ice water and TraPPE-UA methane models in order to correctly predict the behavior of methane hydrates in accordance to their phase equilibria. We also demonstrated that the dispersion corrections applied to short-range interactions lead to artificially induced hydrate growth. We observed that in the confinement of a hydroxylated silica pore, a convex-shaped methane nanobuble forms, and methane hydrate growth primarily takes place in the center of the pore rather than the surfaces where a thin water layer exists. Most importantly, our study showed that in the nanopores methane hydrate growth can indeed take place at pressures which would be too low for the growth of methane hydrates in the bulk.

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.

Phase Diagram of TIP4P/2005 Water at High Pressure

We report a new ice phase that forms spontaneously at the interface between ice VII and liquid water in molecular dynamics simulations of TIP4P/2005 water. The new phase is structurally quite similar to an ice phase originally found to be a precursor in the course of the homogeneous nucleation of ice VII in liquid water. A part of the water molecules in these ice phases can rotate easily because the number of hydrogen bonds in them is less than four, and thus they can be regarded as partial plastic phases. A rough estimate suggests that these phases are thermodynamically more stable than either ice VI or ice VII for 3 GPa < P < 18 GPa at T = 300 K. Although the partial plastic phases would be metastable states at any point in the phase diagram of real water, they might be realized experimentally with the aid of dopants and/or solid surfaces.

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.

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.

A theoretical study of the dissociation of the sI methane hydrate induced by an external electric field

Molecular dynamics simulations in the equilibrium isobaric-isothermal (NPT) ensemble were used to examine the strength of an external electric field required to dissociate the methane hydrate sI structure. The water molecules were modeled using the four-site TIP4P/Ice analytical potential and methane was described as a simple Lennard-Jones interaction site. A series of simulations were performed at T = 260 K with P = 80 bars and at T = 285 K with P = 400 bars with an applied electric field ranging from 1.0 V nm(-1) to 5.0 V nm(-1). For both (T,P) conditions, applying a field greater than 1.5 V nm(-1) resulted in the orientation of the water molecules such that an ice Ih-type structure was formed, from which the methane was segregated. When the simulations were continued without the external field, the ice-like structures became disordered, resulting in two separate phases: gas methane and liquid water.

Molecular dynamics simulation of halogen bonding in Cl2, BrCl, and mixed Cl2/Br2 clathrate hydrates

Clathrate hydrate phases of dihalogen molecules have properties that differ from those of other guest molecules of similar size. The water oxygen–chlorine distances in the structure I (sI) Cl2 hydrate are smaller than the sum of the van der Waals radii of oxygen and chlorine. Bromine hydrate forms a unique clathrate hydrate structure that is not seen in other guest substances. In mixed Cl2/Br2 structure I hydrate, the water oxygen–bromine distances are also smaller than the sum of the oxygen and bromine van der Waals radii. We previously studied the structure of three dihalogen clathrate hydrates using single crystal X-ray diffraction and described these structural features in terms of halogen bonding between the dihalogen and water molecules. In this work, we perform molecular dynamics simulations of cubic sI Cl2, mixed Cl2/Br2, and BrCl clathrate hydrate phases. We perform quantum chemical computations on the dihalogen molecules to determine the nature of σ-hole near the halogen atoms. We fit the electrostatic potential of the molecules to point charge models including dummy atoms that represent σ-holes adjacent to the halogen molecules. Molecular dynamics simulations are used to determine the lattice constants, radial distribution functions, and guest dynamics in these phases. We determine the effect of guest size and difference in halogen bonding on the properties of the clathrate hydrate phase. Simulations for the Cl2, BrCl, and mixed Cl2/Br2 hydrates are performed with small cages of the sI clathrate hydrate phases completely full or filled with experimental occupancies with Cl2 guests.

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.