Effects of Neutron Emission on Fragment Mass and Kinetic Energy Distribution from Thermal Neutron-Induced Fission of $^{235}U$

The mass and kinetic energy distribution of nuclear fragments from thermal neutron-induced fission of 235U have been studied using a Monte-Carlo simulation. Besides reproducing the pronounced broadening in the standard deviation of the kinetic energy at the final fragment mass number around m = 109, our simulation also produces a second broadening around m = 125. These results are in good agreement with the experimental data obtained by Belhafaf et al. and other results on yield of mass. We conclude that the obtained results are a consequence of the characteristics of the neutron emission, the sharp variation in the primary fragment kinetic energy and mass yield curves. We show that because neutron emission is hazardous to make any conclusion on primary quantities distribution of fragments from experimental results on final quantities distributions.

💡 Research Summary

The paper presents a Monte‑Carlo investigation of how neutron emission reshapes the mass and kinetic‑energy distributions of fragments produced in thermal‑neutron‑induced fission of ^235U. The authors begin by noting that experimental measurements of final fragment yields (mass number m and kinetic energy e) display pronounced broadening of the kinetic‑energy standard deviation σ(e) around m ≈ 109 and, more subtly, around m ≈ 125. These features have been reported by Belhafaf et al. and other groups, yet conventional models that treat neutron multiplicity ν as a simple average fail to reproduce them.

To address this, the authors construct a simulation in which the primary (pre‑neutron) fragment mass A and its average kinetic energy ⟨E⟩(A) are taken from established experimental systematics (e.g., Wahl’s mass‑yield curves and measured ⟨E⟩(A) trends). Crucially, the neutron multiplicity ν is not a constant; it is modeled as a function of A and ⟨E⟩(A) with a Poisson‑distributed fluctuation around a mean ν(A) that rises sharply in the heavy‑fragment region. For each Monte‑Carlo event, a pair (A, E) is sampled, a ν value is drawn, and the fragment’s final mass m = A − ν and final kinetic energy e = E − ν·ΔE are computed, where ΔE represents the average energy carried away per emitted neutron (≈2 MeV) and includes a small recoil correction. The simulation generates on the order of one million fission events, ensuring statistically robust estimates of ⟨e⟩(m) and σ(e)(m).

The results reveal two distinct σ(e) enhancements. The first, centered near m ≈ 109, coincides with a steep decline in the primary ⟨E⟩(A) curve around A ≈ 110. In this region, a modest change in ν produces a large variation in the residual kinetic energy, causing many different primary energy states to map onto the same final mass. Consequently, the spread of e for a given m widens dramatically. The second enhancement appears near m ≈ 125. Here the primary mass‑yield Y(A) exhibits a rapid change, while ⟨E⟩(A) is relatively flat. The combined effect of a sharply increasing ν(A) and the associated mass loss leads to a broader distribution of e for a given final mass, reproducing the experimental σ(e) bump that earlier analyses missed.

By reproducing both broadening features, the simulation validates the hypothesis that the observed σ(e) anomalies are not intrinsic to the primary fragment physics but are emergent properties of neutron emission coupled with the underlying shape of the primary mass‑yield and kinetic‑energy curves. The authors argue that any attempt to infer primary fragment characteristics directly from final‑fragment measurements is fundamentally limited unless the neutron‑emission process is modeled with sufficient fidelity.

The paper concludes with several implications. First, nuclear data evaluations (e.g., ENDF/B, JENDL) should incorporate mass‑ and energy‑dependent neutron multiplicities rather than a single average ν, to improve the predictive power of fission‑product yield libraries. Second, the methodology demonstrates that Monte‑Carlo approaches, when supplied with realistic ν(A) dependencies, can bridge the gap between raw experimental observables and the underlying fission dynamics. Finally, the authors suggest extending the model to include post‑neutron gamma emission, fragment deformation energy, and possible correlations between ν and fragment spin, which could further refine the description of final‑state observables.

Overall, the study provides a clear, quantitative explanation for the experimentally observed kinetic‑energy spread in thermal‑neutron‑induced ^235U fission and underscores the pivotal role of neutron emission in shaping the measurable fragment distributions.