Calculations of Di-Hadron Production via Two-Photon Processes in Relativistic Heavy-Ion Collisions

Titel: Berekeningen van Di-hadronproductie via Twee-fotonprocessen in Relativistische Zware-ionenbotsingen

Auteurs: Luobing Wang, Xinbai Li, Zebo Tang, Xin Wu en Wangmei Zha (USTC, China)

1. Probleemstelling en Context

Relativistische zware-ionenbotsingen (bij RHIC en LHC) worden voornamelijk gebruikt om het quark-gluonplasma (QGP) te bestuderen. In ultra-perifere botsingen (UPC), waarbij de impactparameter groter is dan de som van de kernstralen, worden hadronische interacties onderdrukt en domineren elektromagnetische processen. Hoewel twee-fotonprocessen ($\gamma\gamma$) die leiden tot de productie van leptonparen (dileptonen, zoals $e^+e^-$ en $\mu^+\mu^-$) uitgebreid zijn bestudeerd en goed overeenkomen met kwantumelektrodynamica (QED)-voorspellingen, ontbreekt er een kwantitatieve theoretische basis voor de productie van hadronparen (mesonen en baryonen) via deze mechanismen.

Hoewel er recente experimentele observaties zijn van exclusieve proton-antiproton ($p\bar{p}$) productie bij STAR en $K^+K^-$ productie bij LHC, zijn er slechts weinig theoretische voorspellingen beschikbaar om deze metingen te benchmarken. Een specifiek probleem is de extrapolatie van data van $e^+e^-$-colliders (waar de fotonen een virtueelheid $Q^2 \sim 10^{-2}$ GeV$^2$ hebben) naar zware-ionen-UPC (waar de fotonen quasi-reëel zijn met $Q^2 \lesssim 10^{-3}$ GeV$^2$). De auteurs willen deze kennislacune opvullen door een geünificeerde berekening te leveren voor $\pi^+\pi^-$, $K^+K^-$ en $p\bar{p}$ productie.

2. Methodologie

De auteurs hanteren de Equivalent Photon Benadering (EPA) om de productie van di-hadronen in ultra-perifere botsingen te modelleren. De kern van de methode bestaat uit de volgende stappen:

• Fotonfluxberekening: De flux van quasi-reële fotonen wordt berekend door de elektromagnetische velden van de bewegende kernen te behandelen als een bundel virtuele fotonen.
  • Voor afstanden kleiner dan de kernstraal ($r < R$) wordt de ladingsverdeling gemodelleerd met een Woods-Saxon-verdeling (twee-parameter Fermi-verdeling), wat rekening houdt met de nucleaire structuur.
  • Voor grotere afstanden wordt een puntlading-benadering gebruikt.
  • De flux wordt geïntegreerd over de impactparameter $b$ en de transversale coördinaten.
• Hadronische onderdrukking: Om te garanderen dat alleen echte UPC-events worden geteld (zonder kern-kern botsingen), wordt de waarschijnlijkheid $P_{NH}(b)$ berekend dat er geen hadronische interactie optreedt. Dit wordt gedaan met het Glauber-model, waarbij de inelastische nucleon-nucleon doorsnede ($\sigma_{inel}^{NN}$) en de kern-overlappingsfunctie worden gebruikt.
• Mutuele Coulomb-excitatie (MCE): De auteurs houden rekening met de kans dat de kernen worden geëxciteerd (en neutronen uitzenden) door de intense fotonvelden. Dit wordt gemodelleerd via Poisson-statistieken gebaseerd op de Giant Dipole Resonance (GDR).
• Elementaire Doorsnede ($\sigma_{\gamma\gamma \to h\bar{h}}$): De kern van de berekening is het "vouwen" (folding) van de berekende fotonflux met de elementaire doorsnede voor $\gamma\gamma \to \text{hadronen}$. Deze elementaire doorsnedes worden direct afgeleid van bestaande experimentele data uit $e^+e^-$-colliders (zoals Belle, BaBar, LEP) voor de kanalen:
  • $\gamma\gamma \to \pi^+\pi^-$
  • $\gamma\gamma \to K^+K^-$
  • $\gamma\gamma \to p\bar{p}$
• Kinematische Acceptatie: De berekeningen worden gefilterd door de specifieke acceptatiecriteria van de STAR-detector (voor Au+Au bij $\sqrt{s_{NN}} = 200$ GeV) en typische LHC-detectors (voor Pb+Pb bij $\sqrt{s_{NN}} = 5.36$ TeV). Dit omvat snijwaarden voor rapiditeit ($Y$), transversale impuls ($p_T$) en pseudorapidity ($\eta$).
• Monte Carlo Simulatie: Een MC-simulatie wordt gebruikt om de transversale impulsverdelingen te genereren en de acceptatiecorrecties toe te passen, inclusief extrapolaties voor hoekverdelingen buiten het bereik van de $e^+e^-$-data (waar $|\cos \theta^*| > 0.6$).

3. Belangrijkste Bijdragen

• Unificatie van Baselines: Het artikel levert de eerste uitgebreide, geünificeerde theoretische baselines voor de productie van licht-mesonparen ($\pi, K$) en baryonparen ($p$) via twee-fotonprocessen in zowel RHIC- als LHC-omgevingen.
• Validatie van de EPA voor Hadronen: Het toont aan dat de EPA, gecombineerd met $e^+e^-$-data, een robuuste voorspelling kan doen voor hadronische eindtoestanden, ondanks de verschillen in fotonvirtueleheid tussen $e^+e^-$ en zware-ionen-systemen.
• Vergelijking met Bestaande Theorie: De resultaten voor $p\bar{p}$ productie bij RHIC worden vergeleken met eerdere theoretische werken. De auteurs tonen aan dat hun resultaten uitstekend overeenkomen met de voorspellingen van Shao et al., terwijl andere studies (zoals Pu et al.) waarden geven die twee tot drie orden van grootte te hoog zijn.
• Systeematische Onzekerheidsanalyse: Een gedetailleerde analyse van onzekerheden wordt gepresenteerd, waaronder nucleaire parameters, fouten in de $e^+e^-$-data en extrapolaties van hoekverdelingen.

4. Resultaten

De berekende productiedoorsneden tonen duidelijke trends:

• Gradiënt in Doorsnede: Er is een duidelijke hiërarchie in de productie-waarschijnlijkheid, bepaald door de massa van de deeltjes:
  $$\sigma_{\pi^+\pi^-} \gg \sigma_{K^+K^-} \gg \sigma_{p\bar{p}}$$
  Elke volgende fase is ongeveer een orde van grootte onderdrukt ten opzichte van de vorige.
• Energie-afhankelijkheid:
  • STAR (Au+Au, 200 GeV): De doorsneden liggen in de orde van microbarn ($\mu$b). De detectoracceptatie is ongeveer 60% voor $p\bar{p}$ en 85% voor mesonen.
  • LHC (Pb+Pb, 5.36 TeV): De doorsneden liggen in de orde van millibarn (mb), wat een toename van ongeveer drie orden van grootte betekent ten opzichte van RHIC. Dit komt door de hogere botsingsenergie en het ontbreken van de strenge "neutron-trigger" vereisten die bij STAR worden gebruikt.
• Transversale Impuls ($p_T$): De $p_T$-spectra van de geproduceerde paren pieken scherp onder 0.15 GeV/c. Dit is een kenmerkend signaal van coherente twee-fotonproductie, waarbij de transversale impuls wordt bepaald door de fotonen zelf en niet door harde botsingen.
• Onzekerheden: De totale systematische onzekerheid wordt geschat op < 4% voor STAR-kinematica en < 10% voor LHC-kinematica. De dominante bron van onzekerheid is de extrapolatie van hoekverdelingen, vooral bij LHC waar een groter deel van de fase-ruimte ($|\cos \theta^*| > 0.6$) wordt bestreken.

5. Betekenis en Conclusie

Dit werk biedt een cruciale theoretische referentie voor de interpretatie van toekomstige metingen bij STAR en LHC.

• Benchmarking: De voorspellingen dienen als directe benchmarks om de nieuwste experimentele data (zoals de recente STAR $p\bar{p}$ resultaten) te valideren.
• Test van QED in Sterke Velden: Het bevestigt dat de intense elektromagnetische velden van relativistische ionen effectief kunnen worden gebruikt om QED-processen te testen in het regime van sterke velden, zelfs voor hadronische eindtoestanden.
• Fysica van de Fotonvirtueleheid: Door de resultaten te vergelijken met $e^+e^-$-data, biedt het inzicht in de invloed van fotonvirtueleheid op de productie van baryonparen, wat beperkingen oplegt aan elektromagnetische interacties en de structuur van baryonen in het regime van bijna-reële fotonen.

Samenvattend vult dit onderzoek een belangrijke leemte in de theoretische fysica van zware-ionenbotsingen op en legt het de basis voor systematische onderzoeken naar hadronische twee-fotonprocessen in de komende jaren.