Double-real corrections to color singlet decay in a parton-shower inspired scheme

Résumé technique : Corrections de double-réel pour la désintégration de singlets de couleur dans un schéma inspiré des cascades de partons

Énoncé du problème
Les prédictions de haute précision pour les expériences de physique des hautes énergies, particulièrement pour les futures installations comme le FCC-ee, nécessitent des calculs de QCD au niveau Next-to-Next-to-Leading Order (NNLO). Un défi primaire dans ces calculs est le traitement des singularités infrarouges (IR) provenant des divergences de gluons mous et des singularités de collision. Bien que les théorèmes de Bloch-Nordsieck et de Kinoshita-Lee-Nauenberg garantissent l'annulation de ces singularités entre les corrections réelles et virtuelles, le calcul pratique nécessite d'isoler ces dernières pour calculer les restes finis. Les méthodes de soustraction NNLO existantes (par exemple, la décomposition de secteurs, la soustraction d'antennes, CoLoRFulNNLO) se heurtent souvent à la complexité des singularités imbriquées dans les états finaux multi-partoniques et à la difficulté d'associer les résultats d'ordre fixe aux cascades de partons avec une précision Next-to-Leading Logarithmic (NLL). Plus précisément, il existe un besoin pour un schéma de soustraction utilisant des mappages cinématiques simples, compatibles avec les définitions de recul des cascades de partons, et séparant les structures singulières de manière à faciliter l'association avec les calculs de resommation.

Méthodologie
Les auteurs introduisent une méthode de soustraction IR locale pour les corrections de double-réel NNLO dans les désintégrations de singlets de couleur (spécifiquement $e^+e^- \to q\bar{q}$). La méthode est construite sur un algorithme récemment développé pour l'assemblage de fonctions de division (splitting functions) à partir de composantes indépendantes et dépendantes du spin.

• Structure des contre-termes : L'élément de matrice de double-réel carré soustrait ($M_{RR}$) est construit en combinant les termes de soustraction NLO ($S_1$) et NNLO ($S_2$). Le schéma sépare le calcul en termes fortement ordonnés ($S_1 \otimes S_1$), en restes purement NNLO ($S_2 - \sum S_1 \otimes S_1$), et en termes mixtes impliquant des contributions finies à $n+1$ partons.
• Cinématique et choix de jauge : La méthode emploie un mappage de quantité de mouvement basé sur les particules identifiées et les quantités de mouvement de recul, itérant les mappages NLO pour gérer les doubles émissions. Crucialement, les auteurs utilisent un vecteur de jauge de type temps (identifié avec la quantité de mouvement de recul) plutôt qu'une jauge axiale de type lumière. Ce choix s'aligne sur les définitions de recul des algorithmes de cascades de partons sûrs pour le NLL, évitant l'impraticabilité des fonctions de division négatives et de la moyenne azimutale associées aux schémas d'association additifs.
• Désenchevêtrement des singularités imbriquées : Pour traiter les singularités imbriquées dans les radiateurs scalaires (où les gluons sont rayonnés de manière cohérente à partir de multipôles de charge QCD), les auteurs emploient des techniques de fractionnement partiel. Cela permet de désenchevêtrer les composantes accentuées par le régime mou (soft-enhanced), en attribuant les singularités à des configurations de collision uniques.
  • Cas NLO : Le fractionnement partiel est appliqué aux facteurs d'éikonal en utilisant des variables angulaires.
  • Cas NNLO (Abélien et Non-Abélien) : Les radiateurs scalaires pour la double émission de gluons et l'émission de paires $q\bar{q}$ sont décomposés. Les contributions non-abéliennes, qui contiennent des pôles de triple-collision enchevêtrés, sont partitionnées à l'aide de poids angulaires pour les mapper vers des secteurs de collision spécifiques.
• Corrélations de spin : Les auteurs définissent une méthode améliorée pour construire des tenseurs de polarisation à l'aide de vecteurs de polarisation partielle. Cela garantit la cohérence entre les contre-termes sur couche (on-shell) et les termes hors-couche (off-shell), évitant les pôles dépendants de la jauge qui surviennent lors de l'assignation de quantités de mouvement mal adaptée dans les termes de soustraction.

Contributions clés et résultats

• Intégration analytique : Les auteurs ont calculé les intégrales d'espace de phase des contre-termes scalaires analytiquement pour la configuration cinématique dos à dos (back-to-back). En utilisant l'unitarité inversée et le package Kira pour la réduction d'intégrales, ils ont dérivé les structures de pôles ($\epsilon^{-4}$ à $\epsilon^0$) pour les contributions d'émission abélienne, non-abélienne et quark-antiquark.
• Validation numérique :
  • Limites infrarouges : La méthode a été validée à l'aide des générateurs d'événements Pepper et Sherpa. Le reste soustrait ($\Delta$) a été montré comme s'annulant linéairement avec le paramètre d'échelle infrarouge $\lambda$ dans les limites de simple-non-résolu et quadratiquement dans les limites de double-non-résolu, confirmant l'annulation correcte des singularités de premier et second ordre.
  • Intégration numérique : Les auteurs ont vérifié que les intégrales d'espace de phase des radiateurs scalaires peuvent être calculées numériquement via la décomposition de secteurs (pySecDec) combinée à la paramétrisation spécifique de l'espace de phase. Les résultats numériques pour les radiateurs intégrés correspondaient aux résultats analytiques "dos à dos" avec une haute précision.
  • Performance : Une intégration préliminaire de la contribution de double-réel à la section efficace de $e^+e^- \to q\bar{q}$ au pôle Z a produit une incertitude statistique de $654$ fb ($< 2 \cdot 10^{-5}$ par rapport à la section efficace de Born) avec $5,1 \cdot 10^8$ points d'espace de phase.
• Taux de jets : Des taux de jets différentiels ($2 \to 3$ et $3 \to 4$) ont été calculés, démontrant la performance de l'algorithme dans des situations de type NLO et confirmant une convergence numérique suffisante pour la physique de précision.

Signification et affirmations
L'article affirme que ce nouveau schéma fournit une séparation directe des structures singulières et établit une correspondance directe avec les calculs de resommation sous forme de cascades de partons. En utilisant des mappages cinématiques itérés simples qui correspondent aux définitions de recul des cascades de partons sûres pour le NLL, la méthode facilite la combinaison des calculs d'ordre fixe NNLO avec la resommation numérique. Les auteurs postulent que cette approche permettra d'utiliser les contre-termes de double-réel comme corrections locales d'ordre supérieur pour les noyaux de division (splitting kernels) des cascades de partons, aidant ainsi au développement de cascades de partons pleinement différentielles à plus haute précision. Ce travail représente une étape vers la réduction de l'écart de précision entre la théorie et l'expérience anticipé pour les futurs collisionneurs électron-positron, notant spécifiquement les exigences de l'option Tera-Z pour le FCC-ee. Les auteurs restent modestes, notant que les résultats actuels sont physiques sans l'inclusion des corrections réel-virtuelles et double-virtuelles, et que des travaux supplémentaires sont nécessaires pour intégrer les fonctions de radiateur à une boucle et dériver les contre-termes correspondants afin de compléter le calcul complet NNLO pour $e^+e^- \to \text{hadrons}$.