Influence of combining rules on the cavity occupancy of clathrate hydrates by Monte Carlo simulations

Nikolaos I. Papadimitriou, Ioannis N. Tsimpanogiannis, Ioannis Economou, Athanassios K. Stubos

Research output: Contribution to journalArticle

26 Citations (Scopus)

Abstract

Assessing the exact amount of gas stored in clathrate-hydrate structures can be addressed by either molecular-level simulations (e.g. Monte Carlo) or continuum-level modelling (e.g. van der Waals-Platteeuw-theory-based models). In either case, the Lorentz-Berthelot (LB) combining rules are by far the most common approach for the evaluation of the parameters between the different types of atoms that form the hydrate structure. The effect of combining rules on the calculations has not been addressed adequately in the hydrate-related literature. Only recently the use of the LB combining rules in hydrate studies has been questioned.In the current study, we report an extensive series of Grand Canonical Monte Carlo simulations along the three-phase (H-Lw-V) equilibrium curve. The exact geometry of hydrate crystals is known from diffraction experiments and, therefore, the formation of hydrates can be simulated as a process of gas adsorption in a solid porous material. We examine the effect of deviations from the LB combining rules on the cavity occupancy of argon hydrates and work towards quantifying it. The specific system is selected as a result of the characteristic behaviour of argon to form hydrates of different structures depending on the prevailing pressure. In particular, an sII hydrate is formed at lower pressures, while an sI hydrate is formed at intermediate pressures, and finally an sH hydrate is formed at higher pressures.

Original languageEnglish
Pages (from-to)2258-2274
Number of pages17
JournalMolecular Physics
Volume112
Issue number17
DOIs
Publication statusPublished - 2 Sep 2014

Fingerprint

clathrates
Hydrates
hydrates
Pressure
cavities
Argon
simulation
Gases
Adsorption
Monte Carlo simulation
argon
Gas adsorption
porous materials
gases
Porous materials
low pressure
Diffraction
continuums
deviation
Atoms

Keywords

  • argon hydrate
  • cavity occupancy
  • clathrate hydrate
  • combining rules
  • Grand Canonical Monte Carlo
  • parametric study

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Condensed Matter Physics
  • Biophysics
  • Molecular Biology

Cite this

Influence of combining rules on the cavity occupancy of clathrate hydrates by Monte Carlo simulations. / Papadimitriou, Nikolaos I.; Tsimpanogiannis, Ioannis N.; Economou, Ioannis; Stubos, Athanassios K.

In: Molecular Physics, Vol. 112, No. 17, 02.09.2014, p. 2258-2274.

Research output: Contribution to journalArticle

Papadimitriou, Nikolaos I. ; Tsimpanogiannis, Ioannis N. ; Economou, Ioannis ; Stubos, Athanassios K. / Influence of combining rules on the cavity occupancy of clathrate hydrates by Monte Carlo simulations. In: Molecular Physics. 2014 ; Vol. 112, No. 17. pp. 2258-2274.
@article{da1ca041665a4864b870a903f46951fb,
title = "Influence of combining rules on the cavity occupancy of clathrate hydrates by Monte Carlo simulations",
abstract = "Assessing the exact amount of gas stored in clathrate-hydrate structures can be addressed by either molecular-level simulations (e.g. Monte Carlo) or continuum-level modelling (e.g. van der Waals-Platteeuw-theory-based models). In either case, the Lorentz-Berthelot (LB) combining rules are by far the most common approach for the evaluation of the parameters between the different types of atoms that form the hydrate structure. The effect of combining rules on the calculations has not been addressed adequately in the hydrate-related literature. Only recently the use of the LB combining rules in hydrate studies has been questioned.In the current study, we report an extensive series of Grand Canonical Monte Carlo simulations along the three-phase (H-Lw-V) equilibrium curve. The exact geometry of hydrate crystals is known from diffraction experiments and, therefore, the formation of hydrates can be simulated as a process of gas adsorption in a solid porous material. We examine the effect of deviations from the LB combining rules on the cavity occupancy of argon hydrates and work towards quantifying it. The specific system is selected as a result of the characteristic behaviour of argon to form hydrates of different structures depending on the prevailing pressure. In particular, an sII hydrate is formed at lower pressures, while an sI hydrate is formed at intermediate pressures, and finally an sH hydrate is formed at higher pressures.",
keywords = "argon hydrate, cavity occupancy, clathrate hydrate, combining rules, Grand Canonical Monte Carlo, parametric study",
author = "Papadimitriou, {Nikolaos I.} and Tsimpanogiannis, {Ioannis N.} and Ioannis Economou and Stubos, {Athanassios K.}",
year = "2014",
month = "9",
day = "2",
doi = "10.1080/00268976.2014.902136",
language = "English",
volume = "112",
pages = "2258--2274",
journal = "Molecular Physics",
issn = "0026-8976",
publisher = "Taylor and Francis Ltd.",
number = "17",

}

TY - JOUR

T1 - Influence of combining rules on the cavity occupancy of clathrate hydrates by Monte Carlo simulations

AU - Papadimitriou, Nikolaos I.

AU - Tsimpanogiannis, Ioannis N.

AU - Economou, Ioannis

AU - Stubos, Athanassios K.

PY - 2014/9/2

Y1 - 2014/9/2

N2 - Assessing the exact amount of gas stored in clathrate-hydrate structures can be addressed by either molecular-level simulations (e.g. Monte Carlo) or continuum-level modelling (e.g. van der Waals-Platteeuw-theory-based models). In either case, the Lorentz-Berthelot (LB) combining rules are by far the most common approach for the evaluation of the parameters between the different types of atoms that form the hydrate structure. The effect of combining rules on the calculations has not been addressed adequately in the hydrate-related literature. Only recently the use of the LB combining rules in hydrate studies has been questioned.In the current study, we report an extensive series of Grand Canonical Monte Carlo simulations along the three-phase (H-Lw-V) equilibrium curve. The exact geometry of hydrate crystals is known from diffraction experiments and, therefore, the formation of hydrates can be simulated as a process of gas adsorption in a solid porous material. We examine the effect of deviations from the LB combining rules on the cavity occupancy of argon hydrates and work towards quantifying it. The specific system is selected as a result of the characteristic behaviour of argon to form hydrates of different structures depending on the prevailing pressure. In particular, an sII hydrate is formed at lower pressures, while an sI hydrate is formed at intermediate pressures, and finally an sH hydrate is formed at higher pressures.

AB - Assessing the exact amount of gas stored in clathrate-hydrate structures can be addressed by either molecular-level simulations (e.g. Monte Carlo) or continuum-level modelling (e.g. van der Waals-Platteeuw-theory-based models). In either case, the Lorentz-Berthelot (LB) combining rules are by far the most common approach for the evaluation of the parameters between the different types of atoms that form the hydrate structure. The effect of combining rules on the calculations has not been addressed adequately in the hydrate-related literature. Only recently the use of the LB combining rules in hydrate studies has been questioned.In the current study, we report an extensive series of Grand Canonical Monte Carlo simulations along the three-phase (H-Lw-V) equilibrium curve. The exact geometry of hydrate crystals is known from diffraction experiments and, therefore, the formation of hydrates can be simulated as a process of gas adsorption in a solid porous material. We examine the effect of deviations from the LB combining rules on the cavity occupancy of argon hydrates and work towards quantifying it. The specific system is selected as a result of the characteristic behaviour of argon to form hydrates of different structures depending on the prevailing pressure. In particular, an sII hydrate is formed at lower pressures, while an sI hydrate is formed at intermediate pressures, and finally an sH hydrate is formed at higher pressures.

KW - argon hydrate

KW - cavity occupancy

KW - clathrate hydrate

KW - combining rules

KW - Grand Canonical Monte Carlo

KW - parametric study

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

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

U2 - 10.1080/00268976.2014.902136

DO - 10.1080/00268976.2014.902136

M3 - Article

VL - 112

SP - 2258

EP - 2274

JO - Molecular Physics

JF - Molecular Physics

SN - 0026-8976

IS - 17

ER -