Electromagnetic Properties of the N=50 Isotones with the p35-i3 Hamiltonian

This paper presents and compares magnetic moments, quadrupole moments, and transition strengths ($B(M1)$ and $B(E2)$) for N=50 isotones calculated using the newly developed p35-i3 Hamiltonian and related variants, assessing their performance against experimental data to evaluate theoretical uncertainties in ab-initio derived shell model interactions.

The Gist

Imagine the atomic nucleus as a tiny, bustling city. In this city, the "residents" are protons and neutrons. Just like people in a city, these residents have specific neighborhoods (called orbitals) they prefer to live in, and they follow strict rules about how they can move and interact.

This paper is a report card for a specific group of atomic cities: those that have exactly 50 neutrons (the "N=50" group). The scientists wanted to see if their computer models could accurately predict how these cities behave when they are "excited" (like when they dance or spin) or how they react to magnetic and electric fields.

Here is a breakdown of the paper using simple analogies:

1. The Goal: Building a Better Map

For decades, physicists have tried to build a "map" (a mathematical formula called a Hamiltonian) that predicts how these nuclear cities behave.

• The Old Maps: Previous maps were good, but they were like hand-drawn sketches. They worked okay, but they missed some details.
• The New GPS: The authors created new, high-tech maps using a method called VS-IMSRG. Think of this as using a super-advanced satellite system that starts with the fundamental laws of physics (how protons and neutrons talk to each other) and refines the map until it matches reality perfectly.

They created three slightly different versions of this new map (named p35-i2, p35-i3, and p30-i3) to see which one was the most accurate.

2. The Test Drive: Magnetic and Electric Properties

To test these maps, the scientists looked at two main things:

• Magnetic Moments (The Compass): How much does the nucleus act like a tiny magnet?
• Quadrupole Moments (The Shape): Is the nucleus perfectly round like a ball, or is it squashed like a rugby ball or stretched like a football?
• Transitions (The Dance Moves): How much energy does it take for the nucleus to jump from one state to another? This is measured as B(M1) and B(E2).

They compared their new maps against real-world data collected from experiments.

3. The Results: The New Maps Win!

The paper found that the new p35-i3 map is the champion.

• Accuracy: When they checked the "compass" (magnetic moments) and the "shape" (quadrupole moments), the new map matched the real data much better than the old maps.
• The "Sweet Spot": They found that using 35 adjustable knobs (parameters) in their formula gave the best results. Using fewer (30) was okay, but 35 was the "Goldilocks" number—just right.
• The "Outlier" Case: There was one tricky city, Gallium-81, where the map got a little confused about a specific dance move. The scientists realized this was because two different "dance moves" (energy states) were so close together that they were mixing up. It's like trying to tell if a dancer is spinning left or right when they are doing both at the same time!

4. Why Does This Matter?

Think of the nucleus as a giant orchestra.

• The Old Theory: Sometimes the orchestra sounded a bit out of tune because the sheet music (the Hamiltonian) wasn't quite right.
• The New Theory: The new map acts like a conductor who knows exactly how every instrument (proton and neutron) should play.
• The Discovery: The paper confirms that for these specific nuclei (between Nickel-78 and Tin-100), the "music" is actually very simple. The protons are mostly just filling up one specific seat in the orchestra (the 0g9/2 orbital). Because they are all sitting in the same seat, they behave very predictably, almost like a choir singing in perfect unison.

5. The Takeaway

The scientists successfully updated the "rulebook" for how these atomic nuclei work.

• For the future: They have provided a better tool for other scientists to predict the properties of nuclei that haven't been discovered yet (especially near the edges of the map, like near Nickel-78 and Tin-100).
• The "Effective Charge": To make the math work, they had to pretend the protons were slightly "heavier" or more "charged" than they actually are. It's like a video game where you have to tweak the character's stats to make the physics engine work smoothly. They found that a specific tweak (an effective charge of 1.8) made the simulation match the real world perfectly.

In short: The authors built a better, more accurate computer model for how atomic nuclei with 50 neutrons behave. They proved that their new model is superior to older ones, giving us a clearer picture of the tiny, magnetic, and shape-shifting world inside the atom.

Technical Summary

1. Problem Statement

The nuclear shell model for protons in the $N=50$ isotone region (between $^{78}\text{Ni}$ and $^{100}\text{Sn}$) has historically relied on phenomenological Hamiltonians fitted to experimental data. While these models describe energy spectra well, there is a need to assess theoretical uncertainties and improve predictive power for electromagnetic observables (magnetic moments, quadrupole moments, and transition probabilities) using modern ab initio starting points.

Specifically, the authors aim to:

• Evaluate the performance of new effective Hamiltonians derived from Valence-Space In-Medium Renormalization Group (VS-IMSRG) methods.
• Quantify the theoretical uncertainties associated with different Hamiltonian choices (varying the number of Singular Value Decomposition parameters and the ab initio starting point).
• Compare these new models against older phenomenological Hamiltonians ($jj44a$, $n50j$) and experimental data for electromagnetic observables.

2. Methodology

Model Space and Hamiltonians

The study focuses on the proton model space $\pi j^4 = \{0f_{5/2}, 1p_{3/2}, 1p_{1/2}, 0g_{9/2}\}$. The authors utilize three newly derived Hamiltonians based on Singular Value Decomposition (SVD) fits to experimental binding and excitation energies:

• p35-i2: 35 SVD parameters, starting from IMSRG(2) (standard approximation truncating operators at the two-body level).
• p35-i3: 35 SVD parameters, starting from IMSRG(3f2) (includes factorized three-body operators in nested commutators).
• p30-i3: 30 SVD parameters, starting from IMSRG(3f2) (used to assess sensitivity to the number of fit parameters).

The ab initio inputs were generated using the EM 1.8/2.0 NN+3N interaction in a harmonic oscillator basis ($\hbar\omega = 12$ MeV, truncated to 13 shells). The IMSRG(3f2) method was found in previous work [11] to provide a superior starting point with lower RMS deviations for energy data compared to IMSRG(2).

Calculation of Observables

Electromagnetic observables were calculated using the NuShellX code. The reduced matrix elements were expressed as a product of One-Body Transition Densities (OBTD) and Single-Particle Matrix Elements (SPME):

• Magnetic Moments ($\mu$) and B(M1): Calculated using orbital-dependent effective $g$-factors fitted to data to account for configuration mixing outside the model space and mesonic-exchange currents.
• Quadrupole Moments ($Q$) and B(E2): Calculated using realistic, nucleus-dependent radial wavefunctions derived from Energy-Density Functionals (EDF) using the Skx Skyrme functional. An effective proton charge of $e_p = 1.8$ was determined to reproduce experimental E2 data, replacing the standard harmonic oscillator assumption ($b^2=4.481$) and effective charge $e_p=2.0$.

Uncertainty Assessment

Theoretical uncertainties were assessed by comparing results across the different Hamiltonians (p35-i2 vs. p35-i3 vs. p30-i3) and against the older $jj44a$ Hamiltonian. The study separated nuclei into two groups based on mass number $A$:

• $A < 88$: Dominated by $\{0f_{5/2}, 1p_{3/2}, 1p_{1/2}\}$ orbitals.
• $A \ge 88$: Dominated by $\{1p_{1/2}, 0g_{9/2}\}$ orbitals.

3. Key Contributions

1. Validation of IMSRG(3f2) as a Starting Point: The study confirms that starting SVD fits with IMSRG(3f2) (Hamiltonian p35-i3) yields better agreement with experimental data and lower theoretical uncertainties compared to IMSRG(2) (p35-i2), particularly for nuclei where the wavefunctions are less dominated by a single orbital.
2. Refined Radial Wavefunctions: The adoption of Skyrme EDF-based radial wavefunctions with an effective charge $e_p = 1.8$ significantly improved the reproduction of quadrupole moments and B(E2) values compared to traditional harmonic oscillator assumptions.
3. Quantification of Theoretical Uncertainties: The authors provide a systematic breakdown of RMS deviations between different Hamiltonians, establishing a baseline for theoretical error bars in shell-model predictions for this mass region.
4. Identification of Configuration Mixing Effects: The paper highlights specific cases (e.g., $^{81}\text{Ga}$, $^{94}\text{Ru}$, $^{95}\text{Rh}$) where discrepancies between theory and experiment arise from subtle configuration mixing or level interference, offering insights into the limitations of the current model space.

4. Results

Magnetic Moments

• Agreement: The p35-i3 Hamiltonian shows excellent agreement with experimental magnetic moments, with an RMS deviation of 0.50 $\mu_N$.
• Uncertainty: This experimental deviation is roughly double the theoretical Hamiltonian uncertainty (0.24 $\mu_N$) derived from comparing p35-i3 and p30-i3.
• Outliers: $^{81}\text{Ga}$ shows a larger experimental moment than predicted, suggesting potential additional orbital dependence in the effective $g$-factors for the $0f_{5/2}$ orbital.

Quadrupole Moments and B(E2)

• Effective Charge: Using $e_p = 1.8$ with EDF radial wavefunctions resulted in an RMS deviation of 3.0 e-fm$^2$ between theory and experiment, consistent with the theoretical uncertainty of 5 e-fm$^2$.
• Comparison: The p35-i3 Hamiltonian outperforms the older $jj44a$ and $n50j$ Hamiltonians, particularly for the $0^+_1 \to 2^+_1$ transitions in even-even nuclei.
• Seniority Patterns: For nuclei with $A \ge 90$, the $8^+_1$ isomeric states in even-even nuclei exhibit uniform $g$-factors and B(E2) values, consistent with pure $\left[0g_{9/2}\right]^n$ configurations with good seniority ($\nu=2$).

Transition Strengths and Anomalies

• B(M1): Small B(M1) values are highly sensitive to Hamiltonian changes due to cancellations in the transition matrix elements.
• Specific Discrepancies:
  • $^{81}\text{Ga}$: The p30-i3 Hamiltonian correctly predicts the gamma-decay branching ratio for the $7/2^-$ state, whereas p35-i3 does not, suggesting sensitivity to the specific SVD parameters.
  • $^{94}\text{Ru}$ and $^{95}\text{Rh}$: Significant discrepancies were found in B(E2) values for specific transitions ($8^+ \to 6^+$ in Ru; $21/2^+ \to 17/2^+$ in Rh). The authors attribute this to mixing between nearby levels of the same spin-parity (yrare vs. yrast bands). A small mixing matrix element (~25 keV) could resolve these discrepancies, potentially requiring a readjustment of Two-Body Matrix Elements (TBME) or inclusion of three-body interactions.

5. Significance

This work represents a critical step in bridging ab initio nuclear theory with phenomenological shell-model calculations. By demonstrating that IMSRG(3f2) provides a robust starting point for effective Hamiltonians, the authors validate the use of modern many-body methods to constrain shell-model parameters.

The study establishes that:

1. Theoretical uncertainties in electromagnetic observables for N=50 isotones are now quantifiable and relatively small (within ~20-30% for moments, slightly higher for transitions).
2. The p35-i3 Hamiltonian is currently the state-of-the-art tool for predicting properties of nuclei in this region, including those near the limits of stability ($^{78}\text{Ni}$ and $^{100}\text{Sn}$).
3. Future improvements require addressing specific configuration mixing effects and potentially incorporating three-body forces directly into the effective TBME determination, rather than relying solely on SVD fits to energy data.

The paper provides a comprehensive benchmark for future experimental campaigns at facilities like FRIB (Facility for Rare Isotope Beams), offering precise predictions for unmeasured magnetic moments and transition rates.