Controlling artificial surface heating in neutron star simulations: Application to hybrid equations of state

Neutronensterren zijn de dichtste objecten in het universum, de ingestorte kernen van massieve sterren die hun brandstof hebben verbruikt. Stel je een stadsgrote sfeer voor die meer massa bevat dan onze hele zon, zo compact samengeperst dat een enkele theelepel van het materiaal ervan miljarden tonnen zou wegen. Deze kosmische laboratoria stellen wetenschappers in staat om materie te bestudelen onder omstandigheden van druk en dichtheid die niet op aarde kunnen worden nagebootst. Om te begrijpen hoe deze sterren zich gedragen, hoe ze bewegen en hoe ze tegen elkaar aan kunnen botsen, vertrouwen onderzoekers op krachtige computersimulaties. Deze digitale modellen fungeren als virtuele laboratoria die complexe vergelijkingen oplossen om de evolutie van de sterren en de zwaartekrachtgolven die zij uitzenden bij botsingen te voorspellen. Echter, voor deze voorspellingen om nuttig te zijn, moeten de computermodellen uiterst nauwkeurig zijn, vooral aan de uiterste rand van de ster waar de dichte materie de lege vacuüm van de ruimte ontmoet.

De grens tussen een neutronenster en de omringende leegte is een plek van extreme moeilijkheid voor computersimulaties. In de echte wereld is deze overgang scherp, maar in een computer wordt de ster gerepresenteerd door een rooster van punten, zoals een digitaal raster. Terwijl de ster binnen dit raster beweegt of van vorm verandert, worstelt de computer met het verwerken van de plotselinge overgang van dicht materiaal naar nietsheid. Deze strijd creëert een subtiele maar hardnekkige fout: de computer voegt per ongeluk warmte toe aan het oppervlak van de ster waar geen fysieke warmte zou mogen bestaan. Deze kunstmatige verhitting wordt niet veroorzaakt door wrijving of kernreacties; het is een foutje in de wiskunde. Na verloop van tijd verspreidt deze valse warmte zich naar binnen, waardoor de ster licht uitzet en de interne energie verandert. Als dit niet wordt gecontroleerd, kan deze numerieke fout de simulatie vervormen, wat ertoe kan leiden dat wetenschappers onjuiste conclusies trekken over de structuur van de ster of de signalen die zij naar het universum stuurt.

In een recente studie onderzocht een onderzoeker van de Goethe Universiteit Frankfurt hoe deze kunstmatige verhitting te stoppen. Het werk richtte zich op een specifiek type computercode die wordt gebruikt om neutronensterren te simuleren, bekend als de BAM-code. Deze code wordt veelvuldig gebruikt om de gewelddadige botsingen van binaire neutronensterren te modelleren, waarbij twee dergelijke sterren naar elkaar toe spiralen en uiteindelijk samensmelten. De onderzoeker testte een nieuwe wiskundige strategie genaamd entropy-based flux-limiting (entropie-gebaseerde fluxbegrenzing). In eenvoudige termen werkt deze methode als een slim filter voor de berekeningen van de computer. Het monitort de stroom van energie en materie aan het oppervlak van de ster en past de berekeningsmethode automatisch aan wanneer het de soort steile veranderingen detecteert die gewoonlijk fouten veroorzaken. Door dit te doen, voorkomt het dat de computer die ongewenste warmte in de eerste plaats genereert.

Om te zien of deze methode werkte, voerde de onderzoeker twee sets simulaties uit. De eerste set betrof geïsoleerde neutronensterren, die alleen in de ruimte liggen, terwijl de tweede set paren sterren simuleerde die naar elkaar toe spiralen. In beide scenario's kreeg de computer een keuze: de standaard berekeningsmethode gebruiken of de nieuwe entropie-gebaseerde filter toepassen. De resultaten waren duidelijk en consistent. Wanneer de standaardmethode werd gebruikt, nam de interne energie van de sterren in de loop van de tijd gestaag toe en werden de sterren zelf groter, alsof ze opblazen door de valse warmte. Deze expansie was geen fysieke realiteit, maar een symptoom van de numerieke fout. Wanneer de nieuwe filter daarentegen werd toegepast, bleven de sterren stabiel. Hun interne energie bleef dicht bij de beginwaarde en hun grootte veranderde niet kunstmatig. De nieuwe methode slaagde erin de sterren te behoeden voor verhitting en expansie, waardoor een veel getrouwere representatie van hun fysieke staat werd behouden.

De studie onderzocht een breed scala aan neutronenstermodel, elk gebaseerd op verschillende theorieën over hoe materie zich gedraagt bij zulke extreme dichtheden. Sommige modellen gingen ervan uit dat de sterren uit één type deeltje bestonden, terwijl andere complexere mengsels bevatten. Ongeacht welk model werd gebruikt, de nieuwe filtermethode presteerde even goed. Het verminderde de kunstmatige verhitting met een aanzienlijke marge en hield de groei van de interne energie op een fractie van wat met de oude methode werd gezien. In de simulaties van binaire sterren, waar de zwaartekracht en de getijdenkrachten intens zijn, bleef de nieuwe methode de verhitting onderdrukken tijdens de lange dans van de inspiral, tot aan het moment van botsing. De sterren die met de nieuwe methode werden gesimuleerd, bleven kleiner en koeler dan die met de oude methode, wat suggereert dat de digitale artefacten effectief werden verwijderd.

Dit werk beweert niet dat het elk probleem in de simulatie van neutronensterren heeft opgelost, noch suggereert het dat de nieuwe methode een magische oplossing is voor alle numerieke uitdagingen. In plaats daarvan levert het sterk bewijs dat een specifieke bron van fouten — de kunstmatige verhitting aan het oppervlak — kan worden gecontroleerd. De onderzoeker merkt op dat hoewel deze tests werden uitgevoerd met een vereenvoudigde manier om de warmte van de ster te modelleren, het succes van de methode suggereert dat het net zo goed zal werken wanneer het wordt toegepast op complexere, meer realistische modellen die temperatuur in meer detail volgen. De bevindingen bieden een pad vooruit voor het creëren van nauwkeurigere simulaties van botsingen tussen neutronensterren. Door deze digitale verstoringen te verwijderen, kunnen wetenschappers hun modellen dieper vertrouwen, zodat zij ervoor kunnen zorgen dat de signalen die zij voorspellen voor detectoren van zwaartekrachtgolven gebaseerd zijn op de ware fysica van het universum in plaats van op de beperkingen van de computercode.