Bulk and spectroscopic nuclear properties within an ab initio renormalized random-phase approximation framework

This study employs an ab initio renormalized random-phase approximation framework with modern chiral three-body forces to successfully calculate bulk and spectroscopic properties of closed-shell nuclei, demonstrating improved agreement with experiments by eliminating quasiboson approximation instabilities while highlighting the necessity of extending beyond the particle-hole space to resolve remaining discrepancies.

The Gist

Imagine the atomic nucleus not as a solid marble, but as a bustling, chaotic dance floor filled with tiny dancers (protons and neutrons). Physicists have been trying to predict exactly how these dancers move, how they hold hands, and what happens when the music changes (like when energy is added).

For decades, scientists have used a mathematical tool called RPA (Random Phase Approximation) to model this dance. Think of RPA as a "group choreography" calculator. It assumes that if one dancer moves, everyone else moves in a synchronized, predictable wave. It's a great shortcut, but it has a major flaw: it assumes the dancers are standing perfectly still and perfectly organized before the music starts. In reality, the dance floor is already jittery and messy even before the music begins.

This paper introduces a new, upgraded version of this calculator called RRPA (Renormalized RPA). Here is the story of what they did, explained simply:

1. The Problem: The "Perfectly Still" Lie

The old method (RPA) makes a simplifying assumption called the Quasi-Boson Approximation.

• The Analogy: Imagine trying to predict how a crowd reacts to a fire drill. The old method assumes everyone is standing in perfect, rigid rows, waiting for a command. It ignores the fact that people are already chatting, shifting their weight, and fidgeting.
• The Result: Because the old method ignores this "fidgeting" (ground-state correlations), it often predicts that the nucleus is too heavy, too big, or that certain dance moves (energy levels) happen at the wrong time. Sometimes, it even predicts impossible things, like the nucleus collapsing into nothingness.

2. The Solution: Acknowledging the Fidgeting

The authors, led by R. Folprecht and colleagues, decided to fix this by using a modern, realistic set of rules for how the dancers interact. They used a "chiral potential," which is like a very detailed instruction manual on how protons and neutrons actually push and pull on each other, including the fact that they sometimes need a third partner to interact properly (three-body forces).

They then applied their new RRPA method.

• The Analogy: Instead of assuming the dancers are in rigid rows, RRPA looks at the dance floor and says, "Okay, I see that the dancers are already jiggling and shifting. Let's calculate the dance moves based on that jiggling."
• The Magic: By accounting for this pre-existing movement, the math stops breaking. The "instabilities" (where the old math predicted the nucleus would collapse) disappear.

3. The Results: A Much Better Dance Floor

The team tested this new method on a huge variety of atomic nuclei, from the tiny Helium-4 (4 dancers) to the massive Lead-208 (208 dancers).

• Binding Energy (How tightly they hold hands): The old method said the nucleus was too loose. The new method predicted how tightly they hold hands with amazing accuracy, matching real-world experiments almost perfectly.
• Size (How big the dance floor is): The old method thought the nucleus was too big. The new method corrected the size to match reality.
• The "Spooky" Levels: In the old method, some energy levels were predicted to be imaginary or negative (mathematically impossible). The new method fixed these, giving realistic numbers that match what scientists see in labs.

4. Why This Matters

Usually, to get this level of accuracy, scientists have to use super-complex, super-slow computer programs (like Coupled-Cluster or IMSRG) that take weeks to run on the world's fastest supercomputers.

• The Analogy: If the old methods were like trying to solve a Rubik's cube by hand, and the super-complex methods were like using a robot arm that takes an hour to solve it, this new RRPA method is like a smart, high-speed app that solves it in seconds with the same accuracy.

The Bottom Line

This paper shows that you don't need to be a super-complex wizard to get perfect results. By simply acknowledging that the nucleus is "jittery" even before it starts moving, and by using a better set of interaction rules, the authors created a tool that is:

1. Accurate: It matches real experiments better than the old "simplified" models.
2. Efficient: It runs much faster than the most advanced methods currently available.
3. Universal: It works for light nuclei and heavy nuclei alike.

In short, they fixed the "group choreography" calculator so it accounts for the fact that the dancers are already moving, giving us a clearer, faster, and more reliable map of the atomic world. This is a big step forward for understanding everything from how stars burn to how we might use nuclear energy in the future.

Technical Summary

1. Problem Statement

The Random-Phase Approximation (RPA) is a standard tool for studying nuclear spectroscopy and response functions (e.g., giant resonances). However, standard RPA relies on the Quasi-Boson Approximation (QBA), which assumes the ground state is an uncorrelated Hartree-Fock (HF) vacuum. This simplification leads to several critical deficiencies when using modern, realistic chiral potentials (which include strong three-body forces):

• Instabilities: The QBA induces unphysical instabilities at low excitation energies, sometimes causing energy levels to collapse or become imaginary.
• Inaccurate Bulk Properties: Standard RPA tends to overestimate binding energies and charge radii compared to experimental data.
• Ground State Correlations: HF alone accounts for only a fraction of the binding energy in closed-shell nuclei. The lack of ground-state correlations in standard RPA prevents a consistent description of both bulk properties and excited states within a unified ab initio framework.

2. Methodology

The authors propose a Renormalized Random-Phase Approximation (RRPA) framework that restores ground-state correlations without abandoning the computational efficiency of the RPA eigenvalue structure.

• Interaction: The study utilizes the modern chiral potential $\Delta$N2LO$_{GO}$(394), which includes two- and three-body forces.
• Basis: A self-consistent Hartree-Fock (HF) basis is generated within an harmonic oscillator (HO) space (up to $N_{max}=14$). Three-body forces are treated at the normal-ordered two-body (NO2B) level.
• Formalism:
  • The excited states are defined as $|\nu\rangle = Q^\dagger_\nu |0\rangle$, where $|0\rangle$ is a correlated ground state rather than the HF vacuum.
  • The creation operator $B^\dagger_{ph}$ is renormalized by a factor $D_{ph} = n_h - n_p$, where $n_r$ are the ground-state occupation numbers derived from the one-body density matrix (OBDM).
  • Unlike standard RPA where $n_p=0$ and $n_h=1$ (yielding $D_{ph}=1$), RRPA calculates $n_r$ self-consistently, allowing them to deviate from 0 and 1.
• Iterative Solution: The method employs a double-iteration scheme:
  1. Solve the standard RPA eigenvalue problem to get initial amplitudes ($X, Y$).
  2. Use these amplitudes to calculate the OBDM and update occupation numbers ($n_r$) and the single-particle basis.
  3. Re-solve the RPA equations with the new basis and occupation numbers.
  4. Repeat until the OBDM converges to a fixed point.
• Scope: The framework is applied systematically to closed-(sub)shell nuclei ranging from $^4$He to $^{208}$Pb.

3. Key Contributions

• First Ab Initio RRPA with Chiral Forces: This is the first study to implement a renormalized RPA directly using state-of-the-art two- and three-body chiral interactions without intermediate smoothing or phenomenological adjustments.
• Removal of QBA Instabilities: By explicitly accounting for ground-state correlations and renormalizing the operators, the method eliminates the low-energy instabilities (level collapse) observed in standard RPA.
• Unified Framework: It provides a single, consistent approach to calculate both bulk properties (binding energies, radii) and spectroscopic properties (excitation spectra, transition strengths) referred to a correlated ground state.
• Computational Efficiency: The approach achieves accuracy comparable to high-cost methods like Coupled-Cluster (CC) and In-Medium Similarity Renormalization Group (IMSRG) but at a significantly lower computational cost.

4. Results

• Binding Energies and Radii:
  • Standard HF severely underestimates binding energies.
  • Standard RPA overestimates binding energies (over-binding) and charge radii.
  • RRPA restores consistency with experimental values, reproducing binding energies with accuracy comparable to CC and IMSRG(2) methods.
• Excitation Spectra:
  • In standard RPA, low-lying levels in nuclei like $^{16}$O and $^{90}$Zr are pushed down unrealistically, and in $^{40}$Ca, the lowest level collapses.
  • RRPA removes these instabilities, yielding stable spectra that agree well with experiments and CC calculations.
• Transition Strengths and Sum Rules:
  • E1 (Dipole) Transitions: RRPA causes a downshift and damping of the main peak. While the energy-weighted sum rule (EWSR) is preserved in RPA, it is underestimated in RRPA due to the depletion of single-particle states near the Fermi surface.
  • E3 (Octupole) Transitions: Standard RPA systematically overestimates the reduced strength $B(E3)$ for the $3^-$ state. RRPA significantly quenches this strength (due to the damping of backward $Y$ amplitudes), bringing it closer to experimental values, though it still slightly underestimates them.
• Occupation Numbers: The analysis of the one-body density matrix shows that deviations of occupation numbers from 0 or 1 are pronounced near the Fermi surface, validating the necessity of the renormalization.

5. Significance and Conclusion

• Validation of RRPA: The study demonstrates that RRPA is a robust, efficient tool for ab initio nuclear structure calculations, particularly for medium and heavy closed-shell nuclei where other methods (like No-Core Shell Model) become computationally prohibitive.
• Limitations: The authors acknowledge that RRPA is confined to the 1p-1h (particle-hole) configuration space. Consequently, it cannot fully satisfy energy-weighted sum rules or describe the full complexity of open-shell nuclei without further extensions (e.g., including $n$-particle-$n$-hole configurations or moving to a quasiparticle basis).
• Future Outlook: The work suggests that while RRPA is not a complete alternative to methods like Coupled-Cluster, it offers a "manageable" bridge between mean-field approaches and high-level ab initio methods. It opens a route for the broader nuclear physics community to utilize modern chiral potentials directly in mean-field-based frameworks.

In summary, the paper establishes RRPA as a critical advancement in nuclear theory, successfully resolving the instabilities of standard RPA and providing a unified, computationally efficient description of nuclear bulk and spectroscopic properties using modern realistic interactions.