Kinetics of CH<sub>4</sub> and CO<sub>2</sub> Hydrate Dissociation and Gas Bubble Evolution via MD Simulation
Bibliographic record
Abstract
Molecular dynamics simulations of gas hydrate dissociation comparing the behavior of CH4 and CO2 hydrates are presented. These simulations were based on a structurally correct theoretical gas hydrate crystal, coexisting with water. The MD system was first initialized and stabilized via a thorough energy minimization, constant volume-temperature ensemble and constant volume-energy ensemble simulations before proceeding to constant pressure-temperature simulations for targeted dissociation pressure and temperature responses. Gas bubble evolution mechanisms are demonstrated as well as key investigative properties such as system volume, density, energy, mean square displacements of the guest molecules, radial distribution functions, H2O order parameter, and statistics of hydrogen bonds. These simulations have established the essential similarities between CH4 and CO2 hydrate dissociation. The limiting behaviors at lower temperature (no dissociation) and higher temperature (complete melting and formation of a gas bubble) have been illustrated for both hydrates. Due to the shift in the known hydrate stability curves between guest molecules caused by the choice of water model as noted by other authors, the intermediate behavior (e.g., 260 K) showed distinct differences however. Also, because of the more hydrogen-bonding capability of CO2 in water, as reflected in its molecular parameters, higher solubility of dissociated CO2 in water was observed with a consequence of a smaller size of gas bubble formation. Additionally, a novel method for analyzing hydrate dissociation based on H-bond breakage has been proposed and used to quantify the dissociation behaviors of both CH4 and CO2 hydrates. Activation energies Ea values from our MD studies were obtained and evaluated against several other published laboratory and MD values. Intrinsic rate constants were estimated and upscaled. A kinetic reaction model consistent with macroscale fitted kinetic models has been proposed to indicate the macroscopic consequences of this analysis.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.001 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.001 | 0.000 |
| Research integrity | 0.000 | 0.001 |
| Insufficient payload (model declined to judge) | 0.002 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".