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

Resumo Técnico: Uma Medição Independente de Modelo da Constante de Hubble a partir de Sirenes Gravitacionais de Ondas Gravitacionais e Observações Eletromagnéticas

Problema
O artigo aborda a "tensão de Hubble", uma discrepância significativa de aproximadamente $5\sigma - 6\sigma$ entre as medições do universo primordial da constante de Hubble ($H_0$) derivadas do Fundo Cósmico de Micro-ondas (CMB) sob o modelo $\Lambda$CDM ($H_0 = 67,4 \pm 0,5$ km s$^{-1}$ Mpc$^{-1}$) e as medições diretas do universo tardio utilizando Supernovas do Tipo Ia (SNe Ia) e a escada de distâncias local ($H_0 = 73,04 \pm 1,04$ km s$^{-1}$ Mpc$^{-1}$). Os métodos tradicionais frequentemente dependem de modelos cosmológicos específicos ou de suposições relativas à escala do horizonte de som. Os autores visam fornecer uma determinação de $H_0$ totalmente independente de modelo que evite esses vieses teóricos, baseando-se exclusivamente em dados observacionais, incorporando conjuntos de dados eletromagnéticos (EM) e de ondas gravitacionais (GW).

Metodologia
O estudo emprega a regressão por Processo Gaussiano (GP), um método bayesiano não paramétrico, para reconstruir a função de distância de luminosidade, $D_L(z)$, e sua primeira derivada diretamente dos dados observacionais sem assumir uma forma funcional predefinida (ex: $\Lambda$CDM).

1. Estrututura: Os autores utilizam o código Python GaPP. Eles modelam $D_L(z)$ como um processo gaussiano definido por uma função média e um kernel de covariância. O kernel exponencial quadrático foi selecionado em detrimento dos kernels de Matérn devido à sua estabilidade e derivadas bem comportadas.
2. Derivação de $H_0$: A constante de Hubble é derivada analiticamente da função reconstruída usando a relação:
   $$H_0 = c \left[ \frac{d}{dz} \left( \frac{D_L(z)}{1+z} \right) \bigg|_{z=0} \right]^{-1}$$
   Isso requer a reconstrução precisa de tanto $D_L(z)$ quanto de sua derivada em $z=0$.
3. Integração de Dados: A análise combina:
   • Ondas Gravitacionais: Dados do catálogo GWTC-3 (90 candidatos com $p_{astro} > 0,5$), fornecendo distâncias de luminosidade absolutas independentes da escada de distâncias cósmicas.
   • Sondas Eletromagnéticas:
     • Cronômetros Cósmicos (CC): Medições diretas de $H(z)$ convertidas em $D_L(z)$.
     • DESI BAO DR2: Dados de Oscilações Acústicas de Bárions (BAO), onde o horizonte de som $r_d$ é tratado como um parâmetro de confusão (nuisance parameter) e marginalizado para manter a independência do modelo.
     • Supernovas do Tipo Ia: Três compilações (Pantheon+SH0ES, Union 3.0 e DES-SN5YR).
     • Explosões de Raios Gama (GRBs): Utilizadas como indicadores de alto redshift via relação de Amati.
4. Tratamento de Erros: O método contabiliza erros de velocidade peculiar em baixos redshifts ($z < 0,05$) ao adicioná-los em quadratura aos erros observacionais. Crucialmente, o framework de GP calcula e propaga explicitamente a covariância entre a função reconstruída $D_L(z)$ e sua derivada $D'_L(z)$ para garantir uma estimativa de incerteza estatisticamente consistente.

Principais Contribuições

• Abordagem Híbrida Independente de Modelo: Este trabalho é um dos primeiros a integrar dados de sirenes padrão GWTC-3 em um framework de GP junto com diversos conjuntos de dados EM para restringir $H_0$ sem depender da escala do horizonte de som ou de modelos específicos de energia escura.
• Quantificação do Impacto de GW: O estudo isola explicitamente a contribuição dos dados de GW para a precisão de $H_0$. Ao comparar execuções de GP com e sem dados de GW, os autores demonstram que a adição dos dados GWTC-3 reduz a incerteza de $1\sigma$ em $H_0$ em aproximadamente 15–20% (de $1,02$ para $0,58$ km s$^{-1}$ Mpc$^{-1}$ em uma combinação específica).
• Tratamento de Covariância: A implementação da diferenciação analítica dentro do framework de GP permite o tratamento rigoroso das correlações entre a função de distância e sua derivada, um fator crítico para minimizar o viés em $z=0$.

Resultados
A reconstrução independente de modelo produz valores de $H_0$ que situam-se entre os resultados de Planck e SH0ES, mas mostram uma forte preferência pelo valor de SH0ES ($H_0 = 73,04 \pm 1,04$ km s$^{-1}$ Mpc$^{-1}$).

• Para a maioria das combinações de dados, a tensão entre o $H_0$ derivado e o resultado SH0ES é inferior a $2\sigma$, indicando que não há uma discrepância estatisticamente significativa.
• A inclusão de dados de GW torna as restrições significativamente mais estreitas em comparação com as análises apenas de EM.
• Uma combinação específica (GWTC3+CC+DESI+DESY5) favorece ligeiramente a medição de Planck, mas a tendência geral alinha-se com os resultados da escada de distâncias local.

Significância e Alegações
Os autores alegam modestamente que, embora este estudo não resolva definitivamente a tensão de Hubble, ele fornece estimativas de $H_0$ robustas e consistentes que favorecem a determinação da SH0ES usando apenas entradas observacionais. A significância do trabalho reside em:

1. Validação de Dados de GW: Demonstrar que as sirenes padrão de GW, mesmo com as limitações atuais (baixa cobertura de redshift e incertezas estatísticas), fornecem uma calibração absoluta que melhora significativamente a precisão das medições de $H_0$ quando combinadas com dados EM.
2. Verificação de Consistência: Confirmar que as medições de GW e EM são consistentes entre si dentro do framework independente de modelo.
3. Utilidade Futura: Estabelecer uma metodologia que pode ser estendida para futuros catálogos de GW maiores de detectores de próxima geração e conjuntos de dados EM mais precisos, oferecendo um caminho para investigar a origem da tensão de Hubble sem viés teórico.

O artigo reconhece limitações atuais, como a dependência de suposições de nível de população para a inferência de redshift em eventos BBH sem contrapartidas eletromagnéticas, observando que um tratamento rigoroso desses sistemas é uma direção para trabalhos futuros.