Four-Factor Formula Calculator

Compute the infinite multiplication factor k∞ from the four-factor product η·ε·p·f for thermal reactor physics analysis.

⚛️ Four-Factor Formula Calculator
1.03.0
1.001.20
0.101.00
0.100.99
k∞ (Infinite Multiplication Factor)
Factor Product (η × ε × p × f)
Reactivity ρ (Δk/k)
Reactivity in pcm
Criticality Status
Required Factor
Condition
Physical Constraint

⚛️ What is the Four-Factor Formula?

The four-factor formula is the foundational equation of thermal reactor physics, expressing the infinite multiplication factor as k∞ = η × ε × p × f. It was formulated by Enrico Fermi and colleagues during the Manhattan Project to decompose neutron multiplication into four physically distinct processes that occur during one complete neutron lifetime, from birth as a fast fission neutron to absorption and the production of the next neutron generation. A reactor with k∞ greater than 1.000 can sustain a chain reaction in principle; the actual criticality of a finite reactor also depends on how many neutrons leak from the core boundary.

Each factor represents a separate stage of the neutron life cycle in a thermal reactor. First, η (eta, the reproduction factor) counts how many fast neutrons are born per thermal neutron absorbed in fuel. Next, ε (epsilon, the fast fission factor) accounts for additional fast fissions in U-238 before neutrons slow down, boosting the count above unity. Then p (the resonance escape probability) represents the fraction of neutrons that survive the dangerous resonance absorption region in U-238 while being moderated to thermal energies. Finally, f (the thermal utilization factor) is the fraction of thermal neutrons actually absorbed in the fuel rather than in moderator, structural materials, or control poisons. The product of these four quantities gives k∞.

The formula applies specifically to thermal reactors, where the neutron spectrum has a distinct thermal population separated from the fast and epithermal regions. It does not directly apply to fast reactors, which lack a well-defined thermal group. For finite reactors, the six-factor formula extends the four-factor formula by adding a fast non-leakage probability and a thermal non-leakage probability to account for geometric leakage: k_eff = η × ε × p × f × P_FNL × P_TNL. For large power reactors, the non-leakage factor P_NL = P_FNL × P_TNL is typically 0.95 to 0.99, so k_eff is close to k∞.

Understanding and optimizing the four factors is central to reactor core design. Nuclear engineers independently adjust each factor through fuel enrichment (which changes η), lattice geometry and moderator ratio (which affects both p and f), and control poisons (which reduce f). The formula also provides a diagnostic framework: if a reactor's measured k_eff is lower than expected, the four-factor decomposition tells engineers which physical process is responsible so they can target the correct design parameter.

📐 Formula

k∞ = η × ε × p × f
k∞ = infinite multiplication factor (dimensionless)
η (eta) = reproduction factor: fission neutrons per thermal neutron absorbed in fuel = ν × σ_f / σ_a(fuel)
ε (epsilon) = fast fission factor: total fission neutrons / thermal fission neutrons (always ≥ 1)
p = resonance escape probability: fraction surviving resonance capture while moderating to thermal energies (0 < p ≤ 1)
f = thermal utilization factor: fraction of thermal neutrons absorbed in fuel = Σ_a(fuel) / Σ_a(total) (0 < f < 1)
Reactivity: ρ = (k∞ − 1) / k∞
Reactivity in pcm: ρ(pcm) = ρ × 105
For finite reactors: k_eff = k∞ × P_NL
P_NL = non-leakage probability = exp(−B²L²) / (1 + B²M²) depending on the model
= geometric buckling (depends on core dimensions)
= migration area (material property of the moderator-fuel lattice)

📖 How to Use This Calculator

Steps

1
Select calculation mode - use "Calculate k∞" to find the multiplication factor from all four inputs, or "Solve for Missing Factor" to determine one unknown factor given a target k∞.
2
Enter the four factor values - type eta (η), epsilon (ε), resonance escape probability (p), and thermal utilization (f) using the sliders or numeric text boxes. Default values represent a typical light water reactor.
3
Click Calculate - the calculator returns k∞, reactivity ρ in both percentage and pcm, and a criticality status: subcritical (k < 1), critical (k ≈ 1), or supercritical (k > 1).
4
Use reverse mode for design targets - switch to "Solve for Missing Factor", set your target k∞, enter the three known factors, and select which factor to solve for. This is useful for lattice optimization: for example, finding the required thermal utilization f to achieve a target k∞ given fixed η, ε, and p.

💡 Example Calculations

Example 1 - Typical Light Water Reactor (LWR) Lattice

LWR with 3.5% enriched UO2: η=2.07, ε=1.05, p=0.85, f=0.70

1
η × ε = 2.07 × 1.05 = 2.1735
2
2.1735 × p = 2.1735 × 0.85 = 1.8475
3
1.8475 × f = 1.8475 × 0.70 = 1.2932
4
ρ = (1.2932 − 1) / 1.2932 = 0.2268 = 22,680 pcm
k∞ = 1.2932 | Reactivity = 22,680 pcm | Status: Supercritical
Try this example →

Example 2 - CANDU Reactor with Natural Uranium

CANDU natural uranium lattice: η=1.73, ε=1.02, p=0.92, f=0.81

1
η × ε = 1.73 × 1.02 = 1.7646
2
1.7646 × p = 1.7646 × 0.92 = 1.6234
3
1.6234 × f = 1.6234 × 0.81 = 1.3150
4
ρ = (1.3150 − 1) / 1.3150 = 0.2395 = 23,954 pcm
k∞ = 1.3150 | Reactivity = 23,954 pcm | Status: Supercritical
Try this example →

Example 3 - Subcritical Assembly with High Parasitic Absorption

Heavy control rod insertion: η=2.07, ε=1.03, p=0.60, f=0.55

1
η × ε = 2.07 × 1.03 = 2.1321
2
2.1321 × p = 2.1321 × 0.60 = 1.2793
3
1.2793 × f = 1.2793 × 0.55 = 0.7036
4
ρ = (0.7036 − 1) / 0.7036 = −0.4215 = −42,150 pcm (deep subcritical)
k∞ = 0.7036 | Reactivity = −42,150 pcm | Status: Subcritical
Try this example →

❓ Frequently Asked Questions

What is the four-factor formula in nuclear reactor physics?+
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.
What does eta (η) represent in the four-factor formula?+
The reproduction factor η is the number of fast neutrons produced per thermal neutron absorbed in the fuel. It equals ν·σ_f / σ_a, where ν is neutrons released 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 mixed fuel, η is computed as a cross-section-weighted average over all fissile nuclides.
Why is the fast fission factor epsilon always greater than or equal to one?+
The fast fission factor ε accounts for extra neutrons produced when fast neutrons (above the U-238 fission threshold of about 1 MeV) cause fissions before being moderated to thermal energies. Since these additional fast fissions only add neutrons, ε is always ≥ 1. For light water reactors, ε is typically 1.02 to 1.08. For a purely thermal reactor with no fast fissions at all, ε would equal exactly 1.
What physical parameter most strongly controls the resonance escape probability p?+
The resonance escape probability p is most sensitive to the moderator-to-fuel ratio and to the fuel temperature. A high moderator-to-fuel ratio promotes rapid moderation through the resonance region, increasing p. Higher fuel temperature broadens the Doppler width of U-238 resonances, increasing resonance absorption and decreasing p. This Doppler effect on p is the primary source of the negative fuel temperature coefficient of reactivity that gives thermal reactors passive safety.
What is the difference between k∞ and k_eff?+
k∞ is the multiplication factor in a hypothetical infinite reactor with no neutron leakage from boundaries. k_eff applies to a real finite reactor: k_eff = k∞ × P_NL, where P_NL is the non-leakage probability. For a critical reactor, k_eff = 1.000 exactly. Large power reactors have P_NL of 0.95 to 0.99, so k_eff is only slightly below k∞. Small research reactors or subcritical assemblies may have P_NL as low as 0.85.
How is reactivity defined and what are its units?+
Reactivity ρ = (k − 1) / k. It is dimensionless, with ρ = 0 at criticality, positive when supercritical, and negative when subcritical. Common units are pcm (percent-milli = 10^−5), so 1 pcm corresponds to ρ = 0.00001. Another unit is the dollar ($), equal to the delayed neutron fraction β (approximately 650 pcm for U-235 fuel). Prompt criticality, where only prompt neutrons are needed to sustain the chain reaction, occurs at ρ = β.
How does fuel burnup change the four factors over a fuel cycle?+
As U-235 is consumed during burnup, η falls because the fissile inventory decreases. Pu-239 buildup from U-238 neutron capture partially compensates but does not fully restore η. Fission products accumulate and absorb thermal neutrons, significantly reducing f. The resonance escape probability p increases slightly as U-238 inventory drops. The net effect is a gradual decrease in k∞ from beginning-of-life values (often 1.3 to 1.5) toward the end-of-life criticality limit near 1.0.
Can the four-factor formula be applied to fast reactors?+
No. The four-factor formula applies only to thermal reactors where the neutron spectrum has a distinct thermal population. Fast reactors operate entirely in the fast and epithermal energy range, so the concepts of resonance escape probability and thermal utilization are not meaningful as defined. Fast reactor analysis uses multigroup diffusion theory or Monte Carlo transport codes with detailed cross-section libraries across dozens of energy groups.
What is the optimum moderator-to-fuel ratio and why are LWRs under-moderated?+
There is a moderator-to-fuel ratio that maximizes the product p × f and therefore k∞. Above this optimum, more moderator increases p but decreases f; below it, less moderator increases f but decreases p. Light water reactors are deliberately designed slightly below the optimum (under-moderated) so that reducing moderator density (as coolant heats or boils) decreases k∞. This gives a negative moderator temperature coefficient, which is a fundamental passive safety feature.
What values of k∞ are typical for common reactor types?+
A light water reactor with 3-5% enriched UO2 fuel typically has k∞ of 1.25 to 1.45 at beginning-of-life. A CANDU reactor with natural uranium has k∞ of about 1.10 to 1.15. High-temperature gas-cooled reactors with graphite moderator and enriched fuel can reach k∞ of 1.4 to 1.6. Highly enriched uranium research reactor lattices can have k∞ above 1.6, providing large shutdown margins.
How does boron concentration in the coolant affect the four factors?+
Dissolved boron in the coolant (used in pressurized water reactors for long-term reactivity control) primarily reduces the thermal utilization factor f. Boron-10 has a very large thermal neutron absorption cross-section, so it competes with the fuel for thermal neutrons. Increasing boron concentration reduces f and therefore k∞, compensating for the excess reactivity at beginning-of-life. As fuel burns up, boron concentration is slowly reduced to maintain criticality throughout the fuel cycle.
How is the four-factor formula used in reactor lattice design?+
The formula provides a structured framework for lattice optimization. Nuclear engineers calculate or measure each factor independently using neutron cross-section libraries and geometric models, then iterate on fuel enrichment (which changes η), fuel pin pitch and moderator volume (affecting both p and f), and burnable poison loading (reducing f to manage excess reactivity at beginning-of-life). Each factor is a separate design lever, making the four-factor decomposition essential for understanding which physical mechanism limits reactor performance.