Proton-Proton Chain

Every joule of sunlight that has ever warmed the Earth was minted in the core of the Sun by a single, astonishingly slow process: four hydrogen nuclei welded into one helium-4 nucleus, with the mass that goes missing reappearing as energy. In stars of roughly a solar mass and below, that welding runs almost entirely through the proton-proton chain. The net reaction is $4\,{}^1\mathrm{H}\to{}^4\mathrm{He}$ , and because the helium-4 nucleus is bound more tightly than four loose protons, about $0.7\%$ of the reacting mass — some $26.73\,\mathrm{MeV}$ per helium nucleus — is released to hold the star up against gravity. The chain is built from a handful of two-body reactions, but its personality comes from its very first step. Fusing two protons requires not only tunnelling through the Coulomb barrier that like charges throw up, but also converting one proton into a neutron in the fleeting instant of contact — a weak-interaction event of extraordinary rarity. A given proton in the solar core waits on the order of ten billion years to take that step, and that single bottleneck sets the pace of the whole star: it is why the Sun burns slowly and lives for ten billion years rather than blazing out in a geological instant. Downstream of the bottleneck the chain is fast — deuterium survives seconds, helium-3 accumulates and closes the loop through one of three branches (PP I, PP II, PP III). Fold the thermal velocity distribution against the tunnelling probability and the energy-generation rate emerges as a gentle $\varepsilon_{pp}\propto T^{4}$ near solar conditions: a stable stellar thermostat, and the physics that sorts the lower main sequence from its hotter, CNO-powered siblings.

Key takeaways

The proton-proton chain is the dominant energy source in stars with central temperatures below about $1.8\times10^{7}\,\mathrm{K}$ , including the Sun, fusing four protons into one $^4\mathrm{He}$ nucleus.
The rate-limiting first step, $p+p\to{}^2\mathrm{H}+e^{+}+\nu_e$ , proceeds only through the weak interaction, making it fantastically slow and setting the long main-sequence lifetime.
Each completed chain liberates $Q\approx26.73\,\mathrm{MeV}$ from a mass deficit of $0.7\%$ ; roughly $0.5\,\mathrm{MeV}$ of that escapes as neutrinos in the PP I branch.
Three branches (PP I, PP II, PP III) share the deuterium and $^3\mathrm{He}$ intermediates but differ in how $^3\mathrm{He}$ is consumed; their relative weights shift with temperature.
Because the charges are small, the Coulomb barrier is low and the temperature sensitivity is mild, $\varepsilon_{pp}\propto T^{4}$ near solar conditions — a gently sloped, self-stabilizing rate law.
The neutrinos emitted along the chain escape the Sun directly and, once flavour oscillation was understood, confirmed the standard solar model and resolved the solar neutrino problem.
Key idea: Four protons become one helium-4 nucleus in a sequence gated by the weak-interaction fusion of two protons; the resulting energy-generation rate scales only as roughly $T^{4}$ , making the pp chain the stable, dominant power source of low-mass main-sequence stars.

Where this lives today

The pp chain is an active frontier of precision physics. In 2014 the Borexino experiment in Italy achieved the first direct, real-time detection of the low-energy pp neutrinos — from step (3), the very reaction that gates the Sun — measuring a flux of about $6\times10^{10}\,\mathrm{cm^{-2}\,s^{-1}}$ in agreement with the standard solar model to a few percent. Because the pp, $^7\mathrm{Be}$ , and $^8\mathrm{B}$ neutrinos are emitted in fixed, temperature-dependent ratios, their fluxes act as a thermometer for the solar interior, confirming a central temperature near $1.5\times10^{7}\,\mathrm{K}$ ; in 2020 Borexino went further and caught the subdominant CNO neutrinos for the first time. The same nuclear inputs — especially the S-factors of the $p+p$ and $^3\mathrm{He}+^4\mathrm{He}$ reactions, refined at underground accelerators like LUNA — feed directly into Big Bang nucleosynthesis, the ages of the oldest stars, and stellar-evolution models across astrophysics.

The physical picture

Fusing two protons is like trying to glue two beads that fiercely repel. First they must tunnel through the Coulomb barrier — already a rare event at the few-keV energies of a stellar core. But that is not enough: the diproton they briefly form is unbound, and it survives only if one proton weak-decays into a neutron ( $p\to n+e^{+}+\nu_e$ ) so the pair binds as deuterium. Tunnelling is rare; catching a weak decay inside the tunnelling window is astronomically rarer. This is the slow gate, and it throttles the luminosity of the entire star. Everything past the gate is fast. Deuterium, once made, is devoured by an ambient proton in about a second to build $^3\mathrm{He}$ , which accumulates until it finds a partner and the chain closes. The first figure lays out this architecture — a slow trunk feeding a rapid three-way cascade — with the real numbers on each link: timescales, solar branch fractions, and total energy released.

Figure

The proton-proton chain as a slow gate feeding a fast three-way cascade — two protons fuse by the weak interaction once every ~10¹⁰ yr per proton, then deuterium and helium-3 form in seconds and the chain closes through PP I, PP II or PP III with solar branchings near 83%, 17% and 0.02%

The gentle temperature dependence is the second essential fact, and it makes the Sun a thermostat rather than a bomb. Nudge the core hotter and $\varepsilon_{pp}\propto T^{4}$ climbs only modestly; the extra pressure lets the core expand and cool back down. Contrast the CNO cycle, whose $\varepsilon\propto T^{\sim17}$ concentrates burning into a tiny convective core, or degenerate helium ignition, which runs away explosively. The second figure plots the pp and CNO rates against temperature and marks the crossover near $1.8\times10^{7}\,\mathrm{K}$ , with the Sun on the pp side.

Figure

Energy-generation rate versus core temperature on log axes — the pp rate rises as T⁴ while the CNO rate rises as T¹⁷, crossing near 1.8×10⁷ K; the Sun at 1.5×10⁷ K sits on the pp side, hotter massive stars on the CNO side

The core equations

The chain is summarized by one net reaction. Four protons disappear; one helium-4 nucleus, two positrons, two electron neutrinos and a spray of gamma-ray photons appear. The positrons promptly annihilate on ambient electrons, adding their rest energy to the gas.

4\,{}^1\mathrm{H}\;\longrightarrow\;{}^4\mathrm{He}+2e^{+}+2\nu_e+\gamma \qquad (1)

where: ${}^1\mathrm{H}$ — proton; ${}^4\mathrm{He}$ — helium-4 (alpha particle); $e^{+}$ — positron, which annihilates with an electron; $\nu_e$ — electron neutrino, which escapes the star; $\gamma$ — gamma-ray photons carrying the released energy. Reading: two of the four protons must become neutrons, and each such conversion emits a neutrino, so the net chain sheds exactly two. The energy released is set entirely by the reactant-product mass difference through $E=mc^2$ . Using atomic (rather than bare-nuclear) masses is a bookkeeping trick that automatically folds in the annihilation of the two positrons with two atomic electrons.

Q=\left(4m_{{}^1\mathrm{H}}-m_{{}^4\mathrm{He}}\right)c^2=26.73\ \mathrm{MeV} \qquad (2)

where: $Q$ — total nuclear energy per chain, including positron annihilation; $m_{{}^1\mathrm{H}}=1.007825\,\mathrm{u}$ — neutral hydrogen atom; $m_{{}^4\mathrm{He}}=4.002602\,\mathrm{u}$ — neutral helium-4 atom; $c$ — speed of light. Reading: the fractional deficit is $0.028698/4.0313\approx0.71\%$ , the "0.7% efficiency" of hydrogen burning. The chain's first link is the one that matters. Two protons fuse to deuterium, and because a proton must become a neutron, the reaction is weak-mediated and by far the slowest step.

{}^1\mathrm{H}+{}^1\mathrm{H}\to{}^2\mathrm{H}+e^{+}+\nu_e\qquad(Q=1.442\ \mathrm{MeV}) \qquad (3)

where: ${}^2\mathrm{H}$ — deuterium (bound proton plus neutron); $Q=1.442\,\mathrm{MeV}$ includes annihilation — $0.420\,\mathrm{MeV}$ before, plus $1.022\,\mathrm{MeV}$ from the two $m_ec^2$ . Reading: the neutrino carries up to $0.420\,\mathrm{MeV}$ ; averaged over its spectrum, about $0.267\,\mathrm{MeV}$ . Deuterium is then captured by a proton within seconds to make helium-3, an electromagnetic reaction essentially instantaneous on stellar timescales.

{}^2\mathrm{H}+{}^1\mathrm{H}\to{}^3\mathrm{He}+\gamma\qquad(Q=5.493\ \mathrm{MeV}) \qquad (4)

where: ${}^3\mathrm{He}$ — helium-3; $\gamma$ — photon carrying the released energy. Reading: because this step is $\sim10^{18}$ times faster than step (3), deuterium is destroyed as fast as it is born and its equilibrium abundance is minuscule. In the dominant PP I branch the chain closes when two helium-3 nuclei meet, returning two protons so the net cost was four protons for one helium-4.

{}^3\mathrm{He}+{}^3\mathrm{He}\to{}^4\mathrm{He}+2\,{}^1\mathrm{H}\qquad(Q=12.859\ \mathrm{MeV}) \qquad (5)

where: ${}^4\mathrm{He}$ — the final helium-4 product; $2\,{}^1\mathrm{H}$ — two protons handed back to the plasma. Reading: to make one $^4\mathrm{He}$ this way, steps (3) and (4) each run twice, and $2(1.442)+2(5.493)+12.859=26.73\,\mathrm{MeV}$ recovers the total $Q$ . For stellar-structure calculations the microphysics is compressed into an energy-generation rate per unit mass, which near solar conditions takes an approximate power-law form.

\varepsilon_{pp}\simeq\varepsilon_{0,pp}\,\rho\,X^{2}\,T_6^{4}\,f_{pp}\,\psi_{pp} \qquad (6)

where: $\varepsilon_{pp}$ — energy generated per gram per second $(\mathrm{erg\,g^{-1}\,s^{-1}})$ ; $\rho$ — mass density; $X$ — hydrogen mass fraction (the rate $\propto X^{2}$ because two protons react); $T_6=T/10^{6}\,\mathrm{K}$ ; $f_{pp}$ — electron-screening enhancement; $\psi_{pp}$ — PP II/PP III branching correction. Reading: the $\rho X^{2}$ marks a two-body proton reaction, and $T_6^{4}$ is the tell-tale gentle slope. That $T^{4}$ is not put in by hand — it follows from thermally averaging the tunnelling cross-section, which peaks at the Gamow energy $E_0$ .

\langle\sigma v\rangle \propto T^{-2/3}\exp\!\left[-\frac{3E_0}{k_BT}\right],\qquad E_0=\left(\frac{b\,k_BT}{2}\right)^{2/3} \qquad (7)

where: $\langle\sigma v\rangle$ — thermally averaged rate per pair; $E_0$ — Gamow-peak (most effective) energy; $b=\pi\sqrt{2\mu}\,Z_1Z_2e^{2}/\hbar$ — Coulomb-barrier parameter ( $\mu$ the reduced mass); $k_B$ — Boltzmann constant. Reading: reactions happen not at $k_BT\sim1.3\,\mathrm{keV}$ but in a narrow window around $E_0\approx5.9\,\mathrm{keV}$ , high on the Maxwell tail where tunnelling is feasible. The local slope of the rate law follows by differentiating (7) — the number a stellar modeller quotes.

\nu\equiv\frac{\partial\ln\varepsilon}{\partial\ln T}=\frac{\tau-2}{3},\qquad \tau=\frac{3E_0}{k_BT} \qquad (8)

where: $\nu$ — local power-law index of the rate; $\tau$ — Gamow exponent. Reading: at $T_6\approx15$ , $\tau\approx13.7$ gives $\nu\approx3.9\approx4$ ; the small charges ( $Z=1$ ) keep $\tau$ small, which is why the pp slope is so gentle.

Derivation ATotal energy release from the mass deficit

The goal is to pin down the most-quoted number of the pp chain — the $26.73\,\mathrm{MeV}$ released per helium-4 — from nothing but tabulated atomic masses, and to see why using atomic masses (which include electrons) rather than bare nuclear masses is exactly the right bookkeeping once the two positrons annihilate. We compute the mass deficit, convert it, then account for the electron rest energies.

Step 1. Form the mass deficit from atomic masses: $\Delta m=4m_{{}^1\mathrm{H}}-m_{{}^4\mathrm{He}}=4(1.007825)-4.002602=0.028698\,\mathrm{u}$ . Why: the reaction destroys four hydrogen atoms' worth of mass and builds one helium atom, so the released energy is set by this difference and nothing else.

Step 2. Convert with $1\,\mathrm{u}=931.494\,\mathrm{MeV}/c^{2}$ : $\Delta m\,c^{2}=0.028698\times931.494=26.73\,\mathrm{MeV}$ . Why: mass and energy are the same currency through $E=mc^{2}$ , and the atomic mass unit already carries the conversion factor, so the arithmetic is a single multiplication.

Step 3. Check the electron bookkeeping. Four neutral hydrogen atoms bring four electrons; one neutral helium atom carries two. The two "extra" electrons are exactly the ones that annihilate the two positrons emitted in the chain, releasing $2m_ec^{2}=1.022\,\mathrm{MeV}$ . Using atomic masses folds this in automatically, so no separate $2m_ec^{2}$ is added, and the deficit is $0.028698/4.0313\approx0.71\%$ of the reacting mass — the "0.7% efficiency" that sets stellar lifetimes. Why: this is the classic trap — with bare nuclear masses you would add the annihilation energy by hand, and forgetting or double-counting it mis-states $Q$ . *The results $\Delta m=0.028698\,\mathrm{u}$ and $Q=26.73\,\mathrm{MeV}$ are confirmed by the step-by-step sum $2(1.442)+2(5.493)+12.859=26.73\,\mathrm{MeV}$ .*

Go deeper · Derivation B · 1 minThe temperature exponent of the pp rate

Derivation BThe temperature exponent of the pp rate

The aim is to show that the famous $\varepsilon_{pp}\propto T^{4}$ is not an empirical fit but a direct consequence of the Gamow-peak structure of the rate, and to extract the general formula $\nu=(\tau-2)/3$ that reads the temperature sensitivity of any charged-particle reaction off its Gamow exponent. We write the rate as a power of temperature, take a logarithmic derivative, and evaluate at solar conditions.

Step 1. Write the rate from equation (7) as $\langle\sigma v\rangle\propto T^{-2/3}e^{-\tau}$ with $\tau=3E_0/k_BT$ , and note that since $E_0\propto T^{2/3}$ we have $\tau\propto T^{2/3}/T=T^{-1/3}$ . Why: both the algebraic prefactor and the exponential depend on temperature, and the exponential — through $\tau$ — is by far the more sensitive, so we must track how $\tau$ itself varies with $T$ .

Step 2. Differentiate logarithmically: $\dfrac{\mathrm{d}\ln\langle\sigma v\rangle}{\mathrm{d}\ln T}=-\dfrac{2}{3}-\dfrac{\mathrm{d}\tau}{\mathrm{d}\ln T}$ , and with $\tau\propto T^{-1/3}$ we have $\mathrm{d}\tau/\mathrm{d}\ln T=-\tau/3$ . Why: the logarithmic derivative is exactly the local power-law index we want, and the chain rule turns the messy exponential into a clean algebraic term in $\tau$ .

Step 3. Combine to get $\nu=\dfrac{\tau-2}{3}$ . Evaluate at the solar core: with $E_0\approx5.9\,\mathrm{keV}$ and $k_BT\approx1.29\,\mathrm{keV}$ at $T_6=15$ , $\tau=3(5.9)/1.29\approx13.7$ , so $\nu=(13.7-2)/3\approx3.9$ . At lower $T$ the exponent $\tau$ grows and $\nu$ steepens; at higher $T$ it flattens — so $\nu\approx4$ is a solar value, not a universal constant. Why: the small charges of the two protons keep $\tau$ small, giving $\nu\approx4$ , far gentler than the $\nu\approx17$ of CNO or $\nu\approx40$ of the triple-alpha process. *These values — $\tau\approx13.7$ and $\nu\approx3.9$ at $T_6=15$ — follow directly from $E_0=1.22(Z_1^{2}Z_2^{2}A\,T_6^{2})^{1/3}\,\mathrm{keV}$ .*

Go deeper · Derivation C · 1 minNeutrino losses and the heat that reaches the gas

Derivation CNeutrino losses and the heat that reaches the gas

The purpose is to separate the total nuclear yield of $26.73\,\mathrm{MeV}$ from the energy that actually heats the star, since the two electron neutrinos emitted per chain stream out of the Sun without interacting and carry their energy away instantly. We count how many times each reaction runs to build one helium-4 through PP I, tally the neutrino energy, and subtract it to find the deposited heat.

Step 1. Count the reactions in a PP I closure: to assemble one $^4\mathrm{He}$ , the weak step (3) runs twice, the capture (4) runs twice, and the closing step (5) runs once. Why: only step (3) makes neutrinos, so the neutrino budget is fixed by how many times the first step fires — here, twice.

Step 2. Tally the neutrino energy. Each pp neutrino carries an average $\langle E_\nu\rangle\approx0.267\,\mathrm{MeV}$ (spectrum endpoint $0.420\,\mathrm{MeV}$ ), so two remove $E_{\nu,\mathrm{tot}}=2\times0.267\approx0.53\,\mathrm{MeV}$ . Why: interacting only weakly, this energy leaves the star at essentially light speed and is lost to the luminosity budget.

Step 3. Subtract to find the deposited heat: $E_{\mathrm{thermal}}=Q-2\langle E_\nu\rangle\approx26.73-0.53\approx26.2\,\mathrm{MeV}$ . Thus about $2\%$ of the nuclear energy escapes immediately as neutrinos in PP I; the other branches lose more because their neutrinos are more energetic. Why: it is this $\sim26.2\,\mathrm{MeV}$ , not the full $26.73$ , that emerges as sunlight after a photon-diffusion delay of $\sim10^{5}$ yr; using the wrong figure biases luminosity estimates by a couple of percent.

E_{\mathrm{thermal}}=Q-2\langle E_\nu\rangle\approx26.73-2(0.267)\approx26.20\ \mathrm{MeV}

The deposited thermal energy $E_{\mathrm{thermal}}\approx26.2\,\mathrm{MeV}$ is consistent with standard-solar-model energy accounting.

Limiting cases

The behaviour and dominance of the pp chain hinge on temperature, composition, and the competition with other channels.

Case	Result	Meaning
Low core temperature $(T\lesssim1.8\times10^{7}\,\mathrm{K})$	$\varepsilon_{pp}\gg\varepsilon_{\mathrm{CNO}}$	pp chain dominates — the Sun and the lower main sequence
High core temperature $(T\gtrsim1.8\times10^{7}\,\mathrm{K})$	$\varepsilon_{\mathrm{CNO}}\gg\varepsilon_{pp}$	the steep CNO cycle takes over, in stars above roughly $1.2\,M_\odot$
Branch versus temperature	PP I for $T_6\lesssim14$ ; PP II for $14\lesssim T_6\lesssim23$ ; PP III for $T_6\gtrsim23$	higher $T$ shifts $^3\mathrm{He}$ destruction toward the $^7\mathrm{Be}$ and $^8\mathrm{B}$ channels
Effective power-law index	$\nu=(\tau-2)/3\approx4$ at $T_6\approx15$	slope steepens at lower $T$ , flattens at higher $T$ ; $\nu\approx4$ is the solar value

Common trapThe 26.73 MeV is not all available to heat the star

The total $Q=26.73\,\mathrm{MeV}$ already includes the $2m_ec^{2}=1.022\,\mathrm{MeV}$ from positron annihilation — so do not add annihilation energy a second time. But it is equally wrong to treat all $26.73\,\mathrm{MeV}$ as heat: the two neutrinos carry off about $0.53\,\mathrm{MeV}$ in PP I (more in PP II and PP III, where the $^7\mathrm{Be}$ line and the energetic $^8\mathrm{B}$ spectrum steal a larger share). Neutrinos leave without depositing their energy, so the heat available as sunlight is about $26.2\,\mathrm{MeV}$ per helium nucleus in PP I. Keep three quantities distinct: the total nuclear $Q$ , the neutrino loss, and the deposited thermal energy.

Worked examples

1 · How many pp chains power the Sun (fully worked)

Given $L_\odot=3.83\times10^{33}\,\mathrm{erg\,s^{-1}}$ and $26.2\,\mathrm{MeV}$ deposited per chain, find the chains per second, the proton consumption rate, and the pp-neutrino flux at Earth ( $d=1\,\mathrm{AU}$ ).

Step 1. Convert to cgs: $26.2\,\mathrm{MeV}=26.2\times1.602\times10^{-6}\,\mathrm{erg}=4.20\times10^{-5}\,\mathrm{erg}$ , using $1\,\mathrm{MeV}=1.602\times10^{-6}\,\mathrm{erg}$ .

Step 2. Divide luminosity by energy per chain: $R_{\mathrm{chains}}=L_\odot/E_{\mathrm{thermal}}=(3.83\times10^{33})/(4.20\times10^{-5})\approx\mathbf{9.1\times10^{37}\ \mathrm{s^{-1}}}$ — also the rate of helium-4 production in the core.

Step 3. Each chain consumes four protons, so $\dot N_p=4R_{\mathrm{chains}}\approx\mathbf{3.6\times10^{38}\ \mathrm{protons\,s^{-1}}}$ — about $6\times10^{11}\,\mathrm{kg\,s^{-1}}$ of hydrogen fused, of which $\sim4\times10^{9}\,\mathrm{kg\,s^{-1}}$ becomes radiated energy.

Step 4. Two neutrinos leave per chain over a sphere of radius $d=1.496\times10^{13}\,\mathrm{cm}$ : $\Phi_\nu=2R_{\mathrm{chains}}/(4\pi d^{2})\approx\mathbf{6\times10^{10}\ \mathrm{cm^{-2}\,s^{-1}}}$ , matching the measured solar pp-neutrino flux — number for number, the Sun shines by hydrogen fusion.

2 · The Gamow peak for the pp reaction (fully worked)

At the solar centre $(T_6=15)$ , find the Gamow-peak energy $E_0$ , compare it to the thermal energy $k_BT$ , and confirm the temperature exponent $\nu\approx4$ . Use $E_0=1.22\,(Z_1^{2}Z_2^{2}A\,T_6^{2})^{1/3}\ \mathrm{keV}$ with reduced mass $A=0.5$ for two protons.

Step 1. Evaluate the Gamow energy: $E_0=1.22\,(1\cdot1\cdot0.5\cdot15^{2})^{1/3}=1.22\,(112.5)^{1/3}=1.22\times4.83\approx\mathbf{5.9\ \mathrm{keV}}$ .

Step 2. Compare with the thermal scale: $k_BT=8.617\times10^{-5}\,\mathrm{eV\,K^{-1}}\times1.5\times10^{7}\,\mathrm{K}=1.29\,\mathrm{keV}$ . So reactions occur at $E_0/k_BT\approx4.6$ times the mean thermal energy — high on the Maxwell-Boltzmann tail, where only a rare fast pair can tunnel.

Step 3. Form the Gamow exponent: $\tau=3E_0/k_BT=3(5.9)/1.29\approx\mathbf{13.7}$ . The peak has a full width of about $4\sqrt{E_0k_BT/3}\approx6.4\,\mathrm{keV}$ , so fusion draws from a thin band of energies around $6\,\mathrm{keV}$ .

Step 4. Read off the temperature sensitivity: $\nu=(\tau-2)/3=(13.7-2)/3\approx\mathbf{3.9\approx4}$ . The gentle slope traces directly to the small $\tau$ , which traces to the small $Z_1Z_2=1$ — the physical root of the pp chain's stability.

3 · The main-sequence lifetime of the Sun

Estimate how long the Sun can shine on the pp chain, assuming it burns about $10\%$ of its mass at $0.7\%$ efficiency.

Step 1. Nuclear store: $E_{\mathrm{nuc}}\approx0.007\times0.1\times M_\odot c^{2}=0.007\times0.1\times(1.99\times10^{33}\,\mathrm{g})(3.0\times10^{10}\,\mathrm{cm\,s^{-1}})^{2}\approx1.25\times10^{51}\,\mathrm{erg}$ .

Step 2. Divide by luminosity: $t_{\mathrm{MS}}\approx E_{\mathrm{nuc}}/L_\odot\approx(1.25\times10^{51})/(3.83\times10^{33})\approx3.3\times10^{17}\,\mathrm{s}\approx\mathbf{1.0\times10^{10}\ \mathrm{yr}}$ . The ten-billion-year solar lifetime is a direct consequence of the slow weak gate: the tiny cross-section of step (3) stretches a modest energy store over a cosmological span.

Physics of the portrait

The pp chain is where two revolutions of twentieth-century physics converged: the quantum tunnelling that lets sub-barrier protons fuse, and the weak interaction that turns one into a neutron. No one did more to weld these into a working theory of how stars shine than Hans Bethe.

Origins · 1 minHans Bethe and how this was discovered

Hans Bethe

1906–2005 · nuclear astrophysics and stellar energy generation

Hans Bethe was the physicist who finally answered a question that had haunted astronomy since the nineteenth century: what keeps the Sun burning? Trained under Arnold Sommerfeld, Bethe fled Germany in 1933 and settled at Cornell, where in 1938 — spurred by a conference George Gamow organized on stellar energy — he set out to identify the reactions that power stars. With Charles Critchfield he worked out the rate of the proton-proton chain's first step, showing that the weak fusion of two protons into deuterium could supply the Sun's luminosity. In a companion paper that same year he laid out the carbon-nitrogen-oxygen cycle and recognized that it, not the pp chain, dominates in hotter, more massive stars — drawing for the first time the temperature boundary between the two burning modes that organizes the whole main sequence. The work earned him the 1967 Nobel Prize in Physics, "for his discoveries concerning the energy production in stars." He built on Arthur Eddington's conviction that subatomic energy must be the source, and on the Atkinson-Houtermans insight that Gamow's tunnelling made stellar fusion possible, supplying the quantitative machinery that turned those ideas into predictive stellar physics.

History

1920

— Eddington argues that subatomic energy, not gravitational contraction, powers the stars

1929

— Atkinson and Houtermans apply Gamow's quantum tunnelling to fusion in stellar interiors

1938

— Bethe and Critchfield derive the rate of the proton-proton chain

1958

— Cameron and others map the PP II and PP III branches through the ⁷Be and ⁸B channels

2002

— SNO confirms neutrino flavour oscillation, resolving the solar neutrino problem

Going further

Texts. Carroll and Ostlie, An Introduction to Modern Astrophysics, 2nd ed., ch. 10 (nuclear energy generation) · Clayton, Principles of Stellar Evolution and Nucleosynthesis, chs. 4–5 (reaction rates and the pp chain) · Kippenhahn, Weigert and Weiss, Stellar Structure and Evolution, 2nd ed., ch. 18 · Bahcall, Neutrino Astrophysics (1989) · Iliadis, Nuclear Physics of Stars, 2nd ed., ch. 5. The total $Q=26.73\,\mathrm{MeV}$ is computed from atomic masses $m_{{}^1\mathrm{H}}=1.007825\,\mathrm{u}$ , $m_{{}^4\mathrm{He}}=4.002602\,\mathrm{u}$ , and $1\,\mathrm{u}=931.494\,\mathrm{MeV}/c^{2}$ , and confirmed by the step-by-step sum $2(1.442)+2(5.493)+12.859=26.73\,\mathrm{MeV}$ . The average pp-neutrino energy is $\langle E_\nu\rangle\approx0.267\,\mathrm{MeV}$ . The solar burning rate uses $L_\odot=3.83\times10^{33}\,\mathrm{erg\,s^{-1}}$ ; the Gamow energy formula $E_0=1.22(Z_1^{2}Z_2^{2}A\,T_6^{2})^{1/3}\,\mathrm{keV}$ gives $E_0\approx5.9\,\mathrm{keV}$ , $\tau\approx13.7$ , $\nu\approx3.9$ at $T_6=15$ ; the pp/CNO crossover sits near $1.8\times10^{7}\,\mathrm{K}$ . Key anchor results: chain rate $\sim9\times10^{37}\,\mathrm{s^{-1}}$ , neutrino flux $\sim6\times10^{10}\,\mathrm{cm^{-2}\,s^{-1}}$ , and main-sequence lifetime $\sim10^{10}\,\mathrm{yr}$ .

Prerequisites. Nuclear Reaction Rates and the Gamow Peak · The Main Sequence · The Hertzsprung-Russell Diagram. Leads to. The CNO Cycle · Helium Burning and the Triple-Alpha Process · Post-Main-Sequence Evolution and the Red Giant Branch · The Mass-Luminosity Relation.

Ask the tutor: Ask me to walk through why a single weak-interaction reaction can throttle an entire star, or to work the Gamow-peak numbers for a hotter stellar core.

Why is fusing two protons so much harder than the later steps?

Two obstacles stack up. Both protons are positively charged, so they must tunnel through a Coulomb barrier — but the same is true of the later steps, so tunnelling alone is not the bottleneck. The decisive obstacle is that two protons cannot bind: the diproton is unstable. For deuterium to form, one proton must convert to a neutron during the brief instant of contact, and that is a weak-interaction beta decay. Weak processes are enormously slower than the strong and electromagnetic ones driving the rest of the chain, so a proton in the solar core waits on the order of ten billion years before it fuses — which is what makes stars long-lived.

If the reaction is so slow, how does the Sun produce so much power?

By sheer numbers. Any individual proton waits billions of years, but the Sun's core holds something like $10^{56}$ protons, so even a vanishingly small per-proton rate multiplies into roughly $9\times10^{37}$ completed chains every second — converting about $6\times10^{11}$ kilograms of hydrogen to helium each second and radiating $3.8\times10^{26}$ watts. The slowness is a feature: it is precisely because the per-proton rate is so tiny that the enormous energy store is metered out over ten billion years instead of released all at once. A faster first step would make a much shorter-lived, and far less hospitable, star.

What is the difference between PP I, PP II and PP III?

All three share the same opening — two protons make deuterium, which captures a proton to make $^3\mathrm{He}$ — and differ only in how $^3\mathrm{He}$ is turned into $^4\mathrm{He}$ . In PP I two $^3\mathrm{He}$ nuclei fuse directly. In PP II, $^3\mathrm{He}$ captures an existing $^4\mathrm{He}$ to make $^7\mathrm{Be}$ , which grabs an electron to become $^7\mathrm{Li}$ , which captures a proton and splits into two $^4\mathrm{He}$ . PP III is a rarer offshoot where $^7\mathrm{Be}$ captures a proton to make $^8\mathrm{B}$ , which beta-decays and breaks apart. Which route dominates depends on temperature: PP I in cooler cores like the Sun's (about $83\%$ ), PP II and PP III growing as the core heats up.

Where do the solar neutrinos come from, and why do they matter?

Every proton-to-neutron conversion emits an electron neutrino, and different steps produce characteristic energies: low-energy pp neutrinos from the first step, a monoenergetic line from $^7\mathrm{Be}$ capture, and high-energy neutrinos from $^8\mathrm{B}$ decay. Because they interact so weakly, neutrinos stream straight out of the core and reach us in eight minutes, giving a real-time view of the Sun's centre — while the photons we see took $\sim10^{5}$ years to random-walk out. Historically, detectors found only about a third of the expected high-energy neutrinos, the "solar neutrino problem." SNO confirmed the resolution in 2002: neutrinos oscillate between flavours in transit, and the solar model was right all along.

Why does the pp chain make a stable star but the triple-alpha process can run away?

It comes down to the temperature exponent $\nu=\partial\ln\varepsilon/\partial\ln T$ . For the pp chain $\nu\approx4$ , because the low nuclear charges give a small Gamow exponent, and that gentle slope makes the core a thermostat: a small temperature rise boosts generation only modestly, the gas expands and cools, and equilibrium is restored. The triple-alpha reaction, by contrast, involves higher charges and a resonance, giving $\nu\approx40$ ; when helium ignites in a degenerate core — where pressure does not respond to temperature — that ferocious sensitivity produces a runaway, the helium flash. Same Gamow-peak physics, wildly different stability, all encoded in $\nu$ .

How do we know the pp chain is really what powers the Sun?

Three independent lines converge. The energy budget: the measured chain rate times the deposited energy reproduces the solar luminosity, and the $0.7\%$ efficiency yields the ten-billion-year lifetime consistent with the ages of the oldest rocks and stars. The neutrinos: detectors from Homestake to Borexino have measured the pp, $^7\mathrm{Be}$ , and $^8\mathrm{B}$ fluxes, and Borexino's 2014 direct detection of the pp neutrinos confirmed the gating reaction in real time. And helioseismology: the sound-speed profile from solar oscillations matches a pp-powered model to fractions of a percent. No competing source survives all three.