Ab initioDFT study of structural and mechanical properties of methane and carbon dioxide hydrates

Ab Initio Studies of Structural and Mechanical Properties of NH3, NO, and N2O Hydrates

The investigation on the mechanical properties of clathrate hydrate is closely related to the exploitation of hydrates and gas transportation. In this article, the structural and mechanical properties of some nitride gas hydrates were studied using DFT calculations. First, the equilibrium lattice structure is obtained by geometric structure optimization; then, the complete second-order elastic constant is determined by energy-strain analysis, and the polycrystalline elasticity is predicted. It is found that the NH₃, N₂O, and NO hydrates all have high elastic isotropy but are different in shear characteristics. This work may lay a theoretical foundation for studying the structural evolution of clathrate hydrates under the mechanical field.

Atomistic-geometric simulations to investigate the mechanical stability of monocrystalline sI methane hydrates under pressure

Gas hydrate mechanical stability under pressure is critically important in energy supply, global warming, and carbon-neutral technologies. The stability of these polyhedral guest-host crystals under increasing pressure is affected by host cage type and face connectivity as well as guest gas occupancy. The geometry-imposed cage connectivity generates crystal lattices that include inclusion-matrix material composite structures. In this paper, we integrate Density Functional Theory simulations with a polyhedral-inspired composite material model that quantifies stability limits, failure modes, and the impact of the type of cage occupancy. DFT reveals the existence of two failure mechanisms under increasing pressure: (i) a multistep lattice breakdown under total occupancy and under only large cage occupancy and (ii) a single-step breakdown under zero occupancy as well as with only small cage occupancy. The DFT-composite model predicts optimal occupancy pathways to generate strength and critical occupancy pathways to promote decomposition.

Modeling of Structure H Carbon Dioxide Clathrate Hydrates: Guest-Lattice Energies, Crystal Structure, and Pressure Dependencies

We performed first-principles computations to investigate the complex interplay of molecular interaction energies in determining the lattice structure and stability of CO₂@sH clathrate hydrates. Density functional theory computations using periodic boundary conditions were employed to characterize energetics and the key structural properties of the sH clathrate crystal under pressure, such as equilibrium lattice volume and bulk modulus. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. Structural relaxations of the fully CO₂-filled and empty sH unit cells yield crystal structure and lattice energies, while their compressibility parameters were derived by including the pressure dependencies. The present quantum chemistry computations suggest anisotropy in the compressibility of the sH clathrate hydrates, with the crystal being less compressible along the a-axis direction than along the c-axis one, in distinction from nearly isotropic sI and sII structures. The detailed results presented here give insight into the complex nature of the underlying guest-host interactions, checking earlier assumptions, providing critical tests, and improving estimates. Such entries may eventually lead to better predictions of thermodynamic properties and formation conditions, with a direct impact on emerging hydrate-based technologies.

First-Principles Elastic and Anisotropic Characteristics of Structure-H Gas Hydrate under Pressure

Evaluating gas hydrates properties contributes valuably to their large-scale management and utilization in fundamental science and applications. Noteworthy, structure-H (sH) gas hydrate lacks a comprehensive characterization of its structural, mechanical, and anisotropic properties. Anisotropic and pressure dependent properties are crucial for gas hydrates’ detection and recovery studies. The objective of this work is the determination of pressure-dependent elastic constants and mechanical properties and the direction-dependent moduli of sH gas hydrates as a function of guest composition. First-principles DFT computations are used to evaluate the mechanical properties, anisotropy, and angular moduli of different sH gas hydrates under pressure. Some elastic constants and moduli increase more significantly with pressure than others. This introduces variations in sH gas hydrate’s incompressibility, elastic and shear resistance, and moduli anisotropy. Young’s modulus of sH gas hydrate is more anisotropic than its shear modulus. The anisotropy of sH gas hydrates is characterized using the unit cell elastic constants, anisotropy factors, and the angular dependent moduli. Structure-properties composition correlations are established as a function of pressure. It is found that compressing filled sH gas hydrates increases their moduli anisotropy. Differences in atomic bonding across a crystal’s planes can be expected in anisotropic structures. Taken together the DFT-based structure–properties–composition relations for sH gas hydrates provide novel and significant material physics results for technological applications.

Exploring CO2 @sI Clathrate Hydrates as CO2 Storage Agents by Computational Density Functional Approaches

The formation of specific clathrate hydrates and their transformation at given thermodynamic conditions depends on the interactions between the guest molecule/s and the host water lattice. Understanding their structural stability is essential to control structure-property relations involved in different technological applications. Thus, the energetic aspects relative to CO₂ @sI clathrate hydrate are investigated through the computation of the underlying interactions, dominated by hydrogen bonds and van der Waals forces, from first-principles electronic structure approaches. The stability of the CO₂ @sI clathrate is evaluated by combining bottom-up and top-down approaches. Guest-free and CO₂ guest-filled aperiodic cages, up to the gradually CO₂ occupation of the entire sI periodic unit cells were considered. Saturation, cohesive and binding energies for the systems are determined by employing a variety of density functionals and their performance is assessed. The dispersion corrections on the non-covalent interactions are found to be important in the stabilization of the CO₂ @sI energies, with the encapsulation of the CO₂ into guest-free/empty cage/lattice being always an energetically favorable process for most of the functionals studied. The PW86PBE functional with XDM or D3(BJ) dispersion corrections predicts a lattice constant in accord to the experimental values available, and simultaneously provides a reliable description for the guest-host interactions in the periodic CO₂ @sI crystal, as well as the energetics of its progressive single cage occupancy process. It has been found that the preferential orientation of the single CO₂ in the large sI crystal cages has a stabilizing effect on the hydrate, concluding that the CO₂ @sI structure is favored either by considering the individual building block cages or the complete sI unit cell crystal. Such benchmark and methodology cross-check studies benefit new data-driven model research by providing high-quality training information, with new insights that indicate the underlying factors governing their structure-driven stability, and triggering further investigations for controlling the stabilization of these promising long-term CO₂ storage materials.

Computational density-functional approaches on finite-size and guest-lattice effects in CO2@sII clathrate hydrate

We performed first-principles computations to investigate guest-host/host-host effects on the encapsulation of the CO₂ molecule in sII clathrate hydrates from finite-size clusters up to periodic 3D crystal lattice systems. Structural and energetic properties were first computed for the individual and first-neighbors clathrate-like sII cages, where highly accurate ab initio quantum chemical methods are available nowadays, allowing in this way the assessment of the density functional (DFT) theoretical approaches employed. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. On this basis, structural relaxations of the CO₂-filled and empty sII unit cells yield lattice and compressibility parameters comparable to experimental and previous theoretical values available for sII hydrates. According to these data, the CO₂ enclathration in the sII clathrate cages is a stabilizing process, either by considering both guest-host and host-host interactions in the complete unit cell or only the guest-water energies for the individual clathrate-like sII cages. CO₂@sII clathrates are predicted to be stable whatever the dispersion correction applied and in the case of single cage occupancy are found to be more stable than the CO₂@sI structures. Our results reveal that DFT approaches could provide a good reasonable description of the underlying interactions, enabling the investigation of formation and transformation processes as a function of temperature and pressure.

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.

Infrared Spectra of Gas Hydrates from First-Principles

The infrared spectra of sII gas hydrates have been computed using density functional theory for the first time, at equilibrium, and under pressure. It is also the first account of a full vibrational analysis (both guest and host vibrations) for gas hydrates with hydrocarbon guest molecules. Five hydrate structures were investigated: empty, propane, isobutane, ethane-methane, and propane-methane sII hydrates. The computed IR spectra are in good agreement with available experimental and theoretical results. The OH stretching frequencies were found to increase, while the H-bond stretching and H₂O libration frequencies decreased with an increase in guest size and cage occupancy and with a decrease in pressure. The H₂O bending vibrations are relatively independent of guest size, cage occupancy, pressure, temperature, and crystal structure. The guest vibrational modes, especially the bending modes, also have minimal pressure dependence. We have also provided more quantitative evidence that gas hydrate material properties are defined by their hydrogen bond properties, by linking H-bond strength to Young's modulus. The results and ensuing vibrational analysis presented in this paper are a valuable contribution to the ongoing efforts into developing more accurate gas hydrate identification and characterization methods in the laboratory, in industry/nature, and even in outer space.

Effects of the CO₂ Guest Molecule on the sI Clathrate Hydrate Structure

This paper analyzes the structural, energetic and mechanical properties of carbon dioxide hydrate clathrates calculated using finite cluster and periodic ab initio density-functional theory methodologies. Intermolecular interactions are described by the exchange-hole dipole moment method. The stability, gas saturation energetics, guest–host interactions, cage deformations, vibrational frequencies, and equation of state parameters for the low-pressure sI cubic phase of the CO₂@H₂O clathrate hydrate are presented. Our results reveal that: (i) the gas saturation process energetically favors complete filling; (ii) carbon dioxide molecules prefer to occupy the larger of the two cages in the sI structure; (iii) blue shifts occur in both the symmetric and antisymmetric stretching frequencies of CO₂ upon encapsulation; and (iv) free rotation of guest molecules is restricted to a plane parallel to the hexagonal faces of the large cages. In addition, we calculate the librational frequency of the hindered rotation of the guest molecule in the plane perpendicular to the hexagonal faces. Our calculated spectroscopic data can be used as signatures for the detection of clathrate hydrates in planetary environments.

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.

Vlugt TJ. Compressibility, thermal expansion coefficient and heat capacity of CH4 and CO2 hydrate mixtures using molecular dynamics simulations

Understanding the thermal and mechanical properties of CH₄ and CO₂ hydrates is essential for the replacement of CH₄ with CO₂ in natural hydrate deposits as well as for CO₂ sequestration and storage.