Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2003 Dec 9;100(25):14645-50.
doi: 10.1073/pnas.2430913100. Epub 2003 Dec 1.

Thermodynamic stability of hydrogen clathrates

Affiliations

Thermodynamic stability of hydrogen clathrates

Serguei Patchkovskii et al. Proc Natl Acad Sci U S A. .

Abstract

The stability of the recently characterized type II hydrogen clathrate [Mao, W. L., Mao, H.-K., Goncharov, A. F., Struzhkin, V. V., Guo, Q., et al. (2002) Science 297, 2247-2249] with respect to hydrogen occupancy is examined with a statistical mechanical model in conjunction with first-principles quantum chemistry calculations. It is found that the stability of the clathrate is mainly caused by dispersive interactions between H2 molecules and the water forming the cage walls. Theoretical analysis shows that both individual hydrogen molecules and nH2 guest clusters undergo essentially free rotations inside the clathrate cages. Calculations at the experimental conditions--2,000 bar (1 bar = 100 kPa) and 250 K confirm multiple occupancy of the clathrate cages with average occupations of 2.00 and 3.96 H2 molecules per D-5(12) (small) and H-5(12)6(4) (large) cage, respectively. The H2-H2O interactions also are responsible for the experimentally observed softening of the H[bond]H stretching modes. The clathrate is found to be thermodynamically stable at 25 bar and 150 K.

PubMed Disclaimer

Figures

Fig. 1.
Fig. 1.
Structure of the S (D-512) and L (H-51264) cages of the type II ice clathrate. Positions of the hydrogen atoms are omitted for clarity. The coordinate axes correspond to the orientation of the model cages used in the calculations.
Fig. 2.
Fig. 2.
Translational potential energy curve for rigid H2 inside an L (H-51264) cage (a)andS(D-512) cage (b).The H2 molecules are displaced from the center of the cage (R = 0) along the (1, 1, 1) direction (Fig. 1). Calculated average occupancies of the S (H-51264) and L (H-51264) cages at 250 K (c) and 150 K (d).
Fig. 3.
Fig. 3.
Calculated dynamical contributions to enthalpy (Upper) and enthropy (Lower) per mole of clathrated H2 at 150 K. The extreme left point corresponds to H2 in ideal gas. Open and filled symbols denote the results obtained in harmonic and “improved” approximations, respectively. Qualitatively similar dependence (not shown) is found at 250 K.

References

    1. Jeffrey, G. A. (1984) in Inclusion Compounds, eds. Atwood, J. L., Davies, J. E. D. & MacNicol, D. D. (Academic, New York).
    1. Suess, E., Bohrmann, G., Greinert, J. & Lausch, E. (1999) Sci. Am. 281 (5), 52–59.
    1. Sloan, E. D. (2000) Hydrate Engineering Monograph (Soc. Plastic Engineer., Richardson, TX), Vol. 21.
    1. Katz, M. E., Pak, D. K., Dickens, G. R. & Miller, K. G. (1999) Science 286, 1531–1533. - PubMed
    1. Hessello, S. P. P., Grocke, D. R., Jenkyns, K. C., Bjerrum, C. J., Farrimond, P., Bell, H. S. M. & Green, O. R. (2000) Nature 406, 392–395. - PubMed

LinkOut - more resources