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

Die Autoren implementieren die analytischen Ausdrücke für die QCD-Korrekturen nächster Ordnung zur $gg \to ZH$-Prozess in ggxy, erweitern dies um Zerfälle des Z-Bosons mit Spin-Korrelationen und Off-Shell-Effekten sowie eine Schnittstelle zu POWHEG, um Simulationen mit Parton-Shows in Pythia zu ermöglichen.

Das Wesentliche

Das große Puzzle: Wie man Higgs-Bosonen und Z-Bosonen am besten „fängt"

Stellen Sie sich das Large Hadron Collider (LHC) am CERN wie eine riesige, superschnelle Teilchen-Schlagmaschine vor. Physiker lassen dort Protonen (die aus Quarks bestehen) mit fast Lichtgeschwindigkeit kollidieren, um neue Teilchen zu entdecken oder bekannte zu studieren.

Ein besonders spannendes Ziel ist es, ein Higgs-Boson (das Teilchen, das anderen Teilchen Masse gibt) zusammen mit einem Z-Boson (ein Botenteilchen der schwachen Kraft) zu beobachten. Man nennt das die „ZH-Produktion".

Das Problem ist: Die Natur ist kompliziert. Wenn man diese Kollisionen berechnet, gibt es verschiedene Wege, wie sie passieren können. Die Wissenschaftler haben in dieser Arbeit einen sehr speziellen und schwierigen Weg genauer untersucht: den Weg, bei dem zwei Gluonen (die Klebstoff-Teilchen, die Quarks zusammenhalten) direkt zu einem Higgs und einem Z-Boson werden.

Hier ist die Geschichte der Arbeit, übersetzt in einfache Bilder:

1. Der neue, präzise Bauplan (Die Theorie)

Bisher hatten die Physiker für diesen speziellen Weg (Gluon + Gluon → Higgs + Z) nur eine grobe Schätzung oder Berechnungen, die bei sehr hohen Energien gut funktionierten, aber bei anderen nicht so genau waren.

Die Autoren dieser Arbeit haben nun einen neuen, hochpräzisen Bauplan erstellt.

• Die Metapher: Stellen Sie sich vor, Sie wollen ein Haus bauen. Früher hatten Sie nur eine Skizze, die bei gutem Wetter (bestimmte Energien) funktionierte. Jetzt haben Sie einen detaillierten 3D-Plan, der auch bei Sturm und Regen (komplexe physikalische Bedingungen) perfekt funktioniert.
• Was sie getan haben: Sie haben komplexe mathematische Formeln (die „analytischen Ausdrücke") in ein Computerprogramm namens ggxy eingebaut. Dieses Programm ist wie ein flexibler Werkzeugkasten. Man kann damit berechnen, wie wahrscheinlich es ist, dass diese Kollision passiert, und zwar mit einer Genauigkeit, die man „NLO" (Next-to-Leading Order) nennt. Das bedeutet: Sie berücksichtigen nicht nur den einfachsten Weg, sondern auch die kleinen, störenden Nebeneffekte, die in der echten Welt passieren.

2. Der „Zerfall" und die Spin-Verwirrung (Die Z-Boson-Entscheidung)

Das Z-Boson ist sehr unbeständig. Es zerfällt sofort in andere Teilchen, meist in ein Paar von Leptonen (z. B. ein Elektron und ein Positron oder zwei Neutrinos).

• Das Problem: Wenn man nur das Z-Boson als stabiles Teilchen berechnet, verliert man wichtige Informationen. Es ist, als würde man einen Tanz nur aus der Ferne beobachten, ohne zu sehen, wie die Tänzer ihre Arme bewegen.
• Die Lösung: Die Autoren haben in ihr Programm eingebaut, dass das Z-Boson sofort in seine Zerfallsprodukte umgewandelt wird. Wichtig dabei ist die „Spin-Korrelation".
  • Die Analogie: Stellen Sie sich zwei Eiskunstläufer vor, die sich drehen. Wenn sie sich trennen, fliegen sie in bestimmte Richtungen. Wenn man ihre Rotation (den Spin) ignoriert, denkt man, sie könnten überall hinfliegen. Aber in der Realität bestimmen ihre Drehungen, wohin sie genau fliegen.
  • Die Arbeit zeigt: Wenn man diese Drehungen (Spin) ignoriert, bekommt man falsche Vorhersagen für die Winkel, in denen die Elektronen wegfliegen. Das neue Programm berücksichtigt diese Drehungen automatisch.

3. Der Computer-Check (Validierung)

Bevor man einem neuen Bauplan vertraut, muss man ihn testen.

• Die Autoren haben ihr neues Programm ggxy mit anderen bekannten Programmen (wie vh@nnlo) und mit alten Berechnungen verglichen.
• Das Ergebnis: Alles passt perfekt überein. Das ist wie wenn zwei verschiedene Architekten denselben Turm berechnen und beide genau die gleiche Höhe herausbekommen. Das gibt Sicherheit, dass die neuen Formeln stimmen.

4. Die Simulation mit dem „Parton-Shower" (POWHEG & Pythia)

Bis jetzt haben wir nur über die harte Physik der Kollision gesprochen. Aber in einem echten Detektor passiert noch mehr: Die Teilchen spritzen herum, wie Funken von einem Feuerwerk. Man nennt das „Parton-Shower".

• Die Herausforderung: Man muss die präzise Rechnung (ggxy) mit einer Simulation verbinden, die dieses „Herumfliegen" nachahmt. Dafür nutzen sie ein Werkzeug namens POWHEG und einen Simulator namens Pythia.
• Der Test: Sie haben simuliert, was passiert, wenn man die Z-Boson-Zerfälle vor dem Feuerwerk berechnet (mit Spin) versus wenn man das Feuerwerk erst laufen lässt und dann das Z-Boson zerfallen lässt.
• Das Fazit: Wenn man die Drehung (Spin) des Z-Bosons ignoriert, sieht das Endergebnis im Detektor völlig falsch aus – besonders bei den Winkeln, in denen die Elektronen landen. Der Unterschied ist riesig (bis zu 100% in manchen Bereichen!). Das neue Programm stellt sicher, dass diese Details korrekt simuliert werden.

5. Warum ist das wichtig? (Der große Zusammenhang)

Am Ende des Papers zeigen die Autoren, wie sich dieser spezielle Weg (Gluon-Gluon) auf die Gesamtzahl der ZH-Ereignisse auswirkt.

• Bei niedrigen Energien ist dieser Weg weniger wichtig.
• Aber bei sehr hohen Energien (wie in Zukunft am LHC oder einem zukünftigen Collider) wird dieser Weg immer wichtiger. Er kann bis zu 25% aller Ereignisse ausmachen!
• Ohne diese präzise Rechnung würden die Physiker die Daten falsch interpretieren und vielleicht wichtige Hinweise auf neue Physik übersehen.

Zusammenfassung für den Alltag

Stellen Sie sich vor, Sie sind ein Detektiv, der herausfinden will, wie ein Verbrechen (eine Teilchenkollision) passiert ist.

1. Früher: Sie hatten nur eine grobe Skizze des Tatorts.
2. Jetzt: Sie haben einen hochauflösenden 3D-Scan, der jede kleine Bewegung (Spin) und jeden Winkel berücksichtigt.
3. Der Test: Sie haben diesen Scan mit alten Beweisen abgeglichen – er passt perfekt.
4. Die Simulation: Sie haben einen Film davon gedreht, der zeigt, wie die Tat sich in der Realität abspielt, inklusive aller kleinen Details.
5. Das Ergebnis: Ohne diese Details würden Sie den Täter (die Physik dahinter) falsch verstehen. Mit diesem neuen Werkzeug sind die Physiker jetzt viel besser gerüstet, um die Geheimnisse des Higgs-Bosons zu entschlüsseln.

Die Arbeit ist also im Grunde ein Upgrade für den Werkzeugkasten der Teilchenphysiker, damit sie die Welt der kleinsten Teilchen noch genauer verstehen können.

Technische Zusammenfassung

Titel: $gg \to ZH$ auf NLO-QCD-Niveau abgeglichen mit Parton-Shows unter Verwendung von ggxy und POWHEG

Autoren: Joshua Davies, Kay Schönwald, Matthias Steinhauser, Daniel Stremmer
Institut: Universität Liverpool, CERN, KIT (Karlsruhe)
Datum: März 2026 (Vorlage für JHEP)

---

1. Problemstellung und Motivation

Die assoziierte Produktion eines Higgs-Bosons mit einem Eichboson ($ZH$-Produktion) ist ein entscheidender Prozess zur Untersuchung der Higgs-Eigenschaften am Large Hadron Collider (LHC). Während die $WH$-Produktion primär über Drell-Yan-Prozesse erfolgt, gibt es für die $ZH$-Produktion zusätzlich einen numerisch wichtigen, gluon-induzierten Kanal ($gg \to ZH$).

Bisherige öffentliche Codes wie vh@nnlo und HAWK decken zwar viele Aspekte ab, fehlten jedoch eine Implementierung der gluon-induzierten $gg \to ZH$-Prozesse mit vollständigen nächsten-leading-order (NLO) QCD-Korrekturen in einem flexiblen, modularen Framework. Zudem war die Behandlung von off-shell-Effekten und Spin-Korrelationen bei den Zerfällen des Z-Bosons in diesem spezifischen Kanal nicht vollständig integriert. Das Ziel dieser Arbeit ist es, diese Lücke zu schließen und ein Werkzeug für präzise Vorhersagen bereitzustellen, das sowohl für feste Ordnung (fixed-order) als auch für Simulationen mit Parton-Shows geeignet ist.

2. Methodik und Implementierung

Theoretische Grundlage

Die Berechnung basiert auf analytischen Ausdrücken für die Zwei-Schleifen-Amplituden, die in der Literatur (Ref. [14]) unter Verwendung komplementärer Näherungsmethoden entwickelt wurden:

• Vorwärts-Näherung (Forward expansion): Gilt für $|t|, q_Z^2, q_H^2 \ll m_t^2, s$.
• Hochenergie-Näherung (High-energy expansion): Gilt für $q_Z^2, q_H^2 \ll m_t^2 \ll s, |t|$.
• Große Top-Massen-Näherung: Wird für niedrige Schwerpunktsenergien ($\sqrt{s} < 200$ GeV) genutzt.

Diese Methoden decken den gesamten Phasenraum mit hoher Präzision ab. Die Amplituden werden durch sechs unabhängige Formfaktoren ($F_{12}^{\pm}, F_2^-, F_3^{\pm}, F_4$) ausgedrückt.

Software-Implementierung in ggxy

Die Autoren implementieren diese Ausdrücke in der C++-Bibliothek ggxy.

• Modularität: Das Framework erlaubt die flexible Wahl des Renormierungsschemas für die Top-Quark-Masse (On-Shell oder $\overline{\text{MS}}$) und der Renormierungsskala.
• Formfaktoren: Es werden Padé-Approximationen für beide Expansionen bei $p_T > 180$ GeV verwendet, um die Konvergenz zu verbessern.
• Reale Korrekturen: Die $2 \to 3$ Matrixelemente werden mit Recola berechnet. Ein besonderer Fokus liegt auf der korrekten Filterung von Diagrammen: Drell-Yan-artige Beiträge (Higgs-Strahlung von externen Quarklinien) werden für den reinen gluon-induzierten Kanal ausgeschlossen, während Beiträge mit geschlossenen Quarkschleifen (wobei Z und Higgs an die Schleife koppeln oder nur das Higgs) berücksichtigt werden.
• Off-Shell-Effekte: Die Zerfälle $Z \to \ell^-\ell^+$ und $Z \to \nu\bar{\nu}$ werden auf Amplitudenniveau implementiert, inklusive Spin-Korrelationen und Off-Shell-Effekten (bis ca. 150 GeV Off-Shell-Masse).

Schnittstelle zu POWHEG und Parton-Shows

Ein neuer POWHEG-Prozess (ggxy_ggZH) wurde entwickelt, der ggxy mit Parton-Shows (via Pythia 8) koppelt.

• Stabile vs. instabile Z-Bosonen: Es gibt zwei Implementierungen: eine mit stabilem Z-Boson (Zerfall durch den Shower) und eine mit direkter Zerfallsbehandlung im Matrixelement (inkl. Spin-Korrelationen).
• Matching: Das POWHEG-Verfahren wird verwendet, um die NLO-Berechnung nahtlos mit der Emission weiterer Partonen zu verbinden. Der Parameter $h_{\text{damp}}$ wird genutzt, um die Real-Korrekturen zu regulieren und die Übereinstimmung mit der festen Ordnung bei hohen Impulsen zu optimieren.

3. Wichtige Beiträge und Ergebnisse

Validierung und Kreuzchecks

• Die Implementierung wurde rigoros gegen bekannte Ergebnisse (Ref. [12], [13], [34]) validiert.
• Feste Ordnung: Die berechneten Wirkungsquerschnitte stimmen innerhalb der statistischen Unsicherheiten mit Referenzwerten überein (z. B. $\sigma_{\text{NLO}} \approx 130.5$ fb bei 13.6 TeV).
• POWHEG: Die festen Ordnungsergebnisse aus POWHEG (nach Phase 2) stimmen mit denen aus ggxy überein (Abweichung < 1 Standardabweichung).

Ergebnisse für feste Ordnung (Fixed-Order)

• Top-Massen-Schema: Der Wechsel zwischen On-Shell- und $\overline{\text{MS}}$-Schema führt zu signifikanten Verschiebungen (bis zu 40 % bei hohen invarianten Massen), was eine wichtige theoretische Unsicherheitsquelle darstellt.
• Kombination mit Drell-Yan: Die Autoren kombinieren ihre $gg \to ZH$-Ergebnisse mit den Drell-Yan-Beiträgen aus vh@nnlo.
  • Der gluon-induzierte Kanal trägt bei LHC-Energien (13.6 TeV) etwa 15 % zum Gesamtwirkungsquerschnitt bei.
  • Bei höheren Energien (bis 100 TeV) steigt dieser Anteil auf bis zu 25 % an.
  • Die NLO-Korrekturen für den $gg$-Kanal sind groß (K-Faktor $\approx 2$) und reduzieren die Skalenunsicherheiten von ~25 % auf ~15 %.

Ergebnisse mit Parton-Shows (Pythia)

• Spin-Korrelationen: Der Vergleich zwischen der Behandlung des Z-Zerfalls im Matrixelement (mit Spin-Korrelationen) und im Shower (ohne) zeigt massive Unterschiede in Winkelverteilungen (z. B. Azimutalwinkel $\Delta\phi_{e^-e^+}$) und transversen Impulsen der Leptonen ($p_{T,e^+}$). Die Unterschiede liegen bei 40–100 % und übersteigen die Skalenunsicherheiten bei weitem. Für rekonstruierte Observablen wie $M_{ZH}$ ist der Effekt jedoch gering.
• Shower-Effekte: Die Emission zusätzlicher Partonen durch den Shower führt zu signifikanten Formänderungen in den Verteilungen von $p_{T,ZH}$ und $p_{T,H}$.
  • Bei $p_{T,H}$ kann der Shower-Effekt im Schwanz der Verteilung über 80 % betragen.
  • Die Wahl des $h_{\text{damp}}$-Parameters (z. B. 250 vs. $\infty$) beeinflusst die Härte des Spektrums im Schwanz erheblich. Ein Wert von $h_{\text{damp}} = 250$ reduziert die Diskrepanzen zur festen Ordnung um etwa die Hälfte.

4. Signifikanz und Ausblick

Diese Arbeit stellt einen Meilenstein für die Präzisionsphysik der Higgs-Produktion dar:

1. Verfügbarkeit: Sie macht die hochkomplexen NLO-QCD-Korrekturen für den gluon-induzierten $ZH$-Kanal erstmals in einem öffentlich zugänglichen, modularen Code (ggxy) verfügbar.
2. Flexibilität: Das Framework erlaubt detaillierte Studien zu Renormierungsschemata und Off-Shell-Effekten, die für zukünftige Hochenergie-Experimente (HL-LHC, FCC) essenziell sind.
3. Genauigkeit: Die Integration von Spin-Korrelationen und Off-Shell-Effekten direkt im Matrixelement ist notwendig, um die experimentellen Messungen von Lepton-Winkeln und Impulsen korrekt zu interpretieren.
4. Zukunft: Die Ergebnisse motivieren die Berechnung von NNLO-Korrekturen für den $gg \to ZH$-Kanal im Grenzfall großer Top-Masse, da die Skalenunsicherheiten dort dominieren und die Reskalierung mit dem schweren Top-Quark bereits sehr gute Ergebnisse liefert.

Die bereitgestellten Beispiel-Dateien und Schnittstellen ermöglichen es der Gemeinschaft, diesen Prozess direkt in eigene Analysen zu integrieren.