Building the Bridge from Fundamental Forces to Element Formation
11 May 2026, Yanjiang
Ab initio calculations of beta-decay lifetimes at the N=50 magic number connect fundamental nuclear forces to the origin of heavy elements.
Every atom of gold, silver, or uranium in your body was forged in a cataclysm — a neutron-star merger or a supernova explosion that scattered heavy elements across the galaxy. But exactly how that forging happens remains one of the great open questions in astrophysics. The r-process, or rapid neutron-capture process, is responsible for about half of the elements heavier than iron. And its speed — whether it can produce enough gold to explain the universe’s abundance — depends on a single, poorly measured quantity: the beta-decay lifetime of nuclei so neutron-rich they barely exist.
A preprint (arXiv:2509.06812) by a team led by Achim Schwenk at the Max-Planck-Institut für Kernphysik and Zhen Li at the Technische Universität Darmstadt now brings these elusive lifetimes within reach. Using only known fundamental forces — nothing but the strong nuclear interaction derived from chiral effective field theory — the team has computed beta-decay lifetimes for a family of nuclei at the N=50 magic number, with no phenomenological adjustments. The results agree well with existing experimental data, and they open the door to predicting lifetimes for nuclei that no experiment has ever measured.
The bottleneck in the r-process
Think of the r-process as a nuclear production line. Neutrons are added to seed nuclei one by one, building up heavier and heavier isotopes. But not all nuclei absorb neutrons eagerly. At certain “magic numbers” — where a nuclear shell is filled — the addition of neutrons becomes harder, and the production line slows down. The nucleus must wait for a beta decay to convert a neutron into a proton, raising the atomic number and allowing neutron capture to resume. These bottlenecks are called “waiting points,” and the N=50 magic number is one of the most important.
The N=50 waiting-point nuclei — isotopes of nickel, copper, zinc, gallium, and germanium — sit at a critical junction in the r-process. Their decay lifetimes determine how long the r-process tarries at this shell closure, and that delay controls the final abundance pattern of elements from strontium to silver. Yet until now, the lifetimes of these nuclei were known only from a handful of experiments, and for the most neutron-rich species, even those measurements were missing. Like a map with large blank patches, the absence of reliable lifetimes made it nearly impossible to connect nuclear physics to the observed element abundances.
Building from first principles
What makes this work remarkable is its starting point. Previous calculations of beta-decay lifetimes relied on phenomenological models that adjusted parameters to fit known data. That approach works well for nuclei that have been measured, but its predictive power for unknown nuclei is uncertain. The team behind this preprint took a different route: they began with the fundamental theory of the strong nuclear force — chiral effective field theory — and built everything upward from that foundation.
Chiral EFT provides a systematic expansion of nuclear forces in terms of the interactions between nucleons. It is not a free parameterization; it respects the symmetries of quantum chromodynamics and describes how pions and nucleons interact at low energies. From this foundation, the team then applied a many-body method called the in-medium similarity renormalization group, or VS-IMSRG. This technique transforms the computationally intractable problem of solving the Schrödinger equation for dozens of interacting nucleons into a manageable set of equations for an effective Hamiltonian that acts only within a small “valence space” — the handful of orbitals where the nuclear structure changes.
The key insight is that the same transformation must be applied not only to the nuclear Hamiltonian but also to the weak-interaction operator that drives beta decay. This is far from trivial: evolving the decay operator consistently with the Hamiltonian ensures that the resulting beta-decay rates carry no hidden approximations. It is like building a bridge by first surveying the terrain with the same tools that will later be used to construct the road — the two must be fully consistent for the path to hold.
The missing ingredient
When the team performed their calculations, they made a crucial discovery. The standard treatment of beta decay uses only the “one-body” part of the weak current — the simplest coupling between a single nucleon and the emitted electron and neutrino. But the nuclear force itself involves exchanges of pions and other mesons, which generate additional contributions known as two-body currents. These are subtle effects: they arise when the decay probe is shared between two interacting nucleons, rather than acting on a single nucleon in isolation.
Including these two-body currents changed the computed lifetimes significantly. Without them, the predicted decay times were systematically too short — the nuclei seemed to decay too quickly, in disagreement with the few available measurements. With the two-body currents included, the lifetimes lengthened and aligned closely with experiment. This is not a minor correction; it is essential for any ab initio prediction of beta decay. Acknowledging the limits of this analogy: two-body currents are not simply an extra ingredient thrown into a recipe — they emerge naturally from the same chiral forces that describe nuclear binding. The framework is self-consistent.
For the benchmark nucleus nickel-78 — a doubly magic isotope that sits at the heart of the N=50 waiting point — the agreement is particularly striking. The team’s calculation reproduces both the measured ground-state energy and the low-lying excited states, and the beta-decay lifetime matches the experimental value within the theoretical uncertainties. This validation gives confidence that the same method can be extended to even more neutron-rich nuclei that cannot yet be produced in laboratories.
Implications for nucleosynthesis
Why does this matter for the origin of the elements? The r-process network calculations that simulate element formation require input from hundreds of nuclear lifetimes, many of which will never be measured directly. Until now, these networks had to rely on extrapolations from phenomenological models, introducing unknown uncertainties. The ab initio approach demonstrated in this work offers a path to replace extrapolation with calculation. The team’s results for the N=50 waiting points suggest that earlier models may have underestimated the lifetimes, which would in turn affect the predicted abundance pattern of elements around the first r-process peak.
The road ahead is clear: the same methodology can be applied to other magic numbers — N=82, N=126 — that control the production of elements like gold and platinum. The team is already working on these extensions, and the computational framework is general enough to accommodate further refinements, such as including first-forbidden transitions that become important for very neutron-rich nuclei. The preprint does not claim to have answered every question about the r-process, but it establishes a solid foundation on which future work can build. For a field that has long relied on empirical guesswork, this is a step toward genuine predictive power.
Yanjiang is an online editor of LoomSci
References
- Zhen Li et al., Ab initio calculations of beta-decay half-lives for $N=50$ neutron-rich nuclei, arXiv:2509.06812