Fundamentals of Nuclear Energy, Atom Densities & Fission Fuel Breeding

1. Scope of Nuclear Reactor Physics

Nuclear reactor physics is the branch of applied nuclear physics and engineering that governs the distribution, transport, slowing down, and multiplication of neutrons inside a nuclear reactor core. The central objective is to sustain and control a steady-state or time-dependent chain reaction, converting microscopic nuclear binding energy into macroscopic thermal power: $$\text{Nuclear Energy} \longrightarrow \text{Thermal Energy} \longrightarrow \text{Mechanical Turbine Work} \longrightarrow \text{Electrical Power}$$

At the subatomic level, atomic nuclei are held together by the residual strong nuclear force, a phenomenological interaction mediated primarily by virtual meson exchange between nucleons (protons and neutrons). Key physical characteristics include:

• Colossal Magnitude: Approximately $100$ to $1000$ times stronger than the electromagnetic Coulomb repulsion at distances $r \approx 1\text{ fm}$.
• Extremely Short Range: Effective only over nuclear dimensions ($r \sim 1\text{ to }2\text{ fm}$), decaying exponentially as $\sim \frac{e^{-\mu r}}{r}$ for $r > 2\text{ fm}$.
• Saturation Property: A nucleon interacts only with its immediate nearest neighbors, meaning the binding energy per nucleon ($B/A \approx 8\text{ MeV}$) is approximately constant across intermediate and heavy nuclei.
• Hard Repulsive Core: Becomes fiercely repulsive at $r < 0.4\text{ fm}$, preventing the collapse of atomic nuclei into point singularities.

2. Energetics of Fission vs Fusion

The curve of binding energy per nucleon ($B/A$ versus mass number $A$) exhibits a prominent maximum near iron-56 (${}^{56}\text{Fe}$, where $B/A \approx 8.79\text{ MeV/nucleon}$). Because nature favors higher binding energy per nucleon (deeper potential wells), energy can be liberated exothermically via two distinct nuclear pathways:

1. Nuclear Fusion (Light Nuclei, $A < 56$): Fusing two light nuclei (such as deuterium and tritium, ${}^2\text{H} + {}^3\text{H} \to {}^4\text{He} + n + 17.6\text{ MeV}$) moves the products up the steep initial incline toward ${}^4\text{He}$ ($B/A \approx 7.07\text{ MeV}$), yielding $\sim 3.5\text{ MeV/nucleon}$.
2. Nuclear Fission (Heavy Nuclei, $A > 230$): Splitting a massive actinide nucleus (${}^{235}\text{U}$) into two intermediate fragments ($A_1 \sim 95, A_2 \sim 140$) shifts the system from $B/A \approx 7.56\text{ MeV}$ up to $B/A \approx 8.5\text{ MeV}$: $$\Delta(B/A) \approx 8.5\text{ MeV} - 7.56\text{ MeV} \approx 0.94\text{ MeV/nucleon}$$ Multiplying by $A = 236$ nucleons: $$Q_{\text{fission}} \approx 236 \times 0.94\text{ MeV} \approx \mathbf{200\text{ MeV per fission event}}$$

This energy density is astronomical: the complete fission of $1\text{ kg}$ of ${}^{235}\text{U}$ releases $8.2 \times 10^{13}\text{ Joules}$ of energy—equivalent to burning approximately $2,500\text{ metric tons}$ of high-grade coal or $14{,}000\text{ barrels}$ of crude oil!

1. Fissile Isotopes

Fissile materials are nuclides capable of undergoing nuclear fission upon capturing a neutron of any kinetic energy, including thermal neutrons with zero or near-zero kinetic energy ($E \approx 0.0253\text{ eV}$): $$n_{\text{th}} + {}^{A}_Z\text{X} \longrightarrow [{}^{A+1}_Z\text{X}]^* \longrightarrow \text{Fission Fragments} + \nu \, n + Q$$ The three practical fissile fuels of modern nuclear technology are:

• Uranium-235 (${}^{235}_{92}\text{U}$): The only naturally occurring fissile isotope on Earth, comprising $0.7204\%$ of natural uranium.
• Plutonium-239 (${}^{239}_{94}\text{Pu}$): Artificially produced via neutron capture in fertile ${}^{238}\text{U}$.
• Uranium-233 (${}^{233}_{92}\text{U}$): Artificially produced via neutron capture in fertile ${}^{232}\text{Th}$.

2. Fissionable Isotopes

Fissionable materials are nuclides capable of undergoing fission when struck by neutrons, but whose compound nucleus requires a threshold kinetic energy (typically $E_n > 1.0\text{ to }1.5\text{ MeV}$) to overcome the fission activation barrier: $$n_{\text{fast}} (E_n > E_{\text{th}}) + {}^{238}_{92}\text{U} \longrightarrow [{}^{239}_{92}\text{U}]^* \longrightarrow \text{Fission}$$ Examples include even-$N$ actinides such as ${}^{238}_{92}\text{U}$, ${}^{232}_{90}\text{Th}$, and ${}^{240}_{94}\text{Pu}$. All fissile materials are fissionable, but not all fissionable materials are fissile.

3. Fertile Isotopes and Nuclear Transmutation

Fertile materials are nuclides that are not themselves fissile by thermal neutrons, but can be converted into fissile isotopes through neutron radiative capture $(n, \gamma)$ followed by subsequent beta-minus ($\beta^-$) decays:

1. The Uranium-Plutonium Cycle: $${}^{238}_{92}\text{U} + n \longrightarrow {}^{239}_{92}\text{U} \xrightarrow[\beta^-, \, 23.5\text{ min}]{} {}^{239}_{93}\text{Np} \xrightarrow[\beta^-, \, 2.356\text{ days}]{} {}^{239}_{94}\text{Pu} \quad (\text{Fissile!})$$
2. The Thorium-Uranium Cycle: $${}^{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} \quad (\text{Fissile!})$$

1. Conversion Ratio and Breeding Ratio

In any nuclear reactor containing fertile material (${}^{238}\text{U}$ or ${}^{232}\text{Th}$), fissile fuel is simultaneously consumed by fission and produced by transmutation. We quantify this balance through the Conversion Ratio ($CR$), defined as: $$CR \equiv \frac{\text{Average rate of production of new fissile nuclei}}{\text{Average rate of consumption (fission + capture) of fissile nuclei}} = \frac{\dot{N}_{\text{fissile, produced}}}{\dot{N}_{\text{fissile, consumed}}}$$ Classification based on $CR$:

• Burner / Converter Reactor ($CR < 1$): Produces fewer fissile nuclei than it consumes. Commercial Light Water Reactors (LWRs) operate with $CR \approx 0.55\text{ to }0.65$.
• Break-Even Reactor ($CR = 1$): Exactly replaces every fissile atom consumed.
• Breeder Reactor ($CR > 1$): Produces more fissile fuel than it consumes. When $CR > 1$, it is designated the Breeding Ratio ($BR \equiv CR$).

2. The Doubling Time ($T_d$)

The doubling time $T_d$ is the operating time required for a breeder reactor to generate enough excess fissile material to completely fuel an identical second reactor (including out-of-pile reprocessing losses): $$T_d \approx \frac{M_{\text{core}}}{G \cdot \dot{M}_{\text{fissile, consumed}}}$$ In terms of reactor electric power $P_e$, thermal efficiency $\eta_{\text{th}}$, and capacity factor $CF$: $$T_d = \frac{M_{\text{fissile, inventory}}}{(BR - 1) \cdot (1 + \alpha) \cdot \dot{N}_{\text{fiss}} \cdot m_{\text{fiss}}}$$ Fast Breeder Reactors (FBRs) utilizing liquid metal cooling (Sodium or Lead) achieve breeding ratios of $BR \approx 1.20\text{ to }1.35$, yielding practical doubling times of $10\text{ to }20\text{ years}$.

1. Mass Defect $\Delta m$

According to Einstein's mass-energy equivalence relation $E = m c^2$, the ground-state mass of any bound nucleus $M(A, Z)$ is strictly less than the combined rest mass of its constituent free protons and neutrons: $$\Delta m \equiv \left[ Z m_p + (A - Z) m_n \right] - M(A, Z)$$ where $m_p = 1.007276466\text{ u}$ is the proton mass, $m_n = 1.008664916\text{ u}$ is the neutron mass, and $1\text{ u} = 931.494\text{ MeV}/c^2$. In terms of neutral atomic masses $M_{\text{atomic}}(A, Z)$ and hydrogen atom mass $m({}^1\text{H}) = 1.007825\text{ u}$: $$\Delta m = \left[ Z m({}^1\text{H}) + (A - Z) m_n \right] - M_{\text{atomic}}(A, Z)$$

2. Total Binding Energy and Binding Energy per Nucleon

The total nuclear binding energy $B(A, Z)$ is the energy required to disassemble the bound nucleus into isolated, non-interacting nucleons: $$B(A, Z) = \Delta m \cdot c^2 = \left( \left[ Z m({}^1\text{H}) + (A - Z) m_n \right] - M_{\text{atomic}}(A, Z) \right) \times 931.494\text{ MeV}$$ The binding energy per nucleon: $$\frac{B}{A} = \frac{B(A, Z)}{A}$$

1. Neutron Separation Energy ($S_n$)

The neutron separation energy $S_n$ is the minimum energy required to remove a single neutron from a nucleus $(A, Z)$, ejecting it to infinity at rest. It is the nuclear equivalent of the atomic ionization potential: $${}^{A}_{Z}\text{X} + S_n \longrightarrow {}^{A-1}_{Z}\text{X} + n$$ Using binding energies or atomic masses: $$S_n = B(A, Z) - B(A-1, Z) = \left[ M(A-1, Z) + m_n - M(A, Z) \right] c^2$$ Similarly, the proton separation energy is: $$S_p = B(A, Z) - B(A-1, Z-1) = \left[ M(A-1, Z-1) + m({}^1\text{H}) - M(A, Z) \right] c^2$$

2. The Fission Barrier and Why ${}^{235}\text{U}$ is Fissile but ${}^{238}\text{U}$ is Not

When an incident neutron with kinetic energy $E_n$ is absorbed by target nucleus ${}^{A}_Z\text{X}$, the resulting compound nucleus $[{}^{A+1}_Z\text{X}]^*$ is formed in an excited state with excitation energy $E^*$: $$E^* = S_n({}^{A+1}_Z\text{X}) + E_{\text{cm}} \approx S_n({}^{A+1}_Z\text{X}) + E_n \left( \frac{A}{A + 1} \right)$$ For thermal neutrons ($E_n \approx 0$), the entire excitation energy is simply the separation energy: $E^* \approx S_n$. To induce fission, $E^*$ must exceed the critical fission barrier energy $E_{\text{crit}}$ required to deform the nucleus past the saddle point of the liquid drop potential:

The profound physical origin is the pairing term in the semi-empirical mass formula: when a neutron is added to odd-$N$ $^{235}\text{U}$, it forms a neutron pair in even-even $^{236}\text{U}$, releasing an extra $\sim 1.2\text{ MeV}$ of pairing energy ($\delta_{\text{pair}} \approx +12 A^{-1/2}\text{ MeV}$). This extra pairing energy pushes $E^*$ well above the critical fission barrier, making $^{235}\text{U}$ readily fissile by zero-energy thermal neutrons!

1. Thermalization and Thermodynamic Equilibrium

In a thermal reactor, fast fission neutrons undergo repeated elastic collisions with light moderator nuclei until they reach approximate thermodynamic equilibrium with the thermal motion of the moderator atoms at core temperature $T$. The velocity distribution of thermal neutrons is given by the Maxwell-Boltzmann distribution: $$n(v) dv = \frac{4 n_0}{\sqrt{\pi}} \left( \frac{m}{2 k_B T} \right)^{3/2} v^2 \exp\left( -\frac{m v^2}{2 k_B T} \right) dv$$ where $n_0$ is the total thermal neutron number density ($\text{neutrons/cm}^3$) and $k_B = 8.61733 \times 10^{-5}\text{ eV/K} = 1.38065 \times 10^{-23}\text{ J/K}$.

2. Characteristic Thermal Parameters

From the Maxwellian distribution, three characteristic velocities and energies arise:

1. Most Probable Speed ($v_0$): Located at the peak of $n(v)$ where $dn/dv = 0$: $$v_0 = \sqrt{\frac{2 k_B T}{m}}$$ The corresponding kinetic energy is: $$E_0 = \frac{1}{2} m v_0^2 = k_B T$$ At the international standard reference temperature $T_0 = 293.61\text{ K}$ ($20.46^{\circ}\text{C}$): $$E_0 = (8.61733 \times 10^{-5}\text{ eV/K})(293.61\text{ K}) = \mathbf{0.0253\text{ eV}}$$ $$v_0 = \sqrt{\frac{2(0.0253 \times 1.6022 \times 10^{-19}\text{ J})}{1.67493 \times 10^{-27}\text{ kg}}} = \mathbf{2200\text{ m/s}}$$ This defines the universally tabulated $2200\text{ m/s}$ standard cross-section reference state.
2. Average Speed ($\bar{v}$): $$\bar{v} = \int_0^\infty v \frac{n(v)}{n_0} dv = \sqrt{\frac{8 k_B T}{\pi m}} = \frac{2}{\sqrt{\pi}} v_0 \approx 1.128 v_0$$
3. Root-Mean-Square Speed ($v_{\text{rms}}$): $$v_{\text{rms}} = \sqrt{\langle v^2 \rangle} = \sqrt{\frac{3 k_B T}{m}} = \sqrt{\frac{3}{2}} v_0 \approx 1.225 v_0$$ The average kinetic energy of the thermal neutron distribution is: $$\bar{E} = \frac{1}{2} m \langle v^2 \rangle = \frac{3}{2} k_B T = 1.5 E_0$$

1. Atom Density of a Pure Element

In reactor physics, reaction rates are evaluated using the macroscopic cross section $\Sigma = N \sigma$, where $N$ is the atom density (number of target nuclei per unit volume, typically expressed in $\text{atoms/cm}^3$ or $\text{atoms/b}\cdot\text{cm}$, where $1\text{ b}\cdot\text{cm} = 10^{-24}\text{ cm}^3$). For a pure element with mass density $\rho$ ($\text{g/cm}^3$) and atomic weight $M$ ($\text{g/mol}$): $$N = \frac{\rho \cdot N_A}{M}$$ where $N_A = 6.02214076 \times 10^{23}\text{ atoms/mol}$ is Avogadro's number.

2. Isotopic Atom Densities in an Enriched Mixture

For an element composed of multiple isotopes $i$, if $a_i$ denotes the atomic fraction ($\sum a_i = 1$), the atomic weight of the mixture is $M = \sum a_i M_i$, and the individual isotopic density is: $$N_i = a_i N = a_i \frac{\rho N_A}{M}$$ If the enrichment is given as a weight (mass) fraction $w_i \equiv m_i / m_{\text{total}}$ ($\sum w_i = 1$): $$N_i = \frac{w_i \cdot \rho \cdot N_A}{M_i}$$ The corresponding atomic fraction $a_i$ can be converted via: $$a_i = \frac{w_i / M_i}{\sum_j (w_j / M_j)}$$

3. Atom Densities in Chemical Compounds ($\text{A}_x\text{B}_y$)

For a homogeneous chemical compound with chemical formula $\text{A}_x\text{B}_y$, molecular weight $M_{\text{mol}} = x M_A + y M_B$, and mass density $\rho_{\text{comp}}$: The molecular density $N_{\text{mol}}$ is: $$N_{\text{mol}} = \frac{\rho_{\text{comp}} \cdot N_A}{M_{\text{mol}}}$$ The constituent elemental atom densities are: $$N_A = x \cdot N_{\text{mol}} = x \frac{\rho_{\text{comp}} N_A}{M_{\text{mol}}}, \qquad N_B = y \cdot N_{\text{mol}} = y \frac{\rho_{\text{comp}} N_A}{M_{\text{mol}}}$$ For example, in water ($\text{H}_2\text{O}$, $\rho = 1.0\text{ g/cm}^3$, $M_{\text{mol}} = 18.015\text{ g/mol}$): $$N_{\text{mol}} = \frac{1.0 \times 6.022 \times 10^{23}}{18.015} \approx 3.343 \times 10^{22}\text{ molecules/cm}^3$$ $$N_H = 2 N_{\text{mol}} \approx 6.686 \times 10^{22}\text{ atoms/cm}^3, \qquad N_O = N_{\text{mol}} \approx 3.343 \times 10^{22}\text{ atoms/cm}^3$$