← Part III: Nuclear Reactions
Chapter 8

Nuclear Fission

Fission Physics

Nuclear fission is the splitting of a heavy nucleus into two (or rarely three) lighter fragments, releasing approximately 200 MeV of energy per event. The process was discovered by Hahn and Strassmann in 1938, with the theoretical explanation provided by Meitner and Frisch using the liquid drop model.

$$^{235}_{92}\text{U} + n \to ^{236}_{92}\text{U}^* \to \text{fragments} + 2\text{-}3\,n + \sim200\text{ MeV}$$

Energy Budget per Fission

Fission fragment kinetic energy:

~170 MeV

Prompt neutron kinetic energy:

~5 MeV

Prompt gamma rays:

~7 MeV

Beta particles from products:

~8 MeV

Gamma rays from products:

~7 MeV

Neutrinos (lost):

~12 MeV

Total:

~200 MeV

Fission Barrier

The liquid drop model predicts a fission barrier arising from the competition between surface energy (resists deformation) and Coulomb energy (favors deformation). The barrier height depends on the fissility parameter:

$$x = \frac{E_C^{(0)}}{2E_S^{(0)}} = \frac{Z^2/A}{(Z^2/A)_{\text{crit}}} \approx \frac{Z^2/A}{49}$$

For $^{236}$U* (compound nucleus): $Z^2/A \approx 36$, $x \approx 0.73$, barrier height ~6 MeV. The neutron separation energy of $^{236}$U is 6.5 MeV, so thermal neutrons can induce fission of $^{235}$U.

For $^{239}$U* (from $^{238}$U + n): the neutron separation energy is only 4.8 MeV, less than the barrier height (~6 MeV). Hence $^{238}$U requires fast neutrons ($E_n > 1$ MeV) for fission.

Chain Reactions and Criticality

The neutron multiplication factor $k$ determines whether a chain reaction is self-sustaining:

$$k = \frac{\text{neutrons in generation } n+1}{\text{neutrons in generation } n}$$
  • - $k < 1$: Subcritical -- chain reaction dies out
  • - $k = 1$: Critical -- steady-state chain reaction (reactor operation)
  • - $k > 1$: Supercritical -- exponentially growing chain reaction

The four-factor formula for an infinite medium is:

$$k_\infty = \eta\,f\,p\,\epsilon$$

where $\eta$ = neutrons per absorption in fuel, $f$ = thermal utilization factor,$p$ = resonance escape probability, $\epsilon$ = fast fission factor.

Python Simulation: Fission Product Yields

The characteristic double-humped fission product mass distribution and neutron multiplicity as a function of incident neutron energy.

Fission Product Yield Distribution

Python

Double-humped mass distribution for U-235 and Pu-239 thermal fission

script.py97 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran Implementation

Criticality calculation: computes k-effective using the four-factor formula.

Reactor Criticality Calculation

Fortran

Four-factor formula k-effective calculation for various enrichment levels

criticality_calc.f9068 lines

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Derivation: Fissility Parameter from the Liquid Drop Model

The liquid drop model treats the nucleus as an incompressible charged fluid. Fission occurs when the Coulomb repulsion overcomes the surface tension. We derive the critical condition for instability against small deformations.

Step 1: Parameterize the Deformation

Consider an axially symmetric quadrupole deformation of a sphere of radius $R_0$. The nuclear surface is described by:

$$R(\theta) = R_0\left[1 + \alpha_2 P_2(\cos\theta)\right]$$

where $\alpha_2$ is the quadrupole deformation parameter and $P_2(\cos\theta) = (3\cos^2\theta - 1)/2$is the Legendre polynomial. Volume conservation to second order requires:

$$R_0 = r_0 A^{1/3}\left(1 - \frac{\alpha_2^2}{4\pi}\frac{1}{5} + \cdots\right)$$

Step 2: Surface and Coulomb Energies Under Deformation

The surface area of the deformed shape, expanded to second order in $\alpha_2$:

$$S = 4\pi R_0^2\left(1 + \frac{2}{5}\alpha_2^2 + \cdots\right)$$

Therefore the surface energy increases:

$$E_S = E_S^{(0)}\left(1 + \frac{2}{5}\alpha_2^2\right), \quad E_S^{(0)} = a_S A^{2/3}$$

The Coulomb energy of a uniformly charged deformed body decreases:

$$E_C = E_C^{(0)}\left(1 - \frac{1}{5}\alpha_2^2\right), \quad E_C^{(0)} = a_C \frac{Z(Z-1)}{A^{1/3}} \approx a_C\frac{Z^2}{A^{1/3}}$$

Step 3: Deformation Energy and Fissility

The total change in energy upon deformation is:

$$\Delta E = (E_S + E_C) - (E_S^{(0)} + E_C^{(0)}) = \alpha_2^2\left(\frac{2}{5}E_S^{(0)} - \frac{1}{5}E_C^{(0)}\right)$$

The nucleus is unstable against deformation when $\Delta E < 0$:

$$\frac{2}{5}E_S^{(0)} < \frac{1}{5}E_C^{(0)} \implies E_C^{(0)} > 2E_S^{(0)}$$

Defining the fissility parameter:

$$x = \frac{E_C^{(0)}}{2E_S^{(0)}} = \frac{a_C Z^2 / A^{1/3}}{2a_S A^{2/3}} = \frac{a_C}{2a_S}\frac{Z^2}{A}$$

The critical value is $x = 1$, corresponding to:

$$\left(\frac{Z^2}{A}\right)_{\rm crit} = \frac{2a_S}{a_C} = \frac{2 \times 17.23}{0.70} \approx 49.2$$

For $x < 1$, a fission barrier exists (height depends on $x$). For $x \geq 1$, no barrier exists and the nucleus instantly fissions.

Fission Barrier Calculation

Barrier Height in the Liquid Drop Model

For a nucleus with fissility parameter $x < 1$, the fission barrier height in the liquid drop approximation is approximately:

$$B_f^{\rm LD} \approx E_S^{(0)} \cdot 0.38(1 - x)^3 \text{ (for } x > 0.67\text{)}$$

For specific nuclei:

$^{236}$U* (x = 0.74)

$B_f \approx 5.8$ MeV

$^{240}$Pu* (x = 0.76)

$B_f \approx 5.0$ MeV

$^{252}$Cf (x = 0.79)

$B_f \approx 4.0$ MeV

Double-Humped Fission Barrier

The liquid drop model predicts a single barrier, but shell effects create a double-humped barrier structure for actinide nuclei. At intermediate deformations ($\alpha_2 \approx 0.3\text{-}0.6$), shell effects produce a second minimum:

$$V_{\rm total}(\alpha) = V_{\rm LD}(\alpha) + \delta E_{\rm shell}(\alpha)$$

The structure of the double-humped barrier:

  • - First minimum: Ground state (spherical or slightly deformed)
  • - First barrier (A): Height $\sim 6$ MeV above ground state
  • - Second minimum: Shape isomeric state, typically 2-3 MeV above ground state
  • - Second barrier (B): Height $\sim 5$ MeV above ground state
  • - Scission point: Nucleus divides into two fragments

The second minimum gives rise to fission isomers -- metastable states with lifetimes of nanoseconds to milliseconds that decay preferentially by fission rather than gamma emission back to the ground state.

Mass and Energy Distributions of Fission Fragments

Asymmetric Mass Division

Thermal fission of actinides produces an asymmetric mass distribution with two peaks. For $^{235}$U(n,f), the light fragment peaks at $A_L \approx 95$ and the heavy fragment at $A_H \approx 140$:

$$A_L + A_H = A_{\rm CN} = 236, \quad \frac{A_H}{A_L} \approx 1.47$$

The asymmetry is driven by shell effects in the fragments, particularly the proximity of the heavy fragment to the doubly-magic $^{132}$Sn (Z=50, N=82). As the excitation energy increases, the mass distribution becomes more symmetric.

Kinetic Energy Distribution

The total kinetic energy (TKE) of the fragments is primarily Coulomb repulsion at the scission point:

$$\text{TKE} \approx \frac{Z_1 Z_2 e^2}{4\pi\epsilon_0 D} = \frac{Z_1 Z_2 \times 1.44}{r_0(A_1^{1/3} + A_2^{1/3}) + d} \text{ MeV}$$

where $D$ is the center-to-center distance at scission, typically$D \approx r_0(A_1^{1/3} + A_2^{1/3}) + 2$ fm. For $^{235}$U fission:

$$\langle\text{TKE}\rangle \approx \frac{36 \times 56 \times 1.44}{12 + 2} \approx 170 \text{ MeV}$$

The Viola-Seaborg systematics give an empirical formula:$\langle\text{TKE}\rangle = 0.1189\,Z^2/A^{1/3} + 7.3$ MeV, which works well across the actinide region.

Charge Distribution

For a given mass split, the charge distribution of fragments is approximately Gaussian, centered at the unchanged charge distribution (UCD) value:

$$Z_{\rm UCD} = \frac{Z_{\rm CN}}{A_{\rm CN}}\cdot A_{\rm frag}$$

The width is $\sigma_Z \approx 0.6$ charge units. The most probable charge is slightly below the UCD value due to the tendency toward equal neutron-to-proton ratios in each fragment (charge polarization).

Prompt and Delayed Neutrons

Prompt Neutrons

Fission fragments are born highly excited and neutron-rich. They de-excite by emitting prompt neutrons (within $\sim 10^{-14}$ s of scission):

$$\bar{\nu}_p(^{235}\text{U}) = 2.43, \quad \bar{\nu}_p(^{239}\text{Pu}) = 2.88, \quad \bar{\nu}_p(^{252}\text{Cf}) = 3.76$$

The prompt neutron energy spectrum follows the Watt distribution:

$$\chi(E) = C\,e^{-E/a}\sinh\sqrt{bE}$$

with $a \approx 0.988$ MeV and $b \approx 2.249$ MeV$^{-1}$ for $^{235}$U. The average prompt neutron energy is $\langle E \rangle \approx 2$ MeV.

Delayed Neutrons

A small fraction ($\beta_{\rm eff} \approx 0.0065$ for $^{235}$U) of fission neutrons are emitted following beta decay of certain fission products (delayed neutron precursors). These are crucial for reactor control:

$$\text{Example: } ^{87}\text{Br} \xrightarrow{\beta^-}_{55.6\text{s}} ^{87}\text{Kr}^* \xrightarrow{n} ^{86}\text{Kr}$$

The delayed neutron fraction is grouped into six time groups with half-lives ranging from 0.2 to 55 seconds. The effective multiplication factor including delayed neutrons is:

$$k_{\rm eff} = k_p(1 - \beta) + k_d\beta$$

Reactor control is possible only because the delayed neutron fraction provides a time scale of seconds (rather than microseconds) for power changes. Without delayed neutrons, the reactor period for a 0.1% excess reactivity would be $T \approx \ell/\delta k \sim 10^{-4}/10^{-3} = 0.1$ s (too fast to control). With delayed neutrons, $T \approx \beta/(\lambda\delta k) \sim 80$ s.

Spontaneous Fission

Spontaneous fission (SF) is the quantum tunneling of the entire nucleus through the fission barrier, analogous to alpha decay tunneling through the Coulomb barrier.

SF Half-Lives

The SF half-life depends exponentially on the barrier height and width. The action integral through the fission barrier is:

$$S = \frac{2}{\hbar}\int_{\alpha_1}^{\alpha_2}\sqrt{2B(\alpha)[V(\alpha) - E_0]}\,d\alpha$$

where $B(\alpha)$ is the collective inertia (mass parameter) and$V(\alpha)$ is the potential energy along the fission path. SF half-lives span an enormous range:

$^{238}$U

$t_{1/2}^{\rm SF} \approx 10^{16}$ y

$^{252}$Cf

$t_{1/2}^{\rm SF} \approx 85$ y

$^{256}$Fm

$t_{1/2}^{\rm SF} \approx 2.9$ h

Systematics

The SF half-life correlates strongly with $Z^2/A$:

$$\log_{10} t_{1/2}^{\rm SF} \approx a - b\frac{Z^2}{A} + \text{shell corrections}$$

SF becomes the dominant decay mode for very heavy nuclei ($Z \gtrsim 100$). Odd-A and odd-odd nuclei have SF half-lives typically $10^3\text{-}10^5$ times longer than neighboring even-even nuclei, due to the specialization energy (the unpaired nucleon blocks the fission path).

Ternary Fission

In a small fraction of fission events (~0.2-0.4%), a light charged particle is emitted in addition to the two main fragments. This is called ternary fission.

Ternary Particle Spectrum

The most common ternary particle is the alpha particle (~90% of ternary events), emitted roughly perpendicular to the fission axis:

  • - Alpha particles: ~0.2% of all fissions, average energy ~16 MeV, emitted at ~82$^\circ$ to the fission axis
  • - Tritons ($^3$H): ~0.01% of fissions, average energy ~8.5 MeV
  • - Other light ions: $^6$He, $^8$He, $^{6-9}$Li, $^{7-11}$Be at much lower rates

Physical Mechanism

Ternary particles are thought to be formed in the neck region between the two nascent fragments just before or at scission. The emission angle and energy distributions can be understood from a Coulomb trajectory calculation:

$$E_\alpha \approx \frac{Z_\alpha e^2}{4\pi\epsilon_0}\left(\frac{Z_L}{D_L} + \frac{Z_H}{D_H}\right)$$

where $D_L$ and $D_H$ are the distances from the alpha to the light and heavy fragments at scission. The perpendicular emission is a natural consequence of the Coulomb fields from the two fragments, which accelerate the ternary particle away from the fission axis.

Share:XRedditLinkedIn
Rate this chapter: