First look at the evaluation of two-loop Feynman integrals for radiative return processes

Résumé technique : Premier aperçu de l'évaluation des intégrales de Feynman à deux boucles pour les processus de retour radiatif

Énoncé du problème
Les tests de précision du Modèle Standard dans les collisionneurs électron-positron à basse énergie (par exemple, BaBar, Belle II, BESIII, KLOE) reposent sur des prédictions théoriques précises pour les processus de retour radiatif, tels que $e^+e^- \to \pi^+\pi^-\gamma$ et $e^+e^- \to \mu^+\mu^-\gamma$. Ces processus sont critiques pour déterminer la contribution de la polarisation du vide hadronique au moment magnétique anomal de la particule muon ($(g-2)_\mu$). Pour résoudre les tensions persistantes entre les mesures expérimentales et les prédictions théoriques, des corrections QED au niveau NNLO (Next-to-Next-to-Leading Order) sont nécessaires.

Le principal goulot d'étranglement pour atteindre la précision NNLO pour ces processus, tout en conservant la dépendance complète vis-à-vis de la masse de l'électron, est l'évaluation des intégrales de Feynman à quatre points et deux boucles. La présence de propagateurs massifs introduit une complexité algébrique significative dans les réductions par intégration par parties (IBP) et les équations différentielles (DE). De plus, ces intégrales impliquent des fonctions spéciales associées à des géométries de genre supérieur, spécifiquement des courbes elliptiques et, dans certains cas, des racines carrées imbriquées. La construction d'une base canonique pour ces intégrales — où les DE sont $\epsilon$-factorisées et exprimées en termes de formes d'un $\text{dlog}$ — est mathématiquement difficile et n'est souvent pas strictement nécessaire pour la phénoménologie numérique.

Méthodologie
Les auteurs se concentrent sur la contribution du rayonnement de l'état initial (ISR) au processus $e^+e^- \to \gamma\gamma^*$, ce qui réduit le problème à l'évaluation d'amplitudes authentiques à quatre points et deux boucles. Ils considèrent trois familles d'intégrales planaires, notées PL1, PL2 et PL3, qui contiennent un total de 212 intégrales maîtres (MI), incluant celles associées à cinq courbes elliptiques non isomorphes et un secteur avec une racine carrée imbriquée.

Au lieu de tenter de construire une base entièrement canonique pour tous les secteurs (ce qui est difficile pour les géométries elliptiques et les racines imbriquées), les auteurs adoptent une stratégie développée dans des travaux précédents [48, 49]. Leur approche consiste en :

1. Construction de la base : Ils construisent une base de MI telle que le système de DE résultant soit polynomial en $\epsilon$ par rapport au régulateur dimensionnel. Pour les secteurs polylogarithmiques, ils visent des formes canoniques. Pour les secteurs elliptiques, ils identifient les MI associés aux différentielles elliptiques de première, deuxième et troisième espèce.
2. Gestion des géométries elliptiques : Pour les secteurs où la nature elliptique n'est pas manifeste dans l'intégrale de coin (par exemple, les secteurs de type "kite" à quatre points), ils utilisent des super-secteurs réductibles et analysent les opérateurs de Picard-Fuchs pour identifier la géométrie et construire des bases appropriées.
3. Gestion des racines carrées imbriquées : Dans un secteur spécifique de PL1 impliquant une racine carrée imbriquée, ils démontrent que bien qu'une transformation existe pour $\epsilon$-factoriser les DE (impliquant un duplet de racines imbriquées), ils préfèrent travailler avec une base qui maintient les DE linéaires en $\epsilon$ et évite l'introduction de structures transcendantes complexes, en s'appuyant plutôt sur des fonctions rationnelles et des racines carrées simples.
4. Intégration numérique : Les DE sont résolues numériquement. Les conditions limites sont déterminées à l'aide du package AMFlow à un point spécifique de l'espace des phases choisi pour simuler des configurations cinématiques de basse énergie. L'intégration est réalisée à l'aide d'une implémentation interne en langage de programmation Julia, employant des intégrateurs non raides (Tsit5 et Vern).

Contributions clés et résultats

• Cadre systématique : Le papier présente la première construction systématique de DE pour les familles d'intégrales planaires pertinentes pour les processus de retour radiatif NNLO avec une dépendance complète à la masse de l'électron.
• Traitement des secteurs elliptiques : Les auteurs identifient et traitent avec succès les secteurs associés à cinq courbes elliptiques distinctes. Ils fournissent des bases explicites pour ces secteurs, y compris les topologies de type "kite" à quatre points difficiles, où le comportement elliptique provient d'intégrales avec des points (puissances de propagateurs $>1$) plutôt que des intégrales de coin.
• Analyse des racines carrées imbriquées : Le document détaille le traitement d'un secteur de PL1 contenant une racine carrée imbriquée. Il construit une base qui satisfait des DE polynomiales en $\epsilon$ sans nécessiter l'introduction explicite de la structure de la racine imbriquée dans les équations différentielles, facilitant ainsi la stabilité numérique.
• Validation numérique : Les auteurs démontrent des évaluations numériques stables dans la région physique pertinente pour les expériences de basse énergie ($s \in \{1.0404, 16, 100\} \text{ GeV}^2$).
  • Performance : Les temps d'évaluation varient de 40 à 800 ms par point de l'espace des phases.
  • Précision : Les comparaisons avec AMFlow montrent un accord d'au moins huit chiffres significatifs à des points de test spécifiques. Les comparaisons avec le package DiffExp (utilisant la double précision) montrent environ 5 chiffres d'accord pour les familles PL1 et 2 chiffres pour les familles plus complexes PL2 et PL3 dans les pires scénarios.
  • Robustesse : La méthode reste stable malgré de grandes hiérarchies d'échelles ($m_e^2/s \sim 10^{-7} - 10^{-8}$).

Signification
Le papier affirme fournir des « blocs de construction essentiels pour les prédictions NNLO des processus de retour radiatif ». En construisant et en intégrant numériquement avec succès les DE pour les secteurs planaires, les auteurs ont surmonté l'obstacle computationnel majeur pour la composante du rayonnement de l'état initial de l'amplitude NNLO.

Le travail est modeste dans sa portée, précisant explicitement qu'il ne couvre que les familles planaires. L'extension aux familles non planaires et l'assemblage de l'amplitude de diffusion complète sont notés comme des travaux en cours. La signification réside dans la démonstration qu'une approche polynomiale en $\epsilon$, combinée à l'intégration numérique, est une stratégie viable et efficace pour traiter les géométries complexes (courbes elliptiques et racines imbriquées) inhérentes aux calculs à deux boucles massifs, contournant ainsi le besoin de formes canoniques analytiques complètes qui sont actuellement difficiles à obtenir pour de telles géométries.