Limits of atomic nuclei predicted

Novel calculations have enabled the study of nearly 700 isotopes between helium and iron, showing which nuclei can exist and which cannot. In an article published in Physical Review Letters, scientists report how they simulated for the first time using innovative theoretical methods a large region of the chart of nuclides based on the theory of the strong interaction.

Atomic nuclei are held together by the strong interaction between neutrons and protons. About ten percent of all known nuclei are stable. Starting from these stable isotopes, nuclei become increasingly unstable as neutrons are added or removed, until neutrons can no longer bind to the nucleus and “drip” out. This limit of existence, the so-called neutron “dripline”, has so far been discovered experimentally only for light elements up to neon. Understanding the neutron dripline and the structure of neutron-rich nuclei also plays a key role in the research programme for the future accelerator facility FAIR at the GSI Helmholtz Centre for Heavy Ion Research in Darmstadt.

In a new study, “Ab Initio Limits of Nuclei”, published in the journal Physical Review Letters as an Editors’ Suggestion with an accompanying synopsis in APS Physics, Professor Achim Schwenk of TU Darmstadt and a Max Planck Fellow at the MPI for Nuclear Physics in Heidelberg, together with scientists from the University of Washington, TRIUMF and the University of Mainz, succeeded in calculating the limits of atomic nuclei using innovative theoretical methods up to medium-mass nuclei. The results are a treasure trove of information about possible new isotopes and provide a roadmap for nuclear physicists to verify them.

The new study is not the first attempt to theoretically explore the extremely neutron-rich region of the nuclear landscape. Previous studies used density functional theory to predict bound isotopes between helium and the heavy elements. Professor Schwenk and colleagues, on the other hand, explored the chart of nuclides for the first time based on ab initio nuclear theory. Starting from microscopic two- and three-body interactions, they solved the many-particle Schrödinger equation to simulate the properties of atomic nuclei from helium to iron. They accomplished this by using a new ab initio many-body method – the In-Medium Similarity Renormalization Group –, combined with an extension that can handle partially filled orbitals to reliably determine all nuclei.

Starting from two- and three-nucleon interactions based on the strong interaction, quantum chromodynamics, the researchers calculated the ground-state energies of nearly 700 isotopes. The results are consistent with previous measurements and serve as the basis for determining the location of the neutron and proton driplines. Comparisons with experimental mass measurements and a statistical analysis enabled the determination of theoretical uncertainties for their predictions, such as for the separation energies of nuclei and thus also for the probability that an isotope is bound or does not exist (see figure).

The new study is considered a milestone in understanding how the chart of nuclides and the structure of nuclei emerges from the strong interaction. This is a key question of the DFG-funded Collaborative Research Center 1245 "Nuclei: From Fundamental Interactions to Structure and Stars" at the TU Darmstadt, within which this research was conducted. Next, the scientists want to extend their calculations to heavier elements in order to advance the input for the simulation of the synthesis of heavy elements. This proceeds in neutron-rich environments in the direction of the neutron dripline and occurs in nature when neutron stars merge or in extreme supernovae.


Original publication:

Ab initio limits of atomic nuclei, S. R. Stroberg, J. D. Holt, A. Schwenk, and J. Simonis, Phys. Rev. Lett. 126, 022501 (2021), DOI: 10.1103/PhysRevLett.126.022501

Physics Synopsis: Predicting the Limits of Atomic Nuclei, M. Stephens

Group Schwenk at TU Darmstadt

Interview with Achim Schwenk ("Welt der Physik", in German)


Prof. Achim Schwenk
Phone: +49 6151 16-21550

Global ab initio calculations from helium to iron (element He and Fe). The gray area shows all 700 calculated atomic nuclei, while the color (and height) for each isotope with neutron number N and proton number Z corresponds to the probability that it is bound (probability bound). For comparison, the experimentally known neutron and proton driplines and the most recently discovered isotope are shown with filled and empty drop-symbols, respectively. (chart: TUD)