⚛️ nuclear theory

Ab Initio Emergence and Collapse of Nuclear Collectivity near N = Z = 40

This paper presents an \textit{ab initio} study using chiral forces and the in-medium similarity renormalization group to successfully describe the emergence and subsequent collapse of enhanced quadrupole collectivity in nuclei near $N=Z=40$ , specifically capturing the strong deformation in $^{80}$ Zr and its reduction in heavier isotopes without empirical effective charges.

Original authors: X. C. Cao, C. F. Jiao, R. Z. Hu

Published 2026-09-15

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: X. C. Cao, C. F. Jiao, R. Z. Hu

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

At the heart of every atom lies a nucleus, a dense cluster of protons and neutrons held together by the strongest force in nature. For decades, physicists have sought to understand how these tiny particles organize themselves. Sometimes they sit in neat, rigid shells, like electrons around an atom. Other times, they move together in a fluid, collective motion, causing the entire nucleus to stretch and squash into a football shape. This phenomenon, known as nuclear collectivity, is the key to understanding why some atoms are stable and others are not. The challenge has always been to explain this behavior from the ground up, starting only with the fundamental forces between individual particles, without relying on shortcuts or guessed parameters. A particularly puzzling region for scientists is found in the middle of the periodic table, where the number of protons equals the number of neutrons. Here, a specific group of atoms suddenly becomes incredibly deformed, only to lose that shape just a few steps away, a transition that has long resisted a complete explanation from first principles.

A team of researchers has now taken a major step toward solving this mystery by simulating the behavior of these nuclei using only the fundamental laws of physics. They focused on a chain of atoms ranging from krypton to ruthenium, specifically looking at the ones where the number of protons matches the number of neutrons. Their goal was to see if they could reproduce the dramatic rise and fall of nuclear deformation purely through calculation, without tweaking the math to fit experimental data. Using a powerful computational framework that combines two advanced methods, they modeled the interactions between protons and neutrons based on a modern theory of nuclear forces. This approach allowed them to track how the internal structure of the nucleus changes as they added or removed particles, revealing exactly why some nuclei become highly collective while others remain stiff.

The results paint a clear picture of what happens inside these atomic cores. In the middle of the chain, specifically in the nucleus of zirconium-80, the researchers found a strong, football-shaped deformation that matched experimental measurements almost perfectly. The simulation predicted the energy levels and the strength of the nuclear transitions with high precision, confirming that the nucleus is indeed highly collective. As they moved toward neighboring atoms like molybdenum and ruthenium, the simulation showed this collective strength fading away rapidly, just as experiments have observed. Crucially, this entire sequence of events emerged naturally from the calculation. The researchers did not need to introduce any empirical adjustments or guesswork to make the numbers work; the behavior arose directly from the underlying forces between the particles.

To understand why this happens, the team looked closely at the energy gaps between the different orbitals where protons and neutrons reside. They discovered that the key to the strong deformation in zirconium-80 is not a collapse of the shell structure, as some had hoped, but rather a specific narrowing of the energy gap between two particular sets of orbits. This narrowing allows protons and neutrons to mix in a way that creates a strong, cooperative pull, stretching the nucleus. The researchers found that this effect is driven primarily by the interaction between protons and neutrons, which acts to lower the energy of these specific orbits. When they artificially blocked the ability of particles to occupy a specific orbit called 2d5/2, the strong deformation vanished instantly, and the nucleus became nearly spherical. This test confirmed that the presence of particles in this specific orbit is essential for the collective behavior to exist.

The study also clarified why the collectivity disappears as you move toward heavier isotopes like molybdenum-86 and ruthenium-88. In these nuclei, the energy gap between the critical orbits widens again, and the number of particles available to participate in the collective motion decreases. Without this favorable arrangement, the nucleus loses its ability to deform and returns to a more rigid, spherical shape. The researchers showed that this entire process—the sudden onset of deformation and its subsequent collapse—is rooted in the evolution of the nuclear shell structure itself, governed by the fundamental forces between protons and neutrons. By successfully describing this complex behavior without relying on phenomenological models, the work provides a robust, microscopic foundation for understanding heavy, deformed nuclei. It demonstrates that the intricate dance of nuclear collectivity can be traced back to the simple, underlying interactions of its constituent parts, marking a steady advance in the ability to predict the properties of matter from first principles.

Technical Summary: Ab Initio Emergence and Collapse of Nuclear Collectivity near N=Z=40

Problem Statement
The paper addresses a central challenge in nuclear structure theory: explaining how nuclear collectivity emerges from underlying nuclear forces, specifically within self-conjugate ( $N=Z$ ) nuclei near mass number $A=80$ . Experimental data indicates a maximum in quadrupole collectivity (large $B(E2; 2^+_1 \to 0^+_1)$ strengths) in $^{76}\text{Sr}$ and $^{80}\text{Zr}$ , followed by a rapid collapse in $^{86}\text{Mo}$ and $^{88}\text{Ru}$ . While previous studies using energy density functionals (EDF) and large-scale shell models (LSSM) have described these trends, they rely on phenomenological forces, density functionals, or empirical effective charges. Furthermore, existing ab initio approaches, such as angular-momentum-projected coupled-cluster calculations, have struggled to simultaneously describe rotational spectra and transition strengths in this heavy, open-shell region. The authors aim to determine if the observed enhancement and subsequent weakening of collectivity can be reproduced within a unified ab initio framework based solely on chiral nuclear forces.

Methodology
The authors employ a unified framework termed VS-IM-GCM, which combines:

Valence-Space In-Medium Similarity Renormalization Group (VS-IMSRG): They derive nucleus-dependent multi-shell effective Hamiltonians from a common chiral nucleon-nucleon (NN) and three-nucleon (3N) interaction (specifically the 1.8/2.0(EM) interaction). The VS-IMSRG unitary transformation decouples a chosen core ( $^{56}\text{Ni}$ ) and active valence space (spanning orbits from $2p_{3/2}$ to $1h_{11/2}$ ) from excitations outside this space. This is performed non-perturbatively using the Magnus formulation, retaining up to two-body operators (IMSRG(2) truncation).
Projected Generator-Coordinate Method (PGCM): The effective valence-space Hamiltonian is solved using PGCM as a variational approximation. Intrinsic Bogoliubov vacua are generated via Variation After Particle-Number Projection (VAPNP) constrained to quadrupole coordinates ( $Q_{20}, Q_{22}$ ).
Consistent Evolution of Operators: Electric quadrupole ( $E2$ ) transition operators are evolved consistently with the Hamiltonian using the same Magnus transformation. Crucially, no empirical effective charges are introduced; the renormalization of the $E2$ operator from configurations outside the valence space is included via induced one- and two-body terms.
Symmetry Restoration: The intrinsic states are projected onto good angular momentum ( $J$ ) and particle numbers ( $N, Z$ ) to restore symmetries broken in the mean-field approximation.

Key Results

Reproduction of Collectivity Trends: The VS-IM-GCM framework successfully reproduces the emergence of strong quadrupole collectivity in $^{76}\text{Sr}$ and $^{80}\text{Zr}$ and its rapid reduction toward $^{86}\text{Mo}$ and $^{88}\text{Ru}$ . The calculated $B(E2; 2^+_1 \to 0^+_1)$ values for $^{80}\text{Zr}$ ( $2137 \, e^2\text{fm}^4$ ) agree well with experimental data ( $1910(180) \, e^2\text{fm}^4$ ).
Structural Evolution in $^{80}\text{Zr}$ : The calculation reveals that while the VAPNP energy surface exhibits coexisting spherical, triaxial, and prolate minima, symmetry restoration (particle-number and angular-momentum projection) favors the prolate minimum. The low-lying rotational band is built on a single, strongly prolate intrinsic structure rather than a broad mixture of shapes.
Mechanism of Collectivity: The study identifies that the $N, Z=40$ effective single-particle energy (ESPE) gap remains open and does not track the $E2$ strength systematics. Instead, the collectivity is driven by the reduction of the $1g_{9/2}–2d_{5/2}$ spacing in $^{76}\text{Sr}$ and $^{80}\text{Zr}$ . This reduction is governed primarily by the proton-neutron monopole interaction.
Orbital Occupancies: The reduced spacing facilitates the simultaneous occupation of the $1g_{9/2}$ and $2d_{5/2}$ orbits (quasi-SU(3) particles) and the creation of holes in the $pf$ shell (pseudo-SU(3) holes). This configuration allows for strong quadrupole correlations. In $^{86}\text{Mo}$ and $^{88}\text{Ru}$ , the $1g_{9/2}–2d_{5/2}$ spacing widens, and $pf$-shell hole occupancy decreases, leading to the collapse of collectivity.
Role of the $2d_{5/2}$ Orbit: A critical test involved constraining the proton and neutron occupations of the $2d_{5/2}$ orbit to zero in $^{80}\text{Zr}$ . This restriction caused the collective distribution to become nearly spherical and suppressed the $B(E2)$ strength from $2137$ to $13.8 \, e^2\text{fm}^4$ , confirming the essential role of the $2d_{5/2}$ orbit in generating deformation.
Isospin Symmetry: Similar shell evolution and occupation deviations were observed for neutrons, indicating that the mechanism is largely isospin symmetric.

Significance and Claims
The paper claims that these results mark a "steady step toward ab initio computations of heavy open-shell nuclei." By starting from the same chiral NN+3N interaction, the authors demonstrate that the abrupt onset and weakening of quadrupole collectivity near $N=Z=40$ can be reproduced without phenomenological adjustments or empirical effective charges. The work connects the observed structural changes directly to shell evolution rooted in microscopic nuclear forces derived from chiral effective field theory. It validates the VS-IM-GCM framework as a viable approach for describing deformed heavy nuclei characterized by strong collectivity, bridging the gap between microscopic forces and emergent collective phenomena in the $A \approx 80$ region.

More like this