On a re-examination of neutron star cooling in transient sources -- No shallow heating required?

Technische Samenvatting: Een herbeoordeling van de afkoeling van neutronensterren in transiënte bronnen

Probleemstelling
De thermische evolutie van neutronensterren (NSs) in laagmassa-röntgenbinairsystemen (LMXBs) na een accretie-uitbarsting vereiste traditioneel de aanroeping van een extra, ad hoc warmtebron bekend als "ondiepe verhitting" (shallow heating). Dit mechanisme, waarvan de fysieke oorsprong onbekend is, wordt doorgaans gemodelleerd als zijnde proportioneel aan de massetoename van de accretiesnelheid ($\dot{M}$), en is noodzakelijk om de hoge effectieve temperaturen ($T_{\infty}^{\text{eff}}$) te reproduceren die worden waargenomen in de vroege afkoelingsfasen (honderden dagen na de uitbarsting). Standaardmodellering gaat ervan uit dat de envelop van de NS fungeert als een randvoorwaarde die onafhankelijk is van $\dot{M}$ en voortdurende thermonucleaire verbranding, waarbij de mogelijkheid wordt genegeerd dat warmte die in de envelop wordt gegenereerd tijdens accretie in de korst kan lekken.

Methodologie
De auteurs maken gebruik van de volledig relativistische 1D-code NSCool om de thermische evolutie van de korst en kern van de NS te simuleren. De studie introduceert een nieuwe oppervlakte-randvoorwaarde (BC) bij de overgang tussen de envelop en de korst (drempeldichtheid $\rho_b \approx 10^7$ g cm$^{-3}$). In tegenstelling tot eerdere modellen die een statische $T_b - T_{\text{eff}}$-relatie gebruiken, is deze "Nieuwe BC" gevoelig voor de massetoename van de accretiesnelheid ($\dot{M}$) en houdt deze rekening met de lekkage van thermonucleaire verhitting vanuit de envelop naar de korst, consistent met het "watershed-model" (Hanawa & Fujimoto 1984; Fujimoto et al. 1987).

De studie analyseert acht uitbarstingsepisoden van zeven verschillende bronnen: MXB 1659-29 (twee uitbarstingen), XTE J1701-462, EXO 0748-676, IGR J17480-2446 (Ter5 X-2), Swift J174805.3-244637 (Ter5 X-3), 1RXS J180408.9-342058, en KS 1731-260. Deze bronnen zijn geselecteerd vanwege hun diverse duur van de uitbarsting ($t_{\text{acc}} \sim 0,1$ tot 24 jaar) en accretiesnelheden ($10^{16}$ tot $10^{18}$ g s$^{-1}$).

Twee passende scenario's werden getest met behulp van Monte Carlo Markov Chain (MCMC)-methoden:

• Scenario A: Modellen die afkoelen zonder ondiepe verhitting. De vrije parameters zijn $\dot{M}$, de kolomdiepte van lichte elementen ($y_L$), de roodverschoven kerntemperatuur ($\tilde{T}_{\text{core}}$), en twee onzuiverheidsfactoren ($Q_{\text{imp},1}$ voor de buitenste korst, $Q_{\text{imp},2}$ voor de binnenste korst).
• Scenario B: Een nultest inclusief ondiepe verhitting. Dit voegt drie parameters toe: de verhittingssterkte ($Q_{\text{sh}}$) en het dichtheidsbereik $[\rho_{\text{sh,min}}, \rho_{\text{sh,max}}]$ waar deze opereert.

Belangrijkste Resultaten

1. Succes van Scenario A: De nieuwe randvoorwaarde reproduceert succesvol de afkoelingscurves voor zeven van de acht bronnen zonder dat ondiepe verhitting nodig is. Het model verklaart op natuurlijke wijze de hoge $T_{\infty}^{\text{eff}}$ direct na de uitbarstingen.
2. Consistentie van Parameters: In Scenario A komen de passende $\dot{M}$-waarden over het algemeen goed overeen met de observationeel afgeleide gemiddelde accretiesnelheden (met uitzondering van Ter5 X-2, waar de fit een snelheid vereist die 3–4 keer hoger is). De onzuiverheidsfactoren vertonen consistent een dalende trend met de diepte ($Q_{\text{imp},1} > Q_{\text{imp},2}$), met magnituden van $\sim 10$ in de buitenste korst en $\sim 1$ in de binnenste korst. De kolomdiepte van de lichte elementen wordt consistent gevonden rond $y_L \sim 10^8$ g cm$^{-2}$.
3. Falen van Scenario B (Fysieke Plausibiliteit): Hoewel Scenario B de data kan fitten, vereist dit onfysische parameters. Specifiek dwingt het de passende $\dot{M}$ een orde van grootte lager te zijn dan de observationele schattingen, gecompenseerd door een excessief hoge ondiepe verhittingssterkte ($Q_{\text{sh}} \gtrsim 5$ MeV baryon$^{-1}$).
4. Modelselectie: Informatiecriteria (AICc en BIC) geven de voorkeur aan Scenario A voor 5 van de 7 bronnen. Voor EXO 0748-676 is Scenario B statistisch geprefereerd, maar de vereiste ondiepe verhittingssterkte blijft onfysisch groot. De afkoelingscurve voor EXO 0748-676 is het enige geval waarbij Scenario A geen bevredigende fit biedt, wat mogelijk te wijten is aan de hoge massa ($2M_{\odot}$) en de lange duur van de uitbarsting van de bron.

Betekenis en Claims
Het artikel betoogt dat de "ondiepe verhitting" die voorheen werd aangeroepen in NS-afkoelingsmodellen, waarschijnlijk een artefact is van vereenvoudigde oppervlakte-randvoorwaarden die de koppeling tussen de accretiesnelheid, thermonucleaire verbranding in de envelop en de onderliggende korst verwaarlozen. Door een randvoorwaarde te implementeren die energie van de envelop naar de korst laat stromen (het watershed-model), demonstreren de auteurs dat het waargenomen afkoelingsgedrag van transiënte bronnen kan worden verklaard met standaardfysica zonder ad hoc verhittingsmechanismen.

De auteurs concluderen dat thermonucleaire verhitting in de NS-envelop een belangrijkere rol speelt bij het vormgeven van afkoelingscurves dan voorheen werd aangenomen. Hoewel zij erkennen dat verbeteringen mogelijk zijn—zoals het opnemen van tijd-afhankelijke accretiesnelheden, variërende NS-massa/straal, en het modelleren van onzuiverheden als continue functies van de diepte—suggereren hun resultaten sterk dat ondiepe verhitting niet vereist is voor veel, indien niet alle, transiënte bronnen. De uitzondering van EXO 0748-676 suggereert dat aanvullende fysica nog steeds nodig kan zijn voor specifieke hoog-massieve of langdurige systemen.