Shell model description of the $N=82$ isotonic chain with a new effective interaction

This paper presents a systematic shell model study of the $N=82$ isotonic chain ($Z=51$–77) using a newly developed effective interaction derived via principal component analysis, which successfully reproduces experimental nuclear properties and provides predictions for proton-rich nuclei beyond current experimental reach.

The Gist

Imagine the atomic nucleus not as a solid ball, but as a bustling, multi-story apartment building where tiny particles called protons and neutrons live. In this building, there are specific "floors" or energy levels where these particles prefer to hang out. Sometimes, a floor is completely full, creating a very stable and happy neighborhood. In nuclear physics, we call these full floors "magic numbers."

This paper is about a specific neighborhood where the neutron floor is completely full (the magic number 82). The scientists wanted to understand how the protons behave in the floors above this stable base, specifically in a range of elements from Tellurium to Iridium.

Here is a breakdown of what they did and found, using simple analogies:

1. The Problem: The "Map" Wasn't Perfect

Scientists have been trying to draw a perfect map of how these protons interact with each other. Previous maps (called "effective interactions") were okay, but they had some errors. They were like a GPS that sometimes told you to turn left when you should have turned right, or predicted a building would be 10 feet tall when it was actually 12.

Specifically, old maps struggled to predict:

• The exact energy levels of certain excited states (like how high a ball bounces).
• The "spin" or orientation of heavier, odd-numbered nuclei (like predicting which way a spinning top will fall).
• The behavior of very heavy, proton-rich nuclei that are hard to study in the lab.

2. The Solution: A New, Smarter Map

The authors created a brand new, high-quality map using a method called Principal Component Analysis (PCA).

Think of this like tuning a massive, complex musical instrument with 165 different strings (the interaction parameters). Instead of trying to tune every single string perfectly by guessing, they used a smart algorithm to find the 30 most important strings that actually change the sound of the music. They then "tuned" these 30 strings by listening to 204 real-world experimental notes (data from actual nuclei) and adjusting the map until the music matched perfectly.

The result? A map that is incredibly accurate. The difference between their predictions and real-world measurements is tiny—only about the width of a single atomic nucleus (102 keV).

3. What They Discovered

With this new, precise map, they were able to describe the "neighborhood" in great detail:

• The "Sub-Closure" at Z=64: They confirmed that at a specific proton number (64, which is Gadolinium), there is a special "sub-floor" that acts like a mini-wall. This makes the nucleus extra stable and harder to excite, much like a building with a reinforced concrete floor in the middle. Their map showed this perfectly.
• Predicting the Unseen: Because their map is so reliable, they used it to predict the properties of nuclei that are so heavy and unstable that scientists haven't been able to measure them yet. They made specific predictions for nuclei like Tantalum-155, Tungsten-156, Rhenium-157, Osmium-158, and Iridium-159. They predicted things like whether these nuclei would hold together or fall apart (emit protons).
• Solving Mysteries: They solved a long-standing puzzle about the "ground state" (the resting position) of certain heavy nuclei. Old maps got the direction of the spin wrong for some of these; the new map got it right every time.

4. The Takeaway

This paper is essentially about building a better, more reliable "rulebook" for how protons behave in a specific region of the atomic world. By using a smarter mathematical approach to fit the data, they created a tool that not only explains what we already know but also confidently predicts what we haven't seen yet.

They didn't just fix the numbers; they provided a clear picture of the underlying structure of these atoms, showing exactly which "floors" the protons are living on and how they interact with their neighbors. This new rulebook is now available for other scientists to use for future studies of these heavy elements.

Technical Summary

Technical Summary: Shell Model Description of the N = 82 Isotonic Chain with a New Effective Interaction

Problem Statement
The nuclear structure of the $N=82$ isotonic chain, characterized by valence protons in the $g_{7/2}d_{5/2}s_{1/2}d_{3/2}h_{11/2}$ shell above the $^{132}\text{Sn}$ core, has been the subject of extensive study. While previous effective interactions (e.g., JJ56PNA, KH5082, CW5082) have provided insights, they exhibit systematic discrepancies when compared to experimental data. Notably, these interactions tend to overestimate the excitation energies of $3^-$ states and fail to correctly predict the ground-state spins and parities for heavier odd-$A$ isotones ($Z > 64$). Furthermore, reliable theoretical predictions are required for proton-rich nuclei near and beyond the proton drip line (e.g., $^{155}\text{Ta}$, $^{156}\text{W}$), where experimental data is scarce or non-existent.

Methodology
The authors employ the full-configuration-interaction shell model using the BIGSTICK code within a model space comprising all proton orbitals between $Z=50$ and $Z=82$. The core innovation lies in the development of a new, high-quality effective interaction optimized via Principal Component Analysis (PCA).

1. Parameter Space: The Hamiltonian includes 5 single-particle energies (SPEs) and 160 proton-proton two-body matrix elements (TBMEs), totaling 165 input parameters.
2. Optimization Strategy: Instead of fitting all parameters directly, the authors utilize PCA to truncate the parameter space. They perform a least-squares fit to 204 experimental data points, which include:
   • 21 nuclear binding energies.
   • 101 excitation energies in even-even nuclei (focusing on yrast and specific non-yrast states like $2^+_2, 3^-$, etc.).
   • 82 excitation energies in odd-$A$ nuclei.
3. PCA Implementation: The error matrix of the fit is diagonalized. The parameter space is truncated to the 30 principal components (linear combinations of the original parameters) associated with the smallest eigenvalues (smallest uncertainties). This approach mitigates overfitting while retaining the most sensitive parameters.
4. Constraints: Standard mass scaling factors are applied to TBMEs. Coulomb effects are implicitly incorporated through the fitting procedure. Effective charges ($e_\pi = 1.5e$) and quenching factors ($g_s = 5.586 \times 0.7$) are used for electromagnetic properties.

Key Contributions and Results

• Interaction Optimization: The new interaction achieves a root-mean-square deviation (RMSD) of 102 keV across the 204 data points. This represents a significant improvement over previous interactions, which often deviated by MeV scales for binding energies.
• Binding and Separation Energies: The model reproduces binding energies for $^{134}\text{Te}$ to $^{151}\text{Tm}$ with an RMSD of 59 keV. It successfully predicts atomic mass excesses and separation energies for nuclei beyond current experimental reach, including $^{156}\text{W}$ through $^{159}\text{Ir}$. The calculations suggest $^{153}\text{Lu}$ is beyond the proton drip line, while $^{154}\text{Hf}$ remains bound. For $^{155}\text{Ta}$, the model predicts single-proton emission, and for $^{156}\text{W}$, it predicts a positive one-proton separation energy but negative two-proton separation energy, consistent with observed suppression of two-proton decay.
• Low-Lying Spectra:
  • Even-Even Nuclei: The model reproduces the systematic trend of $2^+_1$ excitation energies, correctly identifying the peak at $^{146}\text{Gd}$ ($Z=64$) and the local maximum at $^{140}\text{Ce}$, confirming the $Z=64$ subshell closure. It resolves previous discrepancies by accurately predicting the $3^-$ excitation energies, particularly the minimum observed in $^{146}\text{Gd}$.
  • Odd-A Nuclei: The interaction correctly reproduces ground-state spins and parities across the entire chain, including the challenging $1/2^+$ ground state of $^{147}\text{Tb}$ and $3/2^+$ states of $^{155}\text{Ta}$ and $^{157}\text{Re}$, which were previously difficult to describe.
• Electromagnetic Properties:
  • B(E2) Values: The model generally agrees well with experimental $B(E2)$ values. However, deviations are noted for specific transitions in $^{138}\text{Ba}$ and $^{139}\text{La}$, suggesting a possible reversal in the ordering of calculated $9/2^+$ states in $^{139}\text{La}$.
  • Magnetic Moments: Calculated magnetic dipole moments ($\mu$) agree well with experiment (discrepancies $< 0.5 \mu_N$).
  • B(E3) Strengths: The model underestimates $B(E3; 3^- \to 0^+)$ strengths by a factor of roughly two. The authors attribute this to missing contributions from neutron excitations across the $N=82$ shell gap, noting that increasing the effective charge to 2.0 brings the results into agreement.
• Structural Analysis:
  • Subshell Gaps: Mean-field single-particle energies (MSPEs) reveal a subshell gap of 1.788 MeV at $Z=64$ (between $1d_{5/2}$ and $0h_{11/2}$) and 1.062 MeV at $Z=58$ (between $0g_{7/2}$ and $1d_{5/2}$).
  • Seniority: The study analyzes the seniority structure of yrast states. $6^+_1$ states in lighter isotones are dominated by $(0g_{7/2})^2$ or $0g_{7/2}1d_{5/2}$ configurations, while those in heavier isotones ($Z > 64$) correspond to $(0h_{11/2})^2$.
  • Isomerism: The model predicts several very low-lying excited states in odd-$A$ nuclei (e.g., a $1/2^+$ state in $^{155}\text{Ta}$ at 4 keV), suggesting potential isomerism.

Significance and Claims
The paper claims that the newly optimized interaction provides a systematic and reliable description of the $N=82$ isotonic chain, resolving long-standing discrepancies in excitation energies and ground-state spins that plagued previous models. By utilizing PCA to optimize the interaction against a large dataset, the authors demonstrate that a one-major-shell description with an optimized interaction is sufficient to capture complex structural features, such as the $Z=64$ subshell closure and the evolution of seniority patterns, without needing to invoke large-scale excitations across major shells.

The authors emphasize that this interaction serves as a solid foundation for future studies of nuclei near $N=82$. They explicitly cite its application (and ongoing use) in:

1. Studies of nuclear shape evolution and quantum phase transitions in rare-earth nuclei.
2. Systematic investigations of even-even $N=80$ isotones.
3. The first spectroscopic measurement of $^{150}\text{Yb}$.

The work concludes by providing the optimized interaction as supplemental material to facilitate further research in this mass region.