Hayashi Track

The Hayashi track is a steeply sloped, nearly vertical luminosity–temperature locus on the Hertzsprung-Russell diagram traversed downward in luminosity by fully convective pre-main-sequence protostars (M ≲ 3 M_☉, observed as T Tauri stars) during quasi-hydrostatic Kelvin-Helmholtz gravitational contraction, and traversed upward in luminosity by post-main-sequence stars with deep convective envelopes ascending the Red Giant Branch and asymptotic giant branch. Discovered in 1961 by Japanese astrophysicist Chushiro Hayashi at Kyoto University, the track marks the hot boundary of the Hayashi forbidden zone—a thermodynamic domain on the cool right-hand side of the H–R diagram (T_eff ≲ 3,000 to 4,500 K, depending on stellar mass and metallicity) inside which no stable hydrostatic stellar configuration can exist because the photospheric boundary condition mediated by negative hydrogen ion (H⁻) opacity forces the interior entropy gradient to become superadiabatic.

Prior to Hayashi's 1961 theorem, astrophysicists following Louis Henyey, Robert LeLevier, and R. D. Levée (1955) had assumed that contracting protostars were predominantly radiative throughout their interiors, implying that young stars approached the main sequence slowly along a nearly horizontal low-luminosity radiative trajectory. Hayashi proved instead that young low-mass stars emerge onto the optical H–R diagram at luminosities 10 to 100 times brighter than their final main-sequence values and contract rapidly downward along a vertical convective track (L decreasing as R² at roughly constant T_eff ≈ 3,500–4,500 K), solving the long-standing discrepancy between theoretical models and the high observed luminosities of T Tauri stars in young open clusters such as NGC 2264 and the Orion Nebula.

Fully Convective Polytropes of Index n = 1.5

When a dense molecular cloud core (such as a Bok globule) collapses via the Jeans instability, dissociates molecular hydrogen (H₂ → 2 H at T ≈ 2,000 K), and ionizes atomic hydrogen (H → H⁺ + e⁻ at T ≈ 10⁴ K), it settles into a quasi-hydrostatic pre-main-sequence sphere of radius R ≈ 3 to 15 R_☉ and low interior temperature (T_c ~ 10⁶ K). At these cool interior temperatures, Kramers bound-free and free-free radiative opacity (κ_rad ∝ ρ T^−3.5) is exceedingly high throughout the interior, making radiative heat transport inefficient and driving vigorous turbulent convection from the center to just beneath the photosphere.

Because convective mixing in the high-density interior is highly efficient, the specific entropy s is uniform throughout the star (ds/dr ≈ 0), and pressure P and density ρ obey the adiabatic ideal-gas relation P = K ρ^γ with adiabatic exponent γ = 5/3 = 1 + 1/n, corresponding to an Emden polytrope of index n = 1.5. From homologous polytropic scaling, the central pressure P_c, central density ρ_c, and adiabatic constant K of an n = 1.5 sphere of mass M and radius R satisfy:

P = K ρ^5/3, where K = 0.4242 G M^1/3 R

Photospheric Boundary Condition and H⁻ Bound-Free Opacity

While the interior adiabatic constant K is fixed purely by M and R, the star must match smoothly onto a radiative photosphere at optical depth τ_phot = 2/3, where photons escape freely into space at effective temperature T_eff. Hydrostatic equilibrium across the atmospheric scale height H_p = k_BT_eff / (μ m_u g) down to τ = ∫ κ_phot ρ dr = 2/3 requires that the photospheric gas pressure P_phot satisfy:

P_phot ≈ ( 2 / 3 ) ( g / κ_phot(ρ_phot, T_eff) ) = ( 2 G M ) / ( 3 R² κ_phot )

In cool stellar atmospheres (T_eff ≈ 2,800 to 5,000 K), thermal energy (k_BT ≈ 0.3–0.4 eV) is far too low to ionize neutral hydrogen (ionization potential χ_H = 13.6 eV), so free electrons are supplied almost exclusively by easily ionized alkali and alkaline-earth metals (sodium, potassium, calcium, magnesium, iron, with ionization potentials χ_m ≈ 5–8 eV). These free electrons attach weakly to abundant neutral hydrogen atoms to form the negative hydrogen ion (H⁻), which has a tiny electron affinity (binding energy) of χ(H⁻) = 0.754 eV (corresponding to a photo-detachment threshold wavelength of λ = 1.64 μm in the near-infrared).

Consequently, photospheric opacity in cool stars is governed by H⁻ bound-free and free-free absorption, which follows a power law that increases steeply with both density and temperature:

κ_H⁻ = κ₀ Z ρ^a T_eff^b (with a ≈ 0.5, b ≈ 4 to 9, and metallicity Z)

Equating the interior adiabatic relation P_phot = K ρ_phot^5/3 (with ρ_phot = P_phot μ m_u / k_BT_eff) to the photospheric hydrostatic condition P_phot = 2*GM* / (3R² κ_H⁻) and eliminating radius R via the Stefan-Boltzmann law (L = 4π R² σ_SB T_eff⁴) yields the analytical Hayashi Track Scaling Law:

log₁₀ T_eff ≈ ( +0.04 to +0.08 ) log₁₀( L / L_☉ ) + ( +0.15 to +0.20 ) log₁₀( M / M_☉ ) − 0.08 log₁₀ Z + constant

Because the coefficient multiplying log₁₀(L/L_☉) is close to zero (0.04 to 0.08), a 100-fold decrease in luminosity (Δ log₁₀ L = −2) changes log₁₀ T_eff by less than 0.10 dex, producing a nearly vertical evolutionary track on the H–R diagram.

Physical Meaning of the Hayashi Forbidden Zone

What happens if a hypothetical star of mass M and luminosity L is placed at an effective temperature T_eff cooler than the Hayashi limit T_Hayashi(M, L, Z) (to the right of the Hayashi track)?

• Lowering T_eff decreases the H⁻ photospheric opacity κ_H⁻ ∝ T_eff^b precipitously, which raises the photospheric pressure P_phot ∝ 1/κ_H⁻ and would require an interior polytropic index n < 1.5 (a superadiabatic temperature gradient throughout the deep interior).
• Because deep stellar convection has an enormous heat capacity, any superadiabatic excess in the deep interior immediately triggers supersonic convective plumes that transport heat to the surface on a dynamical free-fall timescale (t_dyn ~ ( G ρ̄ )^−1/2 ~ hours to days), raising the surface temperature T_eff back onto the Hayashi boundary n = 1.5. Thus, the region to the right of the Hayashi track is a forbidden zone for any star in hydrostatic equilibrium.

As a pre-main-sequence star contracts along the Hayashi track on its Kelvin-Helmholtz thermal timescale:

t_KH = ( 3 G M² ) / ( 7 R L )

(where the factor 3/7 applies to a fully convective n = 1.5 polytrope), half of the released gravitational potential energy is radiated away as luminosity L while the other half heats the interior according to the virial theorem (T_c ∝ M / R). During this descent, the star burns primordial deuterium (²H) via ²H(p,γ)³He at T_c ≈ 1.0 × 10⁶ K (briefly pausing contraction along the "deuterium birthline" for ~10⁵–10⁶ yr) and subsequently destroys lithium-7 (⁷Li) via ⁷Li(p,α)⁴He at T_c ≈ 2.5 × 10⁶ K.

How the star completes its approach to the Zero-Age Main Sequence (ZAMS) depends strongly on its stellar mass M:

The Transition to the Radiative Henyey Track (0.5 to 3.0 M_☉)

For a 1.0 M_☉ young solar analog, downward contraction along the Hayashi track proceeds for ~10 million years until its luminosity drops to L ≈ 0.5–0.6 L_☉ and its central temperature reaches T_c ≈ 4 × 10⁶ K. Because Kramers radiative opacity falls rapidly with rising temperature (κ_rad ∝ ρ T^−3.5), the central core becomes transparent enough for radiation to carry the outward heat flux without convection. A radiative core nucleates at the center and grows outward in mass fraction.

Once the radiative core encompasses the majority of the star's mass, the star's luminosity is governed by radiative diffusion across the interior rather than the surface H⁻ valve (L ∝ M^5.5 R^−0.5 μ^7.5). Consequently, as radius R continues to contract, luminosity L stops falling and instead rises slightly by 0.2–0.4 dex while T_eff increases sharply toward the left, tracing the nearly horizontal Henyey track for another ~20–30 million years until central temperatures reach T_c ≈ 1.3 × 10⁷ K and hydrostatic proton–proton fusion stabilizes the star on the Zero-Age Main Sequence.

Post-Main-Sequence Ascent Along the Hayashi Limit

The identical physical boundary operates in reverse at the end of a star's main-sequence lifetime. When a post-main-sequence subgiant with an inert helium core expands its outer envelope eastward across the H–R diagram (crossing the Hertzsprung gap on a thermal timescale), its surface temperature cools until T_eff hits the Hayashi limit (~4,500–5,000 K). Because the star cannot enter the Hayashi forbidden zone to the right, its outer layers become deeply convective, and further expansion of the radius from 3 R_☉ to >150 R_☉ forces the star to climb nearly vertically upward along the Hayashi track as a Red Giant Branch star up to the Tip of the Red-Giant Branch and later as an asymptotic giant branch star.

1. Hayashi, C. (1961). "Stellar evolution in early phases of gravitational contraction." Publications of the Astronomical Society of Japan, 13, 450–452. https://ui.adsabs.harvard.edu/abs/1961PASJ...13..450H
2. Hayashi, C., Hōshi, R., & Sugimoto, D. (1962). "Evolution of the Stars." Progress of Theoretical Physics Supplement, 22, 1–183. https://doi.org/10.1143/PTPS.22.1
3. 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
4. Stahler, S. W. (1988). "Understanding young stars: A history." Publications of the Astronomical Society of the Pacific, 100(634), 1474–1485. https://doi.org/10.1086/132352
5. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. (2015). "New evolutionary models for pre-main sequence and main sequence low-mass stars down to the hydrogen-burning limit." Astronomy & Astrophysics, 577, A42. https://doi.org/10.1051/0004-6361/201425481