Dark siren cross-correlations and the sensitivity of $H_0$ to methodological choices

Resumen Técnico: Correlaciones Cruzadas de Sirena Oscura y la Sensibilidad de $H_0$ a las Elecciones Metodológicas

Planteamiento del Problema
Las fuentes de ondas gravitacionales (OG) actúan como "sirenas estándar", proporcionando mediciones de distancia absolutas que ofrecen una vía independiente prometedora para restringir la constante de Hubble ($H_0$) y resolver la tensión actual en los parámetros cosmológicos. Si bien las "sirenas brillantes" (con contrapartes electromagnéticas) y las "sirenas oscuras" (asociando estadísticamente las localizaciones de OG con catálogos de galaxias) han avanzado significativamente, las metodologías de sirena oscura enfrentan desafíos relacionados con catálogos de galaxias incompletos, sesgos sistemáticos y la complejidad computacional de modelar la estructura a gran escala (EAE).

Este artículo investiga la robustez del método de correlación cruzada angular, que correlaciona eventos de OG con catálogos de galaxias a través de su estructura a gran escala compartida, frente a diversas elecciones metodológicas. Específicamente, los autores examinan cómo el tratamiento de la covarianza, la parametrización del sesgo, la agrupación por distancia/corrimiento al rojo y la incompletitud del catálogo impactan el valor inferido de $H_0$. El objetivo es determinar si los sesgos sistemáticos pueden mitigarse eficazmente mediante un modelado adecuado, validando así el enfoque de correlación cruzada para la cosmología de precisión futura.

Metodología
Los autores emplean un marco de correlación cruzada angular tomográfico, analizando el agrupamiento de eventos de OG y galaxias a través de bins de corrimiento al rojo y distancia de luminosidad.

• Formalismo: El método utiliza el espectro de potencia angular cruzado ($C_{\ell}$) entre los trazadores de OG y galaxias. Para reducir el costo computacional y la dimensionalidad, los autores introducen una estadística comprimida, $S_{ij}$, definida como una suma ponderada del espectro de potencia cruzado sobre los multipolos $\ell$. Esta estadística mide eficazmente la correlación cruzada a lo largo de la línea de visión entre los dos trazadores.
• Simulaciones: El análisis se basa en 600 catálogos simulados log-normales generados bajo una cosmología $\Lambda$CDM plana ($H_0 = 70$ km s$^{-1}$ Mpc$^{-1}$, $\Omega_m = 0.25$). Estos simulacros abarcan el rango de corrimiento al rojo $0 \le z \le 0.5$ e incorporan funciones de selección realistas, incluyendo cortes de magnitud e incertidumbres en la localización en el cielo.
• Procedimiento: Las predicciones teóricas se calculan utilizando CAMB bajo la aproximación de Limber, incorporando el modelado directo de los efectos de selección (sesgos de Malmquist y Eddington) directamente en los núcleos de proyección. La inferencia de parámetros se realiza mediante Cadena de Markov de Monte Carlo (MCMC), ajustando $H_0$ y los parámetros de sesgo ($b_{gal}, b_{GW}$).
• Estimación de Covarianza: El estudio compara matrices de covarianza analíticas (derivadas de suposiciones gaussianas) con matrices de covarianza empíricas estimadas a partir de las 600 realizaciones simuladas, aplicando correcciones de Hartlap y Percival para tener en cuenta tamaños de muestra finitos.

Contribuciones y Resultados Clave

1. Estadísticas Comprimidas vs. No Comprimidas:
   Los autores demuestran que comprimir el espectro de potencia angular completo en la estadística resumida $S_{ij}$ retiene la mayoría de la información cosmológica. Si bien el uso de $S_{ij}$ resulta en una pérdida moderada de poder de restricción (aumentando la incertidumbre de $H_0$ en un factor de $\sim 1.4$ para 200 eventos y $\sim 1.5$ para 1000 eventos en comparación con $C_{\ell}$ completo), ofrece ventajas computacionales significativas y estabilidad numérica al evitar la integración directa de funciones de Bessel.

2. Modelado de Covarianza:
   Un hallazgo crítico es la sensibilidad de las restricciones de parámetros a la elección de la matriz de covarianza.

   • Para muestras grandes ($\sim 10^3$ eventos), las covarianzas analíticas y simuladas producen resultados consistentes si se incluyen los términos fuera de la diagonal.
   • Para muestras más pequeñas (200 eventos), la covarianza analítica subestima las incertidumbres en un factor de $\sim 1.5$ en comparación con las covarianzas derivadas de simulaciones. Además, incluir términos fuera de la diagonal en la covarianza analítica para muestras pequeñas conduce a inestabilidad numérica en la inversión de MCMC, impidiendo la convergencia.
   • El estudio concluye que las estimaciones de covarianza basadas en simulaciones son esenciales, particularmente para conteos de eventos más pequeños, para evitar restricciones excesivamente optimistas.

3. Incertidumbres Sistemáticas:

   • Sesgos de Malmquist y Eddington: El artículo confirma que las incertidumbres en la distancia de luminosidad inducen sesgos de Eddington y Malmquist, que pueden sesgar sistemáticamente $H_0$ hacia valores altos si no se modelan. Sin embargo, estos efectos se mitigan eficazmente modelando directamente la probabilidad de detección y los umbrales de selección dentro del núcleo teórico.
   • Parametrización del Sesgo: El estudio prueba parametrizaciones constantes, lineales y cuadráticas para el sesgo de galaxias y de OG. Si bien la parametrización constante introduce un sesgo debido a su incapacidad para capturar la evolución del corrimiento al rojo, las formas lineales y cuadráticas recuperan estimaciones de $H_0$ sin sesgo. Esto sugiere que un conocimiento preciso de la evolución del sesgo no es estrictamente necesario para las observaciones de futuro cercano, siempre que la parametrización sea suficientemente flexible.
   • Incompletitud del Catálogo: Una ventaja mayor identificada es la robustez del método frente a catálogos incompletos. A diferencia de los enfoques estándar de sirena oscura que requieren modelar la población faltante, el método de correlación cruzada solo requiere que la función de selección se modele consistentemente en la predicción teórica. Esto permite utilizar directamente el campo de sobre-densidad observado, obteniendo restricciones de $H_0$ sin sesgo incluso con efectos de selección radial (por ejemplo, límites de magnitud). Sin embargo, combinar catálogos con variabilidad en la completitud a través del cielo puede degradar las restricciones al reducir el poder para restringir parámetros de sesgo globales.

4. Impacto de la Agrupación (Binning):
   La elección de la agrupación por corrimiento al rojo y distancia de luminosidad implica una compensación entre la resolución de la señal y el ruido de disparo.

   • Una agrupación más fina del corrimiento al rojo mejora significativamente las restricciones de $H_0$ porque el ruido de disparo de las galaxias permanece bajo incluso con alta resolución.
   • Una agrupación más fina de la distancia de luminosidad produce rendimientos decrecientes para los eventos de OG debido a su escasez; subdividir estos bins aumenta el ruido de disparo, lo que puede compensar los beneficios de una mayor resolución.
   • La estrategia óptima implica anchos de bin comparables o mayores que la incertidumbre típica de distancia de OG para absorber errores de medición sin una dilución excesiva de la señal.

Significado y Afirmaciones
El artículo afirma que el método de correlación cruzada angular es una herramienta poderosa y complementaria para la cosmología de sirena oscura, capaz de lograr restricciones de $H_0$ sin sesgo siempre que las elecciones metodológicas se gestionen cuidadosamente. Los autores enfatizan que:

• Los sesgos sistemáticos que surgen de incertidumbres de distancia, evolución del sesgo e incompletitud del catálogo pueden mitigarse eficazmente mediante modelado directo en lugar de la reconstrucción explícita de poblaciones faltantes.
• El uso de matrices de covarianza derivadas de simulaciones es crítico para una estimación fiable de incertidumbres, particularmente para muestras de eventos más pequeñas donde las aproximaciones analíticas fallan.
• La estadística comprimida $S_{ij}$ ofrece una alternativa computacionalmente eficiente al análisis completo del espectro de potencia con una pérdida mínima de información.

En última instancia, el trabajo establece que, con un modelado adecuado y una muestra suficientemente grande de eventos de OG precisos, el método de correlación cruzada puede abordar robustamente la tensión de Hubble, ofreciendo una ventaja distinta sobre los enfoques estándar de catálogos de galaxias al incorporar naturalmente los efectos de selección en el marco teórico.