Atmospheric Evolution Models

• Atmospheric evolution models are computational frameworks that simulate the temporal changes in composition, structure, and mass of stellar and planetary atmospheres.
• They integrate hydrodynamics, radiative transfer, chemistry, and escape mechanisms to predict atmospheric behavior across diverse astrophysical contexts.
• These models guide interpretations of observational data for stars, exoplanets, terrestrial planets, and substellar objects by linking physical processes with measurable outcomes.

Atmospheric evolution models are computational and theoretical frameworks that describe the changing physical, chemical, and structural properties of an atmosphere over evolutionary timescales. These models unify constraints from hydrodynamics, radiative transfer, chemistry, planetary/interior processes, external irradiation, and mass-loss (escape) to predict or interpret the time-dependent states of atmospheres in diverse astrophysical and planetary contexts, from terrestrial worlds and exoplanets to substellar objects and very massive stars.

1. Physical Scope, Goals, and Model Architectures

Atmospheric evolution models aim to capture the temporal progression of bulk atmospheric composition, structure, radiative properties, and mass for a body as a function of time. Formally, the models solve coupled equations for mass and energy conservation in the atmosphere, taking into account relevant sources and sinks (e.g., outgassing, escape, condensation, mixing, surface exchange, winds, photochemistry, and hydrodynamics).

The following generic form representing atmospheric inventory for volatile $i$ is characteristic: $$\frac{d N_i}{dt} = \Phi_{\mathrm{source}, i} - \Phi_{\mathrm{escape}, i} - \Phi_{\mathrm{loss}, i}$$ where $N_i$ is the total number of molecules or mass of species $i$ in the atmosphere, and the $\Phi$ terms represent source (outgassing, delivery, surface exchange) and various loss (escape, sequestration, chemical destruction) processes (Avice et al., 2020, Sinclair et al., 2020).

Model architectures fall into:

• Coupled atmosphere–interior models: Explicit feedbacks between atmospheric escape, outgassing, planetary cooling, redox evolution (e.g., PACMAN, Lichtenberg et al. frameworks) (Krissansen-Totton et al., 2022, Lichtenberg et al., 2021).
• Stellar/planetary structure + non-LTE/static atmosphere hybrids: Stellar evolution codes (e.g., GENEC) matched to detailed atmosphere codes (e.g., PoWR) at various optical depths (Josiek et al., 11 Apr 2025).
• Escape-driven evolutionary grids: 1D or 3D hydrodynamic codes constructing grids of atmospheric escape and contraction histories (often for exoplanets) (Kubyshkina et al., 2018).
• Atmosphere-interior cooling tracks for substellar/giant planets: Cooling, contraction, and spectrum models (e.g., ATMO, Sonora, Diamondback) (Phillips et al., 2020, Marley et al., 2021, Morley et al., 2024).
• Statistical/stochastic frameworks: Atmospheric erosion and growth under bombardment by planetesimals, with volatile delivery and stripping (Sinclair et al., 2020).

2. Governing Equations and Core Physical Processes

Atmospheric evolution models are constructed around a core set of conservation equations and physical processes:

A. Mass and Momentum Conservation (Hydrodynamics)

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$

$$\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v}\mathbf{v}) = -\nabla p + \rho \mathbf{g} + \text{external/stellar forces}$$

as in escape and upper atmosphere models (Kubyshkina et al., 2018, Hazra, 25 Feb 2025).

B. Energy Conservation (Radiative-Convective, Non-LTE)

Models implement radiative transfer (LTE/non-LTE), convection (mixing-length, MLT), and detailed heating/cooling: $$F_{\text{rad}} + F_{\text{conv}} = \sigma T_{\text{eff}}^4$$ Numerical radiative transfer is solved using multi-band correlated-$k$, discrete ordinate, or two-stream methods (Marley et al., 2021, Morley et al., 2024, Phillips et al., 2020). In hot star and hot Jupiter models, non-LTE and line-driven opacities are solved in the comoving frame (Josiek et al., 11 Apr 2025).

C. Atmospheric Escape:

Major loss channels modeled include:

• Thermal (Jeans) escape: Kinetic evaporation from the exobase where

$$\Phi_i = \left(\frac{n_i u_{\text{th}}}{2\sqrt{\pi}}\right) (1 + \lambda) \exp(-\lambda)$$

$\lambda = \frac{GM_pm}{R_{\text{exo}}kT}$

• Hydrodynamic escape (Parker wind, energy-limited flow):

$$\dot M_\mathrm{EL} = \frac{\eta \pi R_p^3 F_{\mathrm{XUV}}}{G M_p K}$$

with $K$ a Roche correction [(Hazra, 25 Feb 2025); (Kurokawa et al., 2013); (Chakrabarty et al., 2023)].

• Non-thermal escape: Impact-driven stripping, sputtering, photochemistry, modeled via efficiency factors or detailed reaction grids (Sinclair et al., 2020, Hazra, 25 Feb 2025).

D. Chemistry and Outgassing:

Chemistry includes equilibrium/rainout chemistry (Gibbs minimization), non-equilibrium kinetics, photochemistry, condensation, and outgassing. Evolution models typically solve for chemical networks with hundreds of reactions and dozens of species/phases (Johnstone et al., 2018, Lichtenberg et al., 2021, Phillips et al., 2020).

E. Interface to Interiors and Feedbacks:

Planetary cooling, internal outgassing, volatile recycling, and thermochemical feedbacks are incorporated either as parameterized functions or dynamically through time-dependent energy and mass fluxes (Krissansen-Totton et al., 2022, Phillips et al., 2020).

3. Model Coupling and Approximations: Atmosphere–Structure Interfaces

A critical issue in atmospheric evolution, particularly for stars and gas planets, is the interface between interior structure and atmosphere:

• Stellar models (VMS): GENEC (hydrostatic, gray atmosphere) is typically coupled to PoWR (non-LTE, expanding, radiation-driven wind) at a selected optical depth or density floor. The boundary choice (e.g., above vs. below iron opacity bump) leads to divergent predictions for $T_{\rm eff}$, radius, and spectra, with differences in $T_{\rm eff}$ reaching up to $20\,{\rm kK}$ (Josiek et al., 11 Apr 2025).
• Substellar objects: Interior adiabat is matched at a defined pressure (e.g. 10 bar) to the atmosphere's $T(P)$, and the correct $T_{\rm eff}$ found by solving for flux agreement iteratively. Updated hydrogen–helium EOS, opacities, and vertical mixing modify evolutionary tracks (Marley et al., 2021, Morley et al., 2024, Phillips et al., 2020).

The table below summarizes interface choices for very massive stars as an example:

Connection Point	Structural Feature	Effect on $T_{\rm eff}$, Radius
Shallow ($\tau=20$)	Above iron opacity bump	Matches inflated GENEC photosphere, cooler $T_{\rm eff}$
Deep ($\rho=10^{-7}$g/cm$^3$)	Below inflation zone	Excludes inflated plateau, hotter and more constant $T_{\rm eff}$

Neither interface fully captures the effects of dynamical instability or turbulent pressure in the iron opacity zone, indicating a need for simulations including 3D turbulence/hydrodynamics (Josiek et al., 11 Apr 2025).

4. Key Applications Across Astrophysical and Planetary Regimes

Atmospheric evolution models underpin quantitative predictions and interpretations in a range of scenarios:

A. Stellar Evolution (Very Massive Stars):

• Wind-driven mass-loss, iron-bump inflation, and the influence of atmosphere–structure coupling on predicted $T_{\rm eff}$ and mass-loss rates (Josiek et al., 11 Apr 2025).

B. Exoplanets (Mass Loss, Water Worlds):

• Atmospheric escape grid models (hydrodynamic, energy-limited, impact stripping) are used to reproduce the planetary radius valley, infer hidden populations of water-rich worlds, and identify observability windows for transmission and RV spectroscopy (Chakrabarty et al., 2023, Hazra, 25 Feb 2025, Kubyshkina et al., 2018).
• Model–data comparisons (e.g., migration models vs. observed period-radius distributions) guide target selection for the detection of bare water worlds.

C. Terrestrial Planets (Earth, Mars, Venus):

• Geological proxies (noble-gas/nitrogen isotopes), stochastic impact models, and coupled climate/escape/outgassing codes constrain the evolution of atmospheric mass and composition since planetary formation (Sinclair et al., 2020, Avice et al., 2020, Johnstone et al., 2018).
• Fully coupled models (e.g., PACMAN) for TRAPPIST-1 planets predict the present-day prevalence of CO$_2$ or O$_2$ atmospheres, hydrodynamic escape-driven erosion, and redox budget trajectories (Krissansen-Totton et al., 2022).

D. Substellar Objects (Brown Dwarfs, Giant Exoplanets):

• Cooling and contraction evolution is governed by atmospheric boundary conditions, opacities, and EOS. Evolutionary model grids (ATMO, Sonora, Diamondback) deliver synthetic spectra, photometry, and $T_{\rm eff}$–age–mass tracks, with updates from non-equilibrium chemistry or cloud physics causing $\sim$100–200 K changes in $T_{\rm eff}$ at a given age (Phillips et al., 2020, Marley et al., 2021, Morley et al., 2024).

E. Hot Jupiters (Inflation, Ohmic Dissipation):

• Ohmic dissipation models include coupling between atmospheric wind-driven currents and evolving internal dynamo fields, reproducing the observed inflated radii and predicting time-variable Ohmic efficiency and self-regulation via dynamo–convection feedback (Viganò et al., 18 Jul 2025).

5. Uncertainties, Benchmarks, and Model Limitations

Although atmospheric evolution models are built on first principles, they are subject to substantial uncertainties:

• Interface ambiguity: For massive-star models, the depth where the atmosphere is cut onto the structure model introduces $\sim$0.2 dex uncertainties in HRD position, systematic offsets in spectroscopic mass-loss estimates, and possibly incorrect predictions of wind regime or spectral type (Josiek et al., 11 Apr 2025).
• Escape efficiency regimes: For exoplanets, uncertainties in XUV flux evolution, energy-limited vs. recombination-limited escape, and geometric/stellar wind corrections propagate into large differences in predicted mass-loss and planet survival [(Hazra, 25 Feb 2025); (Chakrabarty et al., 2023); (Kurokawa et al., 2013)].
• Impact stochasticity: For early Earth, stochastic delivery of volatiles via large asteroids dominates over smooth trends, leading to order-of-magnitude scatter in atmosphere mass and composition (Sinclair et al., 2020).
• Chemical & radiative data: Updates to molecular line lists, cloud models, and EOS cause order-100 K effects in substellar cooling and atmospheric structure (Phillips et al., 2020, Morley et al., 2024).

Models are subject to limitations arising from:

• Dimensionality (1D vs. 3D missing turbulence or instabilities)
• Missing physics (e.g., turbulence, MHD, cloud microphysics, dynamical convection)
• Simplified or uncertain initial conditions (e.g., primordial volatile inventory, planet formation/migration history)
• Incomplete feedbacks (e.g., static or uncoupled wind corrections)
• Observational ambiguities (mass–composition–radius degeneracy)

6. Research Outlook and Future Directions

The ongoing development of atmospheric evolution models is driven by demands for increased realism and predictive power:

• Dynamically consistent coupled models: Incorporating 3D radiative–hydrodynamic turbulence, time-dependent non-LTE, and dynamically evolved interfaces is critical for resolving the discrepancies in $T_{\rm eff}$, radius, and wind structure in very massive stars (Josiek et al., 11 Apr 2025).
• Next-generation escape models: Improved MHD + radiative transfer simulations are required to robustly parameterize atmospheric loss in exoplanet evolution and to interpret observed transits and demographic structures (Hazra, 25 Feb 2025).
• Long-term coupled interior–atmosphere evolution: Fully self-consistent models integrating all key volatile cycling, climate, planetary redox, and escape mechanisms are now being benchmarked against planet samples (e.g., TRAPPIST-1, water worlds) and will be essential for interpreting JWST and ELT observations (Krissansen-Totton et al., 2022, Chakrabarty et al., 2023).
• Expanded opacity, cloud, EOS grids: For substellar objects, expansion of metallicity, C/O, and cloud parameter spaces, as well as improved opacities, will refine $T_{\rm eff}$–age–mass relations and synthetic spectra (Morley et al., 2024, Marley et al., 2021, Phillips et al., 2020).
• Empirical validation: Model predictions are being tested against precision RV, transmission spectra, planetary photometry, and exospheric escape diagnostics (e.g., Ly-$\alpha$, He 10830 Å), providing direct feedback for improving physical assumptions (Chakrabarty et al., 2023, Hazra, 25 Feb 2025).

Emerging observational capabilities and high-fidelity computation are rapidly advancing the use of atmospheric evolution models as indispensable tools for interpreting planetary, stellar, and substellar atmospheres across the galaxy.