Revisiting GW170817 at milliarcsecond scale: high-precision constraints on jet geometry and $H_0$

Resumo Técnico: Revisão de GW170817 em Escala de Milissegundos

Declaração do Problema
A detecção da fusão de estrelas de nêutrons binárias GW170817 forneceu a primeira oportunidade para uma medição de "sirene padrão" da constante de Hubble ($H_0$). No entanto, a precisão dessa medição é fundamentalmente limitada pela forte covariância entre a distância de luminosidade ($D_L$) e o ângulo de inclinação ($\iota$) do sistema binário. Embora os dados de ondas gravitacionais (GW) isoladamente não possam restringir com precisão a inclinação para eventos de face-on como GW170817, o jato relativístico lançado pela fusão oferece uma alça eletromagnética para essa geometria. Tentativas anteriores de restringir o ângulo de visão ($\theta_v$) usando a curva de luz do afterglow foram limitadas por degenerescências geométricas entre o ângulo de visão e o ângulo de abertura do jato ($\theta_{cp}$). Por outro lado, observações de Interferometria de Linha de Base Muito Longa (VLBI) do movimento próprio superluminal do jato fornecem restrições ortogonais, mas historicamente dependeram do ajuste de modelos de jatos discretos e semi-analíticos ou subconjuntos limitados do espaço de parâmetros, frequentemente negligenciando informações espaciais completas (como o tamanho da fonte) ou incertezas sistemáticas de maneira rigorosa.

Metodologia
Este estudo apresenta uma estrutura abrangente de ajuste de modelo no plano de visibilidade Bayesiano, projetada para extrair restrições robustas sobre a geometria de visão de GW170817 e, subsequentemente, $H_0$. A metodologia integra os seguintes componentes:

1. Integração de Dados: A análise incorpora todos os quatro conjuntos de dados astrométricos disponíveis em escala de milissegundos:

   • Duas épocas do VLBI da High Sensitivity Array (HSA) (Dias 75 e 230 pós-fusão) de Mooley et al. (2018).
   • Uma época do Global-VLBI (Dia 207) de Ghirlanda et al. (2019), que os autores reduziram independentemente para garantir consistência com a calibração do HSA.
   • Uma medição do Telescópio Espacial Hubble (HST) da posição da kilonova (Dia 8), servindo como âncora da fusão (Dia 0).

2. Modelagem Hidrodinâmica: Em vez de depender de modelos pré-computados discretos, os autores empregam um formalismo de escala aplicado a um contínuo de novos modelos de afterglow de jatos hidrodinâmicos baseados em simulações relativísticas 2-D (Govreen-Segal & Nakar 2023). Esses modelos cobrem quatro distribuições de energia de jato (top-hat e leis de potência com índices $b=3, 4, 6$). A escala permite que os modelos sejam ajustados para ângulos de visão arbitrários ($\theta_v$), ângulos de abertura ($\theta_{cp}$) e tempos do observador, utilizando a invariância dos modelos sob razões específicas de parâmetros.

3. Inferência Bayesiana: A estrutura opera diretamente nas visibilidades observadas do interferômetro de rádio usando o software acelerado por GPU GALARIO e o amostrador aninhado Dynesty (via Bilby). Essa abordagem evita a perda de informação associada à deconvolução no plano de imagem. O modelo inclui parâmetros de incômodo para lidar rigorosamente com incertezas sistemáticas, incluindo deslocamentos astrométricos ($\delta_{ra}, \delta_{dec}$) e incertezas na escala de fluxo ($\delta_F$) para cada época.

4. Restrições Cosmológicas: Para medir $H_0$, a estrutura incorpora verossimilhanças de GW (via um estimador de densidade de kernel de inclinação e distância) e correções de velocidade peculiar. Os autores consideram 28 combinações de catálogos de grupos de galáxias e mapas de reconstrução de velocidade peculiar, empregando Média de Modelos Bayesianos (BMA) para marginalizar sobre essas incertezas sistemáticas.

Principais Contribuições

• Espaço de Parâmetros Contínuo: O estudo avança além do ajuste de modelos discretos para explorar uma faixa contínua de geometrias de jato ($\theta_v$ e $\theta_{cp}$), evitando assim vieses introduzidos pela seleção de estruturas de jato específicas.
• Ajuste no Plano de Visibilidade: Ao ajustar diretamente às visibilidades em vez de centróides no plano de imagem, o método utiliza a informação espacial completa dos dados, incluindo restrições de tamanho da fonte, o que é crítico para quebrar degenerescências.
• Tratamento de Incertezas Sistemáticas: A inclusão de parâmetros de incômodo para incertezas de calibração astrométrica e de fluxo, juntamente com a abordagem BMA para correções de velocidade peculiar, fornece um orçamento de erro mais robusto do que estudos anteriores que frequentemente assumiam correções fixas.
• Redução Independente do Global-VLBI: Os autores realizaram uma redução independente dos dados do Global-VLBI para garantir consistência de quadro com os dados do HSA, identificando e corrigindo deslocamentos astrométricos relativos a reduções anteriores.

Resultados

• Restrições de Ângulo de Visão:
  • Ao assumir uma distância de luminosidade fixa de $D_L = 40,7$ Mpc (baseada em medições de flutuação de brilho superficial), a análise produz um ângulo de visão de $\theta_v = 18,3^\circ - 20,3^\circ$ (intervalo de credibilidade de 68%) ao incorporar priors informados pela curva de luz sobre a geometria do jato.
  • Ao marginalizar sobre $D_L$ e correções de velocidade peculiar, o ângulo de visão preferido muda ligeiramente para $\theta_v = 16,8^\circ - 19,2^\circ$.
• Distância de Luminosidade: O ajuste direto para a distância produz $D_L = 44,0 \pm 1,6$ Mpc. Isso é aproximadamente $1\sigma$ maior do que o valor de $40,7 \pm 2,4$ Mpc derivado do método de flutuação de brilho superficial, destacando a sensibilidade das restrições de ângulo de visão à distância assumida.
• Constante de Hubble ($H_0$): A medição média de modelos bayesianos através de 28 cenários de correção de velocidade peculiar é $H_0 = 65,5 \pm 4,4$ km s$^{-1}$ Mpc$^{-1}$.
  • Este valor está dentro de $0,5\sigma$ do valor da Colaboração Planck do Universo primordial ($67,4 \pm 0,5$ km s$^{-1}$ Mpc$^{-1}$).
  • É $1,7\sigma$ menor do que a medição do Universo tardio SH0ES ($73,04 \pm 1,04$ km s$^{-1}$ Mpc$^{-1}$).

Significado e Alegações
O artigo afirma que esta medição independente do Universo tardio de $H_0$ é mais consistente com o resultado do Planck do Universo primordial do que com o resultado do SH0ES. Os autores sugerem que, se tomado ao pé da letra, essa descoberta poderia desafiar a necessidade de modelos de cosmologia não padrão para resolver a tensão de Hubble. No entanto, os autores permanecem modestos quanto a essa alegação, notando explicitamente várias ressalvas:

• Especificação Incorreta do Modelo: A evolução em tempos tardios da curva de luz do afterglow ($\gtrsim 350$ dias) não é bem ajustada por modelos de síncrotron padrão, e a fonte da emissão excessiva permanece desconhecida. Embora os autores argumentem que sua análise (focada em épocas anteriores) é dominada pela geometria e não por imprecisões do modelo, eles reconhecem que física não modelada poderia enviesar os resultados.
• Alinhamento do Jato: A análise assume que o jato é lançado ao longo do eixo do momento angular orbital. Se existir um desalinhamento significativo, isso poderia enviesar não linearmente a estimativa de $H_0$.
• Domínio da Velocidade Peculiar: Os autores identificam as correções de velocidade peculiar como a principal fonte de incerteza na medição de $H_0$ (contribuindo com ~5% de incerteza fracionária), em vez da restrição do ângulo de visão em si.

O estudo conclui que, embora observações VLBI de alta precisão de afterglows de jatos permaneçam uma ferramenta única e poderosa para restringir ângulos de inclinação a curto prazo, alcançar uma resolução definitiva para a tensão de Hubble exigirá uma amostra maior de sirenes brilhantes e modelagem aprimorada da física dos jatos e das velocidades peculiares.