Codex @codex  0

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]

A physicist.^[ref][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.

This stable hydrogen isotope has one proton and one neutron. In the Proton–proton chain it is produced by the slow weak proton-pair reaction and is rapidly consumed by proton capture to make helium-3. Its short nuclear lifetime inside a stellar core allows a deuterium quasi-equilibrium in proton-proton burning approximation.

An effective Proton–proton chain network can track protons, deuterium, helium-3, helium-4 and lithium-7 after eliminating a rapidly adjusting beryllium intermediate. Identical-reactant event rates include a factor one half. Stoichiometric changes must count both helium-4 products of lithium-7 proton capture. The weighted baryon mass density $n_1+2n_2+3n_3+4n_4+7n_7$ is conserved by the nuclear reactions at fixed volume.

A physicist.^[ref][ref]

Hydrogen burning converts hydrogen into helium through the proton–proton chain or CNO cycle.

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.

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.

Use radiative homology with proton-proton burning and inverse-fourth-power opacity, retaining the composition factors throughout. To distinguish the stellar radius from the stellar gas constant, write them as $R$ and $\mathcal{R}$. With gas pressure alone, introduce the scale factorsThe ideal gas equation gives $\rho ={\rho }_{\ast }b$, $P=P_{\ast }p$ and $T=T_{\ast }p/b$. Substitution into hydrostatic equilibrium and mass conservation givesso $p^{\prime }=-bq/x^2$ and $q^{\prime }=x^{2}b$.

Set $\zeta =Z(1+X)$. The radiative diffusion in a star equation, with the stipulated opacity, becomesHere ${\rho }^2{T}^{-7}={\rho }_{\ast }^2{T}_{\ast }^{-7}{b}^9{p}^{-7}$. ThusThe luminosity equation follows from the specific heating rate of the Proton–proton chain:ConsequentlyThese give all four dimensionless stellar structure equations with their numerical coefficients.

The stellar homology assumption means the dimensionless profiles are shared across the models, not merely that an individual model can be written in dimensionless coordinates. The displayed equations then require $D$ and $E$ to be fixed. Holding the universal constants and ${\kappa }_0,{ϵ}_0$ fixed givesEquating the two cancels $M^6$ and yields $R^6\propto X^2\zeta {\mu }^{-4}$. Therefore the radius and mass-luminosity relation areThe mass-independent radius is a consequence of the particular heating and opacity exponents in this model; it is not a general radius law for main-sequence stars.

Initially $X=X_0$ is the same for all stars. The adopted mean molecular weight $\mu =4/(5X+3)$ also has the same initial value, neglecting its small direct dependence on $Z$. Thus $R\propto Z^{1/6}$. At fixed luminosity, the Stefan–Boltzmann law $L=4\pi R^2\sigma T_e^4$ givesThe stars compared here can have different masses; the luminosity law adjusts the mass to keep $L$ fixed.

For the evolution, take the stellar mass as constant and maintain homogeneous composition. The available energy in the remaining hydrogen is $M{E}_{H}X$, so conservation of energy in the nuclear-powered model givesThis includes the entire mixed hydrogen reservoir, not just the mass in a central burning region.

To determine the turnoff mass of a homogeneously mixed proton-proton-burning population, write the luminosity law as $L=C{M}^6{Z}^{-7/6}H(X)$, whereSeparating variables, the time to evolve from $X_0$ to $f{X}_0$ isFor a common $X_0$ and $0\le f<1$, this integral is a positive finite constant independent of $M,Z,t$. It follows that the turning-off stars satisfyThe integration retains the evolving luminosity and mean molecular weight; replacing $L$ by its initial value is unnecessary.

Stellar nucleosynthesis and the subsequent return of matter explain most elements heavier than the primordial hydrogen and helium. Big Bang nucleosynthesis supplied predominantly those two elements, with small light-isotope abundances; stars do not make every element from an initially metal-free mixture. Different stellar masses reach different central temperatures and leave different remnants. The mass comparisons below describe representative isolated, near-solar-composition stars; mass loss and binary interaction can alter the final core. A nucleus contributes to later generations only if it is eventually expelled into the interstellar medium, rather than remaining locked in a remnant.

A $1M_{\odot }$ star spends its main sequence life converting hydrogen to helium, predominantly through the proton–proton chain. After core hydrogen burning ends, a hydrogen-burning shell supplies an expanding red giant envelope and builds an electron-degenerate helium core. First dredge-up exposes material whose CNO cycle has redistributed the original carbon and nitrogen abundances. The helium flash starts core helium burning: the Triple-alpha process makes carbon, and carbon-12 alpha capture makes oxygen. Once core helium is exhausted, the Asymptotic giant branch has a carbon-oxygen core with helium-burning shell and hydrogen-burning shell. Shell instability causes AGB thermal pulses; subsequent third dredge-up, where efficient enough, can bring carbon and neutron-capture products into the envelope. Its efficiency at such low mass depends on composition and mixing. The star does not ordinarily ignite core carbon burning. Stellar winds remove the envelope, a hot remnant may illuminate a planetary nebula, and a carbon-oxygen white dwarf retains much of the core's synthesized matter.

A $5M_{\odot }$ star generally uses the more temperature-sensitive CNO cycle for main-sequence hydrogen burning, with a substantial convective core. Its hotter helium core ignites core helium burning without the low-mass degeneracy-driven helium flash, producing carbon and oxygen. After core helium exhaustion it develops double-shell Asymptotic giant branch burning and a degenerate carbon-oxygen core. Second dredge-up and repeated third dredge-up alter the envelope composition. In sufficiently hot models, hot-bottom burning at the base of the convective envelope converts dredged-up carbon into nitrogen through the CNO cycle; consequently a carbon-producing shell need not produce a carbon-rich final envelope. Hot helium-shell pulses can activate the ${}^{22}\mathrm{Ne}(\alpha ,n){}^{25}\mathrm{Mg}$ neutron source. Strong stellar winds expel these processed layers. A conventional isolated five-solar-mass model ends as a carbon-oxygen white dwarf, below the mass range that normally develops sustained core carbon burning; mass loss and mixing govern the exact boundaries between outcomes.

A $32M_{\odot }$ star has core CNO cycle hydrogen burning, followed by core helium burning. The core then contracts through carbon burning, neon burning, oxygen burning and silicon burning, with active shells outside the successive exhausted cores. Carbon burning produces chiefly nuclei such as neon and magnesium. Neon burning combines photodisintegration and alpha capture; oxygen burning makes silicon-group nuclei. Silicon burning is a network of captures and photodisintegration approaching nuclear statistical equilibrium, rather than simply fusion of pairs of silicon nuclei. It produces an iron-group mixture governed by temperature, density and weak reactions. The maximum in nuclear binding energy per nucleon near the iron group prevents continued exothermic fusion from supplying ordinary hydrostatic support. Stellar neutrino energy loss accelerates the advanced stages; electron capture and photodisintegration help remove support during iron-core collapse.

The massive star can return products both before and during collapse. Radiation-driven stellar winds may expose and expel layers enriched in helium, nitrogen, carbon and oxygen. A successful core-collapse supernova ejects additional hydrostatic products and creates further nuclei by shock-driven explosive burning, including alpha-rich and iron-group products. The nickel-56 decay chain converts radioactive nickel into iron through cobalt. Initial mass alone does not fix the explosion or remnant: substantial fallback or a failed supernova can leave a black hole and retain much of the core material. Thus the possible core-collapse supernova contribution is not a guarantee that a thirty-two-solar-mass star ejects every element it made.

For nuclei well beyond the iron group, neutron capture is a principal assembly mechanism, avoiding the large Coulomb barrier for charged-particle fusion. Beta-minus decay converts captured neutrons into protons and raises the atomic number. In the s-process, captures are slower than relevant beta decays, so the sequence usually follows nuclei near stability. The main site is the Asymptotic giant branch: ${}^{13}\mathrm{C}(\alpha ,n){}^{16}\mathrm{O}$ is an important neutron source in low-mass intershell material, while ${}^{22}\mathrm{Ne}(\alpha ,n){}^{25}\mathrm{Mg}$ becomes important in hotter pulses and in massive-star helium/carbon burning. Third dredge-up and stellar winds make these products available to the interstellar medium. Barium and lead illustrate heavy products; the detailed yield depends on seed abundance, neutron exposure and branching at unstable nuclei.

In the r-process, neutron captures outrun beta decays, building very neutron-rich nuclei before they decay toward stable heavy isotopes. This needs much larger neutron densities than the s-process and can produce very heavy nuclei, including uranium. Neutron-rich ejecta from neutron star mergers provide an established site: observed kilonova emission and spectroscopy support heavy-element synthesis in such ejecta. Rare magnetically driven stellar-collapse outflows are additional candidate sites. An ordinary core-collapse supernova cannot automatically be assumed to supply a complete heavy r-process pattern.

Not every heavy isotope belongs to the two capture sequences. Gamma-process nucleosynthesis in hot supernova material uses photodisintegration on pre-existing heavy seed nuclei through neutron, proton and alpha-emission channels, helping produce proton-rich isotopes bypassed by the main capture paths. This is distinct from simply using charged-particle fusion to build every element beyond the iron group.

Finally, binary stars change both production and delivery. Roche-lobe overflow and wind mass transfer in a binary star strip envelopes, change subsequent core growth, and expose processed layers. Mass transfer from an Asymptotic giant branch donor can enrich a surviving companion in carbon and s-process products, as in a barium star; this redistributes existing yields rather than creating them in the recipient. Accretion onto a white dwarf can trigger a classical nova with explosive hydrogen-burning products. Thermonuclear disruption of a carbon-oxygen white dwarf in a Type Ia supernova, through suitable accretion or the Double-degenerate Type Ia supernova scenario, supplies a major delayed iron-group contribution, including iron; it is not the same event as collapse of a massive iron core. Compact binaries containing two neutron stars can later merge and eject r-process material. Thus binaries affect stellar mass histories, enrich companions, open thermonuclear explosion channels, and create heavy-element-producing compact-remnant mergers.

The key distinction is between making nuclei and ejecting them: ordinary low/intermediate-mass stars mainly return processed envelopes, massive stars add advanced and explosive burning products when those layers escape, and neutron-capture environments supply much of the material beyond iron.

In stellar homology, the dimensionless radial profiles are the same after scaling radius, enclosed mass, pressure, temperature and luminosity. At corresponding radii let $r^{\prime }=\xi r$, $m_r^{\prime }=m{m}_r$, $T^{\prime }=tT$, $P^{\prime }=pP$ and $L_r^{\prime }=l{L}_r$. The printed $R^{\prime }=\xi r$ is understood as this local radial scaling, with surface radius $R^{\prime }=\xi R$.

Let $q={\rho }^{\prime }/\rho$ and keep composition, opacity coefficient and nuclear coefficient fixed first. Mass conservation gives $q=m/{\xi }^3$, and hydrostatic equilibrium gives $p=mq/\xi =m^2/{\xi }^4$. The ideal gas law then gives $t=p/q=m/\xi$. Scaling the stellar energy-generation rate equation and the stellar radiative temperature gradient gives respectivelyThe second follows from $dT/dr\propto -\kappa \rho L_r/(r^2{T}^3)$, not from energy production. Equating the two luminosity scalings givesThis assumes the denominator is nonzero and a consistent homologous family exists. If the mean molecular weight differs by $u={\mu }^{\prime }/\mu$, then $t=um/\xi$ andRatios of the opacity and energy-generation coefficients multiply the right-hand side. Thus fixed composition is a real restriction, not an automatic property of every stellar sequence.

Use the Stefan–Boltzmann law $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ to locate the family on a Hertzsprung-Russell diagram. It implies $T_{\mathrm{eff}}\propto M^{(Y-2X)/4}$ andFor proton–proton chain burning take the usual local approximation $\eta =4$, with Kramers' opacity law $\lambda =1$, $\nu =-7/2$. ThenFor the CNO cycle with electron-scattering opacity, $\lambda =\nu =0$, soA conventional local choice $\eta =16$ gives $R\propto M^{15/19}$ and slope $76/9\simeq 8.44$. Choosing $\eta =17$ instead gives slope $60/7\simeq 8.57$. The source specifies no numerical nuclear exponents, and the effective exponent changes with temperature; the general expression is the unambiguous answer.

The Hertzsprung-Russell diagram places hotter stars to the left. Its branches rise toward higher luminosity and mass, and the CNO cycle/electron-scattering opacity branch has the steeper logarithmic slope. Their illustrative joining point and normalization are arbitrary because the proportional opacity and reaction laws do not specify absolute stellar scales. These are the fully radiative ideal-gas homology predictions, not exact observed main sequence relations: convection and increasing radiation support limit those assumptions in real stars.

The mass-radius curve of solar-composition substellar objects reflects the transition from weak compression to pressure ionization and electron degeneracy pressure, followed by sustained hydrogen burning. A schematic joining representative object classes is:

The illustration is not an age-specific numerical evolutionary model. In particular, ice giants contain much more heavy material than a solar-composition giant, so a single uniform-composition equation of state does not describe the entire joined curve.

Planet/brown-dwarf naming conventions and deuterium burning do not define a universal discontinuity in the mass-radius relation. Composition, age and irradiation move the curves; the hydrogen-burning transition changes the long-term energy source more fundamentally.

For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate $m$. Let $u$ be specific thermal internal energy and $D/Dt$ follow a mass element. The stellar energy balance equation isThe outward luminosity is energy transported through the mass shell, ${ϵ}_{\mathrm{nuc}}$ is nuclear rest-mass energy released per unit mass and time, and ${ϵ}_{\nu }$ is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by ${ϵ}_{\mathrm{grav}}=-Du/Dt-PD(1/\rho )/Dt$. At fixed composition it is $-TDS/Dt$ by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in ${ϵ}_{\nu }$.

In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives $U\simeq -\Omega /2$, so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is $G{M}^2/(RL)$. After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.

Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is $Q=(m_{\mathrm{initial}}-m_{\mathrm{final}})c^2$, including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:For nonresonant charged-particle reactions the energy weighting contains $\exp [-E/(kT)-\sqrt{E_G/E}]$; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as $T^4$ or $T^{16}$ are local approximations.

For the proton–proton chain, first produce deuterium and helium-3:The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction $p+e^{-}+p\to {}^2\mathrm{H}+{\nu }_e$ also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:In pp III, the beryllium-7 instead captures a proton:Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately $26.7\mathrm{MeV}$ per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, ${ϵ}_{pp}\propto \rho X_H^2{T}^4$ is a useful local approximation.

The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists ofThe carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow ${}^{14}\mathrm{N}(p,\gamma ){}^{15}\mathrm{O}$ reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through ${}^{15}\mathrm{N}(p,\gamma ){}^{16}\mathrm{O}$, ${}^{16}\mathrm{O}(p,\gamma ){}^{17}\mathrm{F}$, ${}^{17}\mathrm{F}\to {}^{17}\mathrm{O}+e^{+}+{\nu }_e$, and ${}^{17}\mathrm{O}(p,\alpha ){}^{14}\mathrm{N}$, after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is ${ϵ}_{\mathrm{CNO}}\propto \rho X_H{X}_{\mathrm{CNO}}{T}^{\eta }$ with a much steeper exponent than the PP chain; the preceding homology problem specifies $\eta =16$. The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.

Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}$, releasing about $7.27\mathrm{MeV}$. Because three alpha particles are required, the specific rate scales as ${\rho }^2{Y}^3$; a schematic resonant rate isIts local temperature exponent is $-3+44/T_8$, of order forty near $T_8=1$. A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$ reaction determines much of the eventual carbon-oxygen mixture.

If the core becomes hot enough, carbon burning follows, for example through ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{20}\mathrm{Ne}+\alpha$ and ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{23}\mathrm{Na}+p$. Neon burning begins with photodisintegration, ${}^{20}\mathrm{Ne}+\gamma \to {}^{16}\mathrm{O}+\alpha$, followed by alpha capture such as ${}^{20}\mathrm{Ne}+\alpha \to {}^{24}\mathrm{Mg}+\gamma$. The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes ${}^{16}\mathrm{O}+{}^{16}\mathrm{O}\to {}^{28}\mathrm{Si}+\alpha$ and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.

The increasing nuclear binding energy per nucleon supplies progressively less energy as products approach the iron group. Further fusion past that region cannot provide sustained net heat to support a core. Massive stars can therefore end with an unstable iron-group core; lower-mass stars do not attain all these burning stages and instead leave remnants such as carbon-oxygen white dwarfs. Heavy-element neutron captures can synthesize nuclei beyond the iron group without being the principal hydrostatic power source.

Finally, stellar neutrino energy loss competes with every heat source. Reaction neutrinos accompany the proton–proton and CNO chains. Hot or dense material also emits neutrino pairs through $e^{-}+e^{+}\to \nu +\bar{\nu }$, plasmon decay, $e^{-}+\gamma \to e^{-}+\nu +\bar{\nu }$, and electron-ion bremsstrahlung. These neutrinos usually escape, unlike the photons whose transport is diffusive; neutrino trapping requires much more extreme collapse conditions. Late burning has smaller fuel-energy reservoirs and strong neutrino cooling, so its duration is much shorter than core hydrogen burning. The surface luminosity is the escaping photon power, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$, not the sum of photons and unobserved neutrino power. Radiative diffusion in a star, convection and electron conduction redistribute energy, while the energy-balance equation distinguishes true sources, losses and storage. Accretion or mass loss adds boundary energy and mechanical work when present, rather than changing the nuclear reaction budget.

Assume stellar homology: fixed dimensionless profiles, uniform fixed composition, negligible radiation pressure, ideal gas support, and radiative transport throughout the model. Mass conservation gives ${\rho }_c\propto M/R^3$. The hydrostatic pressure support equation gives $P_c\propto G{M}^2/R^4$; combining with $P\propto \rho T$ gives $T_c\propto M/R$.

Integrating specific proton–proton chain energy generation over mass givesFor radiative diffusion in a star, $T_c/R\propto {\kappa }_c{\rho }_{c}L/(R^2{T}_c^3)$, hence $L_{\mathrm{rad}}\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. With the Kramers' opacity law ${\kappa }_c\propto {\rho }_c{T}_c^{-7/2}$,Equating generated and transported luminosity gives $M^{1/2}\propto R^{13/2}$, so the radiative homology with proton-proton burning and Kramers opacity relation isThe proportionality constants are fixed only within the adopted homologous family, not for arbitrary evolving stars.

For a family in stellar homology, dimensionless radial profiles are fixed. Mass conservation and hydrostatic equilibrium then give the central scalingsThe ideal gas equation of state, with fixed mean molecular weight, consequently giveswhere fixed dimensional constants such as $G$ are suppressed in homology relations.

The proton–proton chain law $ϵ\propto \rho T^{7/2}$ gives the nuclear luminosityOn the other hand, radiative diffusion in a star givesUsing the Kramers' opacity law ${\kappa }_c\propto {\rho }_c{T}_c^{-7/2}$,Thermal equilibrium requires $L_{\mathrm{nuc}}=L_{\mathrm{rad}}$. Their common mass factor cancels, leaving $R^{-13/2}\propto R^{-1/2}$, and hence

This model captures the gas-pressure support, pp-chain burning, and strongly mass-dependent luminosity of the lower main sequence, including the Sun approximately. Its fully radiative assumption is an idealization: the Sun has a convective envelope, and sufficiently low-mass red dwarfs become largely or fully convective, so real radii are not exactly constant.