Mesomorphology of clathrate hydrates from molecular ordering.

Clathrate hydrates are crystals formed by guest molecules that stabilize cages of hydrogen-bonded water molecules. Whereas thermodynamic equilibrium is well described via the van der Waals and Platteeuw approach, the increasing concerns with global warming and energy transition require extending the knowledge to non-equilibrium conditions in multiphase, sheared systems, in a multiscale framework. Potential macro-applications concern the storage of carbon dioxide in the form of clathrates, and the reduction of hydrate inhibition additives currently required in hydrocarbon production. We evidence porous mesomorphologies as key to bridging the molecular scales to macro-applications of low solubility guests. We discuss the coupling of molecular ordering with the mesoscales, including (i) the emergence of porous patterns as a combined factor from the walk over the free energy landscape and 3D competitive nucleation and growth and (ii) the role of molecular attachment rates in crystallization-diffusion models that allow predicting the timescale of pore sealing. This is a perspective study that discusses the use of discrete models (molecular dynamics) to build continuum models (phase field models, crystallization laws, and transport phenomena) to predict multiscale manifestations at a feasible computational cost. Several advances in correlated fields (ice, polymers, alloys, and nanoparticles) are discussed in the scenario of clathrate hydrates, as well as the challenges and necessary developments to push the field forward.

Novel approach for designing order parameters of clathrate hydrate structures by graph neural network.

Clathrate hydrates continue to be the focus of active research efforts due to their use in energy resources, transportation, and storage-related applications. Therefore, it is crucial to define their essential characteristics from a molecular standpoint. Understanding molecular structure in particular is crucial because it aids in understanding the mechanisms that lead to the formation or dissociation of clathrate hydrates. In the past, a wide variety of order parameters have been employed to classify and evaluate hydrate structures. An alternative approach to inventing bespoke order parameters is to apply machine learning techniques to automatically generate effective order parameters. In earlier work, we suggested a method for automatically designing novel parameters for ice and liquid water structures with Graph Neural Networks (GNNs). In this work, we use a GNN to implement our method, which can independently produce feature representations of the molecular structures. By using the TeaNet-type model in our method, it is possible to directly learn the molecular geometry and topology. This enables us to build novel parameters without prior knowledge of suitable order parameters for the structure type, discover structural differences, and classify molecular structures with high accuracy. We use this approach to classify the structures of clathrate hydrate structures: sI, sII, and sH. This innovative approach provides an appealing and highly accurate replacement for the traditional order parameters. Furthermore, our method makes clear the process of automatically designing a universal parameter for liquid water, ice, and clathrate hydrate to analyze their structures and phases.

Methane hydrate phase equilibrium considering dissolved methane concentrations and interfacial geometries from molecular simulations.

Natural gas hydrates, mainly existing in permafrost and on the seabed, are expected to be a new energy source with great potential. The exploitation technology of natural gas hydrates is one of the main focuses of hydrate-related studies. In this study, a large-size liquid aqueous solution wrapping a methane hydrate system was established and molecular dynamics simulations were used to investigate the phase equilibrium conditions of methane hydrate at different methane concentrations and interfacial geometries. It is found that the methane concentration of a solution significantly affects the phase equilibrium of methane hydrates. Different methane concentrations at the same temperature and pressure can lead to hydrate formation or decomposition. At the same temperature and pressure, in a system reaching equilibrium, the size of spherical hydrate clusters is coupled to the solution concentration, which is proportional to the Laplace pressure at the solid-liquid interface. Lower solution concentrations reduce the phase equilibrium temperature of methane hydrates at the same pressure; as the concentration increases, the phase equilibrium temperature gradually approaches the actual phase equilibrium temperature. In addition, the interfacial geometry of hydrates affects the thermodynamic stability of hydrates. The spherical hydrate particles have the highest stability for the same volume. Through this study, we provide a stronger foundation to understand the principles driving hydrate formation/dissociation relevant to the exploitation of methane hydrates.

Sixty Years of the van der Waals and Platteeuw Model for Clathrate Hydrates-A Critical Review from Its Statistical Thermodynamic Basis to Its Extensions and Applications.

In light of the 60-year anniversary of the publishing of "Clathrate Solutions" by van der Waals and Platteeuw in 2019, we present a critical review of the famed solid solution model first disclosed in 1959. First, we lay out the groundwork in the 1950s aimed at the development of a phenomenological approach to clathrate modeling. Then we review the statistical thermodynamics fundamentals of the model, considering van der Waals and Platteeuw's earlier works, to obtain a consistent interpretation of the model. We turn our focus to clathrate hydrates and discuss the major contributions that led to the current state-of-the-art of gas hydrate thermodynamic modeling. Finally, we present some of the areas in clathrate thermodynamics that we foresee as the new frontiers in this subject. We expect this review to help newcomers to clathrate science in elucidating some subtle aspects of the model and to intrigue clathrate experts with a fresh look on this well-established solid solution model.

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.

Mechanism for H2 diffusion in sII hydrates by molecular dynamics simulations.

Among the many different types of molecules that form clathrate hydrates, H2 is unique as it can easily diffuse into and out of clathrate cages, a process that involves the physical-chemical interactions between guest (H2) and host (water) molecules, and is unlike any other molecular system. The dynamic and nano-scale process of H2 diffusion into binary structure II hydrates, where the large cages are occupied by larger molecules, was studied using molecular dynamics simulation. As the H2 molecules diffused from one cage to another, two types of diffusion processes were observed: (i) when moving between a pair of large cages, the H2 molecules pass through the central part of the hexagonal rings; (ii) however, when the H2 molecules move from a large cage to a small one, it requires one of the pentagonal rings to partially break, as this allows the H2 molecule to pass through the widened space. While the diffusion of H2 molecules between large cages was found to occur more frequently, the presence of SF6 molecules in the large cages was found to inhibit diffusion. Therefore, in order to attain higher H2 storage capacities in binary hydrates, it is suggested that there is an optimal number of large cages that should be occupied by SF6 molecules.

Cage occupancies, lattice constants, and guest chemical potentials for structure II hydrogen clathrate hydrate from Gibbs ensemble Monte Carlo simulations.

In this paper, equilibrium properties of structure II hydrates of hydrogen were determined from Monte Carlo simulations in the isothermal-isobaric Gibbs ensemble. Water and hydrogen molecules are described by the TIP4P/Ice and Silvera-Goldman models, respectively. The use of the Gibbs ensemble has many key advantages for the simulation of hydrates. By the separation of hydrogen vapor and hydrate phases into their own domains, coupled with transfer moves of hydrogen molecules between domains, cage occupancies were determined. Furthermore, the choice of this ensemble also allows equilibrium lattice constants and guest molecule chemical potentials to be straightforwardly estimated. Results for hydrogen mass fractions indicate reasonable agreement with prior simulation data and theoretical models, while detailed analysis of cage occupancy distributions and neighboring cage pair occupancy combinations gives valuable insight into the behavior of this hydrate at the inter-cage scale. These results will aid in the construction of theoretical models, for which knowledge of the occupancy of neighboring cages is of great importance. In support of previous experimental and theoretical works, we also find evidence of double occupancy of a few small cages inside of the hydrate stability zone, albeit at very high pressures; approximately 0.1% of small cages are doubly occupied at 300 MPa, for temperatures of 225 K and 250 K.

Analysis of three-phase equilibrium conditions for methane hydrate by isometric-isothermal molecular dynamics simulations.

To develop prediction methods of three-phase equilibrium (coexistence) conditions of methane hydrate by molecular simulations, we examined the use of NVT (isometric-isothermal) molecular dynamics (MD) simulations. NVT MD simulations of coexisting solid hydrate, liquid water, and vapor methane phases were performed at four different temperatures, namely, 285, 290, 295, and 300 K. NVT simulations do not require complex pressure control schemes in multi-phase systems, and the growth or dissociation of the hydrate phase can lead to significant pressure changes in the approach toward equilibrium conditions. We found that the calculated equilibrium pressures tended to be higher than those reported by previous NPT (isobaric-isothermal) simulation studies using the same water model. The deviations of equilibrium conditions from previous simulation studies are mainly attributable to the employed calculation methods of pressure and Lennard-Jones interactions. We monitored the pressure in the methane phase, far from the interfaces with other phases, and confirmed that it was higher than the total pressure of the system calculated by previous studies. This fact clearly highlights the difficulties associated with the pressure calculation and control for multi-phase systems. The treatment of Lennard-Jones interactions without tail corrections in MD simulations also contributes to the overestimation of equilibrium pressure. Although improvements are still required to obtain accurate equilibrium conditions, NVT MD simulations exhibit potential for the prediction of equilibrium conditions of multi-phase systems.

Micromechanical cohesion force between gas hydrate particles measured under high pressure and low temperature conditions.

To prevent hydrate plugging conditions in the transportation of oil/gas in multiphase flowlines, one of the key processes to control is the agglomeration/deposition of hydrate particles, which are determined by the cohesive/adhesive forces. Previous studies reporting measurements of the cohesive/adhesive force between hydrate particles used cyclopentane hydrate particles in a low-pressure micromechanical force apparatus. In this study, we report the cohesive forces of particles measured in a new high-pressure micromechanical force (MMF) apparatus for ice particles, mixed (methane/ethane, 74.7:25.3) hydrate particles (Structure II), and carbon dioxide hydrate particles (Structure I). The cohesive forces are measured as a function of the contact time, contact force, temperature, and pressure, and determined from pull-off measurements. For the measurements performed of the gas hydrate particles in the gas phase, the determined cohesive force is about 30-35 mN/m, about 8 times higher than the cohesive force of CyC5 hydrates in the liquid CyC5, which is about 4.3 mN/m. We show from our results that the hydrate structure (sI with CO2 hydrates and sII with CH4/C2H6 hydrates) has no influence on the cohesive force. These results are important in the deposition of a gas-dominated system, where the hydrate particles formed in the liquid phase can then stick to the hydrate deposited in the wall exposed to the gas phase.

Nucleation rate analysis of methane hydrate from molecular dynamics simulations.

Clathrate hydrates are solid crystalline structures most commonly formed from solutions that have nucleated to form a mixed solid composed of water and gas. Understanding the mechanism of clathrate hydrate nucleation is essential to grasp the fundamental chemistry of these complex structures and their applications. Molecular dynamics (MD) simulation is an ideal method to study nucleation at the molecular level because the size of the critical nucleus and formation rate occur on the nano scale. Various analysis methods for nucleation have been developed through MD to analyze nucleation. In particular, the mean first-passage time (MFPT) and survival probability (SP) methods have proven to be effective in procuring the nucleation rate and critical nucleus size for monatomic systems. This study assesses the MFPT and SP methods, previously used for monatomic systems, when applied to analyzing clathrate hydrate nucleation. Because clathrate hydrate nucleation is relatively difficult to observe in MD simulations (due to its high free energy barrier), these methods have yet to be applied to clathrate hydrate systems. In this study, we have analyzed the nucleation rate and critical nucleus size of methane hydrate using MFPT and SP methods from data generated by MD simulations at 255 K and 50 MPa. MFPT was modified for clathrate hydrate from the original version by adding the maximum likelihood estimate and growth effect term. The nucleation rates calculated by MFPT and SP methods are within 5%, and the critical nucleus size estimated by the MFPT method was 50% higher, than values obtained through other more rigorous but computationally expensive estimates. These methods can also be extended to the analysis of other clathrate hydrates.

A molecular dynamics study of guest-host hydrogen bonding in alcohol clathrate hydrates.

Clathrate hydrates are typically stabilized by suitably sized hydrophobic guest molecules. However, it has been experimentally reported that isomers of amyl-alcohol C5H11OH can be enclosed into the 5(12)6(4) cages in structure II (sII) clathrate hydrates, even though the effective radii of the molecules are larger than the van der Waals radii of the cages. To reveal the mechanism of the anomalous enclathration of hydrophilic molecules, we performed ab initio and classical molecular dynamics simulations (MD) and analyzed the structure and dynamics of a guest-host hydrogen bond for sII 3-methyl-1-butanol and structure H (sH) 2-methyl-2-butanol clathrate hydrates. The simulations clearly showed the formation of guest-host hydrogen bonds and the incorporation of the O-H group of 3-methyl-1-butanol guest molecules into the framework of the sII 5(12)6(4) cages, with the remaining hydrophobic part of the amyl-alcohol molecule well accommodated into the cages. The calculated vibrational spectra of alcohol O-H bonds showed large frequency shifts due to the strong guest-host hydrogen bonding. The 2-methyl-2-butanol guests form strong hydrogen bonds with the cage water molecules in the sH clathrate, but are not incorporated into the water framework. By comparing the structures of the alcohols in the hydrate phases, the effect of the location of O-H groups in the butyl chain of the guest molecules on the crystalline structure of the clathrate hydrates is indicated.

Biophysical changes induced by xenon on phospholipid bilayers.

Structural and dynamic changes in cell membrane properties induced by xenon, a volatile anesthetic molecule, may affect the function of membrane-mediated proteins, providing a hypothesis for the mechanism of general anesthetic action. Here, we use molecular dynamics simulation and differential scanning calorimetry to examine the biophysical and thermodynamic effects of xenon on model lipid membranes. Our results indicate that xenon atoms preferentially localize in the hydrophobic core of the lipid bilayer, inducing substantial increases in the area per lipid and bilayer thickness. Xenon depresses the membrane gel-liquid crystalline phase transition temperature, increasing membrane fluidity and lipid head group spacing, while inducing net local ordering effects in a small region of the lipid carbon tails and modulating the bilayer lateral pressure profile. Our results are consistent with a role for nonspecific, lipid bilayer-mediated mechanisms in producing xenon's general anesthetic action.

Quantitative measurement and mechanisms for CH4 production from hydrates with the injection of liquid CO2.

The recovery of gas from natural gas hydrates under the permafrost and in oceanic sediments is of particular interest in energy and environmental fields because of the attractive process to release methane gas through the injection of CO2. The sequestration of CO2, a notorious greenhouse gas, in hydrates has the potential to be used in enhanced gas recovery techniques, while simultaneously releasing CH4 locked within the gas bearing hydrates. In this study, we present quantitative experiments to investigate results of possible CH4-CO2 exchange kinetics from injection of liquid CO2 through CH4 hydrates. The experiments performed use CH4 hydrate formed from ice particles (75-90 or 125-150 microns in diameter) at approximately 10.34 MPa and 263 K. In order to reduce unexpected errors, nearly full conversion (>95%) of ice particles to hydrates is achieved. Liquid CO2 is injected into the pressure cell to sweep the residual CH4 atmosphere, ensuring no free CH4 is left in the gas phase. After soaking the hydrate for several hours, CH4 is produced from the hydrates by injecting liquid CO2. The final composition and analysis of the produced CH4 is measured by using in-line gas chromatography. We also measure the CH4 moles after hydrate dissociation to confirm the closure of the total mass balance of the experiment. From these data, we infer the mechanism for CH4 production, identify the penetration depth of the dissociation/exchange on the hydrate particles, and propose physical models describing the mechanism for CH4 production. These experiments are essential in the quantification of the production of CH4 from CH4 hydrates with the injection of CO2.

Adhesion force interactions between cyclopentane hydrate and physically and chemically modified surfaces.

Interfacial interactions between liquid-solid and solid-solid phases/surfaces are of fundamental importance to the formation of hydrate deposits in oil and gas pipelines. This work establishes the effect of five categories of physical and chemical modification to steel on clathrate hydrate adhesive force: oleamide, graphite, citric acid ester, nonanedithiol, and Rain-X anti-wetting agent. Hydrate adhesive forces were measured using a micromechanical force apparatus, under both dry and water-wet surface conditions. The results show that the graphite coating reduced hydrate-steel adhesion force by 79%, due to an increase in the water wetting angle from 42 ± 8° to 154 ± 7°. Two chemical surface coatings (nonanedithiol and the citric acid ester) induced rapid hydrate growth in the hydrate particles; nonanedithiol increased hydrate adhesive force by 49% from the baseline, while the citric acid ester coating reduced hydrate adhesion force by 98%. This result suggests that crystal growth may enable a strong adhesive pathway between hydrate and other crystalline structures, however this effect may be negated in cases where water-hydrocarbon interfacial tension is minimised. When a liquid water droplet was placed on the modified steel surfaces, the graphite and citric acid ester became less effective at reducing adhesive force. In pipelines containing a free water phase wetting the steel surface, chemical or physical surface modifications alone may be insufficient to eliminate hydrate deposition risk. In further tests, the citric acid ester reduced hydrate cohesive forces by 50%, suggesting mild activity as a hybrid anti-agglomerant suppressing both hydrate deposition and particle agglomeration. These results demonstrate a new capability to develop polyfunctional surfactants, which simultaneously limit the capability for hydrate particles to aggregate and deposit on the pipeline wall.

Two-component order parameter for quantifying clathrate hydrate nucleation and growth.

Methane clathrate hydrate nucleation and growth is investigated via analysis of molecular dynamics simulations using a new order parameter. This order parameter (OP), named the Mutually Coordinated Guest (MCG) OP, quantifies the appearance and connectivity of molecular clusters composed of guests separated by water clusters. It is the first two-component OP used for quantifying hydrate nucleation and growth. The algorithm for calculating the MCG OP is described in detail. Its physical motivation and advantages compared to existing methods are discussed.

Observation of interstitial molecular hydrogen in clathrate hydrates.

The current knowledge and description of guest molecules within clathrate hydrates only accounts for occupancy within regular polyhedral water cages. Experimental measurements and simulations, examining the tert-butylamine + H2 + H2O hydrate system, now suggest that H2 can also be incorporated within hydrate crystal structures by occupying interstitial sites, that is, locations other than the interior of regular polyhedral water cages. Specifically, H2 is found within the shared heptagonal faces of the large (4(3)5(9)6(2)7(3)) cage and in cavities formed from the disruption of smaller (4(4)5(4)) water cages. The ability of H2 to occupy these interstitial sites and fluctuate position in the crystal lattice demonstrates the dynamic behavior of H2 in solids and reveals new insight into guest-guest and guest-host interactions in clathrate hydrates, with potential implications in increasing overall energy storage properties.

Surfactant adsorption and interfacial tension investigations on cyclopentane hydrate.

Gas hydrates represent an unconventional methane resource and a production/safety risk to traditional oil and gas flowlines. In both systems, hydrate may share interfaces with both aqueous and hydrocarbon fluids. To accurately model macroscopic properties, such as relative permeability in unconventional systems or dispersion viscosity in traditional systems, knowledge of hydrate interfacial properties is required. This work presents hydrate cohesive force results measured on a micromechanical force apparatus, and complementary water-hydrocarbon interfacial tension data. By combining a revised cohesive force model with experimental data, two interfacial properties of cyclopentane hydrate were estimated: hydrate-water and hydrate-cyclopentane interfacial tension values at 0.32 ± 0.05 mN/m and 47 ± 5 mN/m, respectively. These fundamental physiochemical properties have not been estimated or measured for cyclopentane hydrate to date. The addition of surfactants in the cyclopentane phase significantly reduced the cyclopentane hydrate cohesive force; we hypothesize this behavior to be the result of surfactant adsorption on the hydrate-oil interface. Surface excess quantities were estimated for hydrate-oil and water-oil interfaces using four carboxylic and sulfonic acids. The results suggest the density of adsorbed surfactant may be 2× larger for the hydrate-oil interface than the water-oil interface. Additionally, hydrate-oil interfacial tension was observed to begin decreasing from the baseline value at significantly lower surfactant concentrations (1-3 orders of magnitude) than those for the water-oil interfacial tension.

Adhesion force between cyclopentane hydrate and mineral surfaces.

Clathrate hydrate adhesion forces play a critical role in describing aggregation and deposition behavior in conventional energy production and transportation. This manuscript uses a unique micromechanical force apparatus to measure the adhesion force between cyclopentane hydrate and heterogeneous quartz and calcite substrates. The latter substrates represent models for coproduced sand and scale often present during conventional energy production and transportation. Micromechanical adhesion force data indicate that clathrate hydrate adhesive forces are 5-10× larger for calcite and quartz minerals than stainless steel. Adhesive forces further increased by 3-15× when increasing surface contact time from 10 to 30 s. In some cases, liquid water from within the hydrate shell contacted the mineral surface and rapidly converted to clathrate hydrate. Further measurements on mineral surfaces with physical control of surface roughness showed a nonlinear dependence of water wetting angle on surface roughness. Existing adhesive force theory correctly predicted the dependence of clathrate hydrate adhesive force on calcite wettability, but did not accurately capture the dependence on quartz wettability. This comparison suggests that the substrate surface may not be inert, and may contribute positively to the strength of the capillary bridge formed between hydrate particles and solid surfaces.

Water proton configurations in structures I, II, and H clathrate hydrate unit cells.

Position and orientation of water protons need to be specified when the molecular simulation studies are performed for clathrate hydrates. Positions of oxygen atoms in water are experimentally determined by X-ray diffraction analysis of clathrate hydrate structures, but positions of water hydrogen atoms in the lattice are disordered. This study reports a determination of the water proton coordinates in unit cell of structure I (sI), II (sII), and H (sH) clathrate hydrates that satisfy the ice rules, have the lowest potential energy configuration for the protons, and give a net zero dipole moment. Possible proton coordinates in the unit cell were chosen by analyzing the symmetry of protons on the hexagonal or pentagonal faces in the hydrate cages and generating all possible proton distributions which satisfy the ice rules. We found that in the sI and sII unit cells, proton distributions with small net dipole moments have fairly narrow potential energy spreads of about 1 kJ∕mol. The total Coulomb potential on a test unit charge placed in the cage center for the minimum energy∕minimum dipole unit cell configurations was calculated. In the sI small cages, the Coulomb potential energy spread in each class of cage is less than 0.1 kJ∕mol, while the potential energy spread increases to values up to 6 kJ∕mol in sH and 15 kJ∕mol in the sII cages. The guest environments inside the cages can therefore be substantially different in the sII case. Cartesian coordinates for oxygen and hydrogen atoms in the sI, sII, and sH unit cells are reported for reference.