Technopedia

The technologist’s encyclopedia

Nuclear physics

Nuclear physics describes the structure and behavior of atomic nuclei: the strong-interaction bound states of protons and neutrons (nucleons), the reactions that transform them, and the radioactive decays that reveal the weak interaction. It is relevant whenever the nucleus itself matters: energy production (fission reactors, fusion), medical imaging and radiotherapy, radiometric dating, and the nucleosynthesis that built the elements (see Cosmology).

The nucleus sits between two levels of description. Below it, nucleons are bound states of quarks and gluons (see Elementary particles); the force between nucleons is a residual of the strong interaction, just as the van der Waals force between molecules is a residual of electromagnetism. Above it, the nucleus is simply the positive charge +Ze+Ze that anchors the electron shells of atomic physics. And because nuclei are quantum objects, the tools come from quantum mechanics: bound states and shell filling from the Pauli principle, and decay through barriers from tunneling.

Authoritative sources

The nuclear force

Definition. The nuclear force is the residual strong interaction between nucleons: a short-range attraction (range ≈ 1–2 fm) with a repulsive core, approximately charge-independent, plus the long-range Coulomb repulsion between protons.

Details

Where the force comes from

Quantum chromodynamics confines quarks into protons and neutrons (see Elementary particles). Each nucleon is color-neutral, so — like two neutral atoms — two nucleons feel almost no force at large distance. But their quark clouds polarize each other, leaving a residual attraction at short range, analogous to the van der Waals force between molecules (see molecular bonds). This residual force is what binds nuclei.

Properties

The potential between two nucleons

For two protons, the effective potential is the sum of the short-range nuclear well and the long-range Coulomb repulsion VC(r)=1.44 MeV\cdotpfm/rV_C(r) = 1.44 \text{ MeV·fm}/r:

V(r)≈−V0 e−r/λr/λ+1.44 MeV\cdotpfmrV(r) \approx -V_0\, \frac{e^{-r/\lambda}}{r/\lambda} + \frac{1.44 \text{ MeV·fm}}{r} V (MeV) r (fm) +35 0 −50 1 2 3 4 5 −2.2 MeV (deuteron) nuclear well (Yukawa): depth ≈ 40–100 MeV, range ≈ 1–2 fm Coulomb repulsion 1.44 MeV·fm / r (small for p–p, but ≈ 30 MeV for α + heavy daughter) repulsive core

The Coulomb tail looks tiny on this scale — but it is the whole story of alpha decay and fission: for an α\alpha particle facing a heavy daughter (Z≈90Z \approx 90), the same 1/r1/r term reaches ≈ 30 MeV at the nuclear surface, forming a barrier that the α\alpha can only escape by quantum tunneling.

Example: why the deuteron is fragile

The deuteron’s binding energy (2.22 MeV) is tiny compared with typical nuclear wells (~40 MeV) because it is bound only by the tensor part of the force in the spin-triplet channel — there is no bound p–p or n–n nucleus. This is direct evidence that the nuclear force depends on spin and isospin, not just position.

Authoritative sources

Example problem

Yukawa proposed that a force of range λ\lambda is carried by a particle of mass m=ℏ/(λc)m = \hbar/(\lambda c). The nuclear force has range ≈ 1.4 fm. Estimate the mass of its carrier, in MeV/c2c^2. (Use ℏc=197.3\hbar c = 197.3 MeV·fm.)

Solution mc2=ℏcλ=197.3 MeV\cdotpfm1.4 fm≈141 MeVm c^2 = \frac{\hbar c}{\lambda} = \frac{197.3 \text{ MeV·fm}}{1.4 \text{ fm}} \approx 141 \text{ MeV}

This matches the pion mass (mπc2≈140m_{\pi} c^2 \approx 140 MeV) — the prediction that led to the discovery of the pion in 1947, and the first concrete picture of a force as particle exchange (see Elementary particles).

Nuclear structure: liquid drop and shell models

Definition. A nucleus with ZZ protons and NN neutrons (mass number A=Z+NA = Z + N) is a quantum bound state of AA nucleons. Its binding energy is described by two complementary models: the liquid drop (collective, like an incompressible fluid) and the shell model (independent nucleons filling quantized orbits, like electrons in atoms).

Details

Step 1: the liquid drop model

Because the nuclear force is short-ranged and saturating, a nucleus behaves like a droplet of incompressible fluid: each nucleon contributes a fixed volume energy, only surface nucleons lose neighbors, protons repel each other via Coulomb, and the Pauli exclusion principle penalizes proton–neutron imbalance. The Bethe–Weizsäcker semi-empirical mass formula:

B(A,Z)=aVA−aSA2/3−aCZ(Z−1)A1/3−aA(A−2Z)2A+δ(A,Z)B(A,Z) = a_V A - a_S A^{2/3} - a_C \frac{Z(Z-1)}{A^{1/3}} - a_A \frac{(A-2Z)^2}{A} + \delta(A,Z)

with aV≈15.8a_V \approx 15.8 MeV (volume), aS≈18.3a_S \approx 18.3 MeV (surface), aC≈0.714a_C \approx 0.714 MeV (Coulomb), aA≈23.2a_A \approx 23.2 MeV (asymmetry), and a pairing term δ=±12/A\delta = \pm 12/\sqrt{A} (positive for even–even nuclei, negative for odd–odd, zero for odd AA).

Example. For iron-56 (A=56A = 56, Z=26Z = 26): volume term 15.8×56=884.815.8 \times 56 = 884.8; surface 18.3×562/3=267.918.3 \times 56^{2/3} = 267.9; Coulomb 0.714×26×25/561/3=121.30.714 \times 26 \times 25 / 56^{1/3} = 121.3; asymmetry 23.2×16/56=6.623.2 \times 16/56 = 6.6; pairing +1.6+1.6. Total B≈490.6B \approx 490.6 MeV, i.e. B/A≈8.76B/A \approx 8.76 MeV — within 0.4% of the measured 8.79 MeV.

Step 2: where the liquid drop fails — magic numbers

The liquid drop formula predicts binding energy varying smoothly with ZZ and NN. But experiment shows certain nuclei are exceptionally stable and abundant: those with

N or Z=2, 8, 20, 28, 50, 82, 126N \text{ or } Z = 2,\ 8,\ 20,\ 28,\ 50,\ 82,\ 126

the magic numbers. Nuclei with both NN and ZZ magic (doubly magic: 4{}^4He, 16{}^{16}O, 40{}^{40}Ca, 48{}^{48}Ca, 208{}^{208}Pb) are spherical, rigid, and tightly bound. A smooth droplet model cannot explain this — just as classical physics could not explain the discrete spectral lines that revealed electron shells.

Step 3: the shell model

The nuclear shell model (Mayer and Jensen, 1949) treats each nucleon as moving independently in the average potential well created by all the others (a Woods–Saxon well, similar in spirit to the electron shells of atomic physics). Because nucleons are fermions, the Pauli exclusion principle forces them to fill discrete levels: closed shells correspond to the magic numbers.

The crucial difference from atoms is the spin–orbit coupling: the nuclear spin–orbit splitting is strong and has the opposite sign to the atomic case (where it is a small relativistic effect). A level with orbital angular momentum ℓ\ell splits into j=ℓ+1/2j = \ell + 1/2 (lower) and j=ℓ−1/2j = \ell - 1/2 (higher), and the large gaps after the lowered j=ℓ+1/2j = \ell + 1/2 levels reproduce the magic numbers 28, 50, 82, and 126.

Spin–orbit splitting of the 1f shell (magic number 28) 1f (no spin–orbit) 10 states (ℓ=3, spin ±½) 1f₅⁄₂ (higher) 6 states 1f₇⁄₂ (lower) 8 states strong splitting, opposite sign to atoms ← filling these 8 states gives the gap at Z, N = 28

The shell model explains why 4{}^{4}He (closed 1s shell), 16{}^{16}O (closed 1p shell), and the calcium isotopes (Z=20Z = 20) are so stable, and why neutron shell closures at N=126N = 126 shape the abundance peaks of the heavy elements formed in nucleosynthesis.

Authoritative sources

Example problem

The bottom of the valley of stability is well approximated by Z=A/(2+aC2aAA2/3)Z = A / \left(2 + \frac{a_C}{2 a_A} A^{2/3}\right). Predict the most stable ZZ for A=120A = 120 (use aC=0.714a_C = 0.714, aA=23.2a_A = 23.2), and identify the element.

Solution

A2/3=1202/3≈24.3A^{2/3} = 120^{2/3} \approx 24.3, so the denominator is 2+0.71446.4×24.3=2+0.374=2.3742 + \frac{0.714}{46.4} \times 24.3 = 2 + 0.374 = 2.374.

Z=1202.374≈50.5⇒Z=50Z = \frac{120}{2.374} \approx 50.5 \Rightarrow Z = 50

Z=50Z = 50 is tin — and 50 is a magic number. Tin has the largest number of stable isotopes (10) of any element, a direct signature of the shell closure.

Binding energy: the energy landscape of fission and fusion

Definition. The binding energy per nucleon B/AB/A measures how tightly a nucleus is bound. Energy is released by any rearrangement — fission or fusion — that moves nuclei up the B/AB/A curve toward its maximum at iron.

Details

Step 1: the curve

Plotting B/AB/A against mass number AA gives the most important single picture in nuclear energy physics:

mass number A B/A (MeV per nucleon) 0 2 4 6 8 10 0 50 100 150 200 250 peak: Fe-56 / Ni-62 ≈ 8.8 MeV fusion: light nuclei climb toward the peak fission: heavy nuclei descend toward the peak

Three features to read off the curve:

  1. Rapid rise up to A≈20A \approx 20: each added nucleon is grabbed by the short-range force (this is why nucleosynthesis in stars releases energy).
  2. Broad peak at iron (B/A≈8.8B/A \approx 8.8 MeV): iron-group nuclei are the most tightly bound — the ash of nuclear burning, and the endpoint of exothermic fusion.
  3. Slow decline beyond A≈60A \approx 60: the Coulomb repulsion of the growing proton population wins out (the Z2/A1/3Z^2/A^{1/3} term in the mass formula), so very heavy nuclei are less bound per nucleon and can release energy by splitting.

Step 2: reading energy release off the curve

The energy released by a rearrangement is the total binding energy gained:

Ereleased=Δ ⁣(BA)×AE_{\text{released}} = \Delta\!\left(\frac{B}{A}\right) \times A

Step 3: energy per kilogram

Because nuclear binding energies (~MeV per nucleon) are ~10⁶ times larger than chemical bond energies (~eV per molecule), nuclear fuels are enormously more energy-dense than chemical ones:

Energy sourceEnergy density (J/kg)Mechanism
Coal (chemical)≈ 3 × 10⁷Electron rearrangement (eV-scale bonds)
U-235 fission≈ 8 × 10¹³Heavy nucleus splits toward the peak of the curve
D–T fusion≈ 3.4 × 10¹⁴Light nuclei climb toward the peak of the curve

Applications. Fission reactors exploit the neutron multiplication of the fission reaction (see Nuclear reactions) to sustain a controlled chain reaction. Fusion powers the stars (see Cosmology) and drives terrestrial fusion-energy research, where the D–T reaction is the leading candidate because it has the largest cross-section at reachable temperatures.

Authoritative sources

Example problem

Verify the energy density of D–T fusion fuel. One D–T reaction releases 17.6 MeV; the reactants (deuteron + triton) have a total mass of 5 u. (1 u = 1.6605 × 10⁻²⁷ kg; 1 MeV = 1.602 × 10⁻¹³ J.)

Solution

Energy per reaction: 17.6×1.602×10−13=2.82×10−1217.6 \times 1.602 \times 10^{-13} = 2.82 \times 10^{-12} J.

Mass per reaction: 5×1.6605×10−27=8.30×10−275 \times 1.6605 \times 10^{-27} = 8.30 \times 10^{-27} kg.

Em=2.82×10−128.30×10−27≈3.4×1014 J/kg\frac{E}{m} = \frac{2.82 \times 10^{-12}}{8.30 \times 10^{-27}} \approx 3.4 \times 10^{14} \text{ J/kg}

Compare coal at 3×1073 \times 10^7 J/kg: fusion fuel is about 10 million times more energy-dense — the direct macroscopic consequence of the MeV-vs-eV scale difference between nuclear and chemical binding.

Nuclear reactions

Definition. A nuclear reaction transforms one set of nuclei into another, conserving charge, baryon number, lepton number, energy, and momentum. Its energy balance is the Q-value: Q=(mi−mf)c2Q = (m_i - m_f)c^2, positive for exothermic reactions.

Details

Step 1: the conservation laws

Whatever happens in a nuclear reaction, the books must balance: total charge ZZ, total baryon number AA (for nucleons), lepton number, energy, and momentum are conserved. Writing a reaction is therefore bookkeeping: match the superscripts (mass numbers) and subscripts (atomic numbers) on both sides.

Step 2: the main reaction types

Step 3: the Q-value

The Q-value is the rest-mass difference converted to kinetic energy:

Q=(∑mi−∑mf)c2Q = \left(\sum m_i - \sum m_f\right) c^2

Example. For the D–T reaction, using atomic masses m(2H)=2.014102m({}^{2}\text{H}) = 2.014102 u, m(3H)=3.016049m({}^{3}\text{H}) = 3.016049 u, m(4He)=4.002603m({}^{4}\text{He}) = 4.002603 u, m(n)=1.008665m(n) = 1.008665 u: mi=5.017522m_i = 5.017522 u, mf=5.011268m_f = 5.011268 u, so Q=0.018883×931.5≈17.6Q = 0.018883 \times 931.5 \approx 17.6 MeV — exothermic, as the binding-energy curve predicts.

Authoritative sources

Example problem

The proton–proton chain in the Sun converts four protons into one helium-4 nucleus, two positrons, two neutrinos, and energy. Using m(4He)=4.002603m({}^{4}\text{He}) = 4.002603 u and m(p)=1.007276m(p) = 1.007276 u, compute the total energy released per helium-4 nucleus formed. (Ignore the neutrino rest masses; 1 u = 931.5 MeV/c2c^2.)

Solution

Initial mass: 4×1.007276=4.0291044 \times 1.007276 = 4.029104 u.

Final rest mass: 4.0026034.002603 u — the helium-4 atomic mass already accounts for the two electrons whose places the two positrons effectively take (the positron masses cancel the missing electron masses in the bookkeeping).

Q=(4.029104−4.002603)×931.5=0.026501×931.5≈24.7 MeVQ = (4.029104 - 4.002603) \times 931.5 = 0.026501 \times 931.5 \approx 24.7 \text{ MeV}

Each positron then annihilates with an electron, releasing 2×0.511=1.0222 \times 0.511 = 1.022 MeV; the two annihilations add 2×1.022=2.042 \times 1.022 = 2.04 MeV, bringing the total to 24.7+2.0≈26.724.7 + 2.0 \approx 26.7 MeV — the figure quoted for the Sun’s energy source (see Cosmology).

Radioactivity

Definition. Radioactivity is the spontaneous transformation of an unstable nucleus: alpha decay (strong + Coulomb, tunneling), beta decay (weak interaction), or gamma decay (electromagnetic). The decay rate follows N(t)=N0e−λtN(t) = N_0 e^{-\lambda t} with half-life t1/2=ln⁡2/λt_{1/2} = \ln 2/\lambda.

Details

Step 1: alpha decay — tunneling through the Coulomb barrier

An α\alpha particle (a helium-4 nucleus) is preformed inside a heavy nucleus, trapped by the Coulomb barrier described in the nuclear force section. Classically it could never escape; quantum mechanically it tunnels through the barrier with a tiny but nonzero probability per collision — which is why alpha half-lives range from microseconds to billions of years for nearly the same decay energy (the Geiger–Nuttall law).

r (fm) V (MeV) 5 10 15 20 Eα = 4.2 MeV α preformed in the well Coulomb barrier ≈ 29 MeV at R ≈ 9 fm tunneling region

Example. 238{}^{238}U →234\to {}^{234}Th +α+4.27+ \alpha + 4.27 MeV, with t1/2=4.5×109t_{1/2} = 4.5 \times 10^9 yr — long because the 4.2 MeV α\alpha sits far below the ≈ 29 MeV barrier top.

Step 2: beta decay — the weak interaction at work

In beta-minus decay a neutron converts to a proton via the weak interaction, emitting an electron and an antineutrino through a W−W^- boson (see Elementary particles):

n→p+e−+νˉen \to p + e^- + \bar{\nu}_e

Beta decay moves a nucleus along the valley of stability (changing NN and ZZ at fixed AA) and is the dominant decay mode of the neutron-rich fission fragments. The neutrino carries away variable energy — the observation of its continuous electron spectrum was the first evidence for the neutrino.

Step 3: gamma decay — photons from nuclear levels

After alpha or beta decay the daughter nucleus is often left in an excited state. It sheds the excitation by emitting a high-energy photon — nuclear electromagnetic radiation at the MeV scale. Gamma decay changes neither ZZ nor AA.

Example. 99m{}^{99m}Tc →99\to {}^{99}Tc +γ+ \gamma (140 keV, t1/2=6.0t_{1/2} = 6.0 h): the workhorse of nuclear medicine, injected into patients for imaging because the 140 keV photon escapes the body and is easily detected.

Step 4: a complete decay scheme

The classic decay scheme of cobalt-60 ties all three ideas together: the beta decay creates the excited daughter, and the gamma cascade de-excites it.

Decay scheme of cobalt-60 ⁶⁰Co 0 (5⁺) t½ = 5.27 y 2505.7 keV (4⁺) 1332.5 keV (2⁺) ⁶⁰Ni 0 (0⁺) β⁻ (Emax = 0.318 MeV) γ 1173 keV γ 1332 keV

Note the division of labor: the weak interaction (β⁻) changes the species of nucleus; the electromagnetic interaction (the two γ’s) then drains the excitation energy of the daughter 60{}^{60}Ni. The two sharp gamma lines (1173 and 1332 keV) are used to calibrate gamma detectors.

Step 5: the decay law and radiometric dating

Each unstable nucleus decays with a fixed probability per unit time λ\lambda, independent of age and chemistry:

N(t)=N0e−λt,t1/2=ln⁡2λN(t) = N_0 e^{-\lambda t}, \qquad t_{1/2} = \frac{\ln 2}{\lambda}

Example. Carbon-14 (t1/2=5730t_{1/2} = 5730 yr) is produced in the atmosphere and taken up by living organisms. When they die, the intake stops and the 14{}^{14}C/12{}^{12}C ratio decays exponentially — measuring it dates organic remains up to ~50 000 years (see the example problem below).

Authoritative sources

Example problem

An archaeological wooden sample has 14C/12C{}^{14}\text{C}/{}^{12}\text{C} activity equal to one quarter of that of living wood. How old is it? (t1/2=5730t_{1/2} = 5730 yr.)

Solution

The activity is proportional to the number of 14{}^{14}C nuclei, so N/N0=1/4=e−λtN/N_0 = 1/4 = e^{-\lambda t}.

t=ln⁡4λ=2 t1/2=2×5730≈11 460 yrt = \frac{\ln 4}{\lambda} = 2\, t_{1/2} = 2 \times 5730 \approx 11\,460 \text{ yr}

Each halving of the activity adds one half-life: two halvings → two half-lives.

Interactive simulation

Explore the binding-energy curve interactively: slide the mass number to see where a nucleus sits on the curve, and read off the energy released by fissioning it or by fusing two copies of it. Near iron, neither direction pays.

Quiz

Example problem

A uranium-235 nucleus absorbs a thermal neutron and fissions into barium-141 and krypton-92, releasing 3 neutrons. Given the atomic masses: m(235U)=235.044m({}^{235}\text{U}) = 235.044 u, m(n)=1.009m(n) = 1.009 u, m(141Ba)=140.914m({}^{141}\text{Ba}) = 140.914 u, m(92Kr)=91.926m({}^{92}\text{Kr}) = 91.926 u, compute the Q-value of the reaction in MeV. (1 u = 931.5 MeV/c2c^2.)

Solution

The reaction is: n+235U→141Ba+92Kr+3nn + {}^{235}\text{U} \to {}^{141}\text{Ba} + {}^{92}\text{Kr} + 3n.

Initial mass: mi=m(n)+m(235U)=1.009+235.044=236.053m_i = m(n) + m({}^{235}\text{U}) = 1.009 + 235.044 = 236.053 u.

Final mass: mf=m(141Ba)+m(92Kr)+3 m(n)=140.914+91.926+3(1.009)=235.867m_f = m({}^{141}\text{Ba}) + m({}^{92}\text{Kr}) + 3\,m(n) = 140.914 + 91.926 + 3(1.009) = 235.867 u.

Mass defect: Δm=mi−mf=236.053−235.867=0.186\Delta m = m_i - m_f = 236.053 - 235.867 = 0.186 u.

Q-value: Q=Δm×931.5=0.186×931.5≈173Q = \Delta m \times 931.5 = 0.186 \times 931.5 \approx 173 MeV.

Two remarks. First, using atomic masses is legitimate here: the 92 electrons of the uranium atom cancel exactly against the 56 + 36 = 92 electrons of the barium and krypton atoms, because total ZZ is conserved. Second, the measured fission energy is ≈ 200 MeV rather than 173 MeV because the neutron-rich fragments 141{}^{141}Ba and 92{}^{92}Kr are radioactive: their subsequent beta decays and gamma emissions add ≈ 30 MeV, and real fission produces a range of fragment splits, not just this symmetric-ish channel.