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.
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:
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 < 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:
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
PythonDouble-humped mass distribution for U-235 and Pu-239 thermal fission
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
FortranFour-factor formula k-effective calculation for various enrichment levels
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:
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:
Step 2: Surface and Coulomb Energies Under Deformation
The surface area of the deformed shape, expanded to second order in $\alpha_2$:
Therefore the surface energy increases:
The Coulomb energy of a uniformly charged deformed body decreases:
Step 3: Deformation Energy and Fissility
The total change in energy upon deformation is:
The nucleus is unstable against deformation when $\Delta E < 0$:
Defining the fissility parameter:
The critical value is $x = 1$, corresponding to:
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:
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:
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$:
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:
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:
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:
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):
The prompt neutron energy spectrum follows the Watt distribution:
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:
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:
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:
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$:
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:
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.