Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

This study employs ab initio many-body methods (IMSRG and CC) alongside RPA to calculate isoscalar monopole sum rules in finite nuclei, demonstrating strong agreement between approaches and deriving nuclear matter incompressibility values consistent with phenomenological ranges.

The Gist

Imagine the atomic nucleus not as a static ball of particles, but as a tiny, vibrating drum. When you hit this drum, it doesn't just make a single "thud"; it rings with a specific pitch and a complex pattern of vibrations. In physics, this "ringing" is called a Giant Monopole Resonance. It's the nucleus breathing in and out, expanding and contracting as a whole.

This paper is like a high-tech acoustic engineering study. The researchers wanted to understand exactly how these nuclear "drums" vibrate and, more importantly, what that tells us about the material they are made of. Specifically, they wanted to measure the stiffness (or "squishiness") of nuclear matter. This stiffness is crucial for understanding how neutron stars behave and how supernovas explode.

Here is a breakdown of their journey, using simple analogies:

1. The Problem: How to Listen to the Drum

To figure out how stiff the nucleus is, scientists need to calculate the "moments" of its vibration. Think of a moment like a weighted average of the sound.

• The Hard Way: You could try to listen to every single possible note the drum can make, from the lowest hum to the highest squeak, and add them all up. This is incredibly difficult because there are infinite notes (states) to check.
• The Smart Way: Instead of listening to every note, you can use a mathematical trick to calculate the "average sound" directly from the drum's resting state.

2. The Tools: Three Different Microphones

The authors used three different "microphones" (mathematical methods) to listen to these nuclear drums:

• RPA (Random Phase Approximation): This is like using a basic, cheap microphone. It's fast and gives a decent idea of the sound, but it misses the subtle, complex harmonics. It works well if the drum is made of "soft" rubber, but it gets confused if the drum is made of "hard" steel.
• CC (Coupled-Cluster Theory): This is a high-end, studio-grade microphone. It captures every nuance and complex interaction between the particles. It's very accurate but computationally expensive (it takes a supercomputer a long time to process).
• IMSRG (In-Medium Similarity Renormalization Group): This is another high-end microphone, but it works differently. Instead of listening to the notes one by one, it transforms the drum itself to make the math easier, then calculates the sound from the ground up.

The Big Discovery: The authors found that the two expensive, high-end microphones (CC and IMSRG) agreed with each other perfectly. They also found that the cheap microphone (RPA) was actually quite good if the nuclear forces were "soft" (like the $\Delta$NNLOGO interaction), but it failed to capture the full picture for "harder" forces (like NNLOsat).

3. The Experiment: Testing Different Nuclei

They tested this on four specific "drums" (nuclei): Oxygen-16, Calcium-40, Nickel-56, and Tin-100. These are like different sizes of drums, from small hand-drums to large bass drums.

• They used two different types of "rubber" for the drum skin (two different nuclear interaction models).
• They found that for the "softer" rubber, all three microphones gave similar results.
• For the "harder" rubber, the cheap microphone (RPA) started to sound off-key, while the two expensive ones stayed in perfect harmony.

4. The Goal: Measuring the Stiffness of the Universe

Why do we care about the pitch of these nuclear drums?

• The Finite Nucleus: They calculated how stiff these specific, small nuclei are.
• The Infinite Matter: They then used a clever mathematical bridge (called a "leptodermous expansion," which is just a fancy way of saying "scaling up from a drop of water to an ocean") to predict how stiff infinite nuclear matter would be. This is the stuff that makes up the cores of neutron stars.

5. The Result: A Surprising Finding

When they scaled up their results to infinite nuclear matter, they found something interesting:

• Their calculated stiffness was lower than what other scientists had calculated using direct methods on infinite matter.
• However, their numbers fell right within the range of what experimentalists (people who actually measure things in labs) have observed.

The Analogy: Imagine trying to guess how stiff a giant block of Jell-O is by poking a tiny cube of it.

• Other scientists poked the giant block directly and said, "It's very stiff!"
• These authors poked the tiny cube, used a sophisticated scaling law, and said, "Based on this tiny cube, the giant block is a bit softer."
• Both answers are close to the "real" range, but the authors' method highlights that the way we scale up from small to big matters a lot.

Summary

This paper is a triumph of consistency. It shows that two very different, high-tech mathematical approaches (CC and IMSRG) agree perfectly when studying how atomic nuclei vibrate. It also suggests that while our current models are good, there are still subtle differences between how we calculate the stiffness of a tiny nucleus versus a giant chunk of nuclear matter.

Ultimately, this helps us refine our understanding of the Equation of State (the rulebook for how matter behaves under extreme pressure), which is essential for understanding the most violent and mysterious objects in the universe: neutron stars and supernovas.

Technical Summary

1. Problem Statement

The paper addresses the challenge of achieving a microscopic, ab initio understanding of the Isoscalar Giant Monopole Resonance (ISGMR) in atomic nuclei. The ISGMR is a collective vibration mode that serves as a critical probe for the nuclear equation of state (EOS), specifically the incompressibility of nuclear matter ($K_\infty$).

Historically, these properties were studied using phenomenological Energy Density Functionals (EDFs). While recent advances allow for calculations using chiral effective field theory ($\chi$EFT) interactions, there remains a need to:

• Validate different ab initio many-body frameworks against each other for collective excitations.
• Understand the impact of ground-state correlations on response functions.
• Bridge the gap between finite-nucleus properties and infinite nuclear matter incompressibility using consistent interactions.

2. Methodology

The authors employ two distinct ab initio many-body approaches to calculate the moments of the monopole strength function ($m_k$) for $N=Z$ closed-shell nuclei ($^{16}\text{O}$, $^{40}\text{Ca}$, $^{56}\text{Ni}$, and $^{100}\text{Sn}$). They utilize two different chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions: NNLOsat and $\Delta$NNLO$_{\text{GO}}$ (394).

A. Many-Body Frameworks

1. In-Medium Similarity Renormalization Group (IMSRG):
   • Approach: Uses a ground-state-based technique.
   • Mechanism: The moments are calculated as ground-state expectation values of moment operators ($M_k$). These operators are defined via nested commutators of the Hamiltonian and the excitation operator.
   • Implementation: The Hamiltonian and operators are evolved via a unitary transformation (IMSRG(2) truncation). The moments are obtained by evaluating $\langle \Phi_0 | e^{\Omega} M_k e^{-\Omega} | \Phi_0 \rangle$.
2. Coupled-Cluster (CC) Theory:
   • Approach: Uses an excited-state formulation.
   • Mechanism: The response is constructed using the Lorentz Integral Transform (LIT) combined with the Equation-of-Motion (EOM) CC method (LIT-CC).
   • Implementation: Excited states are generated by acting with excitation operators on the CC ground state. The strength function is discretized using bound pseudo-states, and moments are extracted by integrating the transform.
3. Random Phase Approximation (RPA):
   • Included for comparison to assess the impact of ground-state correlations. Calculations are performed using the same chiral interactions.

B. Extraction of Incompressibility

• Average Energy ($E_{\text{GMR}}$): Calculated from the ratio of moments: $E_{\text{GMR}} = \sqrt{m_1 / m_{-1}}$.
• Finite-Nucleus Incompressibility ($K_A$): Derived using the relation $K_A = \frac{m}{\hbar^2 \langle r^2 \rangle} E_{\text{GMR}}^2$.
• Infinite Matter Incompressibility ($K_\infty$): Extracted by fitting the calculated $K_A$ values to a leptodermous expansion (liquid-drop-like formula):
  $$K_A = K_{\text{vol}} + K_{\text{surf}} A^{-1/3}$$
  Here, $K_{\text{vol}}$ corresponds to the extrapolated $K_\infty$.

3. Key Contributions

• Methodological Consistency: Demonstrates excellent agreement between two fundamentally different ab initio approaches (IMSRG and CC) for calculating monopole sum rules across a range of nuclei.
• Interaction Dependence: Systematically compares the performance of a "harder" interaction (NNLOsat) versus a "softer" interaction ($\Delta$NNLO$_{\text{GO}}$ (394)).
• Ground-State Correlations: Quantifies the role of ground-state correlations by comparing correlated methods (IMSRG/CC) with RPA.
• Bridging Finite and Infinite Matter: Provides a consistent ab initio pathway to extract $K_\infty$ from finite nuclei, highlighting discrepancies with direct infinite matter calculations.

4. Key Results

A. Convergence and Method Comparison

• Agreement: IMSRG and CC results show good quantitative agreement for all studied nuclei and moments ($m_{-1}, m_0, m_1$).
• Interaction Softness:
  • $\Delta$NNLO$_{\text{GO}}$ (394): This softer interaction leads to rapid convergence with respect to the harmonic oscillator basis size ($e_{\text{max}}$) and frequency ($\hbar\omega$). For this interaction, RPA results are very close to the correlated methods (deviations $\sim 3\%$), suggesting that correlations cancel out effectively for this observable.
  • NNLOsat: This harder interaction shows slower convergence and larger dependence on model-space parameters. Significant differences (up to 20% for $m_0$ in $^{100}\text{Sn}$) appear between CC and IMSRG, attributed to the difficulty of converging even moments with hard interactions. RPA systematically underpredicts the results compared to correlated methods.
• Average Energies ($E_{\text{GMR}}$): The ratio $m_1/m_{-1}$ converges faster than individual moments. Theoretical predictions are generally higher than experimental data, though NNLOsat results are closer to experiment.

B. Incompressibility ($K_A$ and $K_\infty$)

• Finite Nuclei: $K_A$ values derived from IMSRG and CC are consistent. NNLOsat yields lower $K_A$ values than $\Delta$NNLO$_{\text{GO}}$.
• Extrapolated $K_\infty$:
  • IMSRG/CC Extrapolation:
    • NNLOsat: $K_\infty \approx 194 - 208$ MeV.
    • $\Delta$NNLO$_{\text{GO}}$: $K_\infty \approx 218 - 236$ MeV.
  • Comparison with Direct Calculations: These extrapolated values are significantly lower (by $\sim 50-60$ MeV for NNLOsat) than values obtained from direct infinite nuclear matter calculations using the same interactions.
  • Phenomenology: Despite the discrepancy with direct matter calculations, the extrapolated values fall within the phenomenological range ($250 \pm 30$ MeV is often cited, but the paper notes consistency with broader phenomenological ranges depending on the dataset).

5. Significance and Outlook

• Validation of Frameworks: The study confirms that IMSRG and CC are robust and consistent tools for describing collective nuclear responses, validating the newly developed moment-operator approach in IMSRG.
• Role of Correlations: The work highlights that for "soft" interactions, simpler methods like RPA may suffice for monopole sum rules, but "hard" interactions require full correlation treatments (IMSRG/CC) for accuracy.
• The "Finite vs. Infinite" Puzzle: The discrepancy between $K_\infty$ extracted from finite nuclei (via leptodermous expansion) and direct infinite matter calculations suggests that current ab initio methods may still miss subtle many-body effects or that the leptodermous expansion requires more data points (e.g., open-shell nuclei) and inclusion of Coulomb effects for higher precision.
• Future Directions: The authors emphasize the need to extend these calculations to open-shell nuclei and asymmetric matter ($N \neq Z$) to improve the stability of the extrapolation and better constrain the nuclear equation of state for astrophysical applications (neutron stars, supernovae).