The next-to-leading order of the differential cross section of the subprocess of annihilation of quark-antiquark pair of prompt photon production in proton-proton collisions at NICA energies

Technische Samenvatting: Differentieel Cross-sectie op Volgende-na-Leidende Orde voor Kwart-Antikwart Annihilatie in Prompt Foton Productie bij NICA-energieën

Probleemstelling
De productie van prompt fotonen in proton-proton botsingen dient als een kritieke sonde voor het begrijpen van de dynamische structuur van het proton en de vorming van quark-gluonplasma (QGP). Hoewel theoretische onderzoeken historisch gezien gericht zijn geweest op hoge energie-regimes (LHC, Tevatron), opereert de Nuclotron-gebaseerde Ion Collider fAcility (NICA) bij een lagere middelpunt-van-het-massasnelheid ($\sqrt{s} = 10$ GeV). Bij deze energieën is de productie van additionele elementaire deeltjes beperkt, wat de bepaling van effectieve cross-secties vereenvoudigt, maar precieze theoretische instrumenten vereist om verschillende productiemechanismen te onderscheiden. Voorgaand werk door de auteurs stelde Leading Order (LO) berekeningen vast voor Compton-verstrooiing, annihilatie en bremsstrahlung-subprocessen. Echter, LO-benaderingen zijn vaak onvoldoende om experimentele data te beschrijven, met name bij grote transversale impulsen, wat hogere-orde correcties noodzakelijk maakt. Deze studie adresseert de behoefte aan een volledige Next-to-Leading Order (NLO) berekening specifiek voor het kwart-antikwart annihilatie-subproces ($q\bar{q} \to \gamma g$) in zowel niet-gepolariseerde als longitudinaal gepolariseerde proton-proton botsingen bij NICA-energieën.

Methodologie
Het onderzoek hanteert een duaal benadering die analytische fixed-order berekeningen combineert met Monte Carlo (MC) event generatie:

1. Analytische NLO Berekening:

   • Diagram Generatie: Met behulp van FeynCalc werden initieel 773 one-loop diagrammen gegenereerd voor het $q\bar{q} \to \gamma$ subproces. Na het toepassen van selectiecriteria—het uitsluiten van zware deeltjes die niet geregistreerd worden bij NICA, driehoekige loops die verdwijnen door de Furry-stelling, "seagull" diagrammen, uitsluitingen van mass-shell renormalisatie en niet-kleur-conserverende transities—werden 13 bijdragende diagrammen behouden.
   • Regularisatie en Renormalisatie: Berekeningen werden uitgevoerd met dimensionale regularisatie binnen het Minimal Supersymmetric Standard Model (MSSM) schema. De sterke koppelingsconstante ($\alpha_s$) en renormalisatieschalen werden gedefinieerd, waarbij de factorisatieschaal gelijk werd gesteld aan de partonische middelpunt-van-het-massasnelheid ($\hat{s}$).
   • Polarisatie: Longitudinale polarisatie werd geïntegreerd met behulp van gepolariseerde Parton Distribution Functions (PDFs) en heliciteit-afhankelijke matrixelementen. De dubbele spin-asymmetrie ($A_{LL}$) werd berekend op basis van cross-secties voor parallel ($\uparrow\uparrow$) en anti-parallel ($\uparrow\downarrow$) proton-polarisaties.
   • PDFs: Het LHAPDF 6.5.5 framework werd gebruikt om NLO PDFs (specifiek CT14nlo) te integreren voor hadron-niveau cross-sectie berekeningen.

2. Monte Carlo Simulatie:

   • Framework: De PYTHIA 8.316 generator werd geïnterface met de POWHEG BOX (V2) om prompt foton productie te simuleren. Deze NLO+PS (Parton Shower) benadering koppelde NLO matrixelementen aan parton showering en hadronisatie.
   • Configuratie: Simulaties werden uitgevoerd voor $\sqrt{s} = 10$ GeV met een transversale impuls-cut ($p_T^{min} = 0.5$ GeV/c) om infrarood divergenties te reguleren. De CT14nlo PDF set werd gebruikt voor de initiële-toestand hadronstructuur.
   • Event Selectie: Prompt fotonen werden geïdentificeerd via PDG code 22 met status code 23, waarbij fiduciale cuts op $p_T$ en rapiditeit ($y$) werden toegepast.

Belangrijkste Bijdragen en Resultaten
Het artikel presenteert de eerste volledige NLO berekeningen voor het $q\bar{q}$ annihilatie-subproces in de context van NICA-energieën, waarbij deze worden vergeleken met LO resultaten en MC simulaties.

• Grootte van NLO Correcties: De studie toont aan dat NLO bijdragen significant zijn bij hoge botsingsenergieën. De NLO differentiële cross-sectie draagt ongeveer 15% bij relatief aan de LO berekening bij NICA-energieën. Dit staat in contrast met hogere-energie omgevingen (LHC/Tevatron), waar NLO bijdragen 35% kunnen bereiken.
• Kinematische Afhankelijkheid:
  • Energie ($\sqrt{s}$): De totale cross-sectie piekt bij $\sqrt{s} \approx 5.2$ GeV en neemt af bij hogere energieën door het lopen van de QCD koppelingsconstante en het gedrag van de PDF.
  • Transversale Impuls ($p_T$): De differentiële cross-sectie ($d\sigma/dp_T$) neemt scherp af met toenemende $p_T$. De ratio van NLO ten opzichte van LO correcties neemt toe met $p_T$, en bereikt tot wel 30% bij grote transversale impulsen.
  • Angulaire en Rapiditeit Distributies: De angulaire distributie ($d\sigma/d\cos\theta$) vertoont een symmetrische piekvorm met een maximum bij $\cos\theta = 0$ (transversaal) en minima bij $\cos\theta = \pm 1$, gedreven door $t$-kanaal virtuele foton uitwisseling. De rapiditeitsdistributie vertoont een maximum bij $|y| \approx 1.95$.
• Polarisatie Effecten:
  • Longitudinale polarisatie verandert de vorm van de kinematische distributies niet, maar beïnvloedt de grootte van de cross-secties aanzienlijk.
  • Asymmetrie: De dubbele spin-asymmetrie ($A_{LL}$) bereikt een maximaal positieve waarde ($+1$) wanneer kwark- en antikwarks spins parallel zijn en een maximaal negatieve waarde ($-1$) wanneer ze anti-parallel zijn. De studie vindt dat de invloed van longitudinale polarisatie significanter is voor NLO berekeningen dan voor LO, met name bij hoge energieën. Dit wordt toegeschreven aan de gevoeligheid van virtuele loop-correcties en hogere-orde parton verstrooiingsconfiguraties voor initiële spin-oriëntaties.
• Vergelijking Analytisch vs. Simulatie:
  • De ratio van de analytische (FeynCalc) naar de gesimuleerde (PYTHIA) cross-secties varieert met de kinematica. Bij lage $p_T$ onderschat de analytische fixed-order berekening de cross-sectie vergeleken met PYTHIA, vanwege logaritmische divergenties in de soft/collineaire limiet, welke PYTHIA oplost via Sudakov ressommatie en Initial State Radiation (ISR).
  • Bij hoge $p_T$ nadert de ratio eenheid, wat duidt op convergentie tussen de fixed-order theorie en de NLO+PS simulatie.

Betekenis en Claims
De auteurs stellen dat dit werk een noodzakelijke theoretische verfijning biedt voor het analyseren van prompt foton productie bij NICA-energieën. De primaire betekenis ligt in de demonstratie dat NLO correcties niet verwaarloosbaar zijn (ca. 15%) zelfs bij deze gematigde energieën, wat het gebruik ervan in toekomstige experimentele analyses rechtvaardigt. Verder benadrukt de studie dat longitudinale polarisatie effecten worden versterkt in NLO berekeningen, wat suggereert dat polarisatiemetingen als een gevoelige sonde kunnen dienen voor QCD-dynamica in dit energie-regime. De succesvolle reproductie van pQCD signaturen door de PYTHIA 8.316 + POWHEG BOX simulatie valideert het gebruik van NLO+PS tools voor het voorspellen van observables in proton-proton botsingen bij NICA, wat een realistischer beschrijving biedt van exclusieve observables vergeleken met fixed-order berekeningen alleen.