The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models

Résumé technique : L'ensemble de données de 15 ans de NANOGrav : Modèles de bruit chromatique personnalisés

Énoncé du problème
Les réseaux de timing de pulsars (PTA) recherchent des ondes gravitationnelles (OG) à des fréquences de l'ordre du nanoherz en mesurant les temps d'arrivée (TOA) des impulsions avec une précision extrême. L'inférence robuste du fond d'ondes gravitationnelles (GWB) nécessite une modélisation complète des sources de bruit non liées aux OG. Une classe critique de ces sources est le bruit chromatique, issu d'effets de propagation interstellaires et interplanétaires tels que les variations de la mesure de dispersion (DM), la diffusion (scattering) et les délais du vent solaire (SW). Ces effets introduisent des délais dépendants de la fréquence ($\Delta t \propto \nu^{-\chi}$) qui, s'ils ne sont pas modélisés ou caractérisés correctement, peuvent fuir dans les paramètres de bruit rouge achromatique. Cette fuite biaise les estimations des paramètres des OG, pouvant imiter ou masquer le signal du GWB. Les analyses précédentes reposaient souvent sur des prescriptions uniformes (par exemple, des modèles DMX à constante par morceaux) qui peuvent échouer à capturer la nature complexe et spécifique à chaque pulsar des effets de propagation, particulièrement l'interaction entre la variabilité du vent solaire et la turbulence du milieu interstellaire (ISM).

Méthodologie
Les auteurs présentent une personnalisation par pulsar des modèles de bruit chromatique pour les 67 pulsars de l'ensemble de données de 15 ans de NANOGrav (NG15). La méthodologie implique :

1. Suite étendue de processus gaussiens (GP) : L'étude utilise une suite étendue de modèles GP pour caractériser simultanément les délais chromatiques. Celles-ci incluent :

   • Modèles dans le domaine de Fourier : Utilisant des densités spectrales de puissance (PSD) en loi de puissance pour la DM et le bruit « Libre Chromatique » (FC).
   • Modèles dans le domaine temporel (TD) : Employant des bases d'interpolation linéaire avec divers noyaux (Exponentiel Quadratique, Quasi-Périodique, Ridge) pour modéliser les variations corrélées dans le temps.
   • Modèles TD-RF : Modèles multidimensionnels permettant aux variations de DM de se décorréler à travers les fréquences radio.
   • Modèles de Vent Solaire (SW) : Une approche hybride combinant un modèle déterministe global à constante par morceaux pour la densité électronique solaire ($n_E$) avec des perturbations stochastiques par processus gaussien (SWGP) pour capturer la variabilité à courte échelle de temps le long de lignes de visée individuelles.

2. Flux de sélection de modèles : Un processus rigoureux de sélection de modèles bayésiens en plusieurs étapes est employé pour chaque pulsar :

   • Criblage initial : Test de trois classes de modèles de base (Fourier, TD, Ridge) pour évaluer la complexité des noyaux et la significativité des composantes FC et SW via les facteurs de Bayes de Savage-Dickey.
   • Optimisation de la taille de la base : Sélection du nombre optimal d'éléments de base (nombre de fréquences $N_f$ ou espacement d'interpolation $dt$) en maximisant le rapport entre l'évidence bayésienne et le temps d'évaluation de la vraisemblance ($Z/t_L$), équilibrant la qualité de l'ajustement et le coût computationnel.
   • Raffinement : Test itératif des variations chromatiques annuelles, des événements transitoires (ex: creux exponentiels) et des variations de DM dépendantes de la bande.
   • Réajustement global : Après la sélection des modèles individuels par pulsar, les paramètres déterministes globaux du vent solaire sont réajustés en utilisant une approche de distribution postérieure factorisée pour assurer la cohérence entre les 67 pulsars.

3. Cadre de comparaison : Les modèles de bruit personnalisés (CNM) sont comparés au modèle de référence standard de NG15 (DMX) et à un modèle basé sur Fourier non personnalisé. Les métriques incluent les changements dans les hyperparamètres du bruit rouge achromatique (RN), les paramètres de bruit blanc (EFAC, EQUAD, ECORR), l'analyse spectrale libre et les courbes de sensibilité.

Contributions clés et résultats

• Modèles de bruit personnalisés : L'étude dérive avec succès des modèles de bruit chromatique uniques et optimisés pour les 67 pulsars. Environ la moitié des pulsars favorisent les modèles temporels (souvent avec des noyaux Ridge ou SE), tandis que l'autre moitié favorise les modèles de domaine de Fourier. Ce partage est corrélé à la présence de vent solaire complexe ou de variations de DM dépendantes de la fréquence, qui sont mieux capturées par les bases temporelles.
• Caractérisation du vent solaire : Le modèle hybride SW suit avec succès le cycle solaire de 11 ans, les variations de densité électronique inférées étant cohérentes avec les données de flux F10.7 cm. L'inclusion de composantes SWGP à grain fin améliore significativement la récupération des indices spectraux de la DM de l'ISM, les déplaçant vers la valeur de Kolmogorov attendue ($\gamma_{DM} \approx 8/3$) et réduisant les biais causés par la variabilité non modélisée du vent solaire.
• Détection du bruit chromatique libre (FC) : Des preuves de bruit chromatique non dispersif (FC) sont trouvées dans 21 des 67 pulsars. Pour six d'entre eux, il s'agit de la première détection rapportée dans les résidus de PTA. Les indices chromatiques ($\chi$) inférés couvrent une large gamme ; si certains s'alignent sur les attentes de diffusion par écran mince ($\chi \approx 4$), d'autres présentent des valeurs élevées ($\chi > 7$) ou faibles ($\chi < 3$), suggérant une physique complexe ou une potentielle erreur de spécification de modèle.
• Impact sur le bruit rouge achromatique : L'application des CNM modifie significativement le bruit rouge achromatique inféré pour 19 des 67 pulsars. Dans 5 cas, le bruit rouge intrinsèque (IRN) auparavant significatif ne l'est plus, indiquant que le signal était probablement d'origine chromatique (ex: diffusion ou effets du vent solaire). Pour d'autres, l'amplitude du bruit rouge diminue ou l'indice spectral change. Notamment, l'amplitude du bruit rouge pour PSR J1713+0747 est considérablement réduite et son indice spectral augmente, sous l'effet de l'inclusion de modèles d'événements chromatiques.
• Bruit blanc et ECORR : Les CNM réduisent généralement le besoin de grands paramètres de bruit corrélé par époque (ECORR), particulièrement pour les pulsars présentant une forte sensibilité à la diffusion ou au vent solaire. Cependant, pour certains pulsars, la significativité de l'ECORR augmente, suggérant que les modèles SW actuels manquent encore de flexibilité pour capturer tous les extrêmes à courte échelle de temps, notamment lors des conjonctions solaires.
• Améliorations de la sensibilité : Les courbes de sensibilité indiquent que les CNM améliorent la sensibilité des PTA, particulièrement aux fréquences moyennes à hautes, en réduisant le « plancher de bruit » et en améliorant la fonction de transmission. Les résidus blanchis restent généralement cohérents avec une distribution gaussienne sous les modèles CNM et DMX, bien que certains pulsars montrent un meilleur ajustement avec les CNM.

Signification
Cette étude démontre qu'une modélisation raffinée et spécifique à chaque pulsar des bruits chromatiques est essentielle pour améliorer la sensibilité et la précision des recherches d'ondes gravitationnelles à basse fréquence. En atténuant la fuite de la puissance chromatique dans le bruit rouge achromatique, ces modèles personnalisés réduisent les biais dans l'inférence du GWB et améliorent la caractérisation du fond stochastique. Les résultats soulignent qu'une prescription de bruit « taille unique » est insuffisante pour la précision requise dans le régime actuel des PTA. Ce travail fournit un cadre robuste pour séparer la structure induite par la propagation des signaux astrophysiques authentiques, soutenant directement les efforts de la collaboration NANOGrav pour détecter et caractériser le fond d'ondes gravitationnelles de l'ordre du nanoherz et les sources potentielles d'ondes continues. De plus, la décomposition des composantes de bruit offre une nouvelle voie pour l'étude de la turbulence de l'ISM, du milieu interplanétaire et de la variabilité de la magnétosphère des pulsars.