H 2 , HD, and D 2 in the small cage of structure II clathrate hydrate: Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates

David Lauvergnat, Peter Felker, Yohann Scribano, David M. Benoit, Zlatko Bacic

Research output: Contribution to journalArticle

Abstract

We report the first fully coupled quantum six-dimensional (6D) bound-state calculations of the vibration-translation-rotation eigenstates of a flexible H 2 , HD, and D 2 molecule confined inside the small cage of the structure II clathrate hydrate embedded in larger hydrate domains with up to 76 H 2 O molecules, treated as rigid. Our calculations use a pairwise-additive 6D intermolecular potential energy surface for H 2 in the hydrate domain, based on an ab initio 6D H 2 -H 2 O pair potential for flexible H 2 and rigid H 2 O. They extend to the first excited (v = 1) vibrational state of H 2 , along with two isotopologues, providing a direct computation of vibrational frequency shifts. We show that obtaining a converged v = 1 vibrational state of the caged molecule does not require converging the very large number of intermolecular translation-rotation states belonging to the v = 0 manifold up to the energy of the intramolecular stretch fundamental (≈4100 cm −1 for H 2 ). Only a relatively modest-size basis for the intermolecular degrees of freedom is needed to accurately describe the vibrational averaging over the delocalized wave function of the quantum ground state of the system. For the caged H 2 , our computed fundamental translational excitations, rotational j = 0 → 1 transitions, and frequency shifts of the stretch fundamental are in excellent agreement with recent quantum 5D (rigid H 2 ) results [A. Powers et al., J. Chem. Phys. 148, 144304 (2018)]. Our computed frequency shift of −43 cm −1 for H 2 is only 14% away from the experimental value at 20 K.

Original languageEnglish (US)
Article number154303
JournalJournal of Chemical Physics
Volume150
Issue number15
DOIs
StatePublished - Apr 21 2019

Fingerprint

clathrates
Vibrational spectra
Hydrates
hydrates
frequency shift
eigenvectors
vibrational states
vibration
Molecules
molecules
Potential energy surfaces
Wave functions
Ground state
degrees of freedom
potential energy
wave functions
ground state
excitation
energy

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

H 2 , HD, and D 2 in the small cage of structure II clathrate hydrate : Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates. / Lauvergnat, David; Felker, Peter; Scribano, Yohann; Benoit, David M.; Bacic, Zlatko.

In: Journal of Chemical Physics, Vol. 150, No. 15, 154303, 21.04.2019.

Research output: Contribution to journalArticle

@article{98c9f4d7fc3c4751b1c8bc90c8cc5a8c,
title = "H 2 , HD, and D 2 in the small cage of structure II clathrate hydrate: Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates",
abstract = "We report the first fully coupled quantum six-dimensional (6D) bound-state calculations of the vibration-translation-rotation eigenstates of a flexible H 2 , HD, and D 2 molecule confined inside the small cage of the structure II clathrate hydrate embedded in larger hydrate domains with up to 76 H 2 O molecules, treated as rigid. Our calculations use a pairwise-additive 6D intermolecular potential energy surface for H 2 in the hydrate domain, based on an ab initio 6D H 2 -H 2 O pair potential for flexible H 2 and rigid H 2 O. They extend to the first excited (v = 1) vibrational state of H 2 , along with two isotopologues, providing a direct computation of vibrational frequency shifts. We show that obtaining a converged v = 1 vibrational state of the caged molecule does not require converging the very large number of intermolecular translation-rotation states belonging to the v = 0 manifold up to the energy of the intramolecular stretch fundamental (≈4100 cm −1 for H 2 ). Only a relatively modest-size basis for the intermolecular degrees of freedom is needed to accurately describe the vibrational averaging over the delocalized wave function of the quantum ground state of the system. For the caged H 2 , our computed fundamental translational excitations, rotational j = 0 → 1 transitions, and frequency shifts of the stretch fundamental are in excellent agreement with recent quantum 5D (rigid H 2 ) results [A. Powers et al., J. Chem. Phys. 148, 144304 (2018)]. Our computed frequency shift of −43 cm −1 for H 2 is only 14{\%} away from the experimental value at 20 K.",
author = "David Lauvergnat and Peter Felker and Yohann Scribano and Benoit, {David M.} and Zlatko Bacic",
year = "2019",
month = "4",
day = "21",
doi = "10.1063/1.5090573",
language = "English (US)",
volume = "150",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "15",

}

TY - JOUR

T1 - H 2 , HD, and D 2 in the small cage of structure II clathrate hydrate

T2 - Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates

AU - Lauvergnat, David

AU - Felker, Peter

AU - Scribano, Yohann

AU - Benoit, David M.

AU - Bacic, Zlatko

PY - 2019/4/21

Y1 - 2019/4/21

N2 - We report the first fully coupled quantum six-dimensional (6D) bound-state calculations of the vibration-translation-rotation eigenstates of a flexible H 2 , HD, and D 2 molecule confined inside the small cage of the structure II clathrate hydrate embedded in larger hydrate domains with up to 76 H 2 O molecules, treated as rigid. Our calculations use a pairwise-additive 6D intermolecular potential energy surface for H 2 in the hydrate domain, based on an ab initio 6D H 2 -H 2 O pair potential for flexible H 2 and rigid H 2 O. They extend to the first excited (v = 1) vibrational state of H 2 , along with two isotopologues, providing a direct computation of vibrational frequency shifts. We show that obtaining a converged v = 1 vibrational state of the caged molecule does not require converging the very large number of intermolecular translation-rotation states belonging to the v = 0 manifold up to the energy of the intramolecular stretch fundamental (≈4100 cm −1 for H 2 ). Only a relatively modest-size basis for the intermolecular degrees of freedom is needed to accurately describe the vibrational averaging over the delocalized wave function of the quantum ground state of the system. For the caged H 2 , our computed fundamental translational excitations, rotational j = 0 → 1 transitions, and frequency shifts of the stretch fundamental are in excellent agreement with recent quantum 5D (rigid H 2 ) results [A. Powers et al., J. Chem. Phys. 148, 144304 (2018)]. Our computed frequency shift of −43 cm −1 for H 2 is only 14% away from the experimental value at 20 K.

AB - We report the first fully coupled quantum six-dimensional (6D) bound-state calculations of the vibration-translation-rotation eigenstates of a flexible H 2 , HD, and D 2 molecule confined inside the small cage of the structure II clathrate hydrate embedded in larger hydrate domains with up to 76 H 2 O molecules, treated as rigid. Our calculations use a pairwise-additive 6D intermolecular potential energy surface for H 2 in the hydrate domain, based on an ab initio 6D H 2 -H 2 O pair potential for flexible H 2 and rigid H 2 O. They extend to the first excited (v = 1) vibrational state of H 2 , along with two isotopologues, providing a direct computation of vibrational frequency shifts. We show that obtaining a converged v = 1 vibrational state of the caged molecule does not require converging the very large number of intermolecular translation-rotation states belonging to the v = 0 manifold up to the energy of the intramolecular stretch fundamental (≈4100 cm −1 for H 2 ). Only a relatively modest-size basis for the intermolecular degrees of freedom is needed to accurately describe the vibrational averaging over the delocalized wave function of the quantum ground state of the system. For the caged H 2 , our computed fundamental translational excitations, rotational j = 0 → 1 transitions, and frequency shifts of the stretch fundamental are in excellent agreement with recent quantum 5D (rigid H 2 ) results [A. Powers et al., J. Chem. Phys. 148, 144304 (2018)]. Our computed frequency shift of −43 cm −1 for H 2 is only 14% away from the experimental value at 20 K.

UR - http://www.scopus.com/inward/record.url?scp=85064476607&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85064476607&partnerID=8YFLogxK

U2 - 10.1063/1.5090573

DO - 10.1063/1.5090573

M3 - Article

C2 - 31005099

AN - SCOPUS:85064476607

VL - 150

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 15

M1 - 154303

ER -