Theoretical analysis and predictions for the double electron capture of $^{124}$Xe

Diese Arbeit präsentiert eine umfassende theoretische Analyse des Zwei-Neutrino-Doppel-Elektronen-Einfangs in $^{124}$Xe durch die Verbesserung der Kern- und Atomstruktur-Berechnungen, was verfeinerte Kernmatrixelemente liefert, spezifische Einfang-Bruchteile für verschiedene Zerfallskanäle vorhersagt (insbesondere 74 % für den KK-Kanal und 24 % für die kumulativen KL$_1$-KO$_1$-Kanäle) und aktualisierte atomare Relaxationsenergien für die Hintergrundmodellierung in Flüssigxenon-Experimenten bereitstellt.

Das Wesentliche

Stellen Sie sich das Atom als eine winzige, geschäftige Stadt vor. Im Inneren dieser Stadt ist das Zentrum, das Xenon-124, das Rathaus, das mit Protonen und Neutronen vollgepackt ist. Normalerweise sind diese Bürger sehr stabil, aber manchmal beschließen sie, sich neu zu ordnen, um es sich bequemer zu machen.

In diesem Ereignis beschließt das Rathaus, zwei seiner eigenen Bewohner (Elektronen) aus den äußeren Vierteln einzufangen und sie in das Zentrum zu ziehen. Wenn dies geschieht, verwandelt sich die Stadt in eine neue Stadt namens Tellur-124 und stößt zwei winzige, unsichtbare Boten namens Neutrinos aus.

Wissenschaftler nennen dies den Doppel-Elektroneneinfang (speziell die Version mit zwei Neutrinos, oder $2\nu$ECEC). Es ist wie ein doppeltes Eintauchen in einen Pool, aber anstelle von Wasser sind es subatomare Teilchen.

Hier ist das, was die Forscher getan haben, einfach erklärt:

1. Ein besseres Bauplan erstellen (Die Theorie)

In der Vergangenheit versuchten Wissenschaftler vorherzusagen, wie oft dieser „Doppel-Dip“ vorkommt, aber ihre Baupläne waren etwas grob. Sie übersahen einige Details darüber, wie sich die Elektronen bewegen und wie der Kern reagiert.

Die Autoren dieser Arbeit beschlossen, einen viel präziseren Bauplan zu erstellen.

• Die „Taylor-Entwicklung“-Analogie: Stellen Sie sich vor, Sie versuchen, den Pfad eines Autos zu beschreiben. Eine einfache Beschreibung könnte nur sagen: „Es fährt vorwärts.“ Eine bessere Beschreibung fügt hinzu: „Es beschleunigt.“ Die beste Beschreibung fügt hinzu: „Es beschleunigt, biegt dann leicht ab, dann wird es langsamer.“ Die Autoren verwendeten ein mathematisches Werkzeug namens „Taylor-Entwicklung“, um diese zusätzlichen Details (bis zur vierten Potenz der Energie) hinzuzufügen. Dies ermöglichte es ihnen, die „Kurven und Verlangsamungen“ des Zerfallsprozesses zu sehen, die frühere Modelle übersehen hatten.
• Die „neuen Verhältnisse“: Da sie diese zusätzlichen Details hinzufügten, entdeckten sie neue Wege, um verschiedene Teile des Prozesses zu vergleichen (die sogenannten $\xi$-Verhältnisse). Betrachten Sie dies als neue Kontrollpunkte auf einer Rennstrecke, die Wissenschaftlern helfen können, das Rennen später besser zu verstehen.

2. Das Viertel betrachten (Der atomare Teil)

Um zu berechnen, wie wahrscheinlich dieses Ereignis ist, muss man genau wissen, wo die Elektronen leben.

• Die „Pauli-Blockade“-Metapher: Stellen Sie sich einen überfüllten Aufzug vor. Wenn der Aufzug voll ist, können Sie nicht einfach jemanden hineindrücken; Sie müssen warten, bis jemand aussteigt. Im Kern sind die „innersten“ Plätze wie ein voller Aufzug. Die Autoren erkannten, dass die eingefangenen Elektronen nicht einfach irgendwohin gehen können; sie werden durch die bereits dort befindlichen Elektronen blockiert. Sie berücksichtigten diese „Überfüllung“, was die Berechnung verändert.
• Die Suche ausweiten: Frühere Studien betrachteten nur die zwei nächsten Viertel zum Kern (die sogenannten K- und L1-Schalen). Die Autoren sagten: „Lassen Sie uns alle Viertel betrachten, sogar die, die weiter draußen liegen (bis zur O-Schale).“ Sie fanden heraus, dass die äußeren Viertel zwar weniger wahrscheinlich eingefangen werden, aber dennoch zu dem Gesamtereignis beitragen.

3. Das Rathaus simulieren (Der nukleare Teil)

Der Kern ist der schwierigste Teil zu simulieren, da er ein chaotischer Haufen von Teilchen ist. Die Autoren verwendeten zwei verschiedene „Simulations-Engines“, um vorherzusagen, wie der Kern sich verhält:

• Engine A (ISM): Dies ist wie eine detaillierte Raum-für-Raum-Simulation. Sie ließen die Simulation mit unterschiedlichen Regeln (den sogenannten „Hamiltonianen“) laufen, um zu sehen, ob die Ergebnisse Bestand haben. Sie fanden heraus, dass sie, wenn sie alle möglichen „Zwischenschritte“ einbezogen, die die Menge der „Stärke“ des Ereignisses niedriger vorhersagten, als es ältere, einfachere Modelle suggerierten.
• Engine B (pn-QRPA): Dies ist eine andere Art von Simulation. Sie passten die Einstellungen dieser Engine an, bis sie mit den realen Daten übereinstimmten, die wir bereits haben. Sie fanden heraus, dass ihre neue, sorgfältigere Berechnung einen viel kleineren „Stärkewert“ ergab als frühere Versuche unter Verwendung dieser Engine.

4. Die Ergebnisse: Was haben sie herausgefunden?

Durch die Kombination ihres besseren Bauplans, ihrer detaillierten Viertel-Karte und ihrer zwei Simulations-Engines trafen sie mehrere Vorhersagen:

• Das Hauptereignis (KK-Kanal): Sie sagen voraus, dass in etwa 74 % der Fälle die zwei Elektronen aus dem allernächsten Viertel (der K-Schale) eingefangen werden. Dies ist eine kleine, aber wichtige Anpassung gegenüber den 72,4 %, die in früheren Experimenten verwendet wurden.
• Die „nächstbesten“ Ereignisse: Sie sagen voraus, dass in etwa 19 % der Fälle ein Elektron aus dem nächsten Viertel und das andere aus dem nächstfolgenden (KL1) kommt.
• Die „kumulative“ Vorhersage: Wenn man alle etwas selteneren Ereignisse zusammenzählt (von KL1 bis KO1), machen diese etwa 24 % des Gesamtereignisses aus. Dies ist etwa ein Drittel des Hauptereignisses.
• Die „Relaxationsenergie“: Wenn die Elektronen eingefangen werden, ist die neue Stadt (Tellur) angeregt und muss zur Ruhe kommen. Dies geschieht durch die Freisetzung von Energie (wie Röntgenstrahlen). Die Autoren berechneten genau, wie viel Energie für jeden Typ des Einfangs freigesetzt wird. Dies ist, als würde man Wissenschaftlern einen spezifischen „Fingerabdruck“ der Energie geben, nach dem sie in ihren Detektoren suchen können.

Warum ist das wichtig?

Das Papier behauptet nicht, Krankheiten zu heilen oder Städte mit Energie zu versorgen. Stattdessen fungiert es als verfeinerte Karte für Entdecker.

Große Experimente, die flüssiges Xenon verwenden (wie jene, die nach Dunkler Materie suchen), beobachten ständig dieses spezifische „Doppel-Dip“-Ereignis. Dieses Ereignis kann jedoch wie „Hintergrundrauschen“ aussehen, das die Daten verwirrt. Indem sie eine genauere Karte darüber liefern, wie oft dies geschieht, welche Energie es freisetzt und aus welchen „Vierteln“ die Elektronen stammen, helfen die Autoren den Experimentalisten, zwischen einem echten Signal und Hintergrundrauschen zu unterscheiden.

Kurz gesagt: Sie haben ein unscharfes, niedrig auflösendes Foto eines seltenen atomaren Ereignisses genommen und es in ein hochauflösendes 3D-Modell verwandelt, damit Wissenschaftler genau wissen, wonach sie in ihren Detektoren suchen müssen.

Technische Zusammenfassung

Technische Zusammenfassung: Theoretische Analyse und Vorhersagen für den doppelten Elektroneneinfang von $^{124}$Xe

Problemstellung
Der Prozess des doppelten Elektroneneinfangs (ECEC), insbesondere der Zwei-Neutrino-Modus ($2\nu$ECEC), ist ein seltener Kernübergang von großer Bedeutung für die Untersuchung der Kernstruktur und der Neutrinoeigenschaften. Während der $2\nu\beta^-\beta^-\text{-Zerfall}$ intensiv untersucht wurde, erhielten ECEC-Modi aufgrund geringerer Übergangswahrscheinlichkeiten weniger Aufmerksamkeit. $^{124}$Xe stellt jedoch einen einzigartigen Fall dar, da es das größte Q-Wert unter den Isotopen besitzt, die zu ECEC-Übergängen fähig sind, was es zu einem primären Ziel für groß angelegte Flüssigxenon-Experimente (z. B. bei der Suche nach Dunkler Materie und neutrinolosem Doppelbetazerfall) macht.

Frühere theoretische Berechnungen der $2\nu$ECEC-Halbwertszeit von $^{124}$Xe litten unter Einschränkungen sowohl in der Kern- als auch in der Atommodellierung. Speziell verließen sich frühere Arbeiten häufig auf die Closure-Approximation für Kernmatrixelemente (NMEs), betrachteten nur K- und L$_1$-Schalen-Einfänge und verwendeten weniger rigorose Behandlungen der atomaren Abschirmung und des Pauli-Blockings. Da nun experimentelle Messungen der $^{124}$Xe-Halbwertszeit verfügbar sind, ist eine präzisere und vollständigere theoretische Beschreibung erforderlich, um diese Daten zu interpretieren und unentdeckte Zerfallskanäle vorherzusagen.

Methodik
Die Autoren stellen einen umfassenden theoretischen Rahmen bereit, der sowohl die Formalismen als auch die spezifischen Berechnungen für $^{124}$Xe verbessert:

1. Verbesserter Formalismus: Die Zerfallsrate wird mittels einer Taylor-Entwicklungsmethode abgeleitet, wobei Terme bis zur vierten Potenz der Leptonenenergien beibehalten werden. Dieser Ansatz führt höhere Ordnungen ein und definiert neue NME-Verhältnisse, $\xi_{31}^{2\nu\text{ECEC}}$ und $\xi_{51}^{2\nu\text{ECEC}}$, die von den Energiedenominatoren der intermediären $1^+$-Zustände abhängen.
2. Kernstruktur-Berechnungen:
   • Interacting Shell Model (ISM): Die Berechnungen wurden im $jj55$-Modellraum (Orbitale $0g_{7/2}, 1d_{5/2}, 1d_{3/2}, 2s_{1/2}, 0h_{11/2}$) durchgeführt. Im Gegensatz zu früheren Studien, die die Single-State-Dominance (SSD) annahmen, führt diese Arbeit eine vollständige Summation über die intermediären $1^+$-Zustände mittels eines Stärkefunktionsansatzes durch. Es wurden zwei effektive Hamiltonoper (GCN5082 und SVD) mit Quenching-Faktoren ($q_H$) verwendet, die auf die $2\nu\beta^-\beta^-\text{-Daten}$ für $^{136}$Xe kalibriert wurden. Trunkierungen, die bis zu vier Nukleonen-Anregungen aus niedrigeren Orbitalen zulassen, wurden einbezogen.
   • Proton-Neutron-Quasipartikel-Random-Phase-Approximation (pn-QRPA): Die Berechnungen nutzten die Isospin-Restaurierung. Die Teilchen-Teilchen-Wechselwirkungsstärke ($g_{pp}^{T=0}$) wurde angepasst, um die experimentelle Halbwertszeit von $^{124}$Xe zu reproduzieren. Die Studie verwendete realistische Nukleon-Nukleon-Potentiale (Argonne V18 und CD-Bonn) und berücksichtigte den BCS-Vakuum-Überlapp.
3. Atomstruktur und Phasenraumfaktoren (PSFs):
   • Das Dirac-Hartree-Fock-Slater (DHFS)-Selbstkonsistenz-Framework wurde verwendet, um Elektronenwellenfunktionen zu erzeugen.
   • Verbesserungen umfassen die rigorose Behandlung der atomaren Abschirmung, diffuse Kernoberflächenkorrekturen, eine realistische Kernladungsdichte und Elektronen-Austausch-Korrelations-Effekte.
   • Entscheidend ist, dass die Berechnung das Pauli-Blocking der innersten Nukleonzustände berücksichtigt und die Summation über den Einfang auf alle verfügbaren s-Wellen-Elektronen (bis zur O$_1$-Subschale) ausweitet, statt sich wie in früheren Studien auf K- und L$_1$-Schalen zu beschränken.
   • Atomare Relaxationsenergien für das endgültige $^{124}$Te-Atom wurden berechnet, um die Hintergrundmodellierung zu unterstützen.

Wichtigste Ergebnisse

• Kernmatrixelemente (NMEs):
  • Die ISM-Berechnungen liefern NMEs, die gegenüber Änderungen der Trunkierung und des effektiven Hamiltonoperators robust sind. Die vollständige Summation über die intermediären Zustände (ohne SSD) reduzierte das NME für den GCN5082-Hamiltonoperator im Vergleich zu trunkierten Ansätzen um etwa 25 %.
  • Die pn-QRPA-Ergebnisse, wenn sie durch die experimentelle Halbwertszeit eingeschränkt sind, liefern NME-Werte, die signifikant niedriger sind als frühere pn-QRPA-Berechnungen (die oft auf unterschiedlichen Anpassungsstrategien für $g_{pp}$ basierten). Die neuen NMEs für $^{124}$Xe sind vergleichbar mit denen des $2\nu\beta^-\beta^-\text{-Zerfalls}$ von $^{128,130}$Te, was ähnliche Paarungseigenschaften widerspiegelt.
• Halbwertszeit-Vorhersagen:
  • Unter Verwendung des ISM stimmt die berechnete Gesamthalbwertszeit für $^{124}$Xe innerhalb eines Faktors von zwei mit dem experimentellen Wert von $(1,1 \pm 0,2_{\text{stat}} \pm 0,1_{\text{sys}}) \times 10^{22}$ J überein.
  • Die Einbeziehung von Taylor-Entwicklungstermen und die verfeinerten PSFs trugen dazu bei, die theoretischen Vorhersagen näher an die experimentellen Daten zu bringen.
• Einfang-Fraktionen und Kanäle:
  • Die vorhergesagte Einfang-Fraktion für den dominanten KK-Kanal beträgt 74,1 %, was von den 72,4 % abweicht, die in früheren experimentellen Analysen verwendet wurden. Diese Differenz ergibt sich aus der exklusiven Berücksichtigung von s-Wellen-Elektronen im theoretischen Modell gegenüber der Einbeziehung aller Elektronen im experimentellen Signalmodell.
  • Der nächstwahrscheinlichere Kanal ist KL$_1$ mit einer Fraktion von etwa 19 %.
  • Die kumulative Einfang-Fraktion der Kanäle KL$_1$ bis KO$_1$ wird auf etwa 24 % geschätzt.
• Atomare Relaxationsenergien:
  • Präzise Relaxationsenergien wurden für verschiedene Loch-Konfigurationen berechnet (z. B. KK bei 64,62 keV, KL$_1$ bei 37,05 keV). Diese Werte sind für die Verwendung in der Hintergrundmodellierung von Flüssigxenon-Experimenten vorgesehen.

Bedeutung und Ansprüche
Das Paper beansprucht, die erste vollständige theoretische Beschreibung des $2\nu$ECEC in $^{124}$Xe zu liefern, die gleichzeitig die Kernstruktur-Behandlung (durch vollständige Summation und Isospin-Restaurierung) und die Atomphysik (durch DHFS, Pauli-Blocking und erweiterte Orbital-Summation) verbessert.

Die Autoren behaupten, dass ihre ISM-Ergebnisse konsistent mit den experimentellen Daten sind, was die Zuverlässigkeit der berechneten NMEs und der neu vorhergesagten Einfang-Fraktionen validiert. Ein primärer Beitrag ist die Vorhersage, dass die kumulativen KL$_1$-KO$_1$-Kanäle (die etwa 24 % der Zerfälle repräsentieren) im Energieintervall von 31,93–37,05 keV experimentell beobachtbar sein sollten. Darüber hinaus legt die Arbeit nahe, dass das Verhältnis der partiellen Halbwertszeiten (z. B. $T_{KK}^{1/2}/T_{KL_1}^{1/2}$) als experimenteller Constraint für den $\xi_{31}^{2\nu\text{ECEC}}$-Parameter dienen könnte, was einen neuen Weg eröffnet, das Zusammenspiel zwischen niedrig liegenden und höher liegenden intermediären Zuständen zu untersuchen. Die Studie betont, dass frühere theoretische Schätzungen der Halbwertszeiten aufgrund der Verwendung von Closure-Approximationen und unvollständigen Orbital-Summationen oft inkompatibel mit den genauesten PSFs waren.