Cage occupancies, lattice constants, and guest chemical potentials for structure II hydrogen clathrate hydrate from Gibbs ensemble Monte Carlo simulations

Paul E. Brumby, Daisuke Yuhara, Tomohiro Hasegawa, David T. Wu, Amadeu K. Sum, Kenji Yasuoka

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

In this paper, equilibrium properties of structure II hydrates of hydrogen were determined from Monte Carlo simulations in the isothermal-isobaric Gibbs ensemble. Water and hydrogen molecules are described by the TIP4P/Ice and Silvera-Goldman models, respectively. The use of the Gibbs ensemble has many key advantages for the simulation of hydrates. By the separation of hydrogen vapor and hydrate phases into their own domains, coupled with transfer moves of hydrogen molecules between domains, cage occupancies were determined. Furthermore, the choice of this ensemble also allows equilibrium lattice constants and guest molecule chemical potentials to be straightforwardly estimated. Results for hydrogen mass fractions indicate reasonable agreement with prior simulation data and theoretical models, while detailed analysis of cage occupancy distributions and neighboring cage pair occupancy combinations gives valuable insight into the behavior of this hydrate at the inter-cage scale. These results will aid in the construction of theoretical models, for which knowledge of the occupancy of neighboring cages is of great importance. In support of previous experimental and theoretical works, we also find evidence of double occupancy of a few small cages inside of the hydrate stability zone, albeit at very high pressures; approximately 0.1% of small cages are doubly occupied at 300 MPa, for temperatures of 225 K and 250 K.

Original languageEnglish
Article number134503
JournalJournal of Chemical Physics
Volume150
Issue number13
DOIs
Publication statusPublished - 2019 Apr 7

Fingerprint

Chemical potential
clathrates
Hydrates
hydrates
Lattice constants
Hydrogen
hydrogen
simulation
Molecules
molecules
Ice
data simulation
ice
Vapors
Monte Carlo simulation
vapor phases
Water
water
Temperature
temperature

ASJC Scopus subject areas

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

Cite this

Cage occupancies, lattice constants, and guest chemical potentials for structure II hydrogen clathrate hydrate from Gibbs ensemble Monte Carlo simulations. / Brumby, Paul E.; Yuhara, Daisuke; Hasegawa, Tomohiro; Wu, David T.; Sum, Amadeu K.; Yasuoka, Kenji.

In: Journal of Chemical Physics, Vol. 150, No. 13, 134503, 07.04.2019.

Research output: Contribution to journalArticle

@article{7f79c545f8294b6a9f91c4c57d10914f,
title = "Cage occupancies, lattice constants, and guest chemical potentials for structure II hydrogen clathrate hydrate from Gibbs ensemble Monte Carlo simulations",
abstract = "In this paper, equilibrium properties of structure II hydrates of hydrogen were determined from Monte Carlo simulations in the isothermal-isobaric Gibbs ensemble. Water and hydrogen molecules are described by the TIP4P/Ice and Silvera-Goldman models, respectively. The use of the Gibbs ensemble has many key advantages for the simulation of hydrates. By the separation of hydrogen vapor and hydrate phases into their own domains, coupled with transfer moves of hydrogen molecules between domains, cage occupancies were determined. Furthermore, the choice of this ensemble also allows equilibrium lattice constants and guest molecule chemical potentials to be straightforwardly estimated. Results for hydrogen mass fractions indicate reasonable agreement with prior simulation data and theoretical models, while detailed analysis of cage occupancy distributions and neighboring cage pair occupancy combinations gives valuable insight into the behavior of this hydrate at the inter-cage scale. These results will aid in the construction of theoretical models, for which knowledge of the occupancy of neighboring cages is of great importance. In support of previous experimental and theoretical works, we also find evidence of double occupancy of a few small cages inside of the hydrate stability zone, albeit at very high pressures; approximately 0.1{\%} of small cages are doubly occupied at 300 MPa, for temperatures of 225 K and 250 K.",
author = "Brumby, {Paul E.} and Daisuke Yuhara and Tomohiro Hasegawa and Wu, {David T.} and Sum, {Amadeu K.} and Kenji Yasuoka",
year = "2019",
month = "4",
day = "7",
doi = "10.1063/1.5084785",
language = "English",
volume = "150",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "13",

}

TY - JOUR

T1 - Cage occupancies, lattice constants, and guest chemical potentials for structure II hydrogen clathrate hydrate from Gibbs ensemble Monte Carlo simulations

AU - Brumby, Paul E.

AU - Yuhara, Daisuke

AU - Hasegawa, Tomohiro

AU - Wu, David T.

AU - Sum, Amadeu K.

AU - Yasuoka, Kenji

PY - 2019/4/7

Y1 - 2019/4/7

N2 - In this paper, equilibrium properties of structure II hydrates of hydrogen were determined from Monte Carlo simulations in the isothermal-isobaric Gibbs ensemble. Water and hydrogen molecules are described by the TIP4P/Ice and Silvera-Goldman models, respectively. The use of the Gibbs ensemble has many key advantages for the simulation of hydrates. By the separation of hydrogen vapor and hydrate phases into their own domains, coupled with transfer moves of hydrogen molecules between domains, cage occupancies were determined. Furthermore, the choice of this ensemble also allows equilibrium lattice constants and guest molecule chemical potentials to be straightforwardly estimated. Results for hydrogen mass fractions indicate reasonable agreement with prior simulation data and theoretical models, while detailed analysis of cage occupancy distributions and neighboring cage pair occupancy combinations gives valuable insight into the behavior of this hydrate at the inter-cage scale. These results will aid in the construction of theoretical models, for which knowledge of the occupancy of neighboring cages is of great importance. In support of previous experimental and theoretical works, we also find evidence of double occupancy of a few small cages inside of the hydrate stability zone, albeit at very high pressures; approximately 0.1% of small cages are doubly occupied at 300 MPa, for temperatures of 225 K and 250 K.

AB - In this paper, equilibrium properties of structure II hydrates of hydrogen were determined from Monte Carlo simulations in the isothermal-isobaric Gibbs ensemble. Water and hydrogen molecules are described by the TIP4P/Ice and Silvera-Goldman models, respectively. The use of the Gibbs ensemble has many key advantages for the simulation of hydrates. By the separation of hydrogen vapor and hydrate phases into their own domains, coupled with transfer moves of hydrogen molecules between domains, cage occupancies were determined. Furthermore, the choice of this ensemble also allows equilibrium lattice constants and guest molecule chemical potentials to be straightforwardly estimated. Results for hydrogen mass fractions indicate reasonable agreement with prior simulation data and theoretical models, while detailed analysis of cage occupancy distributions and neighboring cage pair occupancy combinations gives valuable insight into the behavior of this hydrate at the inter-cage scale. These results will aid in the construction of theoretical models, for which knowledge of the occupancy of neighboring cages is of great importance. In support of previous experimental and theoretical works, we also find evidence of double occupancy of a few small cages inside of the hydrate stability zone, albeit at very high pressures; approximately 0.1% of small cages are doubly occupied at 300 MPa, for temperatures of 225 K and 250 K.

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

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

U2 - 10.1063/1.5084785

DO - 10.1063/1.5084785

M3 - Article

C2 - 30954046

AN - SCOPUS:85063899358

VL - 150

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 13

M1 - 134503

ER -