How Properties of Solid Surfaces Modulate the Nucleation of Gas Hydrate

Surface modification of mesostructured cellular foam to enhance hydrogen storage in binary THF/H2 clathrate hydrate

This study introduces solid-state tuning of a mesostructured cellular foam (MCF) to enhance hydrogen (H2) storage in clathrate hydrates. Grafting of promoter-like molecules (e.g., tetrahydrofuran) at the internal surface of the MCF resulted in a substantial improvement in the kinetics of formation of binary H2-THF clathrate hydrate. Identification of the confined hydrate as sII clathrate hydrate and enclathration of H2 in its small cages was performed using XRD and high-pressure 1H NMR spectroscopy respectively. Experimental findings show that modified MCF materials exhibit a ∼1.3 times higher H2 storage capacity as compared to non-modified MCF under the same conditions (7 MPa, 265 K, 100% pore volume saturation with a 5.56 mol% THF solution). The enhancement in H2 storage is attributed to the hydrophobicity originating from grafting organic molecules onto pristine MCF, thereby influencing water interactions and fostering an environment conducive to H2 enclathration. Gas uptake curves indicate an optimal tuning point for higher H2 storage, favoring a lower density of carbon per nm2. Furthermore, a direct correlation emerges between higher driving forces and increased H2 storage capacity, culminating at 0.52 wt% (46.77 mmoles of H2 per mole of H2O and 39.78% water-to-hydrate conversions) at 262 K for the modified MCF material with fewer carbons per nm2. Notably, the substantial H2 storage capacity achieved without energy-intensive processes underscores solid-state tuning's potential for H2 storage in the synthesized hydrates. This study evaluated two distinct kinetic models to describe hydrate growth in MCF. The multistage kinetic model showed better predictive capabilities for experimental data and maintained a low average absolute deviation. This research provides valuable insights into augmenting H2 storage capabilities and holds promising implications for future advancements.

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.

Mechanical properties of amorphous CO2 hydrates: insights from molecular simulations

Understanding physicochemical properties of amorphous gas hydrate systems is of great significance to reveal structural stabilities of polycrystalline gas hydrate systems. Furthermore, amorphous gas hydrates can occur ordinarily in the nucleation events of gas hydrate systems. Herein, the mechanical properties of amorphous carbon dioxide hydrates are examined by means of all-atom classical molecular dynamic simulations. Our molecular simulation results reveal that mechanical strengths of amorphous carbon dioxide hydrates are evidently governed by temperatures, confining pressures, and ratios of water to carbon dioxide molecules. Notably, under compressive loads, amorphous carbon dioxide hydrates firstly exhibit monotonic strain hardening, followed by an interesting distinct phenomenon characterized by a steady flow stress at further large deformation strains. Furthermore, structural evolutions of amorphous carbon dioxide hydrates are analyzed on the basis of the N-Hbond DOP order parameter. These important findings can not only contribute to our understanding of the structural stabilities of amorphous gas hydrate systems, but also help to develop fundamental understandings about grain boundaries of gas hydrate systems.

Probing the critical nucleus size in tetrahydrofuran clathrate hydrate formation using surface-anchored nanoparticles

Controlling the formation of clathrate hydrates is crucial for advancing hydrate-based technologies. However, the microscopic mechanism underlying clathrate hydrate formation through nucleation remains poorly elucidated. Specifically, the critical nucleus, theorized as a pivotal step in nucleation, lacks empirical validation. Here, we report uniform nanoparticles, e.g., graphene oxide (GO) nanosheets and gold or silver nanocubes with controlled sizes, as nanoprobes to experimentally measure the size of the critical nucleus of tetrahydrofuran (THF) clathrate hydrate formation. The capability of the nanoparticles in facilitating THF clathrate hydrate nucleation displays generally an abrupt change at a nanoparticle-size-determined specific supercooling. It is revealed that the free-energy barrier shows an abrupt change when the nanoparticles have an approximately the same size as that of the critical nucleus. Thus, it is inferred that THF clathrate hydrate nucleation involves the creation of a critical nucleus with its size being inversely proportional to the supercooling. By proving the existence and determining the supercooling-dependent size of the critical nucleus of the THF clathrate hydrates, this work brings insights in the microscopic pictures of the clathrate hydrate nucleation.

A Review of the Effect of Porous Media on Gas Hydrate Formation

Most gas hydrates on the earth are in sediments and permafrost areas, and porous media are often used industrially as additives to improve gas hydrate formation. For further understanding its exploration and exploitation under natural conditions and its application in industry, it is necessary to study the effect of porous media on hydrate formation. The results show that the stacked porous media affects the phase equilibrium of hydrate formation depending on the competition water activity and large specific surface areas, while integrated porous media, such as metal foam, can transfer the hydration heat rapidly and moderate the hydrate phase equilibrium. A supersaturated metal-organic framework is able to significantly improve gas storage performance and can be used as a new material to promote hydrate formation. However, the critical particle size should be studied further for approaching the best promotion effect. In addition, together with the kinetic accelerators, porous media has a synergistic effect on gas hydrate formation. The carboxyl and hydroxyl groups on the surface of porous media can stabilize hydrate crystals through hydrogen bonding. However, the hydroxyl radicals on the silica surface inhibit the combination of CH₄ and free water, making the phase equilibrium conditions more demanding.

"Nanoreactors" for Boosting Gas Hydrate Formation toward Energy Storage Applications

Hydrogen and methane can be molecularly incorporated in ice-like water structures up to mass fractions of 4.3% and 13.3%, respectively. The resulting solid structures, called gas hydrates, offer great potential for the efficient storage of hydrogen and natural gas. However, slow gas encapsulation by bulk water hinders this application. Porous structures have been shown to effectively promote gas hydrate formation and are a potential enabler for the development of hydrate-based gas storage technologies. Here, we offer an insightful perspective on using porous structures as nanoreactors for achieving fast gas hydrate formation for gas storage applications. We critically discuss and elucidate the working mechanisms of nanoreactors and identify the criteria for efficient nanoreactors. Based on the concepts founded, we propose a theoretical framework for designing next-generation porous materials for delivering better promoting effects on gas 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.

An Integrated Experimental and Computational Platform to Explore Gas Hydrate Promotion, Inhibition, Rheology, and Mechanical Properties at McGill University: A Review

(1) Background: Gas hydrates are historically notable due to their prevalence and influence on operational difficulties in the oil and gas industry. Recently, new technologies involving the formation of gas hydrates to accomplish various applications have been proposed. This has created new motivation for the characterization of rheological and mechanical properties and the study of molecular phenomena in gas hydrates systems, particularly in the absence of oil and under pre-nucleation conditions. (2) Methodology: This work reviews advances in research on the promotion, inhibition, rheology, and mechanical properties of gas hydrates obtained through an integrated material synthesis-property characterization-multi-scale theoretical and computational platform at McGill University. (3) Discussion: This work highlights the findings from previous experimental work by our group and identifies some of their inherent physical limitations. The role of computational research methods in extending experimental results and observations in the context of mechanical properties of gas hydrates is presented. (4) Summary and Future perspective: Experimental limitations due to the length and time scales of physical phenomena associated with gas hydrates were identified, and future steps implementing the integrated experimental-computational platform to address the limitations presented here were outlined.

Evidence for Entropically Controlled Interfacial Hydration in Mesoporous Organosilicas

At aqueous interfaces, the distribution and dynamics of adsorbates are modulated by the behavior of interfacial water. Hydration of a hydrophobic surface can store entropy via the ordering of interfacial water, which contributes to the Gibbs energy of solute binding. However, there is little experimental evidence for the existence of such entropic reservoirs, and virtually no precedent for their rational design in systems involving extended interfaces. In this study, two series of mesoporous silicas were modified in distinct ways: (1) progressively deeper thermal dehydroxylation, via condensation of surface silanols, and (2) increasing incorporation of nonpolar organic linkers into the silica framework. Both approaches result in decreasing average surface polarity, manifested in a blue-shift in the fluorescence of an adsorbed dye. For the inorganic silicas, hydrogen-bonding of water becomes less extensive as the number of surface silanols decreases. Overhauser dynamic nuclear polarization (ODNP) relaxometry indicates enhanced surface water diffusivity, reflecting a loss of enthalpic hydration. In contrast, organosilicas show a monotonic decrease in surface water diffusivity with decreasing polarity, reflecting enhanced hydrophobic hydration. Molecular dynamics simulations predict increased tetrahedrality of interfacial water for the organosilicas, implying increased ordering near the nm-size organic domains (relative to inorganic silicas, which necessarily lack such domains). These findings validate the prediction that hydrophobic hydration at interfaces is controlled by the microscopic length scale of the hydrophobic regions. They further suggest that the hydration thermodynamics of structurally heterogeneous silica surfaces can be tuned to promote adsorption, which in turn tunes the selectivity in catalytic reactions.

Experimental Study of the Growth Kinetics of Natural Gas Hydrates in an Oil-Water Emulsion System

In order to explore the growth kinetics characteristics of NGH (natural gas hydrate) in an oil and gas mixed transportation pipeline and ensure the safe transportation of the pipeline, with the high-pressure hydrate experimental loop, an experimental study on the growth characteristics of NGH in an oil-water emulsion system was carried out, and the effects of pressure, flow rate, and water cut on the hydrate induction time, gas consumption, consumption rate, and hydrate volume fraction were explored, and important experimental rules were obtained. The experiment was divided into three stages: in the rapid formation stage of the hydrate, the temperature and gas consumption rose sharply, and the pressure dropped suddenly. The induction time decreased with the increase of pressure, flow rate, and water cut. The induction time of 6 MPa was 86.13 min, which was shortened by 39.68% compared with the induction time of 142.8 min of 5 MPa. The induction time of 1500 kg/h was 88.27 min, which was shorter by 13.91% than that 102.53 min of 550 kg/h. The induction time of 20% water cut was 58.53 min, which was shorter by 13.99% than that 68.4 min of 15% water cut. The gas consumption and hydrate volume fraction were both increased with the increase of pressure and water cut and decreased with the increase in the flow rate. In the whole process of the formation of NGH, the consumption rate first increased and then decreased. The pressure-drop and apparent viscosity increased with the increase of hydrate volume fraction in a certain range. The sensitivity analysis of hydrate induction time based on the standard regression coefficient method showed that the initial pressure played a major role, followed by the flow rate and the water cut. Based on the sensitivity analysis of hydrate volume fraction by the gray correlation method, it was found that the hydrate volume fraction had the closest relationship with the initial pressure, followed by the flow rate and the water cut. Finally, the empirical formulas of induction time and hydrate volume fraction in an oil-water emulsion system were established.

Interfacial study of clathrates confined in reversed silica pores

Storing methane in clathrates is one of the most promising alternatives for transporting natural gas (NG) as it offers similar gas densities to liquefied and compressed NG while offering lower safety risks. However, the practical use of clathrates is limited given the extremely low temperatures and high pressures necessary to form these structures. Therefore, it has been suggested to confine clathrates in nanoporous materials, as this can facilitate clathrate's formation conditions while preserving its CH4 volumetric storage. Yet, the choice of nanoporous materials to be employed as the clathrate growing platform is still rather arbitrary. Herein, we tackle this challenge in a systematic way by computationally exploring the stability of clathrates confined in alkyl-grafted silica materials with different pore sizes, ligand densities and ligand types. Based on our findings, we are able to propose key design criteria for nanoporous materials favoring the stability of a neighbouring clathrate phase, namely large pore sizes, high ligand densities, and smooth pore walls. We hope that the atomistic insight provided in this work will guide and facilitate the development of new nanomaterials designed to promote the formation of clathrates.

Experiment Investigation of SiO2 Containing Amino Groups as a Kinetic Promoter for CO2 Hydrates

To diminish the greenhouse effect by reducing CO₂ emission into the air based on a capture and sequestration method through hydrates, the thermodynamic and kinetic effects of additives on CO₂ hydrate formation under 1.5 MPa in the presence of 5, 6, 8, 10, and 20 wt % RNS-A (reactive SiO₂ containing amino groups) were studied, and the stirrer speed was set to 800 rpm. This paper calculated the gas consumption and explained the possible mechanisms of RNS-A on CO₂ hydrates. The results showed that RNS-A was a kinetic additive instead of a thermodynamic one. It was found that 5-10 wt % RNS-A all shortened the induction time of hydrates, but only 5 and 6 wt % RNS-A increased the gas consumption of CO₂ hydrates. Although we observed the shortest induction time at a 10 wt % RNS-A system, the lowest gas consumption indicated its weak CO₂ capture and storage ability. In addition, when the concentration was 6 wt %, RNS-A had the highest gas consumption and its reaction time was relatively short. Considering the induction time and gas consumption, 6 wt % RNS-A was the optimal RNS-A concentration for CO₂ capture and sequestration, which was the most suitable for practical applications.

Molecular Understanding of Homogeneous Nucleation of CO2 Hydrates Using Transition Path Sampling

Carbon dioxide hydrate is a solid built from hydrogen-bond stabilized water cages that encapsulate individual CO2 molecules. As potential candidates for reducing greenhouse gases, hydrates have attracted attention from both the industry and scientific community. Under high pressure and low temperature, hydrates are formed spontaneously from a mixture of CO2 and water via nucleation and growth. Yet, for moderate undercooling, i.e., moderate supersaturation, studying hydrate formation with molecular simulations is very challenging due to the high nucleation barriers involved. We investigate the homogeneous nucleation mechanism of CO2 hydrate as a function of temperature using transition path sampling (TPS), which generates ensembles of unbiased dynamical trajectories across the high barrier between the liquid and solid states. The resulting path ensembles reveal that at high driving force (low temperature), amorphous structures are predominantly formed, with 4151062 cages being the most abundant. With increasing temperature, the nucleation mechanism changes, and 51262 becomes the most abundant cage type, giving rise to the crystalline sI structure. Reaction coordinate analysis can reveal the most important collective variable involved in the mechanism. With increasing temperature, we observe a shift from a single feature (size of the nucleus) to a 2-dimensional (size and cage type) variable as the salient ingredient of the reaction coordinate, and then back to only the nucleus size. This finding is in line with the underlying shift from an amorphous to a crystalline nucleation channel. Modeling such complex phase transformations using transition path sampling gives unbiased insight into the molecular mechanisms toward different polymorphs, and how these are determined by thermodynamics and kinetics. This study will be beneficial for researchers aiming to produce such hydrates with different polymorphic forms.

Metal-Organic Framework HKUST-1 Promotes Methane Hydrate Formation for Improved Gas Storage Capacity

The large demand of natural gas consumption requires an effective technology to purify and store methane, the main component of natural gas. Metal-organic frameworks and gas hydrates are highly appealing materials for the efficient storage of industrially relevant gases, including methane. In this study, the methane storage capacity of the combination of methane hydrates and HKUST-1, a copper-based metal-organic framework, was studied using high pressure differential scanning calorimetry. The results show a synergistic effect, as the addition of HKUST-1 promoted hydrate growth, thus increasing the amount of water converted to hydrate from 5.9 to 87.2% and the amount of methane stored, relative to the amount of water present, from 0.55 to 8.1 mmol/g. The success of HKUST-1 as a promoter stems mainly from its large surface area, high thermal conductivity, and hydrophilicity. These distinctive properties led to a kinetically favorable decrease in hydrate growth induction period by 4.4 h upon the addition of HKUST-1. Powder X-ray diffraction and nitrogen isotherm suggests that the hydrate formation occurs primarily on the surface of HKUST-1 rather than within the pores. Remarkably, the HKUST-1 crystals show no significant changes in terms of structural integrity after many cycles of hydrate formation and dissociation, which results in the material having a long life cycle. These results confirm the beneficial role of HKUST-1 as a promoter for gas hydrate formation to increase methane gas storage capacity.

Enhanced methane gas storage in the form of hydrates: role of the confined water molecules in silica powders

Methane hydrates are promising materials for storage and transportation of natural gas; however, the slow kinetics and inefficient water to hydrate conversions impede its broad scale utilisation. The purpose of the present study is to demonstrate rapid (2-3 h) and efficient methane hydrate conversions by utilising the water molecules confined in the intra- and inter-granular space of silica powders. All the experiments were conducted with amorphous silica (10 g) powders of 2-30 μm; 10-20 nm grain size, to mimic the hydrate formations in fine sand and clay dominated environments under moderate methane pressure (7-8 MPa). Encasing of methane molecules in hydrate cages was confirmed by Raman spectroscopic (ex situ) and thermodynamic phase boundary measurements. The present studies reveal that the water to hydrate conversion is relatively slower in 10-20 nm grain size silica, although the nucleation event is rapid in both silicas. The process of hydrate conversion is vastly diffusion-controlled, and this was distinctly observed during the hydrate growth in nanosize silica.

Gas Hydrate Crystallization in Thin Glass Capillaries: Roles of Supercooling and Wettability

We designed and implemented an experimental methodology to investigate gas hydrate formation and growth around a water-guest meniscus in a thin glass capillary, thus mimicking pore-scale processes in sediments. The glass capillary acts as a high-pressure optical cell in a range of supercooling conditions from 0.1 °C, i.e., very close to hydrate dissociation conditions, to ∼35 °C, very near the metastability limit. Liquid or gaseous CO₂ is the guest phase in most of the experiments reported in this paper, and N₂ in a few of them. The setup affords detailed microscopic observation of the roles of the key parameters on hydrate growth and interaction with the substrate: supercooling and substrate wettability. At low supercooling (less than 0.5 °C), a novel hydrate growth process is discovered, which consists of a hollow crystal originating from the meniscus and advancing on the guest side along the glass, fed by a thick water layer sandwiched between the glass and this crystal.

Fundamental Investigation of the Effects of Modified Starch, Carboxymethylcellulose Sodium, and Xanthan Gum on Hydrate Formation under Different Driving Forces

The development of a new drilling fluid system with hydrate inhibition is of great significance for drilling safety in gas hydrate-bearing sediments. Considering the importance of the selection of a suitable thickener for drilling fluid systems under weak and strong driving forces, the hydrate inhibition of 0.1–0.5 wt% modified starch (MS), carboxymethylcellulose sodium (CMC), and xanthan gum (XG) aqueous solutions was studied. The applicability of these three thickeners were investigated through hydrate formation experiments, mesostructure observations, water activity tests, bubble retention observations, and rheological property tests. The results show that (1) under weak driving force, 0.3 wt% or higher concentration CMC and 0.3 wt% XG can almost completely inhibit hydrate formation due to the interactions between relatively small amounts of free water and CH4 molecules. Furthermore, the hydrate inhibition of higher XG concentrations was decreased due to their strong foam stability, leading to good contact between free water and CH4 molecules. Meanwhile, the hydrate inhibition of MS was weaker when compared with that of CMC and XG at the same concentrations. (2) Under strong driving force, the existence of the three 0.1–0.5 wt% thickeners could only slow down the hydrate formation rate, and hydrate inhibition due to XG was slightly better than that of the other two. This result implies that the effects of the different mesostructures on hydrate formation were severely weakened. Finally, (3) the tackifying effect of CMC was found to be stronger than that of XG and MS, and the rheological properties of the CMC solution were shown to be relatively weak compared to those of the XG and MS solutions; the CMC solution showed a more significant increase in viscosity with decreasing temperature, which is related to the differences in the mesostructures. Therefore, when the driving force of hydrate formation is relatively low, CMC is a good choice for the drilling fluid system when there is no requirement for cooling, while XG is more applicable for a system that needs cooling. In the case of a stronger driving force, XG is the optimal choice irrespective of whether the drilling fluid system needs cooling or not.

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.

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.

Formation of Methane Hydrate in the Presence of Natural and Synthetic Nanoparticles

Natural gas hydrates occur widely on the ocean-bed and in permafrost regions, and have potential as an untapped energy resource. Their formation and growth, however, poses major problems for the energy sector due to their tendency to block oil and gas pipelines, whereas their melting is viewed as a potential contributor to climate change. Although recent advances have been made in understanding bulk methane hydrate formation, the effect of impurity particles, which are always present under conditions relevant to industry and the environment, remains an open question. Here we present results from neutron scattering experiments and molecular dynamics simulations that show that the formation of methane hydrate is insensitive to the addition of a wide range of impurity particles. Our analysis shows that this is due to the different chemical natures of methane and water, with methane generally excluded from the volume surrounding the nanoparticles. This has important consequences for our understanding of the mechanism of hydrate nucleation and the design of new inhibitor molecules.

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.

Roles of Wettability and Supercooling in the Spreading of Cyclopentane Hydrate over a Substrate

We use transmission optical microscopy to observe cyclopentane hydrate growth in sub-mm, open glass capillaries, mimicking cylindrical pores. The capillary is initially loaded with water and the guest fluid (cyclopentane) and thus possesses three menisci, that between water and cyclopentane (CP) in the middle and two menisci with the vapors at the ends. At temperatures T below the equilibrium temperature T_eq ≈ 7 °C, the hydrate nucleates on the water-CP meniscus, rapidly coating it with an immobile, polycrystalline crust. Continued movement of the other two menisci provides insights into hydrate growth mechanisms, via the consumption and displacement of the fluids. On water-wet glass, the subsequent growth consists of a hydrate "halo" creeping with an underlying water layer on the glass on the CP side of the meniscus. Symmetrically, on CP-wet glass (silane-treated), a halo and a CP layer grow on the water side of the interface. No halo is observed on intermediate wet glass. The halo consists of an array of large monocrystals, over a thick water layer at low supercooling (ΔT = T_eq - T below 5 K), and a finer, polycrystalline texture over a thinner water layer at higher ΔT. Furthermore, the velocity varies as ΔT^α, with α ≈ 2.7, making the early stages of growth very similar to gas hydrate crusts growing over water-guest interfaces. Beyond a length in the millimeter range, the halo and its water layer abruptly decelerate and thin down to submicron thickness. The halo passes through the meniscus with the vapor without slowing down or change of texture. A model of the mass balance of the fluids helps rationalize all of these observations.

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

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

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.

Size Effect of Silica Shell on Gas Uptake Kinetics in Dry Water

Two kinds of dry water (DW) particles are prepared by mixing water and hydrophobic silica particles with nanometer or micrometer dimensions, and the two DW particles are found to have similar size distributions regardless of the size of the silica shell. The CO2 uptake kinetics of DW with nanometer (nanoshell) and micrometer shells (microshell) are measured, and both uptake rate and capacity show the obvious size effect of the silica shell. The DW with a microshell possesses a larger uptake capacity, whereas the DW with a nanoshell has a faster uptake rate. By comparing the uptake kinetics of soluble NH3 and CO2 further, we found that the microshell enhances the stability and the dispersion degree of DW and the nanoshell offers a shorter path for the transit of guest gas into the water core. Furthermore, molecular dynamics simulation is introduced to illustrate the nanosize effect of the silica shell on the initial step of the gas uptake. It is found that the concentration of gas molecules close to the silica shell is higher than that in the bulk water core. With the increase in the size of the silica shell, the amount of CO2 in the silica shell decreases, and it is easier for the gas uptake to reach steady state.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Microscopic Origin of Strain Hardening in Methane Hydrate

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