Ledoux Criterion

The Ledoux criterion is a thermodynamic and fluid-dynamical stability condition formulated by the Belgian astrophysicist Paul Ledoux in 1947 that governs the onset of macroscopic buoyant convection in chemically stratified stellar and planetary interiors. Extending Karl Schwarzschild's 1906 adiabatic stability criterion—which assumes a uniform chemical composition throughout the fluid—the Ledoux criterion incorporates the stabilizing or destabilizing buoyancy force generated by a radial gradient in mean molecular weight μ. Within astronomy and stellar evolution theory, a layer inside a hydrostatic star is dynamically stable against convection under the Ledoux criterion whenever the radiative temperature gradient ∇_rad remains smaller than the sum of the adiabatic temperature gradient ∇_ad and the composition gradient term (φ / δ) ∇_μ. Because nuclear fusion continuously converts lighter nuclei into heavier reaction products in stellar cores and shells, the choice between the Schwarzschild and Ledoux criteria—together with semiconvection and convective overshooting—controls the size of convective cores on the main sequence, the chemical profile of evolved red supergiants, and the final core compactness of supernova progenitors.

Displaced Fluid Parcel Analysis

Consider a spherically symmetric star in hydrostatic equilibrium, where the radial pressure gradient balances local gravitational acceleration g(r) = G m(r) / r²:

dP / dr = −ρ g

The local pressure scale height is defined as H_P ≡ −dr / dlnP = P / (ρ g). Suppose a small fluid parcel at radius r—characterized by ambient pressure P, temperature T, density ρ, and mean molecular weight μ—is displaced upward by an infinitesimal radial distance Δr > 0 into a new environment with ambient properties P + ΔP_env, T + ΔT_env, ρ + Δρ_env, and μ + Δμ_env.

If the displacement occurs rapidly compared to both the thermal diffusion timescale and the microscopic atomic diffusion timescale, the parcel's motion is adiabatic (entropy S is conserved) and composition-preserving (its internal mean molecular weight remains fixed at μ, so Δμ_parcel = 0). Furthermore, because the parcel adjusts to the surrounding pressure over the acoustic sound-crossing time—which is much shorter than the buoyant turnover time—pressure equilibrium is maintained continuously between the parcel and the ambient medium:

ΔP_parcel = ΔP_env = (dP / dr) Δr = −( P / H_P ) Δr

The parcel's density perturbation relative to the ambient medium after rising by Δr is Δρ = Δρ_parcel − Δρ_env. If Δρ > 0, the displaced parcel is denser than its new surroundings and experiences a downward restoring buoyancy force per unit volume f_buoy = −g Δρ < 0, returning it toward its equilibrium position (dynamical stability). Conversely, if Δρ < 0, the parcel is lighter than its surroundings and accelerates further upward, triggering macroscopic convective overturning.

General Equation of State Coefficients δ and φ

For a general stellar plasma whose equation of state ρ = ρ(P, T, μ) includes ideal gas pressure, radiation pressure, and partial electron degeneracy, a differential change in logarithmic density is expressed as:

dlnρ = α dlnP − δ dlnT + φ dlnμ

where the three dimensionless thermodynamic partial derivatives are defined as:

• Isothermal compressibility coefficient: α ≡ (∂lnρ / ∂lnP)_T,μ
• Isobaric thermal expansion coefficient: δ ≡ −(∂lnρ / ∂lnT)_P,μ
• Compositional density coefficient: φ ≡ (∂lnρ / ∂lnμ)_P,T

For a classical ideal gas without radiation pressure (P = ρ k_B T / [μ m_u]), these coefficients simplify to α = 1, δ = 1, and φ = 1. When radiation pressure P_rad = (1/3)*aT*⁴ is significant—as in high-mass O and B stars where β ≡ P_gas / P < 1—the coefficients become α = 1 / β, δ = (4 − 3β) / β, and φ = 1.

Applying this differential relation separately to the adiabatically displaced parcel (dlnT_parcel = ∇_ad dlnP, dlnμ_parcel = 0) and to the ambient stellar stratification (dlnT_env = ∇ dlnP, dlnμ_env = ∇_μ dlnP) under pressure equilibrium (dlnP_parcel = dlnP_env = −Δr / H_P) yields the net density contrast:

Δρ / ρ = (Δρ_parcel − Δρ_env) / ρ = [ −δ (∇_ad − ∇) + φ ∇_μ ] dlnP = [ δ (∇_ad − ∇) − φ ∇_μ ] (Δr / H_P)

Here the three logarithmic gradients with respect to pressure are:

1. Actual temperature gradient: ∇ ≡ dlnT / dlnP
2. Adiabatic temperature gradient: ∇_ad ≡ (∂lnT / ∂lnP)_S,μ
3. Mean molecular weight gradient: ∇_μ ≡ dlnμ / dlnP

Statement of the Ledoux Stability Criterion

For the layer to be dynamically stable against convection, an upward displacement (Δr > 0) must produce a positive density excess (Δρ > 0). Because H_P > 0 and δ > 0 in normal stellar plasma, dynamical stability requires:

∇ < ∇_ad + (φ / δ) ∇_μ ≡ ∇_L

In a static layer prior to the onset of convection, energy is transported outward solely by photon diffusion, so the actual temperature gradient equals the radiative gradient ∇ = ∇_rad = (3 κ L P) / (16π a c G m T⁴), where κ is the Rosseland mean opacity (often governed by Kramers' opacity law or electron scattering). Thus, a stellar zone is Ledoux-stable if:

∇_rad < ∇_ad + (φ / δ) ∇_μ

and Ledoux-unstable (undergoing full dynamical convection) if:

∇_rad > ∇_ad + (φ / δ) ∇_μ

When the chemical composition is spatially uniform (∇_μ = 0), the Ledoux gradient ∇_L reduces to ∇_ad, recovering the classical Schwarzschild criterion (∇_rad < ∇_ad for stability).

Thermal and Compositional Buoyancy Frequencies

The equation of motion for the radially displaced fluid parcel is ρ (d²Δr / dt²) = −g Δρ = −ρ N² Δr, where N is the Brunt–Väisälä (buoyancy) frequency. Substituting Δρ / ρ from the Ledoux derivation splits N² into a thermal stratification component N_T² and a chemical composition component N_μ²:

N² = N_T² + N_μ² = (g δ / H_P) (∇_ad − ∇) + (g φ / H_P) ∇_μ

When N² > 0 (Ledoux stability), a displaced fluid element oscillates vertically about its equilibrium radius at angular frequency N, supporting the propagation of internal gravity modes (g-modes). When N² < 0 (Ledoux instability), N is imaginary and small perturbations grow in amplitude on the dynamical convective turnover timescale τ_dyn ~ |N|⁻¹.

Classification of Double-Diffusive Regimes

Because thermal diffusivity K_T = 4 a c T³ / (3 κ ρ² c_P) in a stellar radiative interior exceeds microscopic ion diffusivity D_ion by six to ten orders of magnitude (giving an inverse Lewis number τ ≡ D_ion / K_T ~ 10⁻¹⁰ to 10⁻⁶), the interplay between ∇_rad − ∇_ad and ∇_μ defines four distinct hydrodynamic stability regimes:

Core Hydrogen and Helium Burning

During core hydrogen burning on the upper main sequence (M ≳ 1.5 M_☉), the convective core recedes in mass coordinate m(r) as hydrogen (μ ≈ 0.60) is fused into helium (μ ≈ 1.34), leaving behind a positive mean molecular weight gradient (∇_μ > 0, where μ increases inward along pressure P). In stars between roughly 10 and 30 M_☉, radiation pressure lowers ∇_ad while electron-scattering opacity creates a shallow pocket just outside the receding core where ∇_ad < ∇_rad < ∇_L. A stellar evolution code using the pure Schwarzschild criterion mixes this zone instantaneously, whereas a code enforcing the Ledoux criterion inhibits dynamical convection at the μ-barrier, producing a stepped semiconvection zone.

Blue-to-Red Supergiant Loops and Core Compactness

After core hydrogen exhaustion, intermediate- and high-mass stars cross the Hertzsprung–Russell diagram to ignite helium via the triple-alpha process. Under the Ledoux criterion, the sharp composition jump ∇_μ at the hydrogen-burning shell restricts the growth of the intermediate convective zone (ICZ), keeping the envelope more compact and favoring blue supergiant survival (such as Sk −69° 202, the blue supergiant progenitor of SN 1987A in the Large Magellanic Cloud). During advanced carbon, neon, oxygen, and silicon shell burning in M > 8 M_☉ stars, the Ledoux criterion determines whether adjacent burning shells merge or remain separated by composition barriers, shifting the final Chandrasekhar-like iron core compactness parameter ξ_2.5 and altering whether the star explodes as a Type II supernova leaving a neutron star or collapses into a black hole.

Why does a positive mean molecular weight gradient (∇_μ > 0) stabilize a star against convection?

Because pressure P increases inward toward the center of a star, a positive gradient ∇_μ ≡ dlnμ / dlnP > 0 means that heavier nuclei (higher μ, such as helium or carbon) reside deeper in the gravitational potential well beneath lighter nuclei (lower μ, such as hydrogen). When a deep, heavy fluid parcel is displaced upward into a layer of lower mean molecular weight at the same pressure and temperature, its higher μ makes it denser than the surrounding lighter gas, creating a downward restoring buoyancy force that resists convective overturning.

When do the Ledoux and Schwarzschild criteria give identical boundaries?

The two criteria yield identical convective boundaries whenever the medium is chemically homogeneous (∇_μ = 0), such as inside a zero-age main-sequence star prior to nuclear burning or inside a surface ionization convection zone where no nuclear composition gradient exists.

1. Ledoux, P. (1947). "Stellar Models with Convection and with Discontinuity of the Mean Molecular Weight." The Astrophysical Journal, 105, 305–321. https://doi.org/10.1086/144905
2. Schwarzschild, K. (1906). "Über das Gleichgewicht der Sonnenatmosphäre." Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1906, 41–53. https://ui.adsabs.harvard.edu/abs/1906NWGot...1...41S
3. Kippenhahn, R., Weigert, A., & Weiss, A. (2012). Stellar Structure and Evolution (2nd ed.). Springer. https://doi.org/10.1007/978-3-642-30304-3
4. Sukhbold, T., & Woosley, S. E. (2014). "The Compactness of Presupernova Cores." The Astrophysical Journal, 783(1), 10. https://doi.org/10.1088/0004-637X/783/1/10