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]

For a nonresonant charged-particle thermonuclear reaction rate, the thermal high-energy tail favors low $E$ while the Coulomb barrier penetration factor favors high $E$. Their product peaks in an intermediate energy window. Resonant reactions such as the Triple-alpha process require an additional resonance treatment.

For a slowly varying nuclear factor and penetration factor $\exp [-2\sqrt{E_B/E}]$, the thermal integral is ${\int }_0^{\infty }{e}^{-\Phi (E)}dE$, where $\Phi (E)=E/(k_{B}T)+2\sqrt{E_B/E}$. The Gamow peak occurs at $E_0=[E_B(k_{B}T{)}^2{]}^{1/3}$, with $\Phi (E_0)=3E_0/(k_{B}T)$ and ${\Phi }^{\prime \prime }(E_0)=3/(2k_{B}T{E}_0)$. The Gaussian integral gives $\sqrt{4\pi k_{B}T{E}_0/3}{e}^{-3E_0/(k_{B}T)}$. Combined with the Maxwellian prefactor $(k_{B}T{)}^{-3/2}$ this produces the displayed thermonuclear reaction rate at fixed reactant density. The approximation assumes $E_B\gg k_{B}T$ and excludes strong resonances or rapidly varying nuclear factors.

In the pp-I limit, helium-3 relaxation time is ${\tau }^{-1}=n_1\sqrt{2{\lambda }_{11}{\lambda }_{33}}$. At fixed proton number density, each thermonuclear reaction rate coefficient obeys $\lambda (T)/\lambda (T_0)=(T/T_0{)}^{-2/3}\exp [-\eta (T_0)((T_0/T{)}^{1/3}-1)]$. Taking their geometric mean gives the displayed relaxation ratio with $\bar{\eta }=[{\eta }_{11}(T_0)+{\eta }_{33}(T_0)]/2$. A local power-law approximation instead gives $\tau \propto T^{-(\alpha +\beta )/2}$. These are estimates rather than a complete radial model: the proton number density, composition and pp-II branch change along a real stellar composition profile.

The event rate counts unordered pairs of reacting nuclei. For unlike species there are $n_i{n}_j$ pairs per unit-volume normalization; for identical species there are approximately $n_i^2/2$ pairs, since counting ordered pairs counts each collision twice. The Kronecker delta therefore supplies exactly the required divisor $1+{\delta }_{ij}$. It is the event rate, not the destruction rate of identical nuclei: one identical-pair reaction destroys two nuclei.

The temperature dependence combines the high-energy tail of the Maxwell-Boltzmann distribution with quantum tunnelling through the Coulomb barrier. The principal exponential factor in the energy integral is $\exp [-E/(k_{B}T)-b/\sqrt{E}]$. Its stationary point, the Gamow peak, has $E_0\propto T^{2/3}$ and exponent $\eta =3E_0/(k_{B}T)\propto T^{-1/3}$. Integrating around this peak gives a prefactor proportional to $T^{-2/3}$, or ${\eta }^2$, when the nuclear factor varies slowly. Thus the supplied thermonuclear reaction rate is locally proportional to ${\eta }^2{e}^{-\eta }$. The slow weak interaction in the proton-proton reaction also makes its overall coefficient very small; its local thermal slope follows the supplied formula.

At fixed number densities, define the local temperature exponent of a thermonuclear reaction by a logarithmic derivative. Since $d\eta /d\log T=-\eta /3$,For proton-proton reactions, $A=1/2$ and $Z_1=1$. At $T_6=15$, ${\eta }_{11}=13.67$ and $\alpha =3.89$. For two helium-3 nuclei, $A=3/2$ and $Z_3^2{Z}_3^2=16$, giving ${\eta }_{33}=49.68$ and $\beta =15.89$. For helium-3 with helium-4, $A=12/7$ and the same charge product gives ${\eta }_{34}=51.95$ and $\gamma =16.65$. HenceThese are local powers near the specified temperature, not constant exponents over an arbitrarily large range. The greater reduced nuclear mass in the $34$ reaction increases its tunnelling exponent and its thermal sensitivity.

For the effective proton-proton reaction network, denote the event rates byThe short-lived beryllium intermediate can be eliminated, so production through $R_{34}$ feeds $n_7$, representing predominantly lithium-7. This rapid conversion is electron capture. Count the nuclei consumed and produced in each reaction: $11$ consumes two protons and makes one deuterium; $21$ consumes a proton and a deuterium and makes one helium-3; $33$ consumes two helium-3 and makes one helium-4 and two protons; $34$ consumes one of each helium isotope; $17$ consumes a proton and lithium-7 and makes two helium-4. It follows thatFor closure, ${\dot{n}}_7={\lambda }_{34}{n}_3{n}_4-{\lambda }_{17}{n}_1{n}_7$. These nuclear terms conserve $n_1+2n_2+3n_3+4n_4+7n_7$, providing an independent stoichiometric check. They describe a fixed local volume; compression or mixing would add transport terms.

The roughly one-second destruction time of deuterium is negligible compared with the proton-burning and helium-3 evolution times. Its abundance rapidly adjusts to the slowly evolving production rate. Set ${\dot{n}}_2\simeq 0$, giving the deuterium quasi-equilibrium in proton-proton burningSubstitution, without setting ${\dot{n}}_1=0$, givesThe helium-3 relaxation time of about $6\times 10^5$ years is likewise short compared with the central proton-burning time and the age of the Sun. On this relaxation time the proton and helium-4 abundances, and the background temperature, can be treated as fixed. Put $A={\lambda }_{33}$, $B={\lambda }_{34}{n}_4$, $C={\lambda }_{11}{n}_1^2/2$. Then ${\dot{n}}_3=C-A{n}_3^2-B{n}_3$, and its nonnegative equilibrium isThe other root is negative and unphysical. Write $n_3=n_{3e}+x$ and subtract the equilibrium equation. Exactly, for a fixed background,Dropping the quadratic displacement givesThus small displacements decay exponentially. A slowly moving equilibrium introduces an additional $-{\dot{n}}_{3e}$ term, negligible to leading order when $\tau$ is short compared with background evolution.

To estimate the temperature for helium-3 freeze-out, use the cooler pp-I-dominated region and neglect $B$ compared with $2A{n}_{3e}$. Then $n_{3e}\simeq n_1\sqrt{{\lambda }_{11}/(2{\lambda }_{33})}$ and ${\tau }^{-1}\simeq n_1\sqrt{2{\lambda }_{11}{\lambda }_{33}}$. With the local powers just obtained and fixed proton number density, $n_{3e}\propto T^{-6}$ and $\tau \propto T^{-10}$. Adopting $t_{\odot }\simeq 4.6\times 10^9\mathrm{yr}$ givesThis is an order-of-magnitude extrapolation of local exponents. Keeping the supplied exponential factors gives the helium-3 freeze-out estimate with Gamow factorswhich yields about $6.8\times 10^6\mathrm{K}$ with the same normalization and fixed-density approximation. Thus the robust estimate is several million kelvin. The supplied data do not fix the local density, composition or pp-II fraction away from the centre, so they cannot determine an exact transition temperature or radius. In particular, using $\tau \propto T^{-16}$ alone would omit the temperature dependence of the equilibrium helium-3 abundance.

The solar hydrogen and helium-3 abundance profiles follow from burning, finite equilibration time and convection. Central hydrogen has been consumed, so $X_1$ is lowest in the centre and increases toward its much less processed outer value. In the hot centre, rapid destruction keeps helium-3 small. Moving outward initially raises its equilibrium abundance because ${\lambda }_{11}/{\lambda }_{33}$ rises strongly as temperature falls. Further out, equilibration is no longer reached within the solar age, and still farther out even production is slow. Hence $X_3$ has an off-centre maximum and then declines toward its small outer abundance. The outer convective envelope mixes each abundance toward a constant value.

The sketch shows these qualitative shapes; the helium-3 ordinate is scaled to its own maximum. The illustrated amplitudes, peak radius and envelope boundary are schematic, not a numerical solar-model solution.

Thermal energy versus electrostatic repulsion. The mean translational kinetic energy of one nonrelativistic thermal proton isUsing the center-to-center touching separation $2r_p$, the Coulomb barrier estimate isThus $E_C/⟨E⟩\simeq 274$. Using $r_p$ rather than $2r_p$ as the rough interaction distance doubles the estimate but leaves the conclusion unchanged. Typical solar thermal energies are far below the classical barrier. The two essential ideas are the high-energy tail of the thermal distribution and quantum tunnelling through the Coulomb barrier: neither a mean-energy estimate nor a sharp classical cutoff describes the reacting population.

For two equal-mass protons, let $v$ be their relative speed. In the center-of-mass frame each has speed $v/2$, so their total collision energy isHere ${\mu }_r$ is the reduced mass, not the stellar mean molecular weight.

In the low-energy, nonresonant approximation, the geometrical quantum-wave factor in a collision cross-section is proportional to the squared de Broglie wavelength, hence to $1/E$. The slowly varying short-range reaction probability is absorbed into $S_0$; for proton–proton fusion it also includes the small weak-reaction probability. The rapidly varying exponential is the quantum tunnelling probability. For a Coulomb potential $V(r)=e^2/r$, its leading semiclassical exponent isThe integral follows from $r=(e^2/E){\sin }^2\theta$; its dimensionless integral is $\pi /2$. A finite nuclear radius modifies the short-distance factor absorbed into $S_0$. In this convention $E_B={\pi }^2{\mu }_r{e}^4/(2{\hslash }^2)$: it is a penetration-energy scale, distinct from the classical contact energy $E_C$.

Thermal averaging and pair counting. Each proton samples collision partners with number density $N$ and collision flux $v$, so its collision rate is $N⟨\sigma v⟩$. Multiplying by the number of protons counts each identical pair twice. The thermonuclear reaction rate is thereforeThe relative speed of two independent thermal protons has a Maxwell-Boltzmann velocity distribution with reduced mass $m_p/2$. Changing variables with $E=m_p{v}^2/4$ gives the normalized energy densitySubstituting $v=2\sqrt{E/m_p}$ and the collision cross-section causes the energy powers to cancel:

The Gamow peak and its Gaussian width. Set $\Phi (E)=E/(kT)+2\sqrt{E_B/E}$. It diverges at either endpoint and has its unique minimum whereThus the energy called $E_G$ in this question is the Gamow peak energy,In the stellar tunnelling regime $E_B\gg kT$,This establishes the requested hierarchy in its physical regime; it is not an inequality valid at arbitrary temperature.

At the minimum, $\Phi (E_G)=3E_G/(kT)$ and ${\Phi }^{\prime \prime }(E_G)=3/(2kT{E}_G)$. HenceThe Gaussian standard deviation is $\sqrt{2kT{E}_G/3}$, small compared with $E_G$ when $kT\ll E_G$. Extending the lower endpoint to $-\infty$ relative to the peak then givesCombining this with the thermal prefactor yields the Gaussian approximation to a nonresonant thermonuclear rate,At fixed reactant density and fixed $S_0,E_B$, comparison with the requested temperature form givesThe constant $\beta$ here has dimensions of temperature; it is not inverse temperature. Resonances, plasma screening or variation of $S(E)$ would require corrections to this particular approximation.

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.

Collective electric charge rearrangement modifies particle interactions in a plasma. In stellar nuclear fusion, this response can alter the Coulomb barrier used to infer a thermonuclear reaction rate.