Gas hydrates in confined space of nanoporous materials: new frontier in gas storage technology

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

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

A State-Dependent Elasto-Plastic Model for Hydrate-Bearing Cemented Sand Considering Damage and Cementation Effects

In order to optimize the efficiency and safety of gas hydrate extraction, it is essential to develop a credible constitutive model for sands containing hydrates. A model incorporating both cementation and damage was constructed to describe the behavior of hydrate-bearing cemented sand. This model is based on the critical state theory and builds upon previous studies. The damage factor Ds is incorporated to consider soil degradation and the reduction in hydrate cementation, as described by plastic shear strain. A computer program was developed to simulate the mechanisms of cementation and damage evolution, as well as the stress-strain curves of hydrate-bearing cemented sand. The results indicate that the model replicates the mechanical behavior of soil cementation and soil deterioration caused by impairment well. By comparing the theoretical curves with the experimental data, the compliance of the model was calculated to be more than 90 percent. The new state-dependent elasto-plastic constitutive model based on cementation and damage of hydrate-bearing cemented sand could provide vital guidance for the construction of deep-buried tunnels, extraction of hydrocarbon compounds, and development of resources.

Evidence of Spontaneous Formation of Two-Dimensional Amorphous Clathrates on Superhydrophilic Surfaces

Clathrate hydrates reserved in the seabed are often dispersed in the pores of coarse-grained sediments; hence, their formation typically occurs under nanoconfinement. Herein, we show the first molecular dynamics (MD) simulation evidence of the spontaneous formation of two-dimensional (2D) clathrate hydrates on crystal surfaces without conventional nanoconfinement. The kinetic process of 2D clathrate formation is illustrated via simulated single-molecule deposition. 2D amorphous patterns are observed on various superhydrophilic face-centered cubic surfaces. Notably, the formation of 2D amorphous clathrate can occur over a wide range of temperatures, even at room temperature. The strong water-surface interaction, the characteristic properties of guest-gas molecules, and the underlying surface structure dictate the formation of 2D amorphous clathrates. Semiquantitative phase diagrams of 2D clathrates are constructed where representative patterns of 2D clathrates for characteristic gas molecules on prototypical Pd(111) and Pt(111) surfaces are confirmed by independent MD simulations. A tunable pattern of 2D amorphous clathrates is demonstrated by changing the lattice strain of the underlying substrate. Moreover, ab initio MD simulations confirm the stability of 2D amorphous clathrate. The underlining physical mechanism for 2D clathrate formation on superhydrophilic surfaces is elucidated, which offers deeper insight into the crucial role of water-surface interaction.

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.

Metal-Organic Frameworks for Greenhouse Gas Applications

Using petrol to supply energy for a car or burning coal to heat a building generates plenty of greenhouse gas (GHG) emissions, including carbon dioxide (CO₂ ), water vapor (H₂ O), methane (CH₄ ), nitrous oxide (N₂ O), ozone (O₃ ), fluorinated gases. These up-and-coming metal-organic frameworks (MOFs) are structurally endowed with rigid inorganic nodes and versatile organic linkers, which have been extensively used in the GHG-related applications to improve the lives and protect the environment. Porous MOF materials and their derivatives have been demonstrated to be competitive and promising candidates for GHG separation, storage and conversions as they shows facile preparation, large porosity, adjustable nanostructure, abundant topology, and tunable physicochemical property. Enormous progress has been made in GHG storage and separation intrinsically stemmed from the different interaction between guest molecule and host framework from MOF itself in the recent five years. Meanwhile, the use of porous MOF materials to transform GHG and the influence of external conditions on the adsorption performance of MOFs for GHG are also enclosed. In this review, it is also highlighted that the existing challenges and future directions are discussed and envisioned in the rational design, facile synthesis and comprehensive utilization of MOFs and their derivatives for practical applications.

Durable Hydrate-phobic Coating with In Situ Self-Replenishing Hydrocarbon Barrier Films for Low Clathrate Hydrate Adhesion

Clathrate hydrate growth, deposition, and plug formation during oil and gas transportation causes blockage of pipelines. An effective strategy to solve this problem is to mitigate the hydrate formation and reduce its adhesion on pipe walls through a coating process. However, durability failure, corrosion, inability to self-heal, high cost, and strong hydrate adhesion remain unsolved issues. To address these challenges, in this work, we present an in situ self-replenishing nonfluorinated durable hydrate-phobic coating of candle soot particles. The candle soot coating reduces hydrate adhesion by promoting a thick barrier film of hydrocarbons between the hydrate and the soot coated substrate. The hydrocarbons permeating the soot coating display a high contact angle for water and inhibit the formation of water bridges between the hydrate and soot coated substrate. The spherical cyclopentane hydrate slides off easily on the candle soot coating inside the cyclopentane environment. The hydrate former, cyclopentane-water emulsion, and THF-water mixture have high contact angles as well as low hydrate adhesion on soot coating simultaneously. In addition, the coating is flow-induced long-term slippery, durable, low cost, anticorrosion, self-cleaning, and suitable for practical applications.

Delving into guest-free and He-filled sI and sII clathrate hydrates: a first-principles computational study

The dynamics of the formation of a specific clathrate hydrate as well as its thermodynamic transitions depend on the interactions between the trapped molecules and the host water lattice. The molecular-level understanding of the different underlying processes benefits not only the description of the properties of the system, but also allows the development of multiple technological applications such as gas storage, gas separation, energy transport, etc. In this work we investigate the stability of periodic crystalline structures, such as He@sI and He@sII clathrate hydrates by first-principles computations. We consider such host water networks interacting with a guest He atom using selected density functional theory approaches, in order to explore the effects on the encapsulation of a light atom in the sI/sII crystals, by deriving all energy components (guest-water, water-water, guest-guest). Structural properties and energies were first computed by structural relaxations of the He-filled and empty sI/sII unit cells, yielding lattice and compressibility parameters comparable to experimental and theoretical values available for those hydrates. According to the results obtained, the He enclathration in the sI/sII unit cells is a stabilizing process, and both He@sI and He@sII clathrates, considering single cage occupancy, are predicted to be stable whatever the XDM or D4 dispersion correction applied. Our results further reveal that despite the weak underlying interactions the He encapsulation has a rather notable effect on both lattice parameters and energetics, with the He@sII being the most energetically favorable in accord with recent experimental observations.

Assessment of DFT approaches in noble gas clathrate-like clusters: stability and thermodynamics

We have assessed the performance and accuracy of different wavefunction-based electronic structure methods, such as DFMP2 and domain-based local pair-natural orbital (DLPNO-CCSD(T)), as well as a variety of density functional theory (DFT) approaches on He@(H₂O)_N cage systems. We have selected representative clathrate-like structures corresponding to the building blocks present in each of the sI, sII and sH natural gas clathrate hydrates, and we have carefully studied the interaction between a He atom with each of their individual cages. We reported well-converged DFMP2 and DLPNO-CCSD(T) reference data, together with interaction and cohesive energies of four different density functionals (two GGA, revPBE and PW86PBE, and two hybrids, B3LYP and PBE0), including diverse dispersion correction schemes (D3(0), D3(BJ), D4 and XDM) for both He-filled and empty clathrate-like cages. After the analysis of the results, we came to the conclusion that the PW86PBE functional, with both XDM and D4 corrections, and the PBE0-D4 functional present reasonably adequate approaches to describe the guest-host noncovalent interactions that take place in such He clathrate hydrates. Taking into account that the He@sII is the only helium clathrate that scientists have been able to synthesize recently, we have performed a thermodynamic study on the individual 5¹² and 5¹²6⁴ cages present in the sII crystal. We determined the change in enthalpy, ΔH, and in Gibbs free energy, ΔG, at various temperatures and pressures, and we found out that in the range of experimental conditions the reactions associated with the encapsulation of the He atom inside the cages are exothermic and spontaneous. Finally, we highlighted the importance of an accurate description of the interaction in He@water mixtures, as a crucial component in construction of reliable data-driven models.

Clathrate Adhesion Induced by Quasi-Liquid Layer

The adhesive force of clathrates to surfaces is a century-old problem of pipeline blockage for the energy industry. Here, we provide new physical insight into the origin of this force by accounting for the existence of a quasi-liquid layer (QLL) on clathrate surfaces. To gain this insight, we measure the adhesive force between a tetrahydrofuran clathrate and a solid sphere. We detect a strong adhesion, which originates from a capillary bridge that is formed from a nanometer-thick QLL on the clathrate surface. The curvature of this capillary bridge is nanoscaled, causes a large negative Laplace pressure, and leads to a strong capillary attraction. The microscopic capillary bridge expands and consolidates over time. This dynamic behavior explains the time-dependent increase of measured capillary forces. The adhesive force decreases greatly upon increasing the roughness and the hydrophobicity of the sphere, which founds the fundamental basics for reducing clathrate adhesion by using surface coating.