Impact of neutron captures on the r-process: I. Propagation of statistical model uncertainties

Das Universum ist eine riesige chemische Fabrik, doch seit Jahrzehnten versuchen Wissenschaftler zu erklären, wie etwa die Hälfte aller Elemente schwerer als Eisen geschmiedet wurde. Diese schweren Elemente, zu denen alles von Gold bis Uran gehört, entstehen bei den gewaltigsten kosmischen Ereignissen, die man sich vorstellen kann, wie etwa der Kollision zweier Neutronensterne. In diesen extremen Umgebungen werden Atomkerne von einer Flut aus Neutronen bombardiert, was den schnellen Aufbau immer schwererer Atome in einem Prozess ermöglicht, der als schneller Neutroneneinfangprozess bekannt ist. Während Astronomen bestätigt haben, dass diese Kollisionen tatsächlich stattfinden und die richtigen Arten von Elementen erzeugen, bleibt die präzise Physik, die regelt, wie sich die Atome während der Explosion verändern, weitgehend ein Rätsel. Da wir diese Bedingungen in einem Labor nicht nachbilden können, müssen sich Forscher auf komplexe Computermodelle verlassen, um die Reise der Materie von einem einfachen Atom zu einem schweren Element zu simulieren. Die Genauigkeit dieser Simulationen hängt vollständig von den Daten ab, die in sie eingespeist werden – insbesondere davon, wie gut wir die internen Energiezustände instabiler, neutronenreicher Atome verstehen, die nur für einen flüchtigen Moment existieren.

Ein Forschungsteam unter der Leitung von Yu An an der Technischen Universität Darmstadt hat sich intensiv mit einem der kritischsten Teile dieses Puzzles beschäftigt: den statistischen Modellen, die vorhersagen, wie diese instabilen Atome Neutronen einfangen. In ihrer Studie versuchten sie nicht, neue Physik zu erfinden, sondern testeten vielmehr, wie stark sich die Endergebnisse einer kosmischen Explosion ändern, wenn die zugrunde liegenden Annahmen über die Atomstruktur leicht verändert werden. Sie konzentrierten sich auf zwei spezifische theoretische Bestandteile: die Kernniveaudichte, die beschreibt, auf wie viele verschiedene Arten ein Atomkern vibrieren oder Energie speichern kann, und die Gamma-Strahlungsstärke-Funktion, die misst, wie wahrscheinlich es ist, dass der Kern Lichtenergie emittiert, während er sich stabilisiert. Durch das Durchführen tausender Simulationen mit verschiedenen Kombinationen dieser theoretischen Modelle kartierte das Team genau, wo unser Wissensdefizit die größten Lücken in unserem Verständnis der chemischen Geschichte des Universums hinterlässt.

Die Forscher begannen damit, die Nachwirkungen einer Neutronensternverschmelzung zu simulieren – ein Szenario, bei dem zwei tote Sterne zusammenstoßen und eine Wolke aus superdichter Materie ausstoßen. Sie führten diese Simulation unter Verwendung von sechzig verschiedenen Kombinationen der ihnen zur Verfügung stehenden theoretischen Modelle durch, wodurch eine breite Palette möglicher Ergebnisse für die endgültige Mischung der Elemente entstand. Die Ergebnisse zeigten eine klare Hierarchie der Unsicherheit. Die Modelle, die beschreiben, wie Kerne Energie speichern (bekannt als Kernniveaudichten), hatten einen weita-gestaltig größeren Einfluss auf die endgültige chemische Zusammensetzung als die Modelle, die beschreiben, wie sie Licht emittieren. Während die großen Spitzen der schweren Elemente, wie sie beispielsweise Gold und Platin entsprechen, unabhängig vom verwendeten Modell relativ stabil blieben, zeigten die Regionen zwischen diesen Spitzen dramatische Unterschiede. In einigen Simulationen würde die Häufigkeit von Elementen in einem bestimmten Massenbereich um 172 bis 178 atomare Einheiten, je nachdem, welches Energiemodell gewählt wurde, entweder um das Dreifache ansteigen oder völlig verschwinden.

Um zu verstehen, warum dies geschieht, untersuchte das Team den zeitlichen Ablauf der Ereignisse innerhalb der expandierenden Trümmerwolke. Nach der anfänglichen Explosion sinkt die Temperatur rapide ab. Zuerst sind die Atome so heiß, dass Neutronen ständig eingefangen und dann wieder abgeschlagen werden, in einem balancierten Tanz. Schließlich fällt die Temperatur so weit ab, dass dieses Gleichgewicht bricht und die Atome beginnen, sich in ihre endgültigen Formen zu stabilisieren. Die Forscher fanden heraus, dass die Wahl des Energiemodells in der „kalten“ Phase dieses Prozesses entscheidend wird, die eintritt, wenn die Temperatur auf etwa 0,45 Milliarden Grad sinkt. In diesem Zeitfenster gilt: Wenn das theoretische Modell vorhersagt, dass die Atome Neutronen langsam einfangen, verweilt das Material in einer spezifischen Region des Periodensystems und erzeugt so eine ausgeprägte Sub-Spitze von Elementen. Wenn das Modell einen schnelleren Einfang vorhersagt, rast das Material vorbei, um noch schwerere Elemente zu werden, und hinterlässt eine Lücke. Diese Sensitivität ist nicht nur eine Laune einer einzelnen Simulation; das Team stellte fest, dass selbst bei der Mittelung der Ergebnisse über 11.000 verschiedene Explosionsszenarien hinweg, um eine realistische galaktische Mischung nachzubilden, diese spezifischen Lücken und Spitzen bestehen blieben.

Die Studie untersuchte auch ein weites Netz hypothetischer Bedingungen, indem sie die anfängliche Dichte und das Verhältnis von Protonen zu Elektronen im Material variierte, um zu sehen, wie verschiedene Umgebungen darauf reagieren. Sie entdeckten, dass der Einfluss der Energiemodelle nicht über alle kosmischen Ereignisse hinweg einheitlich ist. In Umgebungen, in denen der Neutronenvorrat knapp ist, spielt die Wahl des Modells kaum eine Rolle, da die Atome einfach ihren Brennstoff aufbrauchen, bevor sie signifikant verändert werden können. In Umgebungen mit einem reichen Neutronenvorrat und spezifischen thermodynamischen Bedingungen kann der Unterschied zwischen den Modellen jedoch darüber entscheiden, ob ein bestimmtes Element in großer Häufigkeit produziert wird oder kaum vorhanden ist. Die Forscher identifizierten, dass die signifikantesten Unsicherheiten in „kalten“ Szenarien auftreten, in denen die Atome auch lange nach dem Zusammenbruch des ursprünglichen Gleichgewichts weiterhin Neutronen einfangen. Unter diesen Bedingungen können kleine Änderungen in der vorhergesagten Geschwindigkeit des Neutroneneinfangs den gesamten Materiestrom verschieben und entscheiden, ob der Prozess beim zweiten Peak der schweren Elemente stoppt oder bis zum dritten durchstößt.

Letztendlich dient diese Arbeit als entscheidende Orientierungshilfe für zukünftige Experimente. Sie zeigt der Kernphysik genau auf, welche Regionen des Periodensystems am empfindlichsten auf unsere aktuellen theoretischen Lücken reagieren. Die Studie bestätigt, dass wir zwar ein gutes allgemeines Bild davon haben, wie schwere Elemente entstehen, die feinen Details der chemischen Landschaft jedoch stark davon abhängen, wie gut wir die internen Energiezustände instabiler Kerne verstehen. Die Forscher kommen zu dem Schluss, dass wir zur Verfeinerung unseres Verständnisses der chemischen Entwicklung des Universums bessere experimentelle Daten und präzisere theoretische Beschreibungen der Kernniveaudichten benötigen, insbesondere für die spezifischen Isotope, die zwischen den Hauptspitzen der schweren Elemente liegen. Ohtl dieser Verbesserungen werden unsere Simulationen des Kosmos weiterhin einen Schatten der Ungewissheit tragen und nicht in der Lage sein, die präzise Häufigkeit der schweren Elemente, aus denen unsere Welt besteht, vollständig zu erklären.