Shedding Light on (Anti-)nuclei Production with Pion-Nucleus Femtoscopy

Technical Summary: Shedding Light on (Anti-)nuclei Production with Pion-Nucleus Femtoscopy

Problem Statement
High-energy nuclear collisions create an environment where light nuclei and antinuclei (e.g., deuterons, $^4$He) are synthesized at temperatures ($k_B T \sim 100$ MeV) significantly exceeding their binding energies per nucleon (a few MeV). Despite decades of study, the microscopic mechanism governing the formation and survival of these fragile (anti-)nuclei in such hot, dense hadronic environments remains unresolved. Three primary competing scenarios exist:

1. Statistical Hadronization Model (SHM): Nuclei form at chemical freeze-out from a thermal source.
2. Nucleon Coalescence: Nucleons combine at kinetic freeze-out based on phase-space proximity.
3. Pion-Catalyzed Reactions: Continuous regeneration and disintegration occur in the hadronic phase via reactions like $\pi NN \leftrightarrow \pi d$, potentially mediated by intermediate $\Delta(1232)$ resonances.

Bulk observables (yields, $p_T$ spectra, collective flow) have proven insufficient to distinguish these models, as different mechanisms can reproduce similar integrated yields. While correlation femtoscopy (HBT interferometry) has emerged as a precision tool for probing spatio-temporal structures and strong interactions, its application to hadron-nucleus systems (specifically $\pi-d$) to discriminate production mechanisms requires further theoretical clarification, particularly regarding recent ALICE observations of resonance peaks in $\pi^\pm-p$ and $\pi^\pm-d$ correlations.

Methodology
The authors employ a multi-phase transport (AMPT) model to simulate the partonic evolution and hadronization (via quark coalescence) of high-multiplicity $pp$ collisions at $\sqrt{s} = 13$ TeV. For the subsequent hadronic evolution, they solve relativistic kinetic equations for pion-catalyzed reactions ($\pi NN \leftrightarrow \pi d$) to describe deuteron production and dissociation.

Key methodological components include:

• Kinetic Equations: The phase-space distribution of deuterons is evolved using a Boltzmann-type equation incorporating dissociation ($K_>$) and regeneration ($K_<$) rates. The reaction probability is calculated using the impulse approximation (IA), modified to include a Gaussian smearing function to allow for event-by-event simulations ($N_{test}=1$) necessary for femtoscopy, avoiding artificial suppression of correlations.
• Cross-Sections and Interactions: The $\pi^+ - p$ cross-section is dominated by $p$-wave interactions through the $\Delta$ resonance, calculated using a $\pi N \Delta$ interaction Lagrangian. The $\pi^+ - d$ interaction is treated via the impulse approximation, relating it to the $\pi^+ - p$ cross-section with appropriate statistical factors and Wigner functions.
• Correlation Function Calculation: Two-particle correlation functions $C(k^*)$ are computed using the relative momentum $k^*$ and spatial separation $r^*$ in the pair rest frame. The calculation includes final-state interactions (FSI), specifically $p$-wave scatterings and Coulomb corrections, derived from scattering wave functions.
• Comparative Frameworks: The authors also calculate predictions using the nucleon coalescence model (at kinetic freeze-out) and the statistical hadronization model (thermal emission) to serve as benchmarks against the kinetic approach.

Key Results
The study successfully reproduces the experimental correlation functions measured by the ALICE Collaboration for $\pi^+ - p$ and $\pi^+ - d$ pairs in high-multiplicity $pp$ collisions.

• Resonance Peaks: The theoretical model predicts pronounced peaks in the correlation functions corresponding to the $\Delta$ resonance mass.
• Peak Shift Mechanism: Without final-state $p$-wave scatterings, the theoretical peaks appear at $k^* \approx 0.22$ GeV/c (vacuum $\Delta$ decay position). However, the interplay between the $\Delta$ resonance and $p$-wave scatterings shifts both correlation peaks downward to $k^* \approx 0.15$ GeV/c. This shift aligns precisely with the ALICE data.
• Model Discrimination:
  • The kinetic approach (pion-catalyzed reactions) successfully reproduces the position and shape of the experimental peaks.
  • The nucleon coalescence model significantly underestimates the data and produces an unobserved dip at high $k^*$ due to insufficient compensation from $\Delta$-decay contributions against $p$-wave interference.
  • The statistical hadronization model also significantly underestimates the data and generates a pronounced dip, lacking the resonance contributions inherent in the kinetic mechanism.

Significance and Claims
The paper claims that pion-nucleus femtoscopy possesses unique discriminating power for identifying the production mechanisms of light (anti-)nuclei in high-energy nuclear collisions. The results provide compelling evidence that pion-catalyzed reactions (specifically $\pi NN \leftrightarrow \pi d$) play the dominant role in the production of light nuclei in these collisions, rather than statistical hadronization or simple nucleon coalescence at freeze-out.

The authors assert that their approach elucidates the microscopic origin of the resonance structures observed by ALICE, resolving the previous ambiguity regarding the "unexpectedly low" freeze-out temperatures suggested by the peak shifts. Furthermore, the paper suggests that extending this femtoscopic method to heavier nuclei (e.g., $\pi$-$^3$He, $\pi$-$\alpha$) and hypernuclei, as well as to exotic hadrons like $X(3872)$ and $T_{cc}^+$, could provide crucial insights into their production mechanisms in similar high-energy environments.