One-loop amplitudes for $t\bar{t}j$ and $t\bar{t}\gamma$ productions at the LHC through $\mathcal{O}(\epsilon^2)$

Technische Zusammenfassung: Ein-Schleifen-Amplituden für die Produktion von $t\bar{t}j$ und $t\bar{t}\gamma$ am LHC bis zur Ordnung $O(\epsilon^2)$

Problemstellung
Theoretische Vorhersagen für die Produktion von Top-Quark-Paaren in Assoziation mit einem Jet ($t\bar{t}j$) oder einem Photon ($t\bar{t}\gamma$) am Large Hadron Collider (LHC) sind derzeit bis zur Genauigkeit der Next-to-Leading Order (NLO) verfügbar. Um jedoch der beispiellosen Präzision experimenteller Messungen in den aktuellen und kommenden LHC-Läufen gerecht zu werden, müssen Vorhersagen in der Next-to-Next-to-Leading Order (NNLO) der Quantenchromodynamik (QCD) bereitgestellt werden. Die Erreichung der NNLO-Genauigkeit erfordert die Berechnung von Zwei-Schleifen-Streuamplituden für $2 \to 3$-Prozesse. Eine kritische Voraussetzung für die Konstruktion dieser Zwei-Schleifen-Hartfunktionen ist die Verfügbarkeit von Ein-Schleifen-Amplituden, die bis zur Ordnung $O(\epsilon^2)$ in dem Dimensionsregularisierungsparameter $\epsilon$ entwickelt sind. Zwar existieren automatisierte Pakete zur Berechnung des endlichen Teils ($O(\epsilon^0)$) von Ein-Schleifen-Amplituden, doch liefern sie nicht die höheren $\epsilon$-Terme, die für die Subtraktion ultravioletter (UV) und infraroter (IR) Singularitäten sowie zur Extraktion des endlichen Rests in NNLO-Berechnungen notwendig sind. Dieser Artikel schließt diese Lücke, indem er analytische Ausdrücke für die vollfarbigen Ein-Schleifen-Helizitätsamplituden für die Prozesse $pp \to t\bar{t}j$ und $pp \to t\bar{t}\gamma$ bis zur Ordnung $O(\epsilon^2)$ bereitstellt.

Methodik
Die Autoren verwenden ein mehrstufiges Rechenframework, das endliche-Feld-Arithmetik, Tensorzerlegung und die Methode der Pentagon-Funktionen kombiniert:

1. Konstruktion der Helizitätsamplituden: Die Berechnung wird im 't Hooft-Veltman (tHV)-Schema durchgeführt. Die Autoren zerlegen Partialamplituden in unabhängige 4-dimensionale Tensorstrukturen. Massive Spinoren für die Top-Quarks werden unter Verwendung einer Referenzrichtung definiert, um Helizitätsamplituden in Bezug auf masselose Impulse und Spinorprodukte auszudrücken. Die Helizitäts-Subamplituden werden extrahiert, indem die kontrahierten Partialamplituden bei spezifischen Wahlen von Referenzvektoren ausgewertet und das resultierende lineare System invertiert wird.
2. Rekonstruktion über endliche Felder: Um die algebraische Komplexität der Zwischenausdrücke zu bewältigen, nutzen die Autoren Methoden endlicher Felder. Feynman-Diagramme werden mit QGRAF generiert, gefolgt von einer Farbzerlegung. Die resultierenden kontrahierten Partialamplituden werden mittels Integration-by-Parts (IBP)-Identitäten (generiert via NEATIBP und gelöst via FINITEFLOW) auf eine minimale Menge von Master-Integralen (MIs) reduziert. Die rationalen Koeffizienten, die die MIs multiplizieren, werden analytisch aus ihren numerischen Auswertungen über endliche Felder rekonstruiert.
3. Basis der Pentagon-Funktionen: Die MIs werden in Bezug auf eine Basis algebraisch unabhängiger Komponenten ausgedrückt, die als „Pentagon-Funktionen" bekannt sind. Diese Funktionen entsprechen den in $\epsilon$ entwickelten Komponenten der MIs. Die Autoren leiten kanonische Differentialgleichungen (DEs) für diese Funktionen in $d\log$-Form ab. Das Alphabet der DEs besteht aus 283 unabhängigen Buchstaben, einschließlich 29 algebraischer Quadratwurzeln.
4. Numerische Auswertung: Die DEs werden numerisch unter Verwendung der Methode der verallgemeinerten Potenzreihenentwicklung gelöst. Die Autoren nutzen zwei Softwarepakete: DIFFEXP (Mathematica) und das neu veröffentlichte LINE (C/C++). Randbedingungen werden durch Auswertung der MIs an spezifischen Phasenraum-Punkten unter Verwendung des AMFLOW-Pakets mit 70-stelliger Genauigkeit bestimmt. Die Lösung wird auf beliebige physikalische Phasenraum-Punkte fortgesetzt.
5. Renormierung und endlicher Rest: Die nackten Amplituden werden massenrenormiert und vollständig UV-renormiert, indem geeignete Gegenterme hinzugefügt werden. Der endliche Rest wird durch Subtraktion universeller IR-Singularitäten unter Verwendung des MS-Schemas extrahiert.

Hauptbeiträge

• Analytische Ausdrücke: Der Artikel liefert die ersten analytischen Ausdrücke für die vollfarbigen Ein-Schleifen-Helizitätsamplituden für die Produktion von $t\bar{t}j$ und $t\bar{t}\gamma$, entwickelt bis zur Ordnung $O(\epsilon^2)$.
• Basis der Pentagon-Funktionen: Die Autoren konstruieren eine spezifische Basis von Pentagon-Funktionen für diese Prozesse und drücken die Amplituden als Linearkombinationen dieser Funktionen mit rationalen Koeffizienten aus, die in Impuls-Twistor-Variablen geschrieben sind.
• Differentialgleichungen: Kanonische DEs für die Pentagon-Funktionen werden hergeleitet und numerisch gelöst. Die Autoren stellen die DEs, Randwerte und Permutationsregeln für die Funktionen in ancillären Dateien bereit.
• Software-Implementierung: Die Arbeit demonstriert die Anwendung des Pakets LINE zur Lösung von DEs, die Quadratwurzeln beinhalten, und zeigt, dass es für dieses spezifische Gleichungssystem etwa doppelt so schnell ist wie DIFFEXP.
• Benchmark-Ergebnisse: Numerische Benchmark-Ergebnisse für farb- und helizitätssummierte quadrierte Matrixelemente werden für alle relevanten partonischen Kanäle ($gg \to t\bar{t}g$, $q\bar{q} \to t\bar{t}g$, $gg \to t\bar{t}\gamma$, $q\bar{q} \to t\bar{t}\gamma$) bereitgestellt und gegen OPENLOOPS bis zur Ordnung $O(\epsilon^0)$ validiert.

Ergebnisse
Die Autoren haben erfolgreich die Helizitätsamplituden für Prozesse berechnet, die Top-Quark-Paare und masselose Teilchen (Gluonen, Photonen, leichte Quarks) beinhalten. Die analytischen Ausdrücke für die rationalen Koeffizienten der Pentagon-Funktionen wurden rekonstruiert, wobei die Nennergrade durch univariate Partialbruchzerlegung signifikant reduziert wurden. Die numerische Auswertung der Amplituden an einem spezifischen physikalischen Phasenraum-Punkt bestätigte die Gültigkeit der Ergebnisse gegenüber bestehenden NLO-Werkzeugen. Der Artikel stellt fest, dass zwar die rationalen Koeffizienten für $t\bar{t}j$ und $t\bar{t}\gamma$ einige strukturelle Ähnlichkeiten aufweisen, eine vollständige Extraktion von Beziehungen zwischen den beiden Prozessen auf der Ebene der primitiven Amplituden jedoch nicht durchgeführt wurde.

Bedeutung und Behauptungen
Die Autoren stellen fest, dass ihre Ergebnisse für den Aufbau der Zwei-Schleifen-Hartfunktionen unerlässlich sind, die für NNLO-QCD-Berechnungen der Produktion von $t\bar{t}j$ und $t\bar{t}\gamma$ erforderlich sind. Durch die Bereitstellung von Amplituden bis zur Ordnung $O(\epsilon^2)$ ermöglichen sie die analytische Subtraktion von UV- und IR-Divergenzen sowie die Extraktion endlicher Reste, die mit Werkzeugen, die auf $O(\epsilon^0)$ beschränkt sind, sonst nicht zugänglich wären.

Der Artikel behauptet, dass die Darstellung der Amplituden in Bezug auf eine Basis von Pentagon-Funktionen einen kompakten Rahmen bietet, der die endlichen Ausdrücke vereinfacht und die numerische Effizienz verbessert. Darüber hinaus positionieren die Autoren diese Arbeit als einen Meilenstein für zukünftige Zwei-Schleifen-Berechnungen. Sie stellen fest, dass zwar führende-farbige Zwei-Schleifen-Amplituden für $t\bar{t}j$ elliptische Strukturen beinhalten, die hier konstruierte Basis von Pentagon-Funktionen (die polylogarithmische und nicht-polylogarithmische Beiträge trennt) eine nützliche Strategie für den Umgang mit den komplexeren funktionalen Räumen bieten könnte, die in Zwei-Schleifen-Berechnungen erwartet werden, einschließlich jener für $t\bar{t}\gamma$. Die Autoren schließen bescheiden, dass, obwohl die Basis zukünftige Arbeiten erleichtert, die vollständige Realisierung von Zwei-Schleifen-Berechnungen Gegenstand zukünftiger Forschung bleibt.