Impact of ab initio E1 strength on neutron capture rates along the nickel isotopic chain

Tief im Herzen der Sterne ist eine ständige kosmische Alchemie am Werk, die die schweren Elemente schmiedet, aus denen unsere Welt besteht. Während leichtere Elemente wie Kohlenstoff und Sauerstoff in der stetigen Verbrennung stellarer Öfen entstehen, werden die schwersten Elemente, wie Gold und Platin, in gewaltsamen, extremen Ereignissen erschaffen, bei denen Neutronen mit rasender Geschwindigkeit eingefangen werden. Dieser Prozess, bekannt als der schnelle Neutroneneinfangprozess oder r-Prozess, ist die primäre Fabrik des Universums für diese schweren Materialien. Um genau zu verstehen, wie viel von jedem Element produziert wird, müssen Wissenschaftler die präzisen Raten kennen, mit denen Atomkerne Neutronen einfangen. Viele der instabilen, neutronenreichen Kerne, die an diesem kosmischen Tanz beteiligt sind, existieren jedoch nur für Bruchteile einer Sekunde, was es unmöglich macht, sie direkt im Labor zu untersuchen. Folglich müssen sich Forscher auf theoretische Berechnungen verlassen, um das Verhalten dieser flüchtigen Atome vorherzusagen und die Lücken zu füllen, wo direkte Messungen unmöglich sind.

Ein entscheidendes Puzzleteil dieses Rätsels ist, wie diese Kerne mit Licht interagieren, insbesondere wie sie Energie durch eine Art Schwingung absorbieren, die als elektrischer Dipolstärke bezeichnet wird. Vereinfacht ausgedrückt beschreibt dies, wie die Protonen und Neutronen innerhalb eines Atoms relativ zueinander schwingen, wenn sie von Energie getroffen werden. Jahrzehntelang haben Wissenschaftler breite, generalisierte Modelle verwendet, um dieses Verhalten für instabile Atome zu vermuten, doch diese Vermutungen lassen oft die Präzision vermissen, die für eine genaue Simulation der Entstehung von Elementen im Universum erforderlich ist. Eine neue Studie von Francesca Bonaiti und Joanna E. Sobczyk geht dieser Herausforderung direkt entgegen, indem sie diese Wechselwirkungen von Grund auf berechnet, wobei sie eine Methode verwenden, die mit den fundamentalen Kräften zwischen Teilchen beginnt, anstatt sich auf gebildete Vermutungen zu verlassen.

Die Forscher konzentrierten sich auf eine spezifische Familie von Atomen: die Nickel-Isotope. Nickel spielt eine Schlüsselrolle im r-Prozess, da es als Trittstein fungiert, der die Häufigkeit von Elementen mit Massen um den Wert einhundert beeinflusst. Das Team nutzte einen leistungsstarken computergestützten Ansatz namens „ab initio“, was „von Anfang an“ bedeutet. Anstatt ihre Berechnungen an bestehende Daten anzupassen, begannen sie mit den grundlegenden Gesetzen der Kernphysik, insbesondere den Wechselwirkungen zwischen Protonen und Neutronen, die aus einer Theorie namens chiraler effektiver Feldtheorie abgeleitet sind. Sie wandten dies auf eine Kette von Nickelatomen an, die von dem stabilen Nickel-58 bis zum extrem neutronenreichen Nickel-80 reicht. Um die komplexe Mathematik dieser instabilen Open-Shell-Kerne zu bewältigen, entwickelten sie eine neue Erweiterung einer Methode, die als Lorentz-Integral-Transformations-Cluster-Methode bekannt ist. Dies ermöglichte es ihnen, die Reaktion dieser Atome auf Energie zu simulieren, ohne die unmöglichen Streuzustände konstruieren zu müssen, die solche Berechnungen normalerweise verkomplizieren.

Die Ergebnisse zeigten ein signifikantes Merkmal im Verhalten dieser neutronenreichen Atome: eine ausgeprägte, niedrig liegende Stärke knapp oberhalb des Punktes, an dem ein Neutron den Kern verlassen kann. Dieses Phänomen, das oft mit einer „pygmalischen Dipolresonanz“ assoziiert wird, wirkt wie eine weiche, niederenergetische Schwingung, die deutlicher hervortritt, wenn der Kern mehr Neutronen gewinnt. Im Fall von Nickel-80 ist diese niederenergetische Schwingung so stark und liegt so nah an der Entschwellenenergie, dass sie beginnt, das Verhalten des Atoms bereits bei überraschend niedrigen Temperaturen von etwa 0,5 Milliarden Grad zu beeinflussen. Für die anderen Isotope in der Kette verstärkt sich dieser Effekt mit zunehmender Neutronenzahl und wird zu einem wesentlichen Faktor bei Temperaturen zwischen 3 und 10 Milliarden Grad.

Als die Forscher diese neuen Berechnungen aus erster Hand in die Standard-Computercodes einspeisten, die zur Simulation stellarer Explosionen verwendet werden, war der Einfluss deutlich. Die Einbeziehung dieser spezifischen elektrischen Dipolstärke steigerte die vorhergesagten Raten, mit denen diese Nickelkerne Neutronen einfangen, erheblich. Bei den schwereren Isotopen stieg die Einfangrate je nach Temperatur und spezifischem Isotop um Faktoren zwischen 16 Prozent und dem Dreifachen der bisherigen Schätzungen. Diese Steigerung ist keine geringfügige Anpassung; sie stellt eine fundamentale Verschiebung in unserem Verständnis darüber dar, wie diese Kerne mit ihrer Umgebung interagieren. Die Studie hob auch hervor, dass die größte Unsicherheit in ihren Ergebnissen aus dem mathematischen Prozess der Umkehrung ihrer Berechnungen zur Gewinnung des Endergebnisses resultiert und nicht aus den zugrunde liegenden Physikmodellen selbst.

Durch die Bereitstellung der ersten rigorosen Beschreibung dieser elektrischen Dipolstärke entlang einer Isotopenkette aus erster Hand bietet die Arbeit eine neue, zuverlässigere Grundlage für astrophysikalische Simulationen. Sie legt nahe, dass die Produktion schwerer Elemente im Universum effizienter sein könnte als bisher angenommen, insbesondere für die schwersten Nickelisotope. Während die Forscher einräumen, dass zukünftige Arbeiten den mathematischen Inversionsprozess verfeinern müssen, um die Unsicherheiten zu verringern, schlägt diese Studie erfolgreich die Brücke zwischen fundamentaler Kernphysik und der grandiosen Skala der kosmischen Elemententstehung. Sie zeigt, dass wir selbst für die instabilsten Atome nun das Verhalten mit einem Detailgrad vorhersagen können, der zuvor unerreichbar war, und uns der Lösung des Rätsels näher bringen, woher die schweren Elemente in unserem Universum wirklich stammen.