Lattice determination of the higher-order hadronic vacuum polarization contribution to the muon $g-2$

Diese Studie präsentiert die erste Gitter-QCD-Berechnung des hadronischen Vakuumpolarisationsbeitrags nächster Ordnung zur Anomalie des Myon-g-2 mit einer Präzision unter 1 %, die eine signifikante Spannung zu datengetriebenen Auswertungen aufweist.

Das Wesentliche

Das große Rätsel: Warum das Muon „wackelt"

Stellen Sie sich vor, Sie haben einen winzigen, unsichtbaren Kreisel, der sich im Universum dreht. Das ist ein Muon, ein Teilchen, das dem Elektron sehr ähnlich ist, aber viel schwerer. Physiker haben herausgefunden, dass dieser Kreisel nicht ganz so stabil rotiert, wie die aktuellen Gesetze der Physik (das sogenannte „Standardmodell") vorhersagen. Er wackelt ein winziges bisschen mehr.

Dieses Wackeln nennt man das anomale magnetische Moment ($g-2$). Es ist wie ein extrem empfindlicher Kompass, der uns verrät, ob es im Universum noch unbekannte Kräfte oder Teilchen gibt, die wir noch nicht entdeckt haben.

Das Problem: Ein ungelöster Streit

Um zu berechnen, wie stark dieses Wackeln sein sollte, müssen Physiker eine riesige Rechnung aufstellen. Ein Teil dieser Rechnung ist besonders schwierig: Die Wechselwirkung mit sogenannten Hadronen (Teilchen aus Quarks, wie Protonen und Neutronen).

Bisher gab es zwei Wege, dieses Problem zu lösen:

1. Der „Daten-Weg": Man misst in Teilchenbeschleunigern, wie oft bestimmte Teilchen kollidieren, und nutzt diese Messdaten, um die Rechnung zu machen.
2. Der „Rechen-Weg" (Gitter-QCD): Man baut ein virtuelles Gitter im Computer, auf dem man die Gesetze der starken Kraft direkt simuliert, ohne auf Messdaten angewiesen zu sein.

Das Problem: Die Ergebnisse dieser beiden Wege stimmen nicht überein. Es gibt einen Streit unter den Physikern. Ein neuer, sehr genauer Experiment-Datensatz (CMD-3) hat die Situation noch verworrener gemacht.

Die neue Lösung: Der „Gitter-Rechner" schlägt zu

In dieser Arbeit präsentieren die Autoren (eine Gruppe von Physikern aus Mainz, Darmstadt, CERN und Bonn) einen Meilenstein: Sie haben die zweite, noch schwierigere Stufe dieser Rechnung komplett neu berechnet – und zwar nur mit dem Computer (Gitter-QCD).

Bisher konnten sie nur den Hauptteil (die „erste Stufe") berechnen. Jetzt haben sie auch den feinen, hochkomplexen Rest (die „zweite Stufe" oder NLO) berechnet.

Die magische Zauberei: Das „Wegfallen" der Fehler

Das Geniale an ihrer Methode ist eine Art magisches Aufheben, das man sich wie folgt vorstellen kann:

Stellen Sie sich vor, Sie wollen die genaue Höhe eines Berges messen, aber es gibt einen dichten Nebel (Rauschen), der die Sicht behindert.

• Bei der ersten Berechnung (LO) war der Nebel sehr dick.
• Bei der zweiten Berechnung (NLO) haben die Forscher zwei verschiedene Arten von „Nebel" entdeckt:
  1. Ein Nebel, der die Messung nach unten zieht (NLOa).
  2. Ein Nebel, der die Messung nach oben zieht (NLOb).

Wenn man diese beiden Nebel zusammen betrachtet, heben sie sich fast perfekt gegenseitig auf! Das ist wie bei zwei Personen, die auf einer Wippe sitzen: Wenn eine nach unten drückt und die andere genau gleich stark nach oben, bleibt die Wippe ruhig.

Dadurch verschwindet der störende Nebel (das statistische Rauschen) fast vollständig. Das erlaubt den Physikern, eine unglaublich präzise Messung durchzuführen, die viel genauer ist als alles, was vorher mit Daten aus Experimenten erreicht wurde.

Das Ergebnis: Ein neuer, präziser Wert

Die Forscher haben herausgefunden, dass der korrekte Wert für diesen Teil der Rechnung -101,57 ist (in einer speziellen Einheit).

• Genauigkeit: Sie haben eine Fehlerquote von nur 0,6 %. Das ist doppelt so präzise wie die besten bisherigen Schätzungen, die auf Experimentdaten basierten.
• Der Vergleich: Ihr Ergebnis liegt etwas unter dem Durchschnittswert, den die große Fachgruppe („White Paper") bisher angenommen hat. Es stimmt zwar grob überein, aber der Unterschied ist statistisch signifikant.
• Die Spannung: Ihr Ergebnis steht in starkem Konflikt (4,6 σ) zu den älteren Daten, die den neuen CMD-3-Messungen nicht berücksichtigt haben.

Warum ist das wichtig?

Stellen Sie sich vor, Sie bauen ein Haus (das Standardmodell der Physik). Sie haben einen Riss in der Wand (die Diskrepanz beim Muon-Wackeln).

• Früher sagten die Baumeister: „Der Riss kommt von schlechten Ziegeln (Messdaten)."
• Jetzt sagen diese Forscher: „Nein, wenn wir die Rechnung direkt im Computer machen, sehen wir, dass der Riss noch tiefer geht, als wir dachten."

Ihre Arbeit ist ein unabhängiger Beweis. Sie zeigen, dass das Standardmodell wirklich eine Lücke hat, die wir noch nicht verstanden haben. Da ihre Methode (der Computer-Rechnung) völlig unabhängig von den streitenden Experiment-Daten ist, gibt sie uns eine neue, klare Perspektive.

Zusammenfassend: Diese Forscher haben einen extrem schwierigen mathatischen Knoten gelöst, indem sie zwei sich gegenseitig aufhebende Effekte nutzten. Sie haben damit die Genauigkeit der Vorhersage für das Muon-Wackeln verdoppelt und bestätigen, dass es im Universum wahrscheinlich noch etwas gibt, das wir noch nicht kennen – etwas, das unser Verständnis der Physik revolutionieren könnte.

Technische Zusammenfassung

Titel: Gitterbestimmung des Beitrags höherer Ordnung der hadronischen Vakuumpolarisation zum anomalen magnetischen Moment des Myons ($g-2$)

Autoren: Arnau Beltran et al. (CLS-Kollaboration)
Veröffentlichung: LATTICE2025 (Präsentation), arXiv:2604.08351

---

1. Problemstellung und Motivation

Das anomale magnetische Moment des Myons ($a_\mu$) ist einer der präzisesten Tests des Standardmodells (SM) der Teilchenphysik. Die aktuelle Diskrepanz zwischen Experiment (Fermilab E989) und theoretischer Vorhersage wird maßgeblich durch die Unsicherheit im führenden Ordnung (LO) Beitrag der hadronischen Vakuumpolarisation (HVP) verursacht.

• Herausforderung: Die traditionelle datengetriebene dispersive Methode für den LO-Beitrag ist anfällig für Inkonsistenzen in experimentellen Eingangsdaten (insbesondere $e^+e^- \to \pi^+\pi^-$), wie sie durch die CMD-3-Messungen aufgezeigt wurden.
• Notwendigkeit: Für eine konsistente theoretische Vorhersage müssen sowohl der LO- als auch der NLO-Beitrag (Next-to-Leading Order) innerhalb desselben Rahmens berechnet werden. Bisher basierte die NLO-Schätzung (WP25) noch vollständig auf datengetriebenen Methoden und erbte deren Unsicherheiten.
• Ziel: Die erste Gitter-QCD-Berechnung des NLO-HVP-Beitrags mit einer Präzision im Sub-Prozent-Bereich, um eine unabhängige und konsistente Überprüfung des SM zu ermöglichen.

2. Methodik und Formalismus

Die Studie verwendet die Zeit-Impuls-Darstellung (Time-Momentum Representation, TMR) in Kombination mit räumlich summierten Vektor-Korrelatoren.

• Gitterkonfigurationen: Berechnungen wurden auf 35 CLS-Ensembles (Coordinated Lattice Simulations) durchgeführt.
  • Flavours: $N_f = 2+1$ (up, down, strange) mit $O(a)$-verbesserten Wilson-Fermionen.
  • Bereich: Sechs Gitterabstände ($a \approx 0.039 - 0.097$ fm) und Pionmassen von $\sim 430$ MeV bis zum physikalischen Punkt.
  • Skalierung: Bestimmung über die Gradientenfluss-Observable $\sqrt{t_0}$.
• Diagramm-Klassen: Die Berechnung umfasst drei Klassen von NLO-Diagrammen (siehe Abb. 1 im Paper):
  1. NLOa: Photon-Linien und Myon-Schleifen (dominierender negativer Beitrag).
  2. NLOb: Elektron- und Tau-Schleifen (positiver Beitrag, teilweise Kompensation von NLOa).
  3. NLOc: Zwei QCD-Einfügungen (untergeordneter Beitrag).
• Kernfunktionen: Es werden raumartige NLO-Kernfunktionen verwendet, die über Integraltransformationen aus den zeitartigen Kernen gewonnen werden.
• Analyse-Strategie:
  • Fenster-Zerlegung: Die Integrale werden in kurze (SD), mittlere (ID) und lange (LD) Distanzen unterteilt ($t_1=0.4$ fm, $t_2=1.0$ fm).
  • Kürzungseffekt (Cancellation): Ein entscheidendes Merkmal ist die starke gegenseitige Aufhebung der NLOa- und NLOb-Kernfunktionen bei großen euklidischen Zeiten ($t \gtrsim 1.5$ fm). Dies reduziert die Empfindlichkeit gegenüber statistischem Rauschen und endlichen Volumen-Effekten drastisch im Vergleich zum LO-Fall.
  • Kontinuumslimit: Globale chirale-kontinuierliche Fits unter Verwendung des Akaike-Informationskriteriums (AIC) zur Modellmittelung.
  • Korrekturen: Berücksichtigung von endlichen Volumen-Effekten (Hansen-Patella / Meyer-Lellouch-Lüscher Methoden), Isospin-Bruch (elektromagnetisch und stark) und diskretisierungsbedingten Artefakten (insbesondere logarithmisch verstärkte $O(a^2 \ln^2 a)$ Terme im SD-Bereich).

3. Schlüsselbeiträge und Ergebnisse

Die Autoren präsentieren den ersten Gitter-QCD-Wert für den NLO-HVP-Beitrag mit einer Gesamtunsicherheit von 0,6 %.

Ergebnisse (in Einheiten von $10^{-11}$):

• NLOa: $-216.83(81)(87)(99)(81)(34)(31)[1.81]$
• NLOb: $+111.51(55)(64)(64)(58)(16)(25)[1.24]$
• NLOa & b (Summe): $-105.33(28)(28)(35)(24)(18)(5)[61]$
• NLOc: $+3.75(6)(6)(3)(3)[10]$
• Gesamtergebnis (NLO HVP):
  $$a_\mu^{\text{hvp, nlo}} = (-101.57 \pm 0.26_{\text{stat}} \pm 0.54_{\text{syst}}) \times 10^{-11}$$
  (Gesamtunsicherheit: $\pm 0.59 \times 10^{-11}$)

Wichtige Beobachtungen:

1. Präzision: Das Ergebnis ist doppelt so präzise wie die Schätzung im "White Paper 2025" (WP25).
2. Rolle der LD-Region: Durch die starke Kürzung zwischen NLOa und NLOb trägt der lange Distanzbereich (LD), der normalerweise die größte Unsicherheitsquelle darstellt, nur noch etwa ein Drittel zum Gesamtwert bei (im Gegensatz zu ca. zwei Dritteln beim LO).
3. Isospin-Bruch: Die elektromagnetischen und starken Isospin-Bruch-Korrekturen heben sich teilweise auf, was zu einer kleinen und gut kontrollierten Korrektur von $\Delta_{\text{IB}} = 0.06(27) \times 10^{-11}$ führt.

4. Vergleich und Signifikanz

• Vergleich mit WP25: Der berechnete Wert liegt bei $1.4\sigma$ unter dem WP25-Durchschnitt von $(-99.6 \pm 1.3) \times 10^{-11}$.
• Vergleich mit KNT19: Es besteht eine Spannung von $4.6\sigma$ zur datengetriebenen KNT19-Bestimmung (die keine CMD-3-Daten enthält).
• Konsistenz mit LO: Dieses Muster (Gitterwerte niedriger als datengetriebene Werte ohne CMD-3) spiegelt die bereits bekannte Diskrepanz im LO-Bereich wider.
• Bedeutung:
  • Dies ist die erste vollständig konsistente Gitterbestimmung sowohl für LO als auch für NLO HVP.
  • Es eliminiert die methodische Inkonsistenz, bei der bisher nur der LO-Teil aus Gitterrechnungen stammte.
  • Die hohe Präzision zeigt, dass Gitter-QCD in der Lage ist, auch höhere Ordnungen der HVP mit Sub-Prozent-Genauigkeit zu bestimmen, was für die zukünftige Interpretation der $g-2$-Messungen entscheidend ist.

Fazit

Die Arbeit demonstriert, dass durch die systematische Ausnutzung der strukturellen Kürzungen zwischen den NLOa- und NLOb-Diagrammen eine signifikant höhere Präzision erreicht werden kann als beim LO-Beitrag. Das Ergebnis untermauert die Notwendigkeit, die datengetriebenen Eingangsdaten (insbesondere CMD-3 vs. frühere Messungen) zu klären, da die Gitter-QCD nun eine unabhängige, hochpräzise Referenz für die gesamte HVP-Korrektur (LO + NLO) liefert.