Henyey Track

The Henyey track is a nearly horizontal evolutionary trajectory on the Hertzsprung-Russell diagram along which pre-main-sequence stars with masses exceeding approximately 0.5 M_⊙ undergo Kelvin–Helmholtz gravitational contraction at roughly constant (or slightly rising) bolometric luminosity while their effective surface temperature (T_eff) increases steadily toward the Zero-Age Main Sequence (ZAMS). Named after American astrophysicist Louis G. Henyey, who together with Robert LeLevier and R. D. Levée calculated the first numerical models of radiative pre-main-sequence stellar contraction in 1955 (and later developed the Henyey relaxation method for solving boundary-value differential equations of stellar structure in 1959), the Henyey track connects the bottom of the vertical convective Hayashi track to the onset of hydrostatic core hydrogen fusion.

Whereas low-mass protostars initially descend the Hayashi track with fully convective interiors governed by surface negative hydrogen ion (H⁻) opacity, contraction steadily heats their deep interiors (T_c ∝ M / R). In stars with M ≥ 0.5 M_⊙, rising interior temperatures cause Kramers bound-free and free-free radiative opacity (κ_rad ∝ ρ T^−3.5) to plummet near the center, allowing a stable radiative core to nucleate and expand outward. Once radiative diffusion controls the star's energy transport, the star makes a sharp ~90° bend to the left on the H–R diagram, transitioning onto the Henyey track.

Breakdown of Full Convection via Kramers Opacity Decline

According to the Schwarzschild criterion for convection, a stellar layer is stable against convective overturning and transports heat purely by radiative photon diffusion whenever the radiative temperature gradient ∇_rad is smaller than the adiabatic temperature gradient ∇_ad (≈ 0.4 for an ideal monatomic gas):

∇_rad = ( d ln T / d ln P )_rad = ( 3 / 16π a c G ) × ( κ L(r) P ) / ( M(r) T⁴ ) < ∇_ad

While a pre-main-sequence T Tauri star descends the Hayashi track, its interior is relatively cool (T_c ≈ 1 × 10⁶ to 3 × 10⁶ K) and dense, so Kramers radiative opacity—dominated by bound-free photoionization of heavy elements and free-free bremsstrahlung absorption—is large throughout the star. As gravitational contraction reduces the stellar radius R, the virial theorem dictates that the central temperature scales inversely with radius (T_c ∝ G μ m_u M / k_BR) while central density rises as ρ_c ∝ M / R³.

Substituting Kramers' opacity law (κ = κ₀ ρ T^−3.5) into the homology scaling for ∇_rad reveals how radiative transparency evolves during homologous contraction:

• As T_c climbs above ~3.5 × 10⁶ K, the steep T^−3.5 temperature drop in Kramers opacity overwhelms the ρ ∝ R⁻³ density increase, causing κ_c ∝ R^0.5 to decrease as the star shrinks.
• Consequently, ∇_rad drops below ∇_ad = 0.4 at the center of the star first. A stable radiative core forms at r = 0 and grows progressively outward in mass fraction q_rad = M_rad / M as contraction continues.
• When the radiative core encompasses approximately 50% to 70% of the star's total mass, the surface convective envelope thins sufficiently that the photospheric H⁻ valve no longer dictates the global luminosity. Luminosity ceases to fall, marking the Hayashi–Henyey transition point (the local luminosity minimum on the pre-main-sequence track).

Radiative Diffusion and Homology Relations

Once a contracting pre-main-sequence star is predominantly radiative throughout its interior, its bolometric luminosity L is determined by the rate at which photons diffuse down the interior temperature gradient across the radiative opacity barrier:

L = − ( 16π a c r² T³ ) / ( 3 κ ρ ) × ( dT / dr )

Applying standard stellar homology transformations (T ∝ μ M / R, ρ ∝ M / R³, dT/dr ∝ T / R, and Kramers opacity κ = κ₀ ρ T^−3.5) directly to the radiative diffusion equation yields the Henyey Luminosity Scaling Relation:

L ∝ ( R² T⁴ ) / ( κ ρ R ) ∝ ( R T^7.5 ) / ρ² ∝ M^5.5 R^−0.5 μ^7.5 κ₀⁻¹

This single analytical scaling explains all observational features of the Henyey track in astronomy:

1. Nearly Constant, Slightly Rising Luminosity (L ∝ R^−0.5): Because luminosity scales as the inverse square root of stellar radius (R^−0.5), a contraction of the star's radius by a factor of 2 increases its bolometric luminosity by a factor of 2^0.5 ≈ 1.41 (+0.15 dex), tracing a shallow, upward-sloping horizontal path toward the left on the H–R diagram.
2. Rapid Blueward Heating (T_eff ∝ R^−5/8): Combining L ∝ R^−0.5 with the Stefan–Boltzmann radiation law (L = 4π R² σ_SB T_eff⁴) shows that L ∝ M^22/5 T_eff^4/5, or equivalently log₁₀ L ≈ 0.8 log₁₀ T_eff + 4.4 log₁₀ M + constant. As the star contracts, its effective temperature rises steadily as T_eff ∝ R^−5/8.
3. Steep Mass Dependence (L ∝ M^5.5): Higher-mass pre-main-sequence stars cross the H–R diagram on much brighter Henyey tracks and traverse them in a fraction of the Kelvin–Helmholtz timescale (t_KH ∝ M² / R L ∝ M^−3.5).

Quantitative Comparison Across Stellar Mass Regimes

As a star of 1.0 to 3.0 M_⊙ nears the left end of the Henyey track, its central temperature approaches 1.0 × 10⁷ to 1.5 × 10⁷ K. Just before settling onto the Zero-Age Main Sequence (the stellar locus later bounded at high age by the main-sequence turnoff), the evolutionary track executes a characteristic small downward hook (the pre-ZAMS luminosity dip of Δ log L ≈ −0.10 to −0.15 dex):

1. Out-of-Equilibrium ¹²C(p,γ)¹³N(β⁺ν)¹³C(p,γ)¹⁴N Ignition: Even in a 1.0 M_⊙ star that will ultimately burn hydrogen primarily via the proton–proton chain, primordial carbon-12 (X(¹²C) ≈ 0.002) begins capturing protons to form nitrogen-14 via the first branch of the CNO cycle at T_c ≈ 1.1 × 10⁷ K.
2. Temporary Convective Core Expansion: Because the ¹²C + p reaction has a steep T¹⁶ temperature sensitivity, this brief burst of central nuclear energy creates a transient convective core that expands the central regions slightly.
3. Radiative Envelope Adjustment: The expansion of the deep core absorbs a portion of the outward heat flux, causing the surface luminosity L and radius R to dip downward for a few million years until primordial ¹²C is fully converted into ¹⁴N and steady-state p–p chain or full CNO-cycle equilibrium locks the star onto the Zero-Age Main Sequence.

Observationally, intermediate-mass stars (1.5–8 M_⊙) traversing the Henyey track are identified as Herbig Ae/Be stars—young A- and B-type emission-line stars still surrounded by remnant gas and dust in protoplanetary disks within star-forming regions such as the Orion Nebula, whose atmospheric parameters (T_eff, log g) are routinely measured using Strömgren photometry and high-resolution spectroscopy.

1. Henyey, L. G., LeLevier, R., & Levée, R. D. (1955). "The Early Phases of Stellar Evolution." Publications of the Astronomical Society of the Pacific, 67(396), 154–160. https://doi.org/10.1086/126791
2. Henyey, L. G., Wilets, L., Böhm, K. H., LeLevier, R., & Levée, R. D. (1959). "A Method for Automatic Computation of Stellar Evolution." The Astrophysical Journal, 129, 628–636. https://doi.org/10.1086/146661
3. Iben, I., Jr. (1965). "Stellar Evolution. I. The Approach to the Main Sequence." The Astrophysical Journal, 141, 993–1018. https://doi.org/10.1086/148193
4. Palla, F., & Stahler, S. W. (1993). "The Pre-Main-Sequence Evolution of Intermediate-Mass Stars." The Astrophysical Journal, 418, 414–425. https://doi.org/10.1086/173402