2.4 Stellar Nucleosynthesis Processes

Hydrogen Burning Processes

Proton-Proton Chain Reaction

The pp-chain is the dominant hydrogen-burning mechanism in low-mass stars (roughly 0.08 to 1.5 $M_\odot$). It proceeds through three branches, with PP-I, PP-II, and PP-III contributing differently depending on core temperature.

PP-I branch dominates in the Sun and cooler stars. The sequence is:

1. Two protons fuse to form deuterium ($^2H$), releasing a positron and an electron neutrino. This weak-interaction step is the bottleneck of the entire chain.
2. The deuterium nucleus captures another proton to form $^3He$, releasing a gamma ray.
3. Two $^3He$ nuclei combine to produce $^4He$ and two protons.

The net result: four protons become one $^4He$ nucleus, two positrons, two electron neutrinos, and 26.73 MeV of energy.

PP-II and PP-III branches become more significant at higher core temperatures (above ~17 million K). These branches route through $^7Be$ and $^7Li$ as intermediates. PP-III additionally involves $^8B$ decay and produces the high-energy neutrinos that solar neutrino experiments like Super-Kamiokande are most sensitive to.

Energy release comes primarily as gamma rays and kinetic energy of the produced particles, with a small fraction carried away by neutrinos (which escape the star entirely).

CNO Cycle

The CNO cycle dominates hydrogen burning in stars more massive than about 1.3 $M_\odot$. Rather than fusing protons directly, it uses pre-existing $^{12}C$, $^{14}N$, and $^{16}O$ nuclei as catalysts. The catalyst nuclei are recycled at the end of each cycle, so the net reaction is the same as the pp-chain: four protons yield one $^4He$.

CNO-I cycle (the dominant branch):

1. $^{12}C$ captures a proton → $^{13}N$ + $\gamma$
2. $^{13}N$ undergoes $\beta^+$ decay → $^{13}C$ + $e^+$ + $\nu_e$
3. $^{13}C$ captures a proton → $^{14}N$ + $\gamma$
4. $^{14}N$ captures a proton → $^{15}O$ + $\gamma$ (this is the slowest step, so $^{14}N$ accumulates)
5. $^{15}O$ undergoes $\beta^+$ decay → $^{15}N$ + $e^+$ + $\nu_e$
6. $^{15}N$ captures a proton → $^{12}C$ + $^4He$

CNO-II becomes relevant at still higher temperatures and routes through $^{16}O$, $^{17}F$, and $^{17}O$ as intermediates.

A critical distinction: the CNO cycle's energy generation rate scales as roughly $\epsilon \propto T^{16-20}$, compared to $\epsilon \propto T^4$ for the pp-chain. This extreme temperature sensitivity is why the CNO cycle dominates in more massive (hotter-core) stars and why those stars develop convective cores. The CNO cycle also produces a distinctly different neutrino energy spectrum, which has observational consequences for neutrino astronomy.

Advanced Burning Stages

Each successive burning stage requires higher temperatures to overcome increasing Coulomb barriers, runs on shorter timescales, and produces less energy per unit mass. For a ~25 $M_\odot$ star, hydrogen burning lasts millions of years while silicon burning lasts roughly a day.

Helium Burning

Helium burning ignites when the core reaches ~$10^8$ K and hydrogen is exhausted in the core. The primary mechanism is the triple-alpha process:

1. Two $^4He$ nuclei fuse to form $^8Be$, which is extremely unstable (lifetime ~$10^{-16}$ s).
2. Before the $^8Be$ decays, a third $^4He$ nucleus is captured to form an excited state of $^{12}C$.
3. This excited state is the famous Hoyle state (7.65 MeV resonance), which occasionally de-excites to the ground state of $^{12}C$ rather than decaying back into three alpha particles.

Without the Hoyle state resonance, the triple-alpha rate would be negligibly small and carbon production in the universe would be virtually nonexistent. This is one of the most celebrated examples of nuclear physics shaping cosmic chemistry.

Beyond $^{12}C$ formation, alpha capture continues:

• $^{12}C + ^4He \rightarrow ^{16}O + \gamma$
• $^{16}O + ^4He \rightarrow ^{20}Ne + \gamma$

The $^{12}C(\alpha,\gamma)^{16}O$ reaction rate is one of the most important and still somewhat uncertain quantities in nuclear astrophysics. It determines the carbon-to-oxygen ratio at the end of helium burning, which profoundly affects all subsequent evolution. Helium burning occurs in red giant branch and horizontal branch stars.

Carbon and Oxygen Burning

These stages occur only in massive stars ($> 8 \, M_\odot$) during their late evolutionary phases.

Carbon burning (core $T \sim 6 \times 10^8$ to $9 \times 10^8$ K):

Two $^{12}C$ nuclei fuse, and the compound nucleus $^{24}Mg^*$ decays through several channels:

• $^{12}C + ^{12}C \rightarrow ^{20}Ne + ^4He$ (most probable)
• $^{12}C + ^{12}C \rightarrow ^{23}Na + p$
• $^{12}C + ^{12}C \rightarrow ^{23}Mg + n$

The released protons, neutrons, and alpha particles drive secondary reactions that build up elements through silicon and sulfur.

Oxygen burning (core $T \sim 1.5 \times 10^9$ to $2.6 \times 10^9$ K):

Two $^{16}O$ nuclei fuse with primary channels:

• $^{16}O + ^{16}O \rightarrow ^{28}Si + ^4He$
• $^{16}O + ^{16}O \rightarrow ^{31}P + p$
• $^{16}O + ^{16}O \rightarrow ^{31}S + n$

This stage synthesizes elements up through calcium and argon. Each successive burning stage is dramatically shorter than the last because the energy yield per reaction decreases while neutrino losses (which carry energy out of the star without contributing to pressure support) increase sharply.

Silicon Burning

Silicon burning is the final exothermic burning stage, occurring at core temperatures of ~$2.7$ to $3.5 \times 10^9$ K. It does not proceed by direct fusion of two $^{28}Si$ nuclei (the Coulomb barrier would be prohibitive). Instead, it operates through a quasi-statistical equilibrium (QSE) process:

1. Energetic photons photodisintegrate $^{28}Si$ and other intermediate-mass nuclei, liberating protons, neutrons, and alpha particles.
2. These light particles are recaptured by remaining nuclei, gradually building toward the iron-peak elements.
3. The network of forward (capture) and reverse (photodisintegration) reactions approaches nuclear statistical equilibrium (NSE), driving composition toward the most tightly bound nuclei.

The products are iron-peak elements: $^{52}Cr$, $^{55}Mn$, $^{56}Fe$, $^{59}Co$, and $^{56}Ni$ (which later decays to $^{56}Fe$). In fact, the dominant product is $^{56}Ni$, not $^{56}Fe$ directly; the nickel decays to iron via two successive $\beta^+$ decays after ejection.

Since $^{56}Fe$ sits near the peak of the binding energy per nucleon curve, no further exothermic fusion is possible. The inert iron core grows until it exceeds the Chandrasekhar mass (~1.4 $M_\odot$), electron degeneracy pressure can no longer support it, and the core collapses, triggering a core-collapse supernova.

Neutron Capture Processes

Elements heavier than the iron peak cannot be produced by charged-particle fusion under stellar conditions (the Coulomb barrier is too high and the reactions are endothermic). Instead, heavy elements are built primarily through neutron capture, since neutrons carry no charge and face no Coulomb barrier.

Slow Neutron Capture (s-process)

The s-process operates in low- to intermediate-mass stars (1–8 $M_\odot$) during the thermally pulsing asymptotic giant branch (TP-AGB) phase. The defining characteristic: the timescale between successive neutron captures is long compared to typical $\beta$-decay lifetimes of unstable isotopes.

How it works:

1. Iron-group seed nuclei capture neutrons one at a time.
2. If the resulting isotope is unstable, it has time to $\beta$-decay to the next element before capturing another neutron.
3. The nucleosynthesis path therefore follows closely along the valley of beta stability on the chart of nuclides.

Main neutron sources:

• $^{13}C(\alpha, n)^{16}O$ — the primary source in low-mass AGB stars, operating during interpulse periods
• $^{22}Ne(\alpha, n)^{25}Mg$ — becomes important at higher temperatures ($T > 3 \times 10^8$ K), particularly in more massive AGB stars and during thermal pulses

The s-process accounts for roughly half of all isotopes heavier than iron. Its characteristic products include $^{88}Sr$, $^{138}Ba$, and $^{208}Pb$ (the heaviest stable s-process product). Barium and strontium abundances in stellar spectra are commonly used as s-process diagnostics.

Rapid Neutron Capture (r-process)

The r-process requires extreme neutron densities ($> 10^{20}$ neutrons/cm$^3$) and operates on timescales of seconds. These conditions are found in core-collapse supernovae (specifically the neutrino-driven wind above the proto-neutron star) and neutron star mergers.

How it works:

1. Seed nuclei capture neutrons far faster than unstable isotopes can $\beta$-decay, pushing nuclei to extremely neutron-rich territory.
2. Capture continues until reaching the neutron drip line, where the nuclear binding energy can no longer accommodate additional neutrons.
3. At these "waiting points," nuclei sit until they $\beta$-decay, increasing their proton number by one, which opens up room for further neutron captures.
4. Once the neutron flux subsides, the very neutron-rich nuclei undergo chains of $\beta$-decays back toward stability.

The r-process produces the other ~half of heavy elements beyond iron, including many that the s-process cannot reach. It is responsible for essentially all naturally occurring $^{232}Th$, $^{235}U$, and $^{238}U$.

The 2017 detection of gravitational waves from the neutron star merger GW170817, combined with the associated kilonova optical/infrared transient, provided the first direct observational evidence that neutron star mergers are a major r-process site. The relative contributions of mergers versus core-collapse supernovae remain an active area of research.

Neutron Capture Dynamics

The interplay between neutron capture rates and $\beta$-decay rates determines which path nucleosynthesis follows through the chart of nuclides.

• The neutron capture cross-section depends on neutron energy and the nuclear structure of the target. Nuclei near closed neutron shells (magic numbers $N = 50, 82, 126$) have anomalously small capture cross-sections, meaning neutrons "pile up" at these nuclei.
• s-process abundance peaks occur at stable nuclei with neutron magic numbers (e.g., $^{88}Sr$ at $N=50$, $^{138}Ba$ at $N=82$, $^{208}Pb$ at $N=126$), because these nuclei act as bottlenecks.
• r-process abundance peaks are shifted to lower mass numbers relative to the s-process peaks. This happens because the r-process path runs through very neutron-rich unstable isotopes; the waiting points occur at the magic neutron numbers, but at much lower proton number. After the neutron flux ceases and these nuclei $\beta$-decay to stability, the resulting stable isotopes have lower $A$ than the corresponding s-process peaks.
• Branching points occur along the s-process path where a particular isotope has comparable neutron capture and $\beta$-decay timescales. The branching ratio at these points is sensitive to local conditions (neutron density, temperature) and provides a diagnostic tool for constraining physical conditions in AGB interiors.

Additional Nucleosynthesis Processes

Proton Capture Process (p-process)

A small number of proton-rich stable isotopes (the "p-nuclei") cannot be produced by either the s-process or r-process, since both of those build along or beyond the neutron-rich side of stability. The p-process accounts for these rare isotopes.

The p-process operates primarily through photodisintegration in the hot oxygen/neon-rich layers of massive stars during core-collapse supernovae ($T > 2 \times 10^9$ K):

1. Pre-existing s-process and r-process seed nuclei are bombarded by energetic photons.
2. These photons strip neutrons, protons, or alpha particles from the seed nuclei via $(\gamma, n)$, $(\gamma, p)$, and $(\gamma, \alpha)$ reactions.
3. The resulting nuclei may capture protons or undergo $\beta^+$ decay, populating the proton-rich side of the valley of stability.

Examples of p-nuclei include $^{92}Mo$, $^{96}Ru$, and $^{144}Sm$. These isotopes are rare, accounting for less than 1% of heavy element abundances in the solar system. Some p-nuclei (particularly the lightest ones, like $^{92}Mo$ and $^{94}Mo$) remain difficult to produce in sufficient quantities in current models, which is an ongoing puzzle.

The rp-process (rapid proton capture) is a related but distinct mechanism occurring on the surfaces of accreting neutron stars in X-ray bursts. It builds proton-rich nuclei through successive proton captures and $\beta^+$ decays, but the material generally remains bound to the neutron star and does not contribute to galactic chemical enrichment.

Alpha Process and Photodisintegration

The alpha process refers to the successive capture of $^4He$ nuclei by heavier elements. Because alpha particles are tightly bound (high binding energy per nucleon), they are abundant at high temperatures and readily captured.

• Alpha captures on $^{12}C$, $^{16}O$, $^{20}Ne$, $^{24}Mg$, $^{28}Si$, $^{32}S$, $^{36}Ar$, and $^{40}Ca$ build up the sequence of "alpha elements" with even atomic numbers. This is why even-$Z$ elements up through the iron peak are more abundant than their odd-$Z$ neighbors (the Oddo-Harkins rule).
• The alpha process is most efficient at temperatures above $10^9$ K, in the advanced burning stages of massive stars and during explosive nucleosynthesis in supernovae.

Photodisintegration becomes increasingly important at very high temperatures ($T > 3 \times 10^9$ K), where the photon bath is energetic enough to break apart nuclei:

• $(\gamma, n)$, $(\gamma, p)$, and $(\gamma, \alpha)$ reactions compete with the corresponding capture reactions.
• During silicon burning, photodisintegration and capture reactions reach a dynamic equilibrium (NSE), and the composition is driven toward the iron peak by thermodynamics rather than by any single reaction pathway.
• In supernova shock heating, photodisintegration of iron-peak elements back into free nucleons and alpha particles absorbs enormous energy, which is a key factor in the dynamics of core collapse.

The balance between alpha capture and photodisintegration at a given temperature and density determines the equilibrium composition, and therefore the final elemental yields from massive stars and supernovae.