An Investigation of Systematic Effects from Background Priors on PSR J0740$+$6620 Radius Estimates using Synthetic NICER and XMM-Newton Data

Technische Samenvatting: Een onderzoek naar systematische effecten van achtergrondpriors op de straal-schattingen van PSR J0740+6620 met behulp van synthetische NICER- en XMM-Newton-data

Probleemstelling
Nauwkeurige metingen van neutronensterstralen zijn cruciaal voor het beperken van de toestandsvergelijking (EOS) van koud, dicht materie bij supranucleaire dichtheden. Hoewel analyses van synthetische X-ray puls-golfvormdata van de Neutron Star Interior Composition Explorer (NICER) hebben aangetoond dat ze robuust zijn tegen systematische fouten gerelateerd aan thermische emissiepatronen (bijv. de vorm van vlekken en temperatuurverdelingen), blijft er een potentieel significante maar nog niet verkende bron van systematische fout: de misparameterisering van ongemoduleerde achtergrondcomponenten.

Voor de zware pulsar PSR J0740+6620 ($M \approx 2,1 M_\odot$) wordt de geobserveerde telrate gedomineerd door bronnen die niet de oppervlakte-hotspots zijn, waaronder de diffuse röntgenachtergrond en instrumentele ruis. Het correct schatten van deze ongemoduleerde achtergrond is essentieel omdat een overschatting van de achtergrondflux impliceert dat de bijdrage van de hotspots aan de flux lager is. Aangezien de modulatieamplitude in NICER-data vaststaat, vereist een lagere spot-flux een grotere stellaire straal om de geobserveerde modulatiefractie te behouden. Daarom kunnen fouten in het aangenomen achtergrondmodel de straalschattingen beïnvloeden, met name bij de gezamenlijke analyse van NICER- en XMM-Newton-data.

Methodologie
Om deze effecten te onderzoeken, genereerden de auteurs twee sets synthetische data voor PSR J0740+6620, waarbij de observaties van zowel NICER als XMM-Newton werden gesimuleerd:

1. Dataset 1 (Standaard): Synthetische XMM-Newton-data die alleen de "blank-sky" achtergrond bevatten, consistent met langetermijn-gemiddelde observaties van de missie.
2. Dataset 2 (Verhoogd): Synthetische XMM-Newton-data waarbij de achtergrond kunstmatig met een factor $\sim 5$ is verhoogd ten opzichte van het blank-sky niveau, om scenario's te simuleren met aanzienlijke extra ongemoduleerde flux (bijv. door kosmische straling of nabijgelegen heldere bronnen).

Beide datasets maakten gebruik van de best passende puls-golfvormmodelparameters van Miller et al. (2021), uitgaande van twee cirkelvormige hotspots met een uniforme temperatuur op een niet-magnetische waterstofatmosfeer.

De auteurs analyseerden deze synthetische datasets vervolgens met behulp van twee verschillende puls-golfvormmodellen:

• Standaard Model: Gaat ervan uit dat de XMM-Newton-achtergrond strikt de blank-sky component is, behandeld als een door Poisson beperkte prior.
• Uitgebreid Model (Extended Model): Staat een additionele, ongemoduleerde achtergrondcomponent toe in de XMM-Newton-data, begrensd door een gladgestreken versie van de totale geobserveerde telrate.

De analyse maakte gebruik van een hybride sampling-aanpak, waarbij het nested sampling-algoritme MultiNest werd gebruikt voor de initiële verkenning van de parameterruimte, gevolgd door parallel-tempered Markov Chain Monte Carlo (ptemcee) om robuuste convergentie en nauwkeurige posterior-estimatie te waarborgen. De prestaties van het model werden geëvalueerd met behulp van $\chi^2$-statistieken voor de kwaliteit van de fit en de Bayesiaanse bewijskracht ($\log Z$) voor modelvergelijking.

Belangrijkste Resultaten
De studie onderzocht vier specifieke gevallen die de twee datasets combineerden met de twee analysismodellen:

1. Geval 1: Gegevens met een verhoogde achtergrond geanalyseerd met het Uitgebreide Model (Correct Model).
2. Geval 2: Gegevens met een verhoogde achtergrond geanalyseerd met het Standaard Model (Incorrect Model).
3. Geval 3: Standaard achtergrondgegevens geanalyseerd met het Uitgebreide Model (Over-geparameteriseerd Model).
4. Geval 4: Standaard achtergrondgegevens geanalyseerd met het Standaard Model (Correct Model).

De belangrijkste bevindingen zijn:

• Robuustheid van Straalschattingen: Zelfs wanneer het analysemodel een achtergrond aannam die aanzienlijk lager was dan de werkelijke achtergrond (Geval 2, waar het model de achtergrond met een factor $>5$ onderschatte), was de verschuiving in de geïnferieerde straal-posterior ten opzichte van de ware waarde slechts $\sim 1\sigma$. De ware straal bleef binnen de geloofwaardige regio van de posterior-distributie.
• Discriminatie door Bayesiaanse Bewijskracht: In alle gevallen bevoordeelde de Bayesiaanse bewijskracht het model dat de achtergrondcondities van de synthetische data correct beschreef. Specifiek was de bewijskracht voor het correcte achtergrondmodel groter dan voor het incorrecte model, wat aantoont dat Bayesiaanse bewijskracht effectief kan differentiëren tussen achtergrondmodellen, zelfs wanneer de straal-posterior niet sterk vertekend is.
• Convergentie en Onzekerheid: De studie benadrukte dat analyses met enkel MultiNest de parameteronzekerheden kunnen onderschatten (het produceren van onterecht nauwe betrouwbaarheidsintervallen). De hybride MultiNest+ptemcee-aanpak leverde bredere, meer betrouwbare geloofwaardige regio's op, waarbij de ware geïnjecteerde waarden in de meerderheid van de gevallen binnen de $1\sigma$ geloofwaardige regio's vielen.

Betekenis en Claims
Het artikel concludeert dat analyses van NICER-achtige data, wanneer statistische sampling grondig is en het model de data goed beschrijft, nauwkeurige metingen van neutronensterstralen opleveren, zelfs in de aanwezigheid van significante systematische fouten in de aangenomen achtergrondpriors. De resultaten suggereren dat de straal-inferentie robuust is tegen het specifieke type achtergrond-misparameterisering dat is getest. Bovendien valideert de studie het nut van Bayesiaanse bewijskracht bij het selecteren van geschikte achtergrondmodellen, wat het vertrouwen in de gezamenlijke analyse van NICER- en XMM-Newton-data voor het bepalen van de eigenschappen van zware neutronensterren zoals PSR J0740+6620 versterkt.