Planar master integrals for two-loop NLO electroweak light-fermion contributions to $g g \rightarrow Z H$

Cet article présente un calcul analytique des intégrales maîtresses pour les contributions des fermions légers planaires aux corrections électrofaibles d'ordre NLO à deux boucles de $gg \rightarrow ZH$, en utilisant un cadre d'équations différentielles canoniques pour exprimer la plupart des résultats en termes de polylogarithmes de Goncharov tout en traitant les racines carrées imbriquées restantes par des intégrales à une seule variable.

L'essentiel

Imaginez l'univers comme une machine géante et complexe où de minuscules particules entrent constamment en collision et se transforment. L'une des tâches les plus importantes au Grand collisionneur de hadrons (LHC) consiste à faire entrer en collision des particules pour créer une combinaison spécifique et rare : un boson Z (un porteur lourd de l'interaction faible) et un boson de Higgs (la particule qui donne leur masse aux autres).

Bien que la plupart de ces collisions se produisent de manière directe, il existe un canal secondaire sournois et compliqué où deux « gluons » invisibles (des particules qui maintiennent les noyaux atomiques ensemble) entrent en collision pour créer cette paire Z-Higgs. Ce processus est comparable à une entrée secrète par une porte dérobée dans la machine. Même si cela se produit moins souvent que par la porte principale, c'est suffisamment significatif pour que, si nous l'ignorions, notre carte du fonctionnement de l'univers soit légèrement erronée.

Cet article traite du calcul des « plans » de cette entrée secrète par porte dérobée avec une extrême précision. Voici la décomposition de ce que les auteurs ont réalisé, en utilisant des analogies simples :

1. Le Problème : Un Labyrinthe de Possibilités Infinies

Lorsque les physiciens tentent de calculer ce qui se produit lorsque des particules entrent en collision, ils doivent tenir compte de chaque manière possible dont les particules peuvent osciller, former des boucles et interagir pendant la fraction de seconde de l'impact. Ces interactions sont représentées par des diagrammes de Feynman (pensez-y comme des organigrammes du trafic des particules).

Pour cette collision spécifique ($gg \to ZH$), il existe 132 organigrammes différents (diagrammes) impliquant des particules légères (comme les électrons et les quarks légers) formant des boucles. Tenter de résoudre les mathématiques pour les 132 à la fois, c'est comme essayer de boire à un tuyau d'incendie ; c'est trop désordonné.

2. La Solution : Trouver les « Clés Maîtresses »

Les auteurs ont réalisé que les 132 organigrammes sont en réalité construits à partir d'un ensemble plus restreint de blocs de construction fondamentaux. Ils ont utilisé un outil mathématique appelé Intégration par Parties (IBP) pour décomposer le problème massif.

Pensez-y comme à un château de Lego complexe. Vous n'avez pas besoin de calculer la forme de chaque brique individuellement. À la place, vous identifiez les Intégrales Maîtresses (IM) — les formes de briques uniques et essentielles qui, combinées de différentes manières, permettent de construire l'ensemble du château.

• Ils ont découvert que pour les diagrammes « planaires » (plats, non emmêlés), il existe 62 clés maîtresses uniques pour un type d'interaction et 59 pour un autre.
• Une fois que vous connaissez la valeur de ces clés maîtresses, vous pouvez instantanément déterminer la valeur de l'ensemble du château.

3. La Méthode : La Carte « Canonique »

Pour résoudre ces clés maîtresses, les auteurs ont utilisé une technique appelée Méthode des Équations Différentielles Canoniques.

• L'Analogie : Imaginez que vous êtes perdu dans une forêt brumeuse (le problème mathématique). Vous savez que les arbres (les variables) changent, mais vous ne connaissez pas le chemin. Au lieu de deviner, ils ont construit une carte GPS parfaite (la base canonique) qui vous indique exactement comment le chemin change au fur et à mesure de votre déplacement.
• Ils ont utilisé une astuce mathématique appelée développement de Magnus pour redresser la carte. Cela a transformé un ensemble désordonné et emmêlé d'équations en une liste propre et ordonnée où chaque étape est prévisible.

4. L'Obstacle : Les « Racines Carrées Emboîtées »

Alors qu'ils tentaient d'écrire les réponses finales, ils ont heurté un mur. Les mathématiques impliquaient des racines carrées (comme $\sqrt{2}$ ou $\sqrt{x}$).

• Dans les cas simples, vous pouvez vous débarrasser facilement de ces racines carrées, transformant la réponse en une liste nette de fonctions standard (appelées Polylogarithmes de Goncharov ou GPL). Pensez-y comme des « mots » standard dans le langage de la physique.
• Cependant, dans ce problème spécifique, certaines racines carrées étaient emboîtées à l'intérieur d'autres racines carrées (comme une poupée russe). C'était comme essayer de démêler un nœud où le fil est enroulé autour de lui-même d'une manière qui rend impossible de le tirer droit d'un seul coup.
• Le Résultat : Pour la plupart des clés maîtresses, ils ont trouvé une solution « mot » propre. Mais pour quelques-unes des plus compliquées (celles avec les nœuds emboîtés), ils n'ont pas pu les démêler complètement. Au lieu de cela, ils ont dû les laisser sous forme d'intégrales à une boucle.
  • Analogie : Au lieu de vous donner une phrase terminée, ils vous ont donné une phrase avec un emplacement « à compléter » qui nécessite un petit calcul spécifique pour être achevé. Ce n'est pas un mot complet et propre, mais c'est une instruction précise sur la façon de terminer la phrase.

5. La Vérification : Le « Double Contrôle »

Pour s'assurer qu'ils n'avaient pas fait d'erreur dans leur algèbre complexe, ils ont comparé leurs « plans » manuscrits à une simulation par supercalculateur appelée AMFlow.

• Ils ont choisi un point de test spécifique dans la « région euclidienne » (une zone théorique sûre où les mathématiques sont stables) et ont exécuté les calculs.
• Le Résultat : Leurs formules analytiques correspondaient parfaitement aux résultats numériques de l'ordinateur, jusqu'à 30 décimales. C'est l'équivalent mathématique de deux personnes mesurant une table et s'accordant sur la longueur jusqu'à la largeur d'un atome.

Résumé

Cet article ne nous dit pas comment construire un nouvel accélérateur de particules ou guérir une maladie. Au lieu de cela, il fournit les ingrédients mathématiques essentiels et de haute précision nécessaires pour comprendre une collision de particules spécifique et rare au LHC.

En résolvant les « intégrales maîtresses » pour les contributions des fermions légers, les auteurs ont dissipé le brouillard d'une partie spécifique du Modèle Standard. Ils ont fourni les formules exactes dont les physiciens ont besoin pour prédire ce qui se produit lorsque des gluons créent un boson Z et un boson de Higgs, garantissant que les expériences futures pourront repérer les moindres écarts qui pourraient indiquer une nouvelle physique au-delà de ce que nous connaissons actuellement.

Résumé technique

1. Énoncé du problème

L'article traite du calcul des corrections électrofaibles (EF) à l'Ordre Next-to-Leading (NLO) pour la production associée d'un boson de Higgs ($H$) et d'un boson $Z$ via la fusion de gluons ($gg \to ZH$) au Grand Collisionneur de Hadrons (LHC).

• Contexte : Bien que les corrections QCD à ce processus soient bien étudiées, les corrections EF à l'ordre NLO restent largement non calculées. Ces corrections sont cruciales pour les prédictions de haute précision futures, car elles contribuent à environ 4 % à la section efficace (sur la base de processus analogues) et deviennent significatives dans le régime boosté.
• Défi spécifique : Le calcul implique des diagrammes virtuels à deux boucles induits par des boucles de fermions légers. Ces diagrammes introduisent une nouvelle échelle de masse ($m_W$) aux côtés des invariants cinématiques ($s, t, u$) et d'autres masses ($m_Z, m_H$), conduisant à des intégrales de Feynman complexes.
• Portée : Les auteurs se concentrent spécifiquement sur les topologies planaires parmi les huit topologies distinctes issues des boucles de quarks légers (quatre impliquant le vertex $WWH$ et quatre impliquant le vertex $ZZH$). Les topologies non planaires sont exclues de ce calcul analytique spécifique.

2. Méthodologie

Les auteurs emploient la Méthode des Équations Différentielles Canoniques (ED) pour résoudre les intégrales maîtresses (IM) de manière analytique. Le flux de travail implique :

• Réduction des intégrales : En utilisant le package Kira (implémentant l'algorithme de Laporta), les auteurs réduisent les 132 diagrammes de Feynman irréductibles à deux boucles à un ensemble fini d'Intégrales Maîtresses.
  • Branche B1 (vertex WWH) : 62 IM.
  • Branche B2 (vertex WWH) : 59 IM.
• Construction d'une base canonique : Au lieu de résoudre les ED dans une base générique, ils construisent une base canonique $\vec{g}$ où les équations différentielles prennent la forme factorisée en $\epsilon$ :
  $$d\vec{g} = \epsilon \, dA(\vec{x}) \, \vec{g}$$
  Cela est réalisé en utilisant la méthode de l'expansion de Magnus, à partir d'une base linéaire d'intégrales.
• Alphabet et lettres de symbole : La matrice $dA$ est exprimée en termes de formes $d\log$, $dA = \sum M_a d\log(\omega_a)$. L'ensemble des fonctions $\omega_a$ (lettres de symbole) constitue l'"alphabet" de la solution.
  • Les alphabets impliquent des radicaux non triviaux (racines carrées) découlant de la cinématique et des échelles de masse.
  • B1 : Implique deux racines carrées, $r_1$ et $r_2$.
  • B2 : Implique cinq racines carrées, y compris des racines carrées imbriquées ($r_4, r_5$), rendant la rationalisation simultanée impossible.
• Analyse de la structure radicaire : Pour gérer les radicaux, les auteurs introduisent un indice radical (un vecteur binaire) pour classifier les IM en fonction de leur dépendance à des racines carrées spécifiques. Ils construisent des "arbres tropicaux" pour déterminer les sous-systèmes minimaux requis pour chaque IM.
• Rationalisation et solution :
  • Pour les sous-systèmes à structures radicales simples, des transformations de variables sont appliquées pour rationaliser les racines carrées.
  • Les solutions sont exprimées en termes de Polylogarithmes de Goncharov (GPL) jusqu'à l'ordre $\mathcal{O}(\epsilon^4)$.
  • Pour les sous-systèmes les plus complexes (spécifiquement dans B2 à $\mathcal{O}(\epsilon^3)$ et $\mathcal{O}(\epsilon^4)$), où les radicaux imbriqués empêchent une rationalisation complète, les résultats sont représentés comme des intégrales à une seule boucle sur des GPL.
• Conditions aux limites : Les constantes d'intégration sont fixées en utilisant :
  • Des conditions de régularité aux points singuliers (singularités de Landau et singularités spuriées).
  • Des conditions d'annulation à des limites cinématiques spécifiques.
  • Des symétries de permutation entre les topologies.
  • Des limites analytiques (par exemple, $m_W \to m_Z$) et un ajustement numérique via l'algorithme PSLQ.

3. Contributions clés

1. Calcul analytique : La première évaluation analytique des intégrales maîtresses planaires pour les corrections EF à l'ordre NLO à deux boucles de $gg \to ZH$ impliquant des fermions légers.
2. Gestion des radicaux imbriqués : Un cadre systématique pour traiter les intégrales contenant des racines carrées imbriquées. Les auteurs démontrent que, bien que la rationalisation complète soit impossible pour la branche B2 aux ordres supérieurs, les intégrales peuvent toujours être systématiquement réduites à des intégrales à une seule boucle sur des GPL.
3. Bases canoniques : Construction explicite de bases canoniques pour deux familles d'intégrales complexes (B1 et B2) contenant respectivement 62 et 59 IM.
4. Extension au vertex ZZH : Déduction des IM canoniques pour le vertex $ZZH$ (Branches C1 et C2) en prenant la limite $m_W \to m_Z$ ($z \to -1$) des résultats $WWH$, en traitant analytiquement les divergences résultantes.
5. Résultats de haute précision : Fourniture d'expressions analytiques jusqu'à $\mathcal{O}(\epsilon^4)$, ce qui est nécessaire pour le calcul complet EF à l'ordre NLO (puisque l'amplitude à l'arbre est $\mathcal{O}(\epsilon^0)$, le développement en boucles nécessite des ordres supérieurs en $\epsilon$).

4. Résultats

• Branche B1 (WWH) : Les 62 IM sont classées en quatre sous-ensembles basés sur les structures radicales ($00, 10, 01, 11$). Toutes sont exprimables en termes de GPL après des transformations de rationalisation appropriées.
• Branche B2 (WWH) : Les 59 IM sont classées en trois sous-ensembles.
  • Les sous-ensembles avec des structures radicales simples sont exprimés sous forme de GPL.
  • Le sous-ensemble avec la structure la plus complexe ($S_2(11111)$) contient des IM qui, à $\mathcal{O}(\epsilon^3)$ et $\mathcal{O}(\epsilon^4)$, nécessitent une représentation sous forme d'intégrales à une seule boucle sur des GPL en raison de la présence de radicaux imbriqués ($r_4, r_5$) qui ne peuvent être rationalisés simultanément.
• Branches C1/C2 (ZZH) : Les auteurs déduisent avec succès les expressions analytiques pour le vertex $ZZH$ en prenant la limite de masse, vérifiant la cohérence des parties singulières et la régularité des intégrales restantes.
• Validation numérique : Les résultats analytiques ont été vérifiés par rapport à des calculs numériques indépendants utilisant le package AMFlow (méthode du flux de masse auxiliaire). L'accord a été vérifié à 30 chiffres significatifs en un point euclidien de référence, confirmant la justesse des expressions analytiques.

5. Importance

• Phénoménologie de précision : Ces résultats fournissent les éléments constitutifs essentiels (Intégrales Maîtresses) nécessaires pour calculer les corrections EF complètes à l'ordre NLO pour $gg \to ZH$. Ceci est vital pour réduire les incertitudes théoriques en physique du Higgs au LHC et aux futurs collisionneurs.
• Avancement méthodologique : L'article démontre une stratégie robuste pour aborder les intégrales multi-boucles avec plusieurs échelles de masse et radicaux imbriqués. L'approche des "arbres tropicaux" pour classer les dépendances radicales et l'utilisation d'intégrales à une seule boucle pour les cas non rationalisables offrent un modèle pour des calculs complexes similaires dans le Modèle Standard et au-delà.
• Sondes de nouvelle physique : En améliorant la précision théorique de la section efficace $gg \to ZH$, ces calculs renforcent la capacité de détecter des déviations par rapport au Modèle Standard, en particulier dans les couplages Higgs-boson de jauge et les contributions potentielles de nouvelle physique dans la boucle.

En résumé, ce travail représente une avancée significative dans la description théorique de la production de Higgs, comblant le fossé entre les calculs de précision QCD et électrofaibles pour un processus clé du LHC.