research32 min studyLesson 594 in Physics

Triple-alpha reactions and helium burning

Eliminate a short-lived intermediate carefully and derive the density dependence of an effective three-alpha rate.

Before you begin

On this page

The question for today

What you’ll understand

Explain how sequential reactions bridge the unstable mass-eight intermediate.
Derive the effective cubic alpha-density scaling and state its assumptions.
Track carbon production separately from later alpha capture into oxygen.

Two alpha particles do not make a stable mass-eight nucleus. Helium burning nevertheless crosses this gap because a tiny transient beryllium-eight population can capture another alpha before it breaks apart. A resonant excited state in carbon twelve greatly enhances the route that eventually emits radiation and leaves stable carbon. The familiar effective three-body reaction is therefore a compact description of sequential dynamics with a short-lived intermediate, rather than a literal simultaneous collision rule valid under all conditions.

A short-lived population can carry a finite flux

The first step continually forms and destroys beryllium eight. Its abundance can remain extremely small even while formation and destruction events are rapid. If its relaxation time is much shorter than the timescale on which alpha density and temperature change, the intermediate can be treated as approximately stationary. This is a timescale separation, not a declaration that the intermediate is absent. Setting its abundance to exactly zero would remove the carbon-forming route rather than simplify it correctly.

A minimal kinetic model assigns a two-alpha formation coefficient, a breakup rate per intermediate, and a capture coefficient for intermediate plus alpha. The formation event rate includes the identical-pair factor. Neglect reverse photodisintegration of the final carbon and other channels for the moment. Solving the intermediate balance retains capture as well as breakup in its destruction rate, revealing precisely when the commonly used cubic density law follows and when this simplified law begins to fail.

\dot n_8=\frac{k_{\alpha\alpha}}2n_\alpha^2-(\lambda_d+k_{8\alpha}n_\alpha)n_8,\qquad n_8^{\rm qs}=\frac{k_{\alpha\alpha}n_\alpha^2}{2(\lambda_d+k_{8\alpha}n_\alpha)}

The k coefficients have units volume per time, lambda_d is inverse time, and n is number per volume. This illustrative sequential network uses an intermediate quasistationary approximation with prescribed coefficients, not a complete evaluated triple-alpha rate.↔ Scroll wide equations horizontally.

The density law follows from a branching probability

The carbon production rate is the alpha-capture destruction rate times the intermediate population. Equivalently, multiply the intermediate formation rate by its probability of capturing before breaking apart. When breakup dominates, that probability is proportional to alpha density, giving an additional density factor beyond the two-body formation rate. The resulting volumetric production scales as the cube of alpha density. At fixed composition, division by mass density makes the specific heating scale as density squared.

If capture becomes comparable to or faster than breakup in this minimal model, the denominator changes and the cubic law no longer holds. The capture probability approaches one, and the production rate approaches the intermediate formation rate instead. Actual triple-alpha evaluation involves resonance energies, partial widths, thermal populations, and detailed balance; a simple coefficient model is useful for checking the algebra but is not a substitute for those nuclear inputs. The resonance contribution also creates a strong temperature sensitivity within its relevant regime.

r_{3\alpha}=\frac{k_{\alpha\alpha}k_{8\alpha}n_\alpha^3}{2(\lambda_d+k_{8\alpha}n_\alpha)},\qquad k_{8\alpha}n_\alpha\ll\lambda_d\ \Longrightarrow\ r_{3\alpha}\propto n_\alpha^3,\quad \epsilon_{3\alpha}=\frac{Q_{3\alpha}r_{3\alpha}}{\rho}

r counts carbon-forming events per volume per time. Q is energy per net three-alpha conversion and rho is mass per volume. Specific heating epsilon has units energy per mass per time when radiative products thermalize.↔ Scroll wide equations horizontally.

Carbon production and carbon survival are different questions

Three alpha particles are removed net for every newly produced carbon nucleus in the reduced reaction. Temporary breakup returns the first two alpha particles, so it must not be counted as permanent helium consumption. Once carbon exists, capture of another alpha can produce oxygen sixteen. Carbon abundance therefore results from a competition between production and destruction, coupled to the changing helium reservoir and the thermal history. A larger triple-alpha rate alone does not determine a final stellar carbon-to-oxygen ratio.

Thermal feedback further separates a rate calculation from a stellar outcome. In a nondegenerate region, heating can expand the gas and reduce the conditions favoring burning; an electron-degenerate environment can respond differently. Predicting an instability requires the equation of state, energy transport, and cooling as well as a temperature derivative of the nuclear rate. The exercises below deliberately hold the environment fixed and use rescaled coefficients, allowing the intermediate bookkeeping to be inspected without making a universal stellar-yield or flash prediction.

\dot n_\alpha\big|_{3\alpha}=-3r_{3\alpha},\qquad \dot n_{12}=r_{3\alpha}-n_{12}n_\alpha\langle\sigma v\rangle_{12,\alpha}

The second equation includes only carbon production and alpha capture to oxygen in a fixed-volume illustrative parcel. Other captures, reverse reactions, expansion, and mixing are omitted; each term has units number per volume per time.↔ Scroll wide equations horizontally.

Work it through

A rescaled intermediate with competing exits

For an explicitly rescaled kinetic example, the intermediate formation rate is $1000\,\mathrm{cm}^{-3}\mathrm{s}^{-1}$ , its breakup rate is $99\,\mathrm{s}^{-1}$ , and its capture rate per intermediate is $1\,\mathrm{s}^{-1}$ . These are not measured beryllium-eight rates. Find the quasistationary intermediate density, carbon production, and net alpha loss.

The total intermediate destruction rate is $99+1=100\,\mathrm{s}^{-1}$ . Its relaxation time is therefore $0.01\,\mathrm s$ , so the stationary approximation requires the prescribed surroundings to vary much more slowly than this.
Balance formation and destruction to obtain $n_8=1000/100=10\,\mathrm{cm}^{-3}$ . Only one percent of the formed intermediates capture before breakup, giving carbon production $r_{3\alpha}=1(10)=10\,\mathrm{cm}^{-3}\mathrm{s}^{-1}$ .
Each successful net event consumes three alphas, so $\dot n_\alpha=-30\,\mathrm{cm}^{-3}\mathrm{s}^{-1}$ . The much larger intermediate formation rate does not imply permanent loss of two alphas in every formation event because most intermediates break apart again.

The result

The intermediate density is ten per cubic centimetre, producing ten carbon nuclei and consuming thirty alpha particles per cubic centimetre per second.

The living notebook

Explore the model

Interactive lab

Interactive lab · function

Specific heating scales quadratically in density

At fixed temperature and helium fraction in the breakup-dominated regime, identify position with dimensionless mass density. Reset and set coefficient and position to $(1,1)$ , then independently $(1,2)$ . Predict specific-heating proxies 1 and 4. This quadratic graph illustrates the density-squared specific rate only; it does not calculate the cubic volumetric rate, the Hoyle resonance, an intermediate population, or a stellar flash.

y = ax²TangentSigned integral

y at chosen x: 0.608
Tangent slope · 2ax: 1.35
Signed area · ax³/3: 0.182

Coefficient a0.75Position x0.9

t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forVolumetric and per-mass heating rates differ by one density factor.

Pause & reason

Check your understanding

Interactive checks need JavaScript. You can still read the explanation.

Read the checkpoint explanation

Formation supplies two density factors and a small capture-before-breakup probability supplies a third. The stationary approximation also needs fast intermediate relaxation relative to changes of the surroundings.

Your turn

Try it yourself

Double alpha density in the rescaled worked model while keeping both two-body coefficients and the breakup rate fixed. Recompute the stationary carbon rate and its ratio to the original, rather than assuming exact cubic scaling.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Formation grows by four, but capture per intermediate doubles. Include both destruction channels in the new denominator.

Compare with the solution

The formation rate becomes $4000\,\mathrm{cm}^{-3}\mathrm{s}^{-1}$ , and the capture rate becomes $2\,\mathrm{s}^{-1}$ . The intermediate density is $4000/101\,\mathrm{cm}^{-3}$ and the carbon rate is $8000/101\simeq79.21\,\mathrm{cm}^{-3}\mathrm{s}^{-1}$ . Its ratio to the original ten is $800/101\simeq7.921$ , slightly below eight because capture is a larger fraction of intermediate destruction.

Keep the idea

Take it with you

A tiny short-lived intermediate can sustain a finite carbon-production flux.
Cubic volumetric scaling assumes breakup-dominated intermediate balance; specific heating has one fewer density factor.
Carbon survival and stellar stability require further reactions and environmental feedback.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.