A model-independent measurement of the Hubble constant from gravitational-wave standard sirens and electromagnetic observations

Resumen Técnico: Una Medición Independiente de Modelo de la Constante de Hubble a partir de Sirenas Estándar de Ondas Gravitacionales y Observaciones Electromagnéticas

Planteamiento del Problema
El artículo aborda la "tensión de Hubble", una discrepancia significativa de aproximadamente $5\sigma - 6\sigma$ entre las mediciones de la constante de Hubble ($H_0$) del universo temprano derivadas del Fondo Cósmico de Microondas (CMB) bajo el modelo $\Lambda$CDM ($H_0 = 67.4 \pm 0.5$ km s$^{-1}$ Mpc$^{-1}$) y las mediciones directas del universo tardío utilizando supernovas de Tipo Ia (SNe Ia) y la escala de distancia local ($H_0 = 73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$). Los métodos tradicionales suelen depender de modelos cosmológicos específicos o de suposiciones respecto a la escala del horizonte de sonido. Los autores pretenden proporcionar una determinación de $H_0$ totalmente independiente del modelo que evite estos sesgos teóricos, basándose únicamente en datos observacionales e incorporando conjuntos de datos tanto electromagnéticos (EM) como de ondas gravitacionales (GW).

Metodología
El estudio emplea la regresión por Procesos Gaussianos (GP), un método bayesiano no paramétrico, para reconstruir la función de distancia de luminosidad, $D_L(z)$, y su primera derivada directamente a partir de los datos observacionales sin asumir una forma funcional predefinida (por ejemplo, $\Lambda$CDM).

1. Marco de trabajo: Los autores utilizan el código Python Gapp. Modelan $D_L(z)$ como un proceso gaussiano definido por una función de media y un kernel de covarianza. Se seleccionó el kernel exponencial cuadrático sobre los kernels de Matérn debido a su estabilidad y derivadas bien comportadas.
2. Derivación de $H_0$: La constante de Hubble se deriva analíticamente de la función reconstruida utilizando la relación:
   $$H_0 = c \left[ \frac{d}{dz} \left( \frac{D_L(z)}{1+z} \right) \bigg|_{z=0} \right]^{-1}$$
   Esto requiere una reconstrucción precisa tanto de $D_L(z)$ como de su derivada en $z=0$.
3. Integración de Datos: El análisis combina:
   • Ondas Gravitacionales: Datos del catálogo GWTC-3 (90 candidatos con $p_{astro} > 0.5$), que proporcionan distancias de luminosidad absolutas independientes de la escala de distancia cósmica.
   • Sondas Electromagnéticas:
     • Cronómetros Cósmicos (CC): Mediciones directas de $H(z)$ convertidas a $D_L(z)$.
     • DESI BAO DR2: Datos de Oscilaciones Acústicas de Bariones (BAO), donde el horizonte de sonido $r_d$ se trata como un parámetro de molestia (nuisance parameter) y se marginaliza para mantener la independencia del modelo.
     • Supernovas de Tipo Ia: Tres compilaciones (Pantheon+SH0ES, Union 3.0 y DES-SN5YR).
     • Brotes de Rayos Gamma (GRBs): Utilizados como indicadores de alto desplazamiento al rojo mediante la relación de Amati.
4. Gestión de Errores: El método tiene en cuenta los errores de velocidad peculiar en bajos desplazamientos al rojo ($z < 0.05$) sumándolos en cuadratura a los errores observacionales. Crucialmente, el marco de GP calcula y propaga explícitamente la covarianza entre la función reconstruida $D_L(z)$ y su derivada $D'_L(z)$ para asegurar una estimación de incertidumbre estadísticamente consistente.

Contribuciones Clave

• Enfoque Híbrido Independiente de Modelo: Este trabajo es uno de los primeros en integrar los datos de sirenas estándar de GWTC-3 en un marco de GP junto con diversos conjuntos de datos EM para restringir $H_0$ sin depender de la escala del horizonte de sonido o de modelos específicos de energía oscura.
• Cuantificación del Impacto de GW: El estudio aísla explícitamente la contribución de los datos de GW a la precisión de $H_0$. Al comparar las ejecuciones de GP con y sin datos de GW, los autores demuestran que la adición de datos de GWTC-3 reduce la incertidumbre de $1\sigma$ en $H_0$ aproximadamente un 15–20% (de $1.02$a$0.58$ km s$^{-1}$ Mpc$^{-1}$ en una combinación específica).
• Tratamiento de la Covarianza: La implementación de la diferenciación analítica dentro del marco de GP permite el manejo riguroso de las correlaciones entre la función de distancia y su derivada, un factor crítico para minimizar el sesgo en $z=0$.

Resultados
La reconstrucción independiente del modelo arroja valores de $H_0$ que se sitúan entre los resultados de Planck y SH0ES, pero muestran una fuerte preferencia hacia el valor de SH0ES ($H_0 = 73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$).

• Para la mayoría de las combinaciones de datos, la tensión entre el $H_0$ derivado y el resultado de SH0ES es inferior a $2\sigma$, lo que indica que no hay una discrepancia estadísticamente significativa.
• La inclusión de datos de GW estrecha significativamente las restricciones en comparación con los análisis de solo EM.
• Una combinación específica (GWTC3+CC+DESI+DESY5) favorece ligeramente la medición de Planck, pero la tendencia general se alinea con los resultados de la escala de distancia local.

Significado y Reivindicaciones
Los autores afirman modestamente que, si bien este estudio no resuelve definitivamente la tensión de Hubble, proporciona estimaciones de $H_0$ robustas y consistentes que favorecen la determinación de SH0ES utilizando únicamente entradas observacionales. La importancia del trabajo radica en:

1. Validación de los Datos de GW: Demostrar que las sirenas estándar de GW, incluso con sus limitaciones actuales (baja cobertura de desplazamiento al rojo y incertidumbres estadísticas), proporcionan una calibración absoluta que mejora significativamente la precisión de las mediciones de $H_0$ cuando se combinan con datos EM.
2. Verificación de Consistencia: Confirmar que las mediciones de GW y EM son consistentes entre sí dentro del marco independiente del modelo.
3. Utilidad Futura: Establecer una metodología que pueda extenderse a futuros catálogos de GW más grandes de detectores de próxima generación y conjuntos de datos EM más precisos, ofreciendo una vía para investigar el origen de la tensión de Hubble sin sesgo teórico.

El artículo reconoce las limitaciones actuales, como la dependencia de suposiciones de nivel de población para la inferencia del desplazamiento al rojo en eventos de BBH sin contrapartes electromagnéticas, señalando que un tratamiento riguroso de estos sistemas de error sigue siendo una dirección para trabajos futuros.