Transport-Induced Disequilibrium Chemistry

Updated 2 July 2026

Transport-induced disequilibrium chemistry is the deviation of chemical abundances from local equilibrium due to mass transport processes occurring faster than chemical reactions.
It employs three-dimensional models that couple dynamical transport with chemical kinetics to determine quench levels and predict observable spectral fingerprints.
Observational diagnostics reveal marked enhancements in species like CO and CO₂, influencing interpretations of both planetary atmospheres and battery materials.

Transport-induced disequilibrium chemistry defines the persistent deviation of chemical abundances from local thermochemical equilibrium, driven by atmospheric or solid-state mass transport occurring on timescales shorter than those required for chemical interconversion reactions to re-equilibrate the local composition. The phenomenon is central to understanding the observable signatures, bulk composition, and evolution of planetary atmospheres, cold substellar objects, and even solid-state battery materials. The primary mechanisms—vertical mixing, horizontal advection, and (in solids) diffusion—"freeze in" the chemical profile characteristic of deeper or spatially distinct regions, thereby imparting observable spectral or compositional fingerprints that differ significantly from equilibrium predictions.

1. Governing Equations, Timescales, and Quench Criteria

Transport-induced disequilibrium chemistry is fundamentally characterized by the interplay between chemical and transport timescales. For a key interconversion, such as CH₄ + H₂O ⇌ CO + 3 H₂, the local chemical timescale for species X is

$\tau_{\rm chem}(X) \sim [n\langle k_{\rm forward} n_{\rm reactants} \rangle]^{-1},$

where $k_{\rm forward}$ is the rate constant of the rate-limiting step and $n$ the relevant local density. For hydrogen-rich atmospheres at $T\sim100$ – $1000\,\rm K$ , these timescales can be extremely long— $\sim10^{9}-10^{12}\,$ s above $\sim1\,$ bar depending on the reaction network and temperature, making atmospheric chemistry orders of magnitude slower than dynamical processes (Liu et al., 9 Apr 2026).

Vertical transport is parameterized by an eddy diffusion coefficient $K_{zz}$ (in units $\mathrm{cm}^2\,\mathrm{s}^{-1}$ ), yielding a mixing timescale

$\tau_{\rm mix} = \frac{H^2}{K_{zz}},$

with $k_{\rm forward}$ 0 the local atmospheric scale height ( $k_{\rm forward}$ 1–100 km for sub-Neptunes and gas giants).

In strongly irradiated hot-Jupiter atmospheres, horizontal mixing occurs on an advection timescale

$k_{\rm forward}$ 2

where $k_{\rm forward}$ 3 is a characteristic length scale (often the planetary radius) and $k_{\rm forward}$ 4 the mean zonal wind ( $k_{\rm forward}$ 5–3 km s $k_{\rm forward}$ 6).

The quench level ( $k_{\rm forward}$ 7, $k_{\rm forward}$ 8) is then defined by the balance

$k_{\rm forward}$ 9

depending on whether vertical or horizontal mixing dominates. Above this level, the abundance of each species remains fixed at its equilibrium value at the quench point, independent of the local temperature and pressure (Liu et al., 9 Apr 2026, Zamyatina et al., 2022, Miles et al., 2020).

2. Three-Dimensional Models and Parameter Estimation

Rigorous treatment of transport-induced disequilibrium chemistry requires coupling a three-dimensional (3D) general circulation model (GCM) to chemistry modules that solve the full advection–reaction–diffusion equations for each species. This system incorporates the following:

Full Navier–Stokes dynamics with planetary rotation (for winds and advection)
Radiative transfer with pressure- and temperature-dependent opacities for key absorbers (e.g., H₂O, CH₄, CO, CO₂)
Vertical and horizontal transport through resolved wind fields (upwelling, superrotation, polar vortices)
Treatment of chemistry, either explicitly via stiff-ODE solvers for comprehensive kinetic networks (e.g., over 100 species, nearly 2000 reactions), or through temperature-pressure dependent relaxation timescales calibrated from kinetics (Liu et al., 9 Apr 2026, Zamyatina et al., 2022, Mendonça et al., 2018)

A passive tracer, obeying

$n$ 0

is often introduced to empirically estimate $n$ 1 through flux-gradient or mixing-length methods: $n$ 2 with $n$ 3 the vertical velocity and angle brackets denoting horizontal averaging. These diagnostics provide the critical input parameters for 1D kinetic-transport and retrieval models (Liu et al., 9 Apr 2026).

3. Chemical Pathways, Quench Points, and Abundance Profiles

The dominant reactions affecting quenched species in H2-rich giant planet and sub-Neptune atmospheres are:

$n$ 4

At low enough temperatures and high metallicities, as in temperate sub-Neptunes (e.g., K2-18b), the chemical timescales for CO/CH₄, CO/CO₂, and NH₃/N₂ interconversion become so long ( $n$ 5– $n$ 6 s at 1 bar, $n$ 7 K) that vertical mixing can dramatically enhance the upper-atmosphere levels of CO and CO₂—by over 12 orders of magnitude relative to chemical equilibrium. Typical results for K2-18b show equilibrium abundance of CO₂ and CO at $n$ 8 at $n$ 9 bar, compared to $T\sim100$ 0 under disequilibrium, with quenching occurring near $T\sim100$ 1 bar (Liu et al., 9 Apr 2026).

For cool brown dwarfs and cold giant planets, disequilibrium CO is the principal tracer: CO mixing ratios may be enhanced to $T\sim100$ 2 compared to equilibrium values of $T\sim100$ 3. In these objects, the vertical mixing strengths inferred from observations ( $T\sim100$ 4– $T\sim100$ 5) sharply rise with decreasing temperature, reaching their convective limits in the coldest cases (Miles et al., 2020).

The impact of transport-induced disequilibrium on other minor species depends sensitively on each molecule's underlying reaction network and the temperature/pressure gradients. While CO and CO₂ can be strongly enhanced, NH₃ generally remains near or only slightly above equilibrium unless the N₂/NH₃ quench occurs deep in the atmosphere. Non-detection of predicted species like PH₃ (phosphine) in brown dwarf observations is interpreted as evidence of photochemical destruction, incomplete reaction networks, or sub-solar P/H ratios (Miles et al., 2020, Zhang et al., 2024).

4. Observational Signatures and Retrieval Diagnostics

Transport-induced quenching imprints distinctive features in both transmission and emission spectra. Key observables include:

Enhanced CO and CO₂ opacity in the 4.2–4.6 μm region (Spitzer IRAC 4.5 μm, JWST NIRSpec/MIRI)
CH₄ features at 1.65, 2.3, and 3.3 μm, sensitive to both vertical and horizontal mixing
NH₃ signatures at 10.5 μm (JWST MIRI), especially if quenched abundances are elevated (Liu et al., 9 Apr 2026, Blumenthal et al., 2017, Zamyatina et al., 2022)

The magnitude of disequilibrium effects depends on planetary properties and observing geometry. For example, the effect is maximized ("sweet spot") when the quench level for the target molecule falls within the spectral photosphere (e.g., $T\sim100$ 6–1 bar for CO/CH₄ in warm Jupiters with $T\sim100$ 7–1400 K). In these cases, retrieval models predict emission/transmission differences of up to $T\sim100$ 8 ppm in 4–5 μm bands (detectable with JWST/NIRSpec) (Blumenthal et al., 2017, Baxter et al., 2021).

Empirical $T\sim100$ 9 values and their pressure dependence, retrieved from M-band (CO) spectroscopy in brown dwarfs, or from phase-resolved and multi-band transit/eclipse spectra in exoplanets, are essential for quantifying vertical mixing rates and constraining formation/evolution models (Miles et al., 2020, Bardet et al., 15 Jun 2025, Kawashima et al., 2021).

5. Horizontal vs. Vertical Transport: Three-Dimensional Effects

Beyond vertical mixing, horizontal transport is now observationally established as a driver of disequilibrium chemistry. In hot-Jupiter atmospheres, rapid zonal (day-to-night) winds can move parcels across global temperature gradients on timescales ( $1000\,\rm K$ 0– $1000\,\rm K$ 1 s) much shorter than the local chemical timescale on the nightside ( $1000\,\rm K$ 2– $1000\,\rm K$ 3 s). This "horizontal quenching" leads to day- and nightside chemical homogeneity in species like CO and H₂O, and suppresses the predicted nightside CH₄ that would arise in equilibrium (Parmentier et al., 6 May 2026, Mendonça et al., 2018). In NGTS-10Ab, JWST/NIRSpec observations confirm nightside CO and H₂O abundances identical to the dayside, with no detectable CH₄, uniquely attributable to horizontal transport-induced disequilibrium because vertical mixing and metallicity effects cannot independently explain this depletion (Parmentier et al., 6 May 2026).

Three-dimensional GCMs are required to capture spatially heterogeneous quenching, including jet-driven longitudinal mixing and polar vortex isolation, which can lead to localized abundance anomalies of species (e.g., CH₄ at high latitudes). While vertical and latitudinal mixing are secondary to jet-induced longitudinal effects in many cases, their contribution may be observable in planets with strong latitudinal temperature gradients or seasonal variability (Mendonça et al., 2018).

6. Transport-Induced Disequilibrium in Solid-State Materials

Although the archetype is atmospheric chemistry, transport-induced disequilibrium analogues arise in solid-state systems, such as phase-transforming electrodes or multiphase alloys. Diffusing species (e.g., Li⁺, Na⁺, K⁺) can accumulate nonuniformly, driving the system into local chemical potential disequilibrium, activating phase reactions even in the presence of kinetic barriers imposed by elastic or structural constraints (Salvadori et al., 18 Dec 2025, Cattermull et al., 4 Sep 2025). The relevant mass balance for a species $1000\,\rm K$ 4 is

$1000\,\rm K$ 5

with flux $1000\,\rm K$ 6 driven by gradients in chemical potential, coupled to stress and swelling.

In hybrid battery materials, e.g., K $1000\,\rm K$ 7Mn[Fe(CN) $1000\,\rm K$ 8], coupling between composition, strain, and ionic mobility gives rise to metastable, transport-controlled multiphase transitions inaccessible at equilibrium. Macroscopically, this results in rate-dependent voltage plateaux, compositional inhomogeneities, and phase-boundary trapping effects, underscoring the broader relevance of transport-induced disequilibrium chemical phenomena (Cattermull et al., 4 Sep 2025).

7. Analytical, Numerical, and Retrieval Methodologies

State-of-the-art retrieval and modeling approaches implement either full chemical kinetics (solving the coupled system of continuity, advection, and reaction equations) or adopt "chemical relaxation" approximations calibrated with detailed balance and laboratory rate constants to bypass computationally expensive network integration (Mendonça et al., 2018, Kawashima et al., 2021, Blumenthal et al., 2017). Model choices must be matched to data quality: JWST's wide bandpass and high S/N enable constraints on $1000\,\rm K$ 9, C/O, and [M/H] to levels that expose model assumptions and highlight the necessity of including disequilibrium in all forward/exoplanet retrievals (Bardet et al., 15 Jun 2025, Baxter et al., 2021). Neglecting disequilibrium transport systematically biases inferred atmospheric metallicities, elemental ratios, and hence planetary formation histories.

In the atmospheric context, the degeneracy between vertical and horizontal transport, interior temperature, and photochemistry must be accounted for when interpreting molecular abundances and their spectral signatures. Laboratory kinetics, 3D circulation models, and multi-wavelength retrievals are required for robust inference of the quenching process (Liu et al., 9 Apr 2026, Zamyatina et al., 2022).

References