Related Experiment Video
Updated: May 4, 2026

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry
Published on: April 8, 2020
Carbon dioxide hydrate phase equilibrium and cage occupancy calculations using ab initio intermolecular potentials.
Srinath C Velaga1, Brian J Anderson
1National Energy Technological Laboratory , Morgantown, West Virginia 26506, United States.
Accurate carbon dioxide (CO2) cage occupancy in gas hydrates is crucial for CO2 sequestration. This study developed a new model, the VAS model, using ab initio calculations to accurately predict CO2 hydrate properties and cage occupancies.
Area of Science:
- Geochemistry
- Materials Science
- Computational Chemistry
Background:
- Gas hydrate deposits are promising for CO2 sequestration, involving methane replacement.
- Accurate characterization of CO2 cage occupancy is vital for assessing sequestration potential and methane recovery.
- Existing methods for determining guest-host interaction potentials have limitations and inconsistencies.
Purpose of the Study:
- To develop and validate a reliable method for predicting accurate cage occupancies of CO2 in gas hydrates.
- To establish accurate intermolecular potentials for CO2-H2O interactions using ab initio quantum mechanical calculations.
- To assess the CO2 sequestration potential and methane recoverability in gas hydrate systems.
Main Methods:
- Computed the potential energy surface (PES) for H2O-CO2 interactions at the MP2/aug-cc-pVTZ level, correcting for basis set superposition error (BSSE) using the half counterpoise method.
- Developed the "VAS" model by fitting Exponential-6 and Lennard-Jones 6-12 potentials to the ab initio PES, accounting for many-body interactions.
- Calculated reference parameters for structure I carbon dioxide hydrate and predicted equilibrium pressure and cage occupancies using the VAS model. Validated using molecular dynamics simulations for CO2 density and diffusion.
Main Results:
- The VAS model accurately predicted pure CO2 hydrate equilibrium pressure with <3.2% average absolute deviation from experimental data.
- Predicted small cage occupancy for CO2 hydrate ranged from 32% to 51%.
- Predicted large cage occupancy for CO2 hydrate exceeded 98%.
Conclusions:
- The VAS model, based on ab initio calculations, provides accurate intermolecular potentials for CO2-H2O systems.
- The developed model enables reliable prediction of CO2 cage occupancies in gas hydrates, crucial for sequestration assessments.
- This work advances the understanding of gas hydrate behavior for energy and environmental applications.
Related Concept Videos
Phase Diagrams
Real Gases: Effects of Intermolecular Forces and Molecular Volume Deriving Van der Waals Equation
Calculating the Equilibrium Constant
For example, gaseous nitrogen dioxide forms dinitrogen tetroxide according to this equation:
Polyprotic Acids
Homogeneous Equilibria for Gaseous Reactions
For gas-phase reactions, the equilibrium constant may be expressed in terms of either the molar concentrations (Kc) or partial pressures (Kp) of the reactants and products. A relation between these two K values may be simply derived from the ideal gas equation and the definition of molarity. According to the ideal gas equation:
Van der Waals Equation
First, the attractive forces between molecules, which are stronger at higher densities and reduce the pressure, are considered by adding to the pressure a term equal to the square of the molar density multiplied by a positive coefficient a. Second, the...

