Neutron Star Mass across Binary Pulsar Subpopulations: Mass-Spin Correlation, Mass Distributions, and Moment of Inertia Effects

Technische Samenvatting: Massa van Neutronensterren over Subpopulaties van Binaire Pulsars

Probleemstelling
De vorming en evolutie van neutronensterren (NS) in binaire systemen laten hun sporen na in hun massa's, rotatieperioden en baanparameters. Terwijl het "recycling"-paradigma stelt dat accretie van een begeleider neutronensterren opstuwt naar milliseconde-perioden, blijft de fysieke oorsprong van de waargenomen massaverdeling onderwerp van debat. Specifiek is het onduidelijk of zware gerecyclede pulsars ($\sim 1,8 M_\odot$ ) massief zijn geboren of zijn gegroeid via accretie. Bovendien voorspellen twee concurrerende fysieke mechanismen tegengestelde correlaties tussen de NS-massa ($M$) en de rotatieperiode ($P$) in gerecyclede populaties:

1. Accretie-gestuurde recycling: De toename van impulsmoment is gekoppeld aan massatoename, wat een negatieve correlatie voorspelt ($\rho < 0$; snellere rotaties geassocieerd met hogere massa's).
2. Traagheidsmoment (MoI)-mechanisme: Lagere-massa NS's hebben kleinere traagheidsmomenten en kunnen daardoor efficiënter worden opgestuwd voor een gegeven geaccreteerd impulsmoment, wat een positieve correlatie voorspelt ($\rho > 0$).
   Huidige steekproeven worden beperkt door meetonzekerheden en selectie-effecten, wat het moeilijk maakt om deze mechanismen te onderscheiden of te bepalen of de gerecyclede massaverdeling intrinsiek bimodaal is of een scheef unimodale populatie.

Methodologie
De auteurs presenteren een hiërarchische Bayesiaanse analyse van $\sim 50$ Galactische binaire radiopulsars met gemeten NS-massa's, geclassificeerd in Pulsar–Witte Dwerg (PSR–WD) en Dubbele Neutronenster (DNS) systemen.

• Data Curatie: De steekproef is samengesteld uit gevestigde pulsar-massacatalogi. Herhaalde metingen worden gecombineerd via inverse-variantie weging, met voorzichtigheid toegepast op systemen met discrepante metingen (bijv. J0348+0432, B1855+09).
• Classificatie: Systemen worden geclassificeerd op basis van binair type (PSR–WD primaires, DNS primaires, DNS secondaires) in plaats van rotatieperiode om circulariteit te vermijden. Witte dwerg-begeleiders worden onderverdeeld in Helium (He) of Koolstof-Zuurstof/Zuurstof-Neon (CO/ONe) op basis van een massadrempel van $0,45 M_\odot$.
• Statistisch Kader: Een hiërarchisch bivariaat Gaussisch model wordt ingezet om intrinsieke correlaties tussen massa en recycling-proxies (rotatieperiode, baanperiode, excentriciteit, inclinatie) te schatten. Deze aanpak modelleert meetfouten expliciet als latente variabelen om regressie-dilutie (attenuatiebias) te corrigeren.
• Selectiecorrecties: De analyse bevat correcties voor orbitale-Doppler detecteerbaarheid (met behulp van de $\gamma_1$ efficiëntiefactor) en houdt rekening met inclinatiebiases die inherent zijn aan Shapiro-delay massametingen.
• Modelvergelijking: Om het MoI-mechanisme te testen, vergelijken de auteurs regressiemodellen van de rotatieperiode tegen de massa direct versus tegen het traagheidsmoment (afgeleid van universele relaties), met behulp van het Bayesian Information Criterion (BIC).

Belangrijkste Bijdragen en Resultaten

1. Massa–Rotatieperiode Anti-correlatie:
   Het onderzoek vindt matig bewijs voor een anti-correlatie tussen de NS-massa en de rotatieperiode in de gepoolde gerecyclede populatie ($N=47$ primaires).

• Resultaat: $\rho = -0,26$ met een 90% geloofwaardigheidsinterval van $[-0,47, -0,01]$ en een eenzijdige posterior waarschijnlijkheid $P(\rho < 0) = 0,96$ (ongeveer $1,7\sigma$).
• Robuustheid: Het resultaat blijft standhouden wanneer DNS-kandidaten worden opgenomen en selectiecorrecties worden toegepast.
• Subpopulatie Nuance: De PSR–WD steekproef alleen vertoont een suggestieve trend ($\rho = -0,19$), terwijl de DNS-primaries alleen licht positief neigen ($\rho = +0,13$), hoewel de DNS-steekproef te klein is ($N=10$) om conclusief te zijn. De gepoolde significantie wordt grotendeels gedreven door de inclusie van DNS-primaries, die een gebied van lagere massa en langere periode bezetten, waardoor de negatieve diagonaal in het massa–spin vlak wordt uitgebreid.

2. Degeneratie van Mechanismen:
   Hoewel de negatieve correlatie consistent is met accretie-gestuurde recycling, kunnen de auteurs dit niet statistisch onderscheiden van het MoI-mechanisme.

• Reden: Over het geobserveerde massabereik is het traagheidsmoment van de NS $I(M)$ bijna lineair. Bijgevolg is een regressie van de rotatieperiode op de massa statistisch niet te onderscheiden van een regressie op $I(M)$ (BIC verschil $\Delta \text{BIC} = -0,2$).

3. Massaverdelingen en Bimodaliteit:

• Gerecycled vs. Niet-gerecycled: Gerecyclede primaires zijn significant massiever ($1,50 M_\odot$) dan niet-gerecyclede DNS secondaries ($1,30 M_\odot$), wat consistent is met eerdere bevindingen.
• Bimodaliteit: Binnen de gerecyclede populatie wordt een twee-Gaussische mengeling verkozen boven een enkele symmetrische Gaussische verdeling, maar een enkele rechts-scheve (skew-normal) verdeling biedt de meest parsimonieuze beschrijving (laagste BIC). Dit suggereert dat de schijnbare bimodaliteit een artefact kan zijn van een enkele populatie met een hoog-massa staart in plaats van twee verschillende geboortepopulaties.

4. Begeleider Type en Massa:
   Neutronensterren met Helium witte dwerg-begeleiders zijn marginaal massiever ($1,547 M_\odot$) dan die met CO/ONe-begeleiders ($1,484 M_\odot$), een verschil van $\sim 0,06 M_\odot$. Dit verschil is van marginale statistische significantie ($p=0,11$), maar komt overeen met de verwachting dat He-WD kanalen meer uitgebreide massatransfer omvatten.

5. DNS Begeleider Massa–Excentriciteit Correlatie:
   In een hiërarchisch kader bevestigen de auteurs een sterke positieve correlatie tussen de begeleider (tweede-geboren) massa en orbitale excentriciteit in DNS-systemen ($\rho = +0,82$). Dit ondersteunt de theoretische link tussen lage-kick supernova's (die systemen met lage massa en lage excentriciteit produceren) en de vorming van DNS's.

Significantie en Claims
Het artikel presenteert de bevindingen als matig of suggestief bewijs in plaats van een definitieve detectie, primair vanwege de bescheiden steekproefomvang ($N \approx 50$).

• Interpretatie: De resultaten ondersteunen een "twee-kanaals" beeld: PSR–WD systemen zijn "door accretie gegroeid", waarbij ze de massa–spin anti-correlatie en massa-offsets vertonen, terwijl DNS-systemen "geboortemassa-gedomineerd" zijn, waarbij massa's de vormingscondities reflecteren in plaats van de recyclinggeschiedenis.
• Beperkingen: Het onvermogen om tussen de accretie- en MoI-mechanismen te onderscheiden is een direct gevolg van de bijna lineaire relatie tussen massa en traagheidsmoment in het geobserveerde bereik.
• Toekomstperspectief: De auteurs voorspellen dat het bereiken van een betrouwbare $3\sigma$ detectie van de massa–spin correlatie een steekproefomvang van $N \gtrsim 125–150$ precieze-massa gerecyclede systemen vereist, waarschijnlijk een prestatie in het SKA-tijdperk. Ze merken op dat het doorbreken van de accretie/MoI degeneratie mogelijk vereist dat de spin en massa van niet-gerecyclede NS's (bijv. tweede-geboren pulsars in DNS-systemen) worden gemeten om geboorteneigenschappen te isoleren van recyclingeffecten.

Het werk herkadert het debat over gerecyclede massaverdelingen, waarbij wordt gesuggereerd dat de geobserveerde hoog-massa staart beter beschreven kan worden door een rechts-scheve unimodale verdeling dan door een afzonderlijke bimodale geboortepopulatie.