Light antiproton-nucleus systems at low energies with the ab initio NCSM/RGM method

This paper extends the ab initio No-Core Shell Model combined with the Resonating Group Method (NCSM/RGM) to study low-energy antiproton-nucleus dynamics in light systems (${\bar p}+d$, ${\bar p}+{}^3 \mathrm{H}$, and ${\bar p}+{}^3\mathrm{He}$), validating the approach against exact solutions to isolate uncertainties in the $N\bar{N}$ interaction while addressing the significant numerical challenges posed by the hard short-range components of the interaction.

The Gist

Imagine the nucleus of an atom as a bustling, tightly packed city made of tiny citizens called protons and neutrons. Now, imagine a "ghost" visitor arriving: an antiproton. This ghost is the exact opposite of a proton. When it touches a citizen, they don't just say hello; they annihilate each other in a flash of energy, turning into pure light (mesons/pions).

Scientists at CERN want to use these "ghost visitors" to take a picture of the city's edges. They think the ghost is so sensitive that it will only interact with the outer walls of the city, revealing secrets about how the city is built that normal probes can't see.

This paper is about building a super-accurate map to predict exactly how these ghost visitors behave when they approach small atomic cities (like Deuterium, Tritium, and Helium-3).

Here is the breakdown of their work, using everyday analogies:

1. The Challenge: The "Hard" Interaction

The scientists used a powerful mathematical tool called NCSM/RGM. Think of this tool as a high-resolution camera that takes pictures of the atomic city. Usually, this camera works great for normal visitors (protons).

However, the antiproton is tricky. When it gets close to the city, the interaction is incredibly "hard" and violent (like a sledgehammer hitting a glass window).

• The Problem: When the scientists tried to use their standard camera settings (mathematical grids) to capture this violent crash, the images came out blurry and full of static noise. The math struggled to handle the sudden, intense force.
• The Analogy: Imagine trying to take a photo of a lightning strike with a slow shutter speed. You get a blur. Or, imagine trying to fit a square peg (the hard math) into a round hole (the smooth atomic wave). The math kept "ringing" with false echoes (artifacts) because the interaction was too sharp for the grid they were using.

2. The Solution: The "Noise-Canceling" Filter

To fix the blurry photos, the team invented a special filter (called a regulator).

• How it works: They realized the "noise" only happened far away from the center of the city, where the math was struggling to represent the sharp edge of the interaction. They applied a smooth "fade-out" filter to the math at the edges.
• The Result: This smoothed out the static noise without changing the actual picture of the crash in the center. Suddenly, their calculations became stable and accurate, even with the violent antiproton interaction.

3. The Discovery: The "Perimeter Patrol"

Once they had a clear map, they looked at where the antiproton actually "dies" (annihilates).

• The Finding: They confirmed that the antiproton doesn't crash into the center of the city. Instead, it gets caught in the outer suburbs (the "tail" of the nuclear density).
• Why it matters: This is like a burglar who only picks the lock on the front door and never goes inside. Because the antiproton only interacts with the edge, it is a perfect tool for measuring the thickness of the "skin" of an atom. This is crucial for the PUMA experiment at CERN, which wants to study rare isotopes by seeing how big their "skins" are.

4. The Comparison: The "Exact vs. Approximate" Test

To make sure their new map was right, they compared it against the "Gold Standard" (exact mathematical solutions used for very small systems).

• The Verdict: Their map matched the Gold Standard very well for the smallest cities (Deuterium). However, they noticed small differences.
• The Reason: The Gold Standard accounts for every single citizen in the city rearranging themselves instantly. Their map assumes the city stays mostly in one shape while the ghost visits.
• The Takeaway: For the tiny cities, the "shape-shifting" matters a bit. But for larger cities (which they plan to study next), their method is expected to be even better because the city is too big to rearrange instantly anyway.

Summary

In short, these scientists built a new, super-powerful simulation to predict how antimatter behaves near atomic nuclei.

1. They fixed a math glitch caused by the violent nature of antimatter using a "smoothing filter."
2. They proved that antimatter probes are excellent for measuring the "skin" of atoms because they only interact with the edges.
3. They laid the groundwork for future experiments at CERN that will use these "ghost visitors" to map the structure of rare and exotic atoms, helping us understand the universe's building blocks better.

It's like upgrading from a blurry, shaky video of a ghost to a crystal-clear 4K movie, allowing us to finally see exactly where the ghost walks and how it changes the house it visits.

Technical Summary

1. Problem Statement

The paper addresses the theoretical challenge of describing antiproton-nucleus ($\bar{p}$-nucleus) systems at low energies. While experimental facilities like the CERN Antiproton Decelerator (specifically the PUMA experiment) are generating new data on antimatter interactions with nuclei, theoretical frameworks for these systems remain underdeveloped compared to standard nucleon-nucleus physics.

Key challenges include:

• Annihilation: The $\bar{p}$-nucleon interaction involves strong annihilation, converting mass into mesons (primarily pions). This requires a non-Hermitian (complex) Hamiltonian.
• Multi-scale nature: The systems involve nuclear scales (hundreds of MeV for quasi-bound states) and atomic scales (eV for antiprotonic atom level shifts).
• Interaction Complexity: The Nucleon-Antinucleon ($N\bar{N}$) interaction is "hard" (strong short-range components), leading to slow convergence in standard basis expansions.
• Lack of Ab Initio Methods: Existing calculations often rely on phenomenological models or few-body methods (like Faddeev) limited to very light systems ($A \le 3$). There is a need for a systematic ab initio approach that can eventually scale to heavier nuclei ($A \sim 16$).

2. Methodology: NCSM/RGM Adaptation

The authors extend the No-Core Shell Model combined with the Resonating Group Method (NCSM/RGM) to antiprotonic systems.

• Formalism Adaptation:
  • The standard NCSM/RGM uses an antisymmetrized basis to enforce the Pauli exclusion principle between all nucleons.
  • For $\bar{p}$-nucleus systems, the projectile is an antiproton (not a nucleon). Therefore, antisymmetrization between the target and the projectile is not required. This significantly simplifies the formalism: the inter-cluster antisymmetrizer operator ($\hat{A}$) is omitted, and the norm kernel becomes trivial (diagonal).
• Hamiltonian:
  • The system is treated using a non-relativistic Schrödinger equation with a complex optical potential $V = U + W$.
  • Real Part ($U$): Derived from meson-exchange models (specifically the Kohno-Weise potential) using G-parity transformation from nucleon-nucleon interactions.
  • Imaginary Part ($W$): A Woods-Saxon potential representing the annihilation channel.
  • Coulomb Interaction: Both the point-Coulomb interaction (between $\bar{p}$ and individual protons) and the average Coulomb potential (between $\bar{p}$ and the target center of mass) are included.
• Basis Expansion:
  • The wave function is expanded in a Harmonic Oscillator (HO) basis using Jacobi coordinates.
  • The target nucleus wave functions are obtained from NCSM calculations using chiral EFT interactions (N3LO).
  • The relative motion between the $\bar{p}$ and the target is expanded in an HO basis truncated at $N_{max}$.

3. Key Technical Contributions & Innovations

• Handling "Hard" Potentials: The study identifies that the hard short-range nature of the $N\bar{N}$ interaction causes finite-model-space artifacts. In a truncated HO basis, the potential kernels exhibit spurious oscillations at large distances, leading to unstable phase shifts.
• Regulator Strategy: To mitigate these artifacts without requiring computationally prohibitive model spaces, the authors introduce Woods-Saxon type regulators ($f_s(r)$ and $f_c(r)$). These smooth functions suppress the unphysical long-range tails of the HO basis functions while preserving the short-range interaction physics. This allows for rapid convergence of observables.
• Multi-Channel Coupling: The method incorporates target excited states (pseudo-states) to approximate continuum coupling, although the study finds that for low energies, the ground state dominates.
• Unified Framework: The same wave-function ansatz is used to calculate both nuclear properties (scattering, quasi-bound states) and atomic properties (level shifts, widths), bridging the gap between nuclear and atomic physics scales.

4. Results

The authors applied the method to the lightest systems: $\bar{p} + d$ (deuteron), $\bar{p} + ^3H$ (triton), and $\bar{p} + ^3He$.

• Convergence:
  • With the regulator, phase shifts and scattering lengths converge rapidly with increasing $N_{max}$.
  • The inclusion of target pseudo-states (excited/discretized continuum states) has a minor effect (few percent) on scattering lengths, confirming that low-energy dynamics are dominated by the ground state.
• Scattering Observables:
  • Phase Shifts & Scattering Lengths: Calculated for various partial waves ($S, P, D$). The imaginary parts of the scattering lengths are significant due to annihilation.
  • Cross Sections: Reaction and elastic differential cross sections were computed and compared with experimental data (where available). The results show reasonable agreement but highlight discrepancies likely due to the specific $N\bar{N}$ interaction model used.
• Antiprotonic Atoms:
  • Level Shifts and Widths: Calculated for the $1S$ and $2S$ atomic states.
  • Comparison: Results for $\bar{p}-d$ and $\bar{p}-^3He$ were compared with independent Faddeev calculations (exact few-body solutions).
    • For $\bar{p}-d$, there is a discrepancy with Faddeev results, attributed to the weak binding of the deuteron where the cluster approximation ($\bar{p} + d$) may break down compared to a full three-body treatment.
    • For $\bar{p}-^3He$, the agreement is better, though deviations (2–15%) persist, suggesting limitations in the two-body $N\bar{N}$ interaction model.
• Quasi-Bound States:
  • The study predicts nuclear quasi-bound states (deeply bound states with large imaginary widths due to annihilation).
  • Energies for $\bar{p}-d$, $\bar{p}-^3H$, and $\bar{p}-^3He$ are reported. These states are highly sensitive to the $N\bar{N}$ interaction model.
• Annihilation Density:
  • Calculations of the annihilation density $\gamma_a(r)$ confirm that for these light systems, annihilation occurs predominantly at the nuclear periphery (tail of the density distribution, $r \approx 2$ fm). This validates the use of antiprotons as probes for nuclear surface properties.

5. Significance and Outlook

• Validation of Method: The work demonstrates that the NCSM/RGM is a robust, ab initio tool for antinucleon-nucleus systems. The simplification of removing the inter-cluster antisymmetrizer makes it computationally efficient.
• Foundation for Heavier Nuclei: While currently limited to light nuclei due to the "hard" potential requiring large model spaces, the method is designed to scale. Future work will apply the Slater-Determinant (SD) version of NCSM/RGM to heavier targets ($A \sim 16$) using softer interactions (e.g., via Similarity Renormalization Group).
• Experimental Support: The results provide critical theoretical benchmarks for the PUMA experiment at CERN. By quantifying the uncertainties related to the $N\bar{N}$ interaction and missing configurations, the paper helps interpret experimental data regarding nuclear surface densities and neutron skins.
• Future Directions: The framework can be extended to antideuteron projectiles and systems containing hyperons, offering a unified approach to exotic hadron-nucleus interactions.

In summary, this paper successfully adapts a state-of-the-art nuclear many-body method to antimatter systems, overcoming numerical challenges posed by annihilation and hard potentials, and providing a comprehensive set of predictions for light antiprotonic nuclei.