Automated event generation for S-wave quarkonium and leptonium production in NRQCD and NRQED

Deze paper introduceert een uitgebreide versie van het MadGraph5_aMC@NLO-framework die geautomatiseerde berekeningen van S-golf quarkonium- en leptoniumproductie mogelijk maakt binnen NRQCD en NRQED, waarmee robuuste fenomenologische studies in zowel het Standaardmodel als Beyond-the-Standard-Model-scenario's kunnen worden uitgevoerd.

De kern

Stel je voor dat het universum een gigantische, complexe LEGO-set is. De kleinste stukjes waaruit alles is opgebouwd, zijn deeltjes zoals quarks (die protonen en neutronen maken) en elektronen. Soms, onder de juiste omstandigheden, plakken twee van deze deeltjes aan elkaar en vormen ze een stabiel bouwsel. In de natuurkunde noemen we deze gebonden paren quarkonia (als het zware quarks zijn) en leptonia (als het lichte elektronen of muonen zijn).

Deze paper is als het ware een nieuwe, geavanceerde bouwhandleiding voor wetenschappers. Hier is wat ze hebben gedaan, vertaald naar begrijpelijke taal:

1. Het Probleem: De Bouwplaat was te ingewikkeld

Voorheen was het heel lastig om te voorspellen hoe deze deeltjesparen ontstaan in botsende deeltjesversnellers (zoals de LHC in Zwitserland). De wiskunde erachter is extreem complex. Het is alsof je probeert te voorspellen hoe een toren van LEGO-blokken valt, maar je moet rekening houden met duizenden onzichtbare krachten en mogelijke combinaties.

Bestaande software kon dit wel, maar vaak moest je het handmatig doen, stukje bij beetje, of het kon alleen bepaalde soorten deeltjesparen. Het was als een gereedschapskist die alleen hamers had, maar geen schroevendraaiers of tangen.

2. De Oplossing: Een "Alles-in-Één" Bouwrobot

De auteurs van dit paper hebben een nieuwe module toegevoegd aan een beroemde software genaamd MadGraph5_aMC@NLO. Je kunt dit zien als een slimme robot die automatisch de bouwplannen maakt.

• Wat doet hij? Hij kan nu automatisch berekenen hoe vaak en op welke manier deze deeltjesparen (quarkonia en leptonia) worden gemaakt in verschillende soorten botsingen.
• Het grote voordeel: De robot is "slim" genoeg om te weten dat een deeltje niet altijd exact hetzelfde is. Soms is het een strakke, perfecte kluwen (een 'kleur-zinglet' toestand) en soms een wat losser gebonden kluwen (een 'kleur-octet' toestand). De software telt al deze variaties automatisch mee, zonder dat de gebruiker duizenden regels code hoeft te schrijven.

3. De Twee Soorten "Bouwsels"

De robot is getraind op twee soorten bouwsels:

• Quarkonia (De zware broers): Dit zijn paren van zware quarks (zoals charm en bottom). Denk hieraan als zware, stevige LEGO-blokken die aan elkaar plakken. Ze worden vaak gemaakt in de enorme botsingen van de LHC. De wetenschappers hebben getoond dat hun robot nu kan voorspellen hoe vaak deze worden gemaakt, zelfs als er nog andere deeltjes (zoals straling of jets) bij vrijkomen.
• Leptonia (De lichte broers): Dit zijn paren van elektronen, muonen of tau-deeltjes. Denk hieraan als lichte, fragiele bubbels die aan elkaar plakken. Dit is een nieuwe front voor deze software! Voorheen kon niemand dit makkelijk simuleren. De auteurs hebben nu een manier gevonden om ook deze "elektronische atomen" (zoals positronium, een elektron en een positron) te berekenen.

4. Waarom is dit belangrijk? (De Creatieve Analogie)

Stel je voor dat je een detective bent die probeert te begrijpen hoe een gevaarlijke crimineel (een nieuw deeltje of een onbekende kracht) zich verbergt in een drukke stad (de deeltjesversneller).

• Vroeger: Je had een kaart van de stad, maar je moest zelf handmatig elke straat aflopen om te zien of de crimineel daar was. Als je een nieuwe straat wilde checken, moest je de hele kaart opnieuw tekenen.
• Nu: Je hebt een drone (de nieuwe software) die de hele stad in één keer scant. Hij ziet niet alleen de bekende straten, maar ook de steegjes waar de crimineel zich zou kunnen verstoppen.

Bovendien ontdekten de auteurs iets verrassends: Soms zijn de kleine, onbelangrijk lijkende steegjes (subleading bijdragen) juist waar de crimineel zit. Als je alleen kijkt naar de hoofdstraten (de simpele wiskundige regels), mis je het beeld. Hun software laat zien dat je altijd alle mogelijke hoekjes moet checken, anders krijg je een verkeerd antwoord.

5. Wat betekent dit voor de toekomst?

• Voor de LHC: Wetenschappers kunnen nu sneller en nauwkeuriger voorspellen wat ze moeten zien als ze naar zware deeltjes kijken. Dit helpt bij het zoeken naar nieuwe natuurwetten.
• Voor nieuwe deeltjes: Omdat ze nu ook "leptonia" kunnen simuleren, kunnen ze plannen maken voor experimenten om deze zeldzame deeltjesparen daadwerkelijk te vinden.
• Voor iedereen: De software is gratis beschikbaar gemaakt. Het is alsof ze de bouwhandleiding voor de hele wereld hebben gepubliceerd, zodat elke natuurkundige nu deze complexe bouwwerken kan simuleren zonder jarenlang te hoeven studeren.

Kortom: Dit paper introduceert een krachtige, automatische tool die het heelal een stukje begrijpelijker maakt door te laten zien hoe de kleinste bouwstenen van het universum aan elkaar plakken, of het nu zware quarks zijn of lichte elektronen. Het is een enorme stap voorwaarts in het "LEGO-spel" van de deeltjesfysica.

Technische samenvatting

Probleemstelling

De productie van zware gebonden toestanden, zoals quarkonia (bijv. $J/\psi$, $\Upsilon$) en leptonia (bijv. positronium, true muonium), is een cruciaal testveld voor de kwantumchromodynamica (QCD) en kwantumelektrodynamica (QED). Hoewel er bestaande tools zijn voor het berekenen van secties en het genereren van gebeurtenissen, vertonen deze vaak beperkingen:

• Gebrek aan automatisering: Veel bestaande tools (zoals MadOnia of HELAC-Onia) zijn niet volledig geïntegreerd in de moderne MadGraph5_aMC@NLO (MG5_aMC) framework, wat het moeilijk maakt om complexe processen met meerdere gebonden toestanden of processen buiten het Standaardmodel (BSM) te simuleren.
• Beperkte dekking: Er is een gebrek aan geautomatiseerde tools voor leptonium-productie (gebonden toestanden van leptonen), die theoretisch preciezer voorspelbaar zijn dan quarkonia omdat ze vrij zijn van sterke interactie-complexiteiten.
• Complexiteit van NRQCD/NRQED: De berekening vereist een factorisatieformalisme waarbij korte-afstandscross-secties worden gekoppeld aan lange-afstandsmatrix-elementen (LDMEs). Het handmatig combineren van verschillende Fock-toestanden (kleur-singlet en kleur-octet) is foutgevoelig en tijdrovend.

Methodologie

De auteurs hebben een uitgebreide implementatie ontwikkeld binnen het MadGraph5_aMC@NLO framework dat de productie van S-golf quarkonia en leptonia automatiseert. De kern van de methodologie omvat:

1. Theoretisch Kader:

   • Gebruik van NRQCD (voor quarkonia) en NRQED (voor leptonia) factorisatieformaliteiten.
   • De cross-sectie wordt uitgedrukt als een som over Fock-toestanden ($n = ^{2S+1}L_J^{[C]}$), waarbij de korte-afstandscross-sectie ($\hat{\sigma}$) wordt berekend in perturbatieve theorie en de lange-afstandseffecten worden opgesloten in LDMEs.
   • De implementatie ondersteunt momenteel alleen S-golf toestanden ($L=0$). P-golf toestanden zijn voorbehouden voor toekomstige ontwikkelingen.

2. Technische Implementatie (UFO Interface):

   • Er is een nieuwe Universal Feynman Output (UFO) model genaamd sm_onia ontwikkeld.
   • Gebonden toestanden worden behandeld als een nieuwe klasse Boundstate, analoog aan elementaire deeltjes.
   • Fock-toestanden worden gedefinieerd met specifieke notatie (bijv. Jpsi(1|3S11) voor een $J/\psi$ in de $^3S_1^{[1]}$ toestand).
   • Projectie-methode: De auteurs gebruiken covariante projectoren om de amplitude van de constituenten (bijv. $Q\bar{Q}$) te projecteren op de gewenste kwantumtoestanden (kleur, spin, orbitale impulsmoment). Dit gebeurt automatisch tijdens het genereren van Feynman-diagrammen en matrixelementen.
   • LDME Integratie: De lange-afstandsmatrix-elementen worden direct geïntegreerd in de berekening van de gesquareerde amplitudes, in plaats van als externe factoren, wat de flexibiliteit verhoogt.

3. Validatie:

   • De implementatie is gevalideerd tegen HELAC-Onia (een andere toonaangevende tool) voor een breed scala aan processen, variërend van enkelvoudige tot meervoudige quarkonium-productie.
   • Vergelijkingen zijn gemaakt op het niveau van matrixelementen (standalone modus) en totale cross-secties (LO modus), waarbij uitstekende overeenkomst werd gevonden (verschillen binnen floating-point precisie).

Belangrijkste Bijdragen

• Eerste Geautomatiseerde Tool voor Leptonia: Voor het eerst is er een publiek beschikbare, geautomatiseerde tool voor het genereren van gebeurtenissen voor leptonium-productie (positronium, true muonium, true tauonium) binnen een generiek Monte Carlo framework.
• Naadloze Integratie in MG5_aMC: De tool maakt het mogelijk om quarkonium- en leptonium-processen te genereren in combinatie met:
  • Standaardmodel-deeltjes (jets, bosonen).
  • BSM-scenario's en Effectieve Veldtheorieën (via aangepaste UFO-modellen).
  • Deeltjesdouches en hadronisatie (via interface met Pythia).
• Flexibiliteit in Fock-toestanden: Gebruikers kunnen eenvoudig specifieke Fock-toestanden selecteren of fysische toestanden definiëren als een som van bijdragende Fock-toestanden (bijv. $J/\psi = ^3S_1^{[1]} + ^1S_0^{[8]} + ^3S_1^{[8]}$).
• Ondersteuning voor Meerdere Gebonden Toestanden: Het systeem kan processen genereren met meerdere quarkonia of leptonia in de eindtoestand (bijv. drievoudige $J/\psi$ productie), wat essentieel is voor het bestuderen van multi-parton interacties.

Resultaten

De auteurs presenteren een breed scala aan berekende cross-secties voor verschillende botsingsystemen:

• Proton-Proton (LHC, 13 TeV):
  • Inclusieve productie van enkelvoudige charmonia ($\eta_c, J/\psi$) en bottomonia ($\eta_b, \Upsilon$).
  • Geassocieerde productie met elektroweak bosonen ($W, Z, H$) en jets.
  • Productie van $B_c$-mesonen.
  • Meervoudige productie: Cross-secties voor paren ($J/\psi J/\psi$, $\Upsilon \Upsilon$) en zelfs drievoudige $J/\psi$ productie (relevant voor Triple Parton Scattering).
  • Analyse van de impact van deeltjesdouches (Pythia) op de transversale impulsverdelingen ($p_T$).
• Elektron-Proton (EIC/HERA):
  • Berekeningen voor Deep Inelastic Scattering (DIS) en fotoproductie.
  • Onderzoek naar de rol van kleur-octet mechanismen in vergelijking met kleur-singlet kanalen.
• Elektron-Positron (SuperKEKB, Belle/BaBar):
  • Inclusieve en exclusieve productieprocessen (bijv. $e^+e^- \to J/\psi + \gamma$).
  • Vergelijkingen met analytische uitdrukkingen voor leptonia, waarbij de numerieke stabiliteit van de code wordt aangetoond.
• Leptonia:
  • Berekeningen voor positronium, true muonium en true tauonium (ditauonium) in zowel $pp$ als $e^+e^-$ botsingen.
  • Toepassing van een hybride renormalisatieschema om artefacten in de elektroweak koppeling te vermijden.

Significantie

Deze werken vormt een mijlpaal in de theoretische partonfysica:

1. Toegankelijkheid: Het maakt geavanceerde NRQCD/NRQED berekeningen toegankelijk voor een breder publiek door de gebruiksvriendelijke interface van MG5_aMC te gebruiken.
2. Nauwkeurigheid en Betrouwbaarheid: De validatie tegen bestaande tools en analytische formules bevestigt de nauwkeurigheid van de projectiemethoden en de implementatie van LDMEs.
3. Fenomenologische Impact: De tool biedt een krachtig middel om subleading bijdragen te bestuderen. De auteurs benadrukken dat eenvoudige tellingen op basis van koppelingshierarchyën vaak misleidend zijn; kinematische en dynamische factoren (zoals propagator-versterkingen of behoudswetten) kunnen de relatieve gewichten van Fock-toestanden drastisch veranderen.
4. Toekomstperspectief: Het legt de basis voor toekomstige ontwikkelingen, waaronder de automatisering van P-golf toestanden en berekeningen op Next-to-Leading Order (NLO) nauwkeurigheid, wat essentieel is voor precisiestudies in de huidige en toekomstige experimenten (zoals de LHC, EIC en SuperKEKB).

Kortom, dit artikel introduceert een robuust, flexibel en geautomatiseerd instrument dat de studie van gebonden toestanden in hoge-energie fysica significant vooruit helpt, met name door de brug te slaan tussen complexe theoretische formalismen en praktische simulatie-tools.