The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models

Resumen Técnico: El Conjunto de Datos de 15 años de NANOGrav: Modelos de Ruido Cromático Personalizados

Planteamiento del Problema
Las Matrices de Pulsares de Temporización (PTA, por sus siglas en inglés) buscan ondas gravitacionales (GW) de nanohertz mediante la medición de tiempos de llegada de pulsos (TOAs) ultraprecisos. La inferencia robusta del fondo de ondas gravitacionales (GWB) requiere un modelado exhaustivo de las fuentes de ruido no relacionadas con las GW. Una clase crítica de estas fuentes es el ruido cromático, que surge de efectos de propagación interestelares e interplanetarios, tales como variaciones en la medida de dispersión (DM), la dispersión (scattering) y los retrasos del viento solar (SW). Estos efectos introducen retrasos dependientes de la frecuencia ($\Delta t \propto \nu^{-\chi}$) que, si no se modelan o caracterizan adecuadamente, pueden filtrarse hacia los parámetros de ruido rojo acromático. Esta filtración sesga las estimas de los parámetros de las GW, pudiendo imitar u ocultar la señal del GWB. Los análisis previos a menudo dependían de prescripciones uniformes (por ejemplo, modelos DMX de constante por tramos) que pueden fallar al capturar la naturaleza compleja y específica de cada púlsar de estos efectos de propagación, particularmente la interacción entre la variabilidad del viento solar y la turbulencia del medio interestelar (ISM).

Metodología
Los autores presentan una personalización por púlsar de los modelos de ruido cromático para los 67 púlsares en el conjunto de datos de 15 años de NANOGrav (NG15). La metodología involucra:

1. Suite Expandida de Procesos Gaussianos (GP): El estudio utiliza una suite expandida de modelos GP para caracterizar los retrasos cromáticos de forma simultánea. Estos incluyen:

   • Modelos en el dominio de Fourier: Utilizando densidades espectrales de potencia (PSD) de ley de potencia para la DM y el ruido "Cromático Libre" (FC).
   • Modelos en el dominio del tiempo (TD): Empleando bases de interpolación lineal con varios kernels (Exponencial Cuadrática, Cuasi-Periódico, Ridge) para modelar variaciones correlacionadas en el tiempo.
   • Modelos TD-RF: Modelos multidimensionales que permiten que las variaciones de la DM se desacoplen a través de las frecuencias de radio.
   • Modelos de Viento Solar (SW): Un enfoque híbrido que combina un modelo determinista global de densidad electrónica solar ($n_E$) de constante por tramos con perturbaciones estocásticas de GP (SWGP) para capturar la variabilidad de escala de tiempo corta a lo largo de líneas de visión individuales.

2. Flujo de Trabajo de Selección de Modelos: Se emplea un riguroso proceso de selección de modelos bayesianos de múltiples pasos para cada púlsar:

   • Tamizaje Inicial: Prueba de tres clases de modelos base (Fourier, TD, Ridge) para evaluar la complejidad del kernel y la significancia de los componentes FC y SW utilizando Factores de Bayes de Savage-Dickey.
   • Optimización del Tamaño de la Base: Selección del número óptimo de elementos de base (conteo de frecuencia $N_f$ o espaciado de interpolación $dt$) maximizando la relación entre la evidencia bayesiana y el tiempo de evaluación de la verosimilitud ($Z/t_L$), equilibrando la calidad del ajuste frente al costo computacional.
   • Refinamiento: Pruebas iterativas de variaciones cromáticas anuales, eventos transitorios (por ejemplo, caídas exponenciales) y variaciones de DM dependientes de la banda.
   • Reajuste Global: Tras la selección de los modelos de ruido individuales para cada púlsar, se reajustan los parámetros deterministas globales de SW utilizando un enfoque de posterior factorizado a través de los 67 púlsares para asegurar la consistencia.

3. Marco de Comparación: Los modelos de ruido personalizados (CNM) se comparan contra el estándar de la línea base de NG15 (DMX) y un modelo de Fourier genérico "no personalizado". Las métricas incluyen cambios en los hiperparámetros de ruido rojo (RN) acromático, parámetros de ruido blanco (EFAC, EQUAD, ECORR), análisis espectral libre y curvas de sensibilidad.

Contribuciones Clave y Resultados

• Modelos de Ruido Personalizados: El estudio deriva con éxito modelos de ruido cromático optimizados y únicos para los 67 púlsares. Aproximadamente la mitad de los púlsares favorecen los modelos en el dominio del tiempo (a menudo con kernels Ridge o SE), mientras que la otra mitad favorece los modelos en el dominio de Fourier. Esta división se correlaciona con la presencia de viento solar complejo o variaciones de DM dependientes de la frecuencia, que son mejor capturadas por las bases de TD.
• Caracterización del Viento Solar: El modelo híbrido de SW rastrea con éxito el ciclo solar de 11 años, con variaciones de densidad electrónica inferidas consistentes con los datos de flujo de F10.7 cm. La inclusión de componentes SWGP de grano fino mejora significamente la recuperación de los índices espectrales de la DM del ISM, desplazándolos hacia el valor de Kolmogorov esperado ($\gamma_{DM} \approx 8/3$) y reduciendo los sesgos causados por la variabilidad no modelada del viento solar.
• Detección de Ruido Cromático Libre (FC): Se encuentra evidencia de ruido cromático no dispersivo (FC) en 21 de los 67 púlsares. Para seis de ellos, esta representa la primera detección reportada en los residuales de PTA. Los índices cromáticos ($\chi$) inferidos abarcan un amplio rango; mientras que algunos se alinean con las expectativas de dispersión de pantalla delgada ($\chi \approx 4$), otros exhiben valores grandes ($\chi > 7$) o valores pequeños ($\chi < 3$), lo que sugiere una física compleja o una potencial especificación errónea del modelo.
• Impacto en el Ruido Rojo Acromático: La aplicación de los CNM altera significativamente el ruido rojo acromático inferido para 19 de los 67 púlsares. En 5 casos, el ruido rojo intrínseco (IRN) previamente significativo ya no es significativo, indicando que la señal era probablemente de origen cromático (por ejemplo, efectos de dispersión o viento solar). Para otros, la amplitud del ruido rojo disminuye o el índice espectral cambia. Notablemente, la amplitud del ruido rojo para PSR J1713+0747 se reduce significativamente y su índice espectral aumenta, impulsado por la inclusión de modelos de eventos cromáticos.
• Ruido Blanco y ECORR: Los CNM generalmente reducen la necesidad de grandes parámetros de ruido correlacionado por época (ECORR), particularmente para púlsares con fuerte dispersión o sensibilidad al viento solar. Sin embargo, para algunos púlsares, la significancia de ECORR aumenta, sugiriendo que los modelos actuales de SW aún carecen de la flexibilidad suficiente para capturar todos los valores atípicos de escala de tiempo corta, particularmente cerca de las conjunciones solares.
• Mejoras en la Sensibilidad: Las curvas de sensibilidad indican que los CNM mejoran la sensibilidad de las PTA, particularmente en frecuencias medias a altas, al reducir el "piso de ruido" y mejorar la función de transmisión. Los residuales blanqueados generalmente permanecen consistentes con una distribución gaussiana tanto bajo CNM como bajo DMX, aunque púlsares específicos muestran una mejor calidad de ajuste con los CNM.

Significancia
Este trabajo demuestra que un modelado de ruido cromático refinado y específico para cada púlsar es esencial para mejorar la sensibilidad y la precisión de las búsquedas de ondas gravitacionales de baja frecuencia. Al mitigar la filtración de potencia cromática hacia el ruido rojo acromático, estos modelos personalizados reducen los sesgos en la inferencia del GWB y mejoran la caracterización del fondo estocástico. Los resultados resaltan que una prescripción de ruido de "talla única" es insuficiente para la precisión requerida en el régimen actual de las PTA. El trabajo proporciona un marco robusto para separar la estructura inducida por la propagación de las señales astrofísicas genuinas, apoyando directamente los esfuerzos de la colaboración de la NANOGrav para detectar y caracterizar el fondo de ondas gravitacionales de nanohertz y posibles fuentes de ondas continuas. Además, la descomposición de los componentes de ruido ofrece una nueva vía para estudiar el ISM turbulento, el medio interplanetario y la variabilidad del magnetosfera del púlsar.