Fully differential Higgs boson pair production at N$^3$LO with top quark mass effects

Cet article présente les premières prédictions entièrement différentielles pour la production de paires de bosons de Higgs via la fusion de gluons-gluons à N$^3$LO dans la limite du top lourd, démontrant une réduction par trois des incertitudes d'échelle et incorporant les effets de la masse du quark top afin de fournir des références théoriques de haute précision pour les recherches au LHC.

L'essentiel

Imaginez l'univers comme une machine géante et complexe, et le boson de Higgs comme un engrenage crucial qui donne leur poids aux autres particules. Les physiciens cherchent à comprendre comment cet engrenage fonctionne en faisant s'entrechoquer des particules à des vitesses incroyables dans le Grand Collisionneur de Hadrons (LHC). Plus précisément, ils essaient de voir ce qui se passe lorsque deux bosons de Higgs sont créés simultanément. C'est comme essayer d'attraper deux papillons rares et insaisissables au milieu d'une tempête pour voir comment ils interagissent.

Ce document est une avancée majeure dans le « manuel d'instructions » (prédiction théorique) pour repérer ces paires de papillons. Voici la décomposition en termes simples :

1. Le problème : Une boucle très lourde

Pour créer deux bosons de Higgs, les particules s'entrechoquent et créent une « boucle » temporaire impliquant un quark top (la particule la plus lourde connue).

• L'analogie : Imaginez essayer de prédire la trajectoire d'une balle roulant dans un labyrinthe. Le labyrinthe est constitué du quark top. Parce que le quark top est si lourd, le labyrinthe est incroyablement complexe.
• L'ancienne méthode : Pendant des années, les scientifiques ont utilisé un raccourci appelé la « limite du top lourd ». Ils faisaient comme si le quark top était infiniment lourd, ce qui lissait le labyrinthe pour en faire un sol plat et simple. Cela rendait les calculs plus faciles, mais ce n'était pas parfaitement précis, surtout lorsque les particules se déplaçaient très vite.
• La nouvelle méthode : Ce document calcule la trajectoire à travers le véritable labyrinthe (avec le poids réel du quark top), mais seulement pour la première étape du voyage (Ordre suivant, ou NLO). Cependant, pour les parties principales et complexes du voyage, ils utilisent le raccourci du « sol lisse », mais calculé avec un niveau de détail sans précédent.

2. La percée : Calculer à l'échelle « N3LO »

Le document rapporte les premières prédictions entièrement différentielles à l'échelle N3LO (Ordre suivant, suivant, suivant le premier ordre).

• L'analogie : Pensez à calculer la météo.
  • LO (Ordre de base) : « Il pourrait pleuvoir. » (Une estimation très grossière).
  • NLO : « Il va pleuvoir dans l'après-midi. » (Mieux).
  • NNLO : « Il pleuvra à 15h avec 50 % d'humidité. » (Très bon).
  • N3LO : « Il pleuvra à 15h04, avec 50 % d'humidité, et les gouttes frapperont le sol avec un angle de 45 degrés. » (Extrêmement précis).
• Ce qu'ils ont fait : Ils ont calculé la « météo » de la collision du boson de Higgs avec ce niveau de précision extrême. Ils n'ont pas seulement calculé la quantité totale de pluie (section efficace totale) ; ils ont calculé exactement où et comment elle tombe (distributions différentielles), comme la vitesse et l'angle des bosons de Higgs.

3. Les résultats : Un focus plus net

• Réduction de l'incertitude : Avant ce document, les « prévisions » avaient une marge d'erreur importante (comme dire « il pourrait pleuvoir entre 10h et 18h »). Les nouveaux calculs N3LO réduisent considérablement cette fenêtre, diminuant l'incertitude d'environ trois fois. Désormais, la prédiction est assez précise pour atteindre le « niveau du pourcentage ».
• La forme de la tempête : Ils ont découvert que si la quantité totale de « pluie » ne changeait pas beaucoup, la forme de la tempête, elle, changeait. Les nouveaux calculs modifient la façon dont les bosons de Higgs sont distribués en termes de vitesse et de direction. Cela est crucial car si les « prévisions » (théorie) ne correspondent pas à la « météo réelle » (expérience), cela pourrait signifier qu'une nouvelle physique se cache dans les données.

4. Corriger le raccourci du « Top lourd »

Puisque le raccourci du « Top lourd » n'est pas parfait lorsque les particules se déplacent rapidement, les auteurs ont combiné leurs calculs ultra-précis du « sol lisse » avec un calcul plus exact du « véritable labyrinthe » pour la première étape.

• L'analogie : Imaginez que vous avez une carte extrêmement détaillée d'une ville (N3LO), mais que vous savez que votre carte est légèrement erronée concernant la hauteur des bâtiments. Vous prenez une photo brute, à basse résolution, des bâtiments réels (NLO avec la masse réelle) et vous utilisez cette photo pour corriger les hauteurs sur votre carte ultra-détaillée.
• Le résultat : Cette approche hybride offre l'image la plus précise de la production de paires de bosons de Higgs à ce jour. Ils ont constaté que la « masse réelle » du quark top modifie considérablement les prédictions, en particulier pour les bosons de Higgs se déplaçant à grande vitesse ou dans des directions spécifiques.

5. Pourquoi cela importe (selon le document)

Le document stipule que ce niveau de précision est essentiel pour les expériences en cours au LHC.

• L'objectif : Les scientifiques cherchent des signes indiquant que le potentiel de Higgs (le champ d'énergie qui donne leur masse aux particules) se comporte différemment de ce que prédit le Modèle Standard.
• Le besoin : Pour trouver ces infimes différences, vous avez besoin d'une « règle » (prédiction théorique) qui soit incroyablement précise. Si votre règle est floue, vous ne pouvez pas savoir si l'objet que vous mesurez est légèrement différent ou si c'est simplement votre règle qui est mauvaise. Ce document fournit une règle bien plus nette.

En résumé :
Ce document est une leçon magistrale de précision mathématique. Il prend un problème de physique notoirement difficile (la création de deux bosons de Higgs via une boucle de quark top lourd) et le calcule avec le plus haut degré de précision actuellement disponible. En affinant la « carte » de comportement de ces particules, il permet aux expérimentateurs du LHC de chercher une nouvelle physique avec une vision beaucoup plus nette, réduisant le « brouillard » de l'incertitude théorique qui obscurcissait la vue depuis des années.

Résumé technique

Résumé technique : Production de paires de bosons de Higgs différentielle et complète à l'N3LO avec effets de la masse du quark top

Problème et motivation
La production de paires de bosons de Higgs ($gg \to hh$) par fusion gluon-gluon est le canal dominant pour sonder le potentiel de Higgs et le couplage trilinéaire ($\lambda_{3h}$) au Grand Collisionneur de Hadrons (LHC). Bien que les corrections QCD au niveau NLO (Next-to-Leading Order) soient disponibles avec une dépendance complète à la masse du quark top, elles laissent des incertitudes d'échelle significatives (environ $\pm 13\%$) et introduisent de grandes incertitudes théoriques liées au schéma de renormalisation de la masse du quark top (jusqu'à 30 % pour certaines distributions). Des calculs d'ordres supérieurs dans la limite du quark top lourd (HTL, où $m_t \to \infty$) ont été réalisés jusqu'à l'ordre N3LO (Next-to-Next-to-Next-to-Leading Order) pour les sections efficaces inclusives et la distribution de la masse invariante di-Higgs. Cependant, les prédictions différentielles complètes à l'N3LO dans la limite HTL, particulièrement pour d'autres observables cinématiques, étaient auparavant indisponibles ou seulement approximatives. De plus, l'approximation HTL se dégrade lorsque l'énergie centre de masse partonique approche $2m_t$, nécessitant l'inclusion des effets de la masse finie du top pour une phénoménologie fiable.

Méthodologie et cadre de calcul
Les auteurs réalisent le premier calcul QCD différentiel complet à l'N3LO pour $gg \to hh$ dans la limite HTL, combiné à des prédictions NLO incluant la dépendance complète à la masse du quark top ($N_{LO}^{mt}$).

1. Décomposition HTL : Le calcul dans la limite HTL est organisé en utilisant un lagrangien effectif où les boucles de quarks top sont intégrées en interactions de contact locales. La amplitude au carré est décomposée en trois classes basées sur le nombre d'insertions de sommets effectifs :

   • Classe-a : Deux sommets effectifs (dominante, $O(\alpha_s^2)$ à l'ordre LO).
   • Classe-b : Trois sommets effectifs ($O(\alpha_s^3)$ à l'ordre LO).
   • Classe-c : Quatre sommets effectifs ($O(\alpha_s^4)$ à l'ordre LO).
     Pour atteindre la précision N3LO ($O(\alpha_s^5)$), les auteurs calculent les corrections N3LO pour la Classe-a, les corrections NNLO pour la Classe-b, et les corrections NLO pour la Classe-c.

2. Techniques de calcul :

   • Classe-a : Calculée à l'aide du framework NNLOJET. Les auteurs relient la section efficace di-Higgs à la section efficace du Higgs unique (connue à l'N3LO) via une cartographie cinématique. Ils emploient la méthode de $q_T$-slicing pour gérer les divergences infrarouges, séparant l'espace des phases en une région de faible $q_T$ (calculée via la ressommation $q_T$ dans la théorie de l'effet de collision douce et colinéaire ou SCET) et une région de haute $q_T$ (calculée comme $hh + \text{jet}$ à l'ordre NNLO).
   • Classe-b : Les corrections NNLO sont calculées en utilisant une combinaison de NNLOJET (pour la partie de ressommation à bas $q_T$) et MadGraph5_aMC@NLO (pour la partie de rayonnement réel à haute $q_T$), en utilisant le schéma de soustraction FKS.
   • Classe-c : Calculée entièrement à l'ordre NLO avec MadGraph5_aMC@NLO et la soustraction FKS.
   • Validation : L'implémentation est validée par rapport aux résultats inclusifs connus de iHixs2 et aux résultats différentiels de MadGraph5_aMC@NLO à des ordres inférieurs, montrant un accord parfait. L'indépendance des résultats vis-à-vis du paramètre de découpage $p_T^{\text{veto}}$ est vérifiée, confirmant l'annulation des corrections de puissance.

3. Incorporation des effets de la masse du quark top :
   Pour tenir compte des effets de la masse finie du top, les auteurs combinent les résultats N3LO HTL avec des résultats NLO complets dépendant de $m_t$ (calculés avec le générateur ggxy). Ils évaluent trois schémas de combinaison :

   • Additif ($N_{kLO} \oplus N_{LO}^{mt}$): Ajoute les corrections d'ordre supérieur HTL au résultat NLO complet de $m_t$.
   • Amélioré par le Born ($N_{kLO}^{B-i} \oplus N_{LO}^{mt}$): Repondère les corrections HTL par le ratio de la section efficace LO full-$m_t$ sur la section efficace LO HTL.
   • Multiplicatif ($N_{kLO} \otimes N_{LO}^{mt}$): Repondère les prédictions HTL par le ratio de la section efficace NLO full-$m_t$ sur la section efficace NLO HTL.
     Le papier identifie le schéma multiplicatif comme étant le plus stable perturbativement pour les distributions différentielles.

Résultats clés
Les calculs sont présentés pour des collisions proton-proton à $\sqrt{s} = 14$ TeV sous des coupures fiduciales réalistes ($p_{T,h1} > 30$ GeV, $p_{T,h2} > 20$ GeV, $|y_h| < 2.4$).

1. Sections efficaces fiduciales (HTL) :

   • Les corrections QCD N3LO augmentent la section efficace fiduciale NNLO d'environ 3,8 %.
   • Les incertitudes d'échelle sont considérablement réduites : de $\pm 18\%$ à l'ordre NLO à $\pm 5,2\%$ à l'ordre NNLO, et jusqu'à $+1,4\%/-2,9\%$ à l'ordre N3LO. Cela représente une réduction d'un facteur d'environ trois par rapport au NNLO.
   • Les contributions de la Classe-a dominent ($\sim 103\%$ du total), tandis que les Classes-b et -c sont négligeables ($\sim -3,6\%$ et $\sim 0,03\%$, respectivement).

2. Distributions différentielles (HTL) :

   • Masse invariante ($m_{hh}$) : Les corrections N3LO sont relativement plates ($K \approx 1,04$) sauf près du seuil ($m_{hh} \to 2m_h$), où les incertitudes d'échelle restent importantes.
   • Impulsion transverse ($p_T$) : Les spectres de $p_T$ des bosons de Higgs menant et sous-dominant montrent des modifications de forme significatives à l'N3LO. Dans la région de bas $p_T$, les prédictions à ordre fixe sont instables (corrections négatives, grandes incertitudes), indiquant la nécessité de la ressommation. Dans la région de haut $p_T$ ($>100$ GeV), les corrections N3LO modifient la forme des distributions NNLO.
   • Angle azimutal : La distribution $|\phi_{h1} - \phi_{h2}|/\pi$ montre une excellente convergence entre le NNLO et le N3LO, avec des corrections atteignant $\sim 10\%$ dans certains intervalles.
   • Rapidité : Les distributions de rapidité présentent des facteurs $K$ plats et les bandes N3LO se situent à l'intérieur des bandes d'incertitude NNLO.

3. Effets de la masse finie du top :

   • La combinaison du HTL N3LO avec le NLO complet de $m_t$ (via le schéma multiplicatif) révèle que les effets de la masse finie du top sont essentiels même aux masses invariantes les plus basses ($m_{hh} \approx 250$ GeV).
   • Comparé au HTL, les résultats complets de $m_t$ montrent un enrichissement d'un facteur deux dans la région du seuil et une suppression dans la queue de haute masse.
   • Les spectres de $p_T$ sont adoucis dans les queues lorsque la dépendance complète à la masse est incluse.
   • La distribution de l'angle azimutal est supprimée de plus de 10 % dans la région dos à dos et par un facteur d'environ 1,7 dans la région colinéaire.
   • Les incertitudes d'échelle pour les prédictions combinées sont généralement faibles, bien que la dépendance d'échelle de $N_{LO}^{mt}$ diffère légèrement de celle du HTL.

Signification et affirmations
Le papier affirme fournir les premières prédictions différentielles complètes à l'N3LO pour la production de paires de bosons de Higgs dans la limite HTL. En combinant celles-ci avec les effets de masse du top à l'ordre NLO, les auteurs présentent « l'une des prédictions différentielles à l'ordre partiel la plus précises à ce jour » pour les recherches au LHC.

Les principales affirmations concernant l'impact de ce travail incluent :

• Réduction de l'incertitude : Les corrections QCD N3LO réduisent les incertitudes d'échelle des prédictions fiduciales et différentielles NNLO d'un facteur d'environ trois, ramenant l'incertitude théorique au niveau du pourcentage dans la limite HTL.
• Modifications de forme : Les résultats démontrent que les corrections N3LO ne sont pas simplement un facteur de normalisation ; elles modifient significativement les formes des distributions d'impulsion transverse et d'angle azimutal, soulignant la nécessité d'inclure les corrections N3LO pour une phénoménologie de précision.
• Nécessité de la masse du top : L'étude confirme que les effets de la masse finie du quark top ne peuvent être négligés même aux masses invariantes les plus basses et doivent être inclus à l'ordre le plus élevé disponible (NLO dans ce contexte) pour éviter des distorsions significatives dans les distributions différentielles.

Les auteurs notent modestement que bien que les incertitudes d'échelle soient maîtrisées, d'autres incertitudes théoriques subsistent, notamment les corrections de masse finie du top manquantes à l'ordre NNLO et au-delà, les incertitudes de schéma de la masse du quark top (OS vs $\overline{MS}$), les corrections électrofaibles d'ordres supérieurs manquantes, et les contributions des boucles de quarks bottom. Ils suggèrent que des travaux futurs étendant la dépendance complète à la masse au NNLO pourraient encore réduire ces incertitudes.