Angular and Kinetic Properties of Scission Neutrons… — Plain-Language Explanation

Imagine a giant, unstable balloon (an atomic nucleus) that suddenly snaps in two. This is nuclear fission. For decades, scientists have known that when this happens, the two resulting pieces (called fragments) fly apart at incredible speeds and spit out tiny particles called neutrons.

For a long time, physicists thought all these neutrons were "evaporated" later, like steam rising from a hot cup of coffee after the water has boiled. They assumed the fragments were fully formed and moving steadily before they started letting go of these neutrons.

However, this new paper suggests that some neutrons are actually "kicked out" right at the moment the balloon snaps. These are called scission neutrons. They are born in the chaotic, split-second chaos of the break, not from a calm, cooling fragment later on.

Here is how the researchers found proof of these "snap-time" neutrons, explained simply:

1. The Super-Computer Simulation

To see what happens during the split, the scientists didn't use a microscope; they used a super-computer to run a movie of the event using a theory called Time-Dependent Density Functional Theory (TDDFT).

Think of this like a high-speed, 3D video game where they simulate the atoms dancing and breaking apart. In previous versions of this "game," the virtual world was too small. The neutrons would hit the edge of the screen before the scientists could figure out exactly how fast they were going or which way they were flying.

In this study, they built a much bigger virtual world (about 3 times larger than before). This gave the neutrons enough room to fly out and settle down so the scientists could measure them accurately without the "walls" of the simulation messing up the data.

2. The "Speed Limit" Discovery

Once they had a clear view, they looked at the neutrons flying out at specific angles (mostly sideways and slightly backward relative to the split). They found something surprising:

The "No-Go" Zone: There are no scission neutrons with low energy (below about 1.5 to 2 million electron volts). It's as if there is a speed limit; nothing slow is allowed to be a "snap-time" neutron.
The High-Speed Crowd: Instead, these neutrons are all fast. They cluster around a specific high speed (3–3.5 MeV) and then trail off into a long tail of even faster particles.

It's like a crowd of people jumping off a diving board. The "evaporated" neutrons are like people casually stepping off the pool deck later. The "scission" neutrons are like people who are violently flung off the board the exact second it breaks. The ones flung off the board are always moving fast; you never see a slow one from that specific event.

3. Solving the Mystery of the "Missing" Energy

Scientists have been trying to match their computer models with real-world experiments for years. They had a problem:

The Old Model: If you only count the "steam" (evaporated neutrons), your computer model predicts too few high-energy neutrons. It's like trying to fill a bucket with a small cup, but the bucket keeps needing more water than the cup can provide.
The New Model: When the researchers added the "snap-time" neutrons (the ones they found in their big simulation) to the "steam" neutrons, the math finally worked. The combined model perfectly matched the high-energy data measured in real experiments for Uranium and Californium.

4. Why This Matters

This is a big deal because it's the first time a purely microscopic theory (one that doesn't just guess or assume things exist) has predicted these "snap-time" neutrons and proven they are real.

Before: Scientists had to guess that these neutrons existed because the math didn't add up.
Now: The computer simulation naturally produced them without being told to. It's like predicting a storm by watching the clouds move, rather than just assuming a storm is coming because the weather report says so.

The Bottom Line

The paper concludes that when an atom splits, a small but important chunk of the neutrons (about 6% to 10% of the total) are born in the violent moment of the break. These neutrons are distinct because they are always fast and never slow within certain angles.

By finding this "fingerprint" in the data, the researchers have finally separated the "snap-time" neutrons from the "steam" neutrons, giving us a clearer, more accurate picture of how nuclear fission actually works. This helps refine our understanding of the fundamental forces that hold matter together and tear it apart.

Technical Summary: Angular and Kinetic Properties of Scission Neutrons within Time-dependent Density Functional Theory

Problem Statement
Despite the discovery of nuclear fission nearly a century ago, a complete microscopic description of the process remains elusive, particularly regarding the transition from the outer saddle point to scission and the subsequent formation of primary fission fragments (FFs). Current statistical evaporation codes assume all prompt neutrons are emitted isotropically from fully accelerated fragments. However, this framework cannot account for neutrons emitted during the highly non-equilibrium dynamical process of neck rupture, known as scission neutrons (SNs). While SNs have been conjectured since 1939 and studied via various models, they remain experimentally unresolved because prompt-neutron observables alone cannot distinguish between neutrons emitted at scission and those evaporated later. Previous attempts to predict SNs using time-dependent density functional theory (TDDFT) were limited by small simulation domains, preventing the simultaneous determination of angular and kinetic properties before neutrons reached computational boundaries.

Methodology
This study employs a fully microscopic framework based on TDDFT extended to superfluid systems to investigate SN emission in $^{235}\text{U}(n_{th}, f)$ , $^{239}\text{Pu}(n_{th}, f)$ , and $^{252}\text{Cf}(sf)$ .

Simulation Setup: The authors utilized the axially-symmetric harmonic-oscillator solver AxialHOHFB with the SkM* energy density functional (EDF). A critical advancement was the use of a substantially larger simulation domain (an ellipsoid with semi-axes $Z_{max} = 70$ fm and $R_{max} = 55$ fm), approximately 3.2 times larger than in previous studies.
Two-Step Integration: To manage computational constraints while maintaining accuracy, a two-step time integration procedure was employed. The system was evolved from the static constrained Hartree-Fock-Bogoliubov (HFB) solution to the scission moment in a smaller basis ( $B_1$ ), then transformed to the larger basis ( $B_2$ ) to continue evolution until neutrons approached the boundary.
Extraction of SN Properties: SNs were identified by defining a region $V^{(SN)}_\eta$ where the ratio of collective-flow energy density to intrinsic kinetic energy density exceeded a threshold $\eta = 3$ . This criterion isolates neutrons dominated by outward collective flow.
Angular and Kinetic Analysis: Within this region, local velocity fields and collective-flow energy per neutron were calculated. The study focused on the angular interval $110^\circ \leq \theta \leq 150^\circ$ (measured relative to the light FF direction), where the SN yield saturates before boundary reflections occur.
Comparison with Experiment: The calculated SN kinetic energy spectrum was combined with a Maxwellian evaporation model for the evaporated neutron (EN) component. The evaporation model parameters (fragment temperature $T$ and multiplicity ratio $\bar{\nu}_L/\bar{\nu}_H$ ) were constrained exclusively by low-energy experimental data ( $E \leq E_{thr}$ ). This combined model was then compared against high-energy prompt fission neutron spectrum (PFNS) data for $^{239}\text{Pu}$ and $^{252}\text{Cf}$ .

Key Results

SN Spectrum Characteristics: The extracted SN spectra for all three systems exhibit a distinct threshold behavior. SNs are absent below a specific energy threshold ( $E_{thr}$ ), which is approximately 1.5–2 MeV (specifically $\approx 2$ MeV for induced fission and $\approx 1.5$ MeV for spontaneous fission). The spectra peak near 3–3.5 MeV and feature an exponential high-energy tail extending to roughly 18 MeV.
Yield Estimates: For $^{239}\text{Pu}(n_{th}, f)$ , the total SN yield is calculated as $0.19 \pm 0.01$ (approx. 6.6% of total prompt neutrons). For $^{252}\text{Cf}(sf)$ , the yield is $0.39 \pm 0.02$ (approx. 10.4%).
Model Validation: When the calculated SN component is added to the evaporation model (constrained only by low-energy data), the resulting total PFNS accurately reproduces the measured high-energy yield for both $^{239}\text{Pu}$ and $^{252}\text{Cf}$ . In contrast, the evaporation-only model systematically underestimates the high-energy yield, even within 95% confidence intervals.
Angular Stability: The study confirms that reliable asymptotic properties can only be extracted within the angular range $110^\circ \leq \theta \leq 150^\circ$ . Outside this range, transient interference effects and boundary reflections prevent the yield from saturating within the current simulation domain.

Significance and Claims
The paper claims to provide the first identification of a scission-neutron signature within existing measured prompt fission neutron spectra using a fully microscopic framework. By demonstrating that TDDFT naturally produces SNs without a priori assumptions, and that these SNs resolve the discrepancy between evaporation-only models and high-energy experimental data, the study offers direct evidence for a non-negligible scission-neutron component in prompt fission.

The authors posit that this work establishes a microscopic foundation for interpreting fission neutron spectra beyond the standard statistical-evaporation picture. Furthermore, they identify a specific experimental signature: the separation of SN and EN contributions below the energy threshold ( $E_{thr}$ ). This implies that low-energy measurements in the specified angular range provide the most robust constraints for the evaporation component, while the high-energy tail serves as the primary indicator of the SN contribution. The study concludes that while current simulations are limited to specific angular ranges, the methodology successfully isolates a measurable signature of scission neutrons that has been theoretically conjectured for nearly a century.

Angular and Kinetic Properties of Scission Neutrons within Time-dependent Density Functional Theory