Thermal-Energy-Mediated Photoevaporation (TEMP)

• TEMP is a thermally-driven atmospheric escape process where high-energy stellar photons heat the upper layers of planetary atmospheres and disks above the escape temperature.
• It employs photoionization, photoelectric effects, and detailed cooling mechanisms to launch transonic hydrodynamic winds with mass-loss rates typically ranging from 10⁻¹⁰ to 10⁻⁸ M☉/yr.
• TEMP plays a key role in protoplanetary disk dispersal and the evolution of exoplanet envelopes, offering observable signatures such as blue-shifted forbidden lines and radio free-free emission.

Thermal-Energy-Mediated Photoevaporation (TEMP) is the thermally-driven atmospheric escape process by which high-energy stellar photons (EUV, FUV, X-rays) heat the upper layers of a planetary atmosphere or protoplanetary disk above the local escape temperature, resulting in unbound outflows. TEMP provides the dominant mass-loss channel in the dispersal of protoplanetary disks and the erosion of volatile-rich exoplanet envelopes across multiple astrophysical regimes. The mechanism is governed by the balance between radiative heating (primarily via photoionization and photoelectric effects) and radiative plus advective cooling, setting the thermodynamic structure and mass flux of transonic hydrodynamic winds. TEMP is quantitatively distinct from both classical “energy-limited” escape and “recombination-limited” escape, especially in low-gravity environments and regimes with significant molecular cooling or non-equilibrium chemistry.

1. Fundamental Physics and Formalism

At the core of TEMP is the interaction between incident high-energy photons and surface or atmospheric gas, resulting in photoionization, dissociation, and subsequent rapid thermalization of photoelectrons and secondary electrons. For EUV and X-rays, the local volumetric heating rate is typically cast as

$$\Gamma_X = n_H \sigma_X \frac{L_X}{4\pi r^2}\, e^{-\tau_X} E_{\rm heat}$$

where $n_H$ is the hydrogen nuclei number density, $\sigma_X$ is the mean absorption cross-section, $L_X$ is the X-ray luminosity, $\tau_X$ is the line-of-sight X-ray optical depth, and $E_{\rm heat}$ is the mean deposited energy per photoionization (Rab et al., 2023).

The equilibrium between heating ($\Gamma$) and cooling ($\Lambda$) determines the gas temperature profile. Dominant cooling channels include collisionally excited forbidden-line cooling (O I, C II) by hydrogen and electrons (with rate coefficients $k_{ij}^H(T)$), as well as rovibrational molecular cooling (primarily H$_2$ and H$n_H$0O) (Sellek et al., 2024, Wang et al., 2017). In many protoplanetary disk models, additional adiabatic (PdV) losses and dust–gas thermal accommodation must also be included.

The wind is launched when the thermal pressure gradient becomes sufficient to drive the gas through a sonic transition, analogous to the Parker wind. The critical launching radius for a purely thermal (isothermal) flow is the “gravitational radius”

$n_H$1

where $n_H$2 is the local sound speed set by $n_H$3 (Ercolano et al., 2022, Owen et al., 2013). For disk surface layers, this translates to a launch threshold $n_H$4 (Sellek et al., 2024). The local mass-flux at the sonic point is

$n_H$5

with the integrated mass-loss rate obtained by surface integration over the wind-launching region.

Classical energy-limited and recombination-limited regimes assume either all photon power is converted to gravitational lift (EL) or that strong ionization equilibrium thermostats the wind (RL). In contrast, the TEMP regime arises when a significant fraction of input energy is converted into enthalpy (thermal and kinetic), the outflow is not strongly recombination-limited, and the flow temperature profile is strongly shaped by microphysical and chemical cooling (Tang et al., 2 Oct 2025, Misener et al., 4 May 2026).

2. TEMP across Astrophysical Contexts

Protoplanetary Disks

In protoplanetary disks, TEMP is the principal driver of gas dispersal. EUV (13.6–100 eV) creates a hot (%%%%13$n_H$214%%%% K) ionized layer, while FUV (6–13.6 eV) and X-rays (0.1–10 keV) heat neutral/molecular regions to %%%%15$n_H$216%%%%–10$\sigma_X$0 K via photoelectric and photoionization channels (Ercolano et al., 2022, Nakatani et al., 2024). The penetration depths and corresponding thermal structures result in a multi-layer disk atmosphere: a cold midplane, a warm FUV/X-ray heated molecular layer, and a hot wind (Wang et al., 2017, Nakatani et al., 1 Nov 2025).

Steady-state axisymmetric simulations (e.g., PLUTO+PRIZMO, MOCASSIN, FRIED) show typical wind mass-loss rates of $\sigma_X$1–$\sigma_X$2 for solar-type stars, with X-ray luminosity scaling $\sigma_X$3 (Sellek et al., 2024, Ercolano et al., 2022, Wölfer et al., 2019). The mass-loss surface density profiles peak at 1–10 AU, declining as power laws at larger radii, with total dispersal times constrained by the competition between wind and viscous evolution (Owen et al., 2011, 1904.02752).

Dust physics is critical: the dust component is entrained in the wind and can modulate FUV opacity, creating feedback that can self-limit disk loss rates in regions with significant small-grain content (Gárate et al., 2023). The presence of dust traps/substructures can extend the longevity of the dust reservoir inside the wind truncation radius, but cannot prevent rapid dispersal in strong FUV environments.

Exoplanet Atmospheres

In low-mass exoplanets (sub-Neptunes, super-puffs), the TEMP regime arises where the escape parameter at the sonic point, $\sigma_X$4, drops below unity, i.e., enthalpy dominates gravitational binding at the sonic point (Tang et al., 2 Oct 2025). Here EUV and to a lesser extent X-ray photons heat the wind base, but the outflow is governed by the conversion of energy into enthalpy and then kinetic energy, with mass-loss rates scaling as $\sigma_X$5 (where $\sigma_X$6 is the photoionization base) (Tang et al., 2 Oct 2025).

TEMP governs the evolutionary path of close-in exoplanets, explaining observed bimodal radius distributions and the long-term retention of volatiles in “super-puff” planets at low to intermediate irradiation levels (Tang et al., 2 Oct 2025, Misener et al., 4 May 2026). The transition between photoevaporation and alternative (core-powered) escape regimes depends on the location of the Bondi radius with respect to the XUV penetration depth (Owen et al., 2023, Misener et al., 4 May 2026), with planets potentially transitioning through TEMP as they contract and cool.

3. Multi-band Radiation Hydrodynamics and Key Microphysics

Modern TEMP modeling employs coupled radiation-hydrodynamic and thermochemical solvers:

• Hydrodynamics: The conservative Euler equations are solved, including mass, momentum, and energy conservation.
• Radiative transfer: Multi-band (EUV/FUV/X-ray/IR) radiative transfer is implemented with ray-tracing along 1D (spherical, cylindrical) or multidimensional grids, with attenuation by gas and dust (Nakatani et al., 1 Nov 2025, Wang et al., 2017, Misener et al., 4 May 2026).
• Thermochemistry: Operator-split or fully implicit chemical networks evolve $\sigma_X$725–150 species, accounting for photoionization, photodissociation, recombination, charge exchange, and advection (Wang et al., 2017, Nakatani et al., 1 Nov 2025). Non-equilibrium chemistry is essential, as advected H$\sigma_X$8 and other molecules may survive far into the outflow and dominate cooling (Sellek et al., 2024, Wang et al., 2017).
• Microphysics: Efficient line cooling by O I excited by neutral H, collision partners, and molecular coolants (H$\sigma_X$9, H$L_X$0O) sets the cooling rate. Neglecting atomic (O+H) cooling overestimates wind temperatures and the resulting mass flux (by up to an order of magnitude) (Sellek et al., 2024).
• Boundary conditions and feedback: Disk structure, dust content, and metallicity (especially C and O abundances) affect photon penetration and thermal structure, with C-depletion enabling deeper heating and enhanced wind mass fluxes (Wölfer et al., 2019).

Key parameters are compiled in tabulated "grids" (e.g., FRIED, Aiolos) allowing efficient interpolation during population synthesis or coupled evolution calculations.

4. Observational Diagnostics and Empirical Constraints

TEMP winds are empirically traced via low-velocity, blue-shifted forbidden lines, thermal free-free emission, and far-infrared or mm-wave continuum:

• [O I] 6300 Å: Collisionally excited in the hottest neutral regions of the wind, peaking at T~8,000–10,000 K and tracing the inner few AU of the disk. Observed line centroids are mildly blue-shifted (–0.8 km/s), with FWHM ≈ 10 km/s (Rab et al., 2023, Ercolano et al., 2016, Pascucci et al., 2011). Emissivity is sharply peaked in regions of steep temperature gradients.
• [Ne II] 12.8 μm: Traces higher altitude/ionized regions, broader and more blue-shifted than [O I], allowing spatial and kinematic decomposition of the wind structure (Rab et al., 2023).
• H$L_X$1 Pure Rotational Lines: Imaging with JWST/MIRI reveals characteristic "X-shaped" morphologies, with opening angles ≈ 37–50°, matching predictions from thermal wind models and inferring mass-loss rates $L_X$2 (Nakatani et al., 1 Nov 2025).
• Continuum Emission: Wind-driven free-free emission in the radio can be used to quantitate total gas mass loss, offering a distinction from MHD-wind scenarios (Ercolano et al., 2022).
• Population-level trends: Mass-loss rates correlate with X-ray luminosity, and TEMP predicts inside-out clearing leading to transition disks with a short-lived, non-accreting final phase ("thermal sweeping") (Owen et al., 2013, Owen et al., 2011).

Interpretation of emission lines (e.g., [O I] 6300 Å) as wind diagnostics requires careful consideration of their temperature sensitivity, optical thickness, and local density structure (Ercolano et al., 2016, Pascucci et al., 2011).

5. Limits, Transitions, and Comparative Regimes

TEMP does not universally dominate across all environments. Core-powered mass loss can control atmospheric escape in low-gravity, highly irradiated planets before the XUV-penetration depth drops below the Bondi (or sonic) radius (Misener et al., 4 May 2026, Owen et al., 2023). The TEMP regime typically yields higher mass fluxes than pure core-powered escape, but lower than pure energy-limited photoevaporation, and represents an intermediate state where both heating channels contribute.

In protoplanetary disks, magneto-thermal (MHD) winds can coexist with or dominate thermal winds in regions of high field strength (midplane plasma β ≲ 10$L_X$3) (Rodenkirch et al., 2019). In inner disks (R ≲ 0.3 AU), purely thermally-driven TEMP winds are energetically disfavored: the high gas density and strong cooling preclude both mass-loss rates comparable to observations and significant [O I] emission, requiring alternative driving mechanisms such as magnetic launching (Lin et al., 2024).

The TEMP paradigm is also contextually sensitive to metallicity, dust abundance, and the details of stellar irradiation spectra. For instance, C depletion by a factor ∼10 can increase mass-loss by factors of up to ∼4 (Wölfer et al., 2019). Advection, dust evolution, and out-of-equilibrium chemistry further complicate direct application of analytic formulas.

6. TEMP in Disk Dispersal, Evolutionary Outcomes, and Population Synthesis

TEMP is integral to self-consistent disk evolution models, governing the timescales for disk clearing, the emergence of transition disks, and the atmospheric evolution of sub-Neptunes and super-puffs. Viscous evolution equations coupled to local $L_X$4 profiles yield inside-out disk erosion, rapid thermal sweeping of the outer disk once a critical surface density threshold is breached, and the ultimate dispersal of the planet-forming gas reservoir (Owen et al., 2011, Owen et al., 2013).

In population-level studies, TEMP-driven mass loss successfully explains the paucity of accreting large-hole transition disks, the observed radius valley in exoplanet demographics, and the long-term survival of low-density envelopes in super-puff planets within restricted mass–irradiation parameter space (Misener et al., 4 May 2026, Tang et al., 2 Oct 2025).

7. Open Problems and Theoretical Frontiers

Despite its success, TEMP modeling faces ongoing challenges:

• Precise determination of heating efficiencies $L_X$5 and spectral hardness dependencies, particularly under time-dependent irradiation and for different metallicities (Misener et al., 4 May 2026, Wölfer et al., 2019).
• Consistent treatment of dust evolution, grain entrainment, and consequent feedback on radiative transfer and mass loss (Gárate et al., 2023).
• Robust inclusion of non-equilibrium chemistry and advective processes shaping molecular survival and emission line diagnostics (Sellek et al., 2024, Wang et al., 2017).
• Joint modeling of TEMP and MHD-wind contributions—especially in the inner disk—remains necessary for full reconciliation with observational constraints (Lin et al., 2024, Rodenkirch et al., 2019).
• Extending population synthesis to variable environments, stellar ages, and metallicity demands precomputed grids of mass-loss rates (e.g., FRIED, Aiolos) coupled to evolving disk and planetary structures (Misener et al., 4 May 2026, Gárate et al., 2023).

TEMP thus constitutes a unifying, physically rigorous paradigm in the theory of disk dispersal and planetary atmosphere evolution, underpinned by radiation-hydrodynamic simulation, analytic scaling, and a growing set of empirical diagnostics across multiple astrophysical environments.