Benchmarking projected generator coordinate method for nuclear Gamow-Teller transitions

This paper benchmarks a minimal extension of the quantum-number projected generator coordinate method (PGCM) for describing Gamow-Teller transitions and neutrinoless double-beta decay matrix elements in even-even nuclei, demonstrating its validity against exact shell-model solutions and configuration-interaction calculations in calcium and titanium isotopes.

The Gist

The Big Picture: Predicting the Unpredictable

Imagine you are trying to predict exactly how a complex machine will behave when you press a specific button. In the world of atoms, that "machine" is the atomic nucleus, and the "button" is a process called beta decay.

Scientists need to know exactly how nuclei behave during these decays to answer huge questions: How are elements made in stars? Do neutrinos have mass? To get these answers, they need to calculate something called a Nuclear Matrix Element (NME). Think of the NME as the "instruction manual" for how the nucleus rearranges itself. If the manual is wrong, our predictions about the universe are wrong.

The Problem: The "Exact" Way is Too Hard

To get the perfect manual, you would need to track every single particle inside the nucleus (protons and neutrons) and how they interact with every other particle.

• The Analogy: Imagine trying to predict the weather by tracking the path of every single air molecule in the atmosphere. It's theoretically possible, but the math is so massive that even the world's fastest supercomputers can't solve it for anything bigger than a tiny atom.
• The Reality: For heavier atoms (like Calcium or Titanium), the "Exact Solution" (called the Shell Model) is too computationally expensive to run easily. Scientists need a shortcut—a smart approximation.

The Solution: The "Group Photo" Method (PGCM)

The authors of this paper are testing a shortcut method called the Projected Generator Coordinate Method (PGCM).

• The Analogy: Imagine you want to describe a crowd of people.
  • The Exact Way: You take a high-resolution photo of every single person standing perfectly still.
  • The PGCM Way: You take a few "group photos" where the crowd is slightly shifted or deformed (maybe they are leaning left, leaning right, or stretching). Then, you use a computer to mathematically blend these photos together to create a "super-image" that represents the crowd's average behavior.
• The Twist: The authors extended this method to handle Odd-Odd Nuclei.
  • Most nuclei have an even number of protons and neutrons (like a perfect dance couple).
  • Some have an odd number of both (like a dance couple where one partner is missing a shoe, or they are dancing awkwardly).
  • The paper shows how to use the "group photo" method to describe these awkward, odd-numbered dancers, which is much harder to do.

The Experiment: Testing the Shortcut

The researchers tested their method on Calcium (Ca) and Titanium (Ti) isotopes. They used a known, perfect mathematical model (the "Exact Solution") as a ruler to measure how well their "PGCM shortcut" worked.

The Results:

1. Low Energy (The Calm Dancers): When the nucleus is in a calm, low-energy state, the PGCM method works beautifully. It matches the "Exact Solution" almost perfectly.
2. High Energy (The Wild Dancers): When the nucleus gets excited and starts moving wildly, the PGCM method starts to drift a bit. It's still a good guess, but not perfect.
3. The Double-Beta Decay Test: They calculated the "instruction manual" (NME) for a specific decay (Calcium-48 turning into Titanium-48).
   • The Result: Their method predicted a value that was about 57% too high.
   • Why? The method overestimated how easily the nucleus could jump to its very first excited state. It was like predicting a car would accelerate too fast because it overestimated the engine's power for that specific gear.

The Comparison: PGCM vs. The "Cut-and-Paste" Method

The paper also compared their method to another popular shortcut called Configuration Interaction (CI).

• CI Analogy: Imagine trying to fix a broken machine by only looking at the first two broken parts (1-particle, 1-hole) or the first four (2-particle, 2-hole). It's a "cut-and-paste" approach.
• The Verdict: The PGCM method (the "group photo" blend) performed just as well as, and sometimes better than, the CI method, even though CI is a very standard tool.

The Conclusion: A Promising, But Imperfect, Tool

What did they learn?
The PGCM method is a powerful new tool for understanding how nuclei behave, especially for atoms that aren't too heavy. It's like having a very good map for a city you've never visited.

What's next?
The map has some errors (the 57% overestimation). The authors suggest that to fix this, they need to:

1. Take more photos: Use more "deformed" configurations in their blend (expand the generator coordinates).
2. Use a better camera: Incorporate a technique called IMSRG (which acts like a lens that cleans up the noise in the data).

In Simple Terms:
The paper says, "We built a new, smart calculator to predict how atomic nuclei change. It works really well for simple cases and is better than some older calculators. It's not perfect yet—it gets a bit too excited in complex scenarios—but with a few upgrades, it could become the gold standard for predicting nuclear behavior."

Technical Summary

1. Problem Statement

Nuclear weak processes, such as single-$\beta$ and double-$\beta$ decays, are critical for understanding nuclear stability, nucleosynthesis, and physics beyond the Standard Model. Calculating the Nuclear Matrix Elements (NMEs) for these processes is a formidable challenge due to the complexity of the nuclear many-body problem.

• The Specific Challenge: While the Quantum-Number Projected Generator Coordinate Method (PGCM) has been successfully combined with the In-Medium Similarity Renormalization Group (IMSRG) to study neutrinoless double-$\beta$ ($0\nu\beta\beta$) decay, extending it to Gamow–Teller (GT) transitions and two-neutrino double-$\beta$ ($2\nu\beta\beta$) decay is significantly more difficult.
• Reason: Unlike $0\nu\beta\beta$ decay, which can often use the closure approximation, $2\nu\beta\beta$ decay and GT transitions require an accurate description of not just the ground states, but also a large number of intermediate excited states in odd–odd nuclei. These systems involve numerous nearly degenerate configurations that are difficult to capture with standard mean-field approaches.

2. Methodology

The authors developed a minimal extension of the PGCM framework to handle odd–odd nuclear systems and benchmarked it against exact shell-model solutions.

• Hamiltonian: The study utilizes the GXPF1A shell-model Hamiltonian defined within the $fp$ shell (orbitals $0f_{7/2}, 0f_{5/2}, 1p_{3/2}, 1p_{1/2}$) on top of an $^{40}\text{Ca}$ inert core. This allows for comparison with exact diagonalization results.
• PGCM Framework for Odd-Odd Nuclei:
  • Even-Even Nuclei: Wave functions are constructed as superpositions of symmetry-projected Hartree–Fock–Bogoliubov (HFB) states with varying intrinsic quadrupole deformations ($\beta_2$).
  • Odd-Odd Nuclei (Intermediate States): Wave functions are approximated as superpositions of quantum-number-projected two-quasiparticle configurations. These are built on an HFB reference state where the average neutron and proton numbers are constrained to be odd.
  • Projection: Angular momentum ($J$) and particle numbers ($N, Z$) are restored via projection operators. The basis states for odd–odd nuclei involve a quasiparticle pair creation operator ($A^\dagger_{pn} = \beta^\dagger_p \beta^\dagger_n$) acting on the HFB vacuum.
  • Solving: The Hill-Wheeler-Griffin (HWG) equation is solved by diagonalizing the norm and Hamiltonian kernels. A cutoff ($E_{cut} = 48$ MeV) is applied to quasiparticle configurations to ensure convergence.
• Comparative Methods:
  • Exact Shell Model (SM): Full configuration interaction (CI) within the $fp$ shell.
  • Truncated CI: CI calculations with particle–hole truncations up to 1p1h and 2p2h.
  • IMSRG Evolution: The authors also tested the CI method with Hamiltonians and GT operators evolved using the IMSRG to assess the impact of dynamic correlations.

3. Key Contributions

• Framework Extension: Successfully generalized the PGCM to describe GT transitions from even–even nuclei to the low-lying and giant resonance states of odd–odd nuclei.
• Benchmarking: Provided a rigorous benchmark of the PGCM against exact shell-model solutions for Ca and Ti isotopes ($^{42-48}\text{Ca}$ and $^{42-48}\text{Ti}$).
• $2\nu\beta\beta$ Decay Calculation: Calculated the NME for the $2\nu\beta\beta$ decay of $^{48}\text{Ca} \to {}^{48}\text{Ti}$ without invoking the closure approximation, relying instead on the explicit summation over intermediate $1^+$ states.
• Analysis of Discrepancies: Identified the specific physical origins of deviations between PGCM and exact results, particularly regarding shape mixing and the truncation of excitation configurations.

4. Results

• GT Transition Strengths:
  • The PGCM generally reproduces GT transition strengths to both low-lying and giant resonance states in Ca and Ti isotopes well.
  • Performance Trend: Agreement with exact results is excellent for nuclei near closed shells (e.g., $^{42}\text{Ca}$) but deteriorates slightly as the number of valence nucleons increases (e.g., $^{48}\text{Ca}$), reflecting the growing importance of complex many-body correlations.
  • Comparison with CI: The PGCM performance is comparable to, and in some cases superior to, CI calculations truncated at the 2p2h level. The PGCM captures collective correlations (via shape mixing) better than low-order CI, while CI captures specific 2p2h excitations better than the minimal PGCM.
  • Shape Mixing: Using a GCM wave function (mixing different deformations) for the initial state significantly improves the description of GT strengths compared to using a single spherical or deformed state.
• $2\nu\beta\beta$ Decay NME ($^{48}\text{Ca} \to {}^{48}\text{Ti}$):
  • The calculated NME using PGCM overestimates the exact shell-model value by approximately 57%.
  • Root Cause: The overestimation is primarily driven by an overestimation of the GT transition strength from $^{48}\text{Ti}$ to the first excited state ($1^+_1$) of $^{48}\text{Sc}$.
  • Convergence: The cumulative NME converges slowly, with most positive contributions coming from the lowest few intermediate states below 5 MeV.
• IMSRG Impact:
  • Applying IMSRG evolution to the CI method improved results for 1p1h truncation but remained slightly inferior to PGCM.
  • However, CI(2p2h) + IMSRG showed slightly better agreement with exact results than PGCM, suggesting that combining IMSRG with PGCM is a promising future direction.

5. Significance and Conclusions

• Validation: The study validates the PGCM as a reliable framework for describing $\beta$-decay observables in nuclei near closed shells, offering a computationally efficient alternative to full shell-model diagonalization for medium-mass nuclei.
• Limitations: As the valence space expands, the minimal PGCM (using only the lowest HFB state and limited quasiparticle pairs) struggles to capture all complex correlations, leading to deviations in transition strengths.
• Future Outlook: The authors conclude that the predictive power of the PGCM can be significantly enhanced by:
  1. Extending the Generator Coordinate Space: Including more collective coordinates (e.g., higher-order deformations, pairing gaps).
  2. Incorporating IMSRG: Combining the static correlations of PGCM with the dynamic correlations captured by the in-medium similarity renormalization group.
• Impact: This work lays the groundwork for applying PGCM to more complex nuclear systems and for calculating NMEs for $0\nu\beta\beta$ decay with higher precision, which is crucial for interpreting experimental searches for neutrinoless double-beta decay.