Higher-order hadronic vacuum polarization contribution to the muon $g-2$ from lattice QCD

Cet article présente le premier calcul en QCD sur réseau de la contribution d'ordre suivant au plus bas de l'hadronisation du vide à l'anomalie magnétique du muon avec une précision inférieure au pourcent, aboutissant à un résultat de $-101.69(58)\times10^{-11}$ qui est deux fois plus précis que l'estimation du Livre Blanc de 2025 et en tension de $4.8\sigma$ avec les évaluations basées sur les données expérimentales excluant le résultat récent de CMD-3.

L'essentiel

🧲 Le Mystère du "Petit Aimant" de l'Univers

Imaginez que vous essayez de comprendre pourquoi une balle de tennis rebondit d'une certaine manière. Vous avez une théorie très précise (le Modèle Standard de la physique) qui prédit exactement comment elle devrait rebondir. Mais quand vous la lancez réellement, elle rebondit un tout petit peu différemment.

En physique des particules, cette "balle" est le muon (une sorte d'électron très lourd et instable). Sa propriété de rebondir s'appelle son moment magnétique (ou "g-2").

• Le problème : Les mesures expérimentales (ce qu'on voit dans les labos) ne correspondent pas parfaitement aux calculs théoriques.
• L'enjeu : Si l'écart est réel, cela signifie qu'il existe une nouvelle physique cachée, quelque chose que nous ne connaissons pas encore !

🍰 La Recette du Calcul : Le Gâteau Hadronique

Pour calculer la théorie, les physiciens doivent additionner plusieurs ingrédients. Le plus gros morceau du gâteau s'appelle la polarisation du vide hadronique (HVP). C'est l'effet des particules qui apparaissent et disparaissent constamment dans le vide quantique.

Jusqu'à présent, on calculait la partie principale (le "premier ordre") de deux façons :

1. La méthode "Recette de Cuisine" (Données expérimentales) : On regarde ce qui se passe dans les accélérateurs de particules et on utilise ces mesures pour deviner le résultat.
2. La méthode "Cuisine Pure" (Lattice QCD) : On utilise des supercalculateurs pour simuler l'univers depuis les règles de base, sans utiliser les données expérimentales.

Le hic : Ces deux méthodes donnaient des résultats différents pour la partie principale, créant une confusion.

🚀 La Nouvelle Découverte de ce Papier

Cet article présente une avancée majeure : les chercheurs ont calculé la partie suivante (le "second ordre" ou NLO) de ce gâteau, avec une précision incroyable (moins de 1 % d'erreur).

Voici comment ils ont fait, avec des analogies simples :

1. Le Simulateur de l'Univers (La Grille)

Imaginez que vous voulez simuler l'écoulement de l'eau dans une rivière. Vous ne pouvez pas le faire en continu, alors vous posez une grille (une moustiquaire) sur la rivière et vous regardez comment l'eau bouge à chaque carré de la grille.

• Ce que font les auteurs : Ils ont utilisé 35 grilles différentes (des ensembles de données) avec des mailles de tailles variées (de très fines à plus grossières).
• L'astuce : En regardant comment le résultat change quand on affine la grille, ils peuvent deviner ce qui se passerait si la grille était infiniment fine (la réalité parfaite). C'est ce qu'on appelle l'extrapolation vers le "continu".

2. Le Filtre Temporel (Les Fenêtres)

Le calcul est très difficile car il y a deux types de problèmes :

• Le bruit de fond lointain (Fenêtre Longue Distance) : C'est comme essayer d'entendre un chuchotement dans une tempête. Les signaux sont faibles et noyés dans le bruit.
• Les détails trop proches (Fenêtre Courte Distance) : C'est comme essayer de voir les détails d'une image trop floue parce qu'on est trop près de l'objectif.

La solution magique : Les chercheurs ont divisé le calcul en trois "fenêtres" (Court, Moyen, Long).

• La grande révélation : Ils ont découvert que dans la partie "Longue Distance", deux effets opposés (l'un négatif, l'autre positif) s'annulent presque parfaitement, comme deux personnes poussant une voiture dans des directions opposées avec la même force.
• Le résultat : Ce qui restait était très petit, mais très stable. Cela a permis de réduire énormément le "bruit" et d'obtenir un résultat très précis.

3. La Correction de la "Balance" (Isospin)

Dans leur simulation, les physiciens ont d'abord supposé que les particules étaient parfaitement symétriques (comme si les jumeaux étaient identiques). Mais dans la vraie vie, ils sont légèrement différents (comme un jumeau qui mange plus que l'autre).

• Ils ont dû ajouter une petite correction mathématique pour tenir compte de cette différence réelle entre les particules (proton/neutron, up/down). C'est comme peser un gâteau avec une balance qui n'est pas tout à fait à zéro et ajuster le poids final.

📉 Le Résultat : Qui a raison ?

Après tous ces calculs complexes, voici ce qu'ils ont trouvé :

• Leur résultat : Ils obtiennent une valeur précise pour cette partie du calcul.
• La comparaison :
  • Leur résultat est très proche de celui obtenu par la méthode "Recette de Cuisine" (les données expérimentales) si on ignore une mesure récente et controversée (CMD-3).
  • Cependant, il y a un écart de 4,8 fois l'incertitude avec les anciennes données. C'est énorme en physique ! Cela signifie que si leur calcul est juste, il y a un vrai problème avec les anciennes données expérimentales.

🌟 Pourquoi est-ce important ?

1. Indépendance : C'est la première fois qu'on calcule cette partie avec une précision aussi fine sans utiliser les données expérimentales. C'est une vérification indépendante, comme si deux juges différents arrivaient au même verdict sans se parler.
2. La précision : Leur erreur est si petite (0,6 %) qu'elle ne limite plus la précision de la théorie. Maintenant, le problème principal se déplace vers la partie principale du gâteau (le premier ordre).
3. La confiance : Le fait que leur méthode "pure" (Lattice QCD) donne un résultat cohérent avec certaines données expérimentales, mais différent d'autres, aide à trancher le débat : il faut probablement revoir certaines mesures expérimentales.

En résumé

Imaginez que vous essayez de résoudre un puzzle géant. Jusqu'ici, deux équipes s'arrachaient les cheveux car leurs pièces ne s'assemblaient pas.
Cette équipe a construit une nouvelle pièce, parfaite et ultra-précise, en utilisant uniquement les règles de base du puzzle (sans regarder l'image finale).
Leur pièce s'assemble bien avec une partie de l'image, mais pas avec l'autre. Cela suggère que l'autre partie de l'image (les anciennes données) est peut-être mal coloriée.

C'est une étape cruciale pour comprendre si l'univers contient des secrets cachés (nouvelle physique) ou si c'est juste une question de mauvais calculs sur les données existantes.

Résumé technique

Voici un résumé technique détaillé de l'article MITP-26-006 / CERN-TH-2026-035, intitulé « Higher-order hadronic vacuum polarization contribution to the muon g−2 from lattice QCD ».

1. Le Problème et le Contexte

Le moment magnétique anomal du muon, $a_\mu$, est l'un des tests de précision les plus rigoureux du Modèle Standard (SM). La mesure expérimentale finale de l'expérience E989 au Fermilab a atteint une précision de 124 parties par milliard (ppb). Cependant, la prédiction théorique du SM reste incertaine, principalement en raison de la contribution de la polarisation du vide hadronique (HVP) au premier ordre (LO), $a_\mu^{\text{hvp, LO}}$.

Il existe une tension persistante entre les évaluations de la HVP basées sur les données expérimentales (approche dispersive) et les calculs récents en Chromodynamique Quantique sur réseau (Lattice QCD). Pour assurer la cohérence, il est crucial de calculer non seulement le terme LO, mais aussi les contributions d'ordre supérieur, en particulier la contribution au prochain ordre dominant (NLO) d'ordre $\alpha^3$. Jusqu'à présent, les estimations du NLO reposaient sur des méthodes dispersives. Cet article présente la première calcul en Lattice QCD de cette contribution NLO avec une précision inférieure à 1 %.

2. Méthodologie

L'approche repose sur le calcul non perturbatif des diagrammes de Feynman contribuant à la HVP jusqu'à l'ordre $\alpha^3$.

Représentation Temps-Impulsion (TMR)

Les auteurs utilisent la représentation temps-impulsion (Time-Momentum Representation - TMR) pour exprimer les contributions HVP comme des intégrales de convolution entre des noyaux temporels ($\tilde{f}$) et la fonction de corrélation vectorielle à impulsion nulle ($G(t)$).
$$a_\mu^{\text{hvp}, (i)} = \left(\frac{\alpha}{\pi}\right)^{2+N_i} \int dt_1 \dots dt_{m_i} \, \tilde{f}^{(i)}(\hat{t}_1, \dots) \, G(t_1) \dots$$
où $\hat{t} = m_\mu t$.

Diagrammes Calculés

Le calcul couvre trois ensembles de diagrammes NLO :

1. NLOa : Corrections provenant des lignes de photons et des boucles de muons (14 diagrammes). C'est la contribution négative dominante.
2. NLOb : Corrections provenant des boucles de leptons légers (électrons et taus). La boucle électronique apporte une correction positive importante qui compense partiellement NLOa.
3. NLOc : Un seul diagramme avec deux insertions QCD (sous-dominant).

Les auteurs ont également dérivé des solutions analytiques pour les noyaux temporels de NLOc (Annexe A) et des expansions polynomiales pour NLOa et NLOb.

Configurations de Réseau et Extrapolation

• Ensembles : Utilisation de 35 ensembles de jauge générés par l'effort CLS (Coordinated Lattice Simulations) avec des fermions Wilson améliorés à l'ordre $O(a)$ et $N_f = 2+1$ saveurs.
• Paramètres : Six espacements de réseau ($a$) allant de 0,039 fm à 0,097 fm, couvrant une gamme de masses de pions incluant la valeur physique.
• Extrapolation : Une extrapolation globale chiral-continuum est effectuée vers le point physique en utilisant un moyennage de modèles (AIC) pour gérer les incertitudes systématiques liées à la discrétisation et à la masse des quarks.

Fenêtres Temporelles (Window Observables)

Pour contrôler les incertitudes systématiques, l'intégrale est divisée en trois fenêtres temporelles :

• Courte distance (SD) : Dominée par les effets de grille et la QCD perturbative. Une soustraction du comportement à courte distance est appliquée pour éliminer les effets logarithmiques renforcés ($O(a^2 \ln a)$).
• Distance intermédiaire (ID) : Zone de transition.
• Longue distance (LD) : Dominée par le bruit statistique et les effets de volume fini.
• Stratégie clé : Les auteurs exploitent la forte annulation entre les contributions NLOa et NLOb à grande distance. En calculant la somme combinée NLOa&b, le bruit statistique et la sensibilité aux effets de volume fini et au réglage de l'échelle sont considérablement réduits par rapport au calcul séparé.

Corrections

• Effets de volume fini (FV) : Corrigés via la méthode Hansen-Patella et Meyer-Lellouch-Lüscher (HP&MLL).
• Brisure d'isospin (IB) : Estimée via des modèles phénoménologiques (VMD et dépendance de masse des quarks) car les noyaux CCS pour le NLO ne sont pas disponibles.
• Quarks lourds : Contributions des quarks bottom (calculées perturbativement) et effets de désaccord de la mer de quarks charmés (estimés).

3. Résultats Principaux

Le résultat final, extrapolé à la limite du continu et corrigé pour les effets de volume et de brisure d'isospin, est :

$$a_\mu^{\text{hvp, nlo}} = -101.69(25)_{\text{stat}}(53)_{\text{syst}} \times 10^{-11}$$

• Précision : L'erreur relative totale est de 0,6 %, ce qui est comparable à la précision des évaluations dispersives basées sur les données.
• Décomposition :
  • NLOa : $-217.19 \times 10^{-11}$
  • NLOb : $+111.76 \times 10^{-11}$
  • NLOc : $+3.76 \times 10^{-11}$
  • La somme NLOa+NLOb montre une annulation significative, réduisant l'incertitude totale.

4. Comparaison et Signification

• Comparaison avec le "White Paper 2025" : Le résultat est inférieur de 1,5 $\sigma$ à l'estimation du Muon g-2 Theory Initiative (WP25) basée sur les données, mais avec une précision deux fois meilleure.
• Tension avec les données historiques : Il existe une tension de 4,8 $\sigma$ avec l'évaluation dispersive de Keshavarzi et al. (2019) qui exclut les résultats récents de CMD-3. Cette tension est cohérente avec le schéma de déviation observé pour la contribution LO.
• Signification Physique :
  1. Indépendance des données : Ce résultat fournit une prédiction du SM indépendante des données de section efficace hadronique, contournant les incohérences entre les expériences (comme la tension entre BaBar/KLOE et CMD-3).
  2. Validation de la méthode : La précision atteinte (0,6 %) démontre que le Lattice QCD est capable de calculer des termes d'ordre supérieur avec une fiabilité suffisante pour être utilisé dans la prédiction finale du SM.
  3. Cohérence LO-NLO : En fournissant LO et NLO avec la même méthode (réseau), les auteurs éliminent les biais méthodologiques potentiels dans la comparaison théorie-expérience.

Conclusion

Cet article marque une étape majeure dans la résolution de l'énigme du moment magnétique anomal du muon. En fournissant une détermination précise et ab initio de la contribution NLO de la polarisation du vide hadronique, les auteurs renforcent la crédibilité des calculs en Lattice QCD. Le résultat suggère que si la tension expérimentale persiste, elle ne peut pas être résolue par des incertitudes sur les termes d'ordre supérieur, pointant potentiellement vers une nouvelle physique au-delà du Modèle Standard, ou nécessitant une réévaluation approfondie des données expérimentales hadroniques.