A Molecular Dynamics Study on Xe/Kr Separation Mechanisms Using Crystal Growth Method

The separation of xenon/krypton gas mixtures is a valuable but challenging endeavor in the gas industry due to their similar physical characteristics and closely sized molecules. To address this, we investigated the effectiveness of the hydrate-based gas separation method for mixed Xe-Kr gas via molecular dynamics (MD) simulations. The formation process of hydrates facilitates the encapsulation of guest molecules within hydrate cages, offering a potential strategy for gas separation. Higher temperatures and pressures are advantageous for accelerating the hydrate growth rate. The final occupancy of guest molecules and empty cages within 512, 51264, and all hydrate cages were thoroughly examined. An increase in the pressure and temperature enhanced the occupancy rates of Xe in both 512 and 51264 cages, whereas elevated pressure alone improved the occupancy of Kr in 51264 cages. However, the impact of temperature and pressure on Kr occupancy within 512 cages was found to be minimal. Elevated temperature and pressure resulted in a reduced occupancy of empty cages. Predominantly, 51264 cages were occupied by Xe, whereas Kr showed a propensity to occupy the 512 cages. With increasing simulated pressure, the final occupancy of Xe molecules in all cages rose from 0.37 to 0.41 for simulations at 260 K, while the final occupancy of empty cages decreased from 0.24 to 0.2.

Unusual species of methane hydrate detected in nanoporous media using solid state 13C NMR

Methane is considered to be a cubic structure I (CS-I) clathrate hydrate former, although in a number of instances, small amounts of structure II (CS-II) clathrate hydrate have been transiently observed as well. In this work, solid-state magic angle spinning 13C NMR spectra of methane hydrate formed at low temperatures inside silica-based nanoporous materials with pores in the range of 3.8-20.0 nm (CPG-20, Vycor, and MCM-41) show methane in several different environments. In addition to methane encapsulated in the dodecahedral 512 (D) and tetrakaidecahedral 51262 (T) cages typical of the CS-I clathrate hydrate phase, methane guests in pentakaidecahedral 51263 (P) and hexakaidecahedral 51264 (H) cages are also identified, and these appear to be stabilized for extended periods of time. The ratio of methane guests among the D and T cages determined from the line intensities is significantly different from that of bulk CS-I samples and indicates that both CS-I and CS-II are present as the dominant species. This is the first observation of methane in P cages, and the possible structures in which they could be present are discussed. Broad and relatively strong methane peaks, which are also observed in the spectra, can be related to methane dissolved in an amorphous component of water adjacent to the pore walls. Nanoconfinement and interaction with the pore walls clearly have a strong influence on the hydrate formed and may reflect species present in the early stages of hydrate growth.

Insights into the mechanical stability of tetrahydrofuran hydrates from experimental, machine learning, and molecular dynamics perspectives

Natural gas hydrates (NGHs) hold immense potential as a future energy resource and for sustainable applications such as gas capture and storage. Due to the challenging formation conditions, however, their mechanical properties remain poorly understood. Herein, the mechanical characteristics of tetrahydrofuran (THF) hydrates, a proxy for methane hydrates, were investigated at different ice contents, strain rates, and temperatures using uniaxial compressive experiments. The results unveil a distinct behavior in the peak strength of THF hydrates with a varying ice content, strain rate and temperature, exhibiting an increase as the strain rate and temperature decrease, in contrast to the peak strength-strain rate relationship observed in polycrystalline ice. Based on the experimental data, four machine learning (ML) models including extreme gradient boosting (XGboost), multilayer perceptron (MLP), gradient boosting decision tree (GBDT) and decision tree (DT) were developed to predict the peak strength. The XGboost model demonstrates superior predictive performance, emphasizing the significant influence of ice content and temperature on the peak strength of hydrates. Furthermore, molecular dynamics (MD) simulations were employed to gain insights into the dissociation and formation processes of clathrate cages, as well as phase transitions and amorphization occurring at grain boundaries (GBs) involving diverse unconventional clathrate cages, including 51265, 4151062, 4151064, 425861 and 425862, with 425861 and 425862 cages being predominant. This study enhances our understanding of the mechanical properties and deformation mechanisms of hydrates and provides a ML-based predictive framework for estimating the compressive strength of hydrates under diverse coupling conditions. The findings have significant implications for stability assessments of NGHs and the exploitation of NGH resources.

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.

Understanding why constant energy or constant temperature may affect nucleation behavior in MD simulations: A study of gas hydrate nucleation

Molecular dynamics simulations have been widely used in exploring the nucleation behavior of many systems, including gas hydrates. Gas hydrates are ice-like solids in which gas molecules are trapped in water cages. During hydrate formation, a considerable amount of heat is released, and previous work has reported that the choice of temperature control scheme may affect the behavior of hydrate formation. The origins of this effect have remained an open question. To address this question, extensive NVE simulations and thermostatted (NPT and NVT) simulations with different temperature coupling strengths have been performed and compared for systems where a water nanodroplet is immersed in a H2S liquid. Detailed analysis of the hydrate structures and their mechanisms of formation has been carried out. Slower nucleation rates in NVE simulations in comparison to NPT simulations have been observed in agreement with previous studies. Probability distributions for various temperature measures along with their spatial distributions have been examined. Interestingly, a comparison of these temperature distributions reveals a small yet noticeable difference in the widths of the distributions for water. The somewhat reduced fluctuations in the temperature for the water species in the NVE simulations appear to be responsible for reducing the hydrate nucleation rate. We further conjecture that the NVE-impeded nucleation rate may be the result of the finite size of the surroundings (here the liquid H2S portion of the system). Additionally, a local spatial temperature gradient arising from the heat released during hydrate formation could not be detected.

Resolving Temperature-Dependent Hydrate Nucleation Pathway: The Role of "Transition Layer"

Understanding the nucleation of natural gas hydrate (NGH) at different conditions has important implications to NGH recovery and other industrial applications, such as gas storage and separation. Herein, vast numbers of hydrate nucleation events are traced via molecular dynamics (MD) simulations at different degrees of supercooling (or driving forces). Specifically, to precisely characterize a hydrate nucleus from an aqueous system during the MD simulation, we develop an evolutionary order parameter (OP) to recognize the nucleus size and shape. Subsequently, the free energy landscapes of hydrate during nucleation are explored by using the newly developed OP. The results suggest that at 270 K (or 0.92 Tm supercooling, where Tm is the melting point), the near-rounded nucleus prevails during the nucleation, as described from the classical nucleation theory. In contrast, at relatively strong driving forces of 0.85 and 0.88 Tm, nonclassical nucleation events arise. Specifically, the pathway toward an elongated nucleus becomes as important as the pathway toward a near-rounded nucleus. To explain the distinct nucleation phenomena at different supercoolings, a notion of a "transition layer" (or liquid-blob-like layer) is proposed. Here, the transition layer is to describe the interfacial region between the nucleus and aqueous solution, and this layer entails two functionalities: (1) it tends to retain CH4 depending on the degrees of supercooling and (2) it facilitates collision among CH4, which thus promote the incorporation of CH4 into nucleus. Our simulation indicates that compared to the near-rounded nucleus, the transition layer surrounding the elongated nucleus is more evident with the higher collision rate among CH4 molecules. As such, the transition layer tends to promote the elongated nucleus pathway, while offsetting the cost of larger surface free energy associated with the elongated nucleus. At 0.92 Tm, however, the transition layer gradually disappears, and classical nucleation events dominate. Overall, the notion of "transition layer" offers deeper insight into the NGH nucleation at different degrees of supercooling and could be extended to describe other types of hydrate nucleation.

Predicting mechanical properties of CO2hydrates: machine learning insights from molecular dynamics simulations

Understanding the mechanical properties of CO2hydrate is crucial for its diverse sustainable applications such as CO2geostorage and natural gas hydrate mining. In this work, classic molecular dynamics (MD) simulations are employed to explore the mechanical characteristics of CO2hydrate with varying occupancy rates and occupancy distributions of guest molecules. It is revealed that the mechanical properties, including maximum stress, critical strain, and Young's modulus, are not only affected by the cage occupancy rate in both large 51262and small 512cages, but also by the distribution of guest molecules within the cages. Specifically, the presence of vacancies in the 51262large cages significantly impacts the overall mechanical stability compared to 512small cages. Furthermore, four distinct machine learning (ML) models trained using MD results are developed to predict the mechanical properties of CO2hydrate with different cage occupancy rates and cage occupancy distributions. Through analyzing ML results, as-developed ML models highlight the importance of the distribution of guest molecules within the cages, as crucial contributor to the overall mechanical stability of CO2hydrate. This study contributes new knowledge to the field by providing insights into the mechanical properties of CO2hydrates and their dependence on cage occupancy rates and cage occupancy distributions. The findings have implications for the sustainable applications of CO2hydrate, and as-developed ML models offer a practical framework for predicting the mechanical properties of CO2hydrate in different scenarios.

Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure

Using molecular dynamics, the evolution of a metastable solution for "methane + water" was studied for concentrations of 3.36, 6.5, 9.45, 12.2, and 14.8 mol% methane at 270 K and 1 bar during 100 ns. We have found the intriguing behavior of the system containing over 10,000 water molecules: the formation of hydrate-like structures is observed at 6.5 and 9.45 mol% concentrations throughout the entire solution volume. This formation of "blobs" and the following amorphous hydrate were studied. The creation of a metastable methane solution through supersaturation is the key to triggering the collective process of hydrate formation under low pressure. Even the first stage (0-1 ns), before the first fluctuating cavities appear, is a collective process of H-bond network reorganization. The formation of fluctuation cavities appears before steady hydrate growth begins and is associated with a preceding uniform increase in the water molecule's tetrahedrality. Later, the constantly presented hydrate cavities become the foundation for a few independent hydrate nucleation centers, this evolution is consistent with the labile cluster and local structure hypotheses. This new mechanism of hydrogen-bond network reorganization depends on the entropy of the cavity arrangement of the guest molecules in the hydrate lattice and leads to hydrate growth.

A Clathrate Hydrate Structure Hidden in Plain Sight

Clathrate hydrates have diverse crystal structures, and among them, the three (sI, sII, and sH) most prevalent ones cover nearly all known structures, while the norm is to consider other structures only when specific guest molecules are present. Here we report the observation of a hidden clathrate structure: the tetragonal structure (TS-I) in commonly formed gas hydrates, as evidenced from molecular dynamics simulations. We show that when two (or more) sI crystal grains with different growth directions come into contact or when the growth of a sI crystal encounters geometrical frustration, the TS-I results as a cocrystal. We give evidence that TS-I may also play an important role in the combination and/or transition between sI and sII. These results imply that this previously neglected structure may be commonly present whenever sI or sII is formed. This hidden structure must be identified, experimentally and in simulations; confining the possible structures may hinder an in-depth understanding of clathrate hydrates.

Mixed Hydrate Nucleation: Molecular Mechanisms and Cage Structures

The molecular-level details of the formation of mixed gas hydrates remain elusive despite their significance for a variety of scientific and industrial applications. In this study, extensive molecular simulations have been performed to examine the behavior of CH₄/H₂S mixed hydrate nucleation utilizing two different simulation setups varying in compositions and temperatures. The observed behavior exhibits similar phenomenology across the various systems once differences in nucleation rates and guest uptake are accounted for. We find that CH₄ is always enriched in the hydrate phase while the aqueous phase is enriched in H₂S. Even with H₂S as a minor component (i.e., 10% mole fraction), the system can mirror the overall nucleation kinetics of pure H₂S hydrate systems with CH₄-dominant nuclei. Through analyses of cages and their transitions, nonstandard cages, particularly those with 12 faces (e.g., 5¹⁰6²), have been found to be key intermediate cage types in the early stage of nucleation. Additionally, we present previously unreported cage types comprising heptagonal faces (e.g., 5⁹6²7¹) as having a significant role in the early-stage gas hydrate structural transitions.

Three-body aggregation of guest molecules as a key step in methane hydrate nucleation and growth

Gas hydrates have an important role in environmental and astrochemistry, as well as in energy materials research. Although it is widely accepted that gas accumulation is an important and necessary process during hydrate nucleation, how guest molecules aggregate remains largely unknown. Here, we have performed molecular dynamics simulations to clarify the nucleation path of methane hydrate. We demonstrated that methane gather with a three-body aggregate pattern corresponding to the free energy minimum of three-methane hydrophobic interaction. Methane molecules fluctuate around one methane which later becomes the central gas molecule, and when several methanes move into the region within 0.8 nm of the potential central methane, they act as directional methane molecules. Two neighbor directional methanes and the potential central methane form a three-body aggregate as a regular triangle with a distance of ~6.7 Å which is well within the range of typical methane-methane distances in hydrates or in solution. We further showed that hydrate nucleation and growth is inextricably linked to three-body aggregates. By forming one, two, and three three-body aggregates, the possibility of hydrate nucleation at the aggregate increases from 3/6, 5/6 to 6/6. The results show three-body aggregation of guest molecules is a key step in gas hydrate formation.

Effect of Cage Occupancy on Stability and Decomposition of Methane Hydrate

Gas hydrates usually contain a certain number of empty cages that will both affect the hydrate stability and reduce the gas storage capacity. In this work, by MD simulations, we found that the hydrate stability is related to the cage occupancy, the empty cage types, and especially the distribution of empty cages. With the decrease of overall occupancy, the stability of hydrate becomes worse. Under the same overall occupancy, the more concentrated the empty cages are, the more unstable the hydrate is and hence the faster it decomposes. The methane molecules may migrate between distorted cages during the decomposition, resulting in a temporary increase in the stability of hydrate. Hydrates with different empty cage distributions show different decomposition mechanisms: when empty cages are concentrated, the melting rate is fast first due to the rapid decomposition of empty cages, but the remaining filled cages will reduce the melting rate; when empty cages are separated on the contrary, the early melting is slow because of the high local occupancy, and the following melting will be accelerated because of the high melting surface area. It indicates that the empty cage distribution plays a controlling role in hydrate decomposition kinetics at different stages.

On the size, shape and energetics of the hydration shell around alkanes

A large number of clathrate-like cages have been proposed as the very first hydration shell of alkanes. The cages include canonical structures commonly found in clathrate hydrates and many others, not previously reported, derived from the carbon fullerene cavities. These structures have a rich and variegated form, which can adapt to the shape and conformation of the solute. They avoid "wasting" hydrogen bonds, while minimizing the volume cage and maximizing the solute-solvent van der Waals interactions. DFT/M06-2X and MP2 ab initio calculations give comparable structural and energetic results although the latter predicts slightly larger cages for a given solute. It is shown that the van der Waals interactions are substantial and the large exoenergetic values found for isobutane and cyclopentane provide an explanation for the surprising high melting points of related hydrates at room pressure. The encaging enthalpy for various hydrocarbons is similar to the enthalpy of solution measured at a temperature just above the melting point of aqueous hydrocarbon solutions, thus indicating that water molecules should not deviate too much from the configuration with O-H bonds tangentially oriented with respect to the solute surface. The computed trend differs from the enthalpy of solution measured at room temperature, thus the very first hydration shell departs, up to a certain degree, from the clathrate-like structures.

Molecular Origins of Deformation in Amorphous Methane Hydrates

Water and methane can stay together under low temperature and high pressure in the forms of liquid solutions and crystalline solids. From liquid and gaseous states to crystalline solids or the contrary processes, amorphous methane hydrates can occur in these evolution scenarios. Herein, mechanical properties of amorphous methane hydrates are explored for the first time to bridge the gap between mechanical responses of monocrystalline and polycrystalline methane hydrates. Our results demonstrate that mechanical properties of amorphous methane hydrates are strongly governed by our original proposed order parameter, namely, normalized hydrogen-bond directional order parameter. Followed by this important achievement, a multistep deformation mechanism core is proposed to explain mechanical properties of amorphous methane hydrates. Through an extensive detailed analysis of amorphous methane hydrates, our simulation results not only greatly enlarge our fundamental understanding for mechanical responses of amorphous methane hydrates in geological systems but also offer a fresh perspective in structure-property topics of solid materials in future science and technology.

Iterative Cup Overlapping: An Efficient Identification Algorithm for Cage Structures of Amorphous Phase Hydrates

Molecular dynamics studies have revealed that the nucleation pathway of clathrate hydrates involves the evolution from amorphous to crystalline hydrates. In this study, complete cages are further classified into the standard edge-saturated cages (SECs) and nonstandard edge-saturated cages (non-SECs). Centered on studying the structure and evolution of non-SECs and SECs, we propose a novel and efficient algorithm, iterative cup overlapping (ICO), to monitor hydrate nucleation and growth in molecular simulations by identifying SECs and discuss possible causes of the instability of non-SECs. Manipulation of topological information makes it possible for ICO to avoid the repeated searches for identified cages and deduce all SECs with low time costs, improving the efficiency of identification significantly. The accuracy and efficiency of ICO were verified by comparing the identification results with other well-proven algorithms. Furthermore, it was found that non-SECs have short lifetimes and eventually decompose or reorganize into more stable structures. Some evidence suggests that the instability of non-SECs is closely related to the hydrogen-bonding configuration of water-ring aggregations that they contain. The spontaneous evolution of the hydrogen-bonding network into the tetrahedral network may be the main factor that causes the conversion of QWRAs and the evolution of non-SECs.

Unraveling nucleation pathway in methane clathrate formation

Methane clathrates are widespread on the ocean floor of the Earth. A better understanding of methane clathrate formation has important implications for natural-gas exploitation, storage, and transportation. A key step toward understanding clathrate formation is hydrate nucleation, which has been suggested to involve multiple evolution pathways. Herein, a unique nucleation/growth pathway for methane clathrate formation has been identified by analyzing the trajectories of large-scale molecular dynamics (MD) simulations. In particular, ternary water-ring aggregations (TWRAs) have been identified as fundamental structures for characterizing the nucleation pathway. Based on this nucleation pathway, the critical nucleus size and nucleation timescale can be quantitatively determined. Specifically, a methane hydration layer compression/shedding process is observed to be the critical step in (and driving) the nucleation/growth pathway, which is manifested through overlapping/compression of the surrounding hydration layers of the methane molecules, followed by detachment (shedding) of the hydration layer. As such, an effective way to control methane hydrate nucleation is to alter the hydration layer compression/shedding process during the course of nucleation.

Propane and Water: The Cooperativity of Unlikely Molecules to Form Clathrate Structures

Many unanswered questions still exist at the molecular level to understand the nucleation process and mechanism of clathrate hydrates, especially for larger guest molecules that would result in the structure II crystal. Here, we report on molecular dynamics simulations for propane and water to describe the molecular mechanism leading to a structure II system. Through a large number (30) of long (5 μs) and coupled annealing (20 μs) simulations, we detail the prenucleation, nucleation, growth, and annealing of propane clathrate hydrate structures at 250 K and 1800 bar. The results demonstrate the equal importance of the empty and occupied cages in the nucleation of propane hydrates. The critical nucleus size is identified to be eight cages. While separate distinct clusters may exist during the prenucleation period, only one survives to grow beyond the critical nucleus size, with the others remaining subcritical. From the annealing simulations, it is clear that solid rearrangement is a very slow process, and 20 μs is still not long enough to capture long-range ordering resembling the structure II crystal. These results, along with the developed analysis method, have a significant impact in advancing our understanding of the nucleation process for unlike molecules.

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 Simulation Study on the Growth of Structure II Nitrogen Hydrate

Crystal growth of N₂ hydrate in a three-phase system consisting of N₂ hydrate, liquid water, and gaseous N₂ was performed by molecular dynamics simulation at 260 K. Pressure influence on hydrate growth was evaluated. The kinetic properties including the growth rates and cage occupancies of the newly formed hydrate and the diffusion coefficient and concentration of N₂ molecules in liquid phase were measured. The results showed that the growth of N₂ hydrate could be divided into two stages where N₂ molecules in gas phase had to dissolve in liquid phase and then form hydrate cages at the liquid-hydrate interface. The diffusion coefficient and concentration of N₂ in liquid phase increased linearly with increasing pressure. As the pressure rose from 50 to 100 MPa, the hydrate growth rate kept increasing from 0.11 to 0.62 cages·ns^-1·Å^-2 and then dropped down to around 0.40 cages·ns^-1·Å^-2 once the pressure surpassed 100 MPa. During the hydrate formation, the initial sII N₂ hydrate phase set in the system served as a template for the subsequent growth of N₂ hydrate so that no new crystal structure was found. Analysis on the cage occupancies revealed that the amount of cages occupied by two N₂ molecules increased evidently when the pressure was above 100 MPa, which slowed down the growth rate of hydrate cages. Additionally, a small fraction of defective cages including two N₂ molecules trapped in 5¹²6⁵ cages and three N₂ molecules trapped 5¹²6⁸ cages was observed during the hydrate growth.

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'.

Cage fusion from bi-cages to tri-cages during nucleation of methane hydrate: a DFT-D simulation

Water-cage clusters encapsulating guest molecules are the basic components of hydrate crystal structures. Herein, we investigated the fusion process from bi-cages to tri-cages to probe the nucleation mechanism at the initial stage of CH4 hydrate formation by employing dispersion-corrected density functional theory. We found that tri-cages possess high stability by sharing three, rather than two, polygonal faces. In addition, any mono-cage combined with a nonstandard 4151062 cage could achieve considerable stability regardless of which face is shared; this finding illustrates that 4151062 cages are more likely to appear at the early stages of CH4 hydrate nucleation than other nonstandard cages. We then simulated the Raman spectra of CH4 molecules in water-cage to characterize the spectral characteristics of the CH4 hydrate. The C-H symmetric stretching frequency of encapsulated CH4 molecules red-shifted with increasing mono-cage size, which follows the prediction of the "loose cage-tight cage" model. The symmetric stretching vibrational frequencies of trapped CH4 molecules in the tri-cage revealed a clear red-shift compared with those in the component mono- and bi-cages. The cage fusion process and spectroscopic properties described in this work are expected to provide new atomistic insights into CH4 hydrates at the initial nucleation stage.

Methane Hydrate Nucleation within Elastic Confined Spaces: Suitable Spacing and Elasticity Can Accelerate the Nucleation

Elastic materials are candidates for process intensification of gas storage by forming gas hydrate. In this work, molecular dynamics simulations of hydrate nucleation in elastic silica double layers were performed to study the effect of elastic confined spaces on hydrate formation. It is found that in narrow confined spaces, hexagonal rings dominated the hydrogen bond network of water molecules established rapidly by a multisite nucleation mechanism. With molecules added, a bilayer water structure was formed finally because elastic space can adapt the volume expansion. In medium and wide confined spaces, hydrates were formed from a series of "pseudo cages" which are considered as precursors of complete hydrate cages. Moreover, the induction time for nucleation was a minimum when the elasticity of the silica layer changes: nucleation is fastest in the weak-elastic system. When the elasticity increases, it becomes hard to adapt the volume expansion during nucleation and also difficult to nucleate in very weak-elastic systems because of the fluctuation of the layers.

Hydrophobic Hydration and the Effect of NaCl Salt in the Adsorption of Hydrocarbons and Surfactants on Clathrate Hydrates

Adsorption of functional molecules on the surface of hydrates is key in the understanding of hydrate inhibitors. We investigate the adsorption of a hydrocarbon chain, nonionic and ionic surfactants, and ions at the hydrate-aqueous interface. Our results suggest a strong connection between the water ordering around solutes in bulk and the affinity for the hydrates surface. We distinguish two types of water ordering around solutes: (i) hydrophobic hydration where water molecules form a hydrogen bond network similar to clathrate hydrates, and (ii) ionic hydration where water molecules align according to the polarity of an ionic group. The nonionic surfactant and the hydrocarbon chain induce hydrophobic hydration and are favorably adsorbed on the hydrate surface. Adsorption of ions and the ionic headgroups on the hydrate surface is not favorable because ionic hydration and the hydrogen bond structure of hydrates are incompatible. The nonionic surfactant is adsorbed by the headgroup and tail while adsorption of the ionic surfactants is not favorable through the head. Water ordering is analyzed using the hydrogen bond and tetrahedral density profiles as a function of the distance to the chemical groups. The adsorption of solutes is studied through the free energy profiles as a function of the distance to the hydrate surface. Salt lowers the melting temperature of hydrates, disrupts hydrophobic hydration, reduces the solubility of solutes in the aqueous solution, and increases the propensity of solutes to be adsorbed on hydrate surfaces. Our studies are performed by the unbiased and steered molecular dynamics simulations. The results are in line with experiments on the effect of salt and alkanes in hydrate antiagglomeration.

Free energy landscape and molecular pathways of gas hydrate nucleation

Despite the significance of gas hydrates in diverse areas, a quantitative knowledge of hydrate formation at a molecular level is missing. The impediment to acquiring this understanding is primarily attributed to the stochastic nature and ultra-fine scales of nucleation events, posing a great challenge for both experiment and simulation to explore hydrate nucleation. Here we employ advanced molecular simulation methods, including forward flux sampling (FFS), p_B histogram analysis, and backward flux sampling, to overcome the limit of direct molecular simulation for exploring both the free energy landscape and molecular pathways of hydrate nucleation. First we test the half-cage order parameter (H-COP) which we developed for driving FFS, through conducting the p_B histogram analysis. Our results indeed show that H-COP describes well the reaction coordinates of hydrate nucleation. Through the verified order parameter, we then directly compute the free energy landscape for hydrate nucleation by combining both forward and backward flux sampling. The calculated stationary distribution density, which is obtained independently of nucleation theory, is found to fit well against the classical nucleation theory (CNT). Subsequent analysis of the obtained large ensemble of hydrate nucleation trajectories show that although on average, hydrate formation is facilitated by a two-step like mechanism involving a gradual transition from an amorphous to a crystalline structure, there also exist nucleation pathways where hydrate crystallizes directly, without going through the amorphous stage. The CNT-like free energy profile and the structural diversity suggest the existence of multiple active transition pathways for hydrate nucleation, and possibly also imply the near degeneracy in their free energy profiles among different pathways. Our results thus bring a new perspective to the long standing question of how hydrates crystallize.

Overview: Nucleation of clathrate hydrates

Molecular level knowledge of nucleation and growth of clathrate hydrates is of importance for advancing fundamental understanding on the nature of water and hydrophobic hydrate formers, and their interactions that result in the formation of ice-like solids at temperatures higher than the ice-point. The stochastic nature and the inability to probe the small length and time scales associated with the nucleation process make it very difficult to experimentally determine the molecular level changes that lead to the nucleation event. Conversely, for this reason, there have been increasing efforts to obtain this information using molecular simulations. Accurate knowledge of how and when hydrate structures nucleate will be tremendously beneficial for the development of sustainable hydrate management strategies in oil and gas flowlines, as well as for their application in energy storage and recovery, gas separation, carbon sequestration, seawater desalination, and refrigeration. This article reviews various aspects of hydrate nucleation. First, properties of supercooled water and ice nucleation are reviewed briefly due to their apparent similarity to hydrates. Hydrate nucleation is then reviewed starting from macroscopic observations as obtained from experiments in laboratories and operations in industries, followed by various hydrate nucleation hypotheses and hydrate nucleation driving force calculations based on the classical nucleation theory. Finally, molecular simulations on hydrate nucleation are discussed in detail followed by potential future research directions.

CH4 Hydrate Formation between Silica and Graphite Surfaces: Insights from Microsecond Molecular Dynamics Simulations

Microsecond simulations have been performed to investigate CH₄ hydrate formation from gas/water two-phase systems between silica and graphite surfaces, respectively. The hydrophilic silica and hydrophobic graphite surfaces exhibit substantially different effects on CH₄ hydrate formation. The graphite surface adsorbs CH₄ molecules to form a nanobubble with a flat or negative curvature, resulting in a low aqueous CH₄ concentration, and hydrate nucleation does not occur during 2.5 μs simulation. Moreover, an ordered interfacial water bilayer forms between the nanobubble and graphite surface thus preventing their direct contact. In contrast, the hydroxylated-silica surface prefers to be hydrated by water, with a cylindrical nanobubble formed in the solution, leading to a high aqueous CH₄ concentration and hydrate nucleation in the bulk region; during hydrate growth, the nanobubble is gradually covered by hydrate solid and separated from the water phase, hence slowing growth. The silanol groups on the silica surface can form strong hydrogen bonds with water, and hydrate cages need to match the arrangements of silanols to form more hydrogen bonds. At the end of the simulation, the hydrate solid is separated from the silica surface by liquid water, with only several cages forming hydrogen bonds with the silica surface, mainly due to the low CH₄ aqueous concentrations near the surface. To further explore hydrate formation between graphite surfaces, CH₄/water homogeneous solution systems are also simulated. CH₄ molecules in the solution are adsorbed onto graphite and hydrate nucleation occurs in the bulk region. During hydrate growth, the adsorbed CH₄ molecules are gradually converted into hydrate solid. It is found that the hydrate-like ordering of interfacial water induced by graphite promotes the contact between hydrate solid and graphite. We reveal that the ability of silanol groups on silica to form strong hydrogen bonds to stabilize incipient hydrate solid, as well as the ability of graphite to adsorb CH₄ molecules and induce hydrate-like ordering of the interfacial water, are the key factors to affect CH₄ hydrate formation between silica and graphite surfaces.

The effects of ice on methane hydrate nucleation: a microcanonical molecular dynamics study

The NVE simulations realize the ice shrinking when methane hydrate nucleates both heterogeneously and homogeneously.

Unraveling Mixed Hydrate Formation: Microscopic Insights into Early Stage Behavior

The molecular-level details of mixed hydrate nucleation remain unclear despite the broad implications of this process for a variety of scientific domains. Through analysis of mixed hydrate nucleation in a prototypical CH₄/H₂S/H₂O system, we demonstrate that high-level kinetic similarities between mixed hydrate systems and corresponding pure hydrate systems are not a reliable basis for estimating the composition of early stage mixed hydrate nuclei. Moreover, we show that solution compositions prior to and during nucleation are not necessarily effective proxies for the composition of early stage mixed hydrate nuclei. Rather, microscopic details, (e.g., guest-host interactions and previously neglected cage types) apparently play key roles in determining early stage behavior of mixed hydrates. This work thus provides key foundational concepts and insights for understanding mixed hydrate nucleation.

Evidence from mixed hydrate nucleation for a funnel model of crystallization

The molecular-level details of crystallization remain unclear for many systems. Previous work has speculated on the phenomenological similarities between molecular crystallization and protein folding. Here we demonstrate that molecular crystallization can involve funnel-shaped potential energy landscapes through a detailed analysis of mixed gas hydrate nucleation, a prototypical multicomponent crystallization process. Through this, we contribute both: (i) a powerful conceptual framework for exploring and rationalizing molecular crystallization, and (ii) an explanation of phenomenological similarities between protein folding and crystallization. Such funnel-shaped potential energy landscapes may be typical of broad classes of molecular ordering processes, and can provide a new perspective for both studying and understanding these processes.

Molecular Dynamics Simulation of the Crystal Nucleation and Growth Behavior of Methane Hydrate in the Presence of the Surface and Nanopores of Porous Sediment

The behavior of hydrate formation in porous sediment has been widely studied because of its importance in the investigation of reservoirs and in the drilling of natural gas hydrate. However, it is difficult to understand the hydrate nucleation and growth mechanism on the surface and in the nanopores of porous media by experimental and numerical simulation methods. In this work, molecular dynamics simulations of the nucleation and growth of CH4 hydrate in the presence of the surface and nanopores of clay are carried out. The molecular configurations and microstructure properties are analyzed for systems containing one H2O hydrate layer (System A), three H2O hydrate layers (System B), and six H2O hydrate layers (System C) in both clay and the bulk solution. It is found that hydrate formation is more complex in porous media than in the pure bulk solution and that there is cooperativity between hydrate growth and molecular diffusion in clay nanopores. The hydroxylated edge sites of the clay surface could serve as a source of CH4 molecules to facilitate hydrate nucleation. The diffusion velocity of molecules is influenced by the growth of the hydrate that forms a block in the throats of the clay nanopore. Comparing hydrate growth in different clay pore sizes reveals that the pore size plays an important role in hydrate growth and molecular diffusion in clay. This simulation study provides the microscopic mechanism of hydrate nucleation and growth in porous media, which can be favorable for the investigation of the formation of natural gas hydrate in sediments.

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.

Molecular dynamics study on the nucleation of methane + tetrahydrofuran mixed guest hydrate

The nucleation of methane (CH4), tetrahydrofuran (THF), and CH4 + THF hydrates are investigated by microsecond MD simulations. These three systems exhibit distinct structural developments in the aqueous phase quantified by the formation of cage structures of hydrogen bonded water molecules. The development of a cluster of cages in the CH4 system is limited by the scarce CH4 molecules in the solution, while in the THF system it is limited by the short lifetime of cages. In the CH4 + THF mixed guest system, a small cluster of caged CH4 molecules can be rapidly stabilized by abundant neighboring cages of THF molecules. Therefore, the induction time of the CH4 + THF mixed guest system is found to be significantly shorter than that of the pure CH4 and pure THF systems. Furthermore, the structure of cages found in the initially formed cage clusters are often different from the typical 5(12)6(n) (n = 0, 2, 3, 4) cages observed in clathrate hydrate systems. The cluster of cages may grow or transform into structure I or II clathrate hydrate in the later stages.