QCD corrections for Pseudoscalar Higgs decay to 3 partons at higher orders in dimensional regulator

Cette contribution présente une étude des corrections d'ordre supérieur pour la désintégration du Higgs pseudoscalaire en trois partons, fournissant des résultats essentiels pour prédire la distribution différentielle de ce boson associé à un jet dans les collisionneurs hadroniques.

L'essentiel

🎭 Le Higgs Mystérieux et la Danse des Particules

Imaginez que l'Univers est une immense scène de théâtre. En 2012, les physiciens ont découvert un nouvel acteur sur cette scène : le Boson de Higgs. C'est une étoile, mais il y a un mystère : est-il un homme (symétrique) ou une femme (antisymétrique) ? En physique, on parle de propriétés "CP".

Dans cette pièce, il existe un personnage spécial appelé le Higgs Pseudoscalaire (noté A). C'est le "jumeau sombre" ou la version "antisymétrique" du Higgs que nous connaissons. Les scientifiques veulent savoir si le Higgs découvert en 2012 est ce personnage A ou non. Pour le prouver, ils doivent observer comment ce personnage A se comporte quand il se désintègre.

🎱 Le Problème de la Désintégration à 3 Particules

Le papier dont nous parlons étudie un scénario très précis : comment le Higgs A se transforme en trois autres particules (soit trois gluons, soit un quark, un anti-quark et un gluon).

Imaginez que le Higgs A est un magicien qui lance un sort. Au lieu de faire apparaître un seul objet, il en fait apparaître trois qui partent dans toutes les directions.

• Le défi : Calculer exactement comment ces trois objets s'envolent est extrêmement difficile. C'est comme essayer de prédire la trajectoire de trois balles de billard qui rebondissent les unes sur les autres, tout en tenant compte de l'air, du frottement et de la gravité, mais à l'échelle subatomique.

🔍 Pourquoi aller si loin dans les détails ? (Les "Corrections QCD")

En physique, on ne se contente pas d'une première approximation. On veut une précision chirurgicale.

• Niveau 1 (LO) : C'est le dessin au crayon. On voit la forme générale.
• Niveau 2 (NLO) : On ajoute les couleurs.
• Niveau 3 (NNLO) : On ajoute les ombres et les lumières.
• Ce papier (Au-delà de NNLO) : Les auteurs vont encore plus loin. Ils calculent des corrections de troisième ordre (et même plus).

Pourquoi ? Parce que dans les collisionneurs de particules (comme le LHC), les signaux sont très faibles. Si vous voulez détecter ce Higgs A parmi des milliards de collisions, vous devez savoir exactement à quoi ressemblera le signal "parfait". Si votre calcul théorique est imprécis, vous risquez de confondre le signal avec le bruit de fond.

🛠️ L'Outil Magique : Le "Régulateur Dimensionnel"

Pour faire ces calculs, les physiciens utilisent une astuce mathématique appelée le régulateur dimensionnel.

• L'analogie : Imaginez que vous essayez de mesurer la surface d'un nuage. En 2D, c'est flou. En 3D, c'est mieux. Mais pour faire les maths, les physiciens disent : "Et si on travaillait dans un monde à 4, 5 ou 6 dimensions ?"
• Cela permet de gérer les "infinis" qui apparaissent dans les calculs (comme des trous noirs mathématiques). Ils travaillent dans un monde imaginaire à $d = 4 + \epsilon$ dimensions, font leurs calculs, puis ramènent le résultat dans notre monde à 4 dimensions.
• Dans ce papier, les auteurs poussent cette astuce très loin, en développant leurs formules jusqu'à des puissances très élevées de $\epsilon$ (comme $\epsilon^4$), ce qui est un exploit mathématique.

🧮 Le Travail de Géant : FORM et les Ordinateurs

Les équations générées par ces calculs sont monstrueuses.

• L'analogie : C'est comme si vous deviez résoudre une équation qui prendrait toute la bibliothèque de Babel pour être écrite. Les fichiers de calcul pèsent plusieurs Gigaoctets (des centaines de milliers de pages de formules).
• Pour gérer cela, les auteurs utilisent un logiciel spécial appelé FORM, conçu pour manipuler des symboles mathématiques gigantesques.
• Ils ont dû optimiser ces calculs pour qu'ils ne prennent pas des années à tourner sur un ordinateur. Grâce à des astuces de programmation, ils ont réduit le temps de calcul d'une configuration à quelques secondes ou microsecondes, ce qui permet de les utiliser dans des simulations réalistes.

🚀 À Quoi ça sert ? (Le Résultat)

Le but final n'est pas juste de faire de belles maths.

1. Prédiction précise : Ces résultats permettent de prédire exactement combien de Higgs A devraient être produits avec un "jet" (un groupe de particules) dans les collisionneurs.
2. Chasse au trésor : En comparant ces prédictions ultra-précises avec les données réelles du LHC, les physiciens pourront dire : "Ah ! Le Higgs que nous voyons correspond parfaitement à la théorie du Higgs Pseudoscalaire" ou "Non, ce n'est pas lui".
3. Au-delà du Modèle Standard : Si le Higgs A existe, cela ouvre la porte à de nouvelles physiques (comme la Supersymétrie), qui pourraient expliquer des mystères comme la masse des neutrinos ou l'unification des forces.

En Résumé

Ce papier est un manuel de précision extrême. Les auteurs ont construit une "machine à calculer" capable de prédire le comportement d'un Higgs hypothétique avec une précision inégalée. C'est comme passer d'une carte routière approximative à un GPS satellite capable de vous dire exactement où poser le pied, même sur un terrain glissant. Cela aide les physiciens à ne pas rater la découverte qui pourrait changer notre compréhension de l'Univers.

Résumé technique

Voici un résumé technique détaillé de l'article de Pulak Banerjee et al., intitulé « Corrections QCD pour la désintégration du Higgs pseudoscalaire en 3 partons à des ordres supérieurs dans le régulateur dimensionnel ».

1. Problématique et Contexte

Depuis la découverte du boson de Higgs en 2012 au LHC, la détermination de ses propriétés de parité (CP) et la mesure de ses couplages sont devenues cruciales. Bien que le Modèle Standard (MS) soit performant, il laisse des questions ouvertes (masse des neutrinos, problème de la hiérarchie), ce qui motive l'étude de modèles au-delà du Modèle Standard (BSM), notamment le Modèle Standard Supersymétrique Minimal (MSSM).

Dans le MSSM, le secteur scalaire comprend trois scalaires neutres ($h, H, A$) et deux scalaires chargés ($H^\pm$). Le scalaire $A$ est un pseudoscalaire CP-impair. La production de ce pseudoscalaire au LHC se fait principalement par fusion de gluons, un mécanisme similaire à celui du Higgs du MS, mais impliquant des opérateurs effectifs différents ($O_G$ et $O_J$).

L'objectif de ce travail est de fournir les ingrédients nécessaires pour prédire avec une grande précision la production d'un Higgs pseudoscalaire associé à un jet ($A + \text{jet}$) au-delà de l'ordre NNLO (Next-to-Next-to-Leading Order). Pour atteindre cet ordre de précision, il est nécessaire de calculer les amplitudes virtuelles à deux boucles pour la désintégration $A \to 3$ partons ($A \to ggg$ et $A \to q\bar{q}g$), étendues jusqu'à l'ordre $\epsilon^2$ dans le régulateur dimensionnel ($d = 4 + 2\epsilon$).

2. Méthodologie et Cadre Théorique

Cadre Effectif :
Les auteurs travaillent dans le cadre d'une théorie effective où la masse des quarks lourds est supposée infinie. Le lagrangien effectif fait intervenir deux opérateurs :

• $O_G(x)$ : Opérateur couplé au tenseur de champ de gluons.
• $O_J(x)$ : Opérateur chiral contenant $\gamma_5$ et le tenseur de Levi-Civita $\varepsilon_{\mu\nu\rho\sigma}$.
  Les coefficients de Wilson $C_G$ et $C_J$ sont développés en série de la constante de couplage forte $a_s$.

Traitement de $\gamma_5$ et Renormalisation :
La présence de $\gamma_5$ en dimensions $d \neq 4$ pose un problème de définition. Les auteurs adoptent la définition de 't Hooft et Veltman pour $\gamma_5$. Pour respecter les identités de Ward, une renormalisation spécifique du courant singulet axial est effectuée via une constante de renormalisation finie $Z_5^s$, suivie de la renormalisation MS ($Z_{MS}^s$). Les opérateurs renormalisés sont obtenus par une matrice de renormalisation mixant $O_G$ et $O_J$.

Calcul des Amplitudes :

• Processus : $A \to ggg$ et $A \to q\bar{q}g$.
• Outils : Les manipulations de Dirac et de couleur (SU(N)) sont effectuées avec le code interne FORM.
• Intégrales : Les résultats sont exprimés en termes de coefficients rationnels et d'intégrales maîtresses (Master Integrals), calculées analytiquement jusqu'à l'ordre $\epsilon^2$ en utilisant des Polylogarithmes Généralisés (GPLs).
• Optimisation : En raison de la taille massive des polynômes intermédiaires (fichiers de plusieurs Go), les auteurs utilisent l'outil MultivariateApart pour simplifier les identités d'intégration par parties (IBP).

Régularisation et Renormalisation :

• UV : Renormalisation des constantes de couplage et des opérateurs composites.
• IR : Les singularités infrarouges (IR) sont régulées dimensionnellement. Les auteurs vérifient que les pôles IR de leurs amplitudes correspondent aux prédictions universelles de factorisation connues pour les amplitudes à 2 boucles.

3. Résultats Clés

Calculs Analytiques et Numériques :

• Les auteurs ont calculé pour la première fois les éléments de matrice à deux boucles pour la désintégration $A \to ggg$ et $A \to q\bar{q}g$, étendus jusqu'à l'ordre $\epsilon^2$.
• Une implémentation numérique rapide a été développée en Fortran-95 (dérivée de routines Mathematica) pour être utilisée avec des générateurs de phase Monte Carlo.
• Comparaison : Les résultats finaux (parties finies) ont été comparés avec la référence [26] (qui couvre l'ordre NNLO complet pour le Higgs scalaire, adapté ici) et montrent un accord parfait pour les points de phase testés (voir Tableau 1 de l'article).

Performances et Observations :

• Temps de calcul : Après optimisation, le temps de compilation des fichiers Fortran à deux boucles est d'environ 10 minutes. L'évaluation par point de phase est très rapide :
  • Ordre $\epsilon^0$ (NLO) : ~100 $\mu$s.
  • Ordre $\epsilon^2$ (NNLO) : ~19 s pour les points de phase considérés.
• Dominance des canaux : L'ordre $\epsilon^2$ pour le canal $ggg$ est environ 100 fois plus grand que pour le canal $q\bar{q}g$, soulignant l'importance cruciale du calcul du canal purement gluonique pour les prédictions de haute précision.

4. Contributions Majeures

1. Première réalisation analytique : Fourniture des premiers éléments de matrice à deux boucles, étendus à l'ordre $\epsilon^2$, pour la désintégration d'un Higgs pseudoscalaire en trois partons ($ggg$ et $q\bar{q}g$).
2. Préparation pour le N$^3$LO : Ces résultats constituent un ingrédient indispensable pour calculer les corrections de troisième ordre (N$^3$LO) à la production de Higgs + 1 jet via la fusion de gluons, en utilisant les relations de symétrie de croisement (crossing symmetry).
3. Validation et Code : Validation des résultats par comparaison avec des travaux existants et mise à disposition de codes numériques optimisés pour les études phénoménologiques.

5. Signification et Impact

Ce travail est une étape fondamentale pour la physique de précision au LHC et les futurs collisionneurs.

• Précision Théorique : Il permet de réduire les incertitudes théoriques sur les sections efficaces de production du Higgs pseudoscalaire, essentielles pour distinguer les modèles BSM du Modèle Standard.
• Détection de CP : En améliorant la précision des prédictions pour $A + \text{jet}$, les expériences peuvent mieux contraindre la nature CP (impair ou mélange) du boson de Higgs découvert.
• Méthodologie : La démonstration de la faisabilité du calcul de ces termes $\epsilon^2$ complexes pour des opérateurs chiraux ouvre la voie à des calculs similaires pour d'autres processus BSM complexes.

En résumé, cet article fournit les briques théoriques et numériques nécessaires pour pousser les prédictions de la production de Higgs pseudoscalaire au-delà du niveau NNLO, un saut qualitatif nécessaire pour l'ère de haute luminosité du LHC.