Vibrations of Methane in Structure I Clathrate Hydrate—an ab initio Density Functional Molecular Dynamics Study

Formation of a Low-Density Liquid Phase during the Dissociation of Gas Hydrates in Confined Environments

The large amounts of natural gas in a dense solid phase stored in the confined environment of porous materials have become a new, potential method for storing and transporting natural gas. However, there is no experimental evidence to accurately determine the phase state of water during nanoscale gas hydrate dissociation. The results on the dissociation behavior of methane hydrates confined in a nanosilica gel and the contained water phase state during hydrate dissociation at temperatures below the ice point and under atmospheric pressure are presented. Fourier transform infrared spectroscopy (FTIR) and powder X-ray diffraction (PXRD) were used to trace the dissociation of confined methane hydrate synthesized from pore water confined inside the nanosilica gel. The characterization of the confined methane hydrate was also analyzed by PXRD. It was found that the confined methane hydrates dissociated into ultra viscous low-density liquid water (LDL) and methane gas. The results showed that the mechanism of confined methane hydrate dissociation at temperatures below the ice point depended on the phase state of water during hydrate dissociation.

From Infrared Spectra to Macroscopic Mechanical Properties of sH Gas Hydrates through Atomistic Calculations

The vibrational characteristics of gas hydrates are key identifying molecular features of their structure and chemical composition. Density functional theory (DFT)-based IR spectra are one of the efficient tools that can be used to distinguish the vibrational signatures of gas hydrates. In this work, ab initio DFT-based IR technique is applied to analyze the vibrational and mechanical features of structure-H (sH) gas hydrate. IR spectra of different sH hydrates are obtained at 0 K at equilibrium and under applied pressure. Information about the main vibrational modes of sH hydrates and the factors that affect them such as guest type and pressure are revealed. The obtained IR spectra of sH gas hydrates agree with experimental/computational literature values. Hydrogen bond’s vibrational frequencies are used to determine the hydrate’s Young’s modulus which confirms the role of these bonds in defining sH hydrate’s elasticity. Vibrational frequencies depend on pressure and hydrate’s O···O interatomic distance. OH vibrational frequency shifts are related to the OH covalent bond length and present an indication of sH hydrate’s hydrogen bond strength. This work presents a new route to determine mechanical properties for sH hydrate based on IR spectra and contributes to the relatively small database of gas hydrates’ physical and vibrational properties.

Effect of temperature and large guest molecules on the C-H symmetric stretching vibrational frequencies of methane in structure H and I clathrate hydrates

Large molecules such as 2-methylbutane (C₅H₁₂) or 2,2-dimethylbutane (C₆H₁₄) form structure H (sH) hydrates with methane (CH₄) as a help gas. In this study, the Raman spectra of the C–H symmetric stretch region of CH₄ enclathrated within various sH hydrates and structure I CH₄ hydrates were analyzed in the temperature range 137.7–205.4 K. Thermal expansions of these sH hydrate samples were also measured using powder X-ray diffraction. Symmetric stretch vibrational frequencies of CH₄ in host–water cages increased because of varying temperature, and the sizes of the host–water cages also increased; variation of CH₄ in small cages was less than in larger cages. Comparing the variations of the C–H symmetric stretching frequencies of CH₄ in gas hydrates with varying pressure and temperature, we suggest that the observed trend is caused by thermal vibrations of the CH₄ molecule in water cages.

Temperature effect on C–H symmetric stretching frequencies of CH₄ in water cages of sI and sH clathrate hydrates were clarified.

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.

Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases?

The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.

Retracted Article: Effect of temperature and large guest molecules on the C–H symmetric stretching vibrational frequencies of methane in structure H and I clathrate hydrates

Large molecules such as 2-methylbutane (C₅H₁₂) or 2,2-dimethylbutane (C₆H₁₄) form structure H (sH) hydrates with methane (CH₄) as a help gas. In this study, the Raman spectra of the C–H symmetric stretch region of CH₄ enclathrated within various sH hydrates and structure I CH₄ hydrates were analyzed in the temperature range 83–183 K. Thermal expansions of these sH hydrate samples were also measured using powder X-ray diffraction. Symmetric stretch vibrational frequencies of CH₄ in host water cages increased because of varying temperature, and the sizes of the host water cages also increased; variation of CH₄ in small cages was less than in larger cages. Comparing the variations of the C–H symmetric stretching frequencies of CH₄ in gas hydrates with varying pressure and temperature, we suggest that the observed trend is caused by thermal vibrations of the CH₄ molecule in water cages.

Temperature effects on C–H symmetric stretching frequencies of CH₄ in water cages of sI and sH clathrate hydrates were clarified.

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.

DFT calculation of the potential energy landscape topology and Raman spectra of type I CH4and CO2hydrates

Type I methane and carbon dioxide hydrates analyzed using DFT calculations: Raman spectra and potential energy landscape.

C-C stretching Raman spectra and stabilities of hydrocarbon molecules in natural gas hydrates: a quantum chemical study

The presence of specific hydrocarbon gas molecules in various types of water cavities in natural gas hydrates (NGHs) are governed by the relative stabilities of these encapsulated guest molecule-water cavity combinations. Using molecular quantum chemical dispersion-corrected hybrid density functional computations, the interaction (ΔE(host--guest)) and cohesive energies (ΔE(coh)), enthalpies, and Gibbs free energies for the complexes of host water cages and hydrocarbon guest molecules are calculated at the ωB97X-D/6-311++G(2d,2p) level of theory. The zero-point energy effect of ΔE(host-guest) and ΔE(coh) is found to be quite substantial. The energetically optimal host-guest combinations for seven hydrocarbon gas molecules (CH4, C2H6, C3H6, C3H8, C4H8, i-C4H10, and n-C4H10) and various water cavities (D, ID, T, P, H, and I) in NGHs are found to be CH4@D, C2H6@T, C3H6@T, C3H8@T, C4H8@T/P/H, i-C4H10@H, and n-C4H10@H, as the largest cohesive energy magnitudes will be obtained with these host-guest combinations. The stabilities of various water cavities enclosing hydrocarbon molecules are evaluated from the computed cohesive Gibbs free energies: CH4 prefers to be trapped in a ID cage; C2H6 prefer T cages; C3H6 and C3H8 prefer T and H cages; C4H8 and i-C4H10 prefer H cages; and n-C4H10 prefer I cages. The vibrational frequencies and Raman intensities of the C-C stretching vibrational modes for these seven hydrocarbon molecules enclosed in each water cavity are computed. A blue shift results after the guest molecule is trapped from gas phase into various water cages due to the host-guest interactions between the water cage and hydrocarbon molecule. The frequency shifts to the red as the radius of water cages increases. The model calculations support the view that C-C stretching vibrations of hydrocarbon molecules in the water cavities can be used as a tool to identify the types of crystal phases and guest molecules in NGHs.

Calculations of NMR properties for sI and sII clathrate hydrates of methane, ethane and propane

Calculations of NMR parameters (the absolute shielding constants and the spin-spin coupling constants) for 5(12), 5(12)6(2) and 5(12)6(4) cages enclathrating CH4, C2H6 and C3H8 molecules are presented. The DFT/B3LYP/HuzIII-su3 level of theory was employed. The (13)C shielding constants of guest molecules are close to available experimental data. In two cases (the ethane in 5(12) and the propane in 5(12)6(2) cages) the (13)C shielding constants are reported for the first time. Inversion of the methyl/methylene (13)C and (1)H shielding constants order is found for propane in the 5(12)6(2) cage. Topological criteria are used to interpret the changes of values of NMR parameters of water molecules and they establish a connection between single cages and bulk crystal.

A quantum chemistry study of natural gas hydrates

The structure and properties of natural gas hydrates containing hydrocarbons, CO₂, and N₂ molecules were studied by using computational quantum chemistry methods via the density functional theory approach. All host cages involved in I, II, and H types structures where filled with hydrocarbons up to pentanes, CO₂ and N₂ molecules, depending on their size, and the structures of these host-guest systems optimized. Structural properties, vibrational spectra, and density of states were analyzed together with results from atoms-in-a-molecule and natural bond orbitals methods. The inclusion of dispersion terms in the used functional plays a vital role for obtaining reliable information, and thus, B97D functional was shown to be useful for these systems. Results showed remarkable interaction energies, not strongly affected by the type of host cage, with molecules tending to be placed at the center of the cavities when host cages and guest molecules cavities are of similar size, but with molecules approaching hexagonal faces for larger cages. Vibrational properties show remarkable features in certain regions, with shiftings rising from host-guest interactions, and useful patterns in the terahertz region rising from water surface vibrations strongly coupled with guest molecules. Likewise, calculations on crystal systems for the I and H types were carried out using a pseudopotential approach combined with Grimme's method to take account of dispersion.

On the infrared activation of the breathing mode of methane in ice

The IR activation of the ν₁ mode of CH₄ in CH₄ : H₂O ice mixtures is studied using DFT methods for amorphous solids.

Storage capacity and vibration frequencies of guest molecules in CH4 and CO2 hydrates by first-principles calculations

Using first-principle calculations at B97-D/6-311++G(2d,2p) level, we systematically explore the gas capacity of five standard water cavities (5(12), 4(3)5(6)6(3), 5(12)6(2), 5(12)6(4), and 5(12)6(8)) in clathrate hydrate and study the inclusion complexes to infer general trends in vibrational frequencies of guest molecules as a function of cage size and number of guest molecules. In addition, the Raman spectra of hydrates from CO2/CH4 gases are simulated. From our calculations, the maximum cage occupancy of the five considered cages (5(12), 4(3)5(6)6(3), 5(12)6(2), 5(12)6(4), and 5(12)6(8)) is one, one, two, three, and seven for both CH4 and CO2 guest molecules, respectively. Meanwhile, the optimum cage occupancy are one, one, one, two, and four for CO2 molecules and one, one, two, three, and five for CH4 molecules, respectively. Both the C-H stretching frequency of CH4 and the C-O stretching frequency of CO2 gradually decrease as size of the water cages increases. Meanwhile, the C-H stretching frequency gradually increases as the amount of CH4 molecules in the water cavity (e.g., 5(12)6(8)) increases.

Vibrational modes of methane in the structure H clathrate hydrate from ab initio molecular dynamics simulation

Vibrational spectra of guest molecules in clathrate hydrates are frequently measured to determine the characteristic signatures of the molecular environment and dynamical properties of guest-host interactions. Here, we present results of our study on the vibrational frequencies of methane molecules in structure H clathrate hydrates, namely, in the 5(12) and 4(3)5(6)6(3) cages, as the frequencies of stretching vibrational modes in these environments are still unclear. The vibrational spectra of methane molecules in structure H clathrate hydrate were obtained from ab initio molecular dynamics simulation and computed from Fourier transform of autocorrelation functions for each distinct vibrational mode. The calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules were found to be lower in the 4(3)5(6)6(3) cages than in the 5(12) cages (3.8 cm(-1) for symmetric stretching and 6.0 cm(-1) for asymmetric stretching). The C-H bond length and average distance between methane molecules and host-water molecules in 4(3)5(6)6(3) cages were slightly longer than those in the 5(12) cages.

Molecular vibrations of methane molecules in the structure I clathrate hydrate from ab initio molecular dynamics simulation

Vibrational frequencies of guest molecules in clathrate hydrates reflect the molecular environment and dynamical behavior of molecules. A detailed understanding of the mechanism for the vibrational frequency changes of the guest molecules in the clathrate hydrate cages is still incomplete. In this study, molecular vibrations of methane molecules in a structure I clathrate hydrate are calculated from ab initio molecular dynamics simulation. The vibrational spectra of methane are computed by Fourier transform of autocorrelation functions, which reveal distinct separation of each vibrational mode. Calculated symmetric and asymmetric stretching vibrational frequencies of methane molecules are lower in the large cages than in the small cages (8 and 16 cm(-1) for symmetric and asymmetric stretching, respectively). These changes are closely linked with the C-H bond length. The vibrational frequencies for the bending and rocking vibrational modes nearly overlap in each of the cages.

Stability, adsorption, and diffusion of CH₄, CO₂, and H₂ in clathrate hydrates

We present a study of the adsorption and diffusion of CH₄, CO₂, and H₂ molecules in clathrate hydrates using ab initio van der Waals density functional formalism [M. Dion, Phys. Rev. Lett. 92, 246401 (2004)10.1103/PhysRevLett.92.246401]. We find that the adsorption energy is dominated by van der Waals interactions and that, without them, gas hydrates would not be stable. We calculate the maximum adsorption capacity as well as the maximum hydrocarbon size that can be adsorbed. The relaxation of the host lattice is essential for a good description of the diffusion activation energies, which are estimated to be of the order of 0.2, 0.4, and 1.0 eV for H₂, CO₂, and CH₄, respectively.

Application of the independent molecule model to elucidate the dynamics of structure I methane hydrate

Raman spectroscopy has exhibited the C-H stretch (A1 mode) frequency nu1 of hydrated methanes at 2915 cm(-1) for the 5(12) cage and 2905 cm(-1) for the 5(12)6(2) cage. These values are lower than the frequency of 2916.5 cm(-1) in gaseous methane. In this paper, we theoretically examine the Raman spectra observed in methane hydrate by normal mode analysis using the independent molecule model. By a breakdown of the symmetry, the four frequencies in modes A1, E, T2 and T2 observed in gaseous methane are separated into nine frequencies in the hydrate. It is necessary to consider the anharmonic potential energy within methane and hydrogen bonding between methane hydrogen and water oxygen in order to get a result in qualitative agreement with experiment. The frequency in the 5(12)6(2) cage is shifted downward in comparison to the one in 5(12), and the frequencies in the both cages are also shifted downward compared with the frequencies in gas. Calculations are also reported for the isotopic methane (CD4, 13CH4) hydrates.