When Nuclei Tell Time: A New Way to Calculate Beta-Decay Lifetimes
11 May 2026, Yanjiang
A first-principles calculation of beta-decay lifetimes for neutron-rich N=50 nuclei provides a self-consistent clock for the r-process element formation.
Every second, somewhere in the universe, a neutron-rich nucleus decides it has had enough. It emits an electron, a neutrino, and becomes something new. This process — beta-decay — sounds simple. Yet for the most extreme, neutron-heavy nuclei, its timing remains a mystery.
That timer is crucial for understanding how elements form. In the violent environments where half of the heavy elements are made, beta-decay rates control the entire chain of production. Get them wrong, and our story of where gold, platinum, and uranium come from unravels.
Now, a team led by Achim Schwenk at the Max Planck Institute for Nuclear Physics in Heidelberg has done something new. They calculated beta-decay rates from first principles — starting from the fundamental equations of nuclear physics and carrying them through to a final number, without adjustments. Their work appears in a preprint on arXiv (arXiv:2509.06812).
The waiting points
To understand why this matters, we need to look at the r-process — rapid neutron capture. This is how stars forge the heaviest elements. A seed nucleus absorbs neutrons one by one, climbing up the chart of nuclides. Every so often, it hits a “waiting point” where the next neutron capture is slow. There, the nucleus must beta-decay before it can continue climbing.
These waiting points have a special property: their neutron number is a “magic number.” Nuclei with magic numbers are unusually stable. For the r-process, one of the most important magic numbers is N=50 (50 neutrons). Here, the beta-decay rate determines how long the entire process pauses.
Until now, calculating these rates from scratch has been too hard. Physicists relied on models that used experimental data as input, extrapolating into unknown territory. The problem is that for the most neutron-rich N=50 nuclei, there is no data.
Think of the team’s approach as building a clock from its most basic parts, rather than adjusting an existing one. You start with the springs and gears — the fundamental forces between protons and neutrons — and assemble them step by step, never peeking at the finished product to correct your work.
From forces to rates
Their method uses three layers. First, they describe the strong nuclear force using a framework called chiral effective field theory. This treats the force between nucleons as a series of interactions, ordered by importance, much like an expansion that converges to the truth.
Second, they use a method called the in-medium similarity renormalization group. It simplifies the calculations by transforming the problem into a smaller, more manageable space — a “valence space” where only the relevant nuclear states live. The key is that this transformation is consistent. The same method applies both to the nuclear Hamiltonian and to the operator describing the beta-decay process itself.
Third, they calculate the Gamow-Teller transition strength. This measures how likely the beta-decay is to happen, based on how the spins of the protons and neutrons rearrange themselves.
One of the key findings involves what the team calls two-body currents. For decades, most calculations only included one-body currents — treating the decay as a single nucleon changing its identity. But nature is more subtle. Two-body currents capture the effect of nucleons interacting during the decay itself. They are small corrections, like the shimmer of heat on a road, but they matter.
The team found that including these two-body currents systematically increases the calculated decay time. And this brings the theoretical values into good agreement with the limited experimental data that exists.
For the nucleus nickel-78 — a doubly magic waiting point — the result matches the known decay time well. The same method works across a chain of N=50 nuclei, from nickel-78 to germanium-82. The agreement is consistent, not just a lucky hit.
The team also examined first-forbidden transitions — a more exotic decay mode where the spin change is larger. These contribute a small, important component to the total rate. Their calculations show that for these N=50 nuclei, the main contribution comes from allowed Gamow-Teller transitions, with first-forbidden effects adding a minor but measurable correction.
A new kind of prediction
What makes this work different from earlier attempts is simplicity. The approach contains no adjustable parameters. The nuclear forces come from chiral effective field theory. The transformation to the valence space uses a rigorous framework. The weak currents — the currents that drive beta-decay — are derived from the same fundamental theory. Everything connects.
This means the team can make predictions for nuclei that no experiment has reached. And those predictions carry a self-consistency that phenomenological models lack.
The team tested two different nuclear interactions. One, called 1.8/2.0 (EM), reproduces many nuclear properties across the chart. Another, called ΔN²LO_GO (394), uses a different set of parameters. Yet both converge to similar decay times for the N=50 waiting points. This agreement across different force models gives confidence that the predictions are robust.
The team also checked sensitivity to their method. They varied the basis size, the oscillator energy, the truncation of three-body forces. The answer barely changed. This stability is a hallmark of reliable computation.
Why it matters
The r-process produces half of all elements heavier than iron. Its path passes through the N=50 waiting points. The decay rates there determine how much time the process has before a neutrino wind blows the material away. Longer pauses increase the chance that the chain breaks. Faster transitions mean a smooth flow toward the heaviest elements.
Every gold atom in a wedding ring, every platinum atom in a catalytic converter, every uranium atom in a power plant — each one passed through waiting points like these. The time it spent there, determined by beta-decay rates, shaped its journey from a stellar explosion to our world.
This paper does not solve the entire problem. The r-process involves many more waiting points, from N=50, to N=82, to N=126. Each magic number poses its own challenge. The method needs to be extended to heavier systems, where the computational cost grows rapidly. And experimental data from rare isotope facilities, like FRIB in the US or the RIBF in Japan, will test whether these first-principles predictions hold.
But the direction is clear. For the first time, scientists can calculate beta-decay rates from the fundamental laws of physics. Not by fitting data, but by starting with the equations that describe how protons and neutrons interact, and carrying that description through to a prediction.
The team — Zhen Li at TU Darmstadt, Takayuki Miyagi at the University of Tsukuba, and Achim Schwenk — has built a bridge. It connects the abstract world of nuclear forces to the concrete reality of a nucleus waiting to transform. On one side sits chiral effective field theory, a framework that captures the essence of the strong force. On the other sits an observable number: how long a nucleus lives before it becomes something new.
The bridge itself will be tested and refined. But it now stands. And that changes what we can ask.
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