---
title: Thermal-Energy-Mediated Photoevaporation (TEMP)
url: https://www.emergentmind.com/topics/thermal-energy-mediated-photoevaporation-temp
type: topic
---

# Thermal-Energy-Mediated Photoevaporation (TEMP)

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 [2309.01832].

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$_2$O) [2408.00848, 1706.03155]. 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”
\[
R_g = \frac{GM_*}{c_s^2}
\]
where $c_s$ is the local sound speed set by $T$ [2211.10130, 1309.0508]. For disk surface layers, this translates to a launch threshold $T \gtrsim GM_*/(2kR)$ [2408.00848]. The local mass-flux at the sonic point is
\[
\dot\Sigma = \rho_{\rm base} c_s
\]
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 [2510.02201, 2605.02766].

## 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 ($\sim$10$^4$ K) ionized layer, while FUV (6–13.6 eV) and X-rays (0.1–10 keV) heat neutral/molecular regions to $\sim$10$^2$–10$^4$ K via photoelectric and photoionization channels [2211.10130, 2406.18461]. 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 [1706.03155, 2511.00515].

Steady-state axisymmetric simulations (e.g., PLUTO+PRIZMO, MOCASSIN, FRIED) show typical wind mass-loss rates of $10^{-10}$–$10^{-8}\,M_\odot\,\mathrm{yr}^{-1}$ for solar-type stars, with X-ray luminosity scaling $\dot{M}\propto L_X^{1.0\text{--}1.2}$ [2408.00848, 2211.10130, 1910.08565]. 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 [1112.1087, 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 [2310.20214]. 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, $\lambda_s = (G M_p/R_s)/(k T_s/\mu_s)$, drops below unity, i.e., enthalpy dominates gravitational binding at the sonic point [2510.02201]. 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 $\dot{M}_{\rm TEMP} \propto F_{\rm XUV}^{0.6} R_{\rm base}^2$ (where $R_{\rm base}$ is the photoionization base) [2510.02201].

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 [2510.02201, 2605.02766]. 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 [2308.00020, 2605.02766], 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 [2511.00515, 1706.03155, 2605.02766].
- **Thermochemistry**: Operator-split or fully implicit chemical networks evolve $\sim$25–150 species, accounting for photoionization, photodissociation, recombination, charge exchange, and advection [1706.03155, 2511.00515]. Non-equilibrium chemistry is essential, as advected H$_2$ and other molecules may survive far into the outflow and dominate cooling [2408.00848, 1706.03155].
- **Microphysics**: Efficient line cooling by O I excited by neutral H, collision partners, and molecular coolants (H$_2$, H$_2$O) sets the cooling rate. Neglecting atomic (O+H) cooling overestimates wind temperatures and the resulting mass flux (by up to an order of magnitude) [2408.00848].
- **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 [1910.08565].

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 [2309.01832, 1605.04066, 1105.0045]. 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 [2309.01832].
- **H$_2$ 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 $\dot{M} \sim 10^{-10}-10^{-8} M_\odot {\rm yr}^{-1}$ [2511.00515].
- **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 [2211.10130].
- **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") [1309.0508, 1112.1087].

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 [1605.04066, 1105.0045].

## 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 [2605.02766, 2308.00020]. 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$^7$) [1911.04510]. 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 [2401.15419].

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 [1910.08565]. 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 $\dot\Sigma_{\rm TEMP}(R)$ 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 [1112.1087, 1309.0508].

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 [2605.02766, 2510.02201].

## 7. Open Problems and Theoretical Frontiers

Despite its success, TEMP modeling faces ongoing challenges:

- Precise determination of heating efficiencies $\eta$ and spectral hardness dependencies, particularly under time-dependent irradiation and for different metallicities [2605.02766, 1910.08565].
- Consistent treatment of dust evolution, grain entrainment, and consequent feedback on radiative transfer and mass loss [2310.20214].
- Robust inclusion of non-equilibrium chemistry and advective processes shaping molecular survival and emission line diagnostics [2408.00848, 1706.03155].
- Joint modeling of TEMP and MHD-wind contributions—especially in the inner disk—remains necessary for full reconciliation with observational constraints [2401.15419, 1911.04510].
- 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 [2605.02766, 2310.20214].

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.

Source: https://www.emergentmind.com/topics/thermal-energy-mediated-photoevaporation-temp