Numerical calculations of neutron star mountains supported by crustal lattice pressure

Deze studie berekent zelfconsistent temperatuur-asymmetrieën en massa-'bergen' op accreterende neutronensterren met realistische toestandsvergelijkingen, en concludeert dat deze bergen te klein zijn om de spin-equilibrium van lage-massa Röntgen-dubbelsterren te bepalen.

De kern

De onzichtbare berg op een sterrenwolk: Waarom neutronensterren geen zware 'berg' kunnen vormen

Stel je voor dat je een balletje van honderd kilometer dik en zo zwaar als de zon hebt. Dit is een neutronenster. Deze sterren draaien razendsnel om hun as, soms wel honderden keren per seconde.

Wetenschappers hopen dat deze sterren zwaartekrachtsgolven uitzenden. Dit zijn rimpelingen in de ruimtetijd, zoals golven in een vijver. Om deze golven te maken, moet de ster niet perfect rond zijn; hij moet een klein 'bultje' of een 'berg' hebben. Hoe groter die berg, hoe sterker de rimpeling.

De vraag is: Kan zo'n berg echt bestaan, en is hij groot genoeg om door onze telescopen (zoals LIGO) te worden gezien?

1. Het probleem: De ster moet een berg hebben

In het verleden dachten wetenschappers dat neutronensterren enorme bergen konden hebben, gemaakt door de hitte en druk in de buitenste laag (de korst) van de ster.

• De oude theorie: Ze dachten dat als het aan de ene kant van de ster heter is dan aan de andere kant, de atomen in de korst verschuiven. Denk aan een deken die oneffen wordt als je hem aan één kant verwarmt. Dit zou een enorme berg kunnen vormen.
• Het nieuwe idee: De auteurs van dit artikel (Hutchins en Jones) kijken naar een ander mechanisme. Ze zeggen: "Laten we kijken naar de kristalstructuur zelf." De korst van een neutronenster is niet zomaar stof; het is een gigantisch, superhard kristalrooster van atoomkernen. Als dit kristal ongelijkmatig wordt verwarmd, kan de druk in het kristalrooster veranderen en een berg vormen.

2. De berekening: Een gigantische puzzel

De auteurs hebben een zeer complexe computerberekening gemaakt. Ze hebben drie dingen tegelijkertijd berekend:

1. Hoe warm is de ster? (De temperatuur).
2. Hoe verandert de druk in het kristalrooster door die hitte? (De elastische spanning).
3. Hoe vormt dit een berg? (De massa-asymmetrie).

Ze gebruikten de meest moderne en nauwkeurige formules voor hoe materie zich gedraagt onder extreme druk (de zogenoemde BSk-models). Ze hebben ook gekeken naar het effect van magnetische velden. Magnetische velden in neutronensterren zorgen ervoor dat warmte niet gelijkmatig wordt verspreid, net zoals een deken die aan één kant dikker is dan aan de andere. Dit creëert de temperatuurverschillen die nodig zijn voor de berg.

3. Het verrassende resultaat: De berg is te klein

Het nieuws is helaas niet zo spannend als we hoopten.

• De verwachting: We hoopten dat deze bergen groot genoeg zouden zijn om de rotatie van de ster te vertragen en zwaartekrachtsgolven te maken die we kunnen zien.
• De realiteit: De berekeningen tonen aan dat de bergen die door dit kristal-mechanisme worden gevormd, ontzettend klein zijn.
  • Analogie: Stel je een berg voor op een planeet ter grootte van een tennisbal. De berg is dan niet groter dan een stofdeeltje.
  • De auteurs berekenden dat de 'berg' (in wetenschappelijke termen: de ellipticiteit) ongeveer 100.000 tot 1.000.000 keer kleiner is dan wat nodig zou zijn om de ster te vertragen of om door onze huidige telescopen te worden gezien.

4. Waarom is dit zo?

De reden is dat de druk in het kristalrooster, hoewel hij verandert door de hitte, verwaarloosbaar klein is vergeleken met de enorme druk van de vrije neutronen die de ster bij elkaar houden.

• Analogie: Het is alsof je probeert een berg te bouwen op een berg steen met een paar druppels water. De druppels (de warmte in het kristal) zijn niet sterk genoeg om de steen (de zwaartekracht en druk van de ster) te verplaatsen.

5. Wat betekent dit voor de toekomst?

• Geen teleurstelling, maar een realiteit: Dit betekent dat we waarschijnlijk geen zwaartekrachtsgolven zullen zien van deze specifieke 'thermische bergen' op neutronensterren in lage-massa X-ray binairs (LMXBs).
• Andere opties: Misschien moeten we zoeken naar andere soorten sterren, zoals Ultra-Luminous X-ray sources (ULXs). Dit zijn sterren die nog veel sneller en heter zijn, met veel sterkere magnetische velden. Daar zouden de bergen misschien wel groot genoeg kunnen zijn om te zien.
• Betere telescopen nodig: Om deze kleine bergen ooit te zien, hebben we toekomstige telescopen nodig die veel gevoeliger zijn dan de huidige generatie (zoals de Einstein Telescope of Cosmic Explorer).

Samenvatting in één zin

De auteurs hebben bewezen dat de 'bergen' die door hitte en magnetische velden in de korst van neutronensterren worden gevormd, veel te klein zijn om door onze huidige apparatuur te worden opgepikt als zwaartekrachtsgolven; de sterren zijn te perfect rond.

Het is alsof we hoopten een tsunami te zien veroorzaakt door een steen die in een meer valt, maar het bleek slechts een klein plasje water te zijn. We moeten nu zoeken naar de echte tsunami's (zoals bij ULX's) of wachten op betere apparatuur.

Technische samenvatting

Titel: Numerieke berekeningen van neutronenster-bergen ondersteund door kristalroosterdruk van de korst

Auteurs: T. J. Hutchins & D. I. Jones
Publicatie: MNRAS (2026)

---

1. Het Probleem

Gravitatiegolven (GW) worden geacht de rotatiefrequentie van neutronensterren in lage-massa röntgendubbelsters (LMXB's) te bepalen via een evenwichtstoestand (torque balance). Een mechanisme om dergelijke emissie mogelijk te maken, is de vorming van een massa-asymmetrie, of een "berg", die wordt ondersteund door elastische spanningen in de vaste korst van de ster.

Tot nu toe hebben studies zich vaak geconcentreerd op ofwel de oorsprong van de elastische spanningen, ofwel de temperatuur-asymmetrie, maar zelden op het volledige vormingsproces in één zelfconsistent model.

• Bestaande modellen: De traditionele benadering (het "UCB-model" van Ushomirsky et al., 2000) gaat uit van verschuivingen in elektron-vangstlagen ("capture layers") veroorzaakt door temperatuurgradiënten. Recentere berekeningen van de toestandsvergelijking (EOS) tonen echter aan dat er veel minder elektron-vangstlagen zijn dan eerder gedacht, en dat de drempelenergieën lager zijn. Dit ondermijnt de grootte van de bergen die via dit mechanisme kunnen worden gegenereerd.
• Het doel: De auteurs willen een nieuw, zelfconsistent model ontwikkelen dat de vorming van thermische bergen beschrijft via een alternatief mechanisme: de temperatuurafhankelijke druk van het ionische kristalrooster zelf ("crustal lattice pressure"), gekoppeld aan anisotrope warmtetransport door interne magnetische velden.

2. Methodologie

De auteurs voeren een volledig zelfconsistente berekening uit binnen één raamwerk, waarbij ze de volgende stappen doorlopen:

1. Toestandsvergelijking (EOS):

   • Voor het eerst worden realistische EOS-modellen gebruikt die alle regio's van een accreterende neutronenster thermodynamisch consistent beschrijven: de Brussel-Montreal functionals BSk19, BSk20 en BSk21.
   • Om numerieke convergentieproblemen door dichtheidsdiscontinuïteiten (bij kernreacties) te omzeilen, hebben de auteurs analytische benaderingen (fits) ontwikkeld voor deze EOS-data.

2. Hydrostatische en Thermische Structuur:

   • De achtergrondstructuur van de ster wordt berekend met de Tolman-Oppenheimer-Volkoff (TOV) vergelijkingen (algemene relativiteit), maar voor de thermische en elastische berekeningen wordt een Newtoniaans raamwerk gebruikt met een effectieve zwaartekrachtsversnelling om relativistische effecten benaderd te houden.
   • De thermische structuur wordt berekend voor sterren die steady accreteren, met inachtneming van diepe korstverwarming (door niet-evenwichtsreacties) en ondiepe korstverwarming (een onbekend mechanisme nodig om waarnemingen te verklaren).

3. Temperatuur-Asymmetrie:

   • De asymmetrie wordt gegenereerd door anisotrope warmtetransport veroorzaakt door interne magnetische velden (poloïdaal of toroïdaal). Dit leidt tot temperatuurvariaties ($\delta T$) in de korst.
   • Ze berekenen de perturbaties voor zowel korstvelden ($B \sim 10^{12}$ G) als kernvelden ($B \sim 10^8$ G).

4. Elastische Respons en Bergvorming:

   • In plaats van elektron-vangstlagen, gebruiken ze de thermische roosterdruk ($P_{th}$). Deze druk is temperatuurafhankelijk en wordt veroorzaakt door de interactie van ionen in het kristalrooster.
   • Ze lossen de lineaire perturbatievergelijkingen op voor de Lagrangiaanse verplaatsing ($\xi_i$) van de korstmateriaal als reactie op de temperatuurgradiënten en de resulterende drukonevenwichten.
   • Hieruit wordt de massakwadrupoolmoment ($Q_{22}$) en de ellipticiteit ($\varepsilon$) afgeleid.

3. Belangrijkste Bijdragen

• Zelfconsistentie: Dit is het eerste werk dat de temperatuur-asymmetrie, de perturbaties in de korstroosterdruk en de resulterende massa-asymmetrie berekent binnen één enkel, coherent model.
• Realistische EOS: Gebruik van de geavanceerde BSk19/20/21 EOS-modellen, die een volledig thermodynamisch consistent beeld geven van de accreterende korst, in tegenstelling tot oudere, gefragmenteerde modellen.
• Alternatief Mechanisme: De auteurs testen expliciet het mechanisme van "thermische roosterdruk" en vergelijken dit met het verouderde "elektron-vangstlaag"-mechanisme.
• Analytische Fits: Ze ontwikkelen nauwkeurige analytische fits voor de EOS om numerieke integratieproblemen rondom dichtheidsstappen op te lossen.

4. Resultaten

• Grootte van de Bergen: De berekende ellipticiteit ($\varepsilon$) voor deze "magneto-thermo-elastische" bergen is extreem klein. De maximale waarde die ze vinden is $\varepsilon_{max} \approx 2.32 \times 10^{-11}$.
• Vergelijking met Torque-Balance: Om de rotatie van LMXB's te bepalen (torque balance), zou een ellipticiteit in de orde van $\varepsilon \sim 10^{-9}$ tot $10^{-7}$ nodig zijn. De door de auteurs berekende bergen zijn vele ordes van grootte te klein om de spin-equilibrium van deze sterren te dictëren.
• Invloed van Magnetische Velden en Accretie:
  • De ellipticiteit neemt toe met de accretiegraad en de hoeveelheid ondiepe verwarming, maar bereikt een verzadigingspunt wanneer de accretie zo hoog is dat de korst smelt bij lagere dichtheden.
  • Zelfs met sterke korstmagnetische velden ($10^{12}$ G) blijft het effect ver onder de detectielimieten van huidige GW-detectors.
• Vergelijking met Elektron-Vangstlagen:
  • Als men alleen kijkt naar de elektron-vangstlagen die door moderne EOS-modellen worden voorspeld (zonder de "kunstmatige" hoge-drempel lagen uit oude modellen), is het mechanisme van elektron-vangstlagen ook veel minder effectief dan eerder gedacht.
  • De roosterdruk-mechanisme is vergelijkbaar in grootte met de verschuiving van ondiepe vangstlagen, maar beide zijn ver onder de theoretische maximale berggrootte ($\varepsilon_{max} \sim 10^{-7}$).

5. Betekenis en Conclusie

• Detectie van Continue Gravitatiegolven: De resultaten suggereren dat "thermische bergen" ondersteund door roosterdruk waarschijnlijk niet de bron zijn van de spin-equilibrium in LMXB's, en dat ze te klein zijn om met huidige generatie GW-detectors (zoals LIGO/Virgo/KAGRA) gedetecteerd te worden.
• Toekomstperspectief: Hoewel LMXB's waarschijnlijk geen sterke signalen leveren via dit mechanisme, wijzen de auteurs op Ultra-Luminous X-ray (ULX) bronnen. Deze systemen hebben extreme accretiegraden en sterke magnetische velden, wat de temperatuur-asymmetrie en de roosterdruk aanzienlijk kan verhogen. Hoewel hun rotatiesnelheid lager is, zouden ze potentieel interessante doelen kunnen zijn voor toekomstige, meer gevoelige detectors (zoals Einstein Telescope of Cosmic Explorer) die beter zijn op lage frequenties.
• Modelverbeteringen: De auteurs benadrukken dat toekomstig werk moet focussen op het volledig relativistisch maken van het model en het verwijderen van de Cowling-benadering, wat de berekende ellipticiteit mogelijk nog verder kan verlagen (door zwaartekrachtsperturbaties).

Samenvattend: Dit paper toont aan dat hoewel het mechanisme van thermische roosterdruk fysiek plausibel is, de resulterende neutronenster-bergen te klein zijn om de waargenomen rotatie-eigenschappen van accreterende neutronensterren te verklaren of directe detectie van continue gravitatiegolven te garanderen met huidige technologie.