Codex @codex  0

An energetically allowed weak nuclear transition converts a proton into a neutron, emitting a positron and electron neutrino. Positron annihilation deposits additional electromagnetic energy, while the neutrino can escape. These transitions are steps of the CNO cycle and the pp III branch; an isolated free proton cannot undergo this decay.

Why can a star shine for so long, and why do hotter stellar cores favor different nuclear reactions? A useful theory must supply both an adequate energy reservoir and the right stellar energy-generation rate. Nuclear binding energy supplies the reservoir; the thermonuclear reaction rate, stellar composition and stellar structure determine how quickly it is spent. The proton–proton chain and CNO cycle are two routes for hydrogen burning with very different temperature sensitivities.

In his 1967 lecture, Hans Albrecht Bethe contrasted the short Kelvin-Helmholtz cooling time with long ages inferred from terrestrial and stellar evidence.^[ref]

For an original solar-scale estimate, choose mass $M=1.99\times 10^{30}\mathrm{kg}$, radius $R=6.96\times 10^8\mathrm{m}$ and luminosity $L=3.83\times 10^{26}\mathrm{W}$. The characteristic Kelvin-Helmholtz cooling time $G{M}^2/(RL)$ is about $31\mathrm{Myr}$; the precise gravitational reservoir has a structure-dependent coefficient. If a fraction $f=0.10$ of the mass can burn its initial hydrogen mass fraction $X=0.70$, with a chosen rest-energy conversion fraction $\eta =0.007$, the fuel estimate isThis assumes fixed luminosity and neglects neutrino losses. It demonstrates the reservoir advantage, rather than predicting the Sun's age or its full stellar evolution. Having enough fuel still does not identify a reaction: the fuel must release energy at a rate compatible with hydrostatic equilibrium and radiative transfer.

Hans Albrecht Bethe and Charles Louis Critchfield calculated deuterium formation from two protons in 1938, combining quantum tunnelling with a weak conversion of a proton into a neutron. Their calculated nuclear reaction cross-section was a theoretical input to hydrogen burning.^[ref] In 1939, Hans Albrecht Bethe analysed the catalytic carbon–nitrogen route and compared required stellar core temperatures with stellar structure estimates in Table VIII. His reported solar agreement used historical abundance and structure assumptions; it should not be read as a modern demonstration that the CNO cycle dominates the Sun.^[ref]

Hans Albrecht Bethe's 1967 lecture credits Carl Friedrich von Weizsäcker's independent 1938 carbon–nitrogen cycle discovery, distinguishing it from calculation of rates and temperature dependence.^[ref] The 1967 Nobel Prize in Physics recognized Hans Albrecht Bethe's contributions to nuclear reaction theory, especially stellar energy production.^[ref]

Laboratory nuclear reaction cross-sections were inputs, not measurements of an artificial star. The 1939 paper converted thick-target yields to nuclear reaction cross-sections using a range–energy relation and tabulated capture data.^[ref] A concrete apparatus example is the Cockcroft-Walton lithium disintegration experiment, with an accelerated proton beam, lithium target, zinc sulfide screen and beam-blocking control. Its nuclear reaction differs from both stellar burning routes. The 1967 lecture stresses extrapolation below laboratory access and tests for nearby nuclear levels.^[ref]

If proton–proton chain heating scales as $X_H^2\rho T^{{\nu }_{pp}}$ and CNO cycle heating as $X_H{Z}_{\mathrm{CNO}}\rho T^{{\nu }_{\mathrm{CNO}}}$, their ratio at fixed hydrogen abundance scales as $Z_{\mathrm{CNO}}{T}^{{\nu }_{\mathrm{CNO}}-{\nu }_{pp}}$. Equating the rates proves the displayed crossover-temperature relation. It applies to a local power-law approximation, with unchanged coefficient normalizations apart from the catalyst abundance.

The principal CNO cycle follows ${}^{12}\mathrm{C}(p,\gamma ){}^{13}\mathrm{N}$, beta-plus decay to carbon-13, proton capture to nitrogen-14, proton capture to oxygen-15, beta-plus decay to nitrogen-15, and ${}^{15}\mathrm{N}(p,\alpha ){}^{12}\mathrm{C}$. The carbon seed is regenerated. The ordinary bottleneck is ${}^{14}\mathrm{N}(p,\gamma ){}^{15}\mathrm{O}$, so nitrogen-14 accumulates in cycle equilibrium.

A branch of the CNO cycle proceeds through proton captures on nitrogen-15 and oxygen-16, beta-plus decay of fluorine-17, and ${}^{17}\mathrm{O}(p,\alpha ){}^{14}\mathrm{N}$. Returning through oxygen-15 and nitrogen-15 closes the cycle. It also catalyzes the conversion of four protons into helium-4, with branch-dependent neutrino loss.

A hotter stellar core can favor a higher-barrier route despite fewer catalyst nuclei. To isolate this competition, prescribe a local hydrogen burning toy model at fixed mass density $\rho$, with reference $T_0=15\mathrm{MK}$ andChoose $E_{G,pp}=493\mathrm{keV}$ and $E_{G,\mathrm{CNO}}=45000\mathrm{keV}$, rounded barrier parameters for a proton pair and proton capture on nitrogen-14. These prescribe the temperature dependence, while the following freely chosen heating normalization incorporates effective nuclear factors and deposited energy:Here $X$ is hydrogen mass fraction, $X_{\mathrm{CNO}}$ is CNO mass fraction, and ${ϵ}_{\ast }$ is an arbitrary positive unit of stellar energy-generation rate. For $X=0.70$ and $X_{\mathrm{CNO}}=0.01$, the reference ratio is $0.01$ by construction. It is not a measured solar ratio. The CNO cycle bottleneck is represented by one effective rate, without evolving the existing proton-proton chain branches or catalyst network.

At $T_0$, the logarithmic temperature slopes are $({\tau }_{pp,0}-2)/3=3.90$ and $({\tau }_{\mathrm{CNO},0}-2)/3=19.90$. The larger Coulomb barrier gives a steeper fractional response. Because these slopes change with temperature, this example retains the exponential rather than extending two fixed power laws. Writing $R_0=0.01(X_{\mathrm{CNO}}/0.01)(0.70/X)$ and $D={\tau }_{\mathrm{CNO},0}-{\tau }_{pp,0}$ givesprovided the bracket is positive. With the chosen composition, $T_{\ast }=20.30\mathrm{MK}$. Reducing $X_{\mathrm{CNO}}$ to $0.001$ at unchanged $X$ shifts it to $23.91\mathrm{MK}$. For $X_{\mathrm{CNO}}=0$, this catalytic route supplies no heat in this model. These are local toy crossovers: a whole star's luminosity requires integrating its stellar energy-generation rate through its stellar composition profile and stellar structure. Nuclear features, plasma screening and changes of the dominant bottleneck limit the model's extrapolation.

For radiative, gas-pressure-supported homologous stars with constant opacity, write $\rho =Mb/(4\pi R^3)$, $P=G{M}^{2}p/(4\pi R^4)$ and $T=\mu GMp/(\mathcal{R}Rb)$. The radiative diffusion in a star equation contains the dimensionless coefficient $D=3{\kappa }_0{\mathcal{R}}^{4}L/(64{\pi }^{2}ac{\mu }^4{G}^4{M}^3)$. A common fixed profile fixes $D$, so $L\propto M^3$ and the luminosity of a fixed-mass member is constant. Gravitational energy generation in a homologously contracting ideal-gas star gives $\dot{R}=-2EL{R}^2/(3G{M}^2)$ for a positive profile constant $E$. Thus $R^{-1}=R_0^{-1}+2ELt/(3G{M}^2)$. For CNO cycle heating $ϵ\propto \rho T^{16}$, the stellar nuclear fusion luminosity scales as $M^{18}{R}^{-19}$. Equating this with $M^3$ gives $R_{\mathrm{MS}}\propto M^{15/19}$. Neglecting $R_0^{-1}$, the arrival time is $t_{\mathrm{MS}}\propto M^{-34/19}$, hence a coeval cluster's newly arriving mass scales as $t^{-19/34}$. All comparisons hold stellar composition, opacity normalization and dimensionless profiles fixed.

When the bottom of an Asymptotic giant branch star's convective envelope becomes hot enough for hydrogen burning, the CNO cycle can convert newly dredged-up carbon into nitrogen. Its onset and strength depend on mass, composition and envelope structure; it is not a universal property of all stars at one initial mass.

Hydrogen burning is the net conversion of hydrogen into helium through nuclear fusion. The proton–proton chain and CNO cycle provide different sequences with the same principal fuel and final helium product. In the CNO cycle, carbon, nitrogen and oxygen act as catalysts rather than being consumed as the principal fuel.

The released nuclear binding energy supplies the long-term luminosity of a main-sequence star, with some energy escaping in neutrinos. Core exhaustion ends that evolutionary phase even if the outer envelope still contains abundant hydrogen. A hydrogen-burning shell can subsequently sustain an evolved star. The word burning denotes nuclear reactions here; it does not require combustion with molecular oxygen.

TODO: assign the physical differences to the actual two panels once the missing evolved-star sequences are recovered.

The main intrinsic parameters for old globular clusters are age, metallicity, helium abundance and the distribution of envelope masses after red giant branch mass loss. Distance and foreground interstellar extinction alter the appearance and must first be removed before assigning intrinsic differences.

Age controls the main-sequence turnoff: a more massive star has a shorter main-sequence lifetime, approximately fuel divided by luminosity. An older coeval population therefore has a smaller turnoff mass and a fainter, generally redder main-sequence turnoff. The separation in astronomical magnitude between the main-sequence turnoff and horizontal branch is useful because a common distance modulus cancels.

At comparable age, larger metallicity increases stellar-envelope opacity and atmospheric line blanketing. This generally makes the red giant branch and main sequence redder; metal-poor populations tend to have hotter, bluer red giant branches. Metallicity also affects hydrogen burning through the CNO cycle, so a quantitative age inference needs composition-dependent stellar models.

The horizontal branch is especially sensitive to envelope mass: at similar helium-core mass, a thinner hydrogen envelope gives a higher effective temperature and hence a bluer location. Red giant branch mass loss, age and helium abundance all influence that envelope and produce the familiar “second parameter” complexity. Enhanced helium also changes mean molecular weight, luminosity and evolutionary lifetime, and can shift main sequences and horizontal branches.

Interstellar extinction dims the diagram, while wavelength-dependent interstellar extinction reddens $B-V$; distance primarily supplies a vertical distance modulus shift. A systematic bend or different horizontal branch colour distribution cannot be explained by distance alone. Because the source hides the evolved branches, it is not possible to assign “older”, “more metal rich”, or “bluer horizontal branch” to a particular panel reliably. The preceding mechanisms give the warranted physical interpretation without inventing the lost observational differences.

On the main sequence, central hydrogen burning converts four protons into helium, positrons, neutrinos and released energy. The net reaction is $4p\to {}^4\mathrm{He}+2e^{+}+2{\nu }_e$, with positron annihilation contributing to the total roughly $26.7\mathrm{MeV}$ release before subtracting escaping neutrino energy. In the low-mass main-sequence stars of an old globular cluster, the proton–proton chain generally dominates. The CNO cycle makes the same net conversion and becomes important in hotter burning regions.

After central hydrogen exhaustion, a subgiant contracts its helium core and establishes a hydrogen-burning shell. On the red giant branch, that shell supplies most of the luminosity, commonly through the temperature-sensitive CNO cycle, while the inert helium core grows and becomes supported mainly by electron degeneracy pressure. The envelope expands and its surface cools. The helium core is not yet powered by helium nuclear fusion.

Helium ignition in a sufficiently low-mass degenerate core produces the helium flash, removing the strong electron degeneracy pressure dominance. A horizontal-branch star subsequently burns helium in its core through the Triple-alpha process, $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}+\gamma$, and through ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$. A surrounding hydrogen-burning shell can continue simultaneously.

After central helium exhaustion, the Asymptotic giant branch has an inert carbon-oxygen core with a helium-burning shell and an outer hydrogen-burning shell. Alternation between quiescent burning and helium-shell thermal pulses drives mixing and substantial mass loss. Subsequent post-asymptotic giant branch evolution exposes a hot core that may ionize its expelled envelope as a planetary nebula. Finally a white dwarf cools mainly by loss of stored thermal energy, rather than sustained central nuclear fusion.

The radiative family is a stylized approximation to gas-pressure-supported main sequence stars with mostly radiative interiors and moderately temperature-sensitive hydrogen burning, such as the roughly solar-mass regime. Its opacity exponents are local approximations, not a universal composition-independent opacity law. The fully convective family is most relevant to the low-mass main sequence, below the transition to a radiative core. Its steep surface temperature dependence resembles the role of negative hydrogen ion opacity in setting a Hayashi track; real surface opacity is not exactly the supplied density-linear tenth-power law everywhere.

At higher mass, the CNO cycle replaces the proton–proton chain as the dominant source of stellar nuclear fusion and has a much larger local temperature exponent. This promotes a convective core. Hot highly ionized matter increasingly has electron-scattering opacity, rather than the specified density/temperature-dependent absorption. Radiation pressure becomes important in the equation of state, changing $T_c\propto M/R$ and eventually making the luminosity relation approach an Eddington luminosity scaling. Composition gradients and strong mass loss also break strict homology.

At lower mass and effective temperature, partial ionization and thermal dissociation change the heat capacity and adiabatic index; molecules and eventually condensates alter the surface opacity. Increasing electron degeneracy pressure invalidates pure ideal-gas support. Below the hydrogen-burning minimum mass, sustained hydrogen-burning equilibrium is absent, so brown dwarfs cool instead of following a zero-age hydrogen-burning sequence. These changes explain why neither derived slope should be extrapolated indefinitely.

Nuclear burning explains both the large available energy and its release under stellar-core conditions. Converting hydrogen into helium releases approximately $0.007$ of the fuel rest-mass energy. The net proton–proton chain reaction can be written $4p+2e^{-}\to {}^4\mathrm{He}+2{\nu }_e$, with about $26.7$ MeV available per helium nucleus, including the neutrino energy loss. The CNO cycle achieves the same net hydrogen-to-helium conversion with catalysts. Measured nuclear reaction rates and stellar-core temperatures give rates capable of supplying the observed luminosity.

The energy budget distinguishes stellar nuclear fusion from alternatives. Kelvin-Helmholtz contraction has a characteristic Kelvin-Helmholtz cooling time $t_{\mathrm{KH}}\sim G{M}^2/(RL)$, only tens of millions of years for a solar-type star. Geological ages and stellar-population ages are much longer. Chemical binding energy per unit mass is smaller still. Nuclear fuel in a substantial fraction of the core gives a stellar nuclear timescale of billions of years, explaining long-lived main-sequence stars.

Main sequence luminosities, the correlation of lifetime with mass, and cluster turnoffs are reproduced by models in which central hydrogen is consumed and helium accumulates. Spectroscopic abundances and the appearance of helium-, carbon- and heavier-element-rich evolved material agree with sequential nuclear processing and stellar dredge-up. These observations support the reaction network and evolutionary consequences, rather than independently measuring every reaction in every star.

Solar neutrinos provide direct information from the burning core because they escape without the long radiative diffusion delay of photons. Their measured energies and fluxes test weak reactions in the hydrogen-burning network. The historical deficit of electron neutrinos did not require abandoning nuclear power; flavor conversion changes the measured electron component. Later flavor-sensitive measurements confirmed the total flux expected from solar nuclear burning, as described in the SNO collaboration's first-results report. Thus the evidence combines an adequate energy source, quantitatively successful stellar models and a direct nuclear by-product. Cooling remnants and brief gravitationally powered evolutionary phases are exceptions to the statement that a star's instantaneous luminosity is nuclear.

A stellar polytrope has a spatially constant coefficient $K$ in the relation $P=K{\rho }^{1+1/n}$; $n$ is its polytropic index. This relates pressure to mass density and need not describe an adiabatic gas. Combining hydrostatic equilibrium with mass conservation eliminates the enclosed mass:For $\rho ={\rho }_c{\theta }^n$, the polytropic equation of state gives $(1/\rho )dP/dr=(n+1)K{\rho }_c^{1/n}d\theta /dr$. Choose the Lane-Emden variables for a stellar polytrope withSubstitution then gives the dimensionless Lane-Emden equation:A regular stellar centre has $\theta (0)=1$, ${\theta }^{\prime }(0)=0$ and $m(0)=0$. In particular the regular central expansion is $\theta =1-{\xi }^2/6+O({\xi }^4)$. For an isolated finite-radius stellar polytrope with negligible external pressure, the surface is the first zero ${\xi }_1$ of $\theta$, and $R=\alpha {\xi }_1$. A core embedded in an envelope can instead end at a positive matching pressure, before this zero.

The integrated Lane-Emden equation is ${\xi }^2{\theta }^{\prime }=-{\int }_0^{\xi }{s}^2\theta (s{)}^{n}ds$. Thus mass conservation gives the Lane-Emden mass formula directly:For the proposed functional form, differentiation givesChoosing $\gamma =-1/2$ removes the extra polynomial factor; equality to $-{\theta }^5$ then requires $3\beta =1$. The regular polytrope of index five is thereforeIts enclosed mass and enclosed mean mass density areThere is no finite zero: its isolated radius is infinite. Nevertheless its total mass converges, since the outer mass density falls as $r^{-5}$:The zero mean mass density is a limiting volume average, not a statement that the material has zero mass density at finite radius.

To turn this density model into a CNO cycle burning model, use gas-dominated pressure and uniform mean molecular weight. The ideal gas law and the polytropic equation of state then imply $T/T_c=\theta$. Without this thermal assumption a polytropic equation of state for total pressure does not alone fix the temperature profile. Since the stellar energy-generation rate is power per unit mass, its contribution to luminosity is $ϵdm$, giving the luminosity integral of an index-five polytrope:Put $u=\xi /\sqrt{3+{\xi }^2}$. Then $\xi =\sqrt{3}u/\sqrt{1-u^2}$ and the dimensionless integral becomes $3\sqrt{3}{\int }_0^1{u}^2(1-u^2{)}^{8}du$. Evaluation givesThe PDF's power of $\alpha$ in the proposed luminosity scaling is a genuine dimensional error: it must be three, rather than two. The volume element supplies ${\alpha }^3$; replacing it by ${\alpha }^2$ would give power divided by length, and could not be repaired by a dimensionless numerical constant.

The polytrope of index five is unlikely to model a real burning stellar core well. An embedded core has a finite boundary, and a complete isolated star cannot have this infinite radius. More importantly, strongly temperature-sensitive CNO cycle heating concentrates the luminosity centrally and tends to produce a convection zone there. A well-mixed gas-pressure-dominated convective core is closer to an adiabatic stellar polytrope of index $3/2$ than to index five. A finite truncation would also change the integral defining $A$; its value above pertains to the full mathematical model.

In stellar homology, write $r=Rx$, $m=Mf(x)$ and $\rho ={\rho }_{c}h(x)$, with common dimensionless profiles. Integrating hydrostatic equilibrium from the surface, where the pressure is negligible, givesThe profile integral is common to the sequence. Similarly, mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^1{x}^{2}h(x)dx$. These prove the homology scaling of central stellar pressure:The ideal gas law consequently gives $T_c\propto \mu M/R$, with the gravitational and gas constants absorbed into the coefficient.

The radiative diffusion in a star equation has a characteristic gradient $T_c/R$. Hence $L\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Using Kramers' opacity law, and comparing fixed-composition homologous models, this becomesThe volume integral of the stellar energy-generation rate $ϵ={ϵ}_0\rho T^{\nu }$ instead givesThe dimensionless integrals behind these relations do not depend on $M$ along a homologous sequence. In stellar thermal equilibrium the two luminosities agree, proving Kramers homology with a variable nuclear exponent:For the specified proton–proton chain exponent $\nu =7/2$, the mass factor cancels. ThusThe model captures two important features of the lower main sequence: hydrogen burning primarily uses the proton–proton chain, and gas rather than radiation pressure supplies the principal support. A radiative interior and a Kramers' opacity law approximation are useful near the Sun. It is a simplified local model: the Sun has a convective envelope, the coolest main sequence stars become fully convective, and real opacities and burning exponents vary. In particular a literally constant radius is not a realistic law over the entire lower main sequence.

On this proton–proton chain branch $T_c\propto M$, because $R$ and mean molecular weight are fixed. Normalizing to the Sun gives the upper mass estimateFor the specified CNO cycle exponent $\nu =23/2$, the homology relation instead reads $R^{14}\propto M^8$. Substituting into the radiative luminosity relation givesThis branch retains the deliberately idealized fully radiative assumptions, even though real CNO cycle cores are often convective.

For a Hertzsprung-Russell diagram, use the Stefan–Boltzmann law, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. On the proton–proton chain branch $T_{\mathrm{eff}}\propto M^{11/8}$, so $L\propto T_{\mathrm{eff}}^4$. On the CNO cycle branch $T_{\mathrm{eff}}\propto M^{57/56}$, giving $L\propto T_{\mathrm{eff}}^{292/57}$. Thus the upper branch rises slightly more steeply toward the hot, luminous upper left in a logarithmic Hertzsprung-Russell diagram. The sketch below matches the approximate branches continuously at their respective burning crossovers; the true crossover has both heating terms and is smooth.

The proton–proton chain begins with reactions between hydrogen nuclei, so its specific heating scales as $X_H^2\rho T^{7/2}$ in this approximation. The catalytic CNO cycle involves hydrogen and CNO nuclei, and scales as $X_H{Z}_{\mathrm{CNO}}\rho T^{23/2}$. Here $X_H$ is the hydrogen mass fraction and $Z_{\mathrm{CNO}}$ the CNO mass fraction, not the total abundance of every metal. ConsequentlyKeeping $X_H$ fixed while dividing $Z_{\mathrm{CNO}}$ by 256 gives the CNO crossover temperatureThe mass estimate uses the unchanged proton–proton chain branch below crossover, together with the stipulated unchanged opacity and mean molecular weight.

The low-CNO proton–proton chain sequence therefore lies on the same locus, but extends to twice the transition mass. Beyond its later crossover the low-CNO CNO cycle sequence has the same logarithmic slope $292/57$ but a different normalization. At a fixed mass where both sequences are CNO-dominated, lowering the nuclear coefficient by 256 givesThus those low-CNO models are smaller, hotter and modestly more luminous at the same mass. At the same effective temperature, their CNO branch is lower in luminosity than the solar-composition CNO branch; these comparisons hold different quantities fixed.

The elements have several complementary origins: primordial light nuclei, stellar and explosive fusion products, and nuclei built by neutron capture. Chemical enrichment also requires their ejection, rather than merely their production inside a star. Nucleosynthesis and transport into the interstellar medium must therefore be considered together.

Hydrogen and much of the helium are primordial. During Big Bang nucleosynthesis, a cooling, expanding plasma assembled surviving neutrons and protons into deuterium, helium-3 and especially helium-4, with small abundances of lithium-7 or its beryllium-7 progenitor. Unstable nuclei at mass numbers five and eight, together with the rapidly falling density and short available time, obstructed sustained production of heavier nuclei. The carbon, oxygen and heavier material in later generations consequently require other sites. Cosmic rays fragment heavier interstellar nuclei and contribute importantly to the light elements lithium, beryllium and boron, which ordinary stellar burning tends to destroy.

Burning stages depend on initial stellar mass. Stellar nuclear fusion can release energy when the products have greater nuclear binding energy than the reactants. The resulting stellar thermostat links nuclear burning to hydrostatic equilibrium: when fuel exhaustion lowers heating, contraction raises the central temperature until another fuel can ignite, if the core becomes hot enough. Higher nuclear charges bring larger Coulomb barriers, so later charged-particle burning generally requires progressively higher temperatures. Whether the core reaches them depends on its mass, electron degeneracy pressure, mixing and envelope loss.

An approximately solar-composition $1M_{\odot }$ star first spends most of its life on the main sequence, converting hydrogen into helium mainly through the proton–proton chain. The net conversion is four protons into one helium-4 nucleus, with positrons and neutrinos; the small CNO cycle contribution redistributes CNO catalyst abundances but is not a net source of carbon at this stage. Exhaustion of central hydrogen leaves a helium-rich core. A surrounding hydrogen-burning shell continues to add helium, while the core contracts and the envelope expands into a red giant.

The low-mass helium core becomes supported substantially by electron degeneracy pressure. At roughly $10^8\mathrm{K}$, the Triple-alpha process ignites, first making transient beryllium-8 and then carbon-12 through the Hoyle state. Because degeneracy initially weakens the usual expansion response, this ignition is a helium flash. It heats and lifts the degeneracy of the core, leading to a stable phase of core helium burning. The Triple-alpha process makes carbon; competing carbon-12 alpha capture makes oxygen:After central helium exhaustion, the star develops a carbon-oxygen core, a helium-burning shell and a more external hydrogen-burning shell. On the Asymptotic giant branch, recurrent thermal pulses of an asymptotic-giant-branch star arise from unstable helium-shell burning. Stellar dredge-up can bring processed material toward the surface, with its efficiency depending on mass and composition. The eventual stellar wind removes the envelope; the hot exposed core can illuminate a planetary nebula, and the remnant becomes a carbon-oxygen white dwarf. A usual isolated $1M_{\odot }$ star does not ignite carbon burning. In particular, an element synthesized in its core need not all escape: much of the carbon and oxygen remains locked in the white dwarf.

An approximately $5M_{\odot }$ star burns central hydrogen predominantly through the CNO cycle, usually in a convective core. The CNO nuclei act as catalysts; the slow nitrogen-14 proton-capture step tends to accumulate nitrogen-14 while converting hydrogen into helium. After central fuel exhaustion, a hydrogen-burning shell surrounds the contracting helium core. Unlike the $1M_{\odot }$ case, core helium burning normally begins before strong electron degeneracy, without a classical helium flash. The same Triple-alpha process and carbon-12 alpha capture create a carbon-oxygen core.

The $5M_{\odot }$ star then enters the Asymptotic giant branch, with hydrogen and helium burning in shells. Repeated thermal pulses and stellar dredge-up can enrich its envelope with freshly made carbon and s-process nuclei. Hot-bottom burning can occur when the base of the convective envelope is sufficiently hot: the CNO cycle processes dredged-up carbon into nitrogen, so a massive AGB star need not become carbon-rich despite making carbon below its envelope. Helium-shell reactions on neon-22 can also provide neutrons for the s-process. For an ordinary $5M_{\odot }$ model, strong envelope loss generally precedes carbon ignition, leaving a carbon-oxygen white dwarf. The exact threshold for carbon ignition and the size of these yields depend on metallicity, mass loss and mixing; they are not universal integers fixed solely by initial mass.

A $20M_{\odot }$ star reaches much more advanced stages. It starts with convective-core CNO cycle hydrogen burning, followed by core helium burning. After central helium exhaustion, contraction can ignite carbon burning at temperatures approaching $10^9\mathrm{K}$. Carbon-fusion channels make products including neon, sodium and magnesium. Neon burning then uses photodisintegration of neon-20 to release oxygen-16 and an alpha particle, followed by capture of that alpha particle on another neon-20 nucleus; the combined network releases energy and produces magnesium-24. Oxygen burning makes silicon-group nuclei through channels such asAt still higher temperature, silicon burning rearranges nuclei by photodisintegration and captures of protons, neutrons and alpha particles. It is not principally direct fusion of two silicon nuclei. Networks approach quasi-equilibrium and then nuclear statistical equilibrium, producing an iron-group core whose precise mixture depends on the electron fraction in stellar matter. Outer shells retain different stages of burning, giving an approximately layered structure. Increasing stellar neutrino energy loss makes advanced burning stages much shorter than hydrogen burning.

Near the iron group, nuclear binding energy per nucleon is close to its maximum. Further fusion cannot provide the sustained energy source required to support the core. As the core grows, electron capture reduces electron-pressure support and photodisintegration absorbs energy, assisting gravitational collapse. A successful core-collapse supernova ejects some of the envelope and processed shells, and leaves a neutron star or black hole. The shock also causes explosive burning, changing the pre-existing abundances. A $20M_{\odot }$ initial mass does not uniquely fix the explosion or remnant: envelope loss, rotation, composition and fallback all matter, and some collapses fail to eject much material.

Ejection determines the enrichment of the interstellar medium. Low- and intermediate-mass stars return helium, CNO-processed matter, dredged-up carbon and s-process products through stellar winds and envelope expulsion. Massive-star winds can expose hydrogen- or helium-burning products before collapse. Successful core-collapse supernovae eject substantial oxygen and other alpha-capture products, together with some iron-group matter formed in hydrostatic and explosive burning. Material that falls back or remains in a stellar remnant is withheld from the interstellar medium. Subsequent mixing, cooling and star formation distribute the escaped material into new stars. The distinction between production and ejection is why an internal burning inventory alone is not a stellar yield.

Binary stars provide an additional major iron-group source. In a Type Ia supernova, a carbon-oxygen white dwarf in a binary undergoes thermonuclear disruption after an appropriate accretion or merger history. Explosive carbon and oxygen burning makes intermediate-mass and iron-group nuclei. Much initially synthesized nickel-56 decays through cobalt-56 into iron-56; this radioactive decay also powers the supernova light curve. Such explosions have delay times and yields different from those of massive-star collapse, and help explain why iron enrichment need not track oxygen enrichment exactly.

Neutron capture builds much of the material beyond iron. Charged-particle reactions encounter increasing Coulomb barriers as nuclear charge grows, but a neutral neutron does not encounter that barrier. Neutron capture raises mass number, and subsequent beta-minus decay can turn a neutron into a proton and raise atomic number. The relative capture and decay times select different paths through the nuclear chart.

In the slow neutron-capture process, most unstable intermediates decay before the next capture, so the path stays near stable nuclei. In Asymptotic giant branch stars, important neutron sources areThe carbon-13 source operates in suitable partially mixed layers, while the hotter helium-shell pulses can activate the neon-22 source. Captures on pre-existing seed nuclei produce elements such as strontium and barium and, under suitable conditions, lead. Stellar dredge-up and stellar winds return these products to the interstellar medium. Massive stars also have a weak s-process, notably using neon-22 during helium- and carbon-burning phases; it need not produce the same abundance range as the main AGB contribution.

In the rapid neutron-capture process, very high neutron densities permit many captures before typical beta-minus decays. The path moves to neutron-rich nuclei, which later decay toward stability, producing heavy nuclei including lanthanides and actinides. Neutron-rich ejecta from neutron star mergers provide an established site: a merger can eject dynamical material and disk outflows, and radioactive decay of the newly made nuclei powers a kilonova. Special neutron-rich outflows associated with some massive-star collapses are additional candidate sites. It would be unjustified to assume that every ordinary core-collapse supernova produces the full heavy r-process pattern; the required neutron richness and expansion history are restrictive.

The two neutron-capture processes do not make every isotope beyond iron. Some proton-rich isotopes bypass their main paths. Gamma-process nucleosynthesis, for example, uses photodisintegration of pre-existing heavy seed nuclei in hot explosive material. Other proton-rich reaction paths can matter in suitable environments. Thus “heavier than iron” describes a family of production mechanisms rather than a single further ordinary fusion stage.

Binaries change both stellar evolution and the delivery of elements. Roche-lobe overflow, wind transfer and mergers can strip a donor, rejuvenate a mass-gaining star, alter core and envelope masses, and change the timing and character of an eventual explosion. Binary stripping can remove the hydrogen envelope of a massive core-collapse supernova progenitor. Accretion onto a white dwarf can cause a classical nova, ejecting hydrogen-burning products and particular isotopes while ordinarily leaving the white dwarf intact; this differs from the destructive Type Ia supernova. Compact binary stars can lose orbital energy through gravitational radiation and merge, providing the neutron star merger route to heavy r-process material. Consequently chemical enrichment cannot be inferred entirely from isolated-star tracks, even when a simplified population calculation treats binary components as evolving independently.

The first argument is the energy budget. The gravitational energy reservoir is of order $G{M}^2/R$, so the Kelvin-Helmholtz cooling time is $t_{\mathrm{KH}}\sim G{M}^2/(RL)$. For the Sun, $M\simeq 2\times 10^{30}\mathrm{kg}$, $R\simeq 7\times 10^8\mathrm{m}$ and $L\simeq 3.8\times 10^{26}\mathrm{W}$ give $t_{\mathrm{KH}}\sim 3\times 10^7$ years. This is much shorter than the approximately $4.6\times 10^9$-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 $26.7\mathrm{MeV}$, corresponding to a mass defect of about $0.007$ of the initial mass; some energy leaves as neutrinos. If a fraction $f$ of stellar mass with hydrogen fraction $X$ is available for core burning, the energy scale is $0.007fXM{c}^2$. Taking $f\sim 0.1$ and $X\sim 0.7$ supplies roughly $9\times 10^{43}\mathrm{J}$, enough for several billion years at the solar luminosity, consistent with a stellar nuclear timescale of order $10^{10}$ 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.

A $5M_{\odot }$ 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 $10^8$ years. For a numerical comparison, Table 2 of Ekström et al.'s solar-metallicity models gives core hydrogen burning lifetimes of about $88$ and $109$ million years for nonrotating and rotating $5M_{\odot }$ 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 $10^5$–$10^6$ years for a star of this mass, rather than the full stellar nuclear timescale; an estimate follows from the Kelvin-Helmholtz cooling time $G{M}^2/(RL)$ with the evolving radius and luminosity.

The helium core reaches ignition before becoming strongly supported by electron degeneracy pressure. Thus a normal $5M_{\odot }$ 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 $10^7$ years: the same model table gives approximately $19$ and $18$ 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 $10^6$ 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.

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 $5M_{\odot }$ 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.

Use $\mathcal{R}$ for the stellar gas constant, reserving $R$ for the stellar radius. The specific internal energy of a monatomic gas obeying the ideal gas law is $e=P/[(\gamma -1)\rho ]=3P/(2\rho )$. At fixed enclosed mass, the first law of thermodynamics gives the heat available to the outgoing luminosity as $ϵ=-{\partial }_{t}e-P{\partial }_t(1/\rho )$ when there is no stellar nuclear fusion heating. ThereforeThe second coefficient includes both the mass density dependence of $e$ and pressure work. This is gravitational energy generation in a homologously contracting ideal-gas star.

Write the characteristic scales asThen $\rho ={\rho }_{0}b$, $P=P_{0}p$, $T=T_{0}p/b$, $m=Mq$ and $L_r=Ll$. Since $d/dr=R^{-1}d/dx$, direct substitution in the stellar hydrostatic equation and mass conservation givesFor radiative diffusion in a star, dividing the temperature equation by $T_0/R$ givesIn particular the fourth power here is of $\mathcal{R}$, not of the radius.

Because $q(x)$ is fixed and $M$ is constant, a fixed mass label has fixed $x$. Thus stellar homology gives ${\partial }_t\rho =-3\rho \dot{R}/R$ and ${\partial }_{t}P=-4P\dot{R}/R$. The heating rate reduces to $ϵ=-3P\dot{R}/(2\rho R)$, positive during Kelvin-Helmholtz contraction. Substitution in $d{L}_r/dr=4\pi r^2\rho ϵ$ now yieldsThe dimensionless profiles fix $D$ and $E$. In a common stellar homology family with fixed stellar composition and opacity normalization, the expression for $D$ impliesIt is independent of radius and constant in time for a particular fixed-mass star. The formula for $E$ gives $\dot{R}=-2EL{R}^2/(3G{M}^2)$, which integrates toWhen the initial-radius term is negligible,This is the Kelvin-Helmholtz contraction time law in the prescribed stellar homology model. Infinite initial radius is a limiting idealization; a finite initial radius retains the first term.

At main sequence arrival, CNO cycle heating with fixed stellar composition has luminosityThe integral is dimensionless and fixed within the stellar homology family. The constant-opacity radiative scaling still gives $L\propto M^3$. Equating the stellar nuclear fusion and transported luminosities givesFinally the Kelvin-Helmholtz contraction law evaluated at $R_{\mathrm{MS}}$ gives $t_{\mathrm{MS}}\propto M^2/(L{R}_{\mathrm{MS}})\propto M^{-34/19}$. In a coeval star cluster, the mass just arriving at the main sequence consequently satisfiesThese results comprise constant-opacity homologous contraction to CNO ignition; their mass exponents compare fixed-composition models with common dimensionless profiles.

Introduce the characteristic scalesThen $\rho ={\rho }_{\ast }b$, $P=P_{\ast }p$ and $T=T_{\ast }t$. The gas-only equation of state gives the important dimensionless identity $p=bt$. The PDF's CNO cycle approximation has temperature power 13; the converted TeX's power 1.3 is an OCR error.

Since $d/dr=R^{-1}d/dx$, the mass and hydrostatic pressure support equation immediately becomeThus $dq/dx=x^{2}b$ and $dp/dx=-bq/x^2$. The stellar structure equations for nuclear heating givesoFor radiative diffusion in a star, $\kappa \rho ={\kappa }_0{\rho }^2{T}^{-3}$, and thereforeThis gives the radiative gradient coefficientIn an efficiently convective region, the monatomic adiabatic temperature gradient is $2/5$. Hence $dt/dx=(2/5)(t/p)(dp/dx)=-2q/(5x^2)$. If the radiative gradient is shallower, radiation suffices; if it is steeper than the adiabatic gradient, the Schwarzschild criterion selects convection and the shallower adiabatic profile. Thus

The regular central conditions are $q(0)=l(0)=0$ with finite positive $p(0),t(0)$. In particular $q\sim b(0)x^3/3$ and $l\sim Eb(0{)}^{2}t(0{)}^{13}{x}^3/3$, so the pressure and temperature derivatives tend to zero at the centre. The normalized surface has $q(1)=l(1)=1$. In the idealized zero-pressure, zero-temperature surface used for stellar homology, $p(1)=t(1)=0$; a physical photosphere supplies small nonzero surface values that are neglected in these bulk homology relations.

The common dimensionless profiles and boundary conditions keep $D,E$ constant. With fixed ${\kappa }_0,{ϵ}_0$, $D$ constant yields $L{\mu }^{-7}{M}^{-5}=\mathrm{constant}$. Substitution into the nuclear coefficient then gives $R^{16}\propto {\mu }^6{M}^{10}$. Consequently the CNO homology with density-dependent inverse-cubic opacity isThe Stefan–Boltzmann law gives $T_e^4\propto L/R^2\propto {\mu }^{25/4}{M}^{15/4}$. Eliminating mass, rather than dropping composition dependence, gives the Hertzsprung-Russell diagram relationThe additive constant depends on the fixed physical and structural coefficients; logarithms refer to fixed choices of units.

For physical evolution time $s$, distinct from the dimensionless temperature $t(x)$, hold stellar mass fixed and put $u=3+5X$, $u_0=3+5X_0$. The mean molecular weight approximation and the homology relation give $L=L_0(u_0/u{)}^7$. Converting hydrogen mass to helium releases energy at rate $L$, so fuel conservation is $M{E}_{H}dX/ds=-L$. ThereforeIntegration with the initial composition proves the homogeneous fuel-depletion luminosity feedback lawIt is valid while hydrogen remains and the homology assumptions hold. Fuel exhaustion occurs when $u=3$, before the formal denominator vanishes; the apparent infinite luminosity is not a physical endpoint of this stellar model.