The mobility of water molecules through gas hydrates

Exploring vacancy defects in s-I clathrate hydrates

This study investigates the role of vacancy defects in s-I clathrate hydrate structures, particularly in the presence of ethylene oxide (EO) molecules, through first-principles calculations. The structural properties, formation energies, and guest-host interactions of these vacancy defects were examined in both periodic systems and finite-size clusters. Our findings demonstrate that EO molecules significantly stabilize vacancy defects via hydrogen bonding, especially when forming double hydrogen bonds with dangling hydrogens (d-Hs) arising from the molecular vacancy defect. The encapsulation of EO in defect-free cages and its interaction with dangling oxygens (d-Os) were also analyzed, highlighting the superior stabilizing effect of double hydrogen bonds. These results provide new insights into the behavior of vacancy defects in hydrate structures and the potential role of polar guest molecules in enhancing defect stability and facilitating hydrate formation processes.

Bjerrum defects in s-II gas hydrate

The energy and structure of Bjerrum defects in structure II gas hydrates were investigated by using first-principle calculations for finite-size clusters and periodic 3D lattice systems. The formation energies of these defects were calculated for the first time when the cages of the structure II structure were completely empty and the large cage was filled with a THF molecule. Analogous to findings in ice structures, one of the hydrogen atoms forming the D defect was noted to orient toward the cage. If the excess proton resides in the large cage, it acts as an attraction center for the polar guest molecule, i.e., THF. Therefore, the large cage guest THF molecule stabilizes the D/L defect pair and isolated D/L defect formation energies by forming hydrogen bonds with the D defect. In such cases, the defect structure representing a D/L defect pair containing a THF molecule interacting with one of the hydrogen atoms of the D defect mirrors the guest-induced ones. Notably, the classical Bjerrum defect and the guest-induced Bjerrum defect exhibit a similar phenomenon in defective structures. Contrary to existing literature, it is evident that guest-induced Bjerrum defects involve both the L and D components. The insights gained from this study could potentially offer an alternative perspective to understand various experimental observations, such as those related to dielectric and NMR properties.

CH/CH2 Group Clusters Doping Methane Hydrate Cages

Methane hydrate is a crystalline compound with methane molecules as guest species trapped in host water cages. In this study, we detected methane hydrate with water cages doped by (C_aromatic-H)₅ clusters, (C_aromatic-H)₆ clusters, and (3C_aliphatic-H₂ + 2H₂O) clusters using current spectroscopic techniques and differential scanning calorimetry (DSC). Methane molecules are trapped in the doped cages with type sI forming in nanoscale silica gel pores. The relative quantity ratio of host carbon to guest carbon in the doped hydrate sample reaches approximately 3.58. Methane hydrate doped by CH/CH₂ group clusters greatly improves the ability of the hydrate unit cell to store methane and increases the stability of methane hydrate. Fast proton diffusion in the doped methane hydrate was confirmed. The results of this study will provide efficient and energy saving technical support for disruptive changes in hydrate storage and transportation of methane gas technology with a doped and dense solid phase.

Study on the Growth Kinetics and Morphology of Methane Hydrate Film in a Porous Glass Microfluidic Device

Natural gas hydrates are widely considered one of the most promising green resources with large reserves. Most natural gas hydrates exist in deep-sea porous sediments. In order to achieve highly efficient exploration of natural gas hydrates, a fundamental understanding of hydrate growth becomes highly significant. Most hydrate film growth studies have been carried out on the surface of fluid droplets in in an open space, but some experimental visual works have been performed in a confined porous space. In this work, the growth behavior of methane hydrate film on pore interior surfaces was directly visualized and studied by using a transparent high-pressure glass microfluidic chip with a porous structure. The lateral growth kinetics of methane hydrate film was directly measured on the glass pore interior surface. The dimensionless parameter (−∆G/(RT)) presented by the Gibbs free energy change was used for the expression of driving force to explain the dependence of methane hydrate film growth kinetics and morphology on the driving force in confined pores. The thickening growth phenomenon of the methane hydrate film in micropores was also visualized. The results confirm that the film thickening growth process is mainly determined by water molecule diffusion in the methane hydrate film in glass-confined pores. The findings obtained in this work could help to develop a solid understanding on the formation and growth mechanisms of methane hydrate film in a confined porous space.

Characterizing key features in the formation of ice and gas hydrate systems

Crystallization in liquids is critical to a range of important processes occurring in physics, chemistry and life sciences. In this article, we review our efforts towards understanding the crystallization mechanisms, where we focus on theoretical modelling and molecular simulations applied to ice and gas hydrate systems. We discuss the order parameters used to characterize molecular ordering processes and how different order parameters offer different perspectives of the underlying mechanisms of crystallization. With extensive simulations of water and gas hydrate systems, we have revealed unexpected defective structures and demonstrated their important roles in crystallization processes. Nucleation of gas hydrates can in most cases be characterized to take place in a two-step mechanism where the nucleation occurs via intermediate metastable precursors, which gradually reorganizes to a stable crystalline phase. We have examined the potential energy landscapes explored by systems during nucleation, and have shown that these landscapes are rugged and funnel-shaped. These insights provide a new framework for understanding nucleation phenomena that has not been addressed in classical nucleation theory. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.

Ab Initio Studies on the Clathrate Hydrates of Some Nitrogen- and Sulfur-Containing Gases

Ab initio calculations are performed to investigate the host-guest interactions and multiple occupancies of some sulfur- (H₂S, CS₂) and nitrogen-containing (N₂, NO, and NH₃) molecules in dodecahedral, tetrakaidecahedral, and hexakaidecahedral water cages in this work. Five functionals in the framework of density functional theory are compared, and the M06-2X method appears to be the best to predict the binding energies as well as the geometries. Results show that N₂ and NO molecules are more stable in the 5¹²6⁴ cage, while NH₃ and H₂S prefer to stabilize in the 5¹²6² cage. This suggests that the sI hydrates of NH₃ and H₂S exhibit higher stability than the sII structures and that sII NO hydrate is more stable than sI NO hydrate. N₂ is found to be more stable in type II structure with single occupancy and to form type I hydrate with multiple occupancy, which is consistent with the experimental observations. As to the guest molecule CS₂, it may undergo severe structural deformation in the 5¹² and 5¹²6² cage. For multiple occupancies, the 5¹², 5¹²6², and 5¹²6⁴ water cages can trap up to two N₂ molecules, and the 5¹²6⁴ water cage can accommodate two H₂S molecules. This work is expected to provide new insight into the formation mechanism of clathrate hydrates for atmospherically important molecules.

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.

The nucleation of gas hydrates near silica surfaces

Understanding the nucleation and crystal growth of gas hydrates near mineral surfaces and in confinement are critical to the methane recovery from gas hydrate reservoirs. In this work, through molecular dynamics simulation studies, we present an exploration of the nucleation behavior of methane hydrates near model hydroxylated silica surfaces. Our simulation results indicate that the nucleation of methane hydrates can initiate from the silica surfaces despite of the structural mismatch of the two solid phases. A layer of intermediate half-cage structures was observed between the gas hydrate and silica surfaces, apparently helping to minimize the free energy penalty. These results have important implications to our understanding of the effects of solid surfaces on hydrate nucleation processes.

C-C stretching Raman spectra and stabilities of hydrocarbon molecules in natural gas hydrates: a quantum chemical study

The presence of specific hydrocarbon gas molecules in various types of water cavities in natural gas hydrates (NGHs) are governed by the relative stabilities of these encapsulated guest molecule-water cavity combinations. Using molecular quantum chemical dispersion-corrected hybrid density functional computations, the interaction (ΔE(host--guest)) and cohesive energies (ΔE(coh)), enthalpies, and Gibbs free energies for the complexes of host water cages and hydrocarbon guest molecules are calculated at the ωB97X-D/6-311++G(2d,2p) level of theory. The zero-point energy effect of ΔE(host-guest) and ΔE(coh) is found to be quite substantial. The energetically optimal host-guest combinations for seven hydrocarbon gas molecules (CH4, C2H6, C3H6, C3H8, C4H8, i-C4H10, and n-C4H10) and various water cavities (D, ID, T, P, H, and I) in NGHs are found to be CH4@D, C2H6@T, C3H6@T, C3H8@T, C4H8@T/P/H, i-C4H10@H, and n-C4H10@H, as the largest cohesive energy magnitudes will be obtained with these host-guest combinations. The stabilities of various water cavities enclosing hydrocarbon molecules are evaluated from the computed cohesive Gibbs free energies: CH4 prefers to be trapped in a ID cage; C2H6 prefer T cages; C3H6 and C3H8 prefer T and H cages; C4H8 and i-C4H10 prefer H cages; and n-C4H10 prefer I cages. The vibrational frequencies and Raman intensities of the C-C stretching vibrational modes for these seven hydrocarbon molecules enclosed in each water cavity are computed. A blue shift results after the guest molecule is trapped from gas phase into various water cages due to the host-guest interactions between the water cage and hydrocarbon molecule. The frequency shifts to the red as the radius of water cages increases. The model calculations support the view that C-C stretching vibrations of hydrocarbon molecules in the water cavities can be used as a tool to identify the types of crystal phases and guest molecules in NGHs.

New observations and insights into the morphology and growth kinetics of hydrate films

The kinetics of film growth of hydrates of methane, ethane, and methane-ethane mixtures were studied by exposing a single gas bubble to water. The morphologies, lateral growth rates, and thicknesses of the hydrate films were measured for various gas compositions and degrees of subcooling. A variety of hydrate film textures was revealed. The kinetics of two-dimensional film growth was inferred from the lateral growth rate and initial thickness of the hydrate film. A clear relationship between the morphology and film growth kinetics was observed. The shape of the hydrate crystals was found to favour heat or mass transfer and favour further growth of the hydrate film. The quantitative results on the kinetics of film growth showed that for a given degree of subcooling, the initial film thicknesses of the double hydrates were larger than that of pure methane or ethane hydrate, whereas the thickest hydrate film and the lowest lateral growth rate occurred when the methane mole fraction was approximately 0.6.

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.