⚛️ nuclear theory

Energy-Resolved Spin and Parity Distributions of Fission Fragments

This paper presents the first microscopic, energy-resolved simulation of spin and parity distributions in $^{239}$ Pu fission fragments, revealing their non-monotonic evolution with excitation energy and demonstrating how these distributions predict prompt photon multiplicities without adjustable parameters.

Original authors: Petar Marević, Nicolas Schunck, Antonio Bjelčić, Guillaume Scamps

Published 2026-10-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Petar Marević, Nicolas Schunck, Antonio Bjelčić, Guillaume Scamps

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Nuclear fission is the process where a heavy atom, like uranium or plutonium, splits into two lighter pieces. For decades, scientists have known that when this happens, the two resulting pieces, called fragments, do not just fly apart; they also spin rapidly, like tiny tops. This spinning is a mystery because the original atom usually has very little spin to begin with. Understanding how this spin is generated is crucial for several reasons. It helps explain how heavy elements are created in the violent collisions of neutron stars, it improves the safety and efficiency of nuclear reactors by predicting how much heat is released, and it aids in detecting hidden nuclear materials. However, while scientists have measured the average spin of these fragments, they have not fully understood how that spin changes as the energy of the incoming particle that triggers the split increases, nor have they known the precise rules governing the spin and the symmetry of the fragments.

A team of researchers has now used advanced computer simulations to watch this splitting process unfold from start to finish, revealing how the spin and symmetry of the fragments evolve. By modeling the fission of plutonium-239 triggered by neutrons with energies ranging from very low thermal levels up to 10 million electron volts, they discovered that the spin of the fragments is not random. Instead, it is tightly linked to how stretched the nucleus is at the exact moment it breaks. As the energy of the incoming neutron increases, the nucleus stretches further before snapping, and this extra stretching directly causes the fragments to spin faster. The researchers found that this relationship holds true for individual fragments, meaning that the shape of the nucleus at the moment of rupture dictates the spin of the pieces it creates.

The study also tackled a long-standing assumption about the symmetry of these fragments. In many standard models, scientists assume that the fragments are equally likely to have a "positive" or "negative" symmetry, a property known as parity, much like flipping a coin. The new simulations show that this assumption is incorrect, especially at lower energies. At thermal energy levels, the fragments show a strong preference for one type of symmetry over the other, and this preference remains linked to their spin even as the energy increases. Only at higher energies does the symmetry begin to look more random, but even then, the connection between spin and symmetry persists. This finding suggests that the standard way of predicting how fragments decay is missing a key piece of the puzzle.

To test the importance of these new findings, the researchers plugged their detailed spin and symmetry data into a standard computer code used to predict how fragments release energy. When they used the new, more accurate microscopic data instead of the old, simplified formulas, the predicted number of gamma rays emitted by the fragments changed significantly. The new data predicted a much higher number of gamma rays at thermal energies than the old models did, and it showed a different pattern of how this number grows as the neutron energy increases. This indicates that the old formulas, which were tuned to match past experiments, may need to be recalibrated to account for the true, complex behavior of the fragments.

The researchers achieved this by simulating the entire fission event, from the initial state of the nucleus through the moment it splits and finally to the decay of the fragments. They did not just look at the most common outcomes but tracked thousands of different possible paths the nucleus could take. They found that the number of different types of fragments produced increases as the energy goes up, filling in gaps that were empty at lower energies. Crucially, they observed that the spin of the fragments does not increase in a simple, straight line with the energy of the incoming neutron. Instead, the spin rises quickly at first, then the rate of increase slows down and eventually levels off. This leveling off happens because the nucleus reaches a limit in how much it can stretch before breaking, regardless of how much more energy is added.

This work provides the first detailed, energy-resolved map of how spin and symmetry are distributed among fission fragments. It moves beyond simple averages to show how these properties vary for every specific type of fragment produced. The results suggest that the mechanism driving the spin is the geometry of the nucleus at the moment of rupture, rather than some other hidden factor. By showing that the spin is a direct consequence of the nuclear shape, the study offers a clearer picture of the physics at play. The findings also highlight that the assumption of random symmetry in current models is a significant oversimplification that affects our ability to predict nuclear behavior accurately.

The implications of these results extend to the tools used by nuclear scientists and engineers. The statistical codes that simulate how fragments decay rely on inputs about spin and symmetry. The study demonstrates that using the new, detailed microscopic inputs changes the predictions for observable quantities like the number of emitted photons. This suggests that the parameters in current models, which have been adjusted to fit older data, may need to be updated to reflect the true physics of the fission process. While the simulations are powerful, the researchers note that future experiments will be needed to confirm these specific predictions, particularly regarding how the spin changes for different fragment masses. The study opens a path toward a more precise understanding of nuclear fission, replacing guesswork with a detailed, physics-based description of how these atomic pieces behave.

Technical Summary: Energy-Resolved Spin and Parity Distributions of Fission Fragments

Problem Statement
Nuclear fission transforms a heavy nucleus with little or no angular momentum into two rapidly spinning fragments. While the generation of these large angular momenta has been observed for decades, the underlying mechanism remains a central open question in low-energy nuclear physics. Previous studies have largely been limited to single, low excitation energies (thermal fission) or have relied on statistical models that assume equiprobable parity distributions and phenomenological spin inputs. Consequently, the dependence of fragment spins and parities on the excitation energy of the fissioning nucleus (determined by incident neutron energy) has remained largely unexplored. Furthermore, the impact of microscopic spin-parity correlations on prompt fission observables, such as photon multiplicities, has not been quantified.

Methodology
The authors present a fully microscopic simulation of the fission process for the reaction $^{239}\text{Pu}(n,f)$ , covering incident neutron energies ( $E_n$ ) from thermal up to 10 MeV. The approach combines two advanced theoretical frameworks:

Time-Dependent Hartree-Fock-Bogoliubov (TDHFB): Used to simulate the dynamics from the saddle point to scission and fragment separation. The calculations employ the axialhohfb program with the SkM* Skyrme energy density functional and a mixed pairing interaction. To model excitation, an instantaneous quadrupole boost is applied to initial configurations, transferring collective energy into intrinsic degrees of freedom.
Time-Dependent Generator Coordinate Method (TDGCM): Implemented via the felix program, this is used to determine the population weights of the TDHFB trajectories based on the collective flux on the potential energy surface (PES) prior to the saddle point.

The simulation generates 60 trajectories (10 baseline configurations boosted to 6 different energy levels: thermal, 1, 3, 5, 7, and 10 MeV). For fully separated fragments, a joint seven-dimensional projection is performed on:

Angular momentum ( $J$ ) and its fission-axis component ( $K$ ).
Parity ( $\pi$ ).
Neutron ( $N$ ) and proton ( $Z$ ) numbers for both fragments and the total system.

Projections utilize the Pfaffian method in the rest frame of each fragment. The resulting distributions are weighted by TDGCM fluxes to yield energy-resolved spin and parity distributions for a broad range of fragment masses.

Key Results

Spin Evolution: The mean spin of individual fragments evolves non-monotonically with excitation energy, tracking the deformation of the nucleus at scission. A characteristic "sawtooth" pattern in mean spin versus fragment mass persists up to $E_n = 10$ MeV. The yield-weighted average spin increases by nearly $2\hbar$ between thermal and 7 MeV, with the steepest rise occurring within the first MeV above thermal energy, followed by a linear increase and eventual saturation.
Deformation Correlation: The study establishes a strong correlation between the nuclear elongation at scission and the resulting fragment spin. The average spin is a nearly linear function of the mean quadrupole deformation at scission ( $\langle q_{20} \rangle$ ), serving as a better predictor than the incident energy itself. This suggests that the geometry of the neck rupture governs the spin generation.
Parity Distributions: Contrary to the standard statistical assumption of equiprobable positive and negative parities, the microscopic results show strong deviations at thermal energy, particularly for even-even and odd-odd fragments. These deviations diminish as $E_n$ increases, becoming largely negligible by 10 MeV. However, even at 10 MeV, the correlation between parity and spin (natural vs. unnatural parity) remains significant for even-even fragments, indicating that the factorized form assumed in many decay models is not fully recovered.
Impact on Prompt Observables: When these microscopic distributions are fed into the feta framework (based on the Hauser-Feshbach code yahfc) to simulate fragment decay, the results differ significantly from phenomenological models.
- Using microscopic spin-parity inputs (with parity restoration) predicts a higher average photon multiplicity ( $\bar{N}_\gamma = 8.93$ ) compared to the phenomenological model ( $\bar{N}_\gamma = 7.36$ ) and static HFB inputs ( $\bar{N}_\gamma = 5.81$ ).
- The microscopic models predict a slower, linear increase in total photon multiplicity with incident energy ( $\sim 0.09 \text{ MeV}^{-1}$ ), whereas the phenomenological model predicts a much steeper rise.
- The study notes that the higher spin inputs lead to a stronger neutron-photon competition, resulting in a slight decrease in neutron multiplicity.

Significance and Claims
The paper claims to provide the first microscopic, energy-resolved spin and parity distributions of fission fragments for $^{239}\text{Pu}(n,f)$ . Its primary significance lies in demonstrating that:

Predictive Power: Parameter-free microscopic inputs can predict fission observables (like photon multiplicities) without relying on empirical fits to fission data, challenging the necessity of phenomenological spin formulas in current nuclear data evaluation codes.
Mechanism Insight: The non-monotonic evolution of spin and its correlation with scission geometry provide a microscopic explanation for the spin-energy dependence, linking it directly to the nuclear shape at the moment of rupture.
Model Limitations: The results highlight that standard statistical decay assumptions (specifically equiprobable parity and uncorrelated spin-parity distributions) are invalid at thermal energies and may require recalibration even at higher energies to accurately reproduce experimental observables.

The authors conclude that future experimental campaigns, such as those planned at the NFS facility of GANIL, should test these predictions, particularly the mass-dependent reversal of spin trends and the specific impact of spin-parity correlations on prompt gamma emission.

Technical Summary: Energy-Resolved Spin and Parity Distributions of Fission Fragments

More like this