⚛️ nuclear theory

Projected shell model description of nuclear level density: Collective, pair-breaking, and multiquasiparticle regimes in even-even nuclei

This paper proposes a novel projected shell model method that utilizes multi-quasiparticle configurations to successfully describe the nuclear level density in deformed nuclei like $^{164}$ Dy by characterizing its evolution through distinct collective, pair-breaking, and multi-quasiparticle regimes.

Original authors: Jiaqi Wang, Saumi Dutta, Long-Jun Wang, Yang Sun

Published 2026-08-13

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jiaqi Wang, Saumi Dutta, Long-Jun Wang, Yang Sun

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

The Atomic Dance Floor: Counting the Steps

Imagine the universe as a giant construction site, and the atoms as the bricks. Inside every atom's core, the nucleus, there is a chaotic, high-energy dance floor packed with tiny particles called protons and neutrons. These particles aren't just standing still; they are constantly jostling, spinning, and pairing up. When a nucleus gets excited—perhaps because it's been hit by a cosmic ray or is part of a star's furnace—it starts to vibrate and spin faster, creating a whole new set of "dance moves" or energy levels.

Scientists call the number of these possible dance moves at a specific energy level the "nuclear level density." Think of it like counting how many different ways a crowd of people can arrange themselves on a dance floor as the music gets louder and faster. This isn't just a game of counting; knowing exactly how crowded the dance floor gets is crucial for understanding how stars burn their fuel, how heavy elements like gold are forged in cosmic explosions, and even how we can safely manage nuclear waste. For decades, scientists have tried to predict this density using simple math, assuming the particles behave like a gas. But as we'll see, the nucleus is more like a complex, choreographed ballet than a simple gas, and the rules change depending on how much energy is in the room.

The Projected Shell Model: A New Way to Count the Dancers

In this paper, a team of physicists introduces a fresh, sophisticated way to count these nuclear dance moves, focusing on a specific type of nucleus called "deformed" nuclei (which are shaped like rugby balls rather than perfect spheres). They used a method called the Projected Shell Model (PSM) to simulate the nucleus of Dysprosium-164 (164Dy), a heavy, well-studied atom. Instead of guessing the answer with a simple formula, they built a digital model where the particles are treated as "quasiparticles"—excited dancers that can break away from their usual pairs.

The researchers found that the story of nuclear level density isn't a single, smooth curve. Instead, they discovered it unfolds in three distinct chapters, or "regimes," as the energy increases:

The Collective Regime (The Warm-Up): At low energies (below 1.5 MeV), the nucleus moves as a whole. The dancers are all paired up and moving in perfect unison, like a synchronized swimming team. The level density here is low and easy to predict because there are only a few specific, well-known patterns they can follow.
The Pair-Breaking Regime (The Soloists Emerge): As the energy rises between 1.5 and 4.0 MeV, the music gets too intense, and the pairs start to break. This is where the magic happens. The authors suggest that the density doesn't just rise smoothly; it takes a "step." First, one pair of dancers (a neutron pair) breaks, causing a jump in the number of possible states. Then, as the energy climbs further, a second pair (a proton pair) breaks simultaneously. This creates a second, distinct step in the curve. The paper argues that previous methods missed this second step because they assumed the density would just start growing exponentially immediately after the first pair broke. The authors' simulations suggest that the "exponential growth" actually gets delayed until both types of pairs have had a chance to break.
The Multi-Quasiparticle Regime (The Chaotic Party): Once the energy goes above 4.0 MeV, almost all the pairs are broken. The nucleus enters a state of "chaos" where the individual dancers are running wild, and the specific choreography is lost. Here, the level density finally starts to follow the smooth, exponential curve that older theories predicted, looking like a bell-shaped distribution of spins.

What the Paper Rules Out and What It Suggests

The authors are careful to point out what their work changes. They explicitly argue against the idea that the nuclear level density starts growing exponentially the moment the very first pair of particles breaks. Their simulations suggest that this "exponential behavior" is actually delayed. They propose that there is a hidden "step" in the curve around 2.0 to 4.0 MeV, caused by the simultaneous breaking of both neutron and proton pairs. This step is a structural feature of the nucleus that simple statistical models miss.

Furthermore, the paper challenges the common assumption that the number of "even" and "odd" parity states (a quantum property related to how the wave function flips) should be equal. In their simulations of 164Dy specifically, they found that the number of odd-parity states can be nearly double that of even-parity states in certain energy ranges. This isn't a universal rule for all nuclei, nor is it a constant ratio; rather, the balance fluctuates with excitation energy and is a direct result of the specific shape and arrangement of the energy levels inside this particular nucleus.

How Sure Are They?

It is important to note that these findings come from computer simulations using the Projected Shell Model, not from a new physical experiment. The authors have validated their model by showing it can perfectly reproduce the known, low-energy energy levels of 164Dy that have been measured in real experiments. Because the model gets the known facts right, they are confident that its predictions for the higher-energy, unmeasured regions are reliable.

However, they acknowledge that their current calculations are limited to energies up to about 5 MeV. While they suggest that their model reveals a "missing step" in the level density curve that experimentalists might have overlooked, they present this as a strong theoretical prediction based on their detailed quantum mechanical calculations. They invite future experiments to check if this "third step" is indeed there, which would change how we understand the transition from orderly nuclear motion to chaotic nuclear behavior.

In the end, this paper doesn't just give us a new number; it gives us a new story. It tells us that the nucleus doesn't just heat up and spin faster in a simple way. Instead, it goes through a complex evolution where pairs break one by one, creating distinct "steps" in the energy landscape, before finally dissolving into a chaotic, high-energy party. This detailed understanding helps scientists refine the tools they use to map the history of the universe and the behavior of stars.

Technical Summary: Projected Shell Model Description of Nuclear Level Density

Problem Statement
Nuclear level density (NLD) is a fundamental property of atomic nuclei, serving as a critical input for nuclear reaction theories, astrophysical nucleosynthesis networks (particularly the r-process), and reactor physics. While phenomenological models (e.g., constant temperature, back-shifted Fermi-gas) and combinatorial approaches have been widely used, they often rely on empirical parameter adjustments. Microscopic approaches, such as exact shell-model diagonalization, are computationally prohibitive for mid-mass and heavy deformed nuclei due to the immense dimensionality of the Hamiltonian matrix. Furthermore, experimental determination of NLD via the Oslo method, while advanced, relies on discrete level data that is often incomplete at higher excitation energies and lacks precise spin-parity resolution. There is a specific theoretical need to understand the NLD in the low-to-intermediate excitation region (before the onset of purely statistical exponential growth), where the density is highly sensitive to detailed nuclear structure, specifically the interplay between collective motion and the breaking of nucleon pairs.

Methodology
The authors propose a novel application of the Projected Shell Model (PSM) to calculate NLD systematically for deformed nuclei. The methodology involves:

Basis Construction: Starting from a deformed Nilsson potential with BCS pairing correlations, the model constructs a many-body configuration space based on multi-quasiparticle (qp) states. This space includes configurations up to 6-qp states (and theoretically up to 10-qp states using advanced Pfaffian algorithms) to account for broken nucleon pairs.
Angular Momentum Projection: The deformed basis states are projected onto good angular momentum ( $I$ ) and parity ( $\Pi$ ) quantum numbers to form shell-model configurations in the laboratory frame.
Hamiltonian Diagonalization: A two-body effective Hamiltonian, comprising a spherical single-particle term, quadrupole-quadrupole interaction, and monopole/quadrupole pairing interactions, is diagonalized within the projected space.
Case Study: The method is applied to the well-deformed rare-earth nucleus $^{164}$ Dy. The calculation is validated against known discrete energy levels and compared with NLD curves derived from the Oslo method and neutron resonance spacing data.

Key Contributions and Results
The study divides the NLD curve into three distinct regimes based on the underlying nuclear structure mechanisms:

Collective Regime (0 – 1.5 MeV):
- The model successfully reproduces the low-lying discrete levels, including the ground-state rotational band, $\gamma$ -vibrational bands, and octupole vibrational bands.
- The NLD in this region is dominated by collective excitations with all nucleons paired.
Pair-Breaking Regime (1.5 – 4.0 MeV):
- The authors identify a third step structure in the NLD curve beyond 2.0 MeV, which is absent in the standard Oslo curve (which assumes a linear exponential rise from ~2 MeV).
- This structure is attributed to the simultaneous breaking of a neutron pair and a proton pair (forming 4-qp states). The energy threshold for this corresponds to $2(\Delta_n + \Delta_p)$ , approximately 3 MeV in $^{164}$ Dy.
- The study argues that the exponential growth of NLD is delayed compared to previous assumptions because the system must overcome distinct neutron and proton pairing gaps before entering a fully chaotic multi-qp state.
Multi-Quasiparticle (Chaotic) Regime (> 4.0 MeV):
- At higher excitations, the system evolves into a state of chaotic motion where individual configurations mix heavily.
- The NLD curve begins to follow a power law consistent with the constant temperature model, aligning with the Oslo curve at ~4 MeV.

Specific Structural Findings:

Parity Dependence: Contrary to the common statistical assumption that even- and odd-parity states are equally distributed, the PSM calculation for $^{164}$ Dy reveals a significant imbalance. In the pair-breaking and multi-qp regimes, the odd-parity NLD can be nearly double the even-parity NLD. This is driven by the specific availability of deformed single-particle orbitals near the Fermi surfaces that allow for the formation of odd-parity 2-qp and 4-qp states at low excitation energies.
Spin Distribution:
- In the collective and pair-breaking regimes, the spin distribution is irregular and structure-dependent, showing pronounced odd-even staggering (particularly for even-parity states).
- In the multi-qp regime (> 4 MeV), the spin distribution converges to a Gaussian-like shape, consistent with the statistical model of random angular momentum coupling (Ericson formula).
- The authors successfully fit the high-energy spin distribution to the Ericson formula to extract the dispersion parameter $\sigma$ , demonstrating that $\sigma$ is a structure-dependent quantity that can be determined microscopically.

Significance and Claims
The paper claims to provide the first systematic microscopic description of NLD that connects statistical quantities directly to detailed nuclear structure questions. Its primary significance lies in:

Refining the Low-Energy Anchor: By identifying the "missing" step structure in the Oslo curve (the 4-qp contribution), the work suggests that the transition to exponential behavior occurs at higher energies than previously thought, offering a more accurate anchor point for experimental NLD determinations.
Structure Sensitivity: It demonstrates that NLD is not merely a smooth statistical function but is highly sensitive to the specific deformed single-particle structure, particularly regarding parity distribution and the sequential nature of pair breaking.
Methodological Bridge: The PSM approach offers a practical alternative to exact diagonalization for heavy nuclei, capable of providing spin- and parity-resolved NLDs in the critical transition region where phenomenological models fail and experimental data is sparse.

The authors conclude that while the current calculation is limited to lower excitations ( $E \lesssim 5$ MeV) due to computational constraints, the method provides a robust framework for understanding the structural evolution of NLD, with future work planned to address odd-mass nuclei and extend calculations to the neutron separation energy using configuration counting without full diagonalization.

The Atomic Dance Floor: Counting the Steps

The Projected Shell Model: A New Way to Count the Dancers

What the Paper Rules Out and What It Suggests

How Sure Are They?

More like this