NNLOJET: a parton-level event generator for jet cross sections at NNLO QCD accuracy

Dit artikel introduceert de open-source NNLOJET-deeltjesniveau-evenementengenerator, die de antenne-subtractiemethode implementeert om jetsnede en gerelateerde observabelen te berekenen met NNLO QCD-nauwkeurigheid voor elektron-positron-, lepton-hadron- en hadron-hadronbotsingen.

De kern

Het Grote Geheel: Een Beter Kristallen Bol Bouwen

Stel je voor dat je een natuurkundige bent die probeert te voorspellen wat er gebeurt wanneer twee tinyeltjes deeltjes met bijna de lichtsnelheid op elkaar botsen. Je hebt een set regels (het Standaardmodel) die je vertelt hoe deze deeltjes zouden moeten gedragen. Maar de natuur is rommelig. Wanneer deeltjes botsen, stuiteren ze niet zomaar af; ze spuiten een chaotische stortvloed van nieuwe deeltjes uit, zoals confetti dat uit een kanon explodeert.

Om een werkelijk nauwkeurige voorspelling te doen, kun je niet alleen gokken over de hoofdexplosie. Je moet de kleine, bijna onzichtbare rimpelingen en secundaire spreuken berekenen die gebeuren door de explosie. In de wereld van de deeltjesfysica heten deze berekeningen NNLO (Next-to-Next-to-Leading Order). Ze zijn ongelooflijk complex, zoals proberen het exacte pad van elk enkel zandkorreltje in een orkaan te voorspellen.

Dit paper introduceert NNLOJET, een nieuw, open-source computerprogramma dat is ontworpen om deze super-nauwkeurige berekeningen te doen voor "jets" (spreuken van deeltjes). Het is alsof je wetenschappers een kristallen bol in hoge definitie geeft die details kan zien die de oude modellen misten.

Het Kernprobleem: De "Oneindige" Rommel

Wanneer natuurkundigen proberen deze botsingen te berekenen, lopen ze tegen een mathematische nachtmerrie aan die "infrarood singulariteiten" wordt genoemd.

• De Analogie: Stel je voor dat je probeert het aantal mensen in een kamer te tellen, maar elke keer als iemand fluistert, wordt de kamer oneindig luider. Als je probeert het lawaai op te tellen, gaat het getal naar oneindig.
• De Realiteit: In de deeltjesfysica, wanneer deeltjes zeer langzaam bewegen of zeer dicht bij elkaar komen, loopt de wiskunde uit tot oneindig. Aangezien de echte wereld geen oneindige energie heeft, zijn deze oneindigheden gewoon fouten in de rekenmethode.

Om dit op te lossen, gebruiken de auteurs een methode die Antenna Subtraction (Antenne-subtractie) wordt genoemd.

• De Analogie: Denk eraan als geluidsdempende koptelefoons. De koptelefoons luisteren naar het vervelende achtergrondlawaai (de wiskundige oneindigheden) en genereren een "negatieve" geluidsgolf om het perfect op te heffen.
• Hoe het werkt: Het programma berekent de rommelige, oneindige delen apart en trekt ze vervolgens af, waardoor alleen het schone, eindige, reële antwoord overblijft.

Wat is NNLOJET?

NNLOJET is de "motor" die deze berekeningen draait. Voor dit paper waren veel van deze berekeningen als geheime recepten opgesloten in de keuken van één chef-kok. Als je ze wilde gebruiken, moest je de chef smeken of je eigen keuken van scratch bouwen.

NNLOJET verandert het spel door open-source te zijn.

• De Analogie: Het is alsof je het volledige receptenboek en de keukengerei aan het publiek vrijgeeft. Iedereen met een computer kan het downloaden, het gerecht bereiden en de ingrediënten controleren.
• Het "Jet"-Deel: Het programma is gespecialiseerd in "jets". In de deeltjesfysica, wanneer een quark of gluon (de bouwstenen van materie) wordt uitgeschoten, blijft hij niet alleen. Hij verandert direct in een sprei van deeltjes. We noemen deze sprei een "jet". NNLOJET voorspelt precies hoe groot, hoe snel en in welke richting deze spreuken zullen gaan.

Hoe Gebruik Je Het? (De Runcard)

Je hoeft geen code-wizard te zijn om NNLOJET te gebruiken. Het paper legt uit dat je het programma bedient met een Runcard.

• De Analogie: Denk aan de Runcard als een vluchtplan of een receptkaart. Je hoeft niet te weten hoe de vliegtuigmotor werkt; je moet de piloot alleen vertellen waar hij naartoe moet.
• Wat je schrijft: Je vertelt het programma:
  • Waar: Botsen we protonen (zoals bij de Large Hadron Collider) of elektronen?
  • Wat: Zoeken we naar een Z-boson, een Higgs-boson, of gewoon een sprei van jets?
  • Regels: Hoe groot moeten de jets zijn? Hoe snel?
  • Output: Welke grafieken wil je aan het einde zien?

Het paper biedt een gedetailleerde handleiding (secties 5 en 6) over hoe je dit "vluchtplan" schrijft, waarbij alles wordt behandeld, van het type botsing tot de specifieke wiskundige instellingen.

De Werkstroom: De Assemblagelijn

Het berekenen van deze botsingen is zo zwaar dat één computer het niet alleen kan doen. Het zou jaren duren.

• De Analogie: Stel je voor dat je een enorm muurschildering moet schilderen. Eén persoon zou een leven nodig hebben. In plaats daarvan huur je een team van 100 schilders in.
• Het Hulpmiddel: Het paper beschrijft een werkstroomschript (nnlojet-run) dat fungeert als de voorman.
  1. Opwarmen: De voorman stuurt een paar schilders uit om de muur te testen en de beste manier van schilderen uit te zoeken (dit is de aanpassing van het "VEGAS-rooster").
  2. Productie: Zodra het plan staat, stuurt de voorman honderden schilders uit om het echte werk gelijktijdig te doen.
  3. Afronden: Als iedereen klaar is, verzamelt de voorman alle stukken, controleert op fouten (uitbijters) en plakt ze samen tot één perfect plaatje.

Wat Kan Het Berekenen?

Het paper somt de specifieke "recepten" op die momenteel beschikbaar zijn in het kookboek (Versie 1.0.0):

• Elektron-Positron Botsingen: Zoals bij de oude LEP-collider.
• Proton-Proton Botsingen: Zoals bij de huidige Large Hadron Collider (LHC).
• Specifieke Gebeurtenissen: Het kan de creatie voorspellen van:
  • Z-bosonen (zware neven van het foton).
  • W-bosonen (deeltjes die de zwakke kracht dragen).
  • Higgs-bosonen (het deeltje dat massa geeft aan anderen).
  • Fotonen (licht) en jets.
  • Combinaties van al het bovenstaande (bijvoorbeeld een Higgs-boson plus een jet).

Waarom Is Dit Belangrijk?

Het paper concludeert dat door deze code open en gebruiksvriendelijk te maken, het wetenschappers in staat stelt hun theoretische voorspellingen direct te vergelijken met echte experimentele data.

• De Analogie: Als je een modelauto bouwt, moet je weten of hij daadwerkelijk zal rijden. NNLOJET biedt de meest accurate "windtunneltest" die mogelijk is. Als de echte data overeenkomt met de NNLOJET-voorspelling, is ons begrip van het universum correct. Als ze niet overeenkomen, kan het betekenen dat we iets nieuws hebben gevonden!

Kortom: Dit paper release een krachtig, gratis en gebruiksvriendelijk hulpmiddel dat natuurkundigen helpt de rommelige details van deeltjesbotsingen te berekenen met extreme precisie, zodat onze kaart van het universum zo nauwkeurig mogelijk is.

Technische samenvatting

Technische Samenvatting van NNLOJET: Een Deeltjesniveau-Evenementengenerator voor Jet-dwarsdoorsneden met NNLO QCD-nauwkeurigheid

Probleemstelling
Precisietests van het Standaardmodel (SM) zijn sterk afhankelijk van het vergelijken van experimentele metingen van hadronische eindtoestanden (jets) met theoretische voorspellingen. Deze voorspellingen vereisen de berekening van Next-to-Next-to-Leading Order (NNLO) correcties in de Quantum Chromodynamica (QCD) om de grote sterke koppelingsconstante ($\alpha_s$) en de complexe structuur van infrarood (IR) singulariteiten die voortvloeien uit zachte en collineaire straling te hanteren. Hoewel Next-to-Leading Order (NLO) berekeningen grotendeels geautomatiseerd zijn, blijven NNLO-berekeningen uitdagend vanwege de noodzaak om IR-divergenties te cancelen tussen dubbel-virtuele (VV), gemengd real-virtuele (RV) en dubbel-reële (RR) bijdragen. Bestaande NNLO-codes zijn vaak proces-specifiek, niet breed beschikbaar als open-source, of vereisen aanzienlijke handmatige afstelling voor elke observabele. Er is behoefte aan een flexibele, open-source deeltjesniveau-Evenementengenerator die in staat is om jet-dwarsdoorsneden met NNLO-nauwkeurigheid te berekenen in diverse collideromgevingen (elektron–positron, lepton–hadron en hadron–hadron).

Methodologie: Antenne-subtractie
De NNLOJET-code implementeert de antenne-subtractiemethode om IR-singulariteiten te hanteren. Deze methode leidt subtractietermen af op basis van de universele factorisatie van QCD-reële-stralingsamplitudes in hun IR-grenzen.

• Kernmechanisme: De methode identificeert onopgeloste reële-stralingsconfiguraties met paren van harde radiator-deeltjes. De straling tussen hen wordt beschreven door "antennefuncties".
• Subtractietermen: Niet-geïntegreerde subtractietermen worden geconstrueerd uit producten van antennefuncties en gereduceerde matrixelementen van lagere deeltjesmultipliciteit. Deze termen worden geïntegreerd over de antenne-fasruimte om geïntegreerde subtractietermen ($d\hat{\sigma}_T$ en $d\hat{\sigma}_U$) te produceren die expliciet de $\epsilon$-polen van virtuele correcties cancelen.
• Implementatie: De code maakt gebruik van een fasruimte-mapping om de kinematica van de antennefuncties te factoriseren van de gereduceerde matrixelementen. De implementatie beslaat eind–eind, begin–eind en begin–begin radiatorconfiguraties.
• Berekeningskader: NNLOJET is een deeltjesniveau-evenementengenerator die voornamelijk is geschreven in moderne Fortran, met C++ en Python-modules. Het maakt gebruik van de adaptieve Monte Carlo-routine VEGAS voor fasruimte-integratie. De code bevat bibliotheken voor het evalueren van speciale functies (harmonische polylogaritmen) via HPLOG, TDHPL en CHAPLIN, en koppelt aan LHAPDF voor parton-distributiefuncties (PDF's).

Belangrijkste Bijdragen en Kenmerken
Het artikel presenteert de open-source release (v1.0.0) van de NNLOJET-code, verspreid onder de GNU General Public License v3.0. Belangrijke technische bijdragen omvatten:

1. Brede Procesdekking: De code ondersteunt NNLO-berekeningen voor jet-productie in:

   • Elektron–positron-botsingen: $e^+e^- \to 2j$ en $3j$.
   • Lepton–hadron-botsingen: Neutral-current (NC) en charged-current (CC) DIS-processen ($ep \to \ell + 2j$, $ep \to \nu + 2j$, enz.).
   • Hadron–hadron-botsingen: Inclusieve jet-productie ($pp \to j$), dijet-productie ($pp \to 2j$), vectorboson plus jet ($V+1j$), foton plus jet ($\gamma+1j$), di-foton ($\gamma\gamma$), en Higgs-boson plus jet ($H+1j$) productie.
   • Afgeleide Ordes: Het kader biedt impliciet NLO-correcties voor processen met één extra jet en LO voor twee extra jets (bijvoorbeeld het berekenen van $Z+2j$ op NLO via het $ZJ$-proces met een 2-jet vereiste).

2. Flexibele Gebruikersinterface (Runcard): De code wordt aangestuurd via een plain-ASCII "runcard"-bestand, waarmee gebruikers kunnen definiëren:

   • Proces en Collider: Selectie van collidertype (pp, ep, $e^+e^-$), bundelenergieën en jet-algoritmen (anti-$k_T$, Cambridge-Aachen, $k_T$, Durham, Jade).
   • Observabelen en Selecties: Definitie van fiduciële cuts, object-niveau selectoren (jet $p_T$, rapiditeit) en event-niveau selectoren.
   • Histogrammen: Flexibele binning (lineair, logaritmisch, aangepaste randen) en multi-differentiële verdelingen via histogram-selectoren.
   • Schalen: Simultane evaluatie van meerdere combinaties van renormalisatie- en factorisatieschalen voor onzekerheidsschatting.

3. Geautomatiseerde Werkstroom: De distributie bevat een Python-gebaseerde werkstroomschript (nnlojet-run) om taakindiening, resource-toewijzing en resultaat-merging te beheren. Deze werkstroom behandelt:

   • Warm-up en Productiefasen: Adaptieve grid-generatie gevolgd door productieruns.
   • Resourcebeheer: Dynamische aanpassing van rekenresources om een doelwit relatieve nauwkeurigheid te bereiken.
   • Resultaat-Merging: Een robuust algoritme (involvend trimming, k-scan en gewogen combinatie) om gedeeltelijke resultaten van onafhankelijke taken te samenvoegen terwijl de impact van uitbijters wordt geminimaliseerd.

4. Validatie: De implementatie bevat rigoureuze validatiecontroles:

   • Punttests: Verificatie dat de verhouding van matrixelementen tot subtractietermen de eenheid nadert nabij singuliere grenzen.
   • Technische cut-controles: Waarborgen van de onafhankelijkheid van resultaten van de technische infrarood-cutoff.
   • Pool-tests: Expliciete cancelatie van $\epsilon$-polen in RV- en VV-bijdragen.

Resultaten en Validatie
Het artikel beschrijft de implementatie van diverse processen, met verwijzing naar specifieke publicaties voor validatiedetails. Als concreet voorbeeld presenteren de auteurs een berekening voor $Z + 1j$-productie bij de LHC ($\sqrt{s} = 7$ TeV):

• Opstelling: Anti-$k_T$ jet-algoritme ($R=0.4$), NNPDF3.1 NNLO PDF's, en specifieke fiduciële cuts op leptonen en jets.
• Resultaat: De totale NNLO fiduciële dwarsdoorsnede wordt berekend op 68,5(2) pb.
• Differentiële Verdelingen: Het artikel toont de transversale impulsverdeling van de leidende jet, waarbij de convergentie van LO naar NLO naar NNLO wordt getoond, met onzekerheidsbanden die een 7-punts schaalvariatie vertegenwoordigen.

Betekenis en Claims
De auteurs claimen dat NNLOJET een universeel toepasbaar, open-source kader biedt voor NNLO QCD-berekeningen van jet-observabelen, waarmee het gat wordt gedicht dat door eerdere codes is achtergelaten, die vaak propriëtair waren of beperkt tot specifieke procesklassen.

• Toegankelijkheid: Door de code te distribueren onder GPLv3.0 en een gebruikersvriendelijke runcard-syntaxis en werkstroomschript te bieden, streven de auteurs ernaar om hoog-precisie NNLO-voorspellingen toegankelijk te maken voor een bredere gemeenschap van fenomenologen en experimentalisten.
• Flexibiliteit: Het vermogen om willekeurige fiduciële cuts en observabelen te definiëren zonder de code opnieuw te compileren, maakt directe vergelijking met experimentele data mogelijk.
• Uitbreidbaarheid: Het modulaire ontwerp staat toe om nieuwe processen en functies toe te voegen in toekomstige releases.

Het artikel concludeert dat NNLOJET een volwassen hulpmiddel is voor precisiefysica-studies, het testen van de grenzen van het Standaardmodel, en het verfijnen van de bepaling van SM-parameters en parton-distributiefuncties. De code is beschikbaar op HEPForge, en de auteurs nodigen gebruikersfeedback en bugrapporten uit.