# Nuclear Reactions

The primary source of the Sun's power output is the proton-proton chain:
• 1H + 1H -> 2H + e+ + νe
• 2H + 1H -> 3He + γ
• 3He + 3He -> 4He + 1H + 1H

(γ denotes a photon.) The first two of these reactions occurs twice for each occurrence of the third reaction.

In these nuclear reactions, the difference between the masses of the products and reactants can be used to compute the energy released by using
energy = mass * c2
Most people think of this as "converting matter into energy". Since energy is a quality of matter, and not really a separate entity, it is better to think of this process in terms of the constituent quarks of the protons and neutrons involved. Then we can see that each reactant has an energy content associated with the fact that they are bound states of particles which are electrically as well as strongly interacting. The product has a different binding energy, and the fusion process liberates the difference as gamma rays (with electromagnetic energy) and thermal energy of motion. The energy of any neutrinos which are produced is usually included as part of the energy released, because the neutrino has a very small and not very accurately measured mass in comparison with the other particles.
The mass of a proton (1H) is 1.673 * 10-27 kg; the mass of a Helium nucleus (4He) is 6.647 * 10-27 kg, and the mass of a positron (e+) is 9.109 * 10-31 kg. The energy released in neutrinos (νe) and gamma rays (γ) is then
E = (4 * 1.673 * 10-27 - 6.647 * 10-27 - 2 * 9.109 * 10-31) * (3 * 108)2
= 3.881 * 10-12 J

(= 24.225 MeV)

Using the solar constant (1365 W/m2) as the energy flux and 1 AU as the distance in the equation below, we can compute the luminosity of the Sun, obtaining LS = 3.84 * 1026 W. Dividing the luminosity by the energy per reaction gives the number of reactions per second, 9.92 * 1037. Since each reaction consumes 4 protons, the rate of Hydrogen consumption is
9.92 * 1037 * 4 * 1.673 * 10-27
= 6.64 * 1011 kg / s.
The photons released in these reactions are gamma rays, but the bulk of the electromagnetic radiation leaving the surface of the Sun is in the visible spectrum. The average travel distance (mean free path) between collisions for a photon in the center of the Sun is about 1 cm. On average, each photon leaving the surface has interacted with matter about 1021 times since being emitted in the core. Since many of these interactions involve the absorption and re-emission of the photon, the photons emitted at the surface are really distant ancestors of the ones created during fusion. Most interactions involve some energy loss (they are "inelastic"), and this accounts for the shift in spectrum from gamma to visible. The average time for a photon to travel to the surface from the core is
1021 * .01 m / c = 33 billion seconds
or over a thousand years.

In the following list of nuclear reactions occurring in stars, numbers in parentheses are approximate energy yields in MeV; negative values indicate additional energy is required to make these reactions occur. Note that these are not all of the possible reactions; any reaction which conserves energy, momentum, angular momentum, electric charge, baryon number and lepton number can occur. Each reaction occurs with different frequency based on how much energy is required to overcome the electric repulsion between nuclei before the strong force takes over.

• Hydrogen fusion ignites at 107 K (pp chain)
• 1H + 1H -> 2H + e+ + νe (0.38)
• 2H + 1H -> 3He + γ (5.5)
• 3He + 3He -> 4He + 1H + 1H (12.9)
This process takes place not only in the interior of main sequence stars, but on the surfaces of white dwarfs that have accreted Hydrogen from a binary partner. This is called a nova, and the process usually repeats over time.
or (CNO Cycle, dominant beyond 16 MK)

• 12C + 1H -> 13N + γ (2)
• 13N -> 13C + e+ + νe (1.1)*
• 13C + 1H -> 14N + γ (11.2)
• 14N + 1H -> 15O + γ (3.6)
• 15O -> 15N + e+ + νe (1.8)*
• 15N + 1H -> 12C + 4He (4.9)
• Helium fusion ignites at 108 K (Triple Alpha Process)
• 4He + 4He -> 8Be + γ (-0.09)
• 4He + 8Be ->12C + γ (7.4)
• Carbon fusion ignites at 2-6 * 108 K
• 12C + 4He -> 16O + γ (7.2)
• 12C + 12C -> 20Ne + 4He (4.7)
• 12C + 12C -> 23Na + 1H (2.2)
• 12C + 12C -> 23Mg + 1n (-2.6)
• 12C + 12C -> 24Mg + γ (14)

Again, this process does not just occur during the horizontal branch of a massive star. If a Carbon white dwarf accretes enough material from a binary companion to raise its internal temperature sufficiently, the white dwarf will begin Carbon fusion throughout its mass and explode in a Type Ia supernova (Type I spectra lack Hydrogen lines; Type Ia typically have Silicon lines).

The second to last reaction in this set is included to illustrate reactions which release neutrons. These free neutrons will enter into neutron capture reactions below.

• Oxygen fusion ignites at 1.4-2 * 109 K
• 16O + 4He -> 20Ne + γ (4.8)
• 16O + 16O -> 28Si + 4He (0.47)
• Silicon fusion ignites at 3.5 * 109 K
• 28Si + 4He -> 32S + γ (16)
• 32S + 4He -> 36Ar + γ (6.7)
• 36Ar + 4He -> 40Ca + γ (7)
• 40Ca + 4He -> 44Ti + γ (5.1)
• 44Ti + 4He -> 48Cr + γ (7.7)
• 48Cr + 4He -> 52Fe + γ (7.9)
• 52Fe + 4He -> 56Ni + γ (8)
• 28Si + 28Si -> 56Ni + γ (29.2)
• 56Ni -> 56Co + e+ + νe (1.1)*
• 56Co -> 56Fe + e+ + νe (3.5)*
Note that baryon number conservation accounts for the drop in the number of nuclei available to fuse. For instance, it takes 56 1H to (eventually) make just one 56Fe.
The following masses (in atomic mass units (amu); 1 amu = 1.66055 * 10-27 kg) will enable you to compute more accurate energy yields from these reactions.

 n 1.0087 14N 13.9954 32S 31.9633 1H 1.00725 15N 14.9963 36Ar 35.9576 2H 2.01355 15O 14.9987 40Ca 39.9516 3He 3.0149 16O 15.9905 44Ti 43.9476 4He 4.0015 20Ne 19.9869 48Cr 47.9408 8Be 8.00311 23Na 22.9838 52Fe 51.9338 12C 11.9967 23Mg 22.9875 56Fe 55.9206 13C 13.0001 24Mg 23.9784 56Co 55.925 13N 13.0019 28Si 27.979 56Ni 55.9267

You can find the atomic numbers of these atoms from the Periodic Table of the Elements.

Neutron capture can form atoms with atomic masses larger than Iron, for instance in the series of processes:
56Fe + 1n -> 57Fe

57Fe + 1n -> 58Fe

58Fe + 1n -> 59Fe

59Fe -> 59Co + e- + νe

59Co + 1n -> 60Co

60Co -> 60Ni + e- + νe

Observation of the merger of two neutron stars on August 17, 2017 seems to indicate that such mergers may be the primary source of heavy elements in the universe.
But nuclear fusion ends with Iron production. This is because 56Fe fusion is an endothermic process: it takes more energy than it releases. For instance,
56Fe + 56Fe -> 112Te
requires over 44 MeV to take place. While we have seen several such reactions above, they have proceeded because excess energy is available in the active core of the star. Once an Iron core forms, there are no reactions available to provide the additional energy necessary for Iron fusion to occur.
It is often said that 56Fe is the most stable nucleus. A nucleus is said to be stable if its energy per baryon is a minimum for its atomic number. Unstable nuclei undergo radioactive decay, emitting Helium nuclei (alpha decay, which heavier nuclei often do because they lose energy quicker), electrons or positrons (beta decay, accompanied by neutrinos), or photons (gamma decay, if nothing else is allowed). The reactions marked with an asterisk above are all beta decays. Fusion stops with Iron because 56Fe has the lowest energy per baryon of any nucleus (-8.8 MeV / baryon).

Portfolio Exercise: For each of the fusion reactions above, compute the energy yield.

This table shows the solar system abundances of a number of elements and isotopes (the vertical axis is log base 10, so, for instance, since 4He ≈ 3He + 4, there is about 104 = 10,000 times more 4He than 3He):