Narrowing the Gap Between Theory and Evaluations: Angular Momentum Distributions in Fission Fragments

This paper presents a parameter-free microscopic framework that successfully predicts angular momentum distributions and photon multiplicities for neutron-induced fission of $^{235}$U and $^{239}$Pu, demonstrating that theoretical models are now quantitatively competitive with phenomenological approaches.

The Gist

Imagine a giant, unstable balloon (an atomic nucleus) that suddenly pops, splitting into two smaller, spinning balloons (fission fragments). For a long time, scientists knew these smaller balloons spun, but they didn't have a precise way to predict how fast or in what pattern they would spin.

This paper is like a new, high-definition camera that finally captures the exact spinning motion of these fragments right at the moment the big balloon splits. Here is the story of what the researchers found, explained simply:

The Old Problem: Guessing vs. Knowing

For decades, scientists had two ways to understand this splitting:

1. The "Guessing Game" (Phenomenological Models): They used simple rules and adjusted knobs until their predictions matched what they saw in experiments. It worked well, but it was more like tuning a radio to get a clear signal than understanding how the radio works.
2. The "Deep Dive" (Microscopic Theory): They tried to calculate everything from the very bottom up, using the fundamental laws of physics. This was the "holy grail," but the math was so incredibly complex that the computers of the past couldn't handle it. The results were often too fuzzy to be useful.

The Breakthrough: Thanks to massive leaps in computer power, the authors (Petar Marević, Nicolas Schunck, and Marc Verriere) finally built a "deep dive" model that is now just as accurate as the "guessing game." They didn't need to tweak any knobs; they just let the laws of physics do the work.

How They Did It: The "Splitting Moment"

To predict the spin, the team didn't just look at the final result; they simulated the exact moment the nucleus splits (called "scission").

• The Analogy: Imagine stretching a piece of taffy until it's about to snap. The team calculated thousands of different ways the taffy could stretch and thin out.
• The Calculation: For every possible way the nucleus could split, they calculated the probability of the two resulting pieces having a specific amount of spin (angular momentum). They combined all these possibilities to create a complete map of how the fragments spin.

The Surprising Patterns

When they looked at their new map, they found three cool things:

1. The "Sawtooth" Dance: As the size of the fragments changes, their average spin doesn't go up or down smoothly. Instead, it zig-zags up and down like the teeth of a saw. This pattern was known to exist, but their theory predicted it perfectly without any help.
2. The "Sibling" Effect: Even if two fragments have the same total weight, they don't always spin the same way. If one is made of a specific mix of protons and neutrons (like a specific "sibling" in a family), it might spin wildly, while its "sibling" with a slightly different mix spins slowly. This is called isobaric dependence.
   • The Metaphor: Think of it like two identical-looking spinning tops. If one has a tiny weight hidden inside a specific spot, it spins differently than the other, even if they look the same from the outside.
3. No "Tuning" Required: The most impressive part is that they didn't adjust their model to fit the data. They just ran the simulation, and the results matched real-world measurements of how many photons (light particles) are emitted when the fragments cool down.

Why This Matters

Before this, if scientists wanted to simulate how these fragments decay (cool down) in a computer program, they had to rely on those old "guessing game" models with adjustable knobs.

In this paper, the authors took their new, "no-knob" microscopic predictions and fed them into a standard simulation program (called cgmf).

• The Result: The simulation predicted the number of light particles (photons) emitted almost exactly right.
• The Takeaway: This proves that "deep dive" physics is finally ready to compete with the old "guessing" methods. It's a major step forward because it means we can now trust our fundamental understanding of the universe to predict complex nuclear events, rather than just relying on trial and error.

What They Didn't Do

The paper is very careful to say what they didn't do:

• They did not invent a new medical treatment or a new power plant design.
• They did not claim to solve all nuclear physics problems.
• They noted that their model still has some limitations (like ignoring certain tiny rotational effects), but for the main question of "how much do these fragments spin?", the answer is now solid.

In a nutshell: The authors built a super-accurate, physics-based crystal ball that predicts how atomic fragments spin after a split. It works so well that it matches real experiments without needing any "cheat codes" or adjustments, proving that our deep understanding of nature is finally catching up to our practical needs.

Technical Summary

Technical Summary: Angular Momentum Distributions in Fission Fragments

Problem Statement
Nuclear fission theory has historically struggled to bridge the gap between fundamental microscopic modeling and practical applications. While phenomenological models have long provided superior quantitative predictions for fission observables, microscopic approaches have often been limited to qualitative insights. A critical bottleneck in simulating the decay of fission fragments (FFs) within statistical reaction theory codes (e.g., cgmf, freya, fifrelin) is the lack of a robust microscopic framework to predict angular momentum distributions across the full range of fragment masses and charges. Currently, these simulations rely on phenomenological distributions with adjustable parameters fitted to experimental spectra, which may fail in regions where experimental data is scarce. Furthermore, previous theoretical studies of FF angular momentum were largely restricted to a narrow mass range around the most probable fragmentation, leaving the isobaric dependence and global distribution characteristics underexplored.

Methodology
The authors employ a fully microscopic framework based on Energy Density Functional (EDF) theory, specifically utilizing the hfbtho computational code with Skyrme EDFs (SkM*) and a mixed volume-surface contact pairing force. The approach involves three primary steps:

1. Scission Configuration Generation: Constrained Hartree-Fock-Bogoliubov (HFB) calculations are performed on even-even compound systems ($^{236}$U$^*$ and $^{240}$Pu$^*$) to generate scission configurations. These configurations are defined by collective variables: quadrupole moment ($q_{20}$), octupole moment ($q_{30}$), and neck parameter ($q_N$).
2. Angular Momentum Projection: For each scission configuration, projection techniques are applied to extract the probability distribution of angular momentum ($J_F$), neutron number ($N_F$), and proton number ($Z_F$) for the left and right fragments, conditioned on the total system's nucleon numbers.
3. Dynamical Weighting: The probability of populating each scission configuration, $F(q)$, is estimated using the Time-Dependent Generator Coordinate Method (TDGCM) with the Gaussian Overlap Approximation (GOA). The final angular momentum distribution is obtained by summing the projected distributions weighted by the dynamical probabilities of the scission configurations.

The model was applied to neutron-induced fission of $^{235}$U and $^{239}$Pu at an incident neutron energy of $E_n = 1$ MeV. The resulting distributions were subsequently used as inputs for the cgmf statistical decay code to simulate photon and neutron multiplicities, without adjusting any free parameters.

Key Contributions and Results
The study presents the first microscopic prediction of angular momentum distributions covering the entire range of fission fragment masses and charges for $^{235}$U and $^{239}$Pu. Key findings include:

• Sawtooth Pattern and Shell Effects: The calculated average angular momentum exhibits a pronounced sawtooth pattern as a function of fragment mass, consistent with experimental observations. The distributions reveal clear signatures of shell effects; for instance, fragments near magic neutron numbers ($N=82$) show a maximum at $J=0$, while moving away from closed shells significantly increases angular momentum.
• Substantial Isobaric Dependence: A significant finding is the strong isobaric dependence of angular momentum. Within a single isobaric chain (fixed mass number $A_F$), the average angular momentum can vary by up to $\pm 3\hbar$ depending on the specific neutron and proton composition. This contrasts with standard phenomenological models (e.g., Eq. 5 in the paper), which typically disregard isobaric variations, considering only mass and temperature dependencies.
• Quantitative Agreement with Experiment: When the microscopic distributions were used as inputs for cgmf, the resulting photon multiplicities for $^{240}$Pu$^*$ fission reproduced experimental data with high fidelity, requiring no adjustable parameters. The total photon multiplicity was calculated as 6.08, compared to the experimental value of 6.28 (a difference within typical experimental error margins). Neutron multiplicities were also minimally impacted (dropping from 2.86 to 2.80).

Significance and Claims
The paper asserts that these results demonstrate microscopic theory is becoming quantitatively competitive with phenomenological models. By providing a parameter-free prediction of angular momentum distributions that successfully reproduces experimental photon multiplicities, the work validates the predictive power of EDF-based microscopic frameworks. The authors argue that these microscopic distributions can serve as valuable guidance for refining phenomenological models, particularly in regimes where experimental data is limited.

The study acknowledges current limitations, such as the neglect of intrinsic excitations beyond nuclear deformation and specific rotational components ($K$, $0$). However, the successful reproduction of fission spectra without parameter tuning marks a significant step toward a fully microscopic description of fission fragment decay. Future work is identified as focusing on refining angular momentum modeling (e.g., restoring parity symmetry) and improving the microscopic prediction of other critical inputs for statistical reaction theory, such as fission yields and excitation energies.