Four-Factor Formula Calculator

The four-factor formula expresses the infinite multiplication factor as k∞ = η·ε·p·f. It was developed by Enrico Fermi to decompose neutron multiplication into four distinct physical processes: reproduction, fast fission, resonance escape, and thermal utilization. Each factor quantifies one step in the neutron life cycle from birth to the next generation. A product greater than 1 means the reactor can sustain a chain reaction in an unbounded medium.

The reproduction factor η is the number of fast neutrons produced per thermal neutron absorbed in the fuel. It equals ν·σ_f / σ_a where ν is the neutrons per fission, σ_f is the fission cross-section, and σ_a is the total absorption cross-section of the fuel. For pure U-235 at thermal energies η ≈ 2.065. For a fuel mixture, η is weighted by the absorption cross-sections of each fuel nuclide.

The fast fission factor ε accounts for additional fissions caused by fast neutrons before they are moderated to thermal energies. Primarily, fast neutrons can fission U-238 (which has a threshold at about 1 MeV) while slowing down in the fuel. Since these extra fissions add neutrons beyond those from thermal fissions, ε > 1 always. For a light water reactor, ε is typically 1.02 to 1.08; for a CANDU reactor with natural uranium it is closer to 1.02.

The resonance escape probability p is the fraction of neutrons that are moderated to thermal energies without being captured in resonance absorption bands, primarily in U-238. As neutrons slow down through the eV-to-keV range, U-238 has extremely large resonance absorption cross-sections. A well-moderated lattice with a high moderator-to-fuel ratio lets neutrons slow down rapidly through the resonance region, giving p close to 1. For a typical LWR, p is 0.80 to 0.90.

The thermal utilization factor f is the fraction of thermal neutrons absorbed in the fuel (as opposed to the moderator, structural materials, and control poisons). It equals Sigma_a(fuel) / Sigma_a(total) where Sigma_a is the macroscopic absorption cross-section. A higher fuel-to-moderator ratio increases f but decreases p, so there is an optimum lattice geometry. Typical values are 0.60 to 0.90 for power reactor cores.

k∞ is the infinite multiplication factor: the neutron multiplication in a hypothetical reactor of infinite size with no neutron leakage. k_eff (effective multiplication factor) applies to a real finite reactor and includes neutron leakage losses: k_eff = k∞ x P_NL where P_NL is the non-leakage probability. For a critical reactor k_eff = 1.000. Large power reactors have P_NL of 0.95 to 0.99; small research reactors may have P_NL of 0.85 to 0.95.

Reactivity is defined as ρ = (k - 1) / k. For a supercritical system (k > 1) reactivity is positive; for a subcritical system (k < 1) it is negative; for a critical system (k = 1) ρ = 0. Reactivity is commonly expressed in units of pcm (percent-milli, or 10^-5). For example, if k∞ = 1.05, then ρ = 0.05/1.05 ≈ 0.04762 = 4,762 pcm.

A light water reactor (LWR) with 3-5% enriched UO2 fuel typically has k∞ of 1.25 to 1.45 at beginning-of-life (before burnup reduces the fissile inventory). A CANDU reactor with natural uranium has k∞ of about 1.10 to 1.15. A graphite-moderated reactor with natural uranium (like the first Chicago Pile) was designed just above criticality. Research reactors using highly enriched uranium can have k∞ well above 1.5.

As U-235 is consumed, η decreases because fissile inventory drops. Simultaneously, Pu-239 builds up from U-238 neutron capture, partially compensating for the U-235 loss. The resonance escape probability p increases slightly as U-238 is depleted. Fission products (especially Xe-135 and Sm-149) absorb thermal neutrons, reducing f significantly. Collectively, k∞ drops from its beginning-of-life value toward 1 (criticality limit) by end-of-life.

The four-factor formula in its standard form applies strictly to thermal reactors where the neutron life cycle passes through a distinct thermal energy group. Fast reactors do not have a clearly defined thermal neutron population, so p and f as defined are not meaningful. Fast reactor analysis uses multigroup diffusion or transport methods. A related six-factor formula (adding fast leakage and thermal leakage) extends the concept to finite thermal reactors but still does not apply to fast spectra.

Adding moderator increases p (faster slowing-down through resonances) but decreases f (more neutrons absorbed in moderator). The product p x f passes through a maximum at an optimum ratio. For light water reactor fuel assemblies this optimum falls at a hydrogen-to-uranium ratio near the practical design range. Under-moderated lattices have lower p; over-moderated lattices have lower f. LWRs are deliberately designed slightly under-moderated so that adding water (moderator) at higher temperature decreases k∞, providing a negative moderator temperature coefficient.

The formula provides a structured framework for lattice optimization during the fuel assembly design phase. Nuclear engineers independently calculate or measure each factor from cross-section libraries and geometry, then iterate on fuel enrichment, fuel pin pitch, moderator density, and poison loading to achieve a target k∞ with appropriate shutdown margin. Each factor can be adjusted through separate design levers: fuel enrichment changes η, fuel pin diameter and spacing affect p and f, and burnable poisons alter f without permanently removing reactivity.