Overview: Nucleation of clathrate hydrates

Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation. II. Effect of multiple occupancy

In this work, we determine the dissociation line of the nitrogen (N2) hydrate by computer simulation using the TIP4P/Ice model for water and the TraPPE force field for N2. This work is the natural extension of Paper I, in which the dissociation temperature of the N2 hydrate has been obtained at 500, 1000, and 1500 bar [Algaba et al., J. Chem. Phys. 159, 224707 (2023)] using the solubility method and assuming single occupancy. We extend our previous study and determine the dissociation temperature of the N2 hydrate at different pressures, from 500 to 4500 bar, taking into account the single and double occupancy of the N2 molecules in the hydrate structure. We calculate the solubility of N2 in the aqueous solution as a function of temperature when it is in contact with a N2-rich liquid phase and when in contact with the hydrate phase with single and double occupancy via planar interfaces. Both curves intersect at a certain temperature that determines the dissociation temperature at a given pressure. We observe a negligible effect of occupancy on the dissociation temperature. Our findings are in very good agreement with the experimental data taken from the literature. We have also obtained the driving force for the nucleation of the hydrate as a function of temperature and occupancy at several pressures. As in the case of the dissociation line, the effect of occupancy on the driving force for nucleation is negligible. To the best of our knowledge, this is the first time that the effect of the occupancy on the driving force for nucleation of a hydrate that exhibits sII crystallographic structure is studied from computer simulation.

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.

Solubility of carbon dioxide in water: Some useful results for hydrate nucleation

In this paper, the solubility of carbon dioxide (CO2) in water along the isobar of 400 bar is determined by computer simulations using the well-known TIP4P/Ice force field for water and the TraPPE model for CO2. In particular, the solubility of CO2 in water when in contact with the CO2 liquid phase and the solubility of CO2 in water when in contact with the hydrate have been determined. The solubility of CO2 in a liquid-liquid system decreases as the temperature increases. The solubility of CO2 in a hydrate-liquid system increases with temperature. The two curves intersect at a certain temperature that determines the dissociation temperature of the hydrate at 400 bar (T3). We compare the predictions with T3 obtained using the direct coexistence technique in a previous work. The results of both methods agree, and we suggest 290(2) K as the value of T3 for this system using the same cutoff distance for dispersive interactions. We also propose a novel and alternative route to evaluate the change in chemical potential for the formation of hydrates along the isobar. The new approach is based on the use of the solubility curve of CO2 when the aqueous solution is in contact with the hydrate phase. It considers rigorously the non-ideality of the aqueous solution of CO2, providing reliable values for the driving force for nucleation of hydrates in good agreement with other thermodynamic routes used. It is shown that the driving force for hydrate nucleation at 400 bar is larger for the methane hydrate than for the carbon dioxide hydrate when compared at the same supercooling. We have also analyzed and discussed the effect of the cutoff distance of dispersive interactions and the occupancy of CO2 on the driving force for nucleation of the hydrate.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log₁₀ (J/(m³ s)) = -7(5).

Investigation on the Removal Performances of Heavy Metal Copper (II) Ions from Aqueous Solutions Using Hydrate-Based Method

Since heavy metal ion-contaminated water pollutionis becoming a serious threat to human and aquatic lives, new methods for highly efficient removal of heavy metal ions from wastewater are important to tackle environmental problems and sustainable development. In this work, we investigate the removal performances of heavy metal copper (II) ions from aqueous solutions using a gas hydrate-based method. Efficient removal of heavy metal copper (II) ions from wastewater via a methane hydrate process was demonstrated. The influence of the temperature, hydration time, copper (II) ions concentration, and stirring rate on the removal of heavy metal copper (II) ions were evaluated. The results suggested that a maximum of 75.8% copper (II) ions were removed from aqueous solution and obtained melted water with 70.6% yield with a temperature of -2 °C, stirring speed 800 r/min, and hydration time of 4 h with aninitial copper concentration of 100 mg/L. The initial concentration of copper (II) ions in the aqueous solution could be increased to between 100 and 500 mg/L. Meanwhile, our study also indicated that 65.6% copper (II) ions were removed from aqueous solution and the yield of melted water with 56.7%, even with the initial copper concentration of 500 mg/L. This research work demonstrates great potential for general applicability to heavy metal ion-contaminated wastewater treatment and provides a reference for the application of the gas hydrate method in separation.

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

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

Double Life of Methanol: Experimental Studies and Nonequilibrium Molecular-Dynamics Simulation of Methanol Effects on Methane-Hydrate Nucleation

We have investigated systematically and statistically methanol-concentration effects on methane-hydrate nucleation using both experiment and restrained molecular-dynamics simulation, employing simple observables to achieve an initially homogeneous methane-supersaturated solution particularly favorable for nucleation realization in reasonable simulation times. We observe the pronounced "bifurcated" character of the nucleation rate upon methanol concentration in both experiments and simulation, with promotion at low concentrations and switching to industrially familiar inhibition at higher concentrations. Higher methanol concentrations suppress hydrate growth by in-lattice methanol incorporation, resulting in the formation of "defects", increasing the energy of the nucleus. At low concentrations, on the contrary, the detrimental effect of defects is more than compensated for by the beneficial contribution of CH₃ in easing methane incorporation in the cages or replacing it altogether.

Rate Prediction for Homogeneous Nucleation of Methane Hydrate at Moderate Supersaturation Using Transition Interface Sampling

The crystallization of methane hydrates via homogeneous nucleation under natural, moderate conditions is of both industrial and scientific relevance, yet still poorly understood. Predicting the nucleation rates at such conditions is notoriously difficult due to high nucleation barriers, and requires, besides an accurate molecular model, enhanced sampling. Here, we apply the transition interface sampling technique, which efficiently computes the exact rate of nucleation by generating ensembles of unbiased dynamical trajectories crossing predefined interfaces located between the stable states. Using an accurate atomistic force field and focusing on specific conditions of 280 K and 500 bar, we compute for nucleation directly into the sI crystal phase at a rate of ∼10^-17 nuclei per nanosecond per simulation volume or ∼10² nuclei per second per cm³, in agreement with consensus estimates for nearby conditions. As this is most likely fortuitous, we discuss the causes of the large differences between our results and previous simulation studies. Our work shows that it is now possible to compute rates for methane hydrates at moderate supersaturation, without relying on any assumptions other than the force field.

Effects of Layer-Charge Distribution of 2:1 Clay Minerals on Methane Hydrate Formation: A Molecular Dynamics Simulation Study

Molecular dynamics simulations were used to investigate the effects of the external surface of a 2:1 clay mineral with different charge amounts and charge locations on CH₄ hydrate formation. The results showed that 5¹², 5¹²6², 5¹²6³, and 5¹²6⁴ were formed away from the clay mineral surface. The surface of the clay mineral with high- and low-charge layers was occupied by Na⁺ to form various distributions of outer- and inner-sphere hydration structures, respectively. The adsorbed Na⁺ on the high-charge layer surface reduced the H₂O activity by disturbing the hydrogen bond network, resulting in low tetrahedral arrangement of H₂O molecules near the layer surface, which inhibited CH₄ hydrate formation. However, more CH₄ molecules were adsorbed onto the vacancy in the Si-O rings of a neutral-charge layer to form semicage structures. Thus, the order parameter of H₂O molecules near this surface indicated that the arrangement of H₂O molecules resulted in a more optimal tetrahedral structure for CH₄ hydrate formation than that near the negatively charged layer surface. Different nucleation mechanisms of the CH₄ hydrate on external surfaces of clay mineral models were observed. For clay minerals with negatively charged layers (i.e., high and low charge), the homogeneous nucleation of the CH₄ hydrate occurred away from the surface. For a clay mineral with a neutral-charge layer, the CH₄ hydrate could nucleate either in the bulk-like solution homogeneously or at the clay mineral-H₂O interface heterogeneously.

Gas hydrates in sustainable chemistry

This review includes the current state of the art understanding and advances in technical developments about various fields of gas hydrates, which are combined with expert perspectives and analyses.

Nucleation and dissociation of methane clathrate embryo at the gas-water interface

Among natural energy resources, methane clathrate has attracted tremendous attention because of its strong relevance to current energy and environment issues. Yet little is known about how the clathrate starts to nucleate and disintegrate at the molecular level, because such microscopic processes are difficult to probe experimentally. Using surface-specific sum-frequency vibrational spectroscopy, we have studied in situ the nucleation and disintegration of methane clathrate embryos at the methane-gas-water interface under high pressure and different temperatures. Before appearance of macroscopic methane clathrate, the interfacial structure undergoes 3 stages as temperature varies, namely, dissolution of methane molecules into water interface, formation of cage-like methane-water complexes, and appearance of microscopic methane clathrate, while the bulk water structure remains unchanged. We find spectral features associated with methane-water complexes emerging in the induction time. The complexes are present over a wide temperature window and act as nuclei for clathrate growth. Their existence in the melt of clathrates explains why melted clathrates can be more readily recrystallized at higher temperature, the so-called "memory effect." Our findings here on the nucleation mechanism of clathrates could provide guidance for rational control of formation and disintegration of clathrates.

In-situ observation for natural gas hydrate in porous medium: Water performance and formation characteristic

Extensive efforts have been made regarding gas hydrate sample reconstruction in the laboratory for a better understanding and development of natural gas resources. Magnetic resonance imaging (MRI) is a useful method for directly observing the reconstruction of methane hydrate, yet relevant studies remain limited. In this study, a 9.4-T 400-MHz MRI instrument was employed to investigate CH₄ hydrate formation in porous media involving various initial water saturation levels and sand diameters. Pressure histories and MRI signal variations were monitored to discuss the process of CH₄ hydrate growth, and the three main formation stages of induction, rapid growth, and slow formation were determined. Furthermore, the liquid water performance in MRI micro-images was analyzed to predict the characteristics of CH₄ hydrate formation. The results indicated that CH₄ hydrate formed in a spatially and temporally random manner and that pore plugging occurred owing to the residual water encased in grown hydrate. Additionally, phase saturations, water conversion percentages, and formation rates were defined to evaluate the effect of sand diameter and initial water saturation on CH₄ hydrate formation. With the reduction in the diameter of quartz glass beads from 400 μm to 100 μm, the average hydrate formation rate increased from 0.0010 min^-1 to 0.0034 min^-1, respectively. When the initial water saturation decreased to the optimized value (0.22 in this study), the water conversion percentage and hydrate saturation increased.

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

Probing diffusion dynamics during hydrate formation by high field NMR relaxometry and diffusometry

High-field nuclear magnetic resonance (NMR) relaxometry and diffusometry along with magnetic resonance imaging were used to monitor phase transition molecular dynamics during hydrate formation occurring in water droplets dispersed in liquid cyclopentane. 1D T₂ relaxation measurements indicate the extent of hydrate formation as well as a reduction in water droplet size with progression of hydrate growth. MRI intensity maps and T₂ relaxation maps indicate spatially dependent hydrate formation rates due to the heterogeneity of the system. Spectrally resolved diffusion measurements indicate a reduction in the porosity of the hydrate agglomerate as the hydrate shell increases in thickness. A novel signal rise observed in two dimensional T₁-T₂ relaxation correlation experiments indicates complex diffusion dynamics due to coupling between regions with varying relaxation and diffusion. These results indicate the ability to monitor hydrate growth and phase transition molecular dynamics due to evolution of the porous hydrate agglomerate by means of high-field NMR.

Clathrate hydrates in interstellar environment

Clathrate hydrates (CHs) are ubiquitous in earth under high-pressure conditions, but their existence in the interstellar medium (ISM) remains unknown. Here, we report experimental observations of the formation of methane and carbon dioxide hydrates in an environment analogous to ISM. Thermal treatment of solid methane and carbon dioxide-water mixture in ultrahigh vacuum of the order of 10^-10 mbar for extended periods led to the formation of CHs at 30 and 10 K, respectively. High molecular mobility and H bonding play important roles in the entrapment of gases in the in situ formed 5¹² CH cages. This finding implies that CHs can exist in extreme low-pressure environments present in the ISM. These hydrates in ISM, subjected to various chemical processes, may act as sources for relevant prebiotic molecules.

CH4 Nanobubbles on the Hydrophobic Solid-Water Interface Serving as the Nucleation Sites of Methane Hydrate

The surface hydrophobicity of solid particles plays a critical role in the nucleation of gas hydrate formation, and it was found that the hydrophobic surface will promote this nucleation process, but the underlying mechanism is still unveiled. Herein, we proposed for the first time our new theory that the formation of methane nanoscale gas bubbles on the hydrophobic surface provides the nuclei sites for further formation of methane hydrate. First, we studied the effect of hydrophobicity of particles on the nucleation of hydrate. It was found that the hydrophobic graphite and silica particles would promote the nucleation of hydrate, but the hydrophilic silica particles did not promote the methane hydrate nucleation. Then, we designed the atomic force microscopy experiment to explain this mechanism from a nanometer scale. The results showed that the methane nanobubbles were formed on the hydrophobic highly ordered pyrolytic graphite surface, but they were hard to form on the hydrophilic mica surface. These results indicated that the methane nanobubbles on the hydrophobic surface could provide the gas hydrate nucleation sites and may induce a rapid nucleation of methane hydrate.

Diffusion of gas mixtures in the sI hydrate structure

Replacing methane with carbon dioxide in gas hydrates has been suggested as a way of harvesting methane, while at the same time storing carbon dioxide. Experimental evidence suggests that this process is facilitated if gas mixtures are used instead of pure carbon dioxide. We studied the free energy barriers for diffusion of methane, carbon dioxide, nitrogen, and hydrogen in the sI hydrate structure using molecular simulation techniques. Cage hops between neighboring cages were considered with and without a water vacancy and with a potential inclusion of an additional gas molecule in either the initial or final cage. Our results give little evidence for enhanced methane and carbon dioxide diffusion if nitrogen is present as well. However, the inclusion of hydrogen seems to have a substantial effect as it diffuses rapidly and can easily enter occupied cages, which reduces the barriers of diffusion for the gas molecules that co-occupy a cage with hydrogen.

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.