Resonance Escape Probability

The resonance escape probability p is the fraction of neutrons that slow down from fast energies (several MeV) to thermal energies (around 0.025 eV) without being absorbed in the resonance peaks of fissile and fertile isotopes (primarily U-238). It is a key factor in the four-factor formula (k∞ = η ε p f), directly influencing the neutron multiplication factor. In thermal reactors, p is typically around 0.85–0.95, meaning that 5–15% of the neutrons are lost to resonance absorption, primarily in U-238. Improving p through better moderation or fuel design is essential for achieving criticality and reducing enrichment requirements.

The exponent term is the resonance absorption integral, which quantifies the probability that a neutron will be absorbed in a resonance while slowing down. N_i is the number density of the absorbing nuclei (e.g., U-238), I_i is the resonance integral (a measure of the effective absorption cross-section over the resonance energy range), and the denominator ξ Σ_s is the moderating power (the product of the average logarithmic energy decrement ξ and the macroscopic scattering cross-section Σ_s). A higher moderating power means neutrons slow down faster, spending less time at resonance energies, thereby reducing absorption. The exponential form arises because the absorption probability is proportional to the total number of collisions, which follows Poisson statistics.

U-238 has a series of strong absorption resonances in the keV to eV energy range (particularly at 6.7 eV, 20.9 eV, and many higher resonances). These resonances are due to the capture cross-section spiking at specific neutron energies. While U-235 also has resonances, they are much weaker and less significant. The large abundance of U-238 in the fuel (typically 95% or more) means that even a small capture probability per neutron can lead to a significant loss of neutrons. This is why thermal reactors rely on efficient moderation to minimize resonance absorption and why fuels with lower U-238 content (e.g., highly enriched uranium or MOX) have higher p.

The resonance integral is the integral of the effective capture cross-section over energy, weighted by the neutron flux spectrum in the resonance region. It is typically expressed in barns (1 barn = 10⁻²⁸ m²). The resonance integral depends on the material's nuclear properties and the temperature (Doppler broadening). It is determined experimentally or via nuclear data libraries (e.g., ENDF/B, JEFF, JENDL). For U-238 at room temperature, the resonance integral is about 280 barns, but it increases with temperature due to Doppler broadening, which is why reactors operate with a 'temperature coefficient' that affects p and reactivity.

ξ = 1 + (A-1)²/(2A) ln((A-1)/(A+1)), where A is the atomic mass of the moderator nuclei. It represents the average change in the natural logarithm of neutron energy per scattering collision. For hydrogen (A=1), ξ = 1, meaning neutrons lose their energy most efficiently. For heavier moderators like carbon (A=12), ξ ≈ 0.158, so more collisions are needed to slow down a neutron. A higher ξ means faster energy loss, reducing the time spent at resonance energies, which increases p. This is why hydrogen-rich materials (water, polyethylene) are excellent moderators.

The moderating power ξ Σ_s varies significantly among moderators. For H₂O, ξ ≈ 0.92 and Σ_s is high due to high hydrogen density, giving a high moderating power, which increases p. For D₂O, ξ ≈ 0.51 (since deuterium is heavier) and Σ_s is lower, so the moderating power is about 5-10 times lower than H₂O. This means that in a heavy water reactor, neutrons slow down more slowly, spending more time at resonance energies, leading to a lower p (but D₂O also has lower absorption). Graphite has ξ ≈ 0.158 and moderate Σ_s, giving an intermediate moderating power. The combination of moderating power and neutron economy determines the overall reactivity. Thermal reactors often use H₂O because its high p allows lower enrichment.

Doppler broadening is the phenomenon where as temperature increases, the resonance peaks in the cross-section become wider (due to thermal motion of the nuclei). This increases the resonance integral (I_i) because neutrons are more likely to be absorbed over a broader energy range. As a result, the resonance escape probability p decreases with increasing temperature. This provides a negative temperature coefficient of reactivity (the Doppler effect), which is a safety mechanism: if the core heats up, p decreases, reducing reactivity and providing inherent stability. This is one of the most important feedback mechanisms in LWRs.

In a fast reactor, there is no significant moderation, so neutrons remain at fast energies and do not pass through the resonance region. Therefore, resonance absorption is minimized, and p is effectively 1 (or very close to 1). This is why fast reactors can use fuels with a high U-238 content (e.g., depleted uranium or MOX) without the penalty of resonance absorption. However, fast reactors require higher enrichment or a compact core to achieve criticality because the fission cross-sections are lower at fast energies. The resonance escape probability is primarily a concept for thermal and intermediate reactors where moderation occurs.

The four-factor formula is k∞ = η × ε × p × f. Here, η is the number of neutrons produced per absorption in the fuel (depending on the fuel type), ε is the fast fission factor (neutrons produced by fast fission in U-238), p is the resonance escape probability, and f is the thermal utilization factor (the fraction of thermal neutrons absorbed in the fuel). These factors are multiplicative; a small decrease in p (e.g., from 0.90 to 0.85) reduces k∞ by approximately the same percentage, which can be the difference between critical and subcritical. p is often the most sensitive factor to design changes like moderator-to-fuel ratio, fuel enrichment, and temperature.

In a typical PWR, p is about 0.88–0.92 at the beginning of cycle (BOC) for 4.5% enriched fuel. As burnup progresses, the U-235 concentration decreases, and the concentration of fission products and plutonium (which also have resonances) increases. However, the net effect on p is complex: the reduction in U-238 (via conversion to Pu-239) slightly reduces the resonance integral, while the buildup of Pu-240 and other actinides introduces new resonances. Overall, p may increase slightly with burnup (by a few percent) because the U-238 concentration decreases, but the change is modest. Reactor physics codes track p as part of the depletion calculation, and it is used in fuel cycle analysis.

Burnable poisons like gadolinium have resonances in the thermal and epithermal ranges. They are primarily designed to absorb thermal neutrons to control excess reactivity, but they can also absorb neutrons in the resonance region. The effect on p depends on the poison's placement and concentration. Typically, p is calculated for the whole core; adding a poison that absorbs neutrons (even at resonance energies) increases the total absorption, which could reduce p, but the effect is usually small because the poison's resonance absorption is minor compared to U-238. The dominant effect is on the thermal utilization factor (f), not p. However, in advanced designs, the interplay between p and f is optimized.

They are essentially the same concept, but the resonance absorption integral is the product of the resonance integral (I_i) and the number density (N_i) integrated over the system. In the exponent term Σ N_i I_i, the sum over all absorber nuclides gives the total resonance absorption probability. The resonance integral I_i itself is the energy-weighted cross-section over the resonance region, often expressed as ∫ σ(E) dE/E in the slowing-down theory. It depends on the temperature and the type of absorber. The product N_i I_i (sometimes called the 'resonance absorption cross-section' per unit volume) is used in the formula for p.

This is a simplified approximation used in elementary reactor physics, where β represents the effective resonance absorption probability relative to the scattering probability. The exact formula is p = exp(-Σ N_i I_i / ξ Σ_s). For small absorption (weak resonances), the exponential can be approximated by p ≈ 1 - Σ N_i I_i / ξ Σ_s, and in some simplified models, this is written as p = 1/(1 + β) with β = Σ N_i I_i / ξ Σ_s. This approximation is less accurate for large resonance absorption (e.g., for large fuel loadings) but is useful for educational purposes. The exponential form is more accurate and is used in modern physics codes.

Resonance escape probability cannot be measured directly; it is inferred from the overall reactivity balance. Reactor physicists measure the multiplication factor (k) through criticality experiments or by observing reactor period and control rod worth. Using the four-factor formula, p can be derived if the other factors (η, ε, f) are known from nuclear data and calculations. In practice, p is computed using transport theory codes (e.g., MCNP, SCALE) with ENDF cross-section data, and the results are validated against critical experiments (e.g., the F-1 or TRX benchmarks). The uncertainty in p is typically a few percent, which contributes to the overall uncertainty in core design.

Higher enrichment means more U-235 and less U-238 per unit mass of fuel. Since U-238 is the primary resonance absorber, reducing its number density N_i reduces the resonance absorption term Σ N_i I_i, which increases p. For example, increasing enrichment from 3% to 5% reduces the U-238 concentration by about 2% relative to the total, which can increase p by about 1-2%. This is one reason why higher enrichment allows a longer fuel cycle: the improved p contributes to higher reactivity and better neutron economy. However, the effect is modest because the fuel's overall resonance absorption is also influenced by the fuel geometry and moderator-to-fuel ratio.