The highest tightly bound nuclei and the stability prediction of their isobaric nuclei versus beta decay in the liquid drop model

Technische Zusammenfassung: Die am stärksten gebundenen Kerne und die Stabilitätsprognose im Flüssigkeitstropfenmodell

Problemstellung
Die Studie befasst sich mit der anhaltenden Debatte über die Identifizierung des am stärksten gebundenen Kerns und der Genauigkeit theoretischer Modelle bei der Vorhersage von Kernmassen und der Beta-Zerfallsstabilität. Während in der Literatur häufig $^{56}\text{Fe}$ als der am stärksten gebundene Kern angeführt wird, deuten experimentelle Daten aus der AME2020-Kompilation darauf hin, dass $^{62}\text{Ni}$, $^{58}\text{Fe}$ und andere Kerne mittlerer Masse höhere durchschnittliche Bindungsenergien pro Nukleon aufweisen. Darüber hinaus untersuchen die Autoren die Einschränkungen der standardmäßigen fünfgliedrigen Bethe–Weizsäcker-Massenformel (BWMF-5T) und ob die Erweiterung des Modells durch zusätzliche Terme sowie die Optimierung der Koeffizienten durch regionale gegenüber globalen Anpassungsverfahren eine überlegene Vorhersagegenauigkeit für diese spezifischen Kerne liefert.

Methodik
Die Autoren verwenden einen semi-empirischen Ansatz basierend auf dem Flüssigkeitstropfenmodell (Liquid Drop Model, LDM). Die Methodik umfasst drei primäre Komponenten:

1. Modellformulierung: Die Studie nutzt die standardmäßige fünfgliedrige BWMF (Volumen, Oberfläche, Coulomb, Asymmetrie und Paarung) und erweitert diese zu einer zehngliedrigen Formulierung (BWMF-10T), indem fünf zusätzliche Terme integriert werden: der Wigner-Term, die Coulomb-Austauschkorrektur, die Oberflächen-Asymmetrie-Korrektur, der Krümmungsenergie-Term und ein Schalenkorrektur-Term.
2. Datenanpassung und Optimierung: Die Koeffizienten für die 5-Term- und 10-Term-Formeln werden mittels eines Kleinste-Quadrate-Verfahrens gegen die experimentellen Massendaten der AME2020-Kompilation bestimmt. Die Anpassung erfolgt über mehrere Datensätze, um die Sensitivität zu bewerten:
   • Ein globaler Satz von 2.457 gemessenen Kernen.
   • Sukzessive Teilmengen von jeweils 500 Kernen.
   • Ein regionaler Satz von 366 geradzahligen (even–even) Kernen im Massenbereich von Sauerstoff bis Gadolinium.
3. Stabilitätsanalyse: Unter Verwendung der optimierten Koeffizienten berechnen die Autoren die Bindungsenergie pro Nukleon für die sechs am stärksten gebundenen Kandidaten ($^{62}\text{Ni}$, $^{58}\text{Fe}$, $^{56}\text{Fe}$, $^{60}\text{Ni}$, $^{54}\text{Cr}$ und $^{64}\text{Ni}$). Zusätzlich analysieren sie die Beta-Zerfallsstabilität innerhalb isobarer Familien für die Massenzahlen $A = 54, 56, 58, 60, 62$ und $64$, indem sie die Massenparabel ableiten und das stabilste Isobar ($Z_0$) durch Differentiation des Massenausdrucks bestimmen.

Zentrale Beiträge

• Vergleichende Koeffizientenanalyse: Die Arbeit vergleicht systematisch die aus regionalen Anpassungen (Kerne mittlerer Masse) abgeleiteten Koeffizientenwerte mit globalen Anpassungen (gesamte Nuklidkarte). Sie zeigt, dass die Koeffizienten hochsensibel gegenüber der Größe und Diversität des zur Anpassung verwendeten Datensatzes sind.
• Bewertung des erweiterten Modells: Die Studie evaluiert die Auswirkungen der Hinzufügung von fünf spezifischen Termen zur Standard-Massenformel und analysiert, wie diese Terme die Symmetrieenergie-Koeffizienten und die allgemeine Anpassungsqualität beeinflussen.
• Vorhersage der Beta-Stabilität: Die Autoren wenden die verfeinerten Formeln an, um die Stabilitätstrends isobarer Ketten vorherzusagen, berechnen explizit die $Q$-Werte für $\beta^-$- und $\beta^+$-Übergänge und identifizieren die stabilsten Kerne durch Massenminimierung.

Ergebnisse

• Koeffizientensensitivität: Die angepassten Koeffizienten für die grundlegenden fünf Terme zeigen eine systematische Zunahme ihrer Größenordnung, wenn sich der Datensatz von 500 auf 2.457 Kerne erweitert. Beispielsweise steigt der Koeffizient der Symmetrieenergie signifikant an, wenn man von einer regionalen zu einer globalen Anpassung übergeht.
• Regionale vs. globale Genauigkeit: Regionale Anpassungen (speziell für den Bereich Sauerstoff–Gadolinium) übertreffen globale Anpassungen bei der Anwendung auf Kerne mittlerer Masse deutlich. Globale Anpassungen, die schwere und superschwere Kerne einschließen, führen zu systematischen Abweichungen bei der Vorhersage der Bindungsenergien leichterer, stark gebundener Kerne.
• Bindungsenergie-Vorhersagen:
  • Theoretische Berechnungen unter Verwendung regionaler Koeffizienten zeigen eine exzellente Übereinstimmung mit den experimentellen Bindungsenergien für $^{62}\text{Ni}$, $^{58}\text{Fe}$ und $^{64}\text{Ni}$.
  • Der relative Fehler ($\Delta B/B$) für die am stärksten gebundenen Kerne bleibt bei regionalen Anpassungen unter 0,26 %, während globale Anpassungen größere Abweichungen aufweisen.
  • Die Studie bestätigt $^{62}\text{Ni}$ als den Kern mit der experimentell höchsten beobachteten Bindungsenergie pro Nukleon unter den bekannten Nukliden, was konsistent mit den AME2020-Daten ist.
• Schaleneffekte: Die Analyse legt nahe, dass die hohe Bindungsenergie von $^{62}\text{Ni}$ mit seiner Protonenzahl ($Z=28$, eine magische Zahl) und der Neutronenkonfiguration zusammenhängt, während $^{58}\text{Fe}$ von einem geschlossenen Neutronen-Subshell-Zustand ($N=32$) profitiert.
• Beta-Stabilität: Die abgeleiteten Massenparabeln identifizieren korrekt die stabilen Isobare für die ausgewählten Massenketten, wobei der Paarungsterm die Energieunterschiede zwischen Odd-$Z$- und Even-$Z$-Übergängen erfolgreich berücksichtigt.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die primäre Bedeutung dieser Untersuchung darin liegt, dass die regionale Massenanpassung eine überlegene Beschreibung von Kernen mittlerer Masse bietet im Vergleich zu globalen Anpassungsverfahren. Die Autoren behaupten, dass die standardmäßige Bethe–Weizsäcker-Formel, wenn sie mit regionalen Koeffizienten verfeinert wird, eine hochgenaue Näherung der Atommassen für die am stärksten gebundenen Kerne darstellt.

Die Studie kommt zu dem Schluss, dass die Vorhersagegenauigkeit des Flüssigkeitstropfenmodells für spezifische Regionen des Kernen-Diagramms stark von dem zur Optimierung der Koeffizienten verwendeten Datensatz abhängt. Die Autoren merken bescheiden an, dass die Häufigkeit von $^{62}\text{Ni}$ in der Natur (wie etwa die Verbreitung von $^{56}\text{Fe}$) durch stellare Nukleosyntheseprozesse und nicht allein durch die Bindungsenergie bestimmt wird. Die Arbeit hebt die Notwendigkeit weiterer experimenteller Messungen leichter Nickel-Isotope hervor, um die Schalenlücken bei $Z=28$ zu klären, schlägt jedoch selbst keine neuen experimentellen Aufbauten vor.