The NANOGrav 15 yr and 20 yr Datasets: Timing Events and Pulse Shape Changes

Titel: De NANOGrav 15-jaar en 20-jaar datasets: Timing-evenementen en veranderingen in pulsvorm

Datum: 8 april 2026 (conceptversie)
Auteurs: Ben Jacobson-Bell et al. (NANOGrav Collaboration)

---

1. Het Probleem

Pulsars, en met name millisecondepulsars (MSP's), vormen de ruggengraat van Pulsar Timing Arrays (PTA's) zoals NANOGrav, die worden gebruikt om het stochastic gravitatiegolfachtergrondsignaal (GWB) in het nanohertz-bereik te detecteren. Deze detectie vereist dat de aankomsttijden van pulsen (TOA's) met extreem hoge precisie (< 1 µs) kunnen worden gemodelleerd.

Hoewel de gemiddelde pulsvorm van een pulsar over decennia stabiel wordt verondersteld, treden er soms afwijkingen op. Deze omvatten:

• Timing-evenementen: Sudden afwijkingen in de TOA's die niet door het standaard timingmodel worden verklaard.
• Pulsvormveranderingen: Discrete, plotselinge veranderingen in de morfologie van de puls (vaak binnen < 30 dagen), gevolgd door een geleidelijke terugkeer naar de oorspronkelijke vorm.

Deze gebeurtenissen introduceren niet-stationaire ruis die de gevoeligheid van PTA's voor gravitatiegolven kan verstoren. Tot nu toe zijn er slechts een paar bekende gevallen van dergelijke gebeurtenissen in NANOGrav-pulsars (o.a. PSR J1713+0747 en PSR J1643−1224), maar een systematische zoektocht naar nieuwe gevallen ontbrak. Het is cruciaal om deze te identificeren om ze in timingmodellen te kunnen corrigeren of uit de dataset te verwijderen.

2. Methodologie

De auteurs ontwikkelen en toepassen een nieuwe detectiemethode die Hoofdkomponentenanalyse (PCA) combineert met matched filtering.

• Data Voorbewerking:

  • Er worden pulsprofielen van negen NANOGrav-pulsars geanalyseerd (uit de 15-jaar dataset en uitgebreid met data tot 2022/2023 voor specifieke pulsars).
  • Profielen worden uitgelijnd met template-matching, genormaliseerd en gefilterd op signaal-ruisverhouding (S/N).
  • De data worden onderverdeeld in 100 MHz subbanden om chromaticiteit (frequentie-afhankelijkheid) te onderzoeken.

• Principal Component Analysis (PCA):

  • PCA wordt toegepast op de residuen van de pulsprofielen (het verschil tussen het waargenomen profiel en een frequentie-afhankelijke template).
  • Dit decomposeert de variabiliteit in een basis van eigenvectoren (PC's). De projectie van de residuen op deze PC's ($D_i$) levert tijdsreeksen op die gevoelig zijn voor vormveranderingen.
  • Door de residuen te gebruiken in plaats van de ruwe profielen, wordt de bekende, stabiele frequentie-afhankelijkheid van de puls verwijderd.

• Zoekalgoritme (Matched Filtering):

  • De auteurs zoeken naar specifieke patronen in de $D_i$-tijdsreeksen: een snelle stijging gevolgd door een exponentiële afname (FRED-model: Fast Rise, Exponential Decay).
  • Ze gebruiken een dubbele matched filter:
    1. Een FRED-template (snelle stijging, trage daling).
    2. Een DERF-template (omgekeerd: trage stijging, snelle daling).
  • Uitsluitingscriterium: Echte FRED-evenementen zullen een sterke respons geven op de FRED-template en een zwakkere, verschoven respons op de DERF-template. Single-epoch outliers (instrumentale fouten) zullen vergelijkbare responsen geven op beide.
  • Ranking Metric ($\zeta$): Een score wordt berekend op basis van de tijdsvertraging ($\ell$) en de amplitudeverhouding ($\alpha$) tussen de FRED- en DERF-peak. Een lage $\zeta$-score (dicht bij 0) duidt op een sterke kandidaat voor een echt pulsvormveranderings-evenement.

3. Belangrijkste Bijdragen

1. Ontwikkeling van een robuust detectiemiddel: Een geautomatiseerde methode die PCA en matched filtering combineert om zowel bekende als nieuwe pulsvormveranderingen te vinden, inclusief discriminatie tegen instrumentale outliers.
2. Systematische analyse: Toepassing op een sample van negen NANOGrav-pulsars, waaronder de meest stabiele (null-testen) en die met bekende problemen.
3. Karakterisering van chromaticiteit: Analyse van de frequentie-afhankelijkheid van deze evenementen via subband-profielen.
4. Ontdekking van herhaling: Het identificeren van een opvallende herhaling van een langzame pulsvormverandering in PSR B1937+21 na ongeveer 10 jaar.

4. Resultaten

De analyse leverde de volgende resultaten op:

• Herkenning van bekende gebeurtenissen: De methode slaagde erin alle vier de eerder bekende pulsvormveranderings-evenementen te detecteren:

  • Drie gebeurtenissen in PSR J1713+0747 (2008, 2016, 2021).
  • Eén gebeurtenis in PSR J1643−1224 (2015).
  • De methode bevestigt dat deze gebeurtenissen intrinsieke pulsveranderingen zijn en niet alleen dispersiemaat (DM) variaties.

• Nieuwe kandidaat-evenementen: Vier nieuwe kandidaten met een $\zeta$-score beter dan (lager dan) de slechtste bekende gebeurtenis werden gevonden:

  • PSR B1937+21: Twee gebeurtenissen die corresponderen met langdurige "verstrooiingsvariabiliteit" (scattering variability). Opvallend is dat een gebeurtenis uit ~2011 bijna identiek herhaalde in ~2021, met een periode van ongeveer 10 jaar.
  • PSR J1600−3053: Een lage-S/N kandidaat die samenvalt met een periode die eerder uit de timing-analyse was verwijderd vanwege anomalieën in de dispersiemaat (DMX).
  • PSR J0030+0451: Een lage-S/N kandidaat in de hoofd-puls, met een gerelateerd signaal in de interpuls.

• Karakterisering van PSR B1937+21:

  • De herhalende gebeurtenis in PSR B1937+21 is het sterkst bij ~800 MHz.
  • Instrumentale oorzaken en zonnestelsel-effecten (zoals de zonnewind) werden uitgesloten.
  • De chromaticiteit suggereert interstellaire verstrooiing (scattering), hoewel de bron van deze verstrooiing onbekend blijft (er is geen bekend HII-gebied langs de lijn van zicht).

• Null-testen: Voor de minst "ruisende" pulsars (PSR J1909−3744 en PSR J2234+0611) werden geen significante FRED-achtige gebeurtenissen gevonden, wat de betrouwbaarheid van de methode bevestigt.

5. Betekenis en Conclusie

Deze studie is van groot belang voor de toekomst van de NANOGrav en andere PTA's:

• Verbetering van Timing-nauwkeurigheid: Door deze evenementen te identificeren, kunnen timinganalisten ze beter modelleren of specifieke tijdsintervallen uit de dataset verwijderen. De auteurs schatten dat het verwijderen van de drie bekende gebeurtenissen in PSR J1713+0747 ongeveer 6 jaar aan data zou vereisen, wat de noodzaak onderstreept om ze te modelleren in plaats van te verwijderen.
• Inzicht in Pulsar-fysica: De bevindingen suggereren dat pulsvormveranderingen vaker voorkomen dan gedacht en kunnen worden veroorzaakt door zowel intrinsieke magnetosferische processen als interstellaire verstrooiing.
• Toekomstige Observaties: De herhaling in PSR B1937+21 suggereert een mogelijk periodiek fenomeen. Het hervatten van observaties bij lage frequenties (800 MHz) met nieuwe ontvangers (zoals de UWBR op de GBT) is essentieel om dit te bevestigen.
• Chromaticiteit als hulpmiddel: De analyse toont aan dat het bestuderen van de frequentie-afhankelijkheid van deze gebeurtenissen helpt bij het onderscheiden tussen ISM-effecten (verstrooiing) en intrinsieke pulsveranderingen.

Kortom, de paper levert een krachtig nieuw instrument aan voor het monitoren van pulsar-stabiliteit, wat direct bijdraagt aan het verminderen van ruis in de zoektocht naar het gravitatiegolfachtergrondsignaal.