Elsevier

Chemical Physics Letters

Volume 648, 16 March 2016, Pages 75-80
Chemical Physics Letters

Understanding effect of structure and stability on transformation of CH4 hydrate to CO2 hydrate

https://doi.org/10.1016/j.cplett.2016.02.004Get rights and content

Highlights

  • CH4 molecules are likely to reside in 512 cages rather than 51262 cages.

  • CO2 molecules are likely to reside in 51262 cages rather than 512 cages.

  • CO2 has smaller interaction energy than CH4 in the same cage.

  • The entropy dominates the transformation of CH4 hydrate into CO2 hydrate.

  • In the same cage CO2 molecules move much slowly than CH4 molecules.

Abstract

Understanding the transformation process of CH4 hydrate to CO2 hydrate is crucial to develop the CH4single bondCO2 replacement technique for CH4 production and CO2 sequestration. Ab initio calculations show that the transformation will slightly distort the host lattice and decrease the binding strength of guest molecules, but it is a thermodynamically spontaneous process dominated by the entropic contribution. Moreover, ab initio molecular dynamics simulations suggest that the dynamics of the host lattice is independent on the guest molecules, while CO2 in hydrate exhibits slower translational and rotational motion than CH4 in hydrate.

Introduction

Clathrate hydrates are nonstoichiometric crystalline compounds in which small guest molecules are encaged by hydrogen-bonded water cages [1], [2], [3]. Depending on the size and shape of guest molecules, three common crystalline structures of clathrate hydrates have been uncovered and addressed as structure-I (sI, cubic, space group Pm3¯n), structure-II (sII, cubic, Fd3¯m) and structure-H (sH, hexagonal, P6/mmm) [4]. Clathrate hydrates typically form at the ambient temperature (less than 300 K) and moderate pressures (more than 0.6 MPa). They are abundant in Earth's permafrost and subsea sediments [5], [6], and also have been found in other planetary bodies [7], [8]. The most conservative estimates suggest that the amount of carbon in hydrates is equivalent to twice that of other fossil fuels combined [9], [10], and the energy density in hydrates is approximately as the same as that of a compressed gas, that is, one volume of hydrate will dissociate to 180 volumes of gas. Therefore, clathrate hydrates have been evaluated both technically and economically as a potential energy resource and attract increasing worldwide interest [11], [12], [13].

The replacement of CH4 with CO2 in clathrate hydrates has been regarded as a promising method for natural gas production [14], [15], [16], [17], [18], [19], which also provides a way to storage CO2 for a long term [20], [21]. By injecting CO2 to CH4 hydrate deposit, CH4 release from hydrates can be facilitated while the formation of CO2 hydrates can be largely increased. To analyze this complex replacement phenomenon, it is necessary to understand the stability and physical-chemical properties of CH4 and CO2 hydrates. Actually, CH4 and CO2 both are likely to form sI structure [4], [22], which consists of 46 water molecules in one unit cell, forming two dodecahedron (12 pentagonal faces; 512) and six tetradecahedron (12 pentagonal and 2 hexagonal faces; 51262) cages. In addition, researchers have verified that CH4single bondCO2 mixed hydrates can be formed during the replacement process [23], [24], but little attention is paid to the structure and stability of CH4single bondCO2 mixed hydrates. Hence, investigating the properties of CH4single bondCO2 mixed hydrates could benefit CH4 production as well as CO2 sequestration.

Computational simulation technique has been proved to be one of the most effective methods to explore the microscopic information of the hydrates [25], [26], [27]. By considering the electronic interactions through van der Waals density functional theory, Román-Pérez et al. [28] have suggested that both CO2 and CH4 have larger interaction energies in 51262 cage than that in 512 cage, and CH4 always has a larger value than CO2 in the same cage. Whereas, Srivastava and Sastry [29] tested the impact of different functionals and basis sets on interaction energy and demonstrated that CO2 is energetically favorable to reside in the 512 cages. Yi and coworkers [21] have recently investigated the growth mechanism of CH4single bondCO2 mixed hydrate using molecular dynamics simulations and found that CO2 molecules are preferably encaged into 51262 cages. Moreover, through classical molecular dynamics simulations, Geng et al. [23] have pointed out that CH4single bondCO2 mixed hydrate is more stable than pure CH4 hydrate or CO2 hydrate. Hence, in this work, we firstly performed ab initio calculations to optimize CH4 hydrate, CH4single bondCO2 mixed hydrate and CO2 hydrate. Structural and energetic characteristics of hydrates were compared to explore the possibility of transformation of CH4 hydrate to CH4single bondCO2 mixed hydrate and CO2 hydrate. In addition, we carried out ab initio molecular dynamics simulations to observe the structural and thermal fluctuations of hydrates in a more realistic environment. The translational and rotational motions of molecules in hydrates were analyzed to evaluate the kinetic stability of hydrates. Through this work, we could gain a better insight of the transformation of CH4 hydrate to CO2 hydrate, which would helpful to understand the CH4single bondCO2 replacement mechanism.

Section snippets

Computational details

The crystal structure of sI clathrate hydrate was taken from Lenz's data [30]. In pure CH4 hydrate and pure CO2 hydrate, the 512 and 51262 cages were fully filled with CH4 and CO2 molecules, respectively. For CH4single bondCO2 mixed hydrate, two structures were considered: 512 cages were fully filled with CO2 molecules and 51262 cages were fully filled with CH4 molecules, denoting as Mixed-I; 512 cages were fully filled with CH4 molecules and 51262 cages were fully filled with CO2 molecules, denoting as

Optimized structure and thermodynamic stability

Figure 1 shows the optimized structures of CH4 hydrate, Mixed-I hydrate, Mixed-II hydrate and CO2 hydrate. The host lattice of four studied hydrates is nearly same, and CH4 and CO2 molecules are basically resided in the center of the clathrate cages. Table 1 shows that the averaged hydrogen bond length of CH4 hydrate is ∼2.699 Å, which matches Martos-Villa's results [36]. When CO2 molecules are encaged by 512 or/and 51262 cages, the corresponding hydrates have a slightly increased hydrogen bond

Conclusions

Ab initio calculations and simulations were performed to investigate structure and stability of CH4 hydrate, CH4single bondCO2 mixed hydrate and CO2 hydrate, and their relationship to the transformation of CH4 hydrate into CO2 hydrate. The results show that CO2 molecules encaged by 512 cages can induce the distortion of host lattice, while 51262 cages are nicely suitable to accommodate CO2 molecules. In both 512 and 51262 cages, CO2 molecules always have smaller interaction energies than CH4 molecules

Acknowledgements

This work was supported by grants from Natural Science Foundation of Shandong Province (BS2015NJ007), National Natural Science Foundation of China (11504133 and 1137412), and National Basic Research Program of China (2015CB251200).

References (39)

  • B.K. Chastain et al.

    Planet. Space Sci.

    (2007)
  • S.-Y. Lee et al.

    Fuel Process. Technol.

    (2001)
  • Y.F. Makogon et al.

    J. Nat. Gas Sci. Eng.

    (2013)
  • H. Dashti et al.

    J. Nat. Gas Sci. Eng.

    (2015)
  • T. Uchida et al.

    Energy Convers. Manage.

    (1997)
  • M. Ota et al.

    Energy Convers. Manage.

    (2005)
  • M. Ota et al.

    Fluid Phase Equilib.

    (2005)
  • L. Yi et al.

    J. Energy Chem.

    (2014)
  • L. Xu et al.

    Comput. Theor. Chem.

    (2011)
  • Y. Qi et al.

    Energy Convers. Manage.

    (2011)
  • R. Martos-Villa et al.

    J. Mol. Graph. Model.

    (2013)
  • J. Liu et al.

    Chem. Phys. Lett.

    (2015)
  • T.C.W. Mak et al.

    J. Chem. Phys.

    (1965)
  • R.K. McMullan et al.

    J. Chem. Phys.

    (1965)
  • W.L. Vos et al.

    Phys. Rev. Lett.

    (1993)
  • E.D. Sloan

    Nature

    (2003)
  • J. Marcelle-De Silva et al.

    Energies

    (2011)
  • K.A. Kvenvolden

    Rev. Geophys.

    (1993)
  • J.S. Loveday et al.

    Nature

    (2001)
  • Cited by (28)

    • Molecular guest exchange and subsequent structural transformation in CH<inf>4</inf> – CO<inf>2</inf> replacement occurring in sH hydrates as revealed by <sup>13</sup>C NMR spectroscopy and molecular dynamic simulations

      2023, Chemical Engineering Journal
      Citation Excerpt :

      Recent studies on molecular dynamics (MD) simulations examined guest–guest, guest–host, and host–host interactions in CH4 + CO2 hydrates at a molecular scale, suggesting that host–guest interactions and guest orientation affect the structure of the host frameworks [29,30]. MD simulations also suggested that the replacement process distorted the host lattices in CH4 – CO2 replacement in sI hydrates and that CO2 molecules exhibited a slower rotational and translational motion than did CH4 molecules in hydrate lattices [31]. Moreover, a series of microsecond MD simulations showed that the chemical potentials of guest molecules and the barrier of mass transfer could play important roles in the replacement process [32].

    • CO<inf>2</inf> enhanced gas recovery and sequestration in depleted gas reservoirs: A review

      2021, Journal of Petroleum Science and Engineering
      Citation Excerpt :

      Thus, the guest distribution and the exchange behavior depend on the cage. The transformation from CH4 hydrate into CO2 hydrate has exothermic nature with a slight distortion of the host structure (J. Liu et al., 2016). The CH4/CO2 substitution takes place when the partial pressure of CH4 is less than the equilibrium pressure required to constitute pure CH4 hydrate.

    • Research progress and challenges in hydrate-based carbon dioxide capture applications

      2020, Applied Energy
      Citation Excerpt :

      SII CH4 + CO2 binary hydrates and CH4 + CO2 + THF ternary hydrates can form in THF solution [114]. SH CH4–CO2 binary hydrates are more stable than pure CH4 or CO2 hydrates [172]. The stability of hydrate is gradually enhanced as the number of guest molecules (CH4 and CO2 simultaneously) in 51268 cages increases up to five [163].

    View all citing articles on Scopus
    View full text