Theoretical framework for lattice QCD computations of $B\to K \ell^+ \ell^-$ and $\bar{B}_s\to \ell^+\ell^- \gamma$ decays rates, including contributions from charming penguin diagrams

Dieses Papier entwickelt eine Strategie zur Berechnung der Zerfallsamplituden für $B\to K\ell^+\ell^-$ und $\bar{B}_s\to\gamma\ell^+\ell^-$ mittels Gitter-QCD, wobei komplexe Beiträge wie „charming penguins" und chromomagnetische Operatoren durch spektrale Dichtemethoden und nicht-perturbative Renormierung behandelt werden.

Das Wesentliche

Stellen Sie sich das Universum der subatomaren Teilchen als eine riesige, extrem komplexe Fabrik vor. In dieser Fabrik werden Teilchen wie der B-Meson (ein schweres, instabiles Teilchen) produziert, die dann in leichtere Teilchen zerfallen. Physiker wollen genau verstehen, wie dieser Zerfall abläuft, besonders bei Prozessen, die im Standardmodell der Physik sehr selten sind. Diese Seltenheit macht sie zu perfekten „Fenstern", um nach neuen, unbekannten Kräften oder Teilchen zu suchen, die das Standardmodell erweitern könnten.

Das Problem: Um diese Zerfälle vorherzusagen, müssen wir die Wechselwirkungen der starken Kernkraft berechnen. Das ist wie der Versuch, das Wetter in einer einzelnen Zelle eines Ameisenhaufens zu berechnen, während tausende andere Ameisen herumlaufen. Die Mathematik ist so kompliziert, dass man sie nicht einfach mit einem Taschenrechner lösen kann.

Hier kommt das Gitter-QCD (Quantenchromodynamik auf dem Gitter) ins Spiel. Man stellt sich die Raumzeit wie ein dreidimensionales Schachbrett vor. Die Teilchen sind die Figuren auf diesem Brett. Indem man Computer nutzt, um die Bewegung der Figuren auf diesem Gitter zu simulieren, kann man die starken Kräfte berechnen.

Das große Hindernis: Die „Geister" der Zwischenwelt
Bei bestimmten Zerfällen (wie $B \to K \ell^+ \ell^-$ oder $B_s \to \gamma \ell^+ \ell^-$) passiert etwas Tückisches: Auf dem Weg vom Anfangs- zum Endzustand entstehen kurzlebige Zwischenzustände, die „auf der Schiene" (on-shell) existieren können.

• Die Analogie: Stellen Sie sich vor, Sie werfen einen Ball von Punkt A nach Punkt B. Normalerweise fliegt er direkt. Aber in der Quantenwelt kann der Ball kurzzeitig in eine unsichtbare „Geister-Box" springen, dort tanzen und dann wieder herauskommen. Diese „Tänze" sind komplex und haben sowohl einen realen als auch einen imaginären Teil (in der Mathematik).
• Das Problem für die Computer-Simulation: Unsere Computer arbeiten in einer „Euklidischen Welt" (eine Art mathematischer Spiegelwelt), in der Zeit anders läuft als in der echten Welt. In dieser Spiegelwelt sind die „Geister-Tänze" unsichtbar oder verschwinden. Man kann sie nicht direkt sehen, genau wie man ein Echo nicht hören kann, wenn man die Ohren zuhält. Bisher mussten Physiker diese Effekte schätzen oder mit Modellen raten, was zu großen Unsicherheiten führte.

Die Lösung: Der „Spektral-Dichte"-Trick
Die Autoren dieses Papiers haben einen neuen Weg gefunden, um diese „Geister" sichtbar zu machen.

• Die Analogie: Stellen Sie sich vor, Sie hören ein Lied, das von einem Orchester gespielt wird, aber Sie können nur die Summe aller Töne hören, nicht die einzelnen Instrumente. Um das Lied zu verstehen, müssen Sie das Signal „zerlegen".
• Die Autoren nutzen eine Methode namens SFR (Spectral Function Reconstruction) in Kombination mit HLT. Man kann sich das wie einen sehr cleveren Filter vorstellen. Anstatt zu versuchen, das Echo direkt zu hören, simulieren sie das Orchester mit verschiedenen „Verzerrungen" (einem Parameter namens $\epsilon$). Indem sie die Ergebnisse bei verschiedenen Verzerrungen messen und dann mathematisch zurückrechnen (extrapolieren), können sie das ursprüngliche, reine Signal rekonstruieren, auch wenn es in der Computer-Simulation eigentlich „versteckt" war.

Die „Charming Penguins" (Die Charmanten Pinguine)
Ein besonderer Fokus liegt auf Diagrammen, die als „Charming Penguins" bekannt sind.

• Was ist das? In diesen Zerfällen taucht kurzzeitig ein Paar aus einem Charm-Quark und seinem Antiteilchen auf (ein „Loop").
• Warum Pinguine? Der Name ist eine Hommage an die Entdeckung, dass diese Effekte oft übersehen wurden, aber plötzlich „auftauchen" und das Ergebnis stark beeinflussen, ähnlich wie Pinguine, die unerwartet aus dem Eis kommen.
• Das Risiko: Diese „Pinguine" können die Vorhersagen so stark verändern, dass sie wie ein Signal für „Neue Physik" aussehen, obwohl sie nur ein bekannter, aber schwer zu berechnender Effekt sind. Wenn man sie falsch berechnet, könnte man glauben, eine neue Kraft entdeckt zu haben, wo es nur einen Rechenfehler gibt.

Was haben die Autoren getan?

1. Theorie: Sie haben einen vollständigen Bauplan entwickelt, wie man diese „Geister"-Effekte und die „Charming Penguins" mit Gitter-QCD berechnen kann, ohne auf Schätzungen angewiesen zu sein. Sie haben auch erklärt, wie man mathematische Unendlichkeiten (die bei solchen Berechnungen auftreten) sauber entfernt.
2. Testlauf: Sie haben einen ersten, kleinen Testlauf durchgeführt. Sie haben eine vereinfachte Version des Problems gelöst (mit einem leichteren B-Meson, um die Rechenzeit zu sparen).
3. Ergebnis: Der Test war erfolgreich! Sie konnten zeigen, dass ihre Methode funktioniert. Die Ergebnisse sahen gut aus und stimmten mit einfachen Modellen überein, wo man es erwarten würde. Besonders wichtig: Sie konnten die komplexen, imaginären Teile der Amplitude (die „Geister") tatsächlich berechnen.

Warum ist das wichtig?
Bisher waren die Vorhersagen für diese Zerfälle unscharf, wie ein Foto, das unscharf ist. Mit dieser neuen Methode können die Physiker das Foto scharfstellen.

• Wenn die neuen, präzisen Berechnungen immer noch nicht mit den experimentellen Messungen (z. B. vom LHCb-Experiment am CERN) übereinstimmen, dann wissen wir zu 100 %, dass es Neue Physik gibt.
• Wenn sie übereinstimmen, haben wir unser Verständnis des Standardmodells vertieft.

Zusammenfassung für den Alltag:
Stellen Sie sich vor, Sie versuchen, die genaue Schwingung einer Gitarrensaite zu messen, während ein lauter Wind weht. Bisher haben Sie nur geschätzt, wie stark der Wind die Saite beeinflusst. Diese Forscher haben nun ein neues Mikrofon und einen cleveren Algorithmus entwickelt, der den Wind herausfiltert und Ihnen die reine Schwingung der Saite zeigt. Damit können sie endlich genau sagen, ob die Saite (das Teilchen) so klingt, wie die Theorie es sagt, oder ob ein neues, unbekanntes Instrument (Neue Physik) im Orchester mitspielt.

Dieser Schritt ist ein Meilenstein, um die Grenzen unseres Wissens über das Universum zu erweitern und vielleicht eines Tages eine völlig neue Theorie der Physik zu finden.

Technische Zusammenfassung

Titel und Autoren

Titel: Theoretical framework for lattice QCD computations of $B \to K\ell^+\ell^-$ and $\bar{B}_s \to \ell^+\ell^-\gamma$ decay rates, including contributions from charming penguin diagrams
Autoren: Roberto Frezzotti, Nazario Tantalo, Giuseppe Gagliardi, Vittorio Lubicz, Guido Martinelli, Chris T. Sachrajda, Francesco Sanfilippo, Silvano Simula, Luca Silvestrini.

---

1. Problemstellung und Motivation

Flavour-changing neutral current (FCNC) Zerfälle wie $B \to K\ell^+\ell^-$ und $\bar{B}_s \to \gamma\ell^+\ell^-$ sind im Standardmodell (SM) stark unterdrückt und stellen daher ein sensibles Fenster für die Suche nach neuer Physik (NP) dar. Ein Hauptproblem bei der theoretischen Vorhersage dieser Zerfälle ist die Berechnung der hadronischen Matrixelemente.

• Herausforderung: Viele Beiträge zu diesen Amplituden, insbesondere solche mit komplexen Anteilen (Imaginärteilen), entstehen durch on-shell Zwischenzustände, die zwischen dem schwachen Operator und dem elektromagnetischen Strom propagieren. Dazu gehören sogenannte "Charming Penguins" (Diagramme mit Charm-Quark-Schleifen, z. B. Operatoren $O_1^{(c)}$ und $O_2^{(c)}$) sowie Beiträge vom chromomagnetischen Operator $O_8$.
• Limitierung herkömmlicher Methoden: Standard-Lattice-QCD-Techniken basieren auf euklidischen Korrelationsfunktionen, die reell sind. Sie können die komplexen Anteile der Amplituden nicht direkt berechnen, da diese eine analytische Fortsetzung von Minkowski- in euklidische Raumzeit erfordern, was bei on-shell Zwischenzuständen (unterhalb der äußeren Energien) zu Singularitäten führt.
• Aktueller Stand: Bisher wurden diese Beiträge oft durch phänomenologische Modelle (z. B. Vektor-Meson-Dominanz) abgeschätzt, was zu großen Unsicherheiten führt.

2. Methodik und Theoretischer Rahmen

Das Paper entwickelt einen vollständigen theoretischen Rahmen, um diese komplexen Beiträge aus ersten Prinzipien (First Principles) mittels Gitter-QCD zu berechnen.

• Spektrale Dichtemethoden (Spectral Density Methods): Die Autoren nutzen den Spectral Function Reconstruction (SFR) Ansatz (basierend auf Ref. [24]) in Kombination mit der Hansen-Lupo-Tantalo (HLT) Technik (Ref. [23]).
  • Die Methode nutzt die spektrale Darstellung zeitgeordneter Korrelationsfunktionen.
  • Sie ermöglicht die Trennung von Real- und Imaginärteilen der Amplituden, selbst wenn on-shell Zwischenzustände (wie $J/\psi$-Resonanzen) im Spiel sind.
  • Die Amplitude wird durch eine "geglättete" (smeared) Version berechnet, die als Integral über die spektrale Dichte mit einem Kernel (oft $1/(E - m_B - i\epsilon)$) dargestellt wird.
• Erweiterung auf Mehrpunkt-Operatoren: Das Paper generalisiert die Methode von bilokalen auf trilokale Operatoren (für $\bar{B}_s \to \gamma\ell^+\ell^-$), was sechs mögliche Zeitordnungen erfordert, von denen drei Imaginärteile beinhalten.
• Renormierung und Divergenzen:
  • Kontaktterme: Wenn der elektromagnetische Strom und der schwache Operator kollidieren, treten logarithmische UV-Divergenzen auf, die zu den Operatoren $O_7$ und $O_9$ im effektiven Hamiltonian führen.
  • Potenz-Divergenzen: Die Operatoren $O_1^{(c)}$ und $O_2^{(c)}$ mischen mit Operatoren niedrigerer Dimension (Skalare und Pseudoskalare Dichten), was zu Divergenzen in Form von Inversen der Gitterkonstante ($1/a$, $1/a^2$) führt.
  • Lösung: Die Autoren zeigen, wie diese Potenz-Divergenzen nicht-störungstheoretisch subtrahiert werden können, insbesondere unter Verwendung von Twisted-Mass Fermionen mit maximaler Twist. Durch die Nutzung spezifischer Symmetrien (CS, $S_3$, $R_5^{sp}$) und die Einführung von "Replicas" (zusätzlichen Valenzquark-Sorten) können die Mischungen eliminiert werden, indem Matrixelemente zwischen bestimmten Zuständen auf Null gesetzt werden.

3. Schlüsselbeiträge

1. Vollständiger Framework: Entwicklung einer Strategie zur Berechnung der gesamten Amplitude für $B \to K\ell^+\ell^-$ und $\bar{B}_s \to \gamma\ell^+\ell^-$, einschließlich der bisher unzugänglichen Charming-Penguin-Beiträge und des Operators $O_8$.
2. Generalisierung der SFR-Methode: Demonstration, dass die spektrale Dichtemethode auf Prozesse mit drei lokalen Operatoren (trilokal) und mehreren Zeitordnungen erweitert werden kann.
3. Behandlung von Divergenzen: Präsentation eines Verfahrens zur nicht-störungstheoretischen Subtraktion von Potenz-Divergenzen bei Vier-Quark-Operatoren in Gitter-QCD mit Twisted-Mass Fermionen.
4. Proof-of-Principle Berechnung: Durchführung einer ersten explorativen numerischen Studie für den Charming-Penguin-Beitrag in $B \to K\ell^+\ell^-$.

4. Ergebnisse der explorativen Studie

Die Autoren führten eine Berechnung auf einem einzigen Gitter-Ensemble (ETMC, $N_f = 2+1+1$, $a \approx 0.08$ fm) durch.

• Setup: Verwendung einer leichteren als der physikalischen $b$-Quark-Masse ($m_b = 2m_c$) zur Vermeidung großer Gitterartefakte. Berechnung bei einem Kaon-Impuls von $|\vec{p}_K| \approx 250$ MeV.
• Vergleich mit Modellen: Die berechneten euklidischen Korrelationsfunktionen wurden mit einem Modell basierend auf der Vakuum-Sättigungs-Näherung (VSA) verglichen.
  • Für den Operator $O_1^{(c)}$ zeigte sich eine qualitative Übereinstimmung, wobei ein multiplikativer Faktor von ca. 2.3 benötigt wurde, um die Renormierung zu berücksichtigen.
  • Für $O_2^{(c)}$ gab es qualitative Abweichungen vom VSA-Modell, was auf eine stärkere Unterdrückung durch axiale Beiträge hindeutet.
• Resonanzen: Die Methode konnte die Struktur der spektralen Dichte in der Nähe der $J/\psi$- und $\psi(2S)$-Resonanzen sichtbar machen.
• Herausforderung $\epsilon \to 0$: Die Extrapolation auf den physikalischen Grenzwert $\epsilon \to 0$ (um die scharfe Resonanzstruktur zu erhalten) erwies sich als schwierig, da die verfügbare Statistik noch nicht ausreicht, um in den asymptotischen Bereich ($\epsilon \ll \Delta(E)$) vorzudringen, wo eine polynomielle Extrapolation sicher ist.

5. Bedeutung und Ausblick

• Reduktion von Unsicherheiten: Die vorgestellte Methode ermöglicht es, die größten theoretischen Unsicherheiten in der Vorhersage von FCNC-Zerfällen (die Charming-Penguins) aus ersten Prinzipien zu bestimmen, anstatt sie zu modellieren.
• Test des Standardmodells: Präzise Berechnungen dieser Zerfälle sind entscheidend, um die aktuellen Spannungen zwischen Experiment (LHCb) und Theorie zu klären und neue Physik zu identifizieren.
• Zukunft: Das Paper legt den Grundstein für zukünftige, hochpräzise Berechnungen. Die Autoren planen, die Statistik zu erhöhen, die Renormierung vollständig durchzuführen und die Extrapolation $\epsilon \to 0$ mit hybriden Strategien (Kombination von Gitterdaten und Resonanzmodellen) zu meistern.

Fazit: Das Paper stellt einen Meilenstein dar, da es zeigt, dass komplexe hadronische Beiträge in seltenen B-Zerfällen, die bisher als "nicht berechenbar" galten, nun mit Gitter-QCD und spektralen Methoden zugänglich gemacht werden können.