Gas-induced perturbations on the gravitational wave in-spiral of live post-Newtonian LISA massive black hole binaries: 0.1 disk aspect ratio

Resumo Técnico: Perturbações induzidas por gás na inspiração de ondas gravitacionais de binárias de buracos negros massivos pós-newtonianas vivas para o LISA

Problema
Espera-se que o Laser Interferometer Space Antenna (LISA) detecte binárias de buracos negros massivos (MBHBs) na banda de mili-Hertz, particularmente aquelas próximas da fusão. Uma questão crítica em aberto é se os discos circumbinários (CBDs) que circundam esses sistemas deixam marcas detectáveis na forma de onda de ondas gravitacionais (GW). Estudos anteriores frequentemente modelaram os efeitos gasosos adicionando linearmente torques gasosos derivados semi-analiticamente à evolução impulsionada por GW. Essa abordagem assume um desacoplamento entre a dinâmica do gás e a emissão de GW, potencialmente perdendo acoplamentos não lineares. Além disso, embora discos finos ($H/R = 0,03$) tenham sido estudados, o impacto de discos mais espessos ($H/R = 0,1$), que são computacionalmente mais eficientes de simular devido aos tempos viscosos mais curtos, sobre a evolução acoplada GW-gás permanece menos explorado.

Metodologia
Os autores realizam simulações de hidrodinâmica 3D (HD) usando o código GIZMO para modelar uma MBHB de massa igual ($10^6 M_\odot$ de massa total) quase circular inserida em um CBD progrado e localmente isotérmico com uma razão de aspecto $H/R = 0,1$.

• Evolução da Binária: A binária é tratada como um sistema "vivo" evoluindo sob torques gasosos e emissão de GW. A dinâmica é descrita usando correções pós-newtonianas (PN) de 2,5 (2,5PN), onde o termo 2,5PN impulsiona a inspiração induzida por GW.
• Configurações de Simulação: Quatro configurações são comparadas:
  1. Gás + 2,5PN: Evolução completa acoplada (GW + gás).
  2. Gás + 2PN: Torques gasosos presentes, mas a emissão de GW negligenciada (termo 2,5PN removido).
  3. 2,5PN: Evolução impulsionada por GW no vácuo.
  4. 2PN: Evolução no vácuo sem GWs.
• Condições Iniciais: As simulações começam com uma separação binária de $54,5 r_s$ (raios de Schwarzschild), correspondendo a aproximadamente um ano antes da fusão em um redshift $z \sim 1$. O sistema evolui por 600 períodos orbitais ($P_B$), encolhendo para $46,5 r_s$.
• Resolução: Um estudo de resolução confirma que a execução de resolução média ($\Delta x \approx 2,1 r_s$ em $3a$) é suficientemente convergente para medições de torque e deslocamento de fase.

Principais Contribuições

1. Primeira Medição para $H/R=0,1$: Este trabalho relata a primeira medição de torques gravitacionais e de acreção para uma MBHB viva em um CBD com $H/R = 0,1$, comparando casos com e sem emissão de GW concomitante.
2. Dinâmica Acoplada GW-Gás: Ao contrário das abordagens semi-analíticas anteriores, este estudo simula diretamente o acoplamento não linear entre a inspiração impulsionada por GW e a dinâmica do gás, permitindo uma medição precisa do deslocamento de fase resultante.
3. Morfologia de Acreção: O artigo analisa como a morfologia da série temporal de acreção muda quando a emissão de GW se torna o mecanismo evolutivo dominante.

Resultados

• Medições de Torque: Os parâmetros de torque adimensional ($\xi$) foram medidos ao longo de 600 órbitas.
  • Torque de Acreção ($\xi_{acc}$): $\sim 0,44$ (GW+gás) vs. $\sim 0,43$ (apenas gás).
  • Torque Gravitacional ($\xi_{grav}$): $\sim 0,13$ (GW+gás) vs. $\sim 0,16$ (apenas gás).
  • Torque Total: O torque líquido é positivo em ambos os casos, indicando que o gás retarda a inspiração. Os autores observam que, para $H/R=0,1$, o torque de acreção é aproximadamente duas vezes a magnitude do torque gravitacional, contrastando com resultados de discos mais finos onde o torque gravitacional domina.
  • Desacoplamento: O CBD segue a inspiração da binária, mas o torque gravitacional do disco externo ($R > a$) diminui ligeiramente em magnitude conforme a binária encolhe, consistente com o início do desacoplamento.
• Modulação da Taxa de Acreção: A presença de GWs altera modestamente a morfologia da taxa de acreção. Espectrogramas mostram uma diminuição de potência para a modulação da borda interna da cavidade ($\sim 3\%$) e a modulação do período orbital da binária (aumentando de $15\%$ para $\sim 60\%$ ao longo do tempo) conforme o disco se desacopla progressivamente da binária.
• Deslocamentos de Fase:
  • O deslocamento de fase orbital induzido pelo gás com emissão de GW ($\delta\phi_{orb}^{(GW)}$) é medido em 0,12 rad ao longo de 600 ciclos.
  • O deslocamento de fase sem emissão de GW ($\delta\phi_{orb}^{(NoGW)}$) é 0,08 rad.
  • O deslocamento de fase de GW resultante ($\delta\phi_{GW} = 2\delta\phi_{orb}^{(GW)}$) é 0,24 rad.
  • Comparação Analítica: Um cálculo semi-analítico padrão (somando linearmente as taxas de GW e de gás) produz um deslocamento de fase de $\sim 0,002$ rad. O resultado numérico é aproximadamente 60 vezes maior do que essa expectativa analítica, destacando a significância do acoplamento não linear GW-gás.
• Detectabilidade: Para uma fonte em $z=1$ com uma razão sinal-ruído (SNR) de $\sim 1300$, o deslocamento de fase calculado de 0,24 rad excede o limiar de detecção ($8/\text{SNR} \approx 0,006$ rad), sugerindo que o LISA poderia detectar essa impressão ambiental.

Significância e Alegações
Os autores alegam que seus resultados têm implicações diretas para a astronomia multi-mensageiro. A detecção da modulação da taxa de acreção por levantamentos ópticos (ex: LSST/Roman) combinada com a detecção do deslocamento de fase pelo LISA forneceria informações cruciais sobre o ambiente complexo que envolve as MBHBs.

O artigo enfatiza que a grande discrepância entre o deslocamento de fase numérico e a previsão semi-analítica sugere que simplesmente somar linearmente os torques gasosos à evolução de GW é insuficiente para uma modelagem precisa. O estudo serve como um benchmark demonstrando a necessidade de incluir a emissão de GW concomitante em simulações de hidrodinâmica para prever corretamente os observáveis do LISA. Os autores observam que, embora sua abordagem assuma uma condição inicial de estado estacionário e viscosidade de Shakura-Sunyaev (ignorando a turbulência magnetohidrodinâmica), continua sendo o método mais viável para inferir previsões para fontes do LISA em separações onde centenas de órbitas são observáveis.