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

Resumen Técnico: El conjunto de datos de 15 años de NANOGrav: Impactos de los modelos de ruido cromático personalizados en los análisis de ondas gravitacionales

Planteamiento del problema
Las matrices de temporización de púlsares (PTA) dependen de la estabilidad rotacional de los púlsares de milisegundos para detectar ondas gravitacionales (GW) de nanohercios. Sin embargo, los conjuntos de datos de las PTA están fuertemente contaminados por fuentes de ruido, incluyendo el ruido de rotación intrínseco, variaciones en la medida de dispersión (DM) causadas por el medio interestelar, efectos del viento solar y errores en el modelo de temporización. El modelado incorrecto de estas fuentes de ruido cromático (dependientes de la frecuencia) puede sesgar la inferencia de un proceso de espectro común, potencialmente inflando o suprimiendo la firma de correlación de Hellings-Downs (HD) del fondo de ondas gravitacionales estocásticas (GWB). Los análisis previos del conjunto de datos de 15 años de NANOGrav (NG15) utilizaron modelos de ruido estándar (por ejemplo, DMX para variaciones de DM), pero trabajos recientes sugieren que los modelos de ruido cromático (CNM) personalizados pueden mitigar mejor estos sistemas. Este artículo investiga cómo la sustitución de los modelos de ruido estándar por los CNM impacta en la significación del GWB, la caracterización espectral y las búsquedas de ondas continuas (CW) y polarizaciones gravitacionales alternativas.

Metodología
Los autores reanalizaron el conjunto de datos NG15 utilizando los CNM desarrollados en Larsen, Baier et al. (2026). Los componentes metodológicos clave incluyen:

• Modelado de ruido: El modelo estándar "DMX" para la medida de dispersión fue sustituido por un modelo de Proceso Gaussiano (DMGP). Además, el marco de los CNM pone a prueba y aplica señales cromáticas personalizadas (por ejemplo, variaciones cromáticas libres, viento solar estocástico, eventos transitorios) para cada uno de los 67 púlsares basándose en la evidencia bayesiana.
• Implementación de la verosimilitud: Para gestionar el coste computacional de los modelos cromáticos por púlsar, los autores implementaron una "Verosimilitud Cromática Almacenada" (Cached-Chromatic Likelihood). Esta técnica extiende el proceso de marginalización de dos pasos a un proceso de $N+2$ pasos, permitiendo que los parámetros de ruido cromático (fijados a sus valores de máxima a posteriori de los análisis de un solo púlsar) se marginalicen antes del análisis principal de GW. Esto redujo los tiempos de evaluación de la verosimilitud en factores de $\sim2$ (para ruido rojo común no correlacionado) a $\sim4$ (para correlaciones HD) en comparación con las implementaciones tradicionales.
• Alcance del análisis: El estudio reprodujo análisis clave de GW de Agazie et al. (2023a, 2023d, 2024) y Agarwal et al. (2025), incluyendo:
  • Caracterización espectral del GWB (ley de potencia, ley de potencia rota y espectro libre).
  • Significación de la correlación HD (factores de Bayes y Estadística Óptima).
  • Búsquedas de modos de polarización escalar-transversa (ST).
  • Búsquedas de ondas continuas (CW) de todo el cielo y dirigidas.
  • Interpretaciones astrofísicas y cosmológicas del GWB.

Contribuciones y resultados clave

1. Significación mejorada de la correlación HD:

   • El factor de Bayes para las correlaciones HD sobre un proceso de ruido rojo común no correlacionado (CURN) aumentó aproximadamente $8\times$, de $\sim185$ (en Agazie et al. 2023a) a $1571 \pm 14$ usando los CNM con 14 modos de Fourier.
   • El análisis frecuentista utilizando la Estadística Óptima mostró una mejora en la relación señal-ruido (S/N) de $5.0$a$5.3$.
   • El valor $p$ bayesiano disminuyó de $7.9 \times 10^{-4}$ ($3.16\sigma$) a $4.5 \times 10^{-4}$ ($3.32\sigma$).
   • Los autores atribuyen este aumento a una mejor mitigación del ruido, particularmente en el rango de frecuencia de 10–30 nHz, en lugar de la adición de nuevos datos.

2. Caracterización espectral:

   • La amplitud inferida del GWB disminuyó a $A_{\text{GWB}} = 2.1^{+0.6}_{-0.5} \times 10^{-15}$ (con $\gamma = 13/3$ fijo), una reducción de $\sim12\%$ respecto a las estimaciones previas.
   • El índice espectral aumentó a $\gamma_{\text{GWB}} = 3.5^{+0.7}_{-0.6}$ en los análisis de $\gamma$ variable, acercándose a la expectativa teórica de $13/3$ para binarias de agujeros negros supermasivos (SMBHB).
   • Una excursión significativa reportada anteriormente en el espectro de potencia a 16 nHz (en Agazie et al. 2025a) resultó ser un artefacto dependiente del modelo de ruido vinculado a PSR J1713+0747 y ya no es significativo bajo los CNM.

3. Polarizaciones alternativas y ondas continuas:

   • La evidencia para un modo escalar-transverso (ST) no Einsteiniano se redujo. El factor de Bayes que favorece a HD sobre ST aumentó de $\sim2.3$ a $\sim158$, y el límite superior del 95% para la amplitud ST disminuyó en un $\sim26\%$.
   • En las búsquedas de CW de todo el cielo, la significación de un candidato de baja frecuencia ($\sim4.2$ nHz) siguió siendo estadísticamente insignificante cuando se modeló contra un fondo correlacionado HD ($B_{\text{HD+CW}}^{\text{HD}} = 1.69 \pm 0.06$).
   • Las búsquedas de CW dirigidas para candidatos específicos de AGN (SDSS J0729+4008 y SDSS J1536+0441) mostraron una reducción de significación, con los factores de Bayes cayendo a niveles consistentes con la población general de candidatos.
   • El volumen de detección de la matriz NANOGrav aumentó en un factor de $3.2$ para las búsquedas de CW debido a la reducción de los suelos de ruido blanco en frecuencias medias-altas.

4. Presupuesto de ruido rojo y sensibilidad del detector:

   • El número de púlsares que exhiben ruido rojo intrínseco (IRN) significativo junto al proceso común cayó de 12 (en Agazie et al. 2023b) a 6. Esto sugiere que las detecciones previas de IRN eran probablemente ruido cromático mal modelado.
   • Las curvas de sensibilidad indican que los CNM proporcionan ganancias sustanciales de sensibilidad en frecuencias medias-altas, desplazando la banda de frecuencia más sensible de $\sim5$ nHz (DMX) a $\sim6.5$ nHz.

5. Interpretaciones astrofísicas y cosmológicas:

   • Astrofísica (SMBHB): La actualización de la caracterización espectral es consistente con un origen de SMBHB. Los ajustes a modelos de población fenomenológica mostraron solo cambios marginales, con una ligera preferencia por masas características galácticas y pendientes de endurecimiento interno más bajas.
   • Cosmológica: Las reinterpretaciones relacionadas con la inflación cósmica, transiciones de fase de primer orden y cuerdas cósmicas estables mostraron solo ligeros cambios cuantitativos en los posteriors de los parámetros. Para las cuerdas cósmicas, la tensión de mejor ajuste $G\mu$ disminuyó un $\sim35\%$, lo que implica una escala de energía inferior para la transición de fase que produce las cuerdas.

Significación y afirmaciones
El artículo afirma que el impacto principal de los CNM es la mejora en la mitigación del ruido cromático, lo que conduce a una detección más robusta y estadísticamente significativa de la firma de correlación HD. Los autores declaran que la caracterización espectral del GWB tiene ahora "menos tensión" con los modelos de población galáctica y que la amplitud es ligeramente más silenciosa y empinada.

Crucialmente, los autores enfatizan que, si bien la evidencia estadística para el GWB se ha fortalecido, las interpretaciones astrofísicas y cosmológicas generales permanecen en gran medida sin cambios porque están actualmente "dominadas por la prioridad" (prior dominated). El trabajo concluye que el modelado detallado de los efectos cromáticos es imperativo para futuros conjuntos de datos para optimizar la sensibilidad del detector y asegurar la caracterización precisa del fondo de ondas gravitacionales en la escala de los nanohercios, pero no pretende haber resuelto todas las tensiones entre la señal del GWB y los modelos teóricos. El trabajo sirve como una validación del marco de los CNM para futuros análisis de NANOGrav.