Lattice Calculation of Short-Range Contributions to Neutrinoless Double-Beta Decay $\pi^-\to\pi^+ ee$ at Physical Pion Mass

Diese Arbeit präsentiert eine Gitter-QCD-Berechnung der $\pi^- \to \pi^+ ee$-Matrixelemente bei der physikalischen Pionmasse unter Verwendung von Domain-Wall-Fermionen, wobei eine neuartige Methode zur Subtraktion von „um-die-Welt"-Effekten und eine nichtstörungstheoretische Renormierung eingesetzt werden, um Unsicherheiten zu verringern und Diskrepanzen in früheren Studien zu Beiträgen des Neutrinolosen Doppelbetazerfalls im Kurzreichweitenbereich aufzulösen.

Das Wesentliche

Das große Ganze: Ein kosmisches Rätsel lösen

Stellen Sie sich das Universum als ein riesiges Puzzle vor, und eines der fehlenden Teile ist die Natur der Neutrinos. Das sind winzige, geisterhafte Teilchen, die durch alles hindurchrasen. Wissenschaftler wollen wissen: Sind sie „Dirac"-Fermionen (wie normale Elektronen, bei denen Teilchen und Antiteilchen unterschiedlich sind) oder „Majorana"-Fermionen (bei denen ein Teilchen sein eigenes Antiteilchen ist)?

Der einzige Weg, dieses Rätsel zu lösen, besteht darin, nach einem sehr seltenen Ereignis zu suchen, das neutrinoloser Doppel-Beta-Zerfall genannt wird. Es ist wie das Beobachten, wie zwei Atome spontan in andere Atome verwandeln und zwei Elektronen ausspucken, aber ohne Neutrinos freizusetzen. Wenn wir dies geschehen sehen, beweist es, dass Neutrinos ihre eigenen Antiteilchen sind.

Das Problem: Ein verrauschtes Signal

Um vorherzusagen, ob dieses Ereignis stattfinden wird, müssen Physiker schwere Mathematik betreiben. Sie teilen die Berechnung in zwei Teile auf:

1. Der Langstrecken-Teil: Wie ein Flüstern, das einen Raum durchquert.
2. Der Kurzstrecken-Teil: Wie ein Schrei direkt neben Ihrem Ohr.

Dieses Paper konzentriert sich auf den Kurzstrecken-Teil. Konkret berechnen sie, wie zwei Pionen (Teilchen aus Quarks) wechselwirken, um zwei Elektronen zu erzeugen. Stellen Sie sich dies vor wie das Messen der „Lautstärke" dieses Schreis.

Der Konflikt: Zwei verschiedene Teams von Wissenschaftlern hatten zuvor versucht, diese „Lautstärke" mit Supercomputern (sogenanntes Gitter-QCD) zu messen. Ihre Ergebnisse unterschieden sich jedoch um den Faktor zwei. Es war, als würde ein Team sagen, der Schrei sei 60 Dezibel laut, und das andere, er sei 120 Dezibel laut. Diese enorme Uneinigkeit machte es schwierig, den Vorhersagen für das Neutrino-Rätsel zu vertrauen.

Die Lösung: Ein neuer Weg zur Bereinigung der Daten

Die Autoren dieses Papers beschlossen, ihr eigenes Experiment durchzuführen, um den Streit zu schlichten. Sie nutzten einen massiven Supercomputer, um die subatomare Welt zu simulieren. Doch sie stießen auf ein spezifisches technisches Problem: „Um-die-Welt-Effekte".

Die Analogie: Stellen Sie sich vor, Sie nehmen ein Gespräch in einem kleinen, hallenden Raum mit kreisförmigen Wänden auf. Wenn Sie in die Hände klatschen, breitet sich der Schall vorwärts aus, trifft auf die Wand, wickelt sich um den Raum herum und kommt von hinten zu Ihnen zurück. In der Computersimulation ist der „Raum" das Gitter der Raumzeit. Da das Gitter endlich ist, können die Teilchen den gesamten Loop umrunden und die Messung stören, wodurch ein verwirrendes „Echo" entsteht, das die Daten ruiniert.

Die Innovation: Bisherige Methoden versuchten zu erraten, wie man diese Echos auslöschen kann. Dieses Team erfand eine neue Subtraktionsmethode.

• Anstatt zu raten, isolierten sie das „Echo"-Signal direkt aus den Daten.
• Sie berechneten genau, wie stark das Echo war, und subtrahierten es vom Hauptsignal.
• Das Ergebnis: Das „Rauschen" verschwand und hinterließ ein sauberes, stabiles Signal (ein „Plateau"), dem sie vertrauen konnten.

Die Verifikation: Den Lineal-Check durchführen

Um sicherzustellen, dass ihre neue Methode nicht fehlerhaft war, überprüften sie ihre Arbeit gegen einen bekannten Standard. Sie berechneten einen spezifischen Wert (einen sogenannten „Bag-Parameter"), der zuvor von anderen Teams gemessen worden war.

• Ihr Ergebnis stimmte perfekt mit dem vertrauenswürdigen Standard überein.
• Als sie ihr Ergebnis mit dem Team verglichen, das den „Faktor-zwei"-Unterschied hatte, stellten sie fest, dass ihre Zahlen genau das Doppelte der Zahlen des anderen Teams waren.
• Die Schlussfolgerung: Es stellt sich heraus, dass das andere Team wahrscheinlich eine leicht andere „Messlatte" (Normierungskonvention) für ihre Messungen verwendet hat. Sobald man diesen Unterschied berücksichtigt, passen die Datenpunkte tatsächlich sehr gut zusammen. Die Methode der Autoren bestätigt, dass ihre Berechnung korrekt ist, und klärt die Verwirrung auf.

Das Endergebnis

Das Team hat erfolgreich den „Kurzstrecken"-Beitrag zum Prozess des neutrinolosen Doppel-Beta-Zerfalls mit viel höherer Präzision berechnet als zuvor.

• Sie entfernten die „Echos" (Um-die-Welt-Effekte), die frühere Daten durcheinandergebracht hatten.
• Sie verwendeten zwei verschiedene mathematische „Linsen" (Renormierungsschemata), um sicherzustellen, dass ihre Mathematik solide war.
• Sie lieferten eine definitive, hochpräzise Zahl, die Physikern hilft vorherzusagen, ob wir dieses seltene Zerfallssignal schließlich in realen Experimenten sehen werden.

Kurz gesagt: Sie bauten ein besseres Mikroskop, reinigten das statische Rauschen und bestätigten, dass die vorherige Uneinigkeit nur eine Frage der Verwendung unterschiedlicher Maßbänder war. Nun hat die wissenschaftliche Gemeinschaft eine zuverlässige Zahl, um das Rätsel der Neutrinomasse zu lösen.

Technische Zusammenfassung

1. Problemstellung

Der neutrinolose Doppelbeta-Zerfall ($0\nu\beta\beta$) ist ein entscheidender Prozess zur Bestimmung, ob Neutrinos Majorana-Fermionen sind, und zur Untersuchung der Verletzung der Leptonenzahl (LNV). Die Zerfallsrate hängt von nuklearen Matrixelementen (NMEs) ab. In der effektiven Feldtheorie des Standardmodells (SMEFT) entstehen kurzreichweitige Beiträge zur $0\nu\beta\beta$ aus Operatoren der Dimension 9. Diese Beiträge werden durch den Prozess $\pi^- \to \pi^+ ee$ vermittelt, der über Pionenaustausch zwischen Nukleonen erfolgt.

Das Kernproblem: Bisherige Gitter-QCD-Berechnungen der kurzreichweitigen Matrixelemente für $\pi^- \to \pi^+ ee$ (insbesondere Refs. [30] und [31]) berichteten über signifikant voneinander abweichende Ergebnisse. Für den Bag-Parameter $B_\pi^{MS}$ bei $\mu=3$ GeV berichtete Ref. [30] über $\approx 0.42$, während Ref. [31] über $\approx 0.20$ berichtete – eine Diskrepanz von etwa dem Faktor zwei. Diese Inkonsistenz behindert zuverlässige theoretische Vorhersagen für $0\nu\beta\beta$-Zerfallsraten. Darüber hinaus litten frühere Studien unter:

• Systematischen Unsicherheiten: Signifikante „Around-the-World" (ATW)-Effekte aufgrund der Rückwärtsausbreitung von Pionen in periodischen zeitlichen Gittern, welche die Extraktion von Matrixelementen des Grundzustands verfälschen.
• Chiraler Extrapolation: Einige Studien stützten sich auf unphysikalische Pionmassen, was eine Extrapolation zum physikalischen Punkt erforderte und zusätzliche Modellabhängigkeiten einführte.

2. Methodik

Die Autoren führten eine hochpräzise Gitter-QCD-Berechnung unter Verwendung von Domain-Wall-Fermion (DWF)-Ensembles durch, die von der RBC/UKQCD-Kollaboration generiert wurden.

A. Gitteraufbau

• Ensembles: Es wurden $N_f = 2+1$-Ensembles verwendet, darunter zwei mit physikalischen Pionmassen (48I und 64I) mit unterschiedlichen Gitterabständen ($a^{-1} \approx 1.73$ und $2.36$ GeV), um eine direkte Kontinuumsextrapolation zu ermöglichen. Auch Ensembles mit unphysikalischen Massen wurden für Kreuzprüfungen verwendet.
• Operatoren: Berechnung der Matrixelemente für fünf Vier-Quark-Operatoren der Dimension 9 ($O_1, O_2, O_3, O'_1, O'_2$), die für den Übergang $\pi^- \to \pi^+ ee$ relevant sind.

B. Minderung von „Around-the-World" (ATW)-Effekten

Eine wesentliche Innovation dieser Arbeit ist eine neue Methode zum Umgang mit ATW-Effekten, die auftreten, wenn Pionen über die periodische zeitliche Grenze hinweg propagieren und die Korrelationsfunktionen kontaminieren.

• Analyse: Die Autoren analysierten die zeitliche Struktur von Drei-Punkt-Korrelationsfunktionen und identifizierten Beiträge des gewünschten Diagramms (A) sowie ATW-Diagramme (B, C, D).
• Neue Subtraktionsmethode: Anstatt sich auf Verhältnismethoden (R1 oder R2) zu verlassen, die ATW-Effekte nur teilweise unterdrücken, schlugen die Autoren eine direkte Subtraktionstechnik vor.
  • Sie identifizierten Zeitbereiche, in denen der ATW-Beitrag eines einzelnen Pions, das die Grenze überquert (Diagramm B), dominiert.
  • Sie rekonstruierten den ATW-Beitrag direkt aus den Gitterdaten unter Verwendung spezifischer Zeitabstände ($t_{\pi\pi}, t_\pi$).
  • Formel: Das korrigierte Matrixelement wird durch Subtraktion des rekonstruierten ATW-Terms erhalten:
    $$O^{(sub)}_i(t) = N_\pi^{-2} e^{2m_\pi t} \left[ C^i_3(t,t) - 2 C^i_3(t_{\pi\pi}, T - t_{\pi\pi} - t_\pi) e^{-m_\pi(T-t_\pi-2t_{\pi\pi})} \right]$$
  • Diese Methode stellte erfolgreich stabile Plateaus in den effektiven Massenplots wieder her und reduzierte die systematischen Fehler im Vergleich zu früheren Verhältnismethoden erheblich.

C. Renormierung

• Schema: Eine nichtstörungstheoretische Renormierung wurde mit der RI/SMOM-Methode (Rome-Southampton) unter Verwendung nicht-exzeptioneller Kinematik durchgeführt.
• Schemata: Zwei Schemata wurden eingesetzt, um Fehler durch störungstheoretische Trunkierung abzuschätzen: $(\gamma_\mu, \gamma_\mu)$ und $(\not{q}, \not{q})$.
• Matching: Die Ergebnisse wurden unter Verwendung von Matching-Koeffizienten auf Ein-Schleifen-Niveau in das $\overline{MS}$-Schema bei $\mu = 3$ GeV umgerechnet. Der Unterschied zwischen den beiden RI/SMOM-Schemata wurde zur Abschätzung der systematischen Unsicherheit durch störungstheoretische Trunkierung verwendet.

3. Wichtige Ergebnisse

A. Bloße Matrixelemente

Die Studie extrahierte bloße Matrixelemente $\langle \pi^+ | O_i | \pi^- \rangle$ bei der physikalischen Pionmasse.

• Kreuzprüfung mit Bag-Parametern: Die Autoren überprüften ihre Normierung, indem sie das Matrixelement von $O_3$ (bezogen auf den neutralen Meson-Bag-Parameter $B_K$) mit früheren, von FLAG rezensierten Ergebnissen verglichen. Die Ergebnisse stimmten perfekt mit Ref. [38] überein (welches dieselbe Normierungskonvention verwendete), waren jedoch etwa zweimal so groß wie die Ergebnisse in Ref. [31].
• Normierungsdiskrepanz: Das Verhältnis der Ergebnisse der Autoren zu Ref. [31] beträgt für alle fünf Operatoren konsistent $\approx 2.0$. Dies deutet auf einen globalen Normierungsfaktorfehler in Ref. [31] hin, und nicht auf eine physikalische Diskrepanz.

B. Renormierte Matrixelemente und Bag-Parameter

Die Endergebnisse im $\overline{MS}$-Schema ($\mu = 3$ GeV) lauten:

• Bag-Parameter: $B_\pi^{MS}(3 \text{ GeV}) = 0.3769(24)_{stat}(76)_{PT}(1)_L$.
  • Dieser Wert ist innerhalb der Unsicherheiten mit Ref. [30] ($0.421$) konsistent, unterscheidet sich jedoch signifikant von Ref. [31] ($0.197$).
  • Das Verhältnis zu Ref. [31] beträgt $1.91(18)$, was das Normierungsproblem des Faktors 2 bestätigt.
• Kontinuumsextrapolation: Unter Verwendung der Ensembles mit physikalischer Masse (48I und 64I) führten die Autoren eine direkte Kontinuumsextrapolation durch, wodurch Unsicherheiten im Zusammenhang mit der chiralen Extrapolation eliminiert wurden.
• Reduktion der Unsicherheit: Die Gesamtunsicherheiten in dieser Arbeit sind im Vergleich zu früheren Studien erheblich reduziert aufgrund von:
  1. Direkter Berechnung bei physikalischer Pionmasse.
  2. Verbesserter statistischer Präzision.
  3. Effektiver Beseitigung von ATW-systematischen Fehlern.
  4. Robuster Abschätzung von störungstheoretischen Trunkierungsfehlern.

4. Bedeutung und Schlussfolgerung

• Auflösung der Diskrepanz: Diese Arbeit liefert eine unabhängige Kreuzprüfung, die die in Ref. [30] und im FLAG-Review verwendete Normierung stark unterstützt, und identifiziert gleichzeitig einen wahrscheinlichen Normierungsfehler in Ref. [31]. Die Diskrepanz des Faktors 2 wird als Problem der Normierungskonvention und nicht als fundamentale physikalische Meinungsverschiedenheit aufgelöst.
• Methodischer Fortschritt: Die vorgeschlagene Subtraktionsmethode für ATW-Effekte ist eine robuste Technik, die auf andere Gitterberechnungen angewendet werden kann, die leichte Mesonen in periodischen Volumina beinhalten, und verbessert die Zuverlässigkeit zukünftiger Berechnungen von $0\nu\beta\beta$-Matrixelementen.
• Auswirkung auf die Physik: Durch die Bereitstellung präziser, kontinuumsextrapolierter Matrixelemente bei physikalischen Pionmassen reduziert diese Arbeit die nichtstörungstheoretischen QCD-Unsicherheiten in der theoretischen Vorhersage von $0\nu\beta\beta$-Zerfallsraten. Dies ist entscheidend für die Interpretation experimenteller Ergebnisse aus Experimenten der nächsten Generation (z. B. LEGEND, nEXO, CUPID), um die Natur der Neutrinomasse und die Skala neuer Physik zu bestimmen.

Zusammenfassend etabliert dieses Papier einen neuen Benchmark für kurzreichweitige Beiträge zum neutrinolosen Doppelbeta-Zerfall, indem es frühere Inkonsistenzen durch eine rigorose Methodik, Simulationen mit physikalischen Massen und eine neuartige Behandlung von endlichen-Volumen-Artefakten auflöst.