Mechanical properties of amorphous CO2 hydrates: insights from molecular simulations

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

Nanoscale Prediction of the Thermal, Mechanical, and Transport Properties of Hydrated Clay on 106- and 1015-Fold Larger Length and Time Scales

Coupled thermal, hydraulic, mechanical, and chemical (THMC) processes, such as desiccation-driven cracking or chemically driven fluid flow, significantly impact the performance of composite materials formed by fluid-mediated nanoparticle assembly, including energy storage materials, ordinary Portland cement, bioinorganic nanocomposites, liquid crystals, and engineered clay barriers used in the isolation of hazardous wastes. These couplings are particularly important in the isolation of high-level radioactive waste (HLRW), where heat generated by radioactive decay can drive the temperature up to at least 373 K in the engineered barrier. Here, we use large-scale all-atom molecular dynamics simulations of hydrated smectite clay nanoparticle assemblages to predict the fundamental THMC properties of hydrated compacted clay over a wide range of temperatures (up to 373 K) and dry densities relevant to HLRW management. Equilibrium simulations of clay-water mixtures at different hydration levels are analyzed to quantify material properties, including thermal conductivity, heat capacity, thermal expansion, suction, water and ion self-diffusivity, and hydraulic conductivity. Predictions are validated against experimental results for the properties of compacted bentonite clay. Our results demonstrate the feasibility of using atomistic-level simulations of assemblages of clay nanoparticles on scales of tens of nanometers and nanoseconds to infer the properties of compacted bentonite on scales of centimeters and days, a direct upscaling over 6 orders of magnitude in space and 15 orders of magnitude in time.

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.

Conformer Selection Upon Dilution with Water: The Fascinating Case of Liquid Ethylene Glycol Studied via Molecular Dynamics Simulations

The aqueous solution of ethylene glycol (EG) is a binary liquid mixture that displays rich conformational and structural behaviour, which has not yet been adequately explored through atomistic molecular dynamics simulations. Herein, employing an accurate force field for EG, several physical properties of this solution are calculated to be in quantitative agreement with experimental data. While 79 % of molecules in neat liquid EG exist with their central OCCO dihedral in the gauche state, this fraction increases to 89 % in the dilute aqueous solution, largely in response to the increase in the static dielectric constant of the solution from that of neat liquid EG. The increase in gauche conformers increases the mean dipole moment of EG molecules in the solution which is additionally contributed by specific conformational states of the two terminal HOCC dihedral angles.

Performance-Based Screening of Porous Materials for Carbon Capture

Computational screening methods have changed the way new materials and processes are discovered and designed. For adsorption-based gas separations and carbon capture, recent efforts have been directed toward the development of multiscale and performance-based screening workflows where we can go from the atomistic structure of an adsorbent to its equilibrium and transport properties at different scales, and eventually to its separation performance at the process level. The objective of this work is to review the current status of this new approach, discuss its potential and impact on the field of materials screening, and highlight the challenges that limit its application. We compile and introduce all the elements required for the development, implementation, and operation of multiscale workflows, hence providing a useful practical guide and a comprehensive source of reference to the scientific communities who work in this area. Our review includes information about available materials databases, state-of-the-art molecular simulation and process modeling tools, and a complete catalogue of data and parameters that are required at each stage of the multiscale screening. We thoroughly discuss the challenges associated with data availability, consistency of the models, and reproducibility of the data and, finally, propose new directions for the future of the field.

Molecular Origins of Deformation in Amorphous Methane Hydrates

Water and methane can stay together under low temperature and high pressure in the forms of liquid solutions and crystalline solids. From liquid and gaseous states to crystalline solids or the contrary processes, amorphous methane hydrates can occur in these evolution scenarios. Herein, mechanical properties of amorphous methane hydrates are explored for the first time to bridge the gap between mechanical responses of monocrystalline and polycrystalline methane hydrates. Our results demonstrate that mechanical properties of amorphous methane hydrates are strongly governed by our original proposed order parameter, namely, normalized hydrogen-bond directional order parameter. Followed by this important achievement, a multistep deformation mechanism core is proposed to explain mechanical properties of amorphous methane hydrates. Through an extensive detailed analysis of amorphous methane hydrates, our simulation results not only greatly enlarge our fundamental understanding for mechanical responses of amorphous methane hydrates in geological systems but also offer a fresh perspective in structure-property topics of solid materials in future science and technology.

A novel hybrid method for the calculation of methane hydrate-water interfacial tension along the three-phase (hydrate-liquid water-vapor) equilibrium line

We use a novel hybrid method to explore the temperature dependence of the solid-liquid interfacial tension of a system that consists of solid methane hydrate and liquid water. The calculated values along the three-phase (hydrate-liquid water-vapor) equilibrium line are obtained through the combination of available experimental measurements and computational results that are based on approaches at the atomistic scale, including molecular dynamics and Monte Carlo. An extensive comparison with available experimental and computational studies is performed, and a critical assessment and re-evaluation of previously reported data is presented.

Mechanical creep instability of nanocrystalline methane hydrates

Mechanical creep behaviors of natural gas hydrates are of importance for understanding the mechanical instability of gas hydrate-bearing sediments on Earth. Limited by the experimental challenges, intrinsic creep mechanisms of nanocrystalline methane hydrates remain largely unknown yet at the molecular scale. Herein, using large-scale molecular dynamics simulations, mechanical creep behaviors of nanocrystalline methane hydrates are investigated. It is revealed that mechanical creep responses are greatly dictated by internal microstructures of crystalline grain size and external conditions of temperature and static stress. Interestingly, a long steady-state creep is observed in nanocrystalline methane hydrates, which can be described by a modified constitutive Bird-Dorn-Mukherjee model. Microstructural analysis shows that deformations of crystalline grains, grain boundary diffusion and grain boundary sliding collectively govern the mechanical creep behaviors of nanocrystalline methane hydrates. Furthermore, structural transformation also appears to be important in their mechanical creep behaviors. This study provides new insights into understanding the mechanical creep scenarios of gas hydrates.

Liquid dibromomethane under pressure: a computational study

Molecular dynamics simulations have been performed on liquid dibromomethane at thermodynamic states corresponding to temperature in the range 268-328 K and pressure varying from 1 bar to 3000 bar. The interaction model is a simple effective two-body pair potential with atom-atom Coulomb and Lennard-Jones interactions and molecules are rigid. Thermodynamic properties have been studied, including the isobaric thermal expansion coefficient, the isothermal compressibility, the heat capacities and the speed of sound. The simulation results exhibit a crossing of the isotherms of the isobaric thermal expansion coefficient at about 800 bar in very good agreement with the prediction of an isothermal fluctuation equation of state predicting such a crossing in the pressure range 650-900 bar, though experimental results up to 1000 bar do not find any crossing.

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.

Molecular dynamics simulations of interfacial properties of the CO2-water and CO2-CH4-water systems

Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO₂ and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO₂-water interface has a lower interfacial tension than the CH₄-water interface. A greater number of CO₂ molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO₂ and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO₂. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO₂-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO₂ capture and sequestration processes are discussed.

Molecular dynamics simulation of sI methane hydrate under compression and tension

Molecular dynamics (MD) analysis of methane hydrate is important for the application of methane hydrate technology. This study investigated the microstructure changes of sI methane hydrate and the laws of stress–strain evolution under the condition of compression and tension by using MD simulation. This study further explored the mechanical property and stability of sI methane hydrate under different stress states. Results showed that tensile and compressive failures produced an obvious size effect under a certain condition. At low temperature and high pressure, most of the clathrate hydrate maintained a stable structure in the tensile fracture process, during which only a small amount of unstable methane broke the structure, thereby, presenting a free-motion state. The methane hydrate cracked when the system reached the maximum stress in the loading process, in which the maximum compressive stress is larger than the tensile stress under the same experimental condition. This study provides a basis for understanding the microscopic stress characteristics of methane hydrate.

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.

The effect of surfactants on hydrate particle agglomeration in liquid hydrocarbon continuous systems: a molecular dynamics simulation study

Anti-agglomerants (AAs), both natural and commercial, are currently being considered for gas hydrate risk management of petroleum pipelines in offshore operations. However, the molecular mechanisms of the interaction between the AAs and gas hydrate surfaces and the prevention of hydrate agglomeration remain critical and complex questions that need to be addressed to advance this technology. Here, we use molecular dynamics (MD) simulations to investigate the effect of model surfactant molecules (polynuclear aromatic carboxylic acids) on the agglomeration behaviour of gas hydrate particles and disruption of the capillary liquid bridge between hydrate particles. The results show that the anti-agglomeration pathway can be divided into two processes: the spontaneous adsorption effect of surfactant molecules onto the hydrate surface and the weakening effect of the intensity of the liquid bridge between attracted hydrate particles. The MD simulation results also indicate that the anti-agglomeration effectiveness of surfactants is determined by the intrinsic nature of their molecular functional groups. Additionally, we find that surfactant molecules can affect hydrate growth, which decreases hydrate particle size and correspondingly lower the risk of hydrate agglomeration. This study provides molecular-level insights into the anti-agglomeration mechanism of surfactant molecules, which can aid in the ultimate application of natural or commercial AAs with optimal anti-agglomeration properties.

Schematic of anti-agglomeration effect of surfactants promoting gas hydrate particle dispersion.

Size and shape effects on the thermodynamic properties of nanoscale volumes of water

Small systems are known to deviate from the classical thermodynamic description, among other things due to their large surface area to volume ratio compared to corresponding big systems. As a consequence, extensive thermodynamic properties are no longer proportional to the volume, but are instead higher order functions of size and shape. We investigate such functions for second moments of probability distributions of fluctuating properties in the grand-canonical ensemble, focusing specifically on the volume and surface terms of Hadwiger's theorem, explained in Klain, Mathematika, 1995, 42, 329-339. We resolve the shape dependence of the surface term and show, using Hill's nanothermodynamics [Hill, J. Chem. Phys., 1962, 36, 3182], that the surface satisfies the thermodynamics of a flat surface as described by Gibbs [Gibbs, The Scientific Papers of J. Willard Gibbs, Volume 1, Thermodynamics, Ox Bow Press, Woodbridge, Connecticut, 1993]. The Small System Method (SSM), first derived by Schnell et al. [Schnell et al., J. Phys. Chem. B, 2011, 115, 10911], is extended and used to analyze simulation data on small systems of water. We simulate water as an example to illustrate the method, using TIP4P/2005 and other models, and compute the isothermal compressibility and thermodynamic factor. We are able to retrieve the experimental value of the bulk phase compressibility within 2%, and show that the compressibility of nanosized volumes increases by up to a factor of two as the number of molecules in the volume decreases. The value for a tetrahedron, cube, sphere, polygon, etc. can be predicted from the same scaling law, as long as second order effects (nook and corner effects) are negligible. Lastly, we propose a general formula for finite reservoir correction to fluctuations in subvolumes.

Revealing at the molecular level the role of the surfactant in the enhancement of the thermal properties of the gold nanofluid system used for concentrating solar power

A surfactant acts as a kind of net around gold that enhances at the molecular level the gold nanofluid thermal properties.

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

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

Microscopic Origin of Strain Hardening in Methane Hydrate

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

Ab initio modelling of methane hydrate thermophysical properties

Using density functional theory, the second-order elastic constants, heat capacity, compressibility, and thermal expansion coefficient of methane hydrate were calculated.

Mechanical instability of monocrystalline and polycrystalline methane hydrates

Despite observations of massive methane release and geohazards associated with gas hydrate instability in nature, as well as ductile flow accompanying hydrate dissociation in artificial polycrystalline methane hydrates in the laboratory, the destabilising mechanisms of gas hydrates under deformation and their grain-boundary structures have not yet been elucidated at the molecular level. Here we report direct molecular dynamics simulations of the material instability of monocrystalline and polycrystalline methane hydrates under mechanical loading. The results show dislocation-free brittle failure in monocrystalline hydrates and an unexpected crossover from strengthening to weakening in polycrystals. Upon uniaxial depressurisation, strain-induced hydrate dissociation accompanied by grain-boundary decohesion and sliding destabilises the polycrystals. In contrast, upon compression, appreciable solid-state structural transformation dominates the response. These findings provide molecular insight not only into the metastable structures of grain boundaries, but also into unusual ductile flow with hydrate dissociation as observed during macroscopic compression experiments.

Sediment-hosted gas hydrates may release vast quantities of methane upon failure, but destabilizing mechanisms at the molecular level are poorly understood. Here, the authors study the deformation using simulations and find that failure differs between single crystals and polycrystalline hydrates.