The Mass Function of Neutron Stars from Core-Collapse Supernova Simulations

Diep in het hart van massieve sterren speelt zich een laatste, gewelddadig drama af. Wanneer een ster die veel zwaarder is dan onze Zon zijn brandstof opgebruikt, stort de kern onder zijn eigen gewicht in, wat een supernova-explosie teweegbrengt die een hele melkweg kan overstralen. In de nasleep overleeft de kern vaak als een neutronenster, een object dat zo dicht is dat een enkele theelepel van het materiaal evenveel zou wegen als een berg. Decennialang hebben astronomen geprobeerd de geboorte van deze kosmische overblijfselen te begrijpen. Ze hebben de massa's van neutronensterren gemeten die in binaire systemen worden gevonden, waar twee sterren om elkaar heen draaien, maar de redenen achter de specifieke verdeling van deze gewichten bleven ongrijpbaar. Is er een standaardgewicht voor een pasgeboren neutronenster, of worden ze geboren met een grote verscheidenheid aan massa's? En hoe bepaalt het chaotische proces van de dood van een ster of hij een neutronenster wordt of volledig instort tot een zwart gat?

Een nieuwe studie door een team van onderzoekers brengt de meest geavanceerde computersimulaties van deze stellaire explosies samen met de nieuwste waarnemingen van neutronensterren om deze vragen te beantwoorden. Door de interne structuur van de stervende ster in kaart te brengen aan de hand van het gewicht van het overblijfsel dat het achterlaat, heeft het team een theoretisch beeld geconstrueerd van hoe neutronensterren worden geboren. Hun werk suggereert dat de variëteit aan massa's die we vandaag de dag zien niet alleen een resultaat is van sterren die later in hun leven aan gewicht winnen, maar in het allereerste moment van hun creatie geschreven staat. De bevindingen bieden een coherente verklaring voor waarom de meeste neutronensterren rond een specifiek gewicht clusteren, terwijl een kleinere groep zich naar veel zwaardere massa's uitstrekt, en ze geven een nieuwe schatting voor hoe vaak massieve sterren er helemaal niet in slagen om neutronensterren te worden, maar in plaats daarvan instorten tot zwarte gaten.

De onderzoekers begonnen door te kijken naar de ingewikkelde relatie tussen de initiële grootte van een massieve ster en de structuur van zijn kern vlak voordat hij explodeert. Ze gebruikten geavanceerde driedimensionale computermodellen die de complexe fysica van een supernova simuleren, waarbij ze bijhouden hoe de kern van de ster zich gedraagt in de laatste momenten van zijn leven. Deze simulaties onthulden een duidelijk patroon: de interne "compactheid" van de kern van de ster — een maatstaf voor hoe dicht de materie is samengeperst op een specifiek punt binnenin — bepaalt direct de massa van de neutronenster die wordt gevormd. Een minder compacte kern laat meestal een lichtere neutronenster achter, terwijl een compactere kern resulteert in een zwaardere. Deze connectie stelde het team in staat om de bekende verdeling van massieve sterren in het universum te nemen en de verdeling van neutronenster-massa's op het moment van hun geboorte te voorspellen.

Wanneer zij deze theoretische geboorteverdeling vergeleken met de werkelijke massa's van neutronensterren die in de hemel worden waargenomen, ontstond er een discrepantie. De simulaties voorspelden dat de meest voorkomende pasgeboren neutronensterren iets lichter zouden moeten zijn dan de gemiddelde waargenomen massa, met een secundaire groep zwaardere sterren. De werkelijke gegevens vertoonden echter een piek bij een iets hoger gewicht. De onderzoekers realiseerden zich dat dit verschil verklaard kon worden door twee hoofdfactoren: het feit dat neutronensterren massa winnen nadat ze zijn geboren, en de manier waarop hun geboorte bepaalt of ze in een binair systeem blijven.

Neutronensterren in binaire systemen trekken vaak materiaal aan van hun begeleidende sterren, een proces dat accretie wordt genoemd, wat het gewicht van de neutronenster in de loop van de tijd toevoegt. Het team paste correcties toe op de geobserveerde gegevens om te schatten wat de massa's bij de geboorte zouden zijn geweest, waarbij het gewicht dat door hun begeleiders is gewonnen, werd weggehaald. Deze aanpassing verschoof de geobserveerde piek omlaag, waardoor deze veel beter overeenkwam met de theoretische voorspelling. Bovendien benadrukte de studie de rol van "kicks" (stoten) die neutronensterren bij de geboorte krijgen. Net zoals een raket naar voren wordt voortgestuwd door zijn uitlaatgassen, ontvangt een neutronenster een plotselinge duw in de tegenovergestelde richting van de explosie. De simulaties suggereren dat zwaardere neutronensterren gemiddeld sterkere kicks ontvangen. Deze krachtige kicks kunnen de gravitationele band die het binaire systeem bij elkaar houdt verbreken, waardoor de zwaardere neutronensterren in de leegte worden uitgeworzigd waar ze moeilijker te vinden zijn. Deze selectieve verwijdering betekent dat de binaire systemen die we observeren bevooroordeeld zijn naar lichtere neutronensterren die mildere kicks ontvingen en in hun baan bleven.

Door de effecten van zwart gat-vorming, massaacretie en deze geboortekicks te combineren, creëerde het team een verfijnd model dat de geobserveerde gegevens opmerkelijk goed benadert. Hun model suggereert dat het universum een breed spectrum aan neutronenster-massa's produceert bij de geboorte, variërend van ongeveer 1,2 tot 2,0 keer de massa van onze Zon. De meest voorkomende geboortemassa is rond de 1,35 keer de massa van de Zon, wat de prominente piek verklaart die in observaties wordt gezien. De zwaardere staart van de verdeling, die zich uitstrekt tot 2,0 zonnemassa's, is een natuurlijk resultaat van de variërende structuren van de progenitor-sterren, in plaats van enkel het gevolg van latere accretie.

De studie leverde ook een nieuwe schatting voor het deel van de massieve sterren die hun leven als zwart gat eindigen in plaats van als neutronenster. Door te analyseren welke typen sterren in hun simulaties er niet in slaagden een neutronenster te produceren, berekenden de onderzoekers dat ongeveer 21 procent van de massieve sterren in het relevante massabereik direct instort tot zwarte gaten. Dit getal is geen eenvoudige regel gebaseerd op de initiële grootte van de ster; het hangt eerder af van de complexe, niet-lineaire structuur van de kern van de ster, wat betekent dat sterren van vergelijkbare grootte zeer verschillende loten kunnen hebben.

Dit werk vormt een belangrijke stap voorwaarts in het verbinden van de theorie over hoe sterren sterven met de realiteit van wat we waarnemen. Het gaat voorbij aan het oude idee dat alle neutronensterren met bijna dezelfde massa worden geboren, en vervangt dit door een dynamisch beeld waarbij de initiële structuur van de ster een breed scala aan uitkomsten dicteert. De overeenkomst tussen de complexe simulaties en de geobserveerde gegevens suggereert dat ons begrip van supernova-explosies volwassen wordt. Hoewel er nog steeds onzekerheden zijn, zoals de precieze invloed van de chemische samenstelling van de ster of de details van binaire interacties, demonstreert de studie dat de moderne supernova-theorie nu in staat is om de populatie van compacte objecten die we in de hemel zien te verklaren. Het opent een nieuw tijdperk waarin de geboorte van deze extreme objecten kan worden teruggevoerd naar de specifieke condities van de sterren die hen creëerden, waardoor een verzameling metingen verandert in een coherent verhaal van stellaire dood en wedergeboorte.