Comparison of the hadronic vacuum polarization between hadronic $\tau$-decay data and lattice QCD

Diese Arbeit vergleicht isospinsymmetrische Gitter-QCD-Berechnungen der hadronischen Vakuumpolarisation mit dispersiven Ergebnissen, die aus korrigierten Daten hadronischer $\tau$-Zerfälle abgeleitet wurden, und findet insgesamt eine allgemein gute Übereinstimmung, deckt jedoch signifikante Diskrepanzen im Vier-Pion-Modus $2\pi^-\pi^+\pi^0$ auf, wenn diese im Vergleich zu den Erwartungen aus Pais-Relationen und $e^+e^-$-Wirkungsquerschnitten bewertet werden.

Das Wesentliche

Das große Ganze: Eine winzige Wackelbewegung messen

Stellen Sie sich das Universum als eine riesige, komplexe Maschine vor. Einer ihrer bekanntesten Teile ist das Myon, ein Teilchen, das wie ein winziger, sich drehender Kreisel wirkt. Wissenschaftler haben gemessen, wie dieser Kreisel wackelt (sein „anomales magnetisches Moment") mit unglaublicher Präzision.

Um jedoch genau vorherzusagen, wie stark er sollte wackeln, basierend auf unseren aktuellen physikalischen Regeln (dem Standardmodell), müssen Wissenschaftler einen „Nebel" virtueller Teilchen berücksichtigen, die um das Myon herum ins und aus dem Nichts auftauchen. Dieser Nebel wird Hadronische Vakuumpolarisation (HVP) genannt.

Das Problem ist, dass die Berechnung dieses Nebels unglaublich schwierig ist. Es gibt zwei Hauptmethoden, mit denen Wissenschaftler versuchen, ihn zu messen:

1. Die „Gitter"-Methode: Verwendung von Supercomputern, um die Gesetze der Physik von Grund auf zu simulieren (wie den Bau eines digitalen Modells des Nebels).
2. Die „Daten"-Methode: Betrachtung realer Experimente, bei denen Teilchen zusammenstoßen, um diesen Nebel zu erzeugen, und anschließende Messung der Ergebnisse.

Lange Zeit stimmten diese beiden Methoden nicht überein. Die „Gitter"-Ergebnisse und die „Daten"-Ergebnisse passten nicht zusammen, was ein Rätsel in der Physik schuf.

Das neue Experiment: Eine andere Kamera verwenden

Dieses Papier versucht, das Rätsel zu lösen, indem es eine andere Art von „Kamera" für die Daten-Methode verwendet.

Normalerweise betrachten Wissenschaftler Daten aus Elektron-Positron-Kollisionen (das Zusammenstoßen eines Elektrons und eines Positrons). Aber dieses Papier verwendet Daten aus Tau-Zerfällen.

• Die Analogie: Stellen Sie sich vor, Sie versuchen, die Form einer bestimmten Wolkensorte zu messen.
  • Methode A (Elektron-Kollisionen): Sie schauen durch ein Teleskop, das manchmal ein wenig statische Störung bekommt (genannt „Isospin-Bruch").
  • Methode B (Tau-Zerfälle): Sie schauen durch ein anderes Teleskop, das einen leicht anderen Winkel sieht.
  • Das Ziel: Die Autoren nehmen die „Tau"-Daten, bereinigen sie, um die statischen Störungen zu entfernen (Korrektur der physikalischen Unterschiede zwischen den beiden Methoden), und vergleichen sie mit der „Gitter"-Computersimulation.

Was sie taten

Die Autoren verwendeten eine massive Menge an Daten aus Tau-Teilchen-Zerfällen (ein schwerer Cousin des Elektrons). Sie konzentrierten sich darauf, wie diese Teilchen in kleinere Stücke zerfallen (wie Pionen).

1. Bereinigung der Daten: Die Tau-Daten sind nicht perfekt; sie weisen winzige Unterschiede im Vergleich zur idealen „reinen" physikalischen Welt auf, die in den Computersimulationen verwendet wird. Die Autoren bauten einen mathematischen „Filter", um diese Unterschiede zu korrigieren, und übersetzten im Wesentlichen die Tau-Daten in die Sprache der Computersimulation.
2. Der Vergleich: Sie verglichen diese bereinigten Tau-Daten mit den Ergebnissen der Mainzer und BMW-Supercomputer-Gruppen (den Gitter-Teams).

Die Ergebnisse: Gute Nachrichten und ein seltsamer Fehler

1. Die guten Nachrichten (Allgemeine Übereinstimmung)
Zum größten Teil stimmten die beiden Methoden sehr gut überein.

• Die Analogie: Es ist wie bei zwei verschiedenen Wetterstationen, die die Temperatur messen. Obwohl sie verschiedene Thermometer verwenden, sagen beide, es seien 22 °C.
• Die Erkenntnis: Als sie den gesamten „Nebel" (den Beitrag zum Myon-Wackeln) und die „mittleren Entfernungs"-Teile davon betrachteten, passten die tau-basierten Daten und die Gitter-Computersimulationen gut zusammen. Dies legt nahe, dass die Computersimulationen wahrscheinlich korrekt sind und dass die früheren Uneinigkeiten möglicherweise auf Probleme mit den Elektron-Positron-Daten zurückzuführen waren, nicht auf die Computermodelle.

2. Der seltsame Fehler (Das Vier-Pionen-Problem)
Allerdings fanden sie eine spezifische Stelle, an der die Daten nicht mit den Regeln des Universums übereinstimmten.

• Die Analogie: Stellen Sie sich vor, Sie backen einen Kuchen. Sie haben ein Rezept (die „Pais-Relationen"), das besagt, dass Sie bei Mischung von 4 Eiern und 2 Tassen Mehl ein bestimmtes Ergebnis erhalten.
  • Als sie sich eine bestimmte Kuchensorte ansahen (der 2π−π+π0-Modus, oder eine spezifische Art, wie vier Teilchen zerfallen), sagten die „Tau"-Daten, der Kuchen habe eine Größe, aber die „Elektron"-Daten sagten, er habe eine andere Größe.
  • Die Autoren prüften dies gegen das „Rezept" (theoretische Regeln) und fanden einen signifikanten Unterschied. Die Tau-Daten für diese spezifische Vier-Teilchen-Kombination stimmten nicht mit dem überein, was die Elektron-Daten und die theoretischen Regeln vorhersagten.

Die Schlussfolgerung

• Insgesamt: Das Papier stellt fest, dass bei Verwendung von Tau-Zerfallsdaten (korrekt bereinigt) eine sehr gute Übereinstimmung mit Gitter-QCD (den Supercomputersimulationen) besteht. Dies unterstützt die Idee, dass die Supercomputer-Ergebnisse wahrscheinlich die richtigen sind.
• Der Vorbehalt: Es gibt einen spezifischen, komplexen Teil der Daten (der das Zerfallen von vier Teilchen auf eine bestimmte Weise betrifft), bei dem die Tau-Daten und die Elektron-Daten erheblich voneinander abweichen. Dies deutet darauf hin, dass es möglicherweise ein Problem mit der Messung oder dem Verständnis dieses spezifischen Teils des Teilchenzerfalls gibt, aber es zerstört nicht die allgemeine Übereinstimmung für die Hauptberechnung.

Kurz gesagt: Die Autoren verwendeten eine neue Art von Daten (Tau-Zerfälle), um die Computersimulationen zu überprüfen. Die Prüfung bestand für das große Ganze und bestätigte die Computermodelle, hob jedoch ein spezifisches, verwirrendes Detail in den Daten hervor, das noch geklärt werden muss.

Technische Zusammenfassung

Technische Zusammenfassung: Vergleich der hadronischen Vakuumpolarisation zwischen Daten hadronischer $\tau$-Zerfälle und Gitter-QCD

Problemstellung
Die Bestimmung des Erwartungswerts des Standardmodells (SM) für das anomale magnetische Moment des Myons, $a_\mu$, wird derzeit durch die Unsicherheit im Beitrag der hadronischen Vakuumpolarisation (HVP), $a_\mu^{\text{HVP}}$, begrenzt. Es besteht eine signifikante Diskrepanz zwischen Gitter-QCD-Berechnungen des Beitrags der hadronischen Vakuumpolarisation für die Verbindung mit leichten Quarks (lqc) und dispersiven Schätzungen, die aus dem $R$-Verhältnis gemessen in $e^+e^- \to \text{Hadronen}(\gamma)$-Prozessen abgeleitet werden. Diese Spannung wird primär dem exklusiven Kanal $e^+e^- \to \pi^+\pi^-$ um die $\rho$-Resonanz herum zugeschrieben. Da der isospin-1 ($I=1$)-Anteil der HVP im Isospin-Limit direkt mit dem lqc-Beitrag verknüpft ist und da $I=1$-Spektralfunktionen auch über nicht-strange hadronische $\tau$-Zerfälle zugänglich sind (nach Korrektur für Isospin-Bruch), untersucht diese Arbeit, ob $\tau$-Zerfallsdaten eine unabhängige dispersive Bestimmung liefern können, die mit Gitter-QCD-Ergebnissen übereinstimmt.

Methodik
Die Autoren konstruieren die isospin-1-Vektorstrom-subtrahierte Vakuumpolarisation, $\Pi^{I=1}(Q^2)$, und die zugehörige euklidische Zeit-Korrelationsfunktion, $C(t)$, unter Verwendung einer dispersiven Darstellung. Die Spektralfunktion, $\rho_V^{I=1}(s)$, wird aus drei Komponenten aufgebaut:

1. Datenbereich ($s \leq m_\tau^2$): Die inklusive vektor-isovektore Spektralfunktion wird aus nicht-strangen hadronischen $\tau$-Zerfällen abgeleitet. Die Autoren nutzen eine Kombination von Daten aus den Experimenten ALEPH, OPAL und Belle, wie in Ref. [19] zusammengestellt. Dies umfasst exklusive Beiträge von Zwei-Pion- ($\pi^-\pi^0$), Vier-Pion- ($2\pi^-\pi^+\pi^0$ und $\pi^-3\pi^0$) sowie Restmodi.
2. Störungstheoretischer Bereich ($s > m_\tau^2$): Für Invariantmassen oberhalb der $\tau$-Masse, wo keine experimentellen Daten verfügbar sind, wird die Spektralfunktion unter Verwendung von masseloser störungstheoretischer QCD (PT) im fünf-Schleifen-Formalismus mit einem geschätzten Koeffizienten fünfter Ordnung modelliert.
3. Dualitätsverletzungen (DV): Ein phänomenologisches Modell wird eingesetzt, um Dualitätsverletzungen Quark-Hadron im Bereich oberhalb von $m_\tau^2$ zu berücksichtigen.

Um einen direkten Vergleich mit Gitter-QCD-Ergebnissen zu ermöglichen, die in „isoQCD" (isospinsymmetrische QCD ohne elektromagnetische Korrekturen) berechnet wurden, unterliegen die auf $\tau$ basierenden Spektraldaten einer rigorosen Isospin-Bruch-Korrektur (IB). Die Autoren gehen davon aus, dass IB-Effekte vom Zwei-Pion-Modus dominiert werden. Der Korrekturfaktor $F_{\text{IB}}(s)$ berücksichtigt:

• Unterschiede in kurzreichweitigen elektroschwachen Faktoren ($S_{\text{EW}}$ vs. $S_{\pi\pi}^{\text{EW}}$).
• Strahlungskorrekturen ($G_{\text{EM}}(s)$), die aus dispersiven Analysen abgeleitet wurden.
• Phasenraum-Anpassungen aufgrund von Pion-Massendifferenzen ($m_{\pi^\pm} \neq m_{\pi^0}$).
• Isospin-Bruch-Effekte im Pion-Formfaktor $f_+(s)$, modelliert unter Verwendung eines durch chirale Störungstheorie (ChPT) inspirierten Ansatzes.

Die Studie bewertet mehrere Größen:

• Die subtrahierte Vakuumpolarisation $\Pi(Q^2)$ für $Q^2 \in [0.5, 12] \text{ GeV}^2$.
• Den gesamten lqc-Beitrag zu $a_\mu$, bezeichnet als $a_\mu^{\text{HVP, lqc}}$.
• Die Kurzreichweiten-(SD), Zwischen-(W1) und Langreichweiten-(LD) Fenstergrößen, definiert von RBC/UKQCD.

Zusätzlich führen die Autoren alternative Analysen durch, bei denen $\tau$-Daten nur für spezifische exklusive Modi (2$\pi$ oder 2$\pi$+4$\pi$) verwendet werden, während andere Beiträge aus der KNT19 $e^+e^-$-Zusammenstellung entnommen werden, um den Einfluss spezifischer Kanal-Diskrepanzen zu isolieren.

Hauptbeiträge und Ergebnisse

1. Allgemeine Übereinstimmung: Die Arbeit findet im Allgemeinen eine gute Übereinstimmung zwischen den auf $\tau$ basierenden dispersiven Ergebnissen und Gitter-QCD-Berechnungen für $a_\mu^{\text{HVP, lqc}}$ und die Fenstergrößen. Die auf $\tau$ basierenden Ergebnisse liegen innerhalb von $1\sigma$ der Gitter-Durchschnitte (speziell die WP25-Durchschnitte aus Ref. [3]) für die gesamte HVP und das LD-Fenster, obwohl die SD- und W1-Fenster etwas größere relative Unterschiede zeigen (bis zu $1,0\sigma$).
2. Spannung bei $\Pi(Q^2)$: Für die subtrahierte Vakuumpolarisation $\Pi(Q^2)$ liegen die auf $\tau$ basierenden Ergebnisse systematisch niedriger als die Mainzer Gitter-Ergebnisse (Ref. [14]) um $2,1\sigma$ bis $3,0\sigma$ im Bereich $Q^2 = 0,5$ bis $12 \text{ GeV}^2$. Diese Ergebnisse zeigen jedoch eine bessere Übereinstimmung mit den BMW-Gitter-Ergebnissen (Ref. [15]).
3. Diskrepanzen im Vier-Pion-Modus: Ein detaillierter Vergleich der Vier-Pion-Spektralverteilungen zeigt signifikante Spannungen zwischen $\tau$-basierten Daten und Erwartungen, die aus $e^+e^-$-Daten über Pais-Relationen abgeleitet wurden. Insbesondere zeigt der Modus $2\pi^-\pi^+\pi^0$ Diskrepanzen im Bereich von $3,3\sigma$ (SD-Fenster) bis $4,4\sigma$ (LD-Fenster) beim Vergleich von ALEPH $\tau$-Daten mit KNT19 $e^+e^-$-Kombinationen. Der Modus $\pi^-3\pi^0$ zeigt eine gute Übereinstimmung.
4. Auswirkung des Modus-Austauschs: Wenn die auf $\tau$ basierende 2$\pi$-Spektralverteilung verwendet wird, um die prä-CMD-3 KNT19 $e^+e^-$ 2$\pi$-Eingabe zu ersetzen, stimmen die daraus resultierenden datengesteuerten lqc-Beiträge gut mit Gitter-Ergebnissen überein. Wenn jedoch auch die auf $\tau$ basierenden 4$\pi$-Daten (speziell der Modus $2\pi^-\pi^+\pi^0$) verwendet werden, um die KNT19-Eingabe zu ersetzen, bleiben die Ergebnisse innerhalb der Fehler mit den Gitter-Durchschnitten vereinbar, führt jedoch die starke Korrelation zwischen den 2$\pi$- und 4$\pi$-Beiträgen zu einer $3,2\sigma$-Spannung zwischen den Strategien „nur 2$\pi$-$\tau$" und „2$\pi$+4$\pi$-$\tau$" für das LD-Fenster.
5. Fehlerhaushalt: Die Autoren stellen fest, dass für niedrigere $Q^2$ der datengesteuerte Teil von $\Pi(Q^2)$ den Fehlerhaushalt dominiert (ca. 91 % bei $Q^2=0,5 \text{ GeV}^2$), was die $\tau$-basierten Fehler mit Gitter-Fehlern konkurrenzfähig macht.

Bedeutung
Die Arbeit zeigt, dass die Verwendung inklusiver hadronischer $\tau$-Zerfallsdaten, korrigiert für Isospin-Bruch, eine dispersive Bestimmung der HVP liefert, die weitgehend mit aktuellen Gitter-QCD-Berechnungen konsistent ist. Dies stützt die Hypothese, dass die Diskrepanz zwischen Gitter- und $e^+e^-$-basierten HVP-Schätzungen auf Probleme spezifisch für die $e^+e^-$-Daten (insbesondere den $\pi^+\pi^-$-Kanal) zurückzuführen sein könnte, und nicht auf ein fundamentales Versagen der Gitter-Methode oder des dispersiven Ansatzes selbst.

Die Autoren heben jedoch hervor, dass signifikante Diskrepanzen in den Vier-Pion-Kanälen zwischen $\tau$-Daten und aus Pais-Relationen abgeleiteten $e^+e^-$-Erwartungen bestehen bleiben. Obwohl diese Diskrepanzen die allgemeine Übereinstimmung zwischen den gesamten $\tau$-basierten und Gitter-HVP-Werten derzeit nicht ungültig machen, deuten sie darauf hin, dass der Vier-Pion-Sektor weiterer Untersuchung bedarf. Die Arbeit unterstreicht die Nützlichkeit von $\tau$-Zerfallsdaten als unabhängige Kreuzprüfung für den HVP-Beitrag zu $a_\mu$, insbesondere im Kontext der Auflösung der Anomalie des Myon $g-2$. Die Autoren bleiben bescheiden und stellen fest, dass die systematischen Unsicherheiten bei der Modellierung von Isospin-Bruch-Effekten (insbesondere im Pion-Formfaktor) und die beobachteten Spannungen in den Vier-Pion-Modi eine definitive Klärung aller Diskrepanzen in diesem Stadium verhindern.