Formation of methane nano-bubbles during hydrate decomposition and their effect on hydrate growth

Study of CO2-hydrate formation in contact with bulk nanobubbles: An investigation from experiment and molecular-dynamics simulations

HYPOTHESIS

Nanobubbles (NBs) have been extensively investigated as a sustainable promoter for gas hydrate nucleation, which also contribute to the hydrate memory effect. However, less attention afforded to their effects on the hydrate-growth process, thus lacking a complete perspective of the overall effects from NBs on hydrate formation. We hypothesize that their effect on CO2 hydrate growth may vary depending on the properties of NBs.

EXPERIMENTS AND SIMULATIONS

This study investigates CO2-hydrate nucleation and growth with a dual methodology. Laboratory experiments were conducted using bulk NBs generated either from CO2 hydrate dissociation or electric-field-based electrostriction in virgin water. Simultaneously, molecular dynamics simulations examined hydrate growth in contact with a single NB containing CO2 molecules.

FINDINGS

Experimental results indicate that NBs promote hydrate nucleation, with finer ones from electric fields leading to a slight promotion at the onset of hydrate growth, followed by inhibition. MD simulations reveal that while NBs can serve as a gas source for growth, denser NBs hinder the process due to stronger gas-water interactions and CO2 clustering. These results suggest that optimizing NB size and concentration is critical for maximizing gas-hydrate formation efficiency in industrial applications such as gas storage and carbon capture.

Molecular insights into methane hydrate dissociation: Role of methane nanobubble formation

Understanding the underlying physics of natural gas hydrate dissociation is necessary for efficient CH4 extraction and in the exploration of potential additives in the chemical injection method. Silica being "sand" is already present inside the reservoir, making the silica nanoparticle a potential green additive. Here, molecular dynamics (MD) simulations have been performed to investigate the dissociation of the CH4 hydrate in the presence and absence of ∼1, ∼2, and ∼3 nm diameter hydrophilic silica nanoparticles at 100 bar and 310 K. We find that the formation of a CH4 nanobubble has a strong influence on the dissociation rate. After the initial hydrate dissociation, the rate of dissociation slows down till the formation of a CH4 nanobubble. We find the critical concentration and size limit to form the CH4 nanobubble to be ∼0.04 mole fraction of CH4 and ∼40 to 50 CH4 molecules, respectively. The solubility of CH4 and the chemical potential of H2O and CH4 are determined via Gibbs ensemble Monte Carlo simulations. The liquid phase chemical potential of both H2O and CH4 in the presence and absence of the nanoparticle is nearly the same, indicating that the effect of this additive will not be significant. While the formation of the hydration shell around the nanoparticle via hydrogen bonding confirms the strength of interactions between the water molecules and the nanoparticle in our MD simulations, the contact of the nanoparticle with the interface is infrequent, leading to no explicit effect of the nanoparticle on the dynamics of methane hydrate dissociation.

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.

Design of eco-friendly antifreeze peptides as novel inhibitors of gas-hydration kinetics

In this study, peptides designed using fragments of an antifreeze protein (AFP) from the freeze-tolerant insect Tenebrio molitor, TmAFP, were evaluated as inhibitors of clathrate hydrate formation. It was found that these peptides exhibit inhibitory effects by both direct and indirect mechanisms. The direct mechanism involves the displacement of methane molecules by hydrophobic methyl groups from threonine residues, preventing their diffusion to the hydrate surface. The indirect mechanism is characterized by the formation of cylindrical gas bubbles, the morphology of which reduces the pressure difference at the bubble interface, thereby slowing methane transport. The transfer of methane to the hydrate interface is primarily dominated by gas bubbles in the presence of antifreeze peptides. Spherical bubbles facilitate methane migration and potentially accelerate hydrate formation; conversely, the promotion of a cylindrical bubble morphology by two of the designed systems was found to mitigate this effect, leading to slower methane transport and reduced hydrate growth. These findings provide valuable guidance for the design of effective peptide-based inhibitors of natural-gas hydrate formation with potential applications in the energy and environmental sectors.

Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate

In this work, the effects of finite size on the determination of the three-phase coexistence temperature (T3) of the carbon dioxide (CO2) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved (hydrate-aqueous solution-liquid CO2) are placed together in the same simulation box. By varying the number of molecules of each phase, it is possible to analyze the effect of simulation size and stoichiometry on the T3 determination. In this work, we have determined the T3 value at 8 different pressures (from 100 to 6000 bar) and using 6 different simulation boxes with different numbers of molecules and sizes. In two of these configurations, the ratio of the number of water and CO2 molecules in the aqueous solution and the liquid CO2 phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO2 in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO2 available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above T3, overestimating the T3 value obtained from direct coexistence simulations. The simulation results obtained for the CO2 hydrate confirm the sensitivity of T3 depending on the size and composition of the system, explaining the discrepancies observed in the original work by Míguez et al. [J. Chem Phys. 142, 124505 (2015)]. Non-stoichiometric configurations with larger unit cells show a convergence of T3 values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.

Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate

Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.

Path-dependent morphology of CH4 hydrates and their dissociation studied with high-pressure microfluidics

Methane hydrates (MHs) have been considered a promising future energy source due to their vast resource volume and high energy density. Understanding the behavior of MH formation and dissociation at the pore-scale and the effect of MH distribution on the gas-liquid two phase flow is of critical importance for designing effective production strategies from natural gas hydrate (NGH) reservoirs. In this study, we devised a novel high-pressure microfluidic chip apparatus that is capable of direct observation of MH formation and dissociation behavior at the pore-scale. MH nucleation and growth behavior at 10.0 MPa and dissociation via thermal stimulation with gas bubble generation and evolution were examined. Our experimental results reveal that two different MH formation mechanisms co-exist in pores: (a) porous-type MH with a rough surface formed from CH4 gas bubbles at the gas-liquid interface and (b) crystal-type MH formed from dissolved CH4 gas. The growth and movement of crystal-type MH can trigger the sudden nucleation of porous-type MH. Spatially, MHs preferentially grow along the gas-liquid interface in pores. MH dissociation under thermal stimulation practically generates gas bubbles with diameters of 20.0-200.0 μm. Based on a custom-designed image analysis technique, three distinct stages of gas bubble evolution were identified during MH dissociation via thermal stimulation: (a) single gas bubble growth with an expanding water layer at an initial slow dissociation rate, (b) rapid generation of clusters of gas bubbles at a fast dissociation rate, and (c) gas bubble coalescence with uniform distribution in the pore space. The novel apparatus designed and the image analysis technique developed in this study allow us to directly capture the dynamic evolution of the gas-liquid interface during MH formation and dissociation at the pore-scale. The results provide direct first-hand visual evidence of the growth of MHs in pores and valuable insights into gas-liquid two-phase flow behavior during fluid production from NGHs.

Einstein-Stokes relation for small bubbles at the nanoscale

As the physicochemical properties of ultrafine bubble systems are governed by their size, it is crucial to determine the size and distribution of such bubble systems. At present, the size or size distribution of nanometer-sized bubbles in suspension is often measured by either dynamic light scattering or the nanoparticle tracking analysis. Both techniques determine the bubble size via the Einstein-Stokes equation based on the theory of the Brownian motion. However, it is not yet clear to which extent the Einstein-Stokes equation is applicable for such ultrafine bubbles. In this work, using atomic molecular dynamics simulation, we evaluate the applicability of the Einstein-Stokes equation for gas nanobubbles with a diameter less than 10 nm, and for a comparative analysis, both vacuum nanobubbles and copper nanoparticles are also considered. The simulation results demonstrate that the diffusion coefficient for rigid nanoparticles in water is found to be highly consistent with the Einstein-Stokes equation, with slight deviation only found for nanoparticle with a radius less than 1 nm. For nanobubbles, including both methane and vacuum nanobubbles, however, large deviation from the Einstein-Stokes equation is found for the bubble radius larger than 3 nm. The deviation is attributed to the deformability of large nanobubbles that leads to a cushioning effect for collision-induced bubble diffusion.

Molecular dynamics of the spontaneous generation mechanism of natural gas hydrates during methane nanobubble rupture

Natural gas hydrates have garnered significant attention as a potential new source of alternative energy, and understanding their formation mechanism is of paramount importance for efficient utilization and pipeline transportation. However, there is no consensus among academics on the formation mechanism of natural gas hydrates. In this paper, we propose a method for promoting the rapid formation of natural gas hydrates based on the rupture of methane nanobubbles, which creates local high temperature and pressure to facilitate the mixing of methane and water. The rapid decrease in system temperature and pressure during the process further enhances the formation of gas hydrates. Using molecular dynamics simulations, we theoretically verify the formation of natural gas hydrates. Our results indicate that the instantaneous rupture of methane nanobubbles induced by shock waves leads to a dramatic increase in the local molecular motion velocity around the bubbles. This results in extreme local high temperature and high pressure, leading to complete mixing of methane and water and rapid formation of gas hydrates during the cooling and pressure drop of the mixture. We confirm our findings by analyzing F3-order parameters, F4-order parameters, and water cage statistics.

Fast Formation of Hydrate Induced by Micro-Nano Bubbles: A Review of Current Status

Hydrate-based technologies have excellent application potential in gas separation, gas storage, transportation, and seawater desalination, etc. However, the long induction time and the slow formation rate are critical factors affecting the application of hydrate-based technologies. Micro-nano bubbles (MNBs) can dramatically increase the formation rate of hydrates owing to their advantages of providing more nucleation sites, enhancing mass transfer, and increasing the gas–liquid interface and gas solubility. Initially, the review examines key performance MNBs on hydrate formation and dissociation processes. Specifically, a qualitative and quantitative assembly of the formation and residence characteristics of MNBs during hydrate dissociation is conducted. A review of the MNB characterization techniques to identify bubble size, rising velocity, and bubble stability is also included. Moreover, the advantages of MNBs in reinforcing hydrate formation and their internal relationship with the memory effect are summarized. Finally, combining with the current MNBs to reinforce hydrate formation technology, a new technology of gas hydrate formation by MNBs combined with ultrasound is proposed. It is anticipated that the use of MNBs could be a promising sustainable and low-cost hydrate-based technology.

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.

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

Dynamic Dissociation Behaviors of sII Hydrates in Liquid Water by Heating: A Molecular Dynamics Simulation Approach

An understanding of the dynamic behavior of subtle hydrate dissociation in the liquid water phase is fundamental for gas production from marine hydrate reservoirs. Molecular dynamics simulations are performed in this study to investigate the dissociation kinetics of pure propane and binary propane + methane sII hydrates in a liquid water environment. The results show that faster hydrate dissociation rates are observed at higher initial temperatures. The hydrate phase dissociates from the cluster surface to the inside in a layer-by-layer manner under the simulation temperature conditions, which is similar to the behavior of sI hydrates and is independent of the hydrate crystal type. Compared to the binary sII hydrate, the pure sII hydrate dissociates more easily under the same initial temperature conditions, which can be attributed to the stabilizing effect of guest molecules in the hydrate cages. The empty cages collapse in one step, in contrast to the two-step pathway induced by the guest-host interaction. In addition, a hydrocarbon phase forms in the binary hydrate dissociation system instead of nanobubbles. These results can provide molecular-level insights into the dynamic mechanism of hydrate dissociation and theoretical guidance for gas recovery by thermal injection from marine hydrate reservoirs.

Molecular Insights into the Effect of Nitrogen Bubbles on the Formation of Tetrahydrofuran Hydrates

In this work, a molecular dynamics simulation was conducted to study the microscopic mechanism of how nitrogen bubbles affect the formation of THF hydrates at the molecular level. The results obtained reveal that the nitrogen bubble can promote the formation of THF hydrates. In the system with a nitrogen bubble, more THF-filled cages were generated, and the crystal structure was more orderly. The promotion of nitrogen bubbles on hydrate crystallization comes from the dissolution of nitrogen molecules. Some of dissolved nitrogen molecules can be enclosed in small hydrate cages near the nitrogen bubble, which can serve as stable sites for hydrate crystal growth, resulting in the fact that THF-filled cages connected with N₂-filled cages are much more stable and have a long lifetime. The results in this work can help to understand the promotion effect of micro- and nano-air bubbles on the crystallization of THF hydrates.

Molecular Insights into Factors Affecting the Generation Behaviors, Dynamic Properties, and Interfacial Structures of Methane Gas Bubbles

Molecular dynamics simulations were performed to study the effects of temperatures, pressures, and methane mole fractions on the generation behaviors, dynamic properties, and interfacial structures of methane gas bubbles. Methane gas bubbling can be promoted by high temperatures and high mole fractions of methane, which come from the generation of larger methane clusters in solution. Bubbles were found to be highly dynamic, with more methane molecules exchanging between bubbles and the surrounding solution at high pressures and in systems with high mole fractions of methane. The interfacial structures between bubbles and the surrounding solution were rough at a molecular level, and the roughness of the outermost methane and water molecules was high at high temperatures, low pressures, and in systems with high methane mole fractions. The dissolution of methane molecules depended on the interactions between the outermost methane and water molecules, which would become stronger with decreasing temperatures, increasing pressures, and decreasing methane mole fractions. The results obtained can help in understanding both the generation behaviors of bubbles when gas hydrates decompose and the re-nucleation behaviors of gas hydrates in the presence of bubbles.

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.

Molecular simulations on the stability and dynamics of bulk nanobubbles in aqueous environments

Nanobubbles have attracted significant attention due to their unexpectedly long lifetimes and stabilities in liquid solutions. However, explanations for the unique properties of nanobubbles at the molecular scale are somewhat controversial. Of special interest is the validity of the Young-Laplace equation in predicting the inner pressure of such bubbles. In this work, large-scale molecular dynamics simulations were performed to study the stability and diffusion of nanobubbles of methane in water. Two types of force field, atomistic and coarse-grained, were used to compare the calculated results. In accordance with predictions from the Young-Laplace equation, it was found that the inner pressure of the nanobubbles increased with decreasing nanobubble size. Consequently, a large pressure difference between the nanobubble and its surroundings resulted in the high solubility of methane molecules in water. The solubility was considered to enable nanobubble stability at exceptionally high pressures. Smaller bubbles were observed to be more mobile via Brownian motion. The calculated diffusion coefficient also showed a strong dependence on the nanobubble size. However, this active mobility of small nanobubbles also triggered a mutable nanobubble shape over time. Nanobubbles were also found to coalesce when they were sufficiently close. A critical distance between two nanobubbles was thus identified to avoid coalescence. These results provide insight into the behavior of nanobubbles in solution and the mechanism of their unique stability while withstanding high inner pressures.

Promoting Effect of Ultra-Fine Bubbles on CO2 Hydrate Formation

When gas hydrates dissociate into gas and liquid water, many gas bubbles form in the water. The large bubbles disappear after several minutes due to their buoyancy, while a large number of small bubbles (particularly sub-micron-order bubbles known as ultra-fine bubbles (UFBs)) remain in the water for a long time. In our previous studies, we demonstrated that the existence of UFBs is a major factor promoting gas hydrate formation. We then extended our research on this issue to carbon dioxide (CO2) as it forms structure-I hydrates, similar to methane and ethane hydrates explored in previous studies; however, CO2 saturated solutions present severe conditions for the survival of UFBs. The distribution measurements of CO2 UFBs revealed that their average size was larger and number density was smaller than those of other hydrocarbon UFBs. Despite these conditions, the CO2 hydrate formation tests confirmed that CO2 UFBs played important roles in the expression of the promoting effect. The analysis showed that different UFB preparation processes resulted in different promoting effects. These findings can aid in better understanding the mechanism of the promoting (or memory) effect of gas hydrate formation.

Structure and thermodynamics of empty clathrate hydrates below the freezing point of water

Recently prepared as a new H2O phase, ice XVI was obtained by degassing a Ne-sII clathrate hydrate under vacuum, however very little is known of that crystalline solid under temperatures (T ≤ 220 K) and pressures (p ≤ 5000 bar) relevant for the Earth's environment and geochemistry. In this work, atomically detailed calculations using long time-scale molecular simulations, seldom paralelled before, are employed to probe empty sII clathrate hydrates. It is found that the volumetric response to an applied pressure-temperature gradient is accurately described by the Parsafar and Mason equation of state with an accuracy of at least 99.7%. Structural deformation induced upon the crystals is interpreted by monitoring the unit cell length and isobaric thermal expansivity, whilst benchmarked against previous neutron diffraction measurements of ice XVI and hexagonal ice under room pressure conditions; a critical comparison is established with other sII guest occupied lattices (CH4, CO2 and CnH2n+2 with n = 2, 3, 4), often found in permafrost regions and in the margins of continental shelves. Such an analysis reveals that empty sII frameworks are slightly more stable to thermal deformation than their sI analogues and that hexagonal ice is the structurally most stable of the condensed H2O phases addressed here. Of paramount importance for the oil and natural gas industries, heat capacities obtained from enthalpy profiles are identical for the sI and sII empty clathrates up to 2000 bar and diverge by only ∼7.3% at 5000 bar. The canonical tetrahedral symmetry of water-bonded networks is analysed in terms of an angular and a distance order parameters, which are observed to decrease (increase) as pressure (temperature) increases (decreases). The results now being reported constitute a landmark for future studies dealing with high-pressure and very low-temperature conditions, characteristic of the Earth's permafrost environment and other planetary interiors, whilst contributing to expand our knowledge regarding the recently discovered ice XVI phase.

Contribution of Ultra-Fine Bubbles to Promoting Effect on Propane Hydrate Formation

To investigate experimentally how ultra-fine bubbles (UFBs) may promote hydrate formation, we examined the formation of propane (C₃H₈) hydrate from UFB-infused water solution using two preparation methods. In one method, we used C₃H₈-hydrate dissociated water, and in the other, C₃H₈-UFB-included water prepared with a generator. In both solutions, the initial conditions had a UFB number density of up to 10⁹ mL⁻¹. This number density decreased by only about a half when stored at room temperature for 2 days, indicating that enough amount of UFBs were stably present at least during the formation experiments. Compared to the case without UFBs, the nucleation probabilities within 50 h were ~1.3 times higher with the UFBs, and the induction times, the time period required for the bulk hydrate formation, were significantly shortened. These results confirmed that UFB-containing water promotes C₃H₈-hydrate formation. Combined with the UFB-stability experiments, we conclude that a high number density of UFBs in water contributes to the hydrate promoting effect. Also, consistent with previous research, the present study on C₃H₈ hydrates showed that the promoting effect would occur even in water that had not experienced any hydrate structures. Applying these findings to the debate over the promoting (or “memory”) effect of gas hydrates, we argue that the gas dissolution hypothesis is the more likely explanation for the effect.

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

Acoustic-propagation properties of methane clathrate hydrates from non-equilibrium molecular dynamics

Given methane hydrates' importance in marine sediments, as well as the widespread use of seabed acoustic-signaling methods in oil and gas exploration, the elastic characterization of these materials is particularly relevant. A greater understanding of the properties governing phonon, sound, and acoustic propagation would help to better classify methane-hydrate deposits, aiding in their discovery. Recently, we have published a new nonequilibrium molecular-dynamics (NEMD) methodology to recreate longitudinal and transverse perturbations, observing their propagation through a crystalline lattice by various metrics, to study the underlying S- and P-wave velocities (achieving excellent agreement with experiment) [Melgar et al., J. Phys. Chem. 122(5), 3006-3013 (2018); ibid.150, 084101 (2019)]. Here, we apply these NEMD methods to methane-clathrate systems to study acoustic-propagation characteristics, as well as the lattice elastic behavior. In so doing, we determine S- and P-wave velocities in excellent accord with experiment; we also ascertain the allowable magnitude range of acoustic perturbation and establish a threshold for lattice breakup and hydrate decomposition. Interestingly, upon dissociation, we observe the formation of methane nanobubbles, which agrees with previous studies on the microscopic fundamentals of hydrate dissociation by various means.

Removal of Fluorine from Wet-Process Phosphoric Acid Using a Solvent Extraction Technique with Tributyl Phosphate and Silicon Oil

The deep removal of fluorine from wet-process phosphoric acid is currently a very serious issue. In this paper, an efficient liquid-liquid separation method based on a bubble membrane was developed to solve this problem. Tributyl phosphate (TBP) and silicon oil (SIO) were used as the organic phase. The effects of the component proportion in the organic phase (TBP/SIO v/v), organic to aqueous phase ratio (O/A), pH, temperature, and reaction time on the extraction ratio were investigated. The extraction ratio of fluorine was 98.4% when using only one stage with the following conditions: 90 °C, pH -0.46, volume ratio (TBP/SIO v/v) of 7:3, phase ratio (O/A) of 1:5, stirring speed of 200 rpm, and reaction time of 50 min. Fourier-transform infrared spectroscopy and inverted fluorescence microscopy were used to investigate the reaction mechanism and reaction kinetics. In addition, the scrubbing and stripping process was investigated. When a 2 mol/L sodium hydroxide solution ([NaOH]) was used as the stripping agent with a phase ratio (O/A) of 1:10, a stirring speed of 200 rpm, and a reaction time of 30 min, a maximum stripping ratio of 90.1% was obtained.

New insights into decomposition characteristics of nanoscale methane hydrate below the ice point

In this paper, molecular dynamics simulation was used to study the decomposition process of nanoscale methane hydrate at 1 atm and 227 K. The results predict that methane hydrate decomposes into supercooled water (SCW) and methane gas and the resulting SCW turns into very high density amorphous ice (VHDA). The density of the VHDA is as high as 1.2-1.4 g cm^-3. The X-ray diffraction phase analysis showed that VHDA has a broad peak at 2θ of around 30°. The VHDA encapsulates the methane hydrate crystal so that the crystal can survive for a long time. The dissolved gas from the hydrate melt cannot escape out of the VHDA in a short time. The simulation results reveal new molecular insights into the decomposition behaviour of nanoscale methane hydrate below the ice point.

Formation of a nanobubble and its effect on the structural ordering of water in a CH4-N2-CO2-H2O mixture

The replacement of methane (CH4) from its hydrate by a mixture of nitrogen (N2) and carbon dioxide (CO2) involves the dissociation of methane hydrate leading to the formation of a CH4-N2-CO2-H2O mixture that can significantly influence the subsequent steps of the replacement process. In the present work, we study the evolution of dissolved gas molecules in this mixture by applying classical molecular dynamics simulations. Our study shows that a higher CO2 : N2 ratio in the mixture enhances the formation of nanobubbles composed of N2, CH4 and CO2 molecules. To understand how the CO2 : N2 ratio affects nanobubble nucleation, the distribution of molecules in the bubble formed is examined. It is observed that unlike N2 and CH4, the density of CO2 in the bubble reaches a maximum at the surface of the bubble. The accumulation of CO2 molecules at the surface makes the bubble more stable by decreasing the excess pressure inside the bubble as well as surface tension at its interface with water. It is found that a frequent exchange of gas molecules takes place between the bubble and the surrounding liquid and an increase in concentration of CO2 in the mixture leads to a decrease in the number of such exchanges. The effect of nanobubbles on the structural ordering of water molecules is examined by determining the number of water rings formed per unit volume in the mixture. The role of nanobubbles in water structuring is correlated to the dynamic nature of the bubble arising from the exchange of gas molecules between the bubble and the liquid.

Fast methane diffusion at the interface of two clathrate structures

Methane hydrates naturally form on Earth and in the interiors of some icy bodies of the Universe, and are also expected to play a paramount role in future energy and environmental technologies. Here we report experimental observation of an extremely fast methane diffusion at the interface of the two most common clathrate hydrate structures, namely clathrate structures I and II. Methane translational diffusion—measured by quasielastic neutron scattering at 0.8 GPa—is faster than that expected in pure supercritical methane at comparable pressure and temperature. This phenomenon could be an effect of strong confinement or of methane aggregation in the form of micro-nanobubbles at the interface of the two structures. Our results could have implications for understanding the replacement kinetics during sI–sII conversion in gas exchange experiments and for establishing the methane mobility in methane hydrates embedded in the cryosphere of large icy bodies in the Universe.

Methane dynamics at the interface of ice clathrate structures is expected to play a role in phenomena ranging from gas exchange to methane mobility in planetary cryospheres. Here, the authors observe extremely fast methane diffusion at the interface of the two most common clathrate hydrate structures.

Influence of Sodium Chloride on the Formation and Dissociation Behavior of CO2 Gas Hydrates

We present an experimental study on the formation and dissociation characteristics of carbon dioxide (CO₂) gas hydrates using Raman spectroscopy. The CO₂ hydrates were formed from sodium chloride/water solutions with salinities of 0–10 wt %, which were pressurized with liquid CO₂ in a stirred vessel at 6 MPa and a subcooling of 9.5 K. The formation of the CO₂ hydrate resulted in a hydrate gel where the solid hydrate can be considered as the continuous phase that includes small amounts of a dispersed liquid water-rich phase that has not been converted to hydrate. During the hydrate formation process we quantified the fraction of solid hydrate, x_H, and the fraction of the dispersed liquid water-rich phase, x_L, from the signature of the hydroxyl (OH)-stretching vibration of the hydrate gel. We found that the fraction of hydrate x_H contained in the hydrate gel linearly depends on the salinity of the initial liquid water-rich phase. In addition, the ratio of CO₂ and water was analyzed in the liquid water-rich phase before hydrate formation, in the hydrate gel during growth and dissociation, and after its complete dissociation again in the liquid water-rich phase. We observed a supersaturation of CO₂ in the water-rich phase after complete dissociation of the hydrate gel and were able to show that the excess CO₂ exists as dispersed micro- or nanoscale liquid droplets in the liquid water-rich phase. These residual nano- and microdroplets could be a possible explanation for the so-called memory effect.

A modeling study of methane hydrate decomposition in contact with the external surface of zeolites

Methane hydrate dissociates on the external surface of siliceous zeolites with methane absorbed by the solid and water forming a liquid-like phase.

Polarization response of clathrate hydrates capsulated with guest molecules

Clathrate hydrates are characterized by their water cages encapsulating various guest atoms or molecules. The polarization effect of these guest-cage complexes was studied with combined density functional theory and finite-field calculations. An addition rule was noted for these systems whose total polarizability is approximately equal to the polarizability sum of the guest and the cage. However, their distributional polarizability computed with Hirshfeld partitioning scheme indicates that the guest-cage interaction has considerable influence on their polarization response. The polarization of encapsulated guest is reduced while the polarization of water cage is enhanced. The counteraction of these two opposite effects leads to the almost unchanged total polarizability. Further analysis reveals that the reduced polarizability of encapsulated guest results from the shielding effect of water cage against the external field and the enhanced polarizability of water cage from the enhanced bonding of hydrogen bonds among water molecules. Although the charge transfer through the hydrogen bonds is rather small in the water cage, the polarization response of clathrate hydrates is sensitive to the changes of hydrogen bonding strength. The guest encapsulation strengthens the hydrogen bonding network and leads to enhanced polarizability.

Carbon dioxide induced bubble formation in a CH4–CO2–H2O ternary system: a molecular dynamics simulation study

The role of carbon dioxide in the formation of gas bubbles in a CH₄–CO₂–H₂O ternary system is studied using molecular dynamics simulations.