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

Diep in het hart van sterren is voortdurend een kosmische alchemie aan de gang, die de zware elementen smeedt waaruit onze wereld is opgebouwd. Terwijl lichtere elementen zoals koolstof en zuurstof worden geboren in de gestage verbranding van stellaire ovens, worden de zwaarste elementen, zoals goud en platina, gecreëerd in gewelddadige, extreme gebeurtenissen waarbij neutronen met een razendsnelle snelheid worden gevangen. Dit proces, bekend als het snelle neutronenvangstproces, of r-proces, is de primaire fabriek van het universum voor deze zware materialen. Om precies te begrijpen hoeveel van elk element wordt geproduceerd, moeten wetenschappers de exacte snelheden weten waarmee atoomkernen neutronen vangen. Veel van de onstabiele, neutronrijke kernen die betrokken zijn bij deze kosmische dans bestaan echter slechts gedurende fracties van een seconde, waardoor het onmogelijk is om ze direct in een laboratorium te bestuderen. Daarom moeten onderzoekers vertrouwen op theoretische berekeningen om het gedrag van deze vluchtige atomen te voorspellen, om zo de gaten op te vullen waar directe meting onmogelijk is.

Een cruciaal onderdeel van deze puzzel is hoe deze kernen interageren met licht, specif kind hoe ze energie absorberen via een type vibratie dat elektrische dipoolsterkte wordt genoemd. In eenvoudige termen beschrijft dit hoe de protonen en neutronen binnen een atoom ten opzichte van elkaar verschuiven wanneer ze door energie worden geraakt. Decennialang hebben wetenschappers brede, gegeneraliseerde modellen gebruikt om deze gedragingen voor onstabiele atomen te raden, maar deze gissingen missen vaak de precisie die nodig is om de creatie van elementen in het universum nauwkeurig te simuleren. Een nieuwe studie door Francesca Bonaiti en Joanna E. Sobczyk pakt deze uitdaging rechtstreeks aan door deze interacties vanaf de grond af aan te berekenen, met behulp van een methode die begint bij de fundamentele krachten tussen deeltjes in plaats van te vertrouwen op geïnformeerde gissingen.

De onderzoekers concentreerden zich op een specifieke familie van atomen: de nikkelisotopen. Nikkel is een sleutelspeler in het r-proces, omdat het fungeert als een tussenstap die de overvloed van elementen met massa's rond de honderd beïnvloedt. Het team gebruikte een krachtige computationele benadering genaamd "ab initio", wat "vanaf het begin" betekent. In plaats van hun berekeningen aan te passen aan bestaande gegevens, begonnen ze met de fundamentele wetten van de kernfysica, specifiek de interacties tussen protonen en neutronen afgeleid van een theorie genaamd chirale effectieve veldentheorie. Ze pasten dit toe op een keten van nikkelatomen, variërend van het stabiele nikkel-58 tot het extreem neutronrijke nikkel-80. Om de complexe wiskunde van deze onstabiele, open-schil kernen aan te pakken, ontwikkelden ze een nieuwe uitbreiding van een methode die bekend staat als de Lorentz integral transform coupled-cluster techniek. Dit stelde hen in staat om te simuleren hoe deze atomen reageren op energie zonder dat ze de onmogelijke verstrooiingstoestanden hoefden te construeren die dergelijke berekeningen gewoonlijk compliceren.

De resultaten onthulden een significant kenmerk in het gedrag van deze neutronrijke atomen: een duidelijke, laagliggende sterkte net boven het punt waar een neutron uit de kern kan ontsnappen. Dit fenomeen, vaak geassocieerd met een "pygmee-dipoolresonantie", werkt als een zachte, laagenergetische vibratie die prominenter wordt naarmate de kern meer neutronen krijgt. In het geval van nikkel-80 is deze laagenergetische vibratie zo sterk en ligt deze zo dicht bij de ontsnappingsdrempel dat het de interactie van het atoom begint te beïnvloeden bij verrassend lage temperaturen, rond de 0,5 miljard graden. Voor de andere isotopen in de keten wordt dit effect sterker naarmate het aantal neutronen toeneemt, waarbij het een belangrijke factor wordt bij temperaturen tussen de 3 en 10 miljard graden.

Toen de onderzoekers deze nieuwe, first-principles berekeningen in de standaard computercodes voerden die worden gebruikt om stellaire explosies te simuleren, was de impact duidelijk. De inclusie van deze specifieke elektrische dipoolsterkte verhoogde de voorspelde snelheden waarmee deze nikkelkernen neutronen vangen aanzienlijk. Voor de zwaardere isotopen steeg de vangstsnelheid met factoren variërend van 16 procent tot drie keer de vorige schattingen, afhankelijk van de temperatuur en het specifieke isotoop. Deze verbetering is geen kleine aanpassing; het vertegenwoordigt een fundamentele verschuiving in ons begrip van hoe deze kernen met hun omgeving interageren. De studie benadrukte ook dat de grootste bron van onzekerheid in hun resultaten voortkomt uit het wiskundige proces van het omkeren van hun berekeningen om tot het uiteindelijke antwoord te komen, en niet van de onderliggende fysieke modellen zelf.

Door de eerste rigoureuze, first-principles beschrijving van deze elektrische dipoolsterkte langs een isotopenketen te bieden, biedt het werk een nieuw, betrouwbaarder fundament voor astrofysische simulaties. Het suggereert dat de productie van zware elementen in het universum efficiënter kan zijn dan voorheen gedacht, met name voor de zwaarste nikkelisotopen. Hoewel de onderzoekers erkennen dat toekomstig werk de wiskundige inversie moet verfijnen om de onzekerheden te verminderen, slaagt deze studie erin de kloof te overbruggen tussen fundamentele kernfysica en de grootschalige vorming van kosmische elementen. Het demonstreert dat we zelfs voor de meest onstabiele atomen nu hun gedrag kunnen voorspellen met een niveau van detail dat voorheen onbereikbaar was, waardoor we dichter bij het oplossen van het mysterie komen van waar de zware elementen in ons universum werkelijk vandaan komen.