Shell effects and the neutron emission within the multi-dimensional Langevin model for fission

This paper presents a multi-dimensional Langevin model that dynamically simulates the time evolution of fissioning nuclei with stochastic neutron emission, enabling the calculation and comparison of pre-scission neutron multiplicities, fragment mass distributions, and neutron energy spectra against experimental data.

The Gist

Imagine a heavy atomic nucleus, like a giant, wobbling drop of water made of protons and neutrons. Sometimes, this drop gets so excited (heated up) that it starts to stretch, elongate, and eventually snaps in two. This is nuclear fission.

For decades, physicists have tried to predict exactly how this snap happens. They use complex math called Langevin equations to simulate the nucleus stretching like taffy. But there's a catch: as the nucleus stretches, it doesn't just sit there. It's also sweating. It's losing energy by shooting out tiny particles called neutrons.

This paper is about a new, more realistic way to simulate this process. The authors, Ivanyuk and his team, decided to stop ignoring the "sweat" (neutron emission) and start tracking it step-by-step as the nucleus stretches.

Here is the breakdown of their work using simple analogies:

1. The Setup: The Stretching Taffy

Think of the nucleus as a piece of hot, sticky taffy being pulled apart by two people.

• The Old Way: Previous simulations would pull the taffy apart and say, "Okay, it broke here. Now let's guess how many bits of sugar (neutrons) fell off." They treated the stretching and the sugar-falling as two separate events.
• The New Way: The authors built a simulation where, every single second the taffy stretches, they check: "Is it hot enough right now to spit out a piece of sugar?"
  • If yes: They immediately remove that piece of sugar, cool the taffy down a bit (because losing sugar takes energy), and then keep pulling.
  • If no: They just keep pulling.

2. The "Sweat" Formula (The Neutron Emission Rate)

To know when the nucleus spits out a neutron, the authors had to invent a new rulebook.

• Imagine the nucleus is a crowded dance floor (the Fermi gas). The dancers are moving fast.
• The "exit door" is the edge of the nucleus.
• The authors used a clever math trick (based on the Continuity Equation) to calculate exactly how many dancers are likely to run out the door at any given moment, based on how hot the dance floor is and how wide the door is.
• This allowed them to calculate the "evaporation rate" of neutrons in real-time, rather than guessing at the end.

3. The Results: What Happened When They Added the "Sweat"?

When they ran their simulation for Uranium-236 (a heavy nucleus), they found some fascinating things:

• The "Cooling" Effect: Every time a neutron leaves, the nucleus gets cooler. This is like a runner taking off their heavy jacket; they can't run as fast anymore. Because the nucleus cools down, it sometimes loses the energy needed to break apart completely.
• The "Shell" Magic: Nuclei have "magic numbers" of particles that make them extra stable (like a perfectly balanced stack of blocks). The authors found that by cooling the nucleus down via neutron emission, these "magic" structures stay strong for longer. This helps explain why the two pieces of the broken nucleus come out in specific sizes, matching real-world experiments much better than before.
• When do they leave? They discovered that neutrons don't just leave at the very end (when the nucleus snaps).
  • At lower temperatures, neutrons mostly leave after the nucleus has already started stretching significantly (past the "hump" or barrier).
  • At higher temperatures, some neutrons leave right from the start, while the nucleus is still sitting comfortably in its original shape.

4. Why Does This Matter?

Think of nuclear fission like a car crash.

• Old Simulations: You could predict the speed of the cars before the crash, but you couldn't predict exactly how the metal would crumple or how much debris would fly off.
• This Paper: By tracking every little piece of debris (neutron) as it flies off during the crash, the authors can now predict the final shape of the wreckage (the mass of the fragments) and the amount of debris with much higher accuracy.

The Bottom Line

The authors successfully combined the "stretching" of the nucleus with the "cooling" caused by shooting out neutrons.

• Before: They treated the nucleus as a static object that cooled down only at the end.
• Now: They treat it as a dynamic, cooling object that changes its shape and energy while it is stretching.

This tiny change in the math made their predictions match real-life experimental data almost perfectly. It's like finally adding the "wind resistance" to a video game car physics engine; suddenly, the car behaves exactly like a real car on a real road. This helps scientists better understand nuclear energy, nuclear waste, and the fundamental forces that hold our universe together.

Technical Summary

1. Problem Statement

The paper addresses the challenge of simultaneously modeling the shape evolution of a fissioning nucleus and the emission of pre-scission neutrons using the Langevin approach.

• Context: While the Langevin equation has been successfully used to describe fission dynamics (mass distributions, kinetic energies) in 1–5 dimensions, previous studies often neglected the dynamic coupling between shape evolution and particle evaporation (neutron emission) during the descent from the compound nucleus to the scission point.
• Limitations of Prior Work: Existing models often treated neutron emission as a post-process or used approximations that did not account for the instantaneous reduction of excitation energy ($E^*$) and the change in the potential energy surface (PES) due to the loss of a neutron. Furthermore, standard emission models often provided total probabilities for the whole process rather than instantaneous rates required for step-by-step integration.
• Goal: To develop a self-consistent framework where neutron emission occurs stochastically during the Langevin trajectory, dynamically altering the system's excitation energy and shell structure, and to analyze how this affects fission observables like fragment mass distributions and neutron multiplicity.

2. Methodology

The authors modified the standard multi-dimensional Langevin approach to incorporate instantaneous neutron emission.

A. The Langevin Framework

• Dimensionality: The model uses a four-dimensional shape parametrization based on the Two-Center Shell Model (TCSM). The variables are:
  1. $z_0/R_0$: Distance between centers of left and right oscillator potentials.
  2. $\delta_1, \delta_2$: Deformations of the left and right parts.
  3. $\alpha$: Mass asymmetry.
  • Constraint: The neck parameter $\epsilon$ is fixed at 0.25 (slightly lower than the standard 0.35) to reduce computational cost while maintaining accuracy for $^{236}\text{U}$.
• Equations of Motion: The time evolution of collective coordinates ($q_\mu$) and momenta ($p_\mu$) is governed by Langevin equations including:
  • Temperature-dependent free energy $F(q, T)$.
  • Macroscopic-microscopic potential energy (Folded Yukawa + Strutinsky shell corrections).
  • Friction (Wall-and-window formula) and Inertia (Werner-Wheeler approximation).
  • Random forces satisfying the fluctuation-dissipation theorem with an effective temperature $T^*$.

B. Neutron Emission Model

The core innovation is the derivation of an instantaneous neutron emission rate ($dN/dt$) integrated into the Langevin steps.

• Derivation: Based on the quantum kinetic equation and the Fermi-gas model. The authors derived an analytical approximation for the flux of neutrons escaping the nucleus through a surface integral.
• Emission Rate Formula:
  $$\frac{dN}{dt} \approx S_{int} \frac{\pi R_0^2 T}{\hbar \pi} \left(\frac{9\pi N}{4}\right)^{2/3} \left(1 + \frac{T}{E_F} - \frac{\mu}{E_F}\right) e^{\mu/T}$$
  Where $S_{int}$ is a deformation-dependent surface integral, $\mu$ is the chemical potential (energy of the last occupied neutron state), $E_F$ is the Fermi energy, and $T$ is the nuclear temperature.
• Stochastic Implementation:
  1. At every integration step ($\Delta t = 0.05$ fm/c), the number of emitted neutrons $\Delta N = \Delta t (dN/dt)$ is calculated.
  2. A random number $r_n \in [0,1]$ is generated. If $\Delta N > r_n$, a neutron is emitted.
  3. Dynamic Update: Upon emission:
     • The local excitation energy is reduced by the neutron separation energy ($S_n = -\mu$) plus the average kinetic energy of the emitted neutron ($\langle E_n \rangle \approx 3/2 T$).
     • The system switches to the PES layer corresponding to the new, lower neutron number ($N-1$).
     • The trajectory continues with the updated energy and potential.
  4. Trajectories that lose too much energy fall below the fission barrier and are classified as evaporation residues (not fission events).

3. Key Contributions

1. Self-Consistent Coupling: The paper presents a method to couple the stochastic Langevin dynamics with the stochastic emission of neutrons, allowing the fission trajectory to dynamically "feel" the loss of excitation energy and changes in shell structure.
2. Analytical Emission Rate: Development of a computationally efficient analytical expression for $dN/dt$ based on the continuity equation and Fermi-gas statistics, suitable for integration into time-dependent simulations.
3. Shell Effect Preservation: The model demonstrates how neutron emission preserves shell effects by reducing excitation energy, which prevents the "washing out" of shell corrections that typically occurs at high excitation energies.
4. Detailed Emission Timing: The ability to map exactly when (at which deformation) neutrons are emitted during the fission process, a quantity inaccessible to experiment.

4. Results

The model was applied to the fission of $^{236}\text{U}$ at excitation energies $E^* = 10, 20, 30,$ and $40$ MeV.

• Fission Fragment Mass Distributions:
  • Including neutron emission significantly improves agreement with experimental data compared to "first-chance" fission (no emission) calculations.
  • The emission of neutrons reduces the excitation energy, thereby enhancing the impact of shell effects on the mass yields. This leads to a better reproduction of the double-humped mass distribution structure observed in experiments.
• Neutron Emission Timing (Deformation Dependence):
  • At low $E^*$ (10 MeV): Neutrons are emitted almost exclusively after the system crosses the fission barrier (between the saddle and scission). Emission from the ground state region is suppressed because the energy loss would trap the system in the potential well (preventing fission).
  • At high $E^*$ (>20 MeV): A small fraction of neutrons is emitted from the ground state region, but the majority are still emitted between the saddle and scission points.
• Pre-Scission Neutron Spectra:
  • The calculated neutron energy spectra follow a Watt distribution, consistent with experimental prompt fission neutron spectra. The average energy is approximately 2 MeV.
• Pre-Scission Neutron Multiplicity ($M_{pre}$):
  • The calculated $M_{pre}$ values for $^{236}\text{U}$ show good agreement with the limited available experimental data (from $\alpha + ^{232}\text{Th}$ and $p + ^{236/238}\text{U}$ reactions).
  • $M_{pre}$ increases with excitation energy, as expected.

5. Significance

• Theoretical Advancement: This work bridges the gap between macroscopic fission dynamics and microscopic particle evaporation. It proves that neglecting the dynamic feedback of neutron emission leads to inaccuracies in predicting fission fragment properties, particularly at lower excitation energies where shell effects are dominant.
• Experimental Insight: The model provides a theoretical explanation for the "missing" pre-scission neutrons in low-energy fission, attributing them to trajectories that fail to reach scission due to energy loss (evaporation residues).
• Predictive Power: The ability to calculate the distribution of emitted neutrons with respect to deformation and energy offers new observables for future experimental verification and improves the reliability of nuclear data for reactor physics and astrophysical applications.
• Validation: The successful reproduction of both mass distributions and neutron multiplicities validates the use of the Langevin approach with dynamic shell corrections and particle emission for describing complex fission dynamics.