Pre-Main-Sequence Evolution

A pre-main-sequence star is a star caught in the act of not yet being one. The infalling envelope that built it has cleared, so it is optically visible and in near-hydrostatic equilibrium — but its core is still too cool to burn hydrogen at the rate that will one day support it. In that gap it has one power source: the gravitational energy released as it slowly shrinks. This is the bridge between star formation — the Jeans instability, molecular clouds, protostellar collapse — and the theory of stellar structure and evolution. The evolution is quasi-static, and that word carries the physics. At every instant the star obeys hydrostatic equilibrium, yet its radius creeps inward on the Kelvin–Helmholtz (thermal) timescale — long compared with the free-fall time, short compared with the nuclear one. The virial theorem then does something counterintuitive: of the gravitational energy liberated, only half escapes as light while the other half heats the interior, so the star gets hotter as it radiates. Because opacity in the cool outer layers is enormous, low-mass stars begin fully convective, which pins their surface temperature and produces the near-vertical Hayashi track in the Hertzsprung–Russell diagram — a locus bounding a forbidden region no hydrostatic star can enter. As contraction heats the interior the opacity falls, a radiative core grows, and the star peels off onto the nearly horizontal Henyey track toward higher effective temperature. The journey ends at the zero-age main sequence, when core hydrogen burning first supplies the entire luminosity and contraction stops — unless the star is too light to get there, in which case it stalls forever as a brown dwarf.

Key takeaways

A pre-main-sequence (PMS) star has left its protostellar accretion phase but has not yet ignited stable core hydrogen burning; it shines on gravitational energy released by slow Kelvin–Helmholtz contraction.
The virial theorem splits the released energy exactly in half — half radiated, half retained as heat — so a contracting self-gravitating gas has negative heat capacity and warms toward ignition as it loses energy.
Low-mass PMS stars are fully convective and descend the near-vertical Hayashi track at almost fixed $T_{\mathrm{eff}}\approx 3500\text{–}4500\,\mathrm{K}$ ; a forbidden region lies to its cool side.
Once a radiative core forms the star turns left onto the nearly horizontal Henyey track, contracting at roughly constant luminosity while $T_{\mathrm{eff}}$ climbs to the main sequence.
The contraction clock is $t_{\mathrm{KH}}=GM^2/(RL)\approx 3\times10^{7}\,\mathrm{yr}$ for the Sun, $\sim 10^{5}\,\mathrm{yr}$ for massive stars and $\sim 10^{9}\,\mathrm{yr}$ for the lowest masses; below $\approx 0.075\,M_\odot$ the core never reaches $10^{7}\,\mathrm{K}$ and the object cools forever as a brown dwarf.
Key idea: a pre-main-sequence star is a gravitationally contracting, quasi-hydrostatic gas sphere that converts gravitational potential energy into radiation on the Kelvin–Helmholtz timescale — descending the convective Hayashi track and crossing the radiative Henyey track — until its central temperature reaches $\sim 10^{7}\,\mathrm{K}$ and hydrogen fusion switches on at the main sequence.

Where this lives today

Pre-main-sequence theory is the working backbone of star-formation astrophysics, and it is being rewritten right now by data. Fitting theoretical PMS isochrones to the color–magnitude diagrams of nearby star-forming regions — the Orion Nebula Cluster, Taurus, Upper Scorpius — yields ages of $1$ – $10\,\mathrm{Myr}$ and the star-formation history of molecular clouds; Gaia has sharpened this with precise parallaxes and clean membership, exposing real age spreads and offsets between track sets. The lithium depletion boundary — the mass at which fully convective stars have just failed to burn their surface lithium — gives some of the most robust absolute ages available and calibrates the method. JWST now probes the youngest embedded objects along the birthline and pins the substellar boundary near $0.075\,M_\odot$ . And eclipsing PMS binaries, where masses and radii are measured directly, remain the benchmark — repeatedly showing that magnetic activity and starspots inflate young stars above the radii standard models predict, an open problem current tracks still work to resolve.

The physical picture

Picture a hot balloon deflating under its own weight in slow motion. As it shrinks, gravity does work on the gas; the virial theorem sends half of that work out as radiation and banks the other half as internal heat. This is the strangest fact about self-gravitating systems: they behave as though their heat capacity were negative — take energy away and they warm up. A PMS star radiates for tens of millions of years while its central temperature climbs toward hydrogen ignition, powering its own future. The Hayashi track is best read as a boundary, not a road. To its cool side, convection cannot carry the luminosity the star must radiate, and no static solution exists. A contracting low-mass star, fully convective throughout, is pinned against this boundary at nearly fixed $T_{\mathrm{eff}}$ because its surface is controlled by the fierce temperature sensitivity of H $^-$ opacity. As it shrinks at constant $T_{\mathrm{eff}}$ , its luminosity $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ can only fall — so it slides straight down. The first figure traces that descent for a solar-mass star, and the turn onto the Henyey track when a radiative core forms.

Figure

The pre-main-sequence track of the Sun in the Hertzsprung–Russell diagram: it descends the near-vertical Hayashi track at Teff ≈ 4400 K from about 10 down to 1 solar luminosity, then turns left onto the Henyey track and reaches the zero-age main sequence at Teff ≈ 5800 K after roughly 30 million years — the shaded region to the right is forbidden to any hydrostatic star

Why do low-mass stars stay convective while massive ones go radiative? Opacity again. In cool outer layers the opacity is so high that radiation could only carry the luminosity down an impossibly steep gradient; convection wins instead. As contraction warms the interior, the opacity drops and a radiative core grows outward from the center — the moment the star leaves the Hayashi boundary and turns onto the Henyey track. The PMS phase is then a race between the contraction clock $t_{\mathrm{KH}}$ and the ignition threshold, and mass decides the winner. The second figure shows how brutally that clock depends on mass: massive stars reach the main sequence in $10^5$ years, while the least massive are still contracting after $10^9$ .

Figure

The punchline: pre-main-sequence contraction time t_KH = GM²/RL versus stellar mass — it falls steeply from about 10⁹ yr at 0.1 solar masses through 3×10⁷ yr for the Sun to only 10⁵ yr at 10 solar masses, so massive stars arrive on the main sequence while the least massive are still shrinking

The core equations

Everything begins with mechanical balance. At each radius the outward pressure gradient must support the weight of the overlying gas, or the star would collapse on a free-fall time — the equation it satisfies at every instant of its slow contraction.

\frac{\mathrm{d}P}{\mathrm{d}r}=-\frac{G\,m(r)\,\rho(r)}{r^{2}} \qquad (1)

where: $P$ — pressure at radius $r$ ; $m(r)$ — mass enclosed within $r$ ; $\rho(r)$ — local mass density; $G$ — gravitational constant. The pressure must drop outward fast enough that its gradient carries the local weight $g\rho=Gm\rho/r^2$ . Integrate that balance against volume over the star and the pointwise equation collapses to a single global relation between heat and gravity — the scalar virial theorem, the master key to the PMS phase.

2K+\Omega=0,\qquad E_{\mathrm{tot}}=K+\Omega=\tfrac{1}{2}\Omega \qquad (2)

where: $K$ — total internal (kinetic) energy of the gas; $\Omega$ — gravitational potential energy (negative); $E_{\mathrm{tot}}$ — total energy. The punchline is the middle result: the total energy is exactly half the (negative) gravitational energy, so a bound star sits at $E_{\mathrm{tot}}<0$ with a thermal reserve of $|\Omega|/2$ . Differentiate the virial relation in time and you learn what powers the star: because $E_{\mathrm{tot}}=\tfrac12\Omega$ , the luminosity radiated is half the rate at which gravitational energy is released.

L=-\frac{\mathrm{d}E_{\mathrm{tot}}}{\mathrm{d}t}=-\frac{1}{2}\frac{\mathrm{d}\Omega}{\mathrm{d}t} \qquad (3)

where: $L$ — surface luminosity powered by contraction; $\Omega$ — gravitational energy, growing more negative as $R$ shrinks. As $R$ falls, $-\dot\Omega>0$ and $L>0$ ; simultaneously $K=-E_{\mathrm{tot}}$ rises, so the interior heats even as it radiates. Dividing the available reservoir by the luminosity gives the natural clock of the phase — how long gravity alone can keep the star shining.

t_{\mathrm{KH}}=\frac{G M^{2}}{R\,L} \qquad (4)

where: $M$ — stellar mass; $R$ — radius; $L$ — luminosity. This Kelvin–Helmholtz (thermal) timescale is $\approx 3\times10^{7}\,\mathrm{yr}$ for the Sun — far shorter than its age, which is how we know the main-sequence Sun is not contraction-powered, yet it is the correct duration of the PMS phase. The surface ties luminosity, radius, and temperature together as a blackbody. With $T_{\mathrm{eff}}$ pinned by convective-envelope physics, this relation is what forces the Hayashi descent to be vertical.

L\approx 4\pi R^{2}\,\sigma\,T_{\mathrm{eff}}^{4} \qquad (5)

where: $\sigma$ — Stefan–Boltzmann constant; $T_{\mathrm{eff}}$ — effective temperature; $R$ — photospheric radius. Hold $T_{\mathrm{eff}}$ fixed and $L\propto R^{2}$ : a shrinking star can only grow fainter — straight down the HR diagram. Whether a layer carries heat by radiation or convection is decided by comparing two temperature gradients. In the cool, opaque PMS envelope the radiative gradient wins big, which is why these stars start out convective.

\nabla_{\mathrm{rad}}=\frac{3}{16\pi a c G}\frac{\kappa L P}{m\,T^{4}}>\nabla_{\mathrm{ad}}\ \Rightarrow\ \text{convection} \qquad (6)

where: $\nabla_{\mathrm{rad}}$ — radiative gradient $\mathrm{d}\ln T/\mathrm{d}\ln P$ needed to carry $L$ by radiation; $\nabla_{\mathrm{ad}}$ — adiabatic gradient ( $=0.4$ for an ideal monatomic gas); $\kappa$ — opacity (huge in cool envelopes); $a$ — radiation constant. When $\nabla_{\mathrm{rad}}$ exceeds $\nabla_{\mathrm{ad}}$ the layer is unstable and convection takes over. Finally, the virial theorem applied locally gives the central temperature versus how far the star has contracted — the thermometer that governs when hydrogen ignites.

T_c\sim \frac{G M\,\mu\,m_{\mathrm{H}}}{k_{\mathrm{B}}\,R} \qquad (7)

where: $T_c$ — central temperature; $\mu$ — mean molecular weight; $m_{\mathrm{H}}$ — hydrogen mass; $k_{\mathrm{B}}$ — Boltzmann constant. As $R$ shrinks $T_c$ rises; crossing $\sim 10^{7}\,\mathrm{K}$ ignites the proton–proton chain and halts contraction.

Derivation AThe virial theorem and the negative heat capacity of a contracting star

The goal is to derive, from hydrostatic equilibrium alone, the relation that powers pre-main-sequence evolution: a self-gravitating ideal-gas sphere has $E_{\mathrm{tot}}=\tfrac12\Omega$ , so when it radiates energy away its internal temperature actually rises. This "negative heat capacity" is the engine of contraction, and it follows inexorably from mechanical balance plus the ideal-gas law. We integrate the balance equation against volume, integrate by parts, insert the equation of state, and apply energy conservation.

Step 1. Write hydrostatic balance in the mass coordinate, $\mathrm{d}P/\mathrm{d}m=-Gm/(4\pi r^{4})$ , using $\mathrm{d}m=4\pi r^{2}\rho\,\mathrm{d}r$ . Why: the mass coordinate makes the star a fixed integration domain from $0$ to $M$ even as $r(m)$ contracts — exactly what a global energy statement needs.

Step 2. Multiply by the volume $V=\tfrac{4}{3}\pi r^{3}$ and integrate over the star. On the right this produces $\int_0^M -Gm/r\,\mathrm{d}m=\Omega$ , the gravitational potential energy. Why: pairing the pressure gradient with the enclosed volume converts a force balance into an energy balance.

Step 3. Integrate the left side by parts. The boundary term vanishes because $P\to0$ at the surface, leaving $-3\int_0^M (P/\rho)\,\mathrm{d}m=\Omega$ . Why: this trades the awkward $\mathrm{d}P/\mathrm{d}m$ for a transparent pressure-per-unit-mass integral.

Step 4. Insert the monatomic ideal-gas relation $P/\rho=(\gamma-1)u=\tfrac{2}{3}u$ , where $\int u\,\mathrm{d}m=K$ . This turns the left side into $-2K$ . Why: for $\gamma=5/3$ the pressure integral is exactly $\tfrac23$ of the internal energy — the factor that makes the theorem $2K$ .

Step 5. Combine to get the scalar virial theorem, $2K+\Omega=0$ . Then $E_{\mathrm{tot}}=K+\Omega=\tfrac12\Omega<0$ and $K=-E_{\mathrm{tot}}$ . Why: this is the payoff — total energy is half the gravitational energy, and thermal energy tracks minus the total energy.

Step 6. Differentiate: $L=-\dot E_{\mathrm{tot}}=-\tfrac12\dot\Omega>0$ while $\dot K>0$ . As the star radiates, $E_{\mathrm{tot}}$ falls and $K\propto T$ rises. Why: negative heat capacity made explicit — losing energy heats the gas, driving $T_c$ toward ignition.

Go deeper · Derivation B · 1 minThe Kelvin–Helmholtz contraction time for the Sun

Derivation BThe Kelvin–Helmholtz contraction time for the Sun

The aim is to estimate, from the gravitational reservoir alone, how long contraction can power a star of solar luminosity — the number that both sank the nineteenth-century contraction theory of the Sun and correctly gives the pre-main-sequence duration. We take the binding energy, divide by the luminosity, and reduce to years.

Step 1. The gravitational energy of a bound sphere is $|\Omega|\sim \tfrac{3}{10}\,GM_\odot^{2}/R_\odot$ for uniform density (order unity for a realistic centrally condensed star). Why: every reasonable density profile gives $|\Omega|$ within a factor of a few of $GM^2/R$ — all an order-of-magnitude clock needs.

Step 2. Evaluate $GM_\odot^{2}/R_\odot=(6.67\times10^{-11})(1.99\times10^{30})^{2}/(6.96\times10^{8})\approx 3.8\times10^{41}\,\mathrm{J}$ . Why: this is the total reservoir; by the virial theorem only half is ever free to radiate, the rest going into heat.

Step 3. Divide by $L_\odot=3.83\times10^{26}\,\mathrm{W}$ : $t_{\mathrm{KH}}\sim (3.8\times10^{41})/(3.83\times10^{26})\approx 1.0\times10^{15}\,\mathrm{s}$ . Why: dropping the order-unity structure factor keeps this honest as an order of magnitude.

Step 4. Convert to years: $1.0\times10^{15}\,\mathrm{s}\div 3.16\times10^{7}\,\mathrm{s\,yr^{-1}}\approx 3\times10^{7}\,\mathrm{yr}$ . Why: this "Kelvin–Helmholtz age" is a thousand times too short to be the Sun's true age — proving contraction cannot power the main sequence — yet is exactly the timescale on which the pre-main-sequence Sun contracts. *The direct evaluation gives $GM_\odot^{2}/(R_\odot L_\odot)\approx 9.9\times10^{14}\,\mathrm{s}\approx 3.1\times10^{7}\,\mathrm{yr}$ .*

Go deeper · Derivation C · 1 minWhy a forbidden region lies to the right of the Hayashi track

Derivation CWhy a forbidden region lies to the right of the Hayashi track

The purpose here is to show why the Hayashi track behaves as a wall rather than a path: for a fully convective star of given mass and radius there is a minimum effective temperature below which no hydrostatic solution exists, so contracting low-mass stars are pressed against that limit and slide down it. The argument links the fixed entropy of a convective interior to the photospheric boundary condition and the fierce temperature sensitivity of H $^-$ opacity.

Step 1. A fully convective star mixes to nearly uniform specific entropy, so its interior follows a single adiabat $P\propto T^{5/2}$ (an $n=3/2$ polytrope), fixing $P(T)$ once the entropy is chosen. Why: reducing the interior to a one-parameter family makes the surface matching a sharp constraint.

Step 2. The photosphere sits where optical depth $\tau\approx 2/3$ ; matching the deep adiabat outward to that boundary selects a specific $T_{\mathrm{eff}}$ for each $M$ and $R$ . Why: the interior and atmosphere must join continuously, tying the entropy to an observable surface temperature.

Step 3. In these cool envelopes the dominant opacity is the H $^-$ ion, $\kappa\propto \rho^{0.5}\,T^{9}$ — extraordinarily steep in temperature, so a small drop in $T_{\mathrm{eff}}$ demands a large structural change the adiabatic interior cannot supply. Why: the enormous exponent on $T$ turns a mild preference into a rigid wall in the HR diagram.

Step 4. Below the Hayashi $T_{\mathrm{eff}}$ the required gradient exceeds the adiabatic gradient everywhere and no static structure exists; such a star is dynamically unstable. Hence the cool side is forbidden and a contracting star rides the boundary at nearly constant $T_{\mathrm{eff}}$ . Why: this is why the Hayashi track is a vertical descent — $T_{\mathrm{eff}}$ is pinned while only $R$ , hence $L$ , changes.

Limiting cases

Stellar mass reorganizes every prediction — it sets the contraction time, the balance of convective versus radiative interior, and whether hydrogen ever ignites.

Case	Result	Meaning
Very low mass ( $M\lesssim 0.5\,M_\odot$ )	$t_{\mathrm{KH}}\sim 10^{8}\text{–}10^{9}\,\mathrm{yr}$	Fully convective almost to the ZAMS; a long Hayashi descent with almost no Henyey portion.
Solar mass ( $M\approx 1\,M_\odot$ )	$t_{\mathrm{KH}}\approx 3\times10^{7}\,\mathrm{yr}$	Descends the Hayashi track, grows a radiative core, then runs a short Henyey track before ignition.
Intermediate / high mass ( $M\gtrsim 2\,M_\odot$ )	$t_{\mathrm{KH}}\propto M^{2}/(RL)\sim 10^{5}\text{–}10^{6}\,\mathrm{yr}$	Radiative interior forms early; the star spends most of its brief PMS life on the Henyey track.
Substellar ( $M\lesssim 0.075\,M_\odot$ )	$T_{c,\max}<10^{7}\,\mathrm{K}$	Degeneracy halts contraction before ignition; the object becomes a brown dwarf and cools forever (deuterium limit $\sim 0.013\,M_\odot$ ).

Common trapContraction radiates all of the released gravitational energy

A tempting but wrong picture is that a contracting star converts the full drop in gravitational energy $|\Delta\Omega|$ into radiation. The virial theorem for a non-relativistic self-gravitating gas says otherwise: exactly half of the liberated energy leaves as light and half stays behind as internal heat. That retained half is not a bookkeeping curiosity — it is why the interior gets hotter as the star radiates, why $T_c$ climbs toward $10^{7}\,\mathrm{K}$ , and why contraction is a self-terminating process that ends in ignition rather than a slow fade. Treat the star as radiating all of $|\Delta\Omega|$ and you overpredict the luminosity by a factor of two and lose the reason the star evolves at all.

Worked examples

1 · The Kelvin–Helmholtz age of the Sun (fully worked)

Using only the solar mass, radius, and luminosity, estimate how long gravitational contraction could power the Sun.

Step 1. Assemble $GM_\odot^{2}/R_\odot$ with $M_\odot=1.99\times10^{30}\,\mathrm{kg}$ , $R_\odot=6.96\times10^{8}\,\mathrm{m}$ , $G=6.67\times10^{-11}$ SI: the numerator is $(6.67\times10^{-11})(1.99\times10^{30})^{2}=2.64\times10^{50}$ , and dividing by $R_\odot$ gives $GM_\odot^{2}/R_\odot=\mathbf{3.8\times10^{41}\ J}$ .

Step 2. Divide by the solar luminosity $L_\odot=3.83\times10^{26}\,\mathrm{W}$ to get the timescale: $t_{\mathrm{KH}}=GM_\odot^{2}/(R_\odot L_\odot)=(3.8\times10^{41})/(3.83\times10^{26})=\mathbf{9.9\times10^{14}\ s}$ .

Step 3. Convert to years with $1\,\mathrm{yr}=3.16\times10^{7}\,\mathrm{s}$ : $t_{\mathrm{KH}}=9.9\times10^{14}/3.16\times10^{7}\approx\mathbf{3\times10^{7}\ yr}$ . This is the historic contradiction — thirty million years is a thousand times shorter than the geological age of the Earth, so contraction cannot be the Sun's steady power source. But it is exactly right as the duration of the Sun's pre-main-sequence contraction. The solar Kelvin–Helmholtz time is t_KH ≈ 3×10⁷ yr.

2 · Central temperature at ignition (fully worked)

Estimate the central temperature of the present-day Sun from virial scaling, and check that contraction naturally delivers a star to ignition.

Step 1. Use $T_c\sim GM\mu m_{\mathrm{H}}/(k_{\mathrm{B}}R)$ with fully ionized composition $\mu\approx 0.6$ , $m_{\mathrm{H}}=1.67\times10^{-27}\,\mathrm{kg}$ , $k_{\mathrm{B}}=1.38\times10^{-23}\,\mathrm{J\,K^{-1}}$ . The numerator $GM_\odot\mu m_{\mathrm{H}}=(6.67\times10^{-11})(1.99\times10^{30})(0.6)(1.67\times10^{-27})=1.33\times10^{-7}$ .

Step 2. Divide by $k_{\mathrm{B}}R_\odot=(1.38\times10^{-23})(6.96\times10^{8})=9.6\times10^{-15}$ : $T_c\sim (1.33\times10^{-7})/(9.6\times10^{-15})\approx\mathbf{1.4\times10^{7}\ K}$ , within ten per cent of the true solar value $1.5\times10^{7}\,\mathrm{K}$ despite the crude estimate.

Step 3. Read equation (7) as a process: $T_c\propto M/R$ , so as a fixed-mass star contracts, $T_c$ rises with $1/R$ . A protostar begins at $R\sim 4\,R_\odot$ with $T_c$ a few $\times10^{6}\,\mathrm{K}$ ; by $R\approx R_\odot$ the center has crossed $10^{7}\,\mathrm{K}$ , the proton–proton chain switches on, and contraction halts. On the Hayashi descent itself the same shrinkage from $4\,R_\odot$ to $R_\odot$ at fixed $T_{\mathrm{eff}}$ dims the star by $(1/4)^2=16$ , since $L\propto R^2$ — the near-vertical plunge of the first figure. Contraction to R ≈ R⊙ brings T_c to ≈ 1.4×10⁷ K, igniting hydrogen.

Physics of the portrait

The theory of pre-main-sequence evolution is really the story of two tracks, and the deeper of the two — the vertical wall that low-mass stars descend — bears one name. For a decade after the Second World War, theorists assumed young contracting stars were radiative and computed the roughly horizontal Henyey tracks. It took a careful accounting of convection and cool-envelope opacity to overturn that picture and reveal a forbidden zone in the HR diagram.

Origins · 1 minChushiro Hayashi and how this was discovered

Chushiro Hayashi

1920–2010 · stellar structure and the fully convective pre-main-sequence track

Chushiro Hayashi was a Japanese astrophysicist whose 1961 analysis transformed how we understand young stars. Working at Kyoto University, he asked what surface temperature a fully convective star of a given mass and radius could actually have, and found that below a critical $T_{\mathrm{eff}}$ no hydrostatic solution exists — carving out the "forbidden region" on the cool side of what is now the Hayashi track and showing that low-mass pre-main-sequence stars must descend it nearly vertically. This overturned the prevailing radiative picture, in which contracting stars were thought to evolve almost horizontally, and it explained why the youngest stars in a cluster sit high and to the right in the HR diagram. Hayashi's reach was broad — a decade earlier he had corrected Gamow's treatment of the early universe by including electron–positron pairs, foundational to Big Bang nucleosynthesis. But the line he drew across the HR diagram remains one of the cleanest results in stellar structure: a boundary set entirely by convection and H $^-$ opacity.

History

1861

— Kelvin and Helmholtz propose gravitational contraction as the Sun's power source

1926

— Eddington's "Internal Constitution of the Stars" formalizes the virial and thermal balance of stars

1955

— Henyey, Lelevier and Levée compute the radiative pre-main-sequence contraction track

1961

— Hayashi identifies the fully convective track and the forbidden region

1997

— Baraffe and collaborators release low-mass PMS models widely used for dating young clusters

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., chs. 10 and 12 (Hayashi and Henyey tracks, PMS timescales) · Kippenhahn, Weigert and Weiss, Stellar Structure and Evolution, 2nd ed., chs. 3, 22, 24 (virial theorem, the Hayashi line) · Prialnik, An Introduction to the Theory of Stellar Structure and Evolution, ch. 8 · Stahler and Palla, The Formation of Stars (2004), chs. 16–18. Solar constants used throughout: $M_\odot=1.99\times10^{30}\,\mathrm{kg}$ , $R_\odot=6.96\times10^{8}\,\mathrm{m}$ , $L_\odot=3.83\times10^{26}\,\mathrm{W}$ (IAU standard). The burning limits $\approx 0.075\,M_\odot$ (hydrogen) and $\approx 0.013\,M_\odot$ (deuterium) are standard. Key anchor results: $t_{\mathrm{KH}}=GM_\odot^{2}/(R_\odot L_\odot)\approx 9.9\times10^{14}\,\mathrm{s}\approx 3\times10^{7}\,\mathrm{yr}$ ; the virial theorem gives $2K+\Omega=0$ and $E_{\mathrm{tot}}=\tfrac12\Omega$ ; the central-temperature estimate yields $T_c\sim GM_\odot\mu m_{\mathrm{H}}/(k_{\mathrm{B}}R_\odot)\approx 1.4\times10^{7}\,\mathrm{K}$ .

Prerequisites. Molecular Clouds · The Jeans Instability and Jeans Mass · Virial Theorem for Molecular Clouds. Leads to. Main-Sequence Stars and Hydrogen Burning · Brown Dwarfs and the Substellar Boundary · Stellar Structure Equations.

Ask the tutor: Ask me to walk through why a contracting star heats up, or to work an isochrone age for a young cluster step by step.

What is the difference between a protostar and a pre-main-sequence star?

They are consecutive phases of the same forming star, separated by whether it is still gaining mass. A protostar is deeply embedded and grows by accreting from an infalling envelope and disk; it is visible only in the infrared and submillimeter because dust reprocesses its light. A pre-main-sequence star is what remains once accretion has essentially stopped and the envelope has cleared: optically visible, in near-hydrostatic equilibrium, and powered by contraction rather than accretion. The transition is marked by the birthline, where a star of a given mass first appears as an optically visible object. Classical T Tauri and Herbig Ae/Be stars are the low- and intermediate-mass faces of this phase.

Why can only half of the gravitational energy be radiated?

Because the star must stay in hydrostatic equilibrium, which the virial theorem enforces. When it contracts and releases gravitational energy, $2K+\Omega=0$ requires the internal energy $K$ to rise by half that amount to keep the pressure supporting the more compact structure; only the remaining half is free to radiate. If the star tried to radiate all of it, it would be left under-supported and collapse further until balance returned. So the fifty-fifty split is not an accident — it is the condition for the star to remain a star while it shrinks, and exactly why contraction warms rather than cools the interior.

Do all stars pass through both a Hayashi and a Henyey track?

No — the mix depends on mass. Stars below roughly half a solar mass stay fully convective almost to the main sequence, so they descend a long Hayashi track with essentially no Henyey portion. A solar-mass star runs a substantial Hayashi descent, then a short Henyey segment once its radiative core has grown. Intermediate- and high-mass stars develop radiative interiors very early, so they barely touch the Hayashi track and traverse a long Henyey track. The general rule: higher mass means an earlier, larger radiative core, hence a longer Henyey and shorter Hayashi phase.

Why do brown dwarfs never reach the main sequence?

Because their cores become degenerate before they get hot enough. As any contracting body shrinks, its central temperature rises like $1/R$ — but in a low-mass object, electron degeneracy pressure grows fast enough to halt contraction while the center is still below the $\sim 10^{7}\,\mathrm{K}$ needed for sustained fusion. Once degeneracy supports the body, contraction stops, the temperature peaks and declines, and it simply cools at nearly constant radius. The dividing line sits near $0.075\,M_\odot$ : above it, ignition beats degeneracy; below it, the object becomes a brown dwarf. Such bodies may still fuse deuterium transiently down to about $0.013\,M_\odot$ , the conventional planet–brown dwarf boundary.

How do astronomers actually measure the age of a young cluster this way?

They use the fact that PMS luminosity and temperature evolve predictably with time. All stars in a cluster form at roughly the same epoch, so in a color–magnitude diagram the low-mass members trace out an isochrone — a snapshot of many masses at one common age. Overlaying theoretical PMS tracks and reading off the best-matching age dates the cluster, typically $1$ – $10\,\mathrm{Myr}$ for active regions. The method needs careful correction for extinction and reddening by dust, and its accuracy is limited by disagreements among track sets — which is why the lithium depletion boundary and Gaia parallaxes are used to calibrate and cross-check it.

If contraction cannot power the Sun today, why does the Kelvin–Helmholtz time still matter?

Because it is the correct timescale for a different phase of the Sun's life. When Kelvin and Helmholtz proposed contraction as the Sun's power source, the thirty-million-year lifetime it implied was fatally short against the geological evidence — a contradiction resolved only by nuclear fusion. But the same number is exactly right as the duration of the Sun's pre-main-sequence contraction: for tens of millions of years before hydrogen ignited, contraction really was its only power source. So $t_{\mathrm{KH}}$ was never wrong physics — it was the right answer to the wrong question, and it now tells us how long young stars linger on their contraction tracks.