The NANOGrav 15 yr Data Set: Impacts of Customized Chromatic Noise Models on Gravitational Wave Analyses

Résumé technique : L'ensemble de données de 15 ans de NANOGrav : Impacts des modèles de bruit chromatique personnalisés sur les analyses d'ondes gravitationnelles

Énoncé du problème
Les réseaux de chronométrage de pulsars (PTA) reposent sur la stabilité de rotation des pulsars millisecondes pour détecter les ondes gravitationnelles (OG) de l'ordre du nanohertz. Cependant, les ensembles de données des PTA sont fortement contaminés par des sources de bruit, notamment le bruit de spin intrinsèque, les variations du milieu interstellaire (ISM) causant des fluctuations de la mesure de dispersion (DM), les effets du vent solaire et les erreurs de modèle de chronométrage. Une mauvaise modélisation de ces sources de bruit chromatique (dépendantes de la fréquence) peut biaiser l'inférence d'un processus à spectre commun, augmentant ou supprimant potentiellement la signature de corrélation de Hellings-Downs (HD) du fond stochastique d'ondes gravitationnelles (GWB). Les analyses précédentes de l'ensemble de données de 15 ans de NANOGrav (NG15) ont utilisé des modèles de bruit standards (par exemple, DMX pour les variations de DM), mais des travaux récents suggèrent que des modèles de bruit chromatique (CNM) personnalisés pourraient mieux atténuer ces biais systématiques. Cet article étudie comment le remplacement des modèles de bruit standards par des CNM impacte la significativité du GWB, la caractérisation spectrale, ainsi que les recherches d'ondes continues (CW) et de polarisations de gravité alternatives.

Méthodologie
Les auteurs ont réanalysé l'ensemble de données NG15 en utilisant les CNM développés dans Larsen, Baier et al. (2026). Les composantes méthodologiques clés incluent :

• Modélisation du bruit : Le modèle standard « DMX » pour la mesure de dispersion a été remplacé par un modèle de processus gaussien (DMGP). De plus, le cadre CNM teste et applique des signaux chromatiques personnalisés (par exemple, variations chromatiques libres, vent solaire stochastique, événements transitoires) pour chacun des 67 pulsars sur la base de l'évidence bayésienne.
• Implémentation de la vraisemblance : Pour gérer le coût computationnel des modèles de bruit chromatique par pulsar, les auteurs ont implémenté une « Vraisemblance Chromatique Mise en Cache » (Cached-Chromatic Likelihood). Cette technique étend le processus de marginalisation en deux étapes à un processus en $N+2$ étapes, permettant aux paramètres de bruit chromatique (fixés à leurs valeurs de maximum a posteriori issues des analyses de pulsar unique) d'être marginalisés avant l'analyse principale des OG. Cela a réduit les temps d'évaluation de la vraisemblance d'un facteur $\sim2$ (pour le bruit rouge non corrélé commun) à $\sim4$ (pour les corrélations HD) par rapport aux implémentations traditionnelles.
• Portée de l'analyse : L'étude a reproduit les analyses clés de GWB d'Agazie et al. (2023a, 2023d, 2024) et d'Agarwal et al. (2025), incluant :
  • La caractérisation spectrale du GWB (loi de puissance, loi de puissance brisée et spectre libre).
  • La significativité de la corrélation HD (facteurs de Bayes et statistique optimale).
  • Les recherches pour les modes de polarisation scalaire-transverse (ST).
  • Les recherches d'ondes continues (CW) dans tout le ciel et ciblées.
  • Les interprétations astrophysiques et cosmologiques du GWB.

Contributions et résultats clés

1. Significativité accrue de la corrélation HD :

   • Le facteur de Bayes pour les corrélations HD sur un processus de bruit rouge commun non corrélé (CURN) a augmenté d'environ $8\times$, passant de $\sim185$ (dans Agazie et al. 2023a) à $1571 \pm 14$ en utilisant les CNM avec 14 modes de Fourier.
   • L'analyse fréquentiste utilisant la statistique optimale a montré une amélioration du rapport signal sur bruit (S/N) de $5.0$à$5.3$.
   • La valeur $p$ bayésienne est passée de $7.9 \times 10^{-4}$ ($3.16\sigma$) à $4.5 \times 10^{-4}$ ($3.32\sigma$).
   • Les auteurs attribuent cette augmentation à une meilleure atténuation du bruit, particulièrement dans la plage de fréquences 10–30 nHz, plutôt qu'à l'ajout de nouvelles données.

2. Caractérisation spectrale :

   • L'amplitude du GWB inférée a diminué à $A_{\text{GWB}} = 2.1^{+0.6}_{-0.5} \times 10^{-15}$ (à $\gamma = 13/3$ fixé), soit une réduction d'environ $12\%$ par rapport aux estimations précédentes.
   • L'indice spectral a augmenté à $\gamma_{\text{GWB}} = 3.5^{+0.7}_{-0.6}$ dans les analyses à $\gamma$ variable, se rapprochant de l'attente théorique de $13/3$ pour les binaires de trous noirs supermassifs (SMBHB).
   • Une excursion significative précédemment rapportée dans le spectre de puissance à 16 nHz (dans Agazie et al. 2025a) a été identifiée comme un artefact dépendant du modèle de bruit lié à PSR J1713+0747 et n'est plus significative sous les CNM.

3. Polarisations alternatives et ondes continues :

   • L'évidence pour un mode scalaire-transverse (ST) non einsteinien a été réduite. Le facteur de Bayes favorisant HD sur ST est passé d'environ $\sim2.3$ à $\sim158$, et la limite supérieure à 95 % de l'amplitude ST a diminué d'environ $26\%$.
   • Dans les recherches de CW dans tout le ciel, la significativité d'un candidat à basse fréquence ($\sim4.2$ nHz) est restée statistiquement insignifiante lorsqu'elle est modélisée contre un fond corrélé HD ($B_{\text{HD+CW}}^{\text{HD}} = 1.69 \pm 0.06$).
   • Les recherches de CW ciblées pour des candidats AGN spécifiques (SDSS J0729+4008 et SDSS J1536+0441) ont montré une significativité réduite, avec des facteurs de Bayes tombant à des niveaux cohérents avec la population générale de candidats.
   • Le volume de détection du réseau NANOGrav a augmenté d'un facteur $3.2$ pour les recherches de CW grâce à la réduction des planchers de bruit blanc aux fréquences moyennes à hautes.

4. Budget de bruit rouge et sensibilité du détecteur :

   • Le nombre de pulsars présentant un bruit rouge intrinsèque (IRN) significatif aux côtés du processus commun est passé de 12 (dans Agazie et al. 2023b) à 6. Cela suggère que les détections précédentes d'IRN étaient probablement dues à un bruit chromatique mal modélisé.
   • Les courbes de sensibilité indiquent que les CNM procurent des gains substantiels de sensibilité aux fréquences moyennes à hautes, déplaçant la bande de fréquence la plus sensible de $\sim5$ nHz (DMX) à $\sim6.5$ nHz.

5. Interprétations astrophysiques et cosmologiques :

   • Astrophysique (SMBHB) : La caractérisation spectrale mise à jour est cohérente avec une origine SMBHB. Les ajustements aux modèles de population phénoménologiques ont montré des changements seulement marginaux, avec une légère préférence pour des masses de galaxies caractéristiques et des pentes de durcissement interne plus faibles.
   • Cosmologique : Les réinterprétations impliquant l'inflation cosmique, les transitions de phase du premier ordre et les cordes cosmiques stables ont montré des décalages quantitatifs très légers dans les posteriors des paramètres. Pour les cordes cosmiques, la tension la mieux ajustée $G\mu$ a diminué d'environ $35\%$, impliquant une échelle d'énergie plus basse pour la transition de phase produisant les cordes.

Signification et affirmations
L'article affirme que l'impact principal des CNM est l'amélioration de l'atténuation du bruit chromatique, ce qui conduit à une détection plus robuste et statistiquement significative de la signature de corrélation HD. Les auteurs déclarent que la caractérisation spectrale du GWB présente « moins de tension » avec les modèles de population de galaxies et que l'amplitude est légèrement plus faible et plus raide.

Crucialement, les auteurs soulignent que bien que l'évidence statistique pour le GWB se soit renforcée, les interprétations astrophysiques et cosmologiques globales restent largement inchangées car elles sont actuellement « dominées par les priors ». L'article conclut que la modélisation détaillée des effets chromatiques est impérative pour les futurs ensembles de données afin d'optimiser la sensibilité des détecteurs et d'assurer une caractérisation précise du fond d'ondes gravitationnelles de l'ordre du nanohertz, mais il ne prétend pas avoir résolu toutes les tensions entre le signal GWB et les modèles théoriques. Ce travail sert de validation du cadre CNM pour les futures analyses de NANOGrav.