NNLOJET: a parton-level event generator for jet cross sections at NNLO QCD accuracy

Cet article présente le générateur d'événements au niveau des partons NNLOJET, open source, qui implémente la méthode de soustraction par antennes pour calculer les sections efficaces de jets et les observables connexes avec une précision QCD NNLO dans les collisions électron-positron, lepton-hadron et hadron-hadron.

L'essentiel

La Vue d'Ensemble : Construire une Meilleure Boule de Cristal

Imaginez que vous êtes un physicien essayant de prédire ce qui se passe lorsque deux minuscules particules entrent en collision à une vitesse proche de celle de la lumière. Vous disposez d'un ensemble de règles (le Modèle Standard) qui vous indique comment ces particules devraient se comporter. Mais la nature est désordonnée. Lorsque les particules entrent en collision, elles ne rebondissent pas simplement ; elles projettent une douche chaotique de nouvelles particules, comme des confettis explosant depuis un canon.

Pour faire une prédiction vraiment précise, vous ne pouvez pas simplement deviner l'explosion principale. Vous devez calculer les minuscules ondulations, presque invisibles, et les pulvérisations secondaires qui se produisent à cause de l'explosion. Dans le monde de la physique des particules, ces calculs sont appelés NNLO (Next-to-Next-to-Leading Order, ou ordre suivant-le-suivant-le-principal). Ils sont incroyablement complexes, comme essayer de prédire la trajectoire exacte de chaque grain de sable dans un ouragan.

Cet article présente NNLOJET, un nouveau programme informatique open-source conçu pour effectuer ces calculs ultra-précis pour les « jets » (pulvérisations de particules). C'est comme offrir aux scientifiques une boule de cristal haute définition capable de voir des détails que les anciens modèles manquaient.

Le Problème Central : Le Désordre « Infini »

Lorsque les physiciens tentent de calculer ces collisions, ils se heurtent à un cauchemar mathématique appelé « singularités infrarouges ».

• L'Analogie : Imaginez essayer de compter le nombre de personnes dans une pièce, mais chaque fois que quelqu'un chuchote, la pièce devient infiniment plus bruyante. Si vous essayez d'additionner le bruit, le nombre tend vers l'infini.
• La Réalité : En physique des particules, lorsque les particules se déplacent très lentement ou très près les unes des autres, les mathématiques explosent vers l'infini. Puisque le monde réel n'a pas d'énergie infinie, ces infinis ne sont que des erreurs dans la méthode de calcul.

Pour résoudre ce problème, les auteurs utilisent une méthode appelée Soustraction d'Antenne.

• L'Analogie : Pensez-y comme à des écouteurs à réduction de bruit active. Les écouteurs écoutent le bruit de fond gênant (les infinis mathématiques) et génèrent une onde sonore « négative » pour l'annuler parfaitement.
• Comment cela fonctionne : Le programme calcule les parties désordonnées et infinies séparément, puis les soustrait, ne laissant que la réponse propre, finie et réelle.

Qu'est-ce que NNLOJET ?

NNLOJET est le « moteur » qui exécute ces calculs. Avant cet article, beaucoup de ces calculs étaient comme des recettes secrètes enfermées dans la cuisine d'un seul chef. Si vous vouliez les utiliser, vous deviez supplier le chef ou construire votre propre cuisine à partir de zéro.

NNLOJET change la donne en étant open-source.

• L'Analogie : C'est comme publier le livre de recettes complet et les outils de cuisine au public. N'importe qui avec un ordinateur peut le télécharger, cuisiner le plat et vérifier les ingrédients.
• La Partie « Jet » : Le programme est spécialisé dans les « jets ». En physique des particules, lorsqu'un quark ou un gluon (les blocs de construction de la matière) est éjecté, il ne reste pas seul. Il se transforme instantanément en une pulvérisation de particules. Nous appelons cette pulvérisation un « jet ». NNLOJET prédit exactement la taille, la vitesse et la direction de ces pulvérisations.

Comment l'Utiliser ? (La Runcard)

Vous n'avez pas besoin d'être un wizard du code pour utiliser NNLOJET. L'article explique que vous contrôlez le programme à l'aide d'une Runcard.

• L'Analogie : Considérez la Runcard comme un plan de vol ou une fiche de recette. Vous n'avez pas besoin de savoir comment fonctionne le moteur de l'avion ; vous devez simplement dire au pilote où aller.
• Ce que vous écrivez : Vous indiquez au programme :
  • Où : Sommes-nous en train d'écraser des protons (comme au Grand Collisionneur de Hadrons) ou des électrons ?
  • Quoi : Cherchons-nous un boson Z, un boson de Higgs, ou simplement une pulvérisation de jets ?
  • Règles : Quelle taille doivent avoir les jets ? Quelle vitesse ?
  • Sortie : Quels graphiques voulez-vous voir à la fin ?

L'article fournit un manuel détaillé (Sections 5 et 6) sur la manière d'écrire ce « plan de vol », couvrant tout, du type de collision aux paramètres mathématiques spécifiques.

Le Flux de Travail : La Chaîne de Montage

Calculer ces collisions est si lourd qu'un seul ordinateur ne peut pas le faire seul. Cela prendrait des années.

• L'Analogie : Imaginez que vous devez peindre une fresque massive. Une personne mettrait une vie entière. Au lieu de cela, vous engagez une équipe de 100 peintres.
• L'Outil : L'article décrit un script de flux de travail (nnlojet-run) qui agit comme le contremaître.
  1. Échauffement : Le contremaître envoie quelques peintres tester le mur et déterminer la meilleure façon de peindre (c'est l'adaptation de la « grille VEGAS »).
  2. Production : Une fois le plan établi, le contremaître envoie des centaines de peintres pour effectuer le travail réel simultanément.
  3. Finalisation : Lorsque tout le monde a terminé, le contremaître rassemble toutes les pièces, vérifie les erreurs (valeurs aberrantes) et les assemble en une seule image parfaite.

Que Peut-Il Calculer ?

L'article liste les « recettes » spécifiques actuellement disponibles dans le livre de recettes (Version 1.0.0) :

• Collisions Électron-Positron : Comme l'ancien collisionneur LEP.
• Collisions Proton-Proton : Comme le Grand Collisionneur de Hadrons (LHC) actuel.
• Événements Spécifiques : Il peut prédire la création de :
  • Bosons Z (cousins lourds du photon).
  • Bosons W (particules qui transportent la force faible).
  • Bosons de Higgs (la particule qui donne sa masse aux autres).
  • Photons (lumière) et jets.
  • Combinaisons de tout ce qui précède (par exemple, un boson de Higgs plus un jet).

Pourquoi Cela Compte-t-il ?

L'article conclut qu'en rendant ce code ouvert et convivial, il permet aux scientifiques de comparer directement leurs prédictions théoriques avec les données expérimentales réelles.

• L'Analogie : Si vous construisez une maquette de voiture, vous devez savoir si elle roulera réellement. NNLOJET fournit le test de « soufflerie » le plus précis possible. Si les données réelles correspondent à la prédiction NNLOJET, notre compréhension de l'univers est correcte. Si elles ne correspondent pas, cela pourrait signifier que nous avons trouvé quelque chose de nouveau !

En résumé : Cet article publie un outil puissant, gratuit et convivial qui aide les physiciens à calculer les détails désordonnés des collisions de particules avec une précision extrême, garantissant que notre carte de l'univers est aussi précise que possible.

Résumé technique

Résumé technique de NNLOJET : un générateur d'événements au niveau des partons pour les sections efficaces de jets avec une précision QCD NNLO

Énoncé du problème
Les tests de précision du Modèle Standard (MS) reposent fortement sur la comparaison entre les mesures expérimentales des états finaux hadroniques (jets) et les prédictions théoriques. Ces prédictions nécessitent le calcul des corrections de Chromodynamique Quantique (QCD) à l'ordre suivant-le-prochain-à-le-prochain (NNLO) pour gérer la constante de couplage fort ($\alpha_s$) élevée et la structure complexe des singularités infrarouges (IR) résultant du rayonnement mou et collinéaire. Bien que les calculs à l'ordre suivant-le-prochain (NLO) aient été largement automatisés, les calculs NNLO restent difficiles en raison de la nécessité d'annuler les divergences IR entre les contributions doublement virtuelles (VV), mixtes réel-virtuel (RV) et doublement réelles (RR). Les codes NNLO existants sont souvent spécifiques à un processus, pas largement disponibles en tant que code source ouvert, ou nécessitent un ajustement manuel significatif pour chaque observable. Il existe un besoin d'un générateur d'événements au niveau des partons, flexible et open-source, capable de calculer les sections efficaces de jets avec une précision NNLO dans divers environnements de collisionneurs (électron–positron, lepton–hadron et hadron–hadron).

Méthodologie : soustraction par antennes
Le code NNLOJET implémente la méthode de soustraction par antennes pour gérer les singularités IR. Cette méthode dérive des termes de soustraction basés sur la factorisation universelle des amplitudes de rayonnement réel de la QCD dans leurs limites IR.

• Mécanisme central : La méthode identifie les configurations de rayonnement réel non résolues avec des paires de partons radiateurs durs. Le rayonnement entre eux est décrit par des « fonctions d'antenne ».
• Termes de soustraction : Des termes de soustraction non intégrés sont construits à partir de produits de fonctions d'antenne et d'éléments de matrice réduits de multiplicité partonique inférieure. Ces termes sont intégrés sur l'espace de phase de l'antenne pour produire des termes de soustraction intégrés ($d\hat{\sigma}_T$ et $d\hat{\sigma}_U$) qui annulent explicitement les pôles en $\epsilon$ des corrections virtuelles.
• Implémentation : Le code utilise une cartographie de l'espace de phase pour factoriser la cinématique des fonctions d'antenne des éléments de matrice réduits. L'implémentation couvre les configurations de radiateurs final–final, initial–final et initial–initial.
• Cadre de calcul : NNLOJET est un générateur d'événements au niveau des partons écrit principalement en Fortran moderne, avec des modules C++ et Python. Il utilise la routine Monte Carlo adaptative VEGAS pour l'intégration sur l'espace de phase. Le code inclut des bibliothèques pour l'évaluation de fonctions spéciales (polynômes logarithmiques harmoniques) via HPLOG, TDHPL et CHAPLIN, et s'interface avec LHAPDF pour les fonctions de distribution des partons (PDF).

Contributions et fonctionnalités clés
L'article présente la version open-source (v1.0.0) du code NNLOJET, distribuée sous la licence publique générale GNU v3.0. Les contributions techniques clés incluent :

1. Couverture large des processus : Le code prend en charge les calculs NNLO pour la production de jets dans :

   • Collisions électron–positron : $e^+e^- \to 2j$ et $3j$.
   • Collisions lepton–hadron : Processus DIS à courant neutre (NC) et à courant chargé (CC) ($ep \to \ell + 2j$, $ep \to \nu + 2j$, etc.).
   • Collisions hadron–hadron : Production inclusive de jets ($pp \to j$), production de dijets ($pp \to 2j$), boson vectoriel plus un jet ($V+1j$), photon plus un jet ($\gamma+1j$), di-photon ($\gamma\gamma$), et production de boson de Higgs plus un jet ($H+1j$).
   • Ordres dérivés : Le cadre fournit implicitement des corrections NLO pour les processus avec un jet supplémentaire et LO pour deux jets supplémentaires (par exemple, calculer $Z+2j$ à NLO via le processus $ZJ$ avec une exigence de 2 jets).

2. Interface utilisateur flexible (Runcard) : Le code est contrôlé via un fichier « runcard » ASCII brut, permettant aux utilisateurs de définir :

   • Processus et collisionneur : Sélection du type de collisionneur (pp, ep, $e^+e^-$), des énergies de faisceau et des algorithmes de jets (anti-$k_T$, Cambridge-Aachen, $k_T$, Durham, Jade).
   • Observables et sélections : Définition des coupes de fidélité, des sélecteurs au niveau des objets ($p_T$ du jet, rapidité) et des sélecteurs au niveau des événements.
   • Histogrammes : Granularité flexible (linéaire, logarithmique, bords personnalisés) et distributions multi-différentielles via des sélecteurs d'histogrammes.
   • Échelles : Évaluation simultanée de multiples combinaisons d'échelles de renormalisation et de factorisation pour l'estimation des incertitudes.

3. Flux de travail automatisé : La distribution inclut un script de flux de travail basé sur Python (nnlojet-run) pour gérer la soumission des tâches, l'allocation des ressources et la fusion des résultats. Ce flux de travail gère :

   • Phases de préchauffage et de production : Génération adaptative de la grille suivie de runs de production.
   • Gestion des ressources : Adaptation dynamique des ressources de calcul pour atteindre une précision relative cible.
   • Fusion des résultats : Un algorithme robuste (impliquant l'élagage, le balayage k et la combinaison pondérée) pour fusionner les résultats partiels de tâches indépendantes tout en atténuant l'impact des valeurs aberrantes.

4. Validation : L'implémentation inclut des vérifications de validation rigoureuses :

   • Tests ponctuels : Vérification que le rapport des éléments de matrice aux termes de soustraction tend vers l'unité près des limites singulières.
   • Vérifications des coupes techniques : Assurer l'indépendance des résultats par rapport à la coupure infrarouge technique.
   • Tests de pôles : Annulation explicite des pôles en $\epsilon$ dans les contributions RV et VV.

Résultats et validation
L'article détaille l'implémentation de divers processus, en renvoyant à des publications spécifiques pour les détails de validation. À titre d'exemple concret, les auteurs présentent un calcul pour la production $Z + 1j$ au LHC ($\sqrt{s} = 7$ TeV) :

• Configuration : Algorithme de jets anti-$k_T$ ($R=0.4$), PDFs NNPDF3.1 NNLO, et coupes de fidélité spécifiques sur les leptons et les jets.
• Résultat : La section efficace de fidélité NNLO totale est calculée à 68,5(2) pb.
• Distributions différentielles : L'article affiche la distribution de l'impulsion transverse du jet principal, montrant la convergence du LO au NLO puis au NNLO, avec des bandes d'incertitude représentant une variation d'échelle à 7 points.

Importance et affirmations
Les auteurs affirment que NNLOJET fournit un cadre universellement applicable et open-source pour les calculs QCD NNLO des observables de jets, comblant le vide laissé par les codes précédents qui étaient souvent propriétaires ou limités à des classes de processus spécifiques.

• Accessibilité : En distribuant le code sous GPLv3.0 et en fournissant une syntaxe de runcard conviviale et un script de flux de travail, les auteurs visent à rendre les prédictions NNLO de haute précision accessibles à une communauté plus large de phénoménologistes et d'expérimentateurs.
• Flexibilité : La capacité de définir des coupes de fidélité et des observables arbitraires sans recompiler le code permet une comparaison directe avec les données expérimentales.
• Extensibilité : La conception modulaire permet l'ajout de nouveaux processus et fonctionnalités dans les futures versions.

L'article conclut que NNLOJET est un outil mature pour les études de physique de précision, le test de stress du Modèle Standard, et le raffinement de la détermination des paramètres du MS et des fonctions de distribution des partons. Le code est disponible sur HEPForge, et les auteurs invitent aux retours des utilisateurs et aux rapports de bugs.