The specific intensity of a blackbody is direction-independent over its outward hemisphere. Its outward flux per unit frequency is therefore
Multiplying by the stellar surface area gives the spectral luminosity
Using the Planck law, the hydrogen-ionizing photon production rate follows by dividing each spectral contribution by its photon energy and integrating above the threshold:
With and , this becomes
For the initial hydrogen-only photon balance, assume an ionization-bounded nebula at steady photoionization equilibrium, negligible dust absorption and ionizing-photon leakage, negligible collisional ionization, and local reabsorption of photons from direct ground-state recombinations. A ground-state recombination then causes another photoionization and produces no net removal of a proton. Net removal is counted by Case B recombination, whose coefficient is . Integrating the local balance gives
The stellar temperature sets the photon supply; the gas temperature sets the recombination coefficients. They are not generally equal. For uniform hydrogen and constant gas temperature, the Strömgren radius is
with in fully ionized pure hydrogen.
For the hydrogen-helium mixture, define the blackbody ionizing photon function
The three relevant thresholds are approximately , and , for , and respectively. Higher stellar temperature increases the relative supply of the harder photons. These cumulative photon rates cannot be assigned independently to the three species: a photon above a helium threshold can instead ionize hydrogen, and helium recombination radiation can also ionize hydrogen. The absorption cross-sections and the diffuse radiation determine the effective allocation.
A useful sharp-front model of nested hydrogen-helium ionization zones has an inner sphere of radius , a shell ending at , and an outer hydrogen-ionized region ending at , with . Let be uniform nuclear number densities. The electron number densities in these three regions are , and . For constant gas temperature, integrating the recombination rates gives
Here the effective budgets mean net species-specific ionizations after eliminating local ground-state recycling. Their calculation must include photon competition and any interspecies diffuse contribution; the helium coefficients label the recombining ion's emitted spectrum. Once these budgets are specified, solve the third equation for , the second for , and the first for .
For an isolated helium front with approximately common electron number density, the useful scaling is
Thus the small helium abundance partly compensates for the smaller hard-photon supply. A cooler central star produces little and can have a substantially smaller helium-ionized zone. A sufficiently hot star can ionize helium throughout most of the hydrogen region. If independent radius estimates place a helium front beyond the hydrogen front, the independent photon allocations are inconsistent: the fronts and transfer must be treated together. The radii depend on photon hardness, abundance, electron number density and recombination rates, not on thresholds alone. The sharp-front equations are an approximation; they do not assert three exact independently conserved photon budgets for a real mixed nebula.
For a nonuniform nebula, an observed recombination line with effective recombination coefficient has luminosity
If gas temperature and density permit a known approximately constant ratio , the integrated line gives
The unknown density distribution cancels because the same recombination integral appears in both quantities. For a stellar continuum measurement at a nonionizing frequency ,
When the inferred recombination budget traces the full appropriate ionizing-photon supply, define and obtain the Zanstra method equation
Both and the distance cancel. Hydrogen and helium recombination lines probe different hardness thresholds, so their ratios and the stellar continuum constrain and test the assumed spectrum and trapping. For a mixed nebula, use the corresponding transfer-corrected photon budgets. temperature gradients require emission-weighted coefficients; dust, leakage, density-bounded helium zones or departures from a blackbody can produce discrepant hydrogen and helium Zanstra estimates. Line fluxes alone without an appropriate budget or continuum normalization need not uniquely determine .
Finally, locate the derived and luminosity on the Hertzsprung-Russell diagram. The central stars are hot exposed cores following post-asymptotic giant branch evolution: envelope loss creates the planetary nebula, contraction heats the core, and continued shell burning can support an approximately constant luminosity as the star moves leftward. Later the luminosity falls and the remnant approaches the white dwarf cooling sequence. Comparing central-star positions with stellar evolution tracks constrains core masses and evolutionary ages; comparison with nebular expansion ages tests whether the core heats rapidly enough to ionize the expelled envelope. The central-star locus traces the transition from an envelope-losing giant to a compact stellar remnant.
The first argument is the energy budget. The gravitational energy reservoir is of order , so the Kelvin-Helmholtz cooling time is . For the Sun, , and give years. This is much shorter than the approximately -year solar-system age. Chemical binding energies are only of order electronvolts per atom and provide an even smaller reservoir. Thus neither ordinary chemical energy nor continuing gravitational contraction can supply the long-term solar luminosity.
Stellar nuclear fusion provides an adequate reservoir. Conversion of four hydrogen atoms into one helium-4 atom releases about , corresponding to a mass defect of about of the initial mass; some energy leaves as neutrinos. If a fraction of stellar mass with hydrogen fraction is available for core burning, the energy scale is . Taking and supplies roughly , enough for several billion years at the solar luminosity, consistent with a stellar nuclear timescale of order years. The core temperature and density inferred from stellar structure permit the proton–proton chain; hotter stellar cores favor the CNO cycle.
The most direct evidence is the detection of solar neutrinos. They are created in the weak reactions of the proton–proton chain and can escape from the core, unlike the photons whose transport through the interior is slow. Measurements accounting for neutrino oscillation recover a total neutrino flux consistent with nuclear-powered solar models. Their spectra and reaction-dependent fluxes test the core reactions rather than merely the surface effective temperature.
There is also independent evidence from stellar evolution: the main sequence is a long-lived core-hydrogen-burning phase, more massive luminous stars consume fuel more quickly, and the main-sequence turnoff of a star cluster records its age. Subsequent evolution and the production of heavier elements fit further nuclear reactions. Agreement among these observations, nuclear energy release and stellar-model lifetimes is much stronger than the energy-budget argument alone. The combination of a sufficient nuclear reservoir, measured core neutrinos and the observed evolutionary sequence identifies nuclear reactions as the main long-term power source of ordinary stars.
After complete ionization, each hydrogen nucleus supplies one ion and one Electron, while each helium-4 nucleus supplies one ion and two Electrons. Neglecting the small metal mass fraction and putting , the mean molecular weight therefore satisfies
Use stellar homology to compare models with common dimensionless profiles. The stellar mass conservation equation, stellar hydrostatic equation and ideal gas law give
The radiative diffusion in a star equation then gives . Inserting the specified opacity yields . Independently, integrating the stellar energy-generation rate over the mass gives . Equating these in stellar thermal equilibrium gives the radiative homology with fifth-power hydrogen burning and inverse-cubic opacity:
At fixed solar mass these reduce to and . The stellar homology approximation, rather than merely the local power laws, is what permits these relations between whole stars.
Equal luminosities imply . Thus
so to one significant figure. Retaining the unrounded hydrogen mass fraction when comparing central temperatures gives
The first model has the higher central temperature.
For homogeneous hydrogen burning at fixed , fuel conservation gives . Since and , the homogeneous fuel-depletion luminosity feedback obeys
Integration with the initial mean molecular weight gives
This idealized stellar evolution law stops when ; it does not predict a physically attainable infinite luminosity. Indeed at fuel exhaustion, giving , strictly before the formal singularity.
A star of roughly solar stellar composition begins on the zero-age main sequence with a hot convective core, a mostly radiative envelope, and hydrogen burning dominated by the CNO cycle. The strong temperature dependence of the CNO cycle concentrates the stellar energy-generation rate near the centre, helping produce the convective core. Core mixing replenishes hydrogen within that region; it does not homogenize the entire star. The nuclear fuel supply gives a stellar nuclear timescale of order years. For a numerical comparison, Table 2 of Ekström et al.'s solar-metallicity models gives core hydrogen burning lifetimes of about and million years for nonrotating and rotating models, respectively.
After central hydrogen exhaustion, an inert helium core contracts and heats through Kelvin-Helmholtz contraction, while a hydrogen-burning shell supplies the growing luminosity. The envelope expands and cools: the track crosses the Hertzsprung gap toward the red giant region. Its deepening convective envelope causes first dredge-up, bringing material altered by hydrogen burning to the surface. The rapid structural crossing has a thermal timescale, roughly – years for a star of this mass, rather than the full stellar nuclear timescale; an estimate follows from the Kelvin-Helmholtz cooling time with the evolving radius and luminosity.
The helium core reaches ignition before becoming strongly supported by electron degeneracy pressure. Thus a normal star starts core helium burning without the helium flash characteristic of lower-mass stars. The Triple-alpha process produces carbon-12, and carbon-12 alpha capture produces oxygen-16; these reactions build a carbon-oxygen core. There is a convective core and continued hydrogen-burning shell activity. Core helium burning lasts order years: the same model table gives approximately and million years. The changing core and envelope structure can produce a blue loop in the Hertzsprung-Russell diagram; its extent depends on stellar composition, mixing and stellar rotation.
After central helium exhaustion, the carbon-oxygen core contracts and becomes increasingly supported by electron degeneracy pressure. A helium-burning shell forms above it, beneath the hydrogen-burning shell. The star returns to a cool, luminous Asymptotic giant branch configuration. On this early branch, helium-burning shell activity is an important energy source; the overlying hydrogen-burning shell can weaken during readjustment and resume later. At this mass, a deepening convective envelope can cause second dredge-up. This shell-supported evolution to the first AGB thermal pulse takes order years, with model-dependent duration. As the helium-burning shell becomes sufficiently thin, the thin-shell instability ends steady shell burning and marks the onset of AGB thermal pulses.
Figure 1.
Schematic evolution of a five-solar-mass star to the first AGB thermal pulse
. Details.
The original sketch shows the direction of stellar evolution, with higher effective temperature to the left. It illustrates a possible blue loop; its coordinates and loop shape are schematic, rather than a computed stellar model. The stages before the first AGB thermal pulse are the requested track; the later remnant evolution is described below.
In the usual final evolution, stellar winds and an AGB superwind remove the envelope before the carbon-oxygen core can ignite substantial carbon burning. Once shell fuel is exhausted, the exposed core contracts and heats, potentially ionizing the expelled gas as a planetary nebula, and finally becomes a carbon-oxygen white dwarf that loses its stored thermal energy by cooling. The relevant competition is between the AGB envelope-loss and core-growth timescales, not just the total initial mass.
Much higher late-stage mass loss strips the envelope sooner, shortens the Asymptotic giant branch phase, and allows less core growth, usually producing a smaller carbon-oxygen white dwarf. Much lower mass loss prolongs shell burning and permits greater core growth. In an extreme hypothetical case, retention of the envelope could let a degenerate carbon-oxygen core approach the Chandrasekhar mass and undergo explosive carbon burning: the proposed type 1.5 supernova channel. This is not an inevitable fate for a star. Doherty et al.'s evolutionary calculations emphasize the uncertainty of that channel and find that their adopted mass loss removes carbon-oxygen envelopes before such ignition. If carbon burning instead proceeds without disruption, the possible further evolution depends on the resulting core and continued envelope retention; lowering mass loss alone does not establish an electron-capture supernova outcome.
Treat the stars as point masses for the orbital motion, neglecting their spin angular momenta. Let . Their distances from the centre of mass are and . The sum of their orbital angular momenta is
Thus the circular-binary orbital angular momentum is
For conservative binary mass transfer, both and are fixed. In terms of the binary mass ratio , and , hence
Mass conservation also gives , so
Differentiating at fixed gives . Therefore the stipulated Roche lobe law has the Roche-lobe radius response exponent
Define the overfill by . On a time scale short compared with the slow change of , the donor's response is . Thus
During mass loss . If , mass loss increases the overfill and therefore drives more mass loss: this is positive feedback, giving more rapid Roche-lobe overflow. If , mass loss decreases the overfill and can balance the slow expansion due to stellar evolution. More explicitly, with , contact requires
A slow negative transfer rate can satisfy this only on the stable side . The rapid-transfer condition is therefore
Its time scale depends on the physical donor response represented by ; it need not be a dynamical time scale if is a thermal-equilibrium mass-radius exponent.
For Roche-lobe overflow feedback with a constant wind torque factor, retain mass conservation to the stated approximation but allow the prescribed angular-momentum loss. Since , differentiation gives
The new Roche-lobe radius response exponent is
The same overfill calculation gives
A loss of angular momentum has because , and it lowers the critical ratio by making the orbit contract more strongly.
For the specified solar-mass donor exponent , the answer is
Without the wind torque this becomes , so donors slightly less massive than their companions can already lie on the rapid-transfer side of this model. For the lower threshold is , and for every positive satisfies the strict inequality. There is no finite upper bound on from this criterion. Since the paper does not specify , the wind case has a conditional range rather than a unique numerical cutoff.
Write and for the number densities of nickel-56 and iron-54. Both nuclei contain 28 neutrons, while their proton numbers differ by two. Dividing the positive-exponential nuclear statistical equilibrium formulas printed in the PDF therefore gives the nickel–iron nuclear equilibrium ratio
The neutron-density factors cancel. contains the isotope-dependent statistical and translational mass factors suppressed in the supplied proportionality; it is constant in this simplified calculation. The TeX conversion incorrectly reverses the exponential sign, so the PDF sign is essential here.
In the specified restricted composition, equality of total proton and neutron numbers means
Consequently
Differentiation, with the suppressed prefactors held constant, gives
Using the paper's supplied signed energy difference produces
With that negative value, tends to zero at both temperature endpoints, increases below this temperature and decreases above it. At the stationary point , so it is indeed a maximum for the printed toy data.
There is a genuine physical-data qualification: the stated negative energy difference is inconsistent with interpreting as the conventional positive total nuclear binding energy of these actual isotopes. Using atomic masses, whose electron masses cancel in this reaction,
The masses used are , and , from NIST hydrogen masses, NIST iron masses and NIST nickel masses. Small electronic binding corrections do not change the sign. With a positive , the simplified above is strictly decreasing for and has no finite maximum. Thus the numerical maximum follows from the supplied negative parameter, not from consistent physical binding energies; changing the TeX sign to recover a maximum would not repair the PDF's data.
Nuclear statistical equilibrium is appropriate when the relevant forward and inverse strong or electromagnetic reactions are fast compared with changes of temperature, density and composition. It is useful in sufficiently hot late silicon-burning material and in the hottest explosive burning regions of massive stars. Some burning zones reach only quasi-equilibrium, and cooling or expansion eventually freezes the abundances out; an equilibrium expression cannot by itself predict the full ashes when these timescale conditions fail. Temperature, density, conserved baryon number and the electron fraction in stellar matter must all be specified. Negligible free-particle contributions in bulk bookkeeping are not the same as exactly zero free densities in the equilibrium formula.
Strong and electromagnetic reactions preserve total charge and baryon number, so they conserve the proton-to-neutron ratio. Weak reactions need not: electron capture and beta-plus decay convert protons into neutrons and lower it. Such conversions occur in hydrogen-burning networks and in the nitrogen-to-neon chain during helium burning. At high density in advanced stellar cores and core collapse, electron capture becomes particularly important; neutrino emission also affects energy and lepton loss. Beta-minus decay changes the ratio in the opposite direction and can compete during later relaxation. Nuclear equilibrium for strong reactions does not automatically imply weak equilibrium or a fixed proton-to-neutron ratio throughout all stellar evolution.
Let and let denote present stellar age, rather than time measured from the galaxy's birth. With a fixed initial mass function, constant star formation rate also gives a constant number of births per unit time, as in the birth-age distribution under constant star formation. Counting the supplied stellar evolution model's remnants as part of the population, ages have a uniform distribution on . Integrating the constant birth rate therefore gives
Outside this interval the cumulative fraction is zero or one as appropriate.
Write for the birth-mass parameter and . The normalization cancels when integrating the initial mass function:
Thus the mass-tail fraction is
It is one at and below the lower cutoff. The normalization constant has no effect on any number fraction.
For evolutionary-state selection by stellar lifetimes, a red giant has completed its main sequence lifetime but not its giant lifetime. Its stellar lifetime regions in a mass-age diagram are bounded by
The lower boundary meets at and the upper boundary meets it at . A white dwarf lies above the upper lifetime curve, within the same age interval.
Figure 1.
Red-giant and white-dwarf regions in the birth-mass–age plane and their triangular images in uniform mass-tail–age coordinates
. Details.
The probability integral transform, applied to the decreasing mass-tail coordinate, makes uniform on : . The independent birth age gives uniform on , so number fractions are areas in a unit square. Since , the red giant region becomes , a triangle with vertices , and . Its area is . The white dwarf region is a triangle of area . Hence the present individual-star fractions are
The rest are on the main sequence under the given toy lifetime model.
For the binary stars, draw two independent mass-tail coordinates , but only one shared age coordinate : components born together are coeval. The coeval binary population therefore occupies a uniform unit cube . At a fixed age , the state probabilities for either component are
Their evolutionary states have conditional independence given . They generally do not have unconditional independence, because sharing an age correlates their states. The binary fractions below are slice areas averaged over .
For this stellar population, let denote initial mass in units of the solar mass, and let denote present stellar age. Constant formation of equal numbers of stars per unit time makes the stellar age a uniform distribution on Gyr, so is uniform on . The normalized initial mass function has probability density function for , since . Consequently
These formulae have the stated age and mass domains; outside them the relevant cumulative fractions saturate at zero or one. The probability integral transform makes uniform on : gives . The time-independent initial mass function and constant number formation rate give independence of and . All fractions here count objects, including white dwarfs, using the stipulated stellar evolution law.
A red giant has
Thus the mass boundaries are and for positive . Equivalently, for the red giant region runs from to the age cap , and for it runs from to . There are no red giants with . Boundaries have zero probability and their endpoint convention does not affect the fractions. In the uniform square the red giant region is , with triangle vertices , and . The white dwarf region is the triangle .
Figure 1.
Mass-age regions and their uniform-coordinate images
. The right panel magnifies the evolved part of the unit square. All of the remaining region at larger X is main sequence. Details.
At fixed , define the conditional probabilities of a red giant, white dwarf and main sequence star by , and . Their interval widths are
Integration over the uniform distribution of age gives the individual-star fractions
The systems form a coeval binary population: the two binary star components have the same age. Their masses are independent, so has uniform probability density function on and their states have conditional independence given . Unconditional independence of their states would be incorrect: older systems make both evolved states more likely. The law of total probability now gives
The subtraction removes the double counting of systems containing two red giants.
Assume an isolated star of initial mass five solar masses of approximately solar stellar metallicity. Rotation, convective overshooting and mass loss change numerical ages and the extent of a blue loop, so the ages below are estimates. Take zero age at the zero-age main sequence; the pre-main-sequence star phase adds a much shorter contraction time.
The original schematic Hertzsprung-Russell diagram marks the stages discussed below. The strip denotes the approximate Cepheid instability strip; the track is illustrative, rather than a computed stellar model.
Figure 1.
Schematic stellar evolution of an initially five-solar-mass star
. Details.
The main sequence lasts roughly – years. CNO cycle hydrogen burning is concentrated in a mixed convective core. Its declining hydrogen mass fraction and rising helium mass fraction are nearly uniform there, while the retreating core leaves a composition gradient. The unprocessed envelope remains hydrogen rich. At the terminal-age main sequence, core hydrogen is exhausted and a hydrogen-burning shell takes over.
The inert helium core grows through shell burning. The Schönberg-Chandrasekhar limit is approximately
Beyond this limit an isothermal nondegenerate core cannot remain in thermal equilibrium with its envelope. Core contraction and envelope expansion carry the star across the Hertzsprung gap. An initial slower shell-burning interval can precede the rapid crossing; the crossing itself is governed by Kelvin-Helmholtz contraction, with a representative – years. First dredge-up then mixes hydrogen-processed material into the convective envelope, lowering its hydrogen mass fraction and increasing helium and nitrogen, while leaving a composition discontinuity where the envelope later retreats.
At an age still of order years, core helium burning begins quietly: the core is nondegenerate, so there is no helium flash. The Triple-alpha process and subsequent alpha capture build a carbon-oxygen core. Core helium burning lasts roughly – years in representative solar-composition models. Pols's stellar-evolution notes illustrate the substantial dependence of these lifetimes on convective overshooting.
During a blue loop the envelope contracts and effective temperature increases, before the star returns redward. The loop depends on core and envelope structure and the hydrogen discontinuity left by first dredge-up. A loop reaching the Cepheid instability strip produces two further crossings, blueward and redward, in addition to the earlier rapid crossing of the Hertzsprung gap. Cepheid variables pulsate through the opacity mechanism, involving helium ionization. Neither loop extent nor all three strip crossings are guaranteed for every composition or mixing prescription. The hydrogen-burning shell continues moving outward in enclosed mass while central helium is depleted. Lattanzio's five-solar-mass tutorial shows these composition changes.
After helium exhaustion, at an age roughly – years, the star ascends the Asymptotic giant branch. An inert carbon-oxygen core is surrounded by helium-rich material and a hydrogen-rich envelope. Second dredge-up mixes helium and hydrogen burning products into the envelope and reduces the hydrogen-exhausted core mass. Subsequently a helium-burning shell and hydrogen-burning shell alternate in importance. A thermal pulse of an asymptotic-giant-branch star causes intershell convection, expansion and temporary quenching of hydrogen burning; third dredge-up can carry carbon and slow neutron-capture process products outward. Between pulses hydrogen burning rebuilds the helium layer. These mechanisms are illustrated in Lattanzio's AGB tutorial.
The following original stellar composition profile sketches distinguish hydrogen exhaustion from helium exhaustion; their mass boundaries are illustrative.
Figure 2.
Schematic internal composition profiles during stellar evolution
. Details.
The final giant phases add at most a few million years at this level of accuracy. Strong stellar winds remove the envelope; the hot remnant can illuminate a planetary nebula and then cool as a carbon-oxygen white dwarf. A typical remnant is of order one solar mass. The ordinary isolated case does not reach sustained carbon burning and core collapse. Thus the usual endpoint is a carbon-oxygen white dwarf, at a total age of order years, with its subsequent cooling age added separately.
Stellar population Created 2026-10-06 Updated 2026-10-14
A stellar population is an ensemble of stars described by its formation times, initial masses, composition and present stellar evolution states. Number fractions must be distinguished from mass-weighted fractions. Shared formation times can correlate the states of components of a binary star even when their initial masses have independence.