⚛️ nuclear experiments

r-Process Nucleosynthesis With Ab Initio Nuclear Masses Around The N=82 Shell Closure

This study demonstrates that ab initio nuclear mass calculations using the valence-space in-medium similarity renormalization group can refine r-process nucleosynthesis predictions around the N=82 shell closure, revealing a strengthened waiting point that slows the nucleosynthesis flow and significantly alters the abundance patterns of heavy elements compared to phenomenological models.

Original authors: Jan Kuske, Takayuki Miyagi, Almudena Arcones, Achim Schwenk

Published 2026-07-03

📖 3 min read🧠 Deep dive

CC BY 4.0

Original authors: Jan Kuske, Takayuki Miyagi, Almudena Arcones, Achim Schwenk

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

Imagine the universe as a giant cosmic kitchen where the ingredients for heavy elements (like gold, platinum, and uranium) are cooked up in extreme events, such as crashing neutron stars. This cooking process is called the r-process (rapid neutron capture).

Here is the problem: To understand how this cosmic cooking works, scientists need a recipe book that lists the "weight" (mass) of every single ingredient (nucleus) involved. However, many of these ingredients are so unstable and rare that we can't weigh them in a lab yet. So, scientists have to guess their weights using mathematical models.

The Old Way vs. The New Way
For a long time, scientists used a "global" model (like a rough sketch) to guess the weights of all these missing ingredients. It worked okay, but it wasn't perfect.

In this paper, the authors tried a new, high-tech approach. Instead of a rough sketch, they used a super-precise method called ab initio (which means "from the beginning"). Think of this as building a 3D model of the ingredients from scratch using the fundamental laws of physics, rather than just guessing based on similar-looking items. They focused specifically on a critical "kitchen zone" around a specific group of atoms known as the N = 82 shell closure (a region near the element Tin).

The Experiment: A Traffic Jam in the Cosmos
The researchers took their new, high-precision weights and plugged them into a simulation of a neutron star crash to see how it changed the cooking process.

Here is what they found, using a simple analogy:

The Traffic Jam: Imagine the r-process as a river of water flowing downstream, carrying atoms with it. As the atoms flow, they hit a "waiting point" (a bottleneck) where they pause before moving on.
The Old Model: The old, rough sketch predicted the water would flow through this bottleneck relatively quickly.
The New Model: When the authors used their new, precise weights, they found that the bottleneck became much stronger. It was like a dam forming in the river. The water (nuclei) piled up and got stuck at a specific spot (around atomic number N = 86) for longer.

The Results: A Slower Cook and a Shifted Peak
Because the "traffic jam" lasted longer:

Slower Flow: The overall cooking process slowed down.
More Pile-Up: More atoms accumulated at that specific bottleneck, making the "second peak" of heavy elements heavier and more abundant.
The Third Peak Shifted: Because the process was delayed, the final batch of heavy elements (the "third peak") ended up slightly shifted toward heavier weights than the old models predicted.

Why This Matters
The authors showed that for the first time, these super-precise, "from-scratch" calculations can actually change our understanding of how the universe makes heavy elements. They proved that even though we can't measure these atoms in a lab yet, our new theoretical tools are accurate enough to refine our cosmic recipes.

They also checked their work against different types of cosmic events (like different types of supernovas) and found that while the effect varies depending on the specific event, the new method consistently provides a clearer picture of where the "traffic jams" happen in the universe's kitchen.

In Short:
The paper demonstrates that using advanced, first-principles physics to calculate the weight of unstable atoms reveals that the universe's "heavy element factory" operates with more traffic jams and delays than we previously thought, leading to a slightly different final product than our old models predicted.

Problem Statement
The origin of heavy elements beyond iron is primarily attributed to the rapid neutron capture process (r-process). However, the extreme neutron-rich conditions required for this process involve nuclei that are currently experimentally inaccessible. Consequently, theoretical predictions are essential but suffer from significant disagreements among global mass models away from stability. These models, such as the Finite-Range Liquid Drop Model (FRDM), are phenomenological and may not accurately capture the nuclear structure effects critical for determining the r-process path, particularly around shell closures like $N=82$ .

Methodology
The authors integrate ab initio nuclear mass calculations into global r-process nucleosynthesis simulations to refine predictions.

Nuclear Mass Calculations: The study employs the valence-space in-medium similarity renormalization group (VS-IMSRG), an ab initio many-body method. Calculations were performed for 70 isotopes ( $^{122-134}$ Pd, $^{123-135}$ Ag, $^{124-137}$ Cd, $^{125-139}$ In, and $^{126-140}$ Sn) around the $N=82$ shell closure. These calculations utilize nucleon-nucleon (NN) and three-nucleon (3N) interactions derived from chiral effective field theory.
Uncertainty Quantification: To assess many-body truncation errors, two different valence spaces were used: $Z=28-50, N=82-126$ (VS1) and $Z=28-50, N=64-90$ (VS2). The differences in one-neutron separation energies ( $S_n$ ) between these spaces for $N=83-86$ were used to define uncertainty bands.
Nucleosynthesis Simulations: The publicly available network code WinNet was used. The baseline model combined the FRDM2012 mass model with experimental data from AME2020. The new VS-IMSRG masses were incorporated into this baseline, creating three scenarios: a central scenario, a minimum scenario, and a maximum scenario (uncertainty band).
Astrophysical Context: Simulations were run for trajectories representing four distinct astrophysical environments: neutron star mergers (NSM), long-time NSM disk ejecta (NSM-DISK), magneto-rotational supernovae (MRSN), and neutron star-black hole mergers (NSBH). Both individual representative trajectories and mass-integrated ejecta from NSM simulations were analyzed.

Key Results

Impact on the Second Peak: The incorporation of VS-IMSRG masses strengthens the "waiting point" at the second r-process peak ( $N=82$ ). Specifically, the generally lower separation energies and increased pairing correlations for Pd, Ag, and Cd isotopes beyond $N=86$ cause matter to accumulate at $N=86$ . This creates a secondary "waiting region."
Slower Nucleosynthesis Flow: Due to the accumulation at $N=86$ , the overall nucleosynthesis flow slows down. This results in a later freeze-out of the r-process compared to the standard FRDM2012 model.
Shift of the Third Peak: The delayed freeze-out allows for additional neutron captures after the main flow halts, leading to a shift of the third r-process peak ( $N=126$ ) toward higher mass numbers ( $A$ ). This effect is most pronounced in the minimum VS-IMSRG scenario.
Environmental Dependence: The impact is most significant in NSM and NSM-DISK trajectories where the r-process path crosses the region of the new mass calculations. In MRSN trajectories, the path remains closer to stability, and in NSBH trajectories, the path moves far beyond the calculated region; in these two cases, the impact of the new masses is minor.
Robustness: The results hold across different baseline mass models (HFB24, HFBD1M) and beta-decay rate sets (REACLIB, D3C*), confirming that the observed shifts are driven by the ab initio mass inputs rather than other nuclear physics inputs.
Validation: The VS-IMSRG mass scenarios are consistent with recent experimental measurements of In and Sn masses not included in the AME2020 dataset, falling within the estimated many-body uncertainties.

Significance and Claims
The paper claims to demonstrate, for the first time, that ab initio mass calculations can be effectively used to refine r-process predictions compared to global, phenomenological mass models. By focusing on the $N=82$ shell closure region, the authors show that state-of-the-art theoretical masses can alter the predicted abundance patterns, specifically by strengthening the second peak and shifting the third peak.

The authors emphasize that while these ab initio calculations are not yet feasible for all nuclei, their integration with global models provides new insights into the origin of heavy elements. Furthermore, the study highlights specific regions (such as the pairing gap enhancement in Pd, Ag, and Cd beyond $N=86$ ) where further experimental investigation is crucial. The work suggests that high-fidelity ab initio models, including theoretical uncertainties, are a necessary step forward for understanding the r-process, especially given that many relevant nuclei will remain beyond experimental reach even for next-generation facilities.

More like this