On constraining the initial orbital eccentricity of inspiral-dominated gravitational-wave events with TaylorF2Ecck

Resumo Técnico: Restringindo a Excentricidade Orbital Inicial em Eventos de Inspiração Dominante com o TaylorF2Ecck

Problema
O consórcio LIGO-Virgo-KAGRA (LVK) detectou eventos de ondas gravitacionais (OG) de inspiração dominante, como as fusões de estrelas de nêutrons binárias (BNS) GW170817 e GW190425. Esses eventos passam significativamente mais tempo na janela de frequência do detector do que as fusões de buracos negros binários (BBH), tornando-os ideais para restringir a excentricidade orbital residual ($e_0$). Modelos de forma de onda padrão como o TaylorF2 assumem órbitas quase circulares, embora canais de formação astrofísica (ex: captura dinâmica em aglomerados densos) prevejam uma excentricidade não desprezível. Análises anteriores usando o aproximante TaylorF2Ecc restringiram $e_0$, mas omitiram os efeitos de avanço do periastro, que entram na dinâmica de binárias compactas na ordem 1PN. Além disso, modelos excêntricos existentes frequentemente carecem de consistência na ordem de precisão pós-Newtoniana (PN) entre as contribuições de fase e amplitude, o que pode introduzir vieses sistemáticos na estimativa de parâmetros.

Metodologia
Os autores introduzem o TaylorF2Ecck, um novo aproximante de inspiração no domínio da frequência implementado dentro da estrutura LALSuite. Este modelo descreve analiticamente as OG de binárias compactas não rotacionantes em órbitas excêntricas. As principais características técnicas incluem:

• Precisão PN: O modelo incorpora consistentemente os efeitos orbitais, de avanço do periastro e de emissão de OG totalmente até a ordem 3PN na fase de Fourier.
• Ordem de Excentricidade: Inclui contribuições de $O(e_0^2)$ para a fase com precisão de 3PN e correções de $O(e_0)$ para as amplitudes de ordem quadripolar.
• Harmônicos: Diferente do TaylorF2Ecc, que foca apenas no modo dominante, o TaylorF2Ecck inclui múltiplos harmônicos $(j, n)$, especificamente $(1,0), (1,-2), (2,-2)$ e $(3,-2)$, para contabilizar os sidebands e o desdobramento de modos.
• Modelos Comparativos: O estudo emprega uma versão restrita do TaylorF2Ecck (denotada TaylorF2Ecch) que omite os efeitos de avanço do periastro para isolar seu impacto. Também utiliza variantes do TaylorF2 e TaylorF2Ecc quase circulares com fases de Fourier truncadas nas ordens 3PN, 3.5PN, 4PN e 4.5PN para investigar vieses sistemáticos.
• Estrutura de Análise: A estimativa bayesiana de parâmetros foi realizada utilizando o software Bilby com amostragem aninhada (nested sampling). A análise focou nos eventos GW170817 e GW190425, utilizando dados oficiais da LVK e densidades espectrais de potência do BayesWave. A frequência mínima foi definida em 20 Hz para equilibrar o custo computacional e a sensibilidade à excentricidade.

Principais Contribuições

1. Implementação do TaylorF2Ecck: O artigo fornece uma família de modelos analiticamente consistente e computacionalmente eficiente que inclui as contribuições de excentricidade de ordem líder até a ordem 3PN, incluindo o avanço do periastro e efeitos de múltiplos harmônicos.
2. Verificações de Sanidade e Distinguibilidade: Através de cálculos de match e testes de injeção-recuperação, os autores estabeleceram que o TaylorF2Ecck torna-se distinguível do modelo TaylorF2 quase circular em $e_0 \gtrsim 0.01$ e do modelo que carece de avanço do periastro (TaylorF2Ecch) em $e_0 \gtrsim 0.06$ (para parâmetros do tipo GW170817).
3. Sensibilidade da Ordem PN: O estudo identificou um viés sistemático significativo na distribuição posterior da razão de massas ($q$) ao comparar modelos excêntricos com precisão de 3PN contra modelos de fase de Fourier circular de ordem 3.5PN. Isso destaca a necessidade de estender as contribuições excêntricas para pelo menos a ordem 3.5PN para uma estimativa de parâmetros confiável.

Resultados

• Restrições de Excentricidade: Para ambos GW170817 e GW190425, as distribuições posteriores para a excentricidade inicial $e_0$ (avaliada em 20 Hz) tendem para zero.
  • GW170817: $e_0 \leq 0.016$ (intervalo de credibilidade de 90%) com um prior uniforme.
  • GW190425: $e_0 \leq 0.023$ (intervalo de credibilidade de 90%) com um prior uniforme.
• Preferência de Modelo: Comparações de fator de Bayes ($|\Delta \log_{10} B| < 1$) não mostram evidência forte favorecendo o modelo excêntrico TaylorF2Ecck sobre o modelo quase circular TaylorF2, nem sobre o modelo TaylorF2Ecc (que carece de avanço do periastro). Isso indica que, nas excentricidades ínfimas inferidas, a inclusão do avanço do periastro e de harmônicos extras não altera significativamente os resultados.
• Vieses Sistemáticos: Embora as restrições de excentricidade sejam robustas, o estudo encontrou um deslocamento substancial na posterior da razão de massas ao comparar modelos de 3PN com modelos circulares de 3.5PN. Isso sugere que a conhecida correlação entre massa e excentricidade torna a inclusão de correções de fase circular de ordem superior (3.5PN) crítica, mesmo quando a excentricidade é desprezível.
• Comparação com GW200105: Uma análise suplementar do evento NSBH GW200105 rendeu estimativas de excentricidade consistentes com estudos anteriores ($e_0 \approx 0.12$), embora os fatores de Bayes não tenham favorecido fortemente o modelo excêntrico sobre o quase circular usando a implementação atual do TaylorF2Ecck.

Significância e Alegações
O artigo alega que o TaylorF2Ecck fornece um arcabouço robusto e analiticamente consistente para analisar eventos de inspiração de baixa excentricidade. A principal significância reside na restrição rigorosa da excentricidade inicial para GW170817 e GW190425, confirmando que são consistentes com binárias formadas em campo (excentricidade desprezível) em vez de canais de formação dinâmica.

Crucialmente, os autores argumentam que, embora as restrições atuais de $e_0$ sejam robustas, o viés observado nas posteriors da razão de massas ressalta um imperativo metodológico: modelos de forma de onda excêntricos devem incorporar contribuições de excentricidade inicial até pelo menos a ordem 3.5PN para evitar erros sistemáticos na estimão de parâmetros. O artigo conclui que trabalhos futuros devem focar na extensão da precisão PN dos termos excêntricos para 3.5PN e na incorporação de efeitos de spin, fusão (merger) e ringdown para explorar plenamente o potencial das buscas por OG excêntricas.