In short

Nuclear fission is the process in which a heavy nucleus — almost always uranium-235 or plutonium-239 — absorbs a slow neutron, becomes grossly unstable, and splits into two lighter nuclei, two or three free neutrons, and a burst of energy. A single fission event of ^{235}\text{U} releases about 200 MeV — ten million times the energy released by burning one molecule of petrol.

The 200 MeV is bookkeeping on the binding-energy-per-nucleon curve. Uranium sits at B/A \approx 7.6 MeV per nucleon. The fission fragments (barium, krypton, caesium, strontium) sit near the peak at B/A \approx 8.5 MeV per nucleon. The difference, summed over 236 nucleons, is where the energy comes from.

The free neutrons make a chain reaction possible: each fission releases on average \nu \approx 2.4 neutrons, and if even one of them goes on to cause another fission, the reaction sustains itself. Define the multiplication factor k as the average number of next-generation fissions per current fission. Three regimes:

  • k < 1 (subcritical) — the reaction dies out.
  • k = 1 (critical) — steady, self-sustaining power.
  • k > 1 (supercritical) — exponential growth, a bomb.

A nuclear reactor runs at k = 1 using three ingredients: a moderator (heavy water at Tarapur and Kudankulam's predecessors, graphite elsewhere) that slows fast neutrons down so ^{235}\text{U} can absorb them; control rods (cadmium or boron) that soak up excess neutrons to hold k at exactly 1; and a coolant that carries the heat out to a steam turbine. India's PHWR design at Kaiga, Kakrapar and Rajasthan uses natural uranium and heavy water — a choice with a long history we will come back to.

At 2 a.m. on the morning of 28 October 1969, a control-rod operator at Tarapur Atomic Power Station near Mumbai lifted a bundle of cadmium rods a few centimetres out of a tank of boiling water the size of a cricket pitch. Inside that tank sat 132 tonnes of uranium dioxide fuel, stacked in Zircaloy tubes. The moment the rods cleared the fuel, something changed in the water. Neutrons that had been getting absorbed by the cadmium were now free to wander into uranium nuclei. Each absorption had a small probability of splitting a ^{235}\text{U} nucleus; each split released two-ish more neutrons; each of those had a small probability of splitting another ^{235}\text{U}. Within seconds, the number of fissions per second climbed by a factor of a thousand, then a million, then held steady at 4\times 10^{19} per second. The water started to boil. Steam shot up through the turbines. 160 megawatts of electrical power flowed out onto the Maharashtra grid.

That is what nuclear fission looks like, staged. A process that Enrico Fermi first demonstrated in a Chicago squash court in 1942, scaled up and turned into a power plant. To understand why the operator's action worked — and why a control rod has to move in centimetres rather than metres, and why the moderator has to be heavy water rather than ordinary water, and why uranium-235 is special but uranium-238 is not — you need three things: the binding energy curve, the cross-section of a nucleus for absorbing a neutron, and a single dimensionless number called k. This article builds all three.

The picture — a nucleus shaped like a waterdrop

Hold the uranium-235 nucleus in your head as a liquid drop. 92 protons, 143 neutrons, 236 when you add the captured neutron — all packed into a sphere about 7.4 \times 10^{-15} m across. Inside that sphere, two forces compete.

The nuclear force is attractive, short-ranged (it fades to nothing beyond about one nucleon diameter), and acts equally between protons and neutrons. It is what makes the drop want to stay spherical — surface tension, in the liquid analogy. Every nucleon on the surface of the drop has fewer neighbours than one in the interior, so surface nucleons cost energy, and the drop minimises energy by minimising surface area.

The electrostatic (Coulomb) force is repulsive, long-ranged, and acts only between protons. Every pair of the 92 protons pushes on every other pair, across the full diameter of the nucleus. Unlike the nuclear force, Coulomb repulsion does not fade with nucleon count — it grows as Z^2. For a nucleus as large as uranium, the Coulomb energy is very nearly equal to the surface-tension energy. The drop is on the verge of being unstable.

Now feed it one more neutron. The new nucleus is uranium-236 in an excited state — it has 6.5 MeV of extra energy, delivered as the neutron settled into its bound orbital. That energy has to go somewhere. In most nuclei, it would come off as a gamma ray and the nucleus would settle back to its ground state. But in ^{236}\text{U}, 6.5 MeV is almost exactly enough to deform the drop past the point of no return — past the fission barrier, which for uranium sits at about 6 MeV. The drop elongates into a peanut shape, then a dumbbell, then two separate drops joined by a thinning neck. The neck snaps. Coulomb repulsion, no longer held back by surface tension across the narrow waist, flings the two fragments apart at roughly 7% of the speed of light.

Liquid-drop model of fission: five stages from absorption to scissionFive panels showing uranium-235 absorbing a neutron, then deforming through oblate, dumbbell, and necked shapes, finally scissioning into two fragment nuclei that fly apart. ²³⁵U neutron absorbed ²³⁶U* excited, 6.5 MeV above ground deforms, crosses barrier neck thins Ba Kr fragments + 3n + 200 MeV
The liquid-drop model of fission. A slow neutron (small dark dot) is absorbed by $^{235}\text{U}$, forming excited $^{236}\text{U}^*$. Surface tension tries to keep the drop spherical; Coulomb repulsion tries to tear it apart. With 6.5 MeV of excitation, the drop deforms past the fission barrier, elongates, necks, and scissions into two fragments plus a few neutrons.

The liquid-drop picture was proposed by Bohr and Wheeler in 1939 — months after Lise Meitner and Otto Frisch interpreted Hahn and Strassmann's chemistry experiments as "splitting the uranium atom". It gets the broad physics right: fission is energetically favourable, but requires a small activation push to cross the barrier. Modern nuclear theory adds shell-model corrections (the fragments prefer certain magic-number configurations, making the mass distribution asymmetric), but for the energy balance, the liquid drop is all you need.

Where the 200 MeV comes from — read it off the binding curve

Every fission event of ^{235}\text{U} is different — the fragments come out with a distribution of masses, not fixed identities. But a representative reaction is:

^{235}_{92}\text{U} + n \longrightarrow ^{141}_{56}\text{Ba} + ^{92}_{36}\text{Kr} + 3n

Conservation of nucleon number: 235 + 1 = 141 + 92 + 3 \cdot 1 = 236. Conservation of charge: 92 + 0 = 56 + 36 + 0 = 92. Both sides balance.

Now compute the energy released using the binding-energy-per-nucleon curve. The relevant numbers, read off the standard curve:

Nucleus A B/A (MeV) Total B (MeV)
^{235}\text{U} 235 7.59 1784
^{141}\text{Ba} 141 8.33 1175
^{92}\text{Kr} 92 8.67 798

Why: B is the energy that would be needed to rip the nucleus apart into free nucleons. A higher B/A means the nucleus is more tightly bound, i.e., further down the energy hill. The free neutrons have B = 0 — a free neutron is not bound to anything.

Step 1. Compute the total binding energy before fission.

B_{\text{before}} = B(^{235}\text{U}) + B(n) = 1784 + 0 = 1784 \text{ MeV}

Why: only ^{235}\text{U} is bound. The incoming neutron is free. Its 0.025 eV of thermal kinetic energy is negligible compared to the MeV scale of nuclear binding.

Step 2. Compute the total binding energy after fission.

B_{\text{after}} = B(^{141}\text{Ba}) + B(^{92}\text{Kr}) + 3\cdot B(n) = 1175 + 798 + 0 = 1973 \text{ MeV}

Why: barium and krypton are both bound. The three emitted neutrons are free, so they contribute 0.

Step 3. The energy released is the increase in binding energy.

Q = B_{\text{after}} - B_{\text{before}} = 1973 - 1784 = 189 \text{ MeV}

Why: if the fragments are more tightly bound than the parent, the system has fallen to a lower energy state. The energy difference has to come out somewhere — as kinetic energy of the fragments, neutrons, and gamma photons. The exact number varies from 180 to 210 MeV depending on which fragments are produced, so "200 MeV" is the textbook average.

Result: Each fission event liberates about 200 MeV.

Where does that 200 MeV end up? Roughly:

Channel Energy (MeV)
Kinetic energy of the two heavy fragments 167
Kinetic energy of prompt fission neutrons (2–3 of them) 5
Prompt gamma rays 7
Beta particles from fragment decays 8
Antineutrinos from fragment decays 12
Delayed gamma rays from fragment decays 6

Everything except the antineutrinos stays inside the reactor fuel and moderator, gets converted to heat, and is what you eventually collect as steam-turbine power. The antineutrinos fly out through the reactor walls, through the Earth, through the Milky Way, and are lost. Kamiokande and SNOLAB (and, in India, the proposed India-based Neutrino Observatory at Theni) exist to catch a tiny fraction of them.

One fission = 200 MeV = 3.2 \times 10^{-11} J. One kilogram of ^{235}\text{U} contains N_A/0.235 \approx 2.56 \times 10^{24} atoms. If every one of them fissions, that is 8.2 \times 10^{13} J — the energy released by burning about 2,700 tonnes of coal. This is why the fuel loadings of reactors are so small compared to coal plants. Tarapur's 132-tonne fuel charge, replaced every few years, does the work of a coal mountain.

Why ^{235}\text{U} and not ^{238}\text{U} — the slow-neutron puzzle

Natural uranium is 99.27% ^{238}\text{U} and only 0.72% ^{235}\text{U}. Only the rare isotope fissions easily. To understand why, compare the fission barriers.

For ^{235}\text{U}, the barrier is \sim 5.8 MeV. Absorbing a thermal neutron delivers \sim 6.5 MeV of excitation. Excitation exceeds barrier by 0.7 MeV — fission is prompt and certain.

For ^{238}\text{U}, the barrier is \sim 7.0 MeV. Absorbing a thermal neutron delivers only \sim 4.8 MeV of excitation. Excitation is 2.2 MeV short — the drop sits in its excited state and relaxes by gamma emission, not fission.

Why the asymmetry: ^{235}\text{U} has an odd number of neutrons (143). When it absorbs one, the new neutron pairs up with the lone unpaired neutron, releasing the nucleon pairing energy of \sim 1.2 MeV — an extra energy bonus that is exactly what pushes the excitation over the barrier. ^{238}\text{U} already has paired neutrons (146), so the incoming neutron finds no partner and no pairing bonus. Nuclear pairing — the same physics that stabilises even-even nuclei like ^{4}\text{He} and ^{16}\text{O} — is ultimately why ^{235}\text{U} is the fuel of the nuclear age.

^{238}\text{U} can still fission if you hit it hard — with a "fast" neutron of a few MeV of kinetic energy. That kinetic energy, added to the excitation, does clear the barrier. This is why plutonium-fuelled fast breeder reactors (India's prototype FBR at Kalpakkam, Tamil Nadu) can burn ^{238}\text{U} indirectly: fast neutrons convert ^{238}\text{U} into fissile ^{239}\text{Pu}, which then fissions readily.

The bottom line: thermal (slow) neutrons fission ^{235}\text{U} and ^{239}\text{Pu}. Fast neutrons fission ^{238}\text{U} and ^{232}\text{Th} — but the cross-section is a hundred times smaller, so you need a dense neutron flux to make it work.

Chain reaction — the factor k

One fission creates on average \nu = 2.43 prompt neutrons (for thermal ^{235}\text{U} fission). Call the current generation n fissions per second. Most neutrons do not cause another fission — some escape the fuel, some get absorbed by ^{238}\text{U} without fissioning, some get swallowed by control rods or structural materials. Let k be the average number of next-generation fissions produced by each current fission.

k = \nu \cdot P_{\text{fission}}

where P_{\text{fission}} is the probability that a given prompt neutron survives long enough to cause another fission. Then the number of fissions per second grows as:

n_{g+1} = k \cdot n_g

After g generations, starting from n_0:

n_g = n_0 \cdot k^g

Why: each generation multiplies the population by k. This is the same algebra as compound interest, or bacterial growth, or any other process where the rate is proportional to the amount.

Since neutron generations in a thermal reactor take about \tau \approx 0.1 millisecond on average — time to thermalise in the moderator, diffuse to a fuel nucleus, and get absorbed — the time evolution is:

n(t) = n_0 \cdot k^{t/\tau} = n_0 \exp\!\left(\frac{t \ln k}{\tau}\right)

Three regimes:

Interactive: multiplication factor k controls the chain reactionTwo curves of fissions per second vs generation number, one linear-scale and one log-scale. Drag the horizontal marker to change k between 0.9 (subcritical) and 1.1 (supercritical). At k=1 the reaction is steady; below 1 it dies; above 1 it explodes. generation g n(g) / n(0) 0 5 10 15 20 0.01 0.1 1 10 100 k: 0.90 1.00 1.10 drag the red dot to change k
Drag the red dot to change $k$ between 0.9 and 1.1. The red curve shows $\log_{10}(n/n_0)$ vs generation. At $k = 1$ the curve is flat — steady power. At $k = 0.99$, after 20 generations you are down to $0.82$ of the initial — the reaction dies. At $k = 1.01$, after 20 generations you are up to $1.22$; continue for 1000 generations (0.1 s in a real reactor) and the factor becomes $2\times 10^4$. Reactor control is the art of keeping $k = 1$ to within 0.01%.

A reactor is not trying to grow its power — it is trying to hold it steady. The operator's job, when starting up Tarapur that October morning, was to inch k from 0.99 to 1.000001 — just barely supercritical, so the power rose gently over seconds to the target level, then back to exactly 1 to hold it there. If the control rods slip and k reaches 1.01, the power doubles every \tau / \ln(1.01) \approx 10 ms. If k reaches 1.1, it doubles every 1 ms. No control system can catch that. The reactor would melt in seconds.

This is the margin that nuclear engineers live in. Why it has been practical at all is the subject of the next section.

Delayed neutrons — the grace that makes reactors controllable

Not all fission neutrons are prompt. A small fraction, called \beta \approx 0.0065 for ^{235}\text{U}, are delayed — they are emitted not at the moment of scission but seconds later, from the beta decay of certain neutron-rich fragment nuclei (notably ^{87}\text{Br}^{87}\text{Kr}^*^{86}\text{Kr} + n, with a 56-second half-life).

This 0.65\% sounds small, but it completely transforms the control problem. The effective generation time is no longer 0.1 ms (prompt neutrons only) but roughly

\tau_{\text{eff}} = (1 - \beta)\tau_{\text{prompt}} + \beta \cdot \tau_{\text{delay}} \approx 0.9935 \times 10^{-4} + 0.0065 \times 10 \approx 0.065 \text{ s}

Why: the weighted average is dominated by the delayed contribution. Even though only 0.65% of neutrons are delayed, their 10-second mean delay makes them set the pace. Provided k stays within 1 \pm \beta, a reactor is delayed-critical — the effective time constant for power changes is seconds, not milliseconds, and mechanical control rods can comfortably catch it.

If k > 1 + \beta, the reactor becomes prompt critical — it achieves criticality on prompt neutrons alone, the delayed contribution becomes irrelevant, and the time constant crashes back to 0.1 ms. This is the runaway regime. Every reactor safety system is designed around one rule: never cross the prompt-critical boundary. The Chernobyl accident in 1986 was, in essence, a sudden transition to prompt critical. The Indian PHWR design has multiple independent shutdown systems — a fast-dropping shutoff rod bank and an emergency liquid-poison injection (gadolinium nitrate into the moderator) — specifically to keep this from happening.

Critical mass — geometry enters the story

A lump of pure ^{235}\text{U} at room temperature is not automatically supercritical. Some of the neutrons produced by a fission will escape the lump before they get absorbed, because neutrons travel a mean-free-path of a few centimetres through uranium before being absorbed. If the lump is small, most of the path is near the surface, so most neutrons escape. If the lump is large, most of the path is deep inside, so most neutrons get absorbed in fuel.

Balance these. Call \Sigma_f the macroscopic fission cross-section (fissions per neutron per cm travelled) and \Sigma_a the total absorption cross-section (fission + capture). The neutron-production-to-absorption ratio in an infinite lump is

k_\infty = \nu \cdot \frac{\Sigma_f}{\Sigma_a}

For pure ^{235}\text{U}, k_\infty \approx 2.1 — comfortably above 1. In a finite lump of radius R, some neutrons leak out. A leakage factor L(R) reduces the multiplication:

k_{\text{eff}}(R) = k_\infty \cdot L(R) = k_\infty \cdot \frac{1}{1 + (B_g / B_m)^2}

where B_g \sim \pi/R is a geometric buckling that measures how curved the neutron-density distribution is (tighter for smaller R) and B_m is a material buckling that depends on k_\infty and the diffusion properties of the fuel. Setting k_{\text{eff}} = 1 and solving for R gives the critical radius.

For a bare sphere of pure ^{235}\text{U} metal, the critical radius is about 8.5 cm — a grapefruit. The corresponding mass is about 52 kg. Surround the sphere with a neutron reflector (beryllium, tungsten, or natural uranium), and neutrons that would have escaped get bounced back. Critical mass drops to about 15 kg — a softball.

Why this matters: critical mass is the hard engineering number. For ^{239}\text{Pu}, the bare critical mass is only 10 kg — which is why plutonium, not uranium, became the fuel of choice for compact bomb designs. For power reactors, which operate with dilute fuel (natural uranium is only 0.72% ^{235}\text{U}), the "critical mass" concept generalises to a critical core volume containing tonnes of fuel. Tarapur's core is a cylinder 4 m across and 3 m tall, holding 132 tonnes of \text{UO}_2.

Anatomy of a reactor — moderator, control rods, coolant

A power reactor has to solve three problems at once:

  1. Slow the neutrons down. Prompt fission neutrons come out at \sim 2 MeV — too fast to be efficiently absorbed by ^{235}\text{U} (the thermal capture cross-section is 600 barns; the fast cross-section is 1 barn). You need to thermalise them to 0.025 eV, a factor of 10^8 reduction in kinetic energy. This is the job of the moderator.
  2. Hold k at exactly 1. If k drifts high, the power grows; if low, it dies. You need a knob. This is the job of the control rods.
  3. Carry the heat out. Each fission deposits 180 MeV of heat in the fuel. A 500-MW reactor is releasing 500 million joules per second of heat in a few cubic metres of fuel. Unless you extract it continuously, the fuel melts in seconds. This is the job of the coolant.
Schematic of a pressurised heavy-water reactor (PHWR)Cross-section of a PHWR showing the calandria with fuel channels, heavy-water moderator, control rods entering from above, primary coolant loop going to a steam generator, and secondary loop going to turbine and condenser. calandria (heavy-water moderator) fuel channels (natural UO₂) control rods (cadmium) hot D₂O ~300°C cooled D₂O steam generator steam turbine water return condenser sea water (Arabian / Bay of Bengal) primary loop: D₂O coolant (pressurised, 100 atm) secondary loop: ordinary water → steam
Cross-section of the Indian PHWR design (Kaiga, Kakrapar, Rajasthan). Natural-uranium fuel bundles sit in pressure tubes immersed in heavy-water moderator. Pressurised heavy-water coolant (red) carries heat to a steam generator, where ordinary water (dark) boils into steam to drive the turbine. Cadmium control rods dangle from above; they drop by gravity if power fails, instantly shutting the reactor down.

Moderator. The physics is elastic scattering. A neutron collides with a light nucleus and transfers kinetic energy; after about 20 collisions with hydrogen, or 120 with carbon, the neutron has lost enough energy to be thermal. But hydrogen — abundant and cheap in ordinary water (H₂O) — has a problem: it also absorbs neutrons, forming deuterium. In natural uranium (0.72% ^{235}\text{U}), this absorption is enough to push k_\infty below 1. You cannot build a natural-uranium reactor with ordinary water.

Two solutions: either enrich the uranium to 3–5% ^{235}\text{U} (the Tarapur approach, following the American BWR design imported in 1969) and use ordinary water; or keep natural uranium but use heavy water (D₂O), where the deuterium nucleus absorbs 600 times less efficiently than hydrogen. The price of heavy water is the cost of separating it from ordinary water — about 1 part in 6400 naturally, taking industrial distillation columns to concentrate.

India went heavy-water. The reason was geopolitics: after the 1974 Pokhran test, the Nuclear Suppliers Group embargoed enriched uranium to India. Enrichment technology became proliferation-controlled. But natural uranium and heavy water — the latter producible domestically — were not. The Bhabha Atomic Research Centre and the Department of Atomic Energy built a series of PHWRs at Kalpakkam, Narora, Kakrapar, Kaiga, and Rawatbhata, all running on natural-uranium fuel that India could mine at Jaduguda (Jharkhand) and enrich to D₂O at Vadodara (Gujarat). Kudankulam, built later with Russian cooperation, returned to enriched uranium. These two lineages — heavy-water and light-water — are the spine of the Indian nuclear programme.

Control rods. Cadmium or boron carbide, both of which have enormous neutron-absorption cross-sections (2450 barns for cadmium, 3840 barns for ^{10}\text{B} at thermal energies). Inserting a rod deeper into the core increases \Sigma_a, decreases k, and shuts the reaction down. Withdrawing opens k up. For fine control, a few rods move continuously; for shutdown, a separate bank drops fully in within seconds. The Indian PHWRs also have a secondary shutdown system: emergency injection of gadolinium nitrate poison directly into the moderator tank, guaranteeing shutdown even if every rod fails to drop.

Coolant. In a PHWR, heavy water at 100 atm pressure flows through the fuel channels at ~300 °C, never boiling. It carries heat to a steam generator where ordinary-water steam is raised. In a BWR (Tarapur), the coolant is ordinary water that is allowed to boil inside the reactor core itself — simpler, but the steam contains trace tritium. In India's fast breeder at Kalpakkam, the coolant is liquid sodium — needed because water would moderate the fast neutrons the breeder relies on.

Worked examples

Example 1: Power output of a reactor core

The Kaiga Atomic Power Station's Unit 4 is a 220 MW electrical PHWR. Thermal efficiency is about 30%. Estimate (a) the thermal power, (b) the number of fissions per second, and (c) the daily consumption of ^{235}\text{U} (assuming the fuel is 0.72% natural uranium and every ^{235}\text{U} atom that fissions contributes).

Energy flow from fission to grid at Kaiga Unit 4Energy flow diagram: fission events release 200 MeV each, thermal power 733 MW, steam turbine converts 30% to 220 MW electrical, rest is waste heat. fissions 2.3×10¹⁹ /s 200 MeV each thermal heat 733 MW η=30% electrical 220 MW to grid waste heat 513 MW → sea water
Energy flow at Kaiga Unit 4: 200 MeV per fission, multiplied up to 733 MW thermal, of which 30% becomes electricity.

Step 1. Find the thermal power.

P_{\text{thermal}} = \frac{P_{\text{electrical}}}{\eta} = \frac{220 \text{ MW}}{0.30} \approx 733 \text{ MW}

Why: thermal efficiency is the ratio of electrical output to heat input. A steam-turbine plant at ~300 °C hot side and ~30 °C cold side is Carnot-limited to about 47%; real plants achieve 30–35%. The rest is waste heat that the coolant dumps into the Arabian Sea (for Kaiga, which sits on the Karnataka coast).

Step 2. Convert thermal power to fissions per second.

Energy per fission = 200 \text{ MeV} = 200 \times 10^6 \times 1.602 \times 10^{-19} \text{ J} = 3.20 \times 10^{-11} \text{ J}.

\text{Fissions per second} = \frac{P_{\text{thermal}}}{E_{\text{per fission}}} = \frac{7.33 \times 10^8 \text{ W}}{3.20 \times 10^{-11} \text{ J}} \approx 2.3 \times 10^{19} \text{ per second}

Why: 1 W = 1 J/s, and each fission puts 3.2 × 10⁻¹¹ J into the fuel. Divide.

Step 3. Fissions per day and mass of ^{235}\text{U} consumed.

N_{\text{daily}} = 2.3 \times 10^{19} \times 86400 \approx 2.0 \times 10^{24} \text{ fissions/day}

One mole of ^{235}\text{U} is 0.235 kg. 2.0 \times 10^{24} atoms is 2.0\times 10^{24}/6.022\times 10^{23} = 3.32 mol.

m_{^{235}\text{U burned per day}} = 3.32 \text{ mol} \times 0.235 \text{ kg/mol} \approx 0.78 \text{ kg/day}

Why: Avogadro's number converts atoms to moles, and moles to mass. About three-quarters of a kilogram of ^{235}\text{U} per day generates 220 MW. For comparison, the Raichur coal plant in the same state burns about 6,000 tonnes of coal per day to produce the same power. The factor of about 8 million captures the whole case for nuclear energy.

Step 4. Natural-uranium fuel consumption.

Natural uranium is 0.72% ^{235}\text{U}, so the mass of natural uranium containing 0.78 kg of ^{235}\text{U} is:

m_{\text{natural U}} = \frac{0.78}{0.0072} \approx 108 \text{ kg/day}

Why: the other 99.28% is ^{238}\text{U}, which mostly does not fission (though some converts to plutonium and contributes a few per cent to the total energy). The rest sits as spent fuel.

Result: Kaiga-4 fissions about 2.3 \times 10^{19} atoms of ^{235}\text{U} per second, consumes about 0.78 kg of ^{235}\text{U} per day, and draws on about 108 kg of natural uranium per day to keep going. Over a fuel cycle of about 9 months, one refuelling covers about 30 tonnes of natural uranium — a small fleet of trucks.

What this shows: The 200-MeV-per-fission number you read off the binding-energy curve, combined with Avogadro's number and a given reactor power, gives you the fuel-consumption rate directly. No empirical tables, no reactor-specific corrections — just energy bookkeeping.

Example 2: How long to thermalise a fission neutron

A 2-MeV prompt fission neutron enters a heavy-water moderator. Each elastic collision with a deuteron reduces the neutron's average kinetic energy by a factor \alpha where, for a head-on elastic collision between a neutron (mass 1) and a target nucleus (mass A),

\alpha = \left(\frac{A-1}{A+1}\right)^2

(This is the elastic-collision kinetic-energy transfer derived in Elastic Collisions.) For deuterium, A = 2, so \alpha = (1/3)^2 = 1/9. Averaged over all collision angles, the mean energy after n collisions is

E_n = E_0 \cdot \bar{\alpha}^n

with \bar{\alpha} \approx 0.56 for deuterium. How many collisions are needed to bring a 2-MeV neutron to thermal (0.025 eV)?

Neutron energy vs collision number in heavy-water moderatorA log-scale plot of neutron energy in eV vs collision count, from 2×10^6 eV at collision 0 down to 0.025 eV at about collision 31, with deuterium (slope −0.25 per collision) steeper than carbon (slope −0.16). collisions n log₁₀ E (eV) 6 (MeV) 4 2 0 (eV) −2 (thermal) 0 10 20 30 40 D₂O (≈ 31 collisions) graphite (≈ 115 collisions)
Neutron energy on a log scale as a function of collision count. Deuterium (red) kills a neutron's kinetic energy quickly because its mass is close to the neutron's. Graphite (dark) is slower per collision but has a tiny absorption cross-section, so in total-thermalisation-probability terms it is roughly comparable.

Step 1. Set up the slowing-down equation.

E_n = E_0 \cdot \bar{\alpha}^n

Take the logarithm of both sides:

\ln E_n = \ln E_0 + n \ln \bar{\alpha}
n = \frac{\ln(E_n / E_0)}{\ln \bar{\alpha}}

Why: a fractional loss per collision becomes a geometric sequence; taking logs linearises it. This is the standard slowing-down theory in any reactor physics textbook.

Step 2. Substitute numbers. E_0 = 2\times 10^6 eV, E_n = 0.025 eV, \bar{\alpha} = 0.56.

\frac{E_n}{E_0} = \frac{0.025}{2\times 10^6} = 1.25 \times 10^{-8}
\ln(1.25\times 10^{-8}) = -18.20
\ln(0.56) = -0.58

Step 3. Divide.

n = \frac{-18.20}{-0.58} \approx 31 \text{ collisions}

Why: 31 collisions sounds like a lot, but each takes only about a microsecond in liquid D₂O. Total thermalisation time is of order 30 µs, which is why the prompt neutron generation time \tau \sim 10^{-4} s quoted earlier is reasonable.

Step 4. Compare with graphite (A = 12). \bar{\alpha}_{\text{C}} \approx 0.716, so \ln(0.716) = -0.334, giving n = -18.20/-0.334 \approx 55 collisions.

Why: graphite needs almost twice as many collisions to thermalise a neutron, because each collision only shaves off 28% of the energy. Ordinary water (H, A=1) is the extreme: \bar{\alpha}_{\text{H}} \approx 0.5 (in fact \alpha_{\max} = 1 for A=1 — a single head-on collision can stop the neutron entirely), requiring just n \approx 20 collisions. Pure hydrogen would be the best moderator by collision count — except for its fatal habit of absorbing thermal neutrons.

Result: A fission-energy neutron in heavy water undergoes about 31 elastic collisions before it reaches thermal equilibrium with the moderator. In graphite, about 55. In light water, about 20 — but light water absorbs too many to permit natural-uranium fuel.

What this shows: The choice of moderator is a quantitative trade-off between slowing-down power (how many collisions to thermalise, favouring light nuclei) and absorption cross-section (how many neutrons you lose en route, favouring heavier non-capturing nuclei like deuterium and carbon). The moderating ratio \xi \Sigma_s / \Sigma_a captures both: deuterium wins.

Common confusions

If you came here to understand what fission is, how chain reactions work, and why reactors need moderators and control rods, the material above is enough. What follows is for readers who want the four-factor formula, the Indian thorium cycle, and the physics of prompt critical.

The four-factor formula for k_\infty

The infinite-medium multiplication factor for a thermal reactor decomposes into four dimensionless factors:

k_\infty = \eta \cdot \epsilon \cdot p \cdot f
  • \eta = neutrons produced per thermal neutron absorbed in fuel = \nu \cdot \Sigma_f^{fuel} / \Sigma_a^{fuel}. For ^{235}\text{U} in a thermal spectrum, \eta \approx 2.07. If all absorptions in fuel led to fission, \eta would equal \nu, but some absorptions are radiative capture (^{235}\text{U} + n \to ^{236}\text{U} + \gamma), which wastes the neutron.

  • \epsilon = fast fission factor. Some fast neutrons, before thermalising, cause ^{238}\text{U} to fission. Each such fast fission adds to the neutron count. For natural-uranium fuel, \epsilon \approx 1.03 — a small boost.

  • p = resonance escape probability. As a neutron slows down, it passes through energies where ^{238}\text{U} has huge resonance absorption peaks (at 6.67 eV, 20.9 eV, 36.7 eV, …). The probability of missing every resonance is p. For natural-uranium PHWRs, p \approx 0.87.

  • f = thermal utilisation. Of all thermal neutrons absorbed, the fraction absorbed in fuel (rather than moderator, coolant, or structure) is f. For a well-designed PHWR, f \approx 0.92.

Multiply: k_\infty \approx 2.07 \times 1.03 \times 0.87 \times 0.92 \approx 1.71. This is why natural-uranium PHWRs work. In principle, pure natural-uranium surrounded by heavy water gives k_\infty > 1; in a finite core, leakage brings k_{\text{eff}} down to 1 by design.

The four-factor formula makes clear why enrichment matters for light-water reactors: in H₂O, the absorption in hydrogen is high enough that f falls below 0.7 for natural uranium — killing k_\infty. Enriching ^{235}\text{U} to 3–5% raises \eta and overcomes the f deficit.

The thorium cycle — India's long game

India has the world's largest reserves of thorium (estimated 25% of the global total, mostly in the monazite sands of Kerala). But thorium is not directly fissile — ^{232}\text{Th} is fertile, not fuel. The trick is a breeding cycle:

^{232}\text{Th} + n \longrightarrow ^{233}\text{Th} \xrightarrow{\beta^-, 22 \text{ min}} ^{233}\text{Pa} \xrightarrow{\beta^-, 27 \text{ days}} ^{233}\text{U}

The end product ^{233}\text{U} is fissile — it fissions on thermal neutrons with \eta \approx 2.3, better than both ^{235}\text{U} and ^{239}\text{Pu}.

India's three-stage nuclear programme, designed by Homi Bhabha in the 1950s, is:

  1. Stage 1: Natural-uranium PHWRs (Kaiga, Kakrapar, Rawatbhata, Narora). Run on natural U, breed some Pu.
  2. Stage 2: Fast breeder reactors fuelled by the Pu from Stage 1, with thorium blankets that breed ^{233}\text{U}. The Prototype Fast Breeder Reactor (PFBR) at Kalpakkam, built by IGCAR, is the proof of concept.
  3. Stage 3: Thorium-^{233}\text{U} advanced heavy-water reactors — closing the cycle on India's thorium reserves for centuries of energy.

The Kalpakkam PFBR achieved first criticality after years of delay in BARC's programme. Stage 3 remains a research target: the Advanced Heavy Water Reactor (AHWR) design at BARC is prototype-only so far. The long timescale is the price of the thorium path — 232Th needs the breeder catalyst, not available off-the-shelf. But the physics is sound, and the strategic logic (energy independence on domestic fuel) is compelling.

Prompt critical and the Chernobyl lesson

Define the reactivity:

\rho = \frac{k_{\text{eff}} - 1}{k_{\text{eff}}}

A reactor at \rho = 0 is critical. A reactor at \rho = \beta is prompt critical — it can achieve criticality on prompt neutrons alone, without needing the delayed component.

The reactor period T (time for power to multiply by e) for small positive reactivity is, including the delayed-neutron physics,

T \approx \frac{\tau_{\text{delay}}}{\rho / \beta}

For \rho = 0.5 \beta (half prompt-critical), T \approx 20 seconds — controllable by rods moving at a few cm/s.

For \rho = \beta (prompt critical), the formula diverges — the reactor's effective time constant crashes to \tau_{\text{prompt}} / \rho, on the order of milliseconds.

The Chernobyl RBMK reactor had a positive void coefficient — if the coolant boiled (creating voids in the core), reactivity increased, because the light-water coolant had been partly moderating the neutrons (which RBMKs didn't need; they used graphite). In the April 1986 accident, a combination of test conditions and operator errors caused coolant to flash to steam, reactivity to spike past \beta, and the core to go prompt critical. The power rose by a factor of 100 in less than a second. The steam explosion that followed destroyed the reactor.

Every modern reactor is designed with a negative void coefficient (and a negative fuel-temperature coefficient, the Doppler effect — as fuel heats up, ^{238}\text{U} resonance absorption broadens and steals more neutrons, reducing k). PWRs, BWRs, and PHWRs all achieve this by physics rather than by operator action. The Indian PHWR has a small positive void coefficient in the coolant but a much larger negative fuel-temperature coefficient; the net is slightly negative, and additional shutdown systems provide defence in depth.

The Fukushima accident in 2011 was not a criticality accident — the reactors SCRAMmed (all control rods fully inserted) within seconds of the earthquake. The subsequent meltdowns were due to loss of cooling: decay heat from the fission products (about 6% of full power, falling slowly over hours) overwhelmed the inadequate backup cooling after the tsunami. This is the other physics lesson of reactor safety: shutdown is not enough — you must also remove decay heat for days.

Why fission released energy at all — the deeper answer

Look at the binding-energy-per-nucleon curve one more time. Why does it peak at iron, rise steeply for light nuclei, and fall slowly for heavy nuclei? Three terms in the semi-empirical mass formula explain it:

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)
  • Volume term (+a_V A): bulk binding, constant per nucleon. This dominates.
  • Surface term (-a_S A^{2/3}): surface nucleons have fewer neighbours. For small A, this is a big fraction — that is why light nuclei have low B/A.
  • Coulomb term (-a_C Z(Z-1)/A^{1/3}): grows as Z^2/A^{1/3} for heavy nuclei. This is why heavy nuclei have low B/A.
  • Asymmetry term (-a_A (A-2Z)^2/A): nuclei prefer N \approx Z.

The peak at A \approx 56 is where surface and Coulomb penalties cross over. Fission moves a large nucleus with big Coulomb penalty down to mid-mass nuclei with less Coulomb penalty. The energy freed is precisely the Coulomb energy that was bottled up in the big nucleus.

You can even estimate the fission energy directly from the formula. Split ^{236}\text{U} into two equal halves (A = 118, Z = 46). The Coulomb-term change is roughly

\Delta E_C = a_C \left[\frac{92 \cdot 91}{236^{1/3}} - 2 \cdot \frac{46 \cdot 45}{118^{1/3}}\right]

With a_C \approx 0.72 MeV, you get \Delta E_C \approx 0.72 \times (8372/6.17 - 4140/4.90) \approx 0.72 \times (1357 - 845) \approx 370 MeV.

The surface-term change offsets this: \Delta E_S = a_S(A^{2/3} - 2 (A/2)^{2/3}) \approx 17.2 \times (38.0 - 2 \times 23.9) \approx 17.2 \times (-9.8) \approx -169 MeV.

Net: 370 - 169 = 201 MeV. The 200 MeV you read off the B/A curve is precisely the liquid-drop balance between the Coulomb gain and the surface cost — the same physics that sets the fission barrier height, and the same physics Bohr and Wheeler wrote down in 1939.

Where this leads next