Time-domain phenomenological multipolar waveforms for aligned-spin binary black holes in elliptical orbits

Resumo Técnico: Formas de onda multipolares fenomenológicas no domínio do tempo para buracos negros binários com spin alinhado em órbitas elípticas

Definição do Problema
Detectores de ondas gravitacionais (GW) baseados no solo identificaram inúmeras coalescências de binários compactos, a maioria das quais é consistente com órbitas quase circulares (QC) formadas via evolução binária isolada. No entanto, binários formados através de interações dinâmicas em ambientes estelares densos ou via o mecanismo de Kozai-Lidov podem reter uma excentricidade orbital significativa quando entram na banda do detector. Detectar e caracterizar essa excentricidade é crucial para compreender as origens astrofísicas desses sistemas. Embora existam modelos de forma de onda de última geração para binários de spin alinhado e QC, modelos fenomenológicos precisos no domínio do tempo para sistemas excêntricos com modos de ordem superior e efeitos de spin têm sido limitados. Os modelos excêntricos existentes frequentemente dependem de formalismos de One-Body Efetivo (EOB) ou de surrogados de relatividade numérica (NR), que podem ser computacionalmente caros ou limitados em sua cobertura de espaço de parâmetros (por exemplo, restritos a binários sem spin ou de massas comparáveis). Além disso, simulações de NR públicas para binários excêntricos são escassas, complicando o desenvolvimento e a validação de novos modelos.

Metodologia
Os autores introduzem o IMRPhenomTEHM, um novo modelo fenomenológico no domínio do tempo para buracos negros binários (BBHs) excêntricos e de spin alinhado. O modelo é construído estendendo o modelo preciso de spin alinhado e QC, IMRPhenomTHM, com correções pós-newtonianas (PN) excêntricas.

• Construção do Modelo: A estratégia central envolve a adição de correções PN excêntricas à amplitude e à fase dos modos do IMRPhenomTHM. O modelo inclui o modo dominante $(2, \pm 2)$ de harmônicos esféricos e os modos subdominantes $(2, \pm 1)$, $(3, \pm 3)$, $(4, \pm 4)$ e $(5, \pm 5)$.
• Dinâmica Excêntrica: O modelo incorpora a dinâmica completa de órbita média 3PN, incluindo correções de não-spin e de spin (órbita-spin e spin-spin). A evolução orbital é impulsionada por equações de evolução secular derivadas de equações de balanço que conectam fluxos de energia e momento angular.
• Sistemas de Coordenadas: Os autores implementam as equações de evolução secular tanto em coordenadas PN de harmônicos modificados (MH) quanto em coordenadas de One-Body Efetivo (EOB). Eles descobrem que as coordenadas EOB produzem uma melhor concordância com as simulações de NR e uma correspondência mais próxima com as definições de excentricidade baseadas em forma de onda ($e_{GW}$), adotando EOB como padrão.
• Suposição de Circularização: O modelo assume que o binário circulariza-se até a fusão. Consequentemente, as correções excêntricas desaparecem durante as fases de fusão e ringdown, garantindo uma transição suave de volta ao modelo IMRPhenomTHM subjacente.
• Redefinição de Fase: Para lidar com dependências de gauge, os autores exploram redefinições de fase. No entanto, para manter a consistência com a base QC calibrada por NR, o modelo padrão não aplica um deslocamento de fase específico, assumindo que a calibração de NR do modelo subjacente já contabiliza os ajustes necessários.
• Eficiência Computacional: O modelo utiliza uma abordagem altamente otimizada para resolver as equações diferenciais acopladas da evolução orbital. Ao interpolar a evolução da frequência QC em uma grade específica e separar os componentes de órbita média dos componentes oscilatórios, o modelo alcança acelerações significativas em comparação com os modelos excêntricos baseados em EOB.

Principais Contribuições

1. Primeiro Modelo Excêntrico Fenomenológico no Domínio do Tempo: O IMRPhenomTEHM é apresentado como o primeiro modelo de forma de onda fenomenológica no domínio do tempo para BBHs excêntricos com spin alinhado que inclui multipolos de ordem superior.
2. Inclusão de Multipolos e Spin: O modelo inclui interações de órbita-spin e spin-spin até 3PN (embora utilizando 2PN por padrão para consistência com a base QC) e cobre modos até $l=5$.
3. Validação contra NR: O modelo é validado contra 28 simulações de NR excêntricas públicas do catálogo Simulating eXtreme Spacetimes (SXS). Ele atinge valores de infidelidade (unfaithfulness/mismatch) abaixo de 2% para a maioria dos casos, mesmo sem calibração direta para conjuntos de dados de NR excêntricos.
4. Robustez e Espaço de Parâmetros: O modelo mostra-se robusto em um amplo espaço de parâmetros (razões de massa $q \in [1, 20]$, spins $\chi_i \in [-0.995, 0.995]$ e excentricidades $e < 0.4$ a 10 Hz).
5. Desempenho Computacional: Demonstra-se que o IMRPhenomTEHM é aproximadamente uma ordem de magnitude mais rápido que os modelos excêntricos de estado da arte baseados em EOB (SEOBNRv5EHM e TEOBResumDali) para massas baixas a intermediárias, tornando-o altamente adequado para estudos de estimação de parâmetros (PE).

Resultos

• Limite Quasi-Circular: O modelo recupera fielmente o limite QC, com mismatches contra o IMRPhenomTHM permanecendo abaixo de $10^{-5}$ em todo o espaço de parâmetros, garantindo que não haja vieses em estudos de PE devido à extensão excêntrica.
• Comparação com NR: Quando comparado a 28 simulações de NR do SXS, o IMRPhenomTEHM atinge mismatches geralmente abaixo de 2%. Os mismatches aumentam ligeiramente para excentricidades iniciais mais altas ($e_0 \sim 0.3$), atribuídos às limitações das expressões PN expandidas em excentricidade para excentricidades elevadas.
• Estimação de Parâmetros (PE) em Injeções de NR: Em estudos de injeção de ruído zero usando três simulações do SXS, o modelo recupera com sucesso os parâmetros injetados (massas, spins, excentricidade, anomalia média) dentro de intervalos de credibilidade de 90%. A inclusão de modos de ordem superior melhora a recuperação da distância de luminosidade e inclinação.
• Análise de Eventos Reais: O modelo foi aplicado para reanalisar GW150914 e GW190521.
  • GW150914: Os resultados são consistentes com a literatura, apoiando um binário não excêntrico, sem spin e de massa comparável. O fator de Bayes favorece ligeiramente a hipótese QC.
  • GW190521: O modelo produz resultados consistentes com estudos anteriores usando SEOBNRv5EHM. Embora os valores medianos de excentricidade sejam não-zero, os posteriors são amplamente não informativos, consistentes com a falta de evidência de excentricidade neste evento de alta massa. O estudo destaca que o IMRPhenomTEHM evita cortes artificiais de prior em excentricidade que podem ocorrer em outros modelos devido a restrições de geração de forma de onda em altas massas.
• Benchmarks de Tempo: Para um binário com $M=10 M_\odot$ e $e=0.2$, o IMRPhenomTEHM gera formas de onda em $\sim 1$ segundo (incluindo modos superiores), comparado a $\sim 8$ segundos para o SEOBNRv5EHM e $\sim 10^4$ segundos para o TEOBResumDali.

Significância e Alegações
O artigo afirma que o IMRPhenomTEHM fornece uma descrição confiável e computacionalmente eficiente de binários de buracos negros excêntricos com spins alinhados e excentricidades inferiores a $e=0.4$ a 10 Hz. Sua principal significância reside na capacidade de facilitar estudos de inferência Bayesiana em larga escala para as próximas corridas de observação de ondas gravitacionais, onde se espera aumentar a detecção de binários excêntricos. Ao oferecer uma alternativa rápida, multipolar e precisa aos modelos EOB mais lentos, ele permite a exploração sistemática da excentricidade nos dados de GW sem as proibições computacionais dos modelos atuais de última geração. Os autores observam que, embora a precisão atual seja suficiente para os detectores existentes, trabalhos futuros focarão em melhorar o tratamento das correções excêntricas nas amplitudes dos modos (indo além das expansões $O(e^6)$) e potencialmente incorporando precessão de spin.