Squeezed Coal Reveals Surprising Methane and Carbon Dioxide Behavior Under Stress

Deep underground, coal seams are never at rest. As mining operations push deeper into the Earth, the rock layers above press down on coal beds with enormous force, squeezing and compacting the molecular architecture of the coal itself. That squeezing does more than deform the rock: it fundamentally changes how methane and carbon dioxide—the two most consequential gases in the coalbed methane story—stick to and move through the coal matrix. A new molecular simulation study published in the Journal of Saudi Chemical Society has now mapped this hidden interplay in unprecedented atomic detail, and its findings could reshape how engineers approach gas extraction, mine safety, and carbon sequestration in stress-riddled coal reservoirs.

The research team, led by Xiaoqing Liu of Shanxi Datong University together with colleagues from Liaoning Technical University, constructed four molecular models of coal spanning nearly the entire coalification spectrum: lignite, long-flame coal, coking coal, and anthracite. These were not arbitrary cartoons of coal. Each model was built from real experimental data, including elemental analysis, X-ray photoelectron spectroscopy, and solid-state carbon-13 nuclear magnetic resonance spectroscopy, yielding molecular formulas such as C124H99O39N for the lignite and C150H65O11NS for the anthracite. The models captured how coalification progressively condenses aromatic ring structures, shortens aliphatic side chains, and strips away oxygen-containing functional groups, producing an increasingly ordered and rigid macromolecular framework as rank increases.

With these validated structures in hand, the researchers applied isotropic compressive stress to the model unit cells at levels ranging from 40 to 200 megapascals, chosen to reflect the in situ stresses encountered in deep coal seams and beyond. The compression produced exactly what one would expect geometrically: pore free volume, pore surface area, and porosity all fell, while the density of the models rose. But the gas behavior that followed was anything but straightforward. As stress increased, the adsorption of both methane and carbon dioxide traced a distinctive decrease–increase–decrease pattern, a nonlinear signature born from two competing effects acting simultaneously at the molecular scale.

The explanation lies in a tug-of-war between space and force. Compression shrinks the effective adsorption space available inside the coal, which should reduce how much gas the coal can hold. At the same time, however, squeezing the pore walls closer together enhances the overlap of the potential energy fields emanating from those walls, deepening the adsorption wells and strengthening the attraction between gas molecules and the coal surface. In certain stress ranges, this second effect briefly wins, producing a recovery in adsorption capacity before the loss of pore space ultimately dominates and drives adsorption down again. The net result across the full stress range was inhibitory: when stress climbed from zero to 0.20 gigapascals, saturated methane adsorption fell by as much as 43.75 percent in the lignite model, and carbon dioxide adsorption fell by similar margins.

One of the study’s most striking discoveries is that methane is more sensitive to stress than carbon dioxide. The reason is deceptively simple: methane molecules are physically larger. Their greater kinetic diameter makes them more vulnerable to being blocked in compressed pores, so when the coal matrix tightens, methane suffers disproportionately. Carbon dioxide, with its smaller molecular footprint and stronger electrostatic interactions with the polar functional groups in coal, consistently out-adsorbed methane under every condition tested. This asymmetry matters enormously for CO2-enhanced coalbed methane recovery, the technique in which injected carbon dioxide displaces methane from the coal while simultaneously locking away a greenhouse gas underground.

The energy analysis reinforced this hierarchy. Adsorbed gases showed unimodal energy distributions whose peaks shifted toward more negative values as stress increased, indicating that individual adsorption sites became more energetically favorable under compression. In the lignite model, for example, the peak adsorption energy of methane deepened from −5.45 to −5.95 kilocalories per mole across the stress range, while carbon dioxide peaked at values between −7.45 and −7.95 kilocalories per mole. Decomposing the interactions revealed that methane adsorption is almost entirely governed by van der Waals forces, accounting for more than 99 percent of the total interaction, whereas carbon dioxide’s interaction includes a substantial electrostatic component alongside van der Waals contributions of roughly 64 to 82 percent.

Diffusion told a parallel but distinct story. The self-diffusion coefficients of both gases fell steadily under stress, with methane dropping by 46 to 65 percent and carbon dioxide by 39 to 60 percent across the four coal models. Visualizing the diffusion trajectories made the mechanism vivid: under stress-free conditions, most gas molecules wandered between pores along irregular, streamline-like paths, colliding frequently with pore walls and covering large distances. Under 0.20 gigapascals of compression, however, the trajectories collapsed into tight clustered patterns, showing that molecules had become trapped oscillating within isolated pores. Stress, in effect, imprisons the gas in its nanoscale cages.

Coal rank added another layer of complexity. Adsorption capacity rose steadily with coalification, mirroring the growth in pore surface area and free volume in the higher-rank models. Diffusion, by contrast, followed a U-shaped trend across the rank spectrum, dipping in the middle-rank coals before recovering in anthracite. Carbon dioxide diffused faster than methane in every model despite its stronger binding to the coal, because its smaller kinetic diameter lets it slip through confined micropore channels more easily. The practical implication is that carbon dioxide can penetrate the coal matrix more rapidly at the molecular scale, giving it an edge in the competitive race to occupy adsorption sites and displace methane during enhanced recovery operations.

The simulations also probed the mechanical integrity of the gas-bearing coal itself. Bulk and shear moduli increased with both compressive stress and coal rank, while compressibility declined, reflecting the closure of compressible pore space and the rigid aromatic framework of mature coals. Anthracite emerged as the most mechanically robust and structurally stable of the four. Cohesive energy density, a molecular-scale gauge of how tightly the organic matter holds together, climbed with both stress and rank, and was consistently higher in carbon dioxide-bearing models than in methane-bearing ones, suggesting that carbon dioxide injection and loading both strengthen intermolecular cohesion within the simulated coal systems.

For a field racing to balance energy production with climate goals, these molecular-level insights arrive at an opportune moment. Coalbed methane is a lower-emission fossil fuel than coal itself, and CO2-enhanced recovery promises a dual dividend: boosting methane output while permanently sequestering carbon dioxide in the same seams. By showing precisely how stress reshapes the adsorption and diffusion landscape—and why methane and carbon dioxide respond so differently to the squeeze—the study gives reservoir engineers a firmer physical foundation for designing extraction and sequestration strategies in the deep, high-stress coal seams where the industry’s future increasingly lies.

Subject of Research: Molecular simulation of methane and carbon dioxide adsorption and diffusion in coals of different ranks under compressive stress loading

Article Title: Molecular simulation of adsorption and diffusion behavior of CH4/CO2 in coals of different ranks under stress loading

Article References: Liu, X., Liang, F., Xu, Q., Zhao, D., Bai, W., & Yang, X. (2026). Molecular simulation of adsorption and diffusion behavior of CH4/CO2 in coals of different ranks under stress loading. Journal of Saudi Chemical Society, 30(5), Article 73. https://doi.org/10.1007/s44442-026-00128-z

Image Credits: AI Generated

DOI: 10.1007/s44442-026-00128-z

Keywords: coalbed methane, molecular simulation, adsorption, diffusion, compressive stress, coal rank, carbon dioxide sequestration, CO2-ECBM, GCMC, molecular dynamics, coal macromolecular model, pore structure