Nuclear Fission Mechanics, Energy Release & Product Yields

1. Fission Energetics and the Liquid Drop Analogy

The liquid drop model (LDM), formulated by Niels Bohr and John Archibald Wheeler (1939), models the heavy atomic nucleus as an incompressible, charged macroscopic liquid droplet with constant nuclear matter density $\rho_0 \approx 0.17\text{ nucleons/fm}^3$. The nuclear ground state represents an equilibrium configuration between two opposing forces:

• Short-range attractive nuclear force: Represented macroscopically as surface tension energy $E_s$, which acts to minimize nuclear surface area (favoring a spherical shape).
• Long-range repulsive Coulomb force: Represented as electrostatic self-energy $E_c$, which pushes protons apart and drives the nucleus toward elongated, non-spherical deformations.

2. Ellipsoidal Deformation and Fissility Parameter

Consider a spherical nucleus of radius $R_0 = r_0 A^{1/3}$ undergoing a small quadrupole deformation into a prolate spheroid with semi-major axis $a = R_0 (1 + \epsilon)$ and semi-minor axes $b = c = R_0 (1 - \frac{1}{2}\epsilon)$, preserving volume to second order in eccentricity $\epsilon$: $$R(\theta) = R_0 \left[ 1 + \alpha_2 P_2(\cos\theta) \right]$$ Expanding the surface and Coulomb energy terms in powers of the deformation coordinate $\alpha_2$: $$E_s(\alpha_2) = E_s^{(0)} \left( 1 + \frac{2}{5} \alpha_2^2 + \mathcal{O}(\alpha_2^3) \right)$$ $$E_c(\alpha_2) = E_c^{(0)} \left( 1 - \frac{1}{5} \alpha_2^2 + \mathcal{O}(\alpha_2^3) \right)$$ where the undeformed ground-state values from the semi-empirical mass formula (Weizsäcker) are: $$E_s^{(0)} = a_s A^{2/3} \approx 17.8 \times A^{2/3}\text{ MeV}, \qquad E_c^{(0)} = a_c \frac{Z^2}{A^{1/3}} \approx 0.71 \times \frac{Z^2}{A^{1/3}}\text{ MeV}$$ The change in net nuclear energy $\Delta E = \Delta E_s + \Delta E_c$ as a function of deformation $\alpha_2$ is: $$\Delta E(\alpha_2) = \alpha_2^2 \left( \frac{2}{5} E_s^{(0)} - \frac{1}{5} E_c^{(0)} \right) = \frac{1}{5} \alpha_2^2 \left( 2 E_s^{(0)} - E_c^{(0)} \right)$$

For the spherical shape to remain stable against spontaneous deformation, the coefficient of $\alpha_2^2$ must be positive ($\Delta E > 0$). When electrostatic repulsion overcomes surface tension, the sphere becomes unstable to infinitesimal perturbations ($\Delta E < 0$): $$2 E_s^{(0)} - E_c^{(0)} < 0 \implies \frac{E_c^{(0)}}{2 E_s^{(0)}} > 1$$ Substituting the Bethe-Weizsäcker coefficients: $$\frac{0.71 \times Z^2 / A^{1/3}}{2 \times 17.8 \times A^{2/3}} \approx \frac{Z^2 / A}{50.1} > 1 \implies \mathbf{\frac{Z^2}{A} \gtrsim 49\text{ to }50}$$ The dimensionless ratio $x \equiv \frac{E_c^{(0)}}{2 E_s^{(0)}} \approx \frac{Z^2/A}{(Z^2/A)_{\text{crit}}}$ is known as the Bohr-Wheeler fissility parameter. For heavy actinides such as ${}^{238}\text{U}$ ($Z^2/A \approx 35.6$) and ${}^{235}\text{U}$ ($Z^2/A \approx 36.0$), $x \approx 0.72$, meaning they are stable against spontaneous deformation at the ground state but possess a finite fission activation barrier.

3. The Fission Barrier and Saddle Point

As deformation increases beyond the harmonic quadratic regime, higher multipoles ($\alpha_3, \alpha_4$) develop necking, culminating in the saddle point—the maximum potential energy along the lowest-energy fission valley. The height of this peak above the ground state is the fission barrier $E_B$: $$E_B = E_{\text{saddle}} - E_{\text{ground}}$$ Representative fission barrier heights:

• ${}^{235}\text{U}$: $E_B \approx 5.75\text{ MeV}$
• ${}^{238}\text{U}$: $E_B \approx 6.25\text{ MeV}$
• ${}^{239}\text{Pu}$: $E_B \approx 5.50\text{ MeV}$
• ${}^{232}\text{Th}$: $E_B \approx 6.50\text{ MeV}$

When a thermal neutron ($E_n \approx 0.025\text{ eV}$) is captured by ${}^{235}\text{U}$, the binding energy of the added neutron ($S_n = 6.55\text{ MeV}$) exceeds the barrier ($S_n > E_B$), inducing prompt fission with zero threshold. In contrast, capture by ${}^{238}\text{U}$ yields only $S_n = 4.81\text{ MeV} < E_B = 6.25\text{ MeV}$, leaving a deficit $\Delta E \approx 1.44\text{ MeV}$ that must be supplied as incident neutron kinetic energy.

1. Fuel Classification and Enrichment Regimes

Nuclear fuels in fission reactors are categorized according to isotopic composition, physical phase, and chemical matrix:

• Natural Uranium (NatU): $0.7204\%\ {}^{235}\text{U}$, $99.274\%\ {}^{238}\text{U}$, $0.0055\%\ {}^{234}\text{U}$. Used directly in heavy-water moderated reactors (CANDU, PHWR) and carbon dioxide-cooled graphite reactors (Magnox).
• Low-Enriched Uranium (LEU): Uranium enriched to $3\%\text{ to }5\%\ {}^{235}\text{U}$ by mass. Standard commercial fuel for Light Water Reactors (PWR and BWR).
• High-Assay Low-Enriched Uranium (HALEU): Uranium enriched between $5\%\text{ and }20\%\ {}^{235}\text{U}$. Crucial for advanced Small Modular Reactors (SMRs) and Generation IV fast reactors, offering high power density and extended refueling cycles without crossing proliferation safeguards.
• Highly Enriched Uranium (HEU): Uranium enriched to $\ge 20\%\ {}^{235}\text{U}$ (weapons-grade $\ge 90\%$). Restricted to naval propulsion reactors, research reactors, and space fission systems.

2. Mixed Oxide (MOX) and Plutonium Recycling

Mixed Oxide fuel consists of depleted uranium dioxide ($\text{UO}_2$) blended with recycled reactor-grade plutonium dioxide ($\text{PuO}_2$), typically containing $5\%\text{ to }9\%\ \text{Pu}_{\text{fissile}}$ (${}^{239}\text{Pu} + {}^{241}\text{Pu}$): $$\text{MOX} = (1 - x)\,\text{UO}_2 + x\,\text{PuO}_2$$ Key neutron physical consequences of MOX fuel in thermal cores:

• Hardened Neutron Spectrum: Massive capture and fission resonances of ${}^{239}\text{Pu}$ ($0.296\text{ eV}$) and ${}^{240}\text{Pu}$ ($1.056\text{ eV}$) depress the thermal flux and harden the average energy spectrum.
• Reduced Delayed Neutron Fraction: $\beta_{\text{eff}}$ decreases because $\beta({}^{239}\text{Pu}) \approx 0.00215$ compared to $\beta({}^{235}\text{U}) \approx 0.00650$, requiring tighter control rod margins.
• Increased Control Rod Absorption Requirements: Stronger thermal absorption in fuel requires higher boron concentration or enriched boron carbide ($\text{B}_4\text{C}$) control rods.

3. The Thorium-Uranium Fuel Cycle

The thorium fuel cycle exploits fertile ${}^{232}\text{Th}$ to breed fissile ${}^{233}\text{U}$: $${}^{232}_{90}\text{Th} + n \longrightarrow {}^{233}_{90}\text{Th} \xrightarrow[\beta^-, \, 22.3\text{ min}]{} {}^{233}_{91}\text{Pa} \xrightarrow[\beta^-, \, 26.97\text{ days}]{} {}^{233}_{92}\text{U}$$ The major nuclear advantage is that ${}^{233}\text{U}$ exhibits $\eta > 2.25$ across the entire thermal neutron spectrum (superior to ${}^{235}\text{U}$ at $\eta = 2.07$ and ${}^{239}\text{Pu}$ at $\eta = 2.11$), making thermal breeding feasible (e.g. in Molten Salt Breeder Reactors - MSBR). In addition, thorium fuels produce orders of magnitude fewer long-lived transuranic actinides ($\text{Np, Pu, Am, Cm}$) because six successive neutron captures are needed to reach plutonium.

1. Mass Yield Curve: Bimodal vs Symmetric Cleavage

When low-energy (thermal) neutrons induce fission in ${}^{235}\text{U}$ or ${}^{239}\text{Pu}$, the nucleus splits overwhelmingly into two fragments of unequal mass. The percentage yield of fission fragments as a function of mass number $A$ is the fission product mass yield curve $Y(A)$, normalized such that: $$\sum_{A} Y(A) = 200\%$$ since each fission produces exactly two primary fission fragments.

For thermal fission of ${}^{235}\text{U}$:

• Light Fragment Peak: Centered at mass number $A_L \approx 95$ (peak yield $\sim 6.5\%$ at isotopes like ${}^{95}\text{Mo}, {}^{95}\text{Zr}, {}^{90}\text{Sr}$).
• Heavy Fragment Peak: Centered at mass number $A_H \approx 138\text{ to }140$ (peak yield $\sim 6.5\%$ at ${}^{137}\text{Cs}, {}^{135}\text{Xe}, {}^{140}\text{Ba}$).
• Symmetric Fission Valley: At $A \approx 117$, symmetric splitting is suppressed by a factor of over $600$ ($Y(117) \approx 0.01\%$).

2. Shell Effects and Energy-Dependent Peak-to-Valley Ratio

The asymmetric splitting is explained by single-particle nuclear shell model effects in the deformed nascent fragments:

• The heavy fragment favors magic numbers: $Z = 50$ (closed spherical proton shell) and $N = 82$ (closed spherical neutron shell, $A \approx 132$, e.g. ${}^{132}_{50}\text{Sn}$), which provides exceptional quantum stability.
• As incident neutron energy increases from thermal to fast ($E_n \sim 14\text{ MeV}$ in D-T fusion), excitation washes out shell effects; the central valley fills in, and the distribution shifts toward symmetric fission.

1. Microscopic Energy Partition Breakdown

When a ${}^{235}\text{U}$ nucleus captures a thermal neutron and fissions, approximately $200\text{ MeV}$ of total energy is liberated. The microscopic distribution of this energy among products is summarized below:

Energy Component	Energy (MeV)	Fraction (%)	Recovery Mechanism
Kinetic Energy of Fission Fragments	$168 \pm 5$	$84.0\%$	Thermalized locally within $\sim 10\ \mu\text{m}$ in fuel
Kinetic Energy of Prompt Neutrons ($\nu \approx 2.43$)	$5 \pm 0.5$	$2.5\%$	Thermalized via elastic scattering in moderator
Prompt Fission Gamma Rays ($\gamma_{\text{prompt}}$)	$7 \pm 1$	$3.5\%$	Absorbed in fuel, cladding, coolant, shield
Fission Product Beta Decay ($\beta^-$)	$8 \pm 1$	$4.0\%$	Delayed decay heat inside fuel matrix
Fission Product Delayed Gamma Rays ($\gamma_{\text{delayed}}$)	$7 \pm 1$	$3.5\%$	Delayed decay heat in core structure
Antineutrinos ($\bar{\nu}_e$)	$12 \pm 2$	$6.0\%$	Lost entirely (escape reactor vessel)
Radiative Capture Gamma Rays ($(n, \gamma)$ in non-fission)	$3\text{ to }10$	$\sim 2.5\text{ to }5\%$	Exothermic capture in structural cladding/water

2. Recoverable Energy per Fission ($Q_{\text{rec}}$)

The antineutrinos possess a negligible interaction cross section ($\sigma_\nu \sim 10^{-44}\text{ cm}^2$) and escape into outer space without depositing heat. Conversely, excess neutrons produced by fission are absorbed via radiative capture $(n, \gamma)$ in moderator, coolant, and structural cladding, releasing additional binding gamma energy ($\sim 3\text{ to }8\text{ MeV}$). Hence, the net recoverable thermal energy per fission event is: $$Q_{\text{rec}} = Q_{\text{total}} - E_{\bar{\nu}_e} + E_{(n,\gamma)} \approx 207\text{ MeV} - 12\text{ MeV} + 5\text{ MeV} \approx \mathbf{200 \pm 2\text{ MeV}}$$ In engineering calculations: $$1\text{ fission} \approx 200\text{ MeV} = 3.20435 \times 10^{-11}\text{ Joules} = 3.20435 \times 10^{-17}\text{ MW}\cdot\text{s}$$ $$3.12 \times 10^{10}\text{ fissions/second} \longleftrightarrow \mathbf{1\text{ Watt of thermal power}}$$

1. The Watt Fission Spectrum Formula

Neutrons liberated at the moment of scission (prompt fission neutrons) are emitted isotropically from rapidly moving, highly excited fission fragments. Transforming from the center-of-mass frame of the fragment to the laboratory frame results in the classical Watt distribution $\chi(E)$: $$\chi(E) = c \, e^{-E/a} \sinh\left(\sqrt{b E}\right)$$ For thermal neutron fission of Uranium-235 (${}^{235}\text{U}$), empirical parameters fitted to nuclear data (ENDF/B-VIII) are: $$a = 0.988\text{ MeV}, \qquad b = 2.249\text{ MeV}^{-1}$$ The normalization constant $c$ enforces total probability conservation: $$\int_{0}^{\infty} \chi(E) \, dE = 1 \implies c = \sqrt{\frac{4}{\pi a^3 b}} \, e^{-a b / 4}$$ An alternate, widely used analytical approximation developed by Cranberg is: $$\chi(E) = 0.453 \, e^{-1.036 E} \sinh\left(\sqrt{2.29 E}\right) \quad (E\text{ in MeV})$$ or the Maxwellian fission spectrum approximation: $$\chi_M(E) = \frac{2}{\sqrt{\pi} T^{3/2}} E^{1/2} e^{-E/T} \quad (T \approx 1.29\text{ MeV})$$

2. Mean and Most Probable Prompt Neutron Energy

The energy distribution spans eight orders of magnitude, with:

• Most Probable Energy ($E_p$): Found by solving $\frac{d\chi(E)}{dE} = 0$: $$E_p \approx \mathbf{0.73\text{ MeV}}$$
• Average (Mean) Kinetic Energy ($\bar{E}$): $$\bar{E} = \int_{0}^{\infty} E \, \chi(E) \, dE = \frac{3}{2} a + \frac{1}{4} a^2 b \approx \frac{3}{2}(0.988) + \frac{1}{4}(0.988)^2(2.249) \approx \mathbf{1.98\text{ to }2.00\text{ MeV}}$$

Thus, fission neutrons are born fast, with an average speed $\bar{v} = \sqrt{2\bar{E}/m_n} \approx 1.95 \times 10^7\text{ m/s} \approx 0.065\, c$ (6.5% the speed of light)!

1. Prompt Neutrons and Delayed Precursors

Over $99\%$ of fission neutrons are emitted within $\tau_{\text{prompt}} \sim 10^{-14}\text{ to }10^{-13}\text{ seconds}$ of scission—these are prompt neutrons. However, a tiny but vital fraction ($\beta \approx 0.65\%$ in ${}^{235}\text{U}$) are emitted with delays ranging from milliseconds to minutes: $$\text{Fission} \longrightarrow \text{Precursor Nuclide } ({}^{87}\text{Br}) \xrightarrow[\beta^-]{T_{1/2} = 55.6\text{ s}} \text{Emitter } ({}^{87}\text{Kr}^*) \xrightarrow[\text{prompt } n]{< 10^{-15}\text{ s}} {}^{86}\text{Kr} + n$$ The emission of the neutron itself is prompt from an excited daughter nucleus whose excitation energy exceeds the neutron separation energy ($E^* > S_n$), but the emission rate is governed strictly by the preceding beta decay half-life of the parent precursor nuclide.

2. The Six Delayed Neutron Precursor Groups

In reactor physics and point kinetics, the hundreds of distinct delayed neutron emitting fission products are grouped into six standard precursor groups according to decay constant $\lambda_i = \ln 2 / T_{1/2, i}$:

Group $i$	Representative Precursor	Half-Life $T_{1/2}$ (s)	Decay Const $\lambda_i\ (\text{s}^{-1})$	Rel. Abundance $\beta_i / \beta$	Average Energy (keV)
1	${}^{87}\text{Br}$	$55.72$	$0.0124$	$0.033$	$250$
2	${}^{137}\text{I}$	$22.72$	$0.0305$	$0.219$	$560$
3	${}^{138}\text{I}, {}^{89}\text{Br}$	$6.22$	$0.111$	$0.196$	$430$
4	${}^{93}\text{Rb}, {}^{139}\text{I}$	$2.30$	$0.301$	$0.395$	$620$
5	${}^{94}\text{Rb}, {}^{140}\text{I}$	$0.61$	$1.14$	$0.115$	$420$
6	${}^{95}\text{Rb}$	$0.23$	$3.01$	$0.042$	$510$

3. Physical Role in Reactor Control

Without delayed neutrons, the average generation time between fission generations would be dictated entirely by the prompt neutron lifetime $l \sim 10^{-4}\text{ s}$ in thermal reactors ($l \sim 10^{-7}\text{ s}$ in fast reactors). For a tiny reactivity insertion $\Delta k = 0.001$, reactor power would escalate as $e^{t / (l/\Delta k)} = e^{10 t}$, multiplying by $22{,}000$ every second—rendering mechanical control rods completely ineffective! Delayed neutrons stretch the effective generation time to: $$\bar{l}_{\text{eff}} = (1 - \beta) l + \sum_{i=1}^6 \beta_i \tau_i = (1 - \beta) l + \sum_{i=1}^6 \frac{\beta_i}{\lambda_i} \approx 0.1\text{ seconds}$$ This slows reactor power transients by three orders of magnitude, allowing mechanical control rods and electronic safety circuits ample time to maintain safe steady-state operation.

1. Relation Between Thermal Power and Fission Rate

If a reactor operates at steady thermal power $P$ (in Watts), the core-wide total fission rate $\dot{F}$ (fissions per second) is: $$\dot{F} = \frac{P}{Q_{\text{rec}}} = \frac{P\ [\text{W}]}{3.20435 \times 10^{-11}\ [\text{J/fission}]} \approx 3.1208 \times 10^{10} \times P\ [\text{W}]$$ For a typical $3000\text{ MW}_{\text{th}}$ commercial power plant: $$\dot{F} = 3.1208 \times 10^{10} \times (3.0 \times 10^9\text{ W}) \approx \mathbf{9.36 \times 10^{19}\text{ fissions/second}}$$

2. Fissile Fuel Consumption Rate

Each fission consumes one fissile nucleus. In addition, radiative capture without fission $(n, \gamma)$ consumes additional nuclei in ratio $\alpha \equiv \sigma_c / \sigma_f$: $$\dot{N}_{\text{consumed}} = \dot{F} (1 + \alpha) = \dot{F} \frac{\sigma_a}{\sigma_f}$$ For Uranium-235 at thermal energies, $\sigma_f \approx 585\text{ b}$, $\sigma_c \approx 99\text{ b}$, giving capture-to-fission ratio: $$\alpha \approx \frac{99}{585} \approx 0.169 \implies 1 + \alpha \approx 1.169$$ The mass consumption rate of ${}^{235}\text{U}$ per Megawatt-day ($\text{MWd}$) of thermal energy is: $$\Delta m_{\text{fission}} = \frac{1\text{ MWd} \times 86400\text{ s/d} \times 235.044\text{ g/mol}}{(3.20435 \times 10^{-11}\text{ J}) \times (6.02214 \times 10^{23}\text{ atoms/mol})} \approx \mathbf{1.05\text{ g } {}^{235}\text{U} / \text{MWd}}$$ Including radiative capture losses ($1 + \alpha$): $$\Delta m_{\text{total}} = 1.05 \times 1.169 \approx \mathbf{1.23\text{ g } {}^{235}\text{U} / \text{MWd}}$$ A $3000\text{ MW}_{\text{th}}$ plant consumes: $$\Delta m = 3000\text{ MW} \times 1.23\text{ g/MWd} \approx \mathbf{3.69\text{ kg of } {}^{235}\text{U} \text{ per day}}$$

3. Fuel Burnup Formalism

Core Burnup measures the total thermal energy extracted per unit initial mass of heavy metal fuel loaded into the reactor, conventionally expressed in Megawatt-days per Metric Ton of Heavy Metal ($\text{MWd/MTU}$ or $\text{GWd/MTU}$): $$B = \frac{\int_0^T P(t) \, dt}{M_{\text{HM}, 0}} = \frac{E_{\text{th}}}{M_{\text{HM}, 0}}$$ where $M_{\text{HM}, 0}$ is the initial mass of uranium (or uranium + plutonium) metal in metric tons ($1\text{ MT} = 1000\text{ kg}$). Commercial LWRs routinely achieve burnups of $45\text{ to }55\text{ GWd/MTU}$ ($45{,}000\text{ to }55{,}000\text{ MWd/MTU}$), corresponding to the fission of approximately $5\%$ of all initial heavy metal atoms.