Improving low-latency multi-messenger follow-up of neutron star-black hole mergers with mode-by-mode filtering

Resumo Técnico: Melhorando o Acompanhamento de Multi-Mensageiros de Baixa Latência de Fusões de Estrela de Nêutrons–Buraco Negro com Filtragem Modo a Modo

Definição do Problema
A estimativa de parâmetros (PE) rápida para fusões de estrela de nêutrons–buraco negro (NSBH) é crítica para orientar as observações de acompanhamento eletromagnético (EM). Os atuais pipelines de baixa latência, como o BAYESTAR, tipicamente modelam apenas o harmônico quadrúpolo dominante de ondas gravitacionais (GW), $(\ell, |m|) = (2, 2)$. Para binários assimétricos como NSBHs (onde as razões de massa $q \gtrsim 5\text{--}10$), essa abordagem deixa degenerescências fortes entre a distância de luminosidade ($d_L$), inclinação ($\iota$) e parâmetros intrínsecos do binário. Especificamente, em formas de onda compostas apenas por quadrúpolos, $d_L$ correlaciona-se fortemente com $\iota$, e a razão de massa correlaciona-se com o spin efetivo. Essas degenerescências ampliam os volumes de localização, enfraquecem as restrições sobre a natureza do objeto secundário (NS vs. BH) e degradam as estimativas de inclinação, o que dificulta as decisões relativas à detectabilidade EM e associação com galáxias hospedeiras. Embora os modos de ordem superior (HMs) possam quebrar essas degenerescências, a amostragem estocástica completa incluindo HMs é computacionalmente proibitiva para aplicações de baixa latência, muitas vezes exigindo horas ou dias.

Metodologia
Os autores propõem um método de PE de baixa latência que realiza a filtragem modo a modo da série temporal da razão sinal-ruído (SNR) para os harmônicos $(2, 2)$, $(3, 3)$ e $(4, 4)$. A abordagem utiliza a biblioteca flywheel, uma adaptação de baixa latência do código cogwheel, para marginalizar analiticamente ou semi-analiticamente sobre parâmetros que entram simplesmente na verossimilhança (por exemplo, $d_L$ e fase de referência) enquanto utiliza amostragem de importância adaptativa para os parâmetros extrínsecos restantes.

O algoritmo recebe três entradas principais:

1. $\langle d | h_{\ell m} \rangle$: A série temporal de SNR para cada modelo de modo.
2. $\langle h_{\ell m} | h_{\ell' m'} \rangle$: A matriz de covariância de modo.
3. Amostras a priori para razões de SNR de modo ($R_{\ell m}$) derivadas de uma distribuição astrofísica assumida.

O método marginaliza sobre parâmetros extrínsecos do binário e um subconjunto de parâmetros intrínsecos que afetam a razão de amplitude dos modos. A integral da distância de luminosidade é avaliada via interpolação de uma tabela pré-computada, a integral da fase de referência via quadratura do trapézio, e as demais integrais (inclinação, localização no céu, polarização) via amostragem de importância adaptativa. Amostras da posterior são geradas por reamostragem de importância de amostras de Monte Carlo quase-aleatórias usadas na integral.

Resultados Principais
O método foi validado usando dois conjuntos de dados: detecções de NSBH simuladas em uma rede LIGO-Virgo com sensibilidade de design de O5 e dados públicos de eventos NSBH detectados anteriormente.

• Detecções Simuladas de O5: Aplicado a $\sim 1,6 \times 10^4$ eventos simulados com SNR $\ge 8$, a inclusão de HMs melhorou significativamente as restrições comparado a análises apenas de $(2, 2)$:

  • Distância e Volume: O método reduziu a largura do intervalo de credibilidade de 90% para a distância de luminosidade e o volume comóvel de localização por fatores que frequentemente excedem 2. Para um evento do tipo GW190814, o volume comóvel foi reduzido por um fator de $\sim 4,2$.
  • Inclinação e Ângulo de Visão: Os HMs quebraram a degenerescência $d_L$–$\iota$, melhorando as restrições do ângulo de visão em até um fator de 5 para sistemas face-on/off. Isso impacta diretamente a probabilidade de detecção de surtos gama de raios-X (GRBs) de jatos relativísticos.
  • Massa Secundária: As restrições sobre a massa secundária no referencial da fonte ($m_{src,2}$) melhoraram, particularmente para sistemas mais pesados. Para sistemas de alta massa, a probabilidade de o secundário ser uma estrela de nêutrons ($m_{src,2} \le 2,7 M_\odot$) aumentou em até 50%, auxiliando estratégias de acompanhamento de kilonovas.
  • Custo Computacional: A marginalização de modos superiores leva aproximadamente $0,13 \pm 0,03$ segundos por evento em um único núcleo de CPU, com um adicional de $0,11$ segundos para desenhar 1000 amostras da posterior. Isso é comparável aos atuais pipelines de apenas quadrúpolo.

• Validação com Dados Reais:

  • GW190814: Como um evento de alto SNR ($\sim 25$) e alta razão de massa ($q \sim 9$), o GW190814 mostrou uma assinatura clara de modo superior. Aplicar o método aos dados reais estreitou a posterior de $d_L$ por um fator de 1,6, a inclinação por um fator de 2 e o volume de localização por um fator de 2,1, trazendo a posterior de baixa latência para uma concordância mais próxima com os resultados de PE Bayesiano completo.
  • Outros Eventos NSBH: Para as outras três detecções confiáveis de NSBH (GW200105, GW200115, GW230529), que possuem SNRs mais baixos ($\sim 11\text{--}14$) e massas menores, nenhuma melhoria significativa foi observada, o que é consistente com seu conteúdo de modo subdominante mais fraco.

Significância e Alegações
O artigo alega que incluir modos de ordem superior em PE rápida produz medições de localização e de massa secundária substancialmente melhoradas com um custo computacional comparável aos métodos atuais de apenas $(2, 2)$. Os autores argumentam que essas melhorias são cruciais para estratégias de acompanhamento EM, permitindo:

1. Identificação mais eficiente de galáxias hospedeiras a partir de catálogos.
2. Melhor determinação de ângulos de visão para avaliar a detectabilidade de emissão de jato (GRBs/raios-X).
3. Melhor avaliação da detectabilidade de kilonova baseada no potencial de ruptura de maré.
4. Classificação de fonte mais confiante (NS vs. BH) para acionar protocolos de acompanhamento apropriados.

Os autores concluem que seu método está pronto para implantação em futuras corridas de observação (especificamente O5) para generalizar as atuais buscas de baixa latência de apenas quadrúpolo. Eles observam que, embora o método atualmente assuma binários não-precessantes, incorporar precessão é uma extensão natural para trabalhos futuros. O método foi projetado para ser agnóstico ao pipeline, exigindo apenas as séries temporais de SNR modo a modo e sua covariância, que podem ser geradas por qualquer pipeline de busca.