Bridging reaction theory and nuclear structure in $π^\pm$-${}^{48}$Ca scattering

This paper extends the pion-nucleus multiple-scattering framework to include second-order rescattering dynamics and nuclear structure details derived from chiral effective field theory, demonstrating that these corrections are essential for accurately reproducing differential cross sections in $\pi^\pm$-${}^{48}$Ca elastic scattering within the $\Delta(1232)$-resonance region.

The Gist

Imagine the atomic nucleus not as a solid marble, but as a bustling, crowded dance floor filled with tiny dancers (protons and neutrons). Now, imagine shooting a fast-moving pi-meson (a type of subatomic particle) at this dance floor. What happens? The pi-meson doesn't just bounce off the edge; it dives into the crowd, bumps into dancers, gets knocked around, and might even swap partners before finally exiting.

This paper is about building a better map to predict exactly how that pi-meson bounces off a specific, crowded dance floor: the Calcium-48 nucleus.

Here is the story of their work, broken down into simple concepts:

1. The Problem: The "Crowded Dance Floor" is Different

Scientists have been studying how particles bounce off nuclei for a long time. They were very good at predicting what happens when the dance floor is perfectly balanced (equal numbers of protons and neutrons). But Calcium-48 is unbalanced; it has more neutrons than protons. It's like a dance floor where one group of dancers is much larger than the other.

Previous maps (theories) worked well for balanced floors but struggled with unbalanced ones because they didn't account for the specific "charge-swapping" moves that happen when the extra neutrons get involved.

2. The New Map: Adding "Second-Order" Moves

The authors created a new, more detailed map. They realized that to get the prediction right, you can't just look at the first bump. You have to look at the second bump.

• First-Order (The Simple Bounce): The pi-meson hits one dancer and bounces off.
• Second-Order (The Complex Shuffle): The pi-meson hits one dancer, which excites the whole floor. Then, before the pi-meson leaves, it hits a second dancer. Crucially, during this time, the two dancers might swap roles (a proton becomes a neutron and vice versa) or flip their spins.

The authors built a mathematical "potential" (a set of rules for how the pi-meson moves) that includes these complex, two-step shuffles. They found that ignoring these second-order moves is like trying to predict a dance by only watching the first step; you miss the most important part of the routine.

3. The Ingredients: How They Built the Map

To make this map accurate, they needed two specific ingredients:

• The Dancers' Positions (One-Body Density): They used a super-advanced computer method called "Coupled-Cluster Theory" to figure out exactly where the protons and neutrons are sitting in the Calcium-48 nucleus. Think of this as a high-resolution 3D scan of the dance floor.
• The Dancers' Relationships (Two-Body Correlations): They needed to know how the dancers relate to each other. If one moves, how does the neighbor react? They used a slightly simpler method called "Hartree-Fock" to map these relationships.

They tested their map using two different sets of "rules of physics" (called Chiral Effective Field Theory interactions). It's like testing a navigation app with two different map providers. They found that while the details of the dance floor changed slightly depending on which provider they used, the final prediction for how the pi-meson bounces remained surprisingly stable.

4. The Results: The Map Works

They tested their new map against real-world data collected from experiments where scientists actually shot pions at Calcium-48.

• The "Delta" Zone: They focused on a specific energy range (the $\Delta(1232)$ resonance) where the pi-meson and the nucleus interact most strongly, like a dance move that gets everyone excited.
• The Verdict: When they included the "second-order" complex shuffles, their predictions matched the experimental data almost perfectly.
  • If they only used the simple "first-order" bounce, the prediction was off.
  • Once they added the complex, two-step interactions, the curve of the data fit beautifully.

5. Why It Matters (According to the Paper)

The paper claims this work is a bridge. It connects the theory of how nuclei are built (nuclear structure) with the theory of how particles crash into them (reaction theory).

They also noted that while their model works great for Calcium-48, there are still some small mismatches with data from a specific experiment at 130 MeV for negative pions. However, they suggest this might be an issue with the experimental data itself rather than their theory, especially since their model works well for other energies and for a similar nucleus (Calcium-40).

In a nutshell: The authors built a sophisticated, two-step simulation of how a particle bounces off an unbalanced atomic nucleus. By accounting for the complex "dance" between pairs of protons and neutrons, they created a model that accurately predicts real-world experimental results, proving that you can't understand the bounce without understanding the shuffle.

Technical Summary

Problem Statement
The study of pion-nucleus interactions remains critical for modern nuclear physics, particularly for modeling final-state interactions in neutrino-nucleus scattering experiments (e.g., DUNE, MicroBooNE) and for extracting neutron skin thicknesses to constrain the nuclear equation of state. While previous theoretical frameworks successfully described pion scattering off spin-isospin-zero nuclei (such as $^{12}\text{C}$, $^{16}\text{O}$, and $^{40}\text{Ca}$) using a multiple-scattering approach with second-order corrections, these models relied on isospin invariance. This assumption breaks down for nuclei with non-zero isospin, such as $^{48}\text{Ca}$. In such systems, intermediate charge-exchange processes and nucleon spin-flip effects become significant, requiring a more rigorous treatment of the scattering potential that explicitly accounts for the distinct distributions of protons and neutrons and their correlations.

Methodology
The authors extend the Kerman-McManus-Thaler multiple-scattering formalism to derive a second-order pion-nucleus scattering potential suitable for nuclei with non-zero isospin. The methodology involves several key components:

1. Scattering Potential Derivation: The elastic scattering amplitude is treated as a solution to a Lippmann-Schwinger-type integral equation. The scattering potential $V$ is approximated by the sum of first-order ($V^{[1]}$) and second-order ($V^{[2]}$) terms.

   • The first-order potential depends on the one-body densities (form factors) of protons and neutrons.
   • The second-order potential incorporates intermediate nuclear excitations, specifically accounting for pion rescattering on nucleon pairs. For non-zero isospin nuclei, this derivation explicitly separates contributions from proton-proton ($pp$), neutron-neutron ($nn$), and neutron-proton ($np$) pairs, including an exchange correlation function ($C^{(ex)}$) that describes intermediate charge-exchange processes.

2. Nuclear Structure Inputs:

   • One-body densities: Calculated using Coupled-Cluster (CC) theory at the CCSDT-1 level (including singles, doubles, and leading triples excitations).
   • Two-body densities: Computed within the Hartree-Fock (HF) approximation. The authors justify this simplification by noting that HF reasonably describes one-body form factors and is computationally less demanding than full CC calculations for two-body matrices, while the resulting errors are expected to be small.
   • Hamiltonians: The calculations utilize modern nuclear Hamiltonians derived from chiral effective field theory ($\chi$EFT), specifically the $\text{NNLO}_{\text{sat}}(450)$ and $\Delta\text{NNLO}_{\text{GO}}(394)$ interactions, to estimate theoretical uncertainties related to the choice of nuclear force.

3. Scattering Amplitudes: The elementary pion-nucleon scattering amplitudes are taken from the SAID phase-shift analysis (WI08). Their in-medium modifications are characterized by three energy-independent parameters (real and imaginary parts of the $\Delta(1232)$ self-energy and the slope of the $s$-wave isoscalar amplitude), which were previously fitted to $\pi^\pm$-$^{12}\text{C}$ data and are applied here without further tuning.

4. Electromagnetic Effects: The model includes short- and long-range Coulomb interactions using a model-independent Fourier-Bessel series for the $^{48}\text{Ca}$ charge distribution.

Key Contributions

• Generalization of Reaction Theory: The paper provides the first derivation of the second-order pion-nucleus scattering potential for targets with non-zero isospin, explicitly handling the spin-isospin structure of intermediate charge-exchange and spin-flip mechanisms.
• Integration of Ab Initio Structure: The work bridges reaction theory and nuclear structure by employing state-of-the-art ab initio many-body methods (CC and HF) to generate the necessary one- and two-body density inputs, rather than relying on fitted phenomenological densities.
• Uncertainty Quantification: By comparing results from two distinct chiral interactions ($\text{NNLO}_{\text{sat}}$ and $\Delta\text{NNLO}_{\text{GO}}$), the authors quantify the sensitivity of scattering observables to the underlying nuclear structure and the inclusion of explicit $\Delta$-isobars in the chiral expansion.

Results
The model was applied to $\pi^\pm$-$^{48}\text{Ca}$ elastic scattering in the $\Delta(1232)$-resonance region (laboratory kinetic energies 116–230 MeV).

• Second-Order Effects: The inclusion of the second-order potential ($V^{[2]}$) is found to be essential. It significantly improves the reproduction of differential cross-section magnitudes, particularly at lower energies, though it causes only minor shifts in the positions of diffraction minima.
• Comparison with Data: The theoretical predictions show good agreement with experimental data from SIN and LAMPF. The model successfully reproduces the general features of angular distributions for both $\pi^+$ and $\pi^-$.
• Discrepancies: Notable deviations persist for $\pi^-$ scattering at 130 MeV. However, the authors note that data at 116 MeV and 180 MeV are well-described, and $\pi^\pm$ scattering on $^{40}\text{Ca}$ at 130 MeV is also well-reproduced. This suggests the discrepancy at 130 MeV for $\pi^-$ may stem from inconsistencies in the experimental data rather than a failure of the theoretical framework.
• Sensitivity: The scattering observables show mild sensitivity to the choice of nuclear interaction (the difference between $\text{NNLO}_{\text{sat}}$ and $\Delta\text{NNLO}_{\text{GO}}$). The uncertainty associated with the second-order potential arises primarily from differences in the predicted two-body correlation functions, which become more pronounced at larger scattering angles (higher momentum transfer).
• Normalization: When systematic normalization uncertainties in the experimental data are accounted for, the agreement between the model and data improves further, with normalization shifts remaining within the reported experimental uncertainties.

Significance
The paper claims that this work establishes a complete modeling framework for coherent pion photoproduction off $^{48}\text{Ca}$ by setting the stage for accurate descriptions of pion scattering on non-zero isospin nuclei. It demonstrates that second-order rescattering corrections are sizeable and necessary for accurate cross-section predictions in the $\Delta$-resonance region. Furthermore, the successful application of this formalism to $^{48}\text{Ca}$ validates its potential utility for analyzing coherent pion photoproduction and scattering on other relevant targets, such as $^{40}\text{Ar}$, which is crucial for interpreting data from long-baseline neutrino oscillation experiments. The study reinforces the consistency of the approach used for zero-spin/isospin nuclei when extended to more complex nuclear systems.