$gg \to ZH$ at NLO matched to parton showers with ggxy and POWHEG

Dit artikel beschrijft de implementatie van de NLO QCD-correcties voor het proces $gg \to ZH$ in ggxy, inclusief een interface naar POWHEG voor de koppeling met parton showers in Pythia, wat flexibele simulaties met verschillende top-quark-massaschema's en volledige $Z$-bosonvervallingen mogelijk maakt.

De kern

De Zelfrijdende Auto en de Grote Ontdekkingsreis: Een Uitleg van het Wetenschappelijke Papier

Stel je voor dat deeltjesversnellers, zoals de Large Hadron Collider (LHC) in Zwitserland, enorme racecircuits zijn waar atomen met bijna de lichtsnelheid tegen elkaar worden gebotst. Wetenschappers hopen hierbij nieuwe deeltjes te vinden of bestaande te begrijpen, zoals de Higgs-boson (het "Goddeeltje") en de Z-boson.

Dit wetenschappelijke papier gaat over een heel specifieke botsing: twee gluonen (de lijm die atoomkernen bij elkaar houdt) die botsen en een Higgs-boson en een Z-boson produceren. In het jargon heet dit proces gg → ZH.

Hier is wat de auteurs hebben gedaan, vertaald naar alledaagse taal:

1. Het Bouwplan (De Theorie)

Stel je voor dat je een heel complexe machine wilt bouwen. Je hebt een ontwerp nodig. In de natuurkunde zijn dit de wiskundige formules die voorspellen hoe vaak zo'n botsing gebeurt en welke deeltjes eruit komen.

De auteurs hebben een nieuwe, zeer geavanceerde rekenmachine (een computerprogramma genaamd ggxy) gebouwd.

• Het probleem: De oude formules waren te complex of niet beschikbaar voor deze specifieke botsing.
• De oplossing: Ze hebben de wiskundige "blauwdrukken" (die eerder door andere wetenschappers zijn berekend) in hun nieuwe software ingebouwd.
• De analogie: Het is alsof ze een nieuwe, super-snelle GPS hebben gemaakt voor een route die voorheen alleen op papier bestond. Nu kunnen ze de reis in detail plannen.

2. De Twee Manieren van Rekenen (Stabiel vs. Instabiel)

In de natuurkunde zijn deeltjes soms "stabiel" (ze blijven bestaan zoals ze zijn) en soms "instabiel" (ze vervallen direct in andere deeltjes, zoals een vuurwerk dat ontploft).

• Stabiele Z-boson: De auteurs laten zien hoe het programma werkt als de Z-boson gewoon blijft bestaan. Dit is de basisversie.
• Instabiele Z-boson (met vervalling): In werkelijkheid vervalt de Z-boson vaak direct in een paar elektronen of neutrino's. De auteurs hebben hun programma zo aangepast dat het deze "ontploffing" en de draaiing (spin) van de deeltjes meerekent.
• De analogie: Stel je voor dat je een bal gooit.
  • Versie 1: Je rekent uit waar de bal landt als hij heel blijft.
  • Versie 2: Je rekent uit wat er gebeurt als de bal in de lucht ontploft in twee stukjes die allebei een eigen baan volgen. Dit is veel moeilijker te berekenen, maar geeft een realistischer beeld.

3. De "Twee Sporen" (De Topquark)

Een belangrijk detail in de berekening is de massa van het zwaarste deeltje, de topquark. In de wiskunde kun je deze massa op twee manieren behandelen: als een exacte waarde (de "pool" methode) of als een gemiddelde waarde (de "MS" methode).

• De auteurs hebben hun software zo gemaakt dat je kunt kiezen welke methode je wilt gebruiken, net zoals je kunt kiezen of je een auto in kilometers of mijlen wilt laten rijden. Ze ontdekten dat deze keuze het eindresultaat flink kan beïnvloeden, vooral bij hoge energieën.

4. De Realistische Simulatie (POWHEG en Pythia)

Dit is misschien wel het coolste deel. Tot nu toe hebben we het gehad over de "theoretische botsing" op het moment van impact. Maar in een echte detector zie je niet alleen dat ene moment; je ziet een hele cascade van deeltjes die eruit spatten.

• POWHEG: Dit is de "meesterchef" die de perfecte botsing berekent.
• Pythia: Dit is de "assistent" die de rest van het feest regelt. Hij zorgt voor de "parton showers" (de regen van deeltjes) die na de botsing ontstaan.
• De Analogie:
  • ggxy is de architect die het perfecte ontwerp van het huis maakt.
  • POWHEG is de aannemer die het huis bouwt volgens dat ontwerp.
  • Pythia is de interieurontwerper die meubels, gordijnen en verlichting toevoegt om het huis bewoonbaar en realistisch te maken.

De auteurs hebben hun nieuwe architectuur (ggxy) gekoppeld aan de aannemer (POWHEG) en de interieurontwerper (Pythia). Hierdoor kunnen ze nu niet alleen zeggen "dit huis wordt gebouwd", maar ook "dit is hoe het eruitziet als je erin loopt, inclusief alle kleine details".

5. Waarom is dit belangrijk?

De LHC wordt steeds preciezer. Om de resultaten van de machines te begrijpen, moeten de theorieën net zo precies zijn.

• Vroeger was deze specifieke berekening (gluon-gluon botsing naar Higgs en Z) niet beschikbaar in publieke software.
• Nu hebben de auteurs een gratis, openbaar gereedschap gemaakt (de code is online te downloaden).
• Dit stelt andere wetenschappers in staat om hun eigen analyses te doen, net zoals iemand die een nieuwe auto wil bouwen nu gebruik kan maken van een open-source bouwplan in plaats van alles zelf opnieuw uit te vinden.

Samenvattend:
De auteurs hebben een nieuwe, flexibele rekenmachine gebouwd die het heel precies kan voorspellen wat er gebeurt als twee gluonen botsen en een Higgs- en Z-deeltje maken. Ze hebben dit zo gemaakt dat het rekening houdt met de "ontploffing" van deeltjes en dat het gekoppeld kan worden aan simulaties van de rest van het botsingsproces. Hierdoor kunnen natuurkundigen de data van de LHC veel beter begrijpen en nieuwe ontdekkingen doen.

Technische samenvatting

Titel en Context

Het paper beschrijft de implementatie van de recent berekende analytische uitdrukkingen voor de Next-to-Leading Order (NLO) QCD-correcties voor het proces $gg \to ZH$ (gluon-fusie naar een Z-boson en een Higgs-boson) in de C++-bibliotheek ggxy. Daarnaast wordt een interface naar POWHEG ontwikkeld om deze berekeningen te koppelen aan parton-shower-simulaties (met Pythia).

Het Probleem

De associatieve productie van een Higgs-boson met een gauge-boson ($ZH$) is cruciaal voor het bestuderen van de eigenschappen van het Higgs-boson. Hoewel de productie via het Drell-Yan-type proces ($q\bar{q} \to ZH$) goed begrepen is, is het gluon-geïnduceerde kanaal ($gg \to ZH$) numeriek significant, vooral bij hoge energieën.

• Aanwezige tekortkomingen: Tot nu toe waren de volledige NLO QCD-correcties voor $gg \to ZH$ niet beschikbaar in publiek toegankelijke computercodes. Bestaande tools zoals vh@nnlo en HAWK behandelen voornamelijk het Drell-Yan-deel of gebruiken benaderingen.
• Complexiteit: De berekening vereist tweeringsamplitudes (two-loop amplitudes) en de behandeling van complexe fenomenen zoals off-shell effecten van het Z-boson, spin-correlaties bij verval, en de keuze van het renormalisatieschema voor de top-quark massa.

Methodologie

1. Implementatie in ggxy:
De auteurs hebben de analytische uitdrukkingen geïmplementeerd die zijn verkregen via twee complementaire expansiemethoden (zoals beschreven in Ref. [14]):

• Forward-expansie: Geldig wanneer $|t|, q_Z^2, q_H^2 \ll m_t^2, s$.
• High-energy-expansie: Geldig wanneer $q_Z^2, q_H^2 \ll m_t^2 \ll s, |t|$.
• Schakelpunten: Er worden Padé-benaderingen gebruikt en dynamische schakelpunten (tussen 200 en 260 GeV voor $p_T$) om de nauwkeurigheid over de volledige faseruimte te garanderen. Voor lage energieën ($\sqrt{s} < 200$ GeV) wordt overgeschakeld op de expansie in de grote top-quark massa.
• Formulering: De amplitude wordt uitgedrukt in zes onafhankelijke formfactoren ($F_{12}^\pm, F_{2}^-, F_{3}^\pm, F_4$). De bibliotheek biedt functies voor het berekenen van deze formfactoren en de virtuele correcties ($V_{fin}$).

2. Behandeling van Real Correcties:

• De real correcties worden berekend met Recola (in combinatie met Collier en CutTools) voor de 2→3 matrixelementen.
• Een belangrijke nuance is het filteren van "Higgs-Strahlung" diagrammen (Drell-Yan type). De implementatie in ggxy focust uitsluitend op bijdragen met één gesloten quark-lus, waarbij het Z-boson en/of Higgs-boson aan deze lus gekoppeld zijn. Diagrammen waar het Z-boson van een externe quarklijn wordt uitgestraald, worden uitgesloten om duplicatie met het Drell-Yan-proces te voorkomen.

3. Off-shell en Verval:

• De code ondersteunt zowel stabiele Z-bosons als het verval in leptonen ($Z \to \ell^-\ell^+$) of neutrino's ($Z \to \nu\bar{\nu}$).
• Off-shell effecten en spin-correlaties worden meegenomen door de helicity-amplitudes van $gg \to Z^*H$ te combineren met die van $Z^* \to \ell\ell$.
• De maximale "off-shellness" van het Z-boson is beperkt tot ongeveer 150 GeV vanwege de gebruikte expansies.

4. Interface naar POWHEG:

• Een nieuwe module ggxy_ggZH is ontwikkeld voor POWHEG.
• Dit koppelt de NLO-berekeningen van ggxy (virtuele correcties) en Recola (real correcties) aan parton showers (Pythia 8).
• De implementatie maakt het mogelijk om de Z-vervalprocessen ofwel op het niveau van de matrixelementen (met spin-correlaties) ofwel via de parton shower te genereren.

Belangrijkste Resultaten

1. Validatie en NLO Cross Sections:

• De resultaten van ggxy voor de totale en differentieel cross sections werden gevalideerd tegen bestaande literatuur (Ref. [13, 34]) en interne Mathematica-codes.
• Voor stabiele Z-bosons bij $\sqrt{s} = 13.6$ TeV werd een NLO cross section van 130.50 fb gevonden, wat uitstekend overeenkomt met referentiewaarden.
• De NLO correcties zijn groot (ongeveer +85% ten opzichte van LO).

2. Invloed van Renormalisatieschema:

• Er is een vergelijking gemaakt tussen het On-Shell (OS) en $\overline{MS}$ schema voor de top-quark massa.
• Het kiezen van het $\overline{MS}$ schema leidt tot een neerwaartse verschuiving van de cross section bij hoge invariantmassa's (10-40%), wat een onzekerheid introduceert die vergelijkbaar is met de schaalvariatie.

3. Combinatie met Drell-Yan (vh@nnlo):

• De gluon-fusie bijdrage ($gg \to ZH$) kan veilig worden gecombineerd met de Drell-Yan bijdragen uit vh@nnlo.
• Bij LHC-energieën ($\sqrt{s} \approx 13-14$ TeV) draagt het $gg$-kanaal ongeveer 15% bij aan de totale $ZH$-productie. Bij hogere energieën (tot 100 TeV) neemt dit aandeel toe tot wel 25%.

4. Parton Shower Effecten (Pythia):

• Spin-correlaties: Het is cruciaal om spin-correlaties op het niveau van de matrixelementen te behouden. Simulaties zonder spin-correlaties (waarbij de Z als stabiel wordt behandeld en later verval) leiden tot fouten van 40-100% in hoekverdelingen (zoals de azimutale hoek tussen leptonen).
• hdamp-parameter: De parameter $h_{damp}$ in POWHEG wordt gebruikt om de real correcties te splitsen. Een waarde van $h_{damp} = 250$ (in plaats van $\infty$) zorgt ervoor dat de parton shower-resultaten beter overeenkomen met de vaste-orde resultaten in de staart van de verdeling (hoge $p_T$), waardoor de onzekerheid in de staart van de $p_T$-verdeling van het Higgs-boson wordt verminderd.

Significantie en Conclusie

Dit werk vult een belangrijke lacune in de theoretische voorspellingen voor Higgs-fysica op.

1. Beschikbaarheid: Het biedt voor het eerst een publiek beschikbare, flexibele tool (ggxy) voor nauwkeurige NLO-berekeningen van $gg \to ZH$, inclusief de mogelijkheid om verschillende top-quark massaschema's en Z-ontvallen te testen.
2. Precisie: Door de implementatie in POWHEG kunnen experimentele analyses nu rekening houden met parton-shower-effecten en spin-correlaties op een consistente manier, wat essentieel is voor de interpretatie van data van de Large Hadron Collider (LHC) en toekomstige colliders.
3. Fenomenologie: De studie benadrukt dat het negeren van spin-correlaties en off-shell effecten leidt tot significante fouten in specifieke observabelen, wat onderstreept dat deze geavanceerde berekeningen noodzakelijk zijn voor toekomstige precisie-metingen.

De broncode en voorbeeldbestanden zijn beschikbaar via GitLab, wat onderzoekers in staat stelt om deze processen direct in hun eigen analyse-omgevingen te integreren.