Reassessing the gallium anomaly using self-consistent electron wave functions

Diese Arbeit bewertet das langjährige Gallium-Anomalie neu, indem sie die Neutrino-Einfangquerschnitte unter Verwendung selbstkonsistenter Elektronenwellenfunktionen berechnet, die aus der Dirac-Coulomb-Gleichung abgeleitet wurden, und dadurch die globale Signifikanz der Diskrepanz sowie deren Interpretation hinsichtlich steriler Neutrinos aktualisiert.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen zu zählen, wie viele Regentropfen in einen bestimmten Eimer fallen. Sie haben eine sehr präzise mathematische Formel, die genau vorhersagt, wie viele Tropfen basierend auf der Größe des Regens und der Größe des Eimers fallen sollten. Doch jedes Mal, wenn Sie die Tropfen in der Realität tatsächlich zählen, finden Sie weniger Tropfen, als Ihre Formel vorhergesagt hat. Dieser fehlende Regen ist das, was Physiker die „Gallium-Anomalie“ nennen.

Seit über 30 Jahren zeigen Experimente, die Gallium-71 (eine Art von Metall) als „Eimer“ verwenden, um Neutrinos (winzige, geisterhafte Teilchen aus der Sonne oder radioaktiven Quellen) aufzufangen, konsistent weniger Neutrinos als erwartet. Die Lücke zwischen der Vorhersage und der Realität ist so groß geworden, dass sie nun als ein großes Rätsel der Physik gilt.

Diese Arbeit von Cadeddu und Kollegen ist wie ein Team von Mechatronikern, das beschließt, den Motor dieser Vorhersageformel von Grund auf neu aufzubauen, um zu sehen, ob sie einen Rechenfehler gemacht haben.

Der alte Weg vs. der neue Weg

Der alte Motor (Die Näherung):
Zuvor berechneten Wissenschaftler, wie Neutrinos mit den Gallium-Atomen interagieren, unter Verwendung einer „Rohfassung“ der Mathematik. Sie behandelten die Elektronen (die winzigen Teilchen, die den Atomkern umkreisen) so, als wären sie einfache, glatte Wellen, die sich innerhalb des Atoms nicht viel verändern würden. Es war, als würde man die Form einer holprigen Straße abschätzen, indem man nur eine flache Karte betrachtet. Sie nahmen an, dass die Elektronenwelle überall im winzigen Kern gleich war.

Der neue Motor (Die exakte Lösung):
In dieser neuen Studie entschieden sich die Autoren dazu, die flache Karte beiseite zu legen. Stattdessen verwendeten sie ein hochpräzises GPS, um die exakten Gleichungen (die Dirac-Coulomb-Gleichung) zu lösen, die beschreiben, wie sich Elektronen tatsächlich verhalten.

• Die Analogie: Stellen Sie sich vor, der Kern ist eine belebte Tanzfläche. Die alte Methode ging davon aus, dass alle auf der Tanzfläche stillstehen, in einem perfekten Kreis. Die neue Methode zählt tatsächlich jeden einzelnen Tänzer und berücksichtigt dabei, wie sie gegeneinanderstoßen und sich im dichten Raum bewegen. Sie lösten die Mathematik sowohl für die Elektronen, die am Atom festgebunden sind, als auch für diejenigen, die herausfliegen, unter Verwendung eines spezialisierten Computerprogramms, um die exakte Form der Elektronen-„Welle“ zu erhalten.

Der „Averaging“-Trick (Der Mittelwert-Trick)

Ein weiterer entscheidender Unterschied in dieser Arbeit ist die Art und Weise, wie sie die Größe des Kerns handhaben.

• Der alte Weg: Sie wählten einen einzigen Punkt im Zentrum des Kerns (wie die Temperatur eines Raumes zu messen, indem man ein Thermometer genau in die Mitte steckt) und nahmen an, dass dieser den Rest repräsentiert.
• Der neue Weg: Sie erkannten, dass der Kern eine Ausdehnung hat, und daher „mittelten“ sie das Verhalten der Elektronen über das gesamte Volumen des Kerns. Es ist, als würde man die Temperatur an jedem Ort im Raum messen und den wahren Durchschnitt ermitteln, anstatt nur basierend auf der Mitte zu raten.

Was haben sie herausgefunden?

Als sie ihre neuen, präziseren Berechnungen durchführten:

1. Die Vorhersage änderte sich: Ihre neue, genauere Formel sagte voraus, dass weniger Neutrinos aufgefangen werden sollten als die alte Formel es tat.
2. Die Lücke vergrößerte sich: Da ihre neue Vorhersage niedriger ist, wurde der Unterschied zwischen dem, was wir zu sehen erwarten, und dem, was die Experimente tatsächlich sahen, noch größer.
3. Das Ergebnis: Die „fehlenden Neutrinos“ sind nun ein 5,5-Sigma-Problem. In der Wissenschaft ist „Sigma“ ein Maß für das Vertrauen. Ein 5-Sigma-Ergebnis ist der Goldstandard für eine Entdeckung, was bedeutet, dass die Wahrscheinlichkeit, dass es sich bei dieser Diskrepanz nur um einen Zufallsfund handelt, bei weniger als 1 zu 3,5 Millionen liegt.

Die „Sterile Neutrino“-Hypothese

Physiker haben eine Lieblingstheorie, um diesen fehlenden Regen zu erklären: Sterile Neutrinos.

• Die Metapher: Stellen Sie sich vor, die Neutrinos sind wie ein Vogelschwarm, der auf den Eimer zufliegt. Die Theorie besagt, dass einige dieser Vögel „unsichtbar“ (steril) sind. Sie interagieren überhaupt nicht mit dem Eimer; sie fliegen einfach direkt hindurch. Wenn diese unsichtbaren Vögel existieren, könnten sie erklären, warum der Eimer leerer ist als erwartet.

Die Autoren aktualisierten die Mathematik, um zu sehen, ob diese „unsichtbare Vogel“-Theorie immer noch passt.

• Die gute Nachricht: Die Mathematik lässt weiterhin zu, dass diese unsichtbaren Vögel existieren können. Die Gallium-Daten deuten nach wie vor stark auf deren Präsenz hin.
• Die schlechte Nachricht: Andere Experimente (die auf Reaktorneutrinos, Sonnenneutrinos und Massenmessungen blicken) haben „Zäune“ errichtet, die besagen, dass diese unsichtbaren Vögel nicht in der Weise fliegen können, wie es die Gallium-Daten nahelegen. Die Gallium-Daten wollen, dass die Vögel sehr aktiv sind, aber die anderen Zäune sagen, dass sie sehr scheu sein müssen.

Das Fazit

Die Autoren haben das Rätsel nicht gelöst; sie haben es eigentlich noch mysteriöser gemacht. Durch die Verwendung besserer Mathematik und präziserer Elektronenmodelle haben sie bestätigt, dass die „fehlenden Neutrinos“ ein reales, robustes Problem sind und kein Rechenfehler.

Sie kommen zu dem Schluss, dass die Idee des „Sterilen Neutrinos“ zwar nach wie vor der Hauptverdächtige ist, diese jedoch derzeit in einem Patt mit anderen experimentellen Belegen steht. Das Rätsel bleibt ungelöst, und die Autoren schlagen vor, dass ein neues Experiment mit einem anderen Typ von Detektor erforderlich sein könnte, um den Übeltäter endlich zu fassen.

Technische Zusammenfassung

Technische Zusammenfassung: Neubewertung der Gallium-Anomalie unter Verwendung selbstkonsistenter Elektron-Wellenfunktionen

Problemstellung
Die vorliegende Arbeit befasst sich mit der „Gallium-Anomalie“, einer anhaltenden Diskrepanz, die in Neutrino-Einfang-Experimenten unter Verwendung von $^{71}\text{Ga}$ als Target beobachtet wurde. Experimente mit radioaktiven Quellen wie $^{51}\text{Cr}$ und $^{37}\text{Ar}$ (GALLEX, SAGE und jüngst BEST) haben konsistent Neutrino-induzierte inverse Betazerfall-Raten (IBD) gemessen, die systematisch niedriger sind als die theoretischen Vorhersagen. Dieses Defizit, das nun eine Signifikanz von über $4\sigma$ erreicht, stellt die Robustheit theoretischer Wirkungsquerschnittsberechnungen infrage und hat Hypothesen über Physik jenseits des Standardmodells, wie etwa sterile Neutrinos, befeuert. Die Autoren identifizieren, dass frühere theoretische Auswertungen auf Näherungen der führenden Ordnung (LO) für elektronische Wellenfunktionen und inkonsistente Behandlungen von nuklearen und leptonischen Matrixelementen beruhten.

Methodik
Die Autoren unternehmen eine umfassende Revision des theoretischen Neutrino-Einfang-Wirkungsquerschnitts auf $^{71}\text{Ga}$, indem sie konventionelle LO-Näherungen zugunsten eines selbstkonsistenten numerischen Ansatzes aufgeben. Die wesentlichen methodischen Komponenten umfassen:

1. Exakte Lösungen der Wellenfunktion: Anstatt näherungsweise Wellenfunktionen zu verwenden, lösen die Autoren die Dirac-Coulomb-Gleichung sowohl für gebundene (Elektronen-Capture, EC) als auch für Kontinuums-Elektronenzustände (IBD) numerisch. Dies wird durch ein maßgeschneidertes Softwarepaket auf Basis der RADIAL-Bibliothek erreicht, welches die Dirac-Hartree-Fock-Slater-Gleichungen (DHFS) verwendet.
2. Selbstkonsistente Potentiale: Die Berechnungen nutzen ein DHFS-Potential, das nukleare, elektronische und Austausch-Potentiale einschließt. Das nukleare Potential wird mittels einer Zwei-Parameter-Fermi-Verteilung modelliert, welche die erste direkte Messung des $^{71}\text{Ge}$-Nuklearladungsradius ($R_{ch} = 4,032 \pm 0,002$ fm) integriert. Der Austauschpotential-Parameter ($C_{ex}$) wurde auf einen Mittelwert von $1,25$ festgelegt, um Modellvariationen zwischen den Slater- ($3/2$) und Kohn-Sham- ($1$) Approximationen zu berücksichtigen.
3. Mittelungsverfahren: Über die Konvention hinaus, die Wellenfunktionen an einem einzelnen Punkt auszuwerten (z. B. $r_0 = 0$ oder $r_0 = R_{box}$), schlagen die Autoren vor, die Elektron-Wellenfunktionen und die Fermi-Funktion über die nukleare Ladungsdichte ($\rho_{ch}(r)$) zu mitteln. Dies wird sowohl auf die IBD-Wirkungsquerschnittsberechnung als auch auf die Elektronen-Capture-Rate angewandt, die zur Extraktion des nuklearen Matrixelements mittels des Prinzip der detaillierten Bilanz verwendet wird.
4. Aktualisierte Eingabewerte: Die Analyse bezieht aktuelle experimentelle Daten ein, darunter:
   • Einen gewichteten Mittelwert aus vier unabhängigen Messungen der $^{71}\text{Ge}$ EC-Halbwertszeit ($t_{1/2} = 11,465 \pm 0,003$ d).
   • Einen aktualisierten Q-Wert für den EC-Prozess ($232,47 \pm 0,09$ keV).
   • Verfeinerte Wahrscheinlichkeitsverhältnisse für den Elektronen-Capture ($P_L/P_K$ und $P_M/P_K$) sowie Bindungsenergien.
   • Gamow-Teller-Stärken für angeregte Zustände, die aus $(p,n)$- und $(^3\text{He}, ^3\text{H})$-Streudaten abgeleitet wurden.

Wesentliche Beiträge und Ergebnisse

• Revidierte Wirkungsquerschnitte: Die neue Berechnung liefert einen Grundzustands-Wirkungsquerschnitt, der aufgrund des Mittelungsverfahrens etwa 2,8 % niedriger ist als bisherige Schätzungen. Unter Berücksichtigung der Beiträge angeregter Zustände sind die Gesamt-Wirkungsquerschnitte für $^{51}\text{Cr}$- und $^{37}\text{Ar}$-Quellen etwa 2 % kleiner als die in der jüngsten Literatur (Ref. [27]) berichteten Werte.
• Signifikanz der Anomalie: Durch den Vergleich der aktualisierten theoretischen Vorhersagen mit den experimentellen Daten von GALLEX, SAGE und BEST bewerten die Autoren die Signifikanz der Anomalie neu:
  • Unter Verwendung von $(p,n)$-Streudaten für die Beiträge angeregter Zustände erreicht die Diskrepanz $5,5\sigma$ ($R = 0,835^{+0,030}_{-0,048}$).
  • Unter Verwendung von $(^3\text{He}, ^3\text{H})$-Streudaten beträgt die Diskrepanz $4,7\sigma$ ($R = 0,811^{+0,037}_{-0,035}$).
  • Eine konservative Schätzung unter Verwendung nur des Grundzustands-Wirkungsquerschnitts ergibt eine $3,8\sigma$ Anomalie.
• Selbstkonsistenz: Die Autoren betonen, dass durch die Verwendung derselben Software und der gleichen zugrunde liegenden Annahmen für sowohl die Fermi-Funktion (IBD) als auch die Elektronendichte am Kern (EC) systematische Verschiebungen im Verhältnis kompensiert werden, was das Ergebnis robuster macht.

Bedeutung und Interpretation
Das Paper behauptet, dass die Gallium-Anomalie nun signifikanter ist als zuvor berichtet, mit Werten von bis zu $5,5\sigma$. Die Autoren untersuchen die Interpretation dieser Anomalie im Hinblick auf kurzreichweitige aktive-sterile Neutrino-Oszillationen (3+1-Modell) erneut. Ihre Ergebnisse bestätigen, dass die Gallium-Daten, einschließlich der präzisen BEST-Messungen, einen großen aktiven-sterilen Mischungswinkel erfordern ($0,16 \lesssim \sin^2 2\vartheta_{ee} \lesssim 0,54$ bei $2\sigma$ für $\Delta m^2_{41} \gtrsim 0,98 \text{ eV}^2$).

Die Autoren stellen jedoch fest, dass dieser erforderliche große Mischungswinkel im Konflikt mit den Einschränkungen durch Reaktor-Antineutrino-Daten, solare Neutrino-Grenzen und die KATRIN-Neutrinomasse-Messung steht. Infolgedessen kommt das Paper zu dem Schluss, dass die sterile Neutrino-Erklärung zwar die statistische Signifikanz der Anomalie erhöht, aber weiterhin in Spannung mit anderen experimentellen Daten steht. Die Autoren erklären, dass die Lösung der Gallium-Anomalie ein offenes Problem bleibt, das potenziell ein neues Quelleneperiment mit einem anderen Detektor erfordert, um die Diskrepanz aufzuklären.