Neutron tagging of heavy lepton pair production at RHIC and LHC energies

Technische Samenvatting: Neutron Tagging van Zware Leptonenpaarproductie bij RHIC- en LHC-energieën

Probleem en Motivatie
De productie van lepton-antilutonenparen via foton-foton fusie ($\gamma\gamma \to \ell^+\ell^-$) in ultraperifere botsingen (UPC's) van zware ionen dient als een strenge test voor Kwantumelektrodynamica (QED) in sterke elektromagnetische velden. Hoewel elektron-positron paarproductie uitgebreid is bestudeerd, zijn de theoretische beschrijvingen van zwaardere leptonenparen (muonen en tau's) in combinatie met specifieke neutronenemissietopologieën minder ontwikkeld in de literatuur. In UPC's, waarbij kernen elkaar passeren op afstanden van de impactparameter die groter zijn dan de som van hun radii, fungeren de intense Lorentz-gebuste Coulombvelden als bronnen van virtuele fotonen. Deze fotonen kunnen wederzijdse elektromagnetische excitatie van de kernen induceren, wat leidt tot de emissie van neutronen via de Reusachtige Dipoolresonantie (GDR) en andere mechanismen. Experimentele samenwerkingen (STAR, ALICE, ATLAS, CMS) categoriseren gebeurtenissen op basis van neutrondetectie in Zero Degree Calorimeters (ZDC's) in klassen zoals $0n0n$ (geen neutronen), $0nXn$ (neutronen in één ZDC) en $XnXn$ (neutronen in beide ZDC's). Een consistent theoretisch kader is vereist om cross-secties voor deze specifieke neutron-getagde klassen te berekenen om directe vergelijking met experimentele data over verschillende centrum-van-massa-energieën ($\sqrt{s_{NN}} = 200$ GeV, 2,76 TeV en 5,02 TeV) mogelijk te maken.

Methodologie
De auteurs voeren berekeningen uit voor de laagste-orde QED (LO QED) voor $\gamma\gamma \to \ell^+\ell^-$ ($\ell = e, \mu, \tau$) productie in Au+Au en Pb+Pb botsingen. De methodologie bestaat uit drie primaire componenten:

1. Vrije Paarproductiekans: Het proces wordt behandeld door de ionentrajecten te beschouwen als klassieke externe elektromagnetische bronnen. De cross-sectie wordt afgeleid van S-matrix elementen die betrokken zijn bij directe en gekruiste Feynman-diagrammen. Om de meerdimensionale integralen met betrekking tot Bessel-functies, die snelle oscillaties vertonen, te behandelen, gebruiken de auteurs een Monte Carlo-techniek. Ze berekenen een functie $F(q)$ (gerelateerd aan de momentumoverdracht) en passen een enkele exponentiële vorm toe, $F(q) = F(0)e^{-aq}$, met behulp van niet-lineaire kleinste-kwadraten-fitting. Dit maakt de analytische integratie van de impactparameter ($b$) afhankelijkheid mogelijk om de kans op paarproductie, $P_{\ell^+\ell^-}(b)$, te bepalen.

2. Neutronemissie en Nucleaire Excitatie: De kans op nucleaire excitatie wordt gemodelleerd door de fotonflux te integreren met de totale fotoabsorptie cross-sectie ($\sigma_{tot}$). De fotoabsorptie cross-sectie wordt geconstrueerd door een Lorentziaanse fit voor het GDR-gebied (met recente experimentele data voor $^{197}$Au en $^{208}$Pb) te combineren met het Kossov-model voor niet-GDR bijdragen (quasi-deuteron, $\Delta$-resonantie, Roper-resonantie en Regge-Pomeron). Het gemiddelde aantal geabsorbeerde fotonen wordt gebruikt binnen een Poisson-statistiek kader om de kansen te bepalen voor het emitteren van nul ($0n$), één ($1n$) of ten minste één ($Xn$) neutron. Dit maakt de categorisering van gebeurtenissen in $0n0n$, $0n1n$, $1n1n$, $0nXn$ en $XnXn$ klassen mogelijk.

3. Geen Hadronische Interactiekans: Om te garanderen dat de botsingen ultraperifeer blijven, wordt de kans op geen hadronische interactie, $P_{nohad}(b)$, berekend met behulp van het Glauber-model met Woods-Saxon nucleaire dichtheidsverdelingen. Deze kans wordt vermenigvuldigd met de paarproductie- en excitatiekansen om de uiteindelijke getagde cross-secties te verkrijgen.

De berekeningen worden uitgevoerd binnen de fiduciale fase-ruimte restricties gedefinieerd door de STAR, ALICE, ATLAS en CMS experimenten, inclusief cuts op pseudorapiditeit ($|\eta|$), transversale momentum ($p_T$), invariante massa ($M_{\ell\ell}$) en paar-transversale momentum ($p_{T,\ell\ell}$).

Belangrijkste Resultaten
De studie presenteert differentiële en geïntegreerde cross-secties voor elektron-, muon- en tau-paarproductie, uitgesplitst naar neutronemissieklassen:

• Elektronenparen (STAR, $\sqrt{s_{NN}}=200$ GeV): De berekende inclusieve cross-sectie is 5850 $\mu$b. De $XnXn$ getagde cross-sectie wordt berekend als 274,5 $\mu$b, wat ongeveer 4-5% hoger is dan de STAR-meting (261 $\mu$b) en 24,5% hoger dan de STARLight-voorspelling. De differentiële distributies voor invariante massa en transversale momentum vertonen een goede overeenstemming met de STAR-data, hoewel het model de hoog-$p_T$ regio licht onderschat.
• Elektronenparen (ALICE, $\sqrt{s_{NN}}=2,76$ TeV): Voor lage en hoge massaregio's liggen de berekende inclusieve cross-secties (147 $\mu$b en 87,7 $\mu$b) 3-5% onder de ALICE-metingen maar 13-15% boven de STARLight-voorspellingen. Het $0n0n$ kanaal domineert de cross-sectie.
• Elektron- en Muonenparen (ATLAS, $\sqrt{s_{NN}}=5,02$ TeV):
  • Voor elektronen is de inclusieve cross-sectie 233 $\mu$b, wat 8,3% boven de ATLAS-meting (215 $\mu$b) ligt.
  • Voor muonen is de inclusieve cross-sectie 36,5 $\mu$b, wat 7,0% boven de ATLAS-meting (34,1 $\mu$b) ligt.
  • In beide gevallen levert het $0n0n$ kanaal de dominante bijdrage. De differentiële distributies in invariante massa en rapiditeit komen over het algemeen goed overeen met de ATLAS-data, met ratio's nabij eenheid, hoewel de afwijkingen toenemen bij hogere invariante massa's.
• Tau-paren (CMS, $\sqrt{s_{NN}}=5,02$ TeV): De auteurs berekenen een inclusieve cross-sectie van 1043 $\mu$b voor $\tau^+\tau^-$ productie. Dit resultaat komt nauw overeen met de Shao en d'Enterria LO-berekening met de nucleaire ladingvormfactor (1060 $\mu$b), maar is aanzienlijk hoger dan de geëxtrapoleerde CMS fiduciale resultaten (580–850 $\mu$b). Het $0n0n$ kanaal draagt 802 $\mu$b bij.

Betekenis en Claims
Het artikel beweert een consistent theoretisch kader te bieden voor het berekenen van neutron-getagde leptonenpaarproductie over een breed scala aan leptonmassa's en botsingsenergieën. Een primaire bijdrage is de uitbreiding van de neutron-tagging analyse van elektron-positron paren naar zwaardere muon- en tau-paren, die minder aandacht hebben gekregen in de literatuur. De auteurs demonstreren dat hun LO QED-benadering, gecombineerd met een Poisson-gebaseerd neutronemissiemodel en een Glauber-gebaseerde hadronische exclusie, resultaten oplevert die over het algemeen consistent zijn met experimentele data van STAR, ALICE, ATLAS en CMS binnen de verwachte theoretische onzekerheden.

De auteurs merken bescheiden op dat zij geen hogere-orde QED-correcties (zoals final state radiation) opnemen, hoewel zij literatuur citeren die suggereert dat deze effecten minimaal zijn voor totale cross-secties. Zij erkennen ook dat hun behandeling van neutronemissiefracties als constant over kinematische bins een vereenvoudiging is, wat de lichte discrepanties in de vorm van differentiële distributies bij hogere invariante massa's kan verklaren. Desalniettemin valideert de studie het nut van deze vereenvoudigde benadering voor het beschrijven van de bulk van de leptonpaarproductie in UPC's en biedt het een baseline voor toekomstige precisie-studies.