⚛️ nuclear theory

The effect of the coupling between the neutron multiplicity and the fission competition on the evaporation residue cross sections of the $^{48}$ Ca $+^{244}$ Pu reaction

This study demonstrates that a coupled treatment of neutron multiplicity and fission competition, as opposed to a conventional factorized approach, significantly improves the accuracy of survival probability predictions for the $^{48}$ Ca+ $^{244}$ Pu reaction, particularly for the 3n channel, while also highlighting the extreme sensitivity of evaporation residue cross sections to small variations in the fission barrier height.

Original authors: Bakhodir Kayumov, Sarvar Amirov

Published 2026-10-09

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Bakhodir Kayumov, Sarvar Amirov

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

In the quest to understand the very limits of matter, nuclear physicists are constantly trying to build the heaviest possible atoms. These superheavy elements do not exist in nature; they must be created in a laboratory by smashing two lighter atomic nuclei together at incredible speeds. When they collide, they sometimes fuse into a single, massive, and highly excited nucleus. This new creation is unstable, like a spinning top wobbling on a table, and it must quickly cool down to survive. It does this by shooting out tiny particles, specifically neutrons, to shed its excess energy. If it manages to cool down without splitting apart, it becomes a new, heavier element. However, the odds of this happening are incredibly slim. The nucleus is constantly fighting a battle between two forces: the desire to cool down by ejecting neutrons, and the overwhelming urge to split in two, a process called fission. The difference between creating a new element and watching it vanish in a split second often comes down to a fraction of a second and a tiny amount of energy.

A team of researchers in Uzbekistan has taken a fresh look at how this cooling battle plays out, specifically for a reaction that creates an isotope of the element flerovium. They focused on a specific collision: a beam of calcium-48 atoms striking a target of plutonium-244. This reaction is one of the most reliable ways scientists currently have to make superheavy nuclei, yet predicting exactly how many of these new atoms will be created has remained a difficult challenge. The researchers found that the standard way scientists have calculated these odds for decades was missing a crucial piece of the puzzle. By changing how they modeled the cooling process, they discovered that the old methods were significantly underestimating the chances of survival for certain outcomes, revealing that the relationship between how many neutrons are ejected and the risk of the nucleus splitting is far more intertwined than previously thought.

The study centers on a specific moment in the life of a superheavy nucleus. After the initial collision fuses the calcium and plutonium, the resulting nucleus is hot and spinning. It must shed energy to survive. It does this by emitting neutrons one by one. At every single step of this cooling process, the nucleus faces a choice: it can emit another neutron and live to see the next step, or it can split apart and die. In the past, scientists calculated the probability of this survival by treating these two choices as separate events. They would first estimate how many neutrons the nucleus would likely throw away on average, and then, separately, calculate the chance of it splitting at each step along that average path. It was a method that assumed the nucleus followed a single, predictable track.

The new research challenges this separation. The authors argue that the nucleus does not follow a single average track. Instead, every nucleus in the experiment follows its own unique path, losing different amounts of energy at different times. Some nuclei lose energy quickly and become cold early; others stay hot for longer. The risk of splitting is highest when the nucleus is hot. Therefore, the nuclei that survive long enough to emit several neutrons are not the "average" hot nuclei; they are the ones that happened to stay cooler than the average. By treating the number of neutrons and the risk of splitting as independent, the old method missed this subtle filtering effect. It assumed the survivors were just as hot as the average, when in reality, the survivors were the lucky, cooler ones.

To test this idea, the researchers built a new computer model that tracks the entire population of nuclei, not just an average. They simulated the collision and followed thousands of individual nuclei as they cooled down, allowing each one to make its own decisions about whether to emit a neutron or split. This "coupled" approach meant that the number of neutrons emitted and the competition with fission were determined simultaneously for each nucleus based on its actual energy at that moment. The results were striking. For the specific reaction they studied, the new method showed that the chance of creating a specific type of flerovium atom (formed by emitting three neutrons) was much higher than the old method predicted. At certain energy levels, the new calculation showed a survival probability that was seventeen to twenty-four times larger than the traditional estimate.

This massive difference did not mean that the total number of atoms created was wrong, but rather that the distribution of different types of atoms was. The old method was not just off by a little bit; it was misjudging the relative strength of different outcomes. It underestimated the likelihood of the three-neutron outcome while having a much smaller effect on the four-neutron outcome. This suggests that the coupling between cooling and splitting is the key to understanding why some reactions produce more of one isotope than another. The researchers also found that the exact height of the energy barrier that holds the nucleus together is critical. They tested two different ways of calculating this barrier, which differed by less than one million electron volts—a tiny amount in the world of nuclear physics. Yet, this tiny difference changed the predicted survival probability by more than ten times. This highlights that while the new method of tracking the cooling process is a major improvement, the biggest remaining uncertainty lies in knowing the precise strength of the forces holding these superheavy nuclei together.

The team compared their new calculations against real experimental data collected over the last two decades from laboratories in Russia, Germany, and the United States. Their new model, which accounts for the linked nature of cooling and splitting, successfully reproduced eight out of the fourteen measured data points within a factor of two. This is a significant improvement in accuracy for such a complex system. While the model still struggled to perfectly match the data on the rising edge of the energy curve and at very high energies, it provided a much clearer picture of the process than previous attempts. The work demonstrates that to truly predict the creation of superheavy elements, scientists must stop treating the cooling process as a simple, linear average and start accounting for the complex, individual journeys of every single nucleus.

Ultimately, this study does not just refine a number; it changes the understanding of the mechanism itself. It shows that the survival of a superheavy nucleus is not a matter of simple statistics applied to an average path, but a complex filtering process where the survivors are a specific, cooler subset of the original population. The researchers conclude that while their new method is a powerful tool, the path forward requires even more precise knowledge of the nuclear forces at play. Until the exact strength of the fission barrier is known with greater certainty, predictions of how many new elements can be created will remain limited, no matter how sophisticated the computer model becomes. The work stands as a reminder that in the subatomic world, the details of how things happen are just as important as the final result.

Technical Summary: Coupling Neutron Multiplicity and Fission Competition in Superheavy Nucleus Synthesis

Problem Statement
The synthesis of superheavy elements (SHEs) relies on the accurate prediction of evaporation-residue (ER) cross sections, which are determined by the product of three factors: the capture cross section, the probability of compound nucleus (CN) formation ( $P_{CN}$ ), and the survival probability ( $W_{sur}$ ) of the excited CN against fission during neutron emission. While entrance-channel dynamics are often treated with established models, the de-excitation cascade presents significant challenges. Specifically, the conventional statistical treatment of $W_{sur}$ often employs a "factorized" approach. In this method, the realization probability of a specific $xn$ channel is calculated independently of the neutron–fission competition, assuming a mean excitation-energy trajectory. The authors argue that this factorization is an approximation that fails to account for the correlation between neutron multiplicity and fission survival. Because fission preferentially removes nuclei with high excitation energy, the surviving population is "colder" than the mean trajectory suggests. This discrepancy is critical because the fission width depends exponentially on excitation energy; thus, treating the competition as independent of the actual energy distribution of the surviving nuclei may lead to substantial errors in predicting ER cross sections, particularly for the 3n and 4n channels in reactions like $^{48}\text{Ca} + ^{244}\text{Pu}$ .

Methodology
The study investigates the $^{48}\text{Ca} + ^{244}\text{Pu} \to ^{292}\text{Fl}$ reaction, focusing on the de-excitation of the $^{292}\text{Fl}$ compound nucleus. The methodology integrates two distinct frameworks:

Entrance Channel: The capture and fusion probabilities are calculated using the Dinuclear System (DNS) model. This includes solving classical dynamical equations for relative motion with friction and tangential friction terms, accounting for the deformation of the $^{244}\text{Pu}$ target via orientation averaging. The fusion probability $P_{CN}$ is derived from the competition between fusion, quasifission, and drift to symmetry.
De-excitation Cascade: The neutron–fission competition is treated using the Vandenbosch–Huizenga formalism, incorporating Kramers and Strutinsky corrections for fission widths. Level densities are based on the FRDM(2012) tables with Ignatyuk energy dependence.
Coupled vs. Factorized Treatment: The core methodological innovation is the comparison of two treatments for the cascade:
- Factorized Approach: The standard method where the branching ratios are evaluated along a single mean excitation-energy trajectory, and the neutron multiplicity is determined independently.
- Coupled Approach: A rigorous statistical propagation where the full population of the cascade is tracked over the excitation energy. At each step, the probability of neutron emission versus fission is evaluated based on the actual excitation energy of the specific nuclei in the population, rather than a mean value. This accounts for the fact that fission removes the "hot" tail of the distribution, leaving a colder population to emit subsequent neutrons.
Fission Barrier Prescriptions: Two prescriptions for the cold fission barrier ( $B_f^{(0)}$ ) are compared: the tabulated microscopic–macroscopic barrier from Möller et al. and a prescription based on the Myers–Świȩtecki topographical theorem (macroscopic barrier plus ground-state shell correction).

Key Contributions

Implementation of Coupled Cascade: The paper provides a reproducible statistical treatment that removes the factorization assumption, propagating the full excitation-energy population through the evaporation–fission cascade.
Quantification of Factorization Error: The study demonstrates that the factorized approach significantly underestimates survival probabilities for channels lying on the high-energy falling flank of the excitation function.
Sensitivity Analysis: The work isolates the impact of the fission barrier prescription and the cascade treatment, showing that small differences in barrier height (<1 MeV) lead to order-of-magnitude changes in predicted cross sections.
Benchmarking: The method is applied to the $^{48}\text{Ca} + ^{244}\text{Pu}$ reaction, comparing results against a comprehensive set of experimental data from FLNR, GSI, and LBNL.

Results

Impact of Coupling: The coupled treatment yields survival probabilities ( $W_{sur}$ ) for the 3n channel that are 17 times larger at $E^* = 39$ MeV and 24 times larger at $E^* = 47$ MeV compared to the factorized approach. For the 4n channel, the coupled result is raised by factors of 1.7 and 7.5 at the same two energies, respectively. The coupling primarily affects the relative strengths and shapes of neighboring channels (e.g., the ratio $\sigma_{4n}/\sigma_{3n}$ ) rather than the overall normalization.
Fission Barrier Sensitivity: A difference of less than 1 MeV between the two barrier prescriptions (Möller table vs. topographical theorem) results in a change in $W_{sur}$ by more than an order of magnitude (factors of 14 to 24). This confirms that the fission barrier is the dominant source of uncertainty in quantitative predictions.
Comparison with Experiment: Using the coupled cascade and the Möller barrier, the calculated excitation functions reproduce eight of the fourteen measured cross sections within a factor of two. The geometric mean of the calculated-to-experimental ratio is 1.75.
- The 4n channel is reproduced on average (mean ratio 1.04).
- The 3n channel is overestimated by about a factor of three (mean ratio 2.94).
- Significant discrepancies remain on the rising flank (e.g., overestimation by a factor of 29 at $E^* = 34.8$ MeV for 3n) and in the high-energy tail (underestimation by a factor of 7 at $E^* = 53$ MeV for 4n).
- The model fails to satisfy the non-observation limit for the 4n channel at $E^* = 35.2$ MeV, exceeding it by a factor of 19.

Significance and Claims
The paper concludes that the coupling between neutron multiplicity and fission competition is essential for accurately modeling the survival of superheavy nuclei. The factorized treatment, while computationally simpler, is shown to be a poor approximation for these systems, particularly for channels where the mean trajectory does not accurately represent the surviving population.

The authors emphasize that the accuracy of ER cross-section predictions is currently limited more by the uncertainty in the fission barrier (which must be known to within a few hundred keV) than by the statistical model itself. Even with a rigorous coupled treatment, the inability to precisely determine the fission barrier prevents the prediction of absolute cross sections with high precision. The study establishes that the fission barrier is a major source of uncertainty, capable of altering predictions by more than an order of magnitude, and that the coupled treatment is a necessary step toward reducing theoretical uncertainties in the synthesis of superheavy elements.

More like this