Pseudoscalar Higgs boson decay to three parton amplitudes at NNLO to higher orders in the dimensional regulator

Technische Samenvatting: Pseudo-scalaire Higgs-decay naar drie parton-amplitudes op NNLO naar hogere orden in de dimensionale regulator

Probleemstelling
De precieze voorspelling van de productie van het pseudo-scalaire Higgs-boson ($A$) met een jet bij hadronenversnellers vereist berekeningen voorbij de Next-to-Next-to-Leading Order (NNLO) in Kwantumchromodynamica (QCD). Specifiek vereist de berekening van de Next-to-Next-to-Next-to-Leading Order (N3LO) dwarsdoorsnede de onderliggende lagere-orde virtuele correcties, geëxpandeerd naar hogere machten van de dimensionale regulator $\epsilon$ (waarbij $d = 4 - 2\epsilon$). Hoewel de twee-lus resultaten voor het Standaardmodel Higgs + 1 jet proces, geëxpandeerd tot $\mathcal{O}(\epsilon^2)$, recentelijk zijn berekend, ontbraken de overeenkomstige bijdragen voor het pseudo-scalaire Higgs-boson. Dit gat belemmert de volledige N3LO-berekening voor pseudo-scalaire productie, wat cruciaal is voor het verminderen van schaalonzekerheden en het onderzoeken van de CP-eigenschappen van de Higgs-sector in scenario's buiten het Standaardmodel (BSM), zoals het Minimal Supersymmetric Standard Model (MSSM).

Methodologie
De auteurs voeren een twee-lus berekening uit voor het verval van een pseudo-scalaire Higgs-boson naar drie partonen ($A \to ggg$ en $A \to q\bar{q}g$) binnen een Effectieve Veldtheorie (EFT) raamwerk, waarbij de zware topquark is geïntegreerd. De berekening verloopt als volgt:

1. Theoretisch Raamwerk: De interactie wordt beschreven door een effectief Lagrangiaan met twee operatoren, $O_G$ en $O_J$, gekoppeld via Wilson-coëfficiënten $C_G$ en $C_J$. De berekening maakt gebruik van de 't Hooft-Veltman definitie van $\gamma_5$ om chirale grootheden te behandelen in $d \neq 4$ dimensies.
2. Diagramgeneratie en Algebraïsche Manipulatie: Feynman-diagrammen worden gegenereerd met QGRAF. De resulterende amplitudes worden verwerkt met in-house FORM-codes om Dirac- en $SU(N)$ kleuralgebra uit te voeren.
3. Integraalreductie: Het enorme aantal scalaire Feynman-integralen wordt gereduceerd tot een minimale set van Master Integrals (MIs) met behulp van de publieke package KIRA, die het Laporta-algoritme implementeert via Integration-by-Parts (IBP) en Lorentz-invariantie identiteiten. De auteurs maken gebruik van MIs die analytisch zijn berekend in recent werk [68], uitgedrukt in termen van Generalised Polylogarithms (GPLs), wat de eerdere resultaten die beperkt waren tot Harmonic Polylogarithms (HPLs) bij lagere orden in $\epsilon$, uitbreidt.
4. Renormalisatie:
   • UV Renormalisatie: De koppelingsconstante en samengestelde operatoren worden gerenormaliseerd. Er wordt speciale aandacht besteed aan de renormalisatie van de axiale singlet-stroom ($O_J$), wat een eindige renormalisatieconstante $Z_5^s$ vereist om Ward-identiteiten te herstellen en de menging tussen operatoren $O_G$ en $O_J$ te behandelen.
   • IR Regularisatie: De infrarood (IR) singulariteiten worden geverifieerd tegen de universele subtractie-operatoren voorspeld door Catani [69] en gegeneraliseerd door Becher en Neubert [95].
5. Numerieke Implementatie: De eindelijke delen van de amplitudes worden vereenvoudigd met Mathematica en geoptimaliseerd met FORM (gebruikmakend van Horner-schema's en Monte-Carlo boomzoekopdrachten) om efficiënte FORTRAN-95 routines te genereren. Deze routines zijn ontworpen om te koppelen met elke Monte-Carlo fase-ruimte generator.

Belangrijkste Bijdragen

• Eerste Berekening van Hogere-Orde $\epsilon$ Termen: Het artikel presenteert voor het eerst de $\mathcal{O}(\epsilon^1)$ en $\mathcal{O}(\epsilon^2)$ termen van de matrixelementen op NNLO, evenals de $\mathcal{O}(\epsilon^3)$ en $\mathcal{O}(\epsilon^4)$ termen op NLO, voor de processen $A \to ggg$ en $A \to q\bar{q}g$.
• Analytische Identiteiten: De auteurs bieden een uitgebreide set analytische identiteiten die Harmonic Polylogarithms (HPLs) mappen naar Generalised Polylogarithms (GPLs), wat de vergelijking met eerdere resultaten en de evaluatie van de nieuwe amplitudes vergemakkelijkt.
• Validatie: De resultaten worden gevalideerd door de IR-polstructuur te controleren tegen de voorspellingen van Catani en door de eindelijke delen analytisch te mappen naar eerdere berekeningen [76], waarbij een perfecte overeenkomst werd gevonden.
• Publieke Code: De eindige delen van de amplitudes zijn geïmplementeerd in geoptimaliseerde FORTRAN-95 routines en beschikbaar gesteld als bijlagen, klaar voor gebruik in fenomenologische studies.

Resultaten

• Numerieke Prestaties: De auteurs rapporteren de CPU-tijd die nodig is om de geoptimaliseerde matrixelementen per fase-ruimte punt te evalueren. Op NNLO-niveau is de evaluatietijd voor $\mathcal{O}(\epsilon^2)$ termen ongeveer 18–25 seconden per punt, wat aanzienlijk langer is dan de $\mathcal{O}(\epsilon^4)$ termen op NLO (die $\sim$1 seconde duren). Het optimalisatieproces voor de $\mathcal{O}(\epsilon^2)$ NNLO-resultaten vereiste ongeveer 16 uur vanwege de verhoogde complexiteit en het verschijnen van hogere-gewicht GPLs.
• Grootte van de Bijdragen: Het reële deel van de $\mathcal{O}(\epsilon^2)$ bijdrage voor het gluonische kanaal ($S_{ggg}$) wordt waargenomen als ongeveer 100 keer groter dan die van het quark-gluon kanaal ($S_{q\bar{q}g}$).
• Cross-kanaal Relaties: De resultaten zijn afgeleid met oog voor crossing-symmetrie relaties, waardoor ze direct toepasbaar zijn op de berekening van de productie van het pseudo-scalaire Higgs-boson geassocieerd met een jet ($pp \to A + \text{jet}$).

Betekenis
Het artikel stelt dat deze resultaten een noodzakelijk ingrediënt zijn voor de volledige drie-lus (N3LO) dwarsdoorsnede-berekening voor de productie van het pseudo-scalaire Higgs-boson. Door de amplitudes te leveren die geëxpandeerd zijn naar hogere orden in de dimensionale regulator, maakt dit werk de annulatie van IR-singulariteiten op N3LO-niveau mogelijk. De auteurs benadrukken dat hoewel het SM Higgs-boson als CP-even is vastgesteld, de zoektocht naar CP-danse scalaire deeltjes cruciaal blijft. De beschikbaarheid van deze hogere-orde correcties is essentieel voor het verminderen van theoretische onzekerheden in de differentiële distributies (zoals transversale impuls) van pseudo-scalaire Higgs-bosonen bij de LHC, wat de precisie van tests voor BSM-fysica verbetert. Het werk stelt geen nieuwe experimentele opstellingen voor, maar biedt de theoretische infrastructuur die nodig is om toekomstige hoog-precieze collider-data te interpreteren.