Nucleon-pair truncation of the shell model for medium-heavy nuclei

This paper proposes and validates an efficient nucleon-pair truncation scheme for the configuration-interaction shell model, which combines variationally optimized pair condensates with angular momentum projection to accurately describe low-lying states and shape coexistence in medium-heavy nuclei where full calculations are computationally prohibitive.

The Gist

Imagine the atomic nucleus as a bustling, crowded dance floor filled with protons and neutrons. Physicists want to understand how these dancers move, pair up, and spin together to create the shapes and energy levels of different atoms. The most accurate way to do this is the "Shell Model," which tries to track every single dancer's move. However, for medium-to-heavy nuclei, the number of possible dance combinations is so huge (like trying to count every grain of sand on a beach) that even the world's fastest supercomputers get stuck. They simply can't calculate all the possibilities in a reasonable amount of time.

This paper proposes a clever shortcut, a new "truncation scheme" (a way to cut down the work without losing the important details), called PNBCS.

Here is how the authors' method works, broken down into simple concepts:

1. The Problem: Too Many Dancers

The "Full Shell Model" is like trying to write a script for a play where every single actor has to improvise every line and movement simultaneously. It's perfect, but the script is too long to ever finish writing. For heavy nuclei, the "script" (the mathematical calculation) becomes too big for computers to handle.

2. The Solution: The "Pairing" Shortcut

The authors realized that in these nuclear dances, particles often move in pairs. Instead of tracking every individual dancer, they decided to focus on the pairs.

• The Setup: First, they use a standard method (Hartree-Fock) to find the best "dance floor" layout. This gives them a starting shape for the nucleus.
• The Pairing: They then use a method called NBCS (Number Conserved Bardeen-Cooper-Schrieffer). Think of this as organizing the dancers into specific couples that move in sync. Unlike older methods that might lose track of the total number of dancers, this method is strict: it ensures the exact number of protons and neutrons is preserved, just like a bouncer checking IDs at the door.
• The Spin: The initial pairing creates a shape that might be tilted or spinning in a messy way. To fix this, they use a mathematical "filter" called Linear Algebra Projection (LAP). Imagine taking a blurry, spinning photo of the dance floor and using a filter to snap a crystal-clear picture of the dance from a specific angle (good angular momentum). This step is very fast, unlike older methods that required slow, heavy calculations.

3. The Results: A Clearer Picture

The authors tested this new "PNBCS" method on a variety of nuclei, from Titanium to Xenon and beyond.

• The Test: They compared their shortcut method against the "Full Shell Model" (the gold standard) where possible.
• The Outcome: For nuclei that are somewhat round (spherical) to those that are stretched out like rugby balls (deformed), their method produced results that matched the expensive, full calculations almost perfectly.
• The "Shape Coexistence" Discovery: Some nuclei are like chameleons; they can exist in two different shapes at the same time (like a ball that is both round and squashed). The paper found that to describe these tricky nuclei correctly, you need two things: the pairing of the dancers and the ability to mix different "dance routines" (configurations). Their method captures both of these effects well.

4. Predicting the Unseen

Because their method is so fast and accurate, they used it to predict the behavior of nuclei that are currently too difficult for supercomputers to study, such as certain isotopes of Barium and Cerium. They provided a "map" of what these nuclei's energy levels likely look like, filling in gaps that were previously unreachable.

The Bottom Line

The paper introduces a fast, efficient way to study the complex dance of atomic nuclei. By focusing on how particles pair up and using a quick mathematical filter to clean up the results, they can study heavy, complex atoms that were previously too computationally expensive to analyze. It's like finding a way to predict the outcome of a massive, chaotic dance party by focusing on the key couples and their rhythm, rather than trying to track every single footstep.

Technical Summary

Technical Summary: Nucleon-pair truncation of the shell model for medium-heavy nuclei

Problem Statement
The configuration-interaction shell model (SM) is a versatile framework for nuclear structure, capable of generating low-lying states and complex multi-particle correlations. However, for medium- and heavy-mass rotational nuclei, the full configuration-interaction (FCI) dimensions often exceed current computational limits (typically $\sim 2-3 \times 10^{10}$). While existing truncation schemes like the generalized seniority scheme or the nucleon-pair approximation (NPA) reduce the basis size, they face challenges: seniority schemes work best for spherical nuclei, and NPA calculations involving high-spin pairs (G and I pairs) remain computationally burdensome for well-deformed nuclei. Conversely, mean-field approaches like Hartree-Fock-Bogoliubov (HFB) or Bardeen-Cooper-Schrieffer (BCS) handle deformation and pairing efficiently but break rotational symmetry and particle number conservation. Restoring these symmetries via projection is often computationally prohibitive, particularly when performing "variation after projection."

Methodology
The authors propose a computationally efficient truncation scheme termed Projected Number-conserved BCS (PNBCS). The methodology proceeds in three main stages:

1. Hartree-Fock (HF) Basis Generation: The authors generate single-particle states using an HF calculation within a shell model basis that enforces Kramers degeneracy. This ensures degenerate time-reversal partners without imposing additional constraints on shape or parity.
2. Particle-Number Conserved BCS (NBCS): Building on the HF states, the authors construct a variational pair condensate ansatz for open-shell nuclei (containing both valence protons and neutrons). Unlike standard BCS, which violates particle number, this approach utilizes a particle-number conserved formalism (NBCS) where the ground state is an $N$-pair condensate. The pair structure coefficients are optimized by minimizing the energy using a conjugate gradient method and iterative formulas derived from the variational principle. This step incorporates pairing correlations explicitly.
3. Angular Momentum Projection (LAP): Since the HF and NBCS states break rotational symmetry, the authors apply a linear algebraic projection (LAP) to restore good angular momentum. Unlike traditional numerical integration over Euler angles, LAP solves a set of linear equations to extract Hamiltonian and norm matrix elements. This approach is reported to be more than 10 times faster than quadrature methods. The final spectrum is obtained by diagonalizing the Hamiltonian in the space spanned by the projected states.

Key Results
The PNBCS scheme was applied to nuclei across the $pf$, $pf_{5/2}g_{9/2}$, and $sdg_{7/2}h_{11/2}$ shells, covering transitional and deformed regions:

• Shape Evolution: In a schematic model varying the ratio of pairing to quadrupole-quadrupole interactions, PNBCS successfully reproduced the evolution from spherical to rotational limits. It provided superior descriptions of excitation energies and $B(E2)$ transition rates compared to angular-momentum projected HF (PHF), particularly in transitional regions ($2.2 < R_{4/2} < 2.8$).
• Medium-Mass Nuclei ($pf$ and $pf_{5/2}g_{9/2}$ shells): For nuclei such as $^{44,46,48}$Ti, $^{48,50}$Cr, $^{52}$Fe, $^{60,62,64}$Zn, $^{66,68}$Ge, and $^{68}$Se, PNBCS results showed good agreement with full SM calculations and experimental data.
  • The inclusion of pair correlations was found to reduce the moment of inertia (improving level spacing) while deformation (encoded in HF) primarily dictated electric quadrupole transition strengths.
  • Shape Coexistence: In $^{66}$Ge and $^{68}$Ge, where shape coexistence (oblate/triaxial or triaxial minima) is observed, the authors demonstrated that both pair correlations and configuration mixing between different intrinsic HF states are crucial for reproducing side-band energies and $B(E2)$ values. In contrast, for $^{68}$Se, configuration mixing was found to be less critical.
• Heavy-Mass Predictions: The authors applied PNBCS to nuclei where full SM calculations are currently infeasible due to dimensionality ($2 \times 10^{10}$ to $5 \times 10^{13}$). They provided predictions for low-lying states of $^{108,110}$Xe, $^{112,114}$Ba, and $^{116,118,120}$Ce. The results for $^{108,110}$Xe agreed well with available data and SM benchmarks.

Significance and Claims
The paper claims that PNBCS offers a practical and efficient truncation scheme for the shell model, specifically targeting transitional and deformed nuclei in medium- and heavy-mass regions. The primary significance lies in its ability to:

1. Balance Accuracy and Efficiency: It captures essential pairing correlations and configuration mixing while remaining computationally fast, avoiding the heavy multidimensional integrations required by projected HFB or general pair condensate methods.
2. Extend Reach: It enables the study of nuclei (e.g., $^{112-120}$Ba/Ce) that are currently beyond the reach of large-scale configuration-interaction SM calculations.
3. Reproduce Collectivity: The study demonstrates that both pair correlations and the mixing of different intrinsic states are key to reproducing collective features and shape coexistence.

The authors note that while the current implementation uses HF single-particle states as a fixed basis for the NBCS variation, future work could explore variation after projection or the inclusion of broken-pair configurations to further improve results for phenomena like backbending.