Joint Electromagnetic and Gravitational Wave Inference of Binary Neutron Star Merger GW170817 Using Forward-Modeling Ejecta Predictions

Resumo Técnico: Inference Conjunta de Ondas Gravitacionais e Eletromagnéticas para GW170817

1. O Problema

O evento GW170817 (fusão de estrelas de nêutrons binárias) e sua contraparte eletromagnética AT2017gfo (kilonova) inauguraram a era da astronomia multimensageira. O objetivo central é utilizar dados combinados de ondas gravitacionais (GW) e luz (EM) para inferir propriedades fundamentais do sistema, como as massas das estrelas, seus raios e a Equação de Estado (EOS) da matéria densa.

No entanto, existe um elo crítico e incerto nessa cadeia de inferência: a mapeamento entre as propriedades do sistema binário (inferidas das GW) e as propriedades do material ejetado (que geram a kilonova).

• Estudos anteriores utilizaram "fórmulas de ajuste" (fitting formulae) derivadas de simulações de relatividade numérica para prever a massa e a velocidade do material ejetado com base nos parâmetros binários.
• O problema identificado pelos autores é que essas fórmulas possuem incertezas sistemáticas significativas e são frequentemente aplicadas fora de suas regiões de calibração (extrapolação).
• Muitas análises anteriores trataram a velocidade do ejetado como um parâmetro livre sem restrições físicas ou ignoraram as incertezas sistemáticas desses modelos de ajuste, levando a conclusões potencialmente enviesadas sobre a natureza das estrelas de nêutrons.

2. Metodologia

Os autores realizaram uma análise bayesiana conjunta rigorosa, combinando dados de GW e EM, utilizando um modelo de substituição (surrogate model) avançado para a kilonova.

• Modelo de Kilonova (EM): Utilizaram um modelo de dois componentes (ejetado dinâmico rico em lantanídeos e ejetado de vento pobre em lantanídeos) simulado pelo código SuperNu. Para acelerar a inferência, empregaram um surrogate model baseado em Redes Neurais Perceptron Multicamada (MLP), treinado em ~250 simulações de curvas de luz, capaz de interpolar espectros e curvas de luz para filtros observacionais (Rubin e 2MASS).
• Modelos de Ajuste de Ejetado (Forward Models): Testaram três fórmulas de ajuste distintas baseadas em diferentes conjuntos de simulações de relatividade numérica:
  1. KruFo20: Baseado em 52 simulações.
  2. DiCo20: Baseado em 73 simulações.
  3. Nedora21: Baseado em 324 simulações (o maior conjunto).
• Inferência Bayesiana: Utilizaram o framework RIFT para calcular a verossimilhança conjunta ($L_{GW} \times L_{EM}$).
  • Parâmetros Livres: Além das massas, spins e raios das estrelas, introduziram parâmetros de incerteza sistemática para quantificar as falhas nos modelos de ajuste:
    • $\alpha_{dyn}$: Fator de escala para a massa do ejetado dinâmico.
    • $f_{disk}$: Fração do disco de acreção que é ejetada como vento.
    • $\beta_{\phi}$: Parâmetro de incerteza para a velocidade do ejetado.
    • $\sigma_{sys}$: Incerteza sistemática geral nas curvas de luz.
• Abordagem Comparativa: Realizaram três cenários de análise:
  1. Todas as incertezas sistemáticas livres (cenário realista).
  2. Sem incerteza na velocidade ($\beta_{\phi} = 0$).
  3. Sem incerteza na velocidade e fração de disco fixa ($f_{disk} = 30\%$).

3. Principais Contribuições e Resultados

A. Dependência Crítica do Modelo de Ajuste

• Ao incluir as incertezas sistemáticas, os três modelos de ajuste (KruFo20, DiCo20, Nedora21) produzem conclusões distintas e incompatíveis sobre os parâmetros binários e a EOS.
• Enquanto o modelo de luz (EM) sozinho restringe bem as propriedades do ejetado, a adição dos dados de GW não resolve a tensão entre os modelos de ajuste. Cada modelo de ajuste favorece uma região diferente do espaço de parâmetros (ex: raios de estrelas de nêutrons variando de 8 km a 40 km dependendo do modelo escolhido).
• Isso demonstra que a inferência sobre o raio da estrela de nêutrons e a EOS é inteiramente dependente da escolha do modelo de mapeamento (forward model) e não apenas dos dados observacionais.

B. O Papel da Velocidade do Ejetado

• A velocidade do ejetado é um fator determinante. Quando a incerteza na velocidade é ignorada ($\beta_{\phi} = 0$), o modelo força soluções extremas para reconciliar os dados, favorecendo binários altamente assimétricos com raios de estrelas de nêutrons muito grandes ($R_{1.4} > 15$ km) ou simétricos com raios extremos.
• Isso contradiz inferências anteriores que marginalizavam a velocidade sem considerar a tensão física imposta pelos modelos de ajuste.

C. Fração de Ejeção do Disco

• Os resultados indicam que apenas uma pequena fração ($\sim 4-20\%$) da massa do disco de acreção é ejetada como vento para explicar AT2017gfo.
• Isso está em tensão com suposições comuns na literatura (que frequentemente assumem até 40% de ejeção), mas é consistente com simulações de magnetohidrodinâmica relativística (GRMHD) de curta duração com transporte de neutrinos completo.

D. Inconsistência dos Modelos Atuais

• Os autores mostram que os modelos de ajuste atuais são "superconstrangedores" (fazem mais previsões do que têm parâmetros livres) e estão em tensão uns com os outros.
• A aplicação desses modelos fora de suas regiões de calibração (necessária para cobrir o espaço de parâmetros de GW) introduz vieses significativos que não podem ser corrigidos apenas aumentando a incerteza sistemática.

4. Significado e Implicações

• Alerta sobre Inferências Multimensageiras: O estudo alerta que conclusões robustas sobre a Equação de Estado da matéria densa e propriedades de estrelas de nêutrons a partir de eventos como GW170817 são atualmente limitadas pela falta de um mapeamento físico preciso e unificado entre as simulações de relatividade numérica e os modelos de kilonova.
• Necessidade de Novas Simulações: Há uma necessidade urgente de mais simulações de relatividade numérica que cubram um espaço de parâmetros mais amplo (incluindo spins e assimetrias extremas) e que forneçam previsões de ejetado mais detalhadas, especialmente para a evolução temporal e dependência angular.
• Tratamento de Incertezas: A análise demonstra que ignorar as incertezas sistemáticas nos modelos de ejetado (especialmente a velocidade) leva a resultados enviesados. Futuras análises devem tratar a massa e a velocidade do ejetado com parâmetros de incerteza explícitos e separados.
• Futuro da Astronomia Multimensageira: Para que a inferência conjunta seja confiável, é necessário desenvolver mapeamentos físicos mais realistas (não apenas ajustes empíricos) e melhorar a fidelidade das simulações de transferência radiativa. Sem isso, a interpretação de futuros eventos de fusão de estrelas de nêutrons continuará sujeita a grandes incertezas modelares.

Em resumo, o paper conclui que, embora os dados observacionais de GW170817 sejam ricos, nossa capacidade de extrair física fundamental deles está atualmente travada pela incerteza na física de como o material é ejetado durante a fusão, e não pela qualidade dos dados observacionais.