Impact of neutron captures on the r-process: I. Propagation of statistical model uncertainties

The universe is a vast chemical factory, but for decades, scientists have struggled to explain how roughly half of all elements heavier than iron were forged. These heavy elements, which include everything from gold to uranium, are created in the most violent cosmic events imaginable, such as the collision of two neutron stars. In these extreme environments, atomic nuclei are bombarded by a flood of neutrons, rapidly building up heavier and heavier atoms in a process known as the rapid neutron-capture process. While astronomers have confirmed that these collisions do happen and produce the right kinds of elements, the precise physics governing how the atoms change during the explosion remains largely a mystery. Because we cannot recreate these conditions in a laboratory, researchers must rely on complex computer models to simulate the journey of matter from a simple atom to a heavy element. The accuracy of these simulations depends entirely on the data fed into them, specifically how we understand the internal energy states of unstable, neutron-rich atoms that exist only for a fleeting moment.

A team of researchers led by Yu An at the Technical University of Darmstadt has taken a deep dive into one of the most critical pieces of this puzzle: the statistical models used to predict how these unstable atoms capture neutrons. In their study, they did not try to invent new physics but rather tested how much the final results of a cosmic explosion change when the underlying assumptions about atomic structure are slightly altered. They focused on two specific theoretical ingredients: nuclear level density, which describes how many different ways an atomic nucleus can vibrate or store energy, and the gamma-ray strength function, which measures how likely the nucleus is to emit light energy as it settles down. By running thousands of simulations with different combinations of these theoretical models, the team mapped out exactly where our lack of knowledge creates the biggest gaps in our understanding of the universe's chemical history.

The researchers began by simulating the aftermath of a neutron star merger, a scenario where two dead stars smash together and eject a cloud of super-dense material. They ran this simulation using sixty different combinations of the theoretical models available to them, creating a wide range of possible outcomes for the final mix of elements. The results revealed a clear hierarchy of uncertainty. The models describing how nuclei store energy, known as nuclear level densities, had a far greater impact on the final chemical makeup than the models describing how they emit light. While the major peaks of heavy elements, such as those corresponding to gold and platinum, remained relatively stable regardless of the model used, the regions between these peaks showed dramatic differences. In some simulations, the abundance of elements in a specific mass range around 172 to 178 atomic units would triple or disappear entirely depending on which energy model was chosen.

To understand why this happens, the team looked closely at the timing of events inside the expanding cloud of debris. After the initial explosion, the temperature drops rapidly. At first, the atoms are so hot that neutrons are constantly being captured and then knocked off again in a balanced dance. Eventually, the temperature falls enough that this balance breaks, and the atoms begin to settle into their final forms. The researchers found that the choice of energy model becomes critical in the "cold" phase of this process, which occurs when the temperature drops to about 0.45 billion degrees. In this window, if the theoretical model predicts that the atoms capture neutrons slowly, the material lingers in a specific region of the periodic table, creating a distinct sub-peak of elements. If the model predicts faster capture, that material rushes past to become even heavier elements, leaving a gap behind. This sensitivity is not just a quirk of a single simulation; the team found that even when they averaged the results of over 11,000 different explosion scenarios to mimic a realistic galactic mixture, these specific gaps and peaks persisted.

The study also explored a vast grid of hypothetical conditions, varying the initial density and the ratio of protons to electrons in the material to see how different environments react to these uncertainties. They discovered that the impact of the energy models is not uniform across all cosmic events. In environments where the supply of neutrons is scarce, the choice of model matters very little because the atoms simply run out of fuel before they can be significantly altered. However, in environments with a rich supply of neutrons and specific thermodynamic conditions, the difference between models can determine whether a particular element is produced in abundance or barely at all. The researchers identified that the most significant uncertainties arise in "cold" scenarios where the atoms continue to capture neutrons long after the initial equilibrium has broken down. In these conditions, small changes in the predicted speed of neutron capture can shift the entire flow of matter, deciding whether the process stops at the second peak of heavy elements or pushes through to the third.

Ultimately, this work serves as a crucial map for future experiments. It tells nuclear physicists exactly which regions of the atomic chart are most sensitive to our current theoretical gaps. The study confirms that while we have a good general picture of how heavy elements are made, the fine details of the chemical landscape depend heavily on how well we understand the internal energy states of unstable nuclei. The researchers conclude that to refine our understanding of the universe's chemical evolution, we need better experimental data and more precise theoretical descriptions of nuclear level densities, particularly for the specific isotopes that sit between the major peaks of heavy elements. Without these improvements, our simulations of the cosmos will continue to carry a shadow of uncertainty, unable to fully explain the precise abundance of the heavy elements that make up our world.