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

Resumen Técnico: Perturbaciones inducidas por gas en la aproximación de espiral de ondas gravitacionales de binarias de agujeros negros masivos de post-Newtonianos vivos para LISA

Planteamiento del Problema
Se espera que el Interferómetro Espacial de Láser LISA detecte binarias de agujeros negros masivos (MBHB) en la banda de mili-Hertz, particularmente aquellas cerca de la fusión. Una cuestión crítica abierta es si los discos circunbinarios (CBD) que rodean estos sistemas dejan huellas detectables en la forma de onda de las ondas gravitacionales (GW). Estudios previos han modelado los efectos del gas añadiendo linealmente torques gaseosos derivados semi-analíticamente a la evolución impulsada por GW. Este enfoque asume un desacoplamiento entre la dinámica del gas y la emisión de GW, lo que podría omitir acoplamientos no lineales. Además, mientras que se han estudiado los discos delgados ($H/R = 0.03$), el impacto de los discos más gruesos ($H/R = 0.1$), que son computacionalmente más eficientes de simular debido a tiempos viscosos más cortos, sobre la evolución acoplada GW-gas sigue siendo menos explorado.

Metodología
Los autores realizan simulaciones de hidrodinámica 3D (HD) utilizando el código GIZMO para modelar una MBHB de masa igual ($10^6 M_\odot$ de masa total) de órbita cuasi-circular embebida en un CBD progrado, localmente isotérmico, con una relación de aspecto $H/R = 0.1$.

• Evolución de la Binaria: La binaria es tratada como un sistema "vivo" que evoluciona bajo torques gaseosos y emisión de GW. La dinámica se describe mediante correcciones post-Newtonianas (PN) de orden 2.5, donde el término 2.5PN impulsa la aproximación de espiral inducida por GW.
• Configuraciones de Simulación: Se comparan cuatro configuraciones:
  1. Gas + 2.5PN: Evolución acoplada completa (GW + gas).
  2. Gas + 2PN: Torques gaseosos presentes, pero la emisión de GW es negada (término 2.5PN eliminado).
  3. 2.5PN: Evolución impulsada por GW en el vacío.
  4. 2PN: Evolución en el vacío sin GWs.
• Condiciones Iniciales: Las simulaciones comienzan con una separación de la binaria de $54.5 r_s$ (radios de Schwarzschild), lo que corresponde aproximadamente a un año antes de la fusión a un desplazamiento al rojo $z \sim 1$. El sistema evoluciona durante 600 periodos orbitales ($P_B$), reduciéndose a $46.5 r_s$.
• Resolución: Un estudio de resolución confirma que la ejecución de resolución media ($\Delta x \approx 2.1 r_s$ en $3a$) es suficientemente convergente para las mediciones de torque y desfase.

Contribuciones Clave

1. Primera Medición para $H/R=0.1$: Este trabajo reporta la primera medición de torques gravitacionales y de acreción para una MBHB viva en un CBD con $H/R = 0.1$, comparando casos con y sin emisión concurrente de GW.
2. Dinámica Acoplada GW-Gas: A diferencia de los enfoques semi-analíticos previos, este estudio simula directamente el acoplamiento no lineal entre la aproximación de espiral impulsada por GW y la dinámica del gas, permitiendo una medición precisa del desfase orbital resultante.
3. Morfología de la Acreción: El artículo analiza cómo cambia la morfología de la serie temporal de acreción cuando la emisión de GW se convierte en el mecanismo evolutivo dominante.

Resultados

• Mediciones de Torque: Los parámetros de torque adimensionales ($\xi$) se midieron durante 600 órbitas.
  • Torque de Acreción ($\xi_{acc}$): $\sim 0.44$ (GW+gas) frente a $\sim 0.43$ (solo gas).
  • Torque Gravitacional ($\xi_{grav}$): $\sim 0.13$ (GW+gas) frente a $\sim 0.16$ (solo gas).
  • Torque Total: El torque neto es positivo en ambos casos, lo que indica que el gas ralentiza la aproximación de espiral. Los autores señalan que para $H/R=0.1$, el torque de acreción es aproximadamente el doble de la magnitud del torque gravitacional, contrastando con los resultados de discos más delgados donde el torque gravitacional domina.
  • Desacoplamiento: El CBD sigue la aproximación de espiral de la binaria, pero el torque gravitacional del disco exterior ($R > a$) disminuye ligeramente en magnitud a medida que la binaria se reduce, lo cual es consistente con el inicio del desacoplamiento.
• Modulación de la Tasa de Acreción: La presencia de GWs altera modestamente la morfología de la tasa de acreción. Los espectrogramas muestran una disminución de potencia para la modulación del borde interno de la cavidad ($\sim 3\%$) y la modulación del periodo orbital de la binaria (aumentando del $15\%$ a $\sim 60\%$ con el tiempo) a medida que el disco se desacopla progresivamente de la binaria.
• Desfases:
  • El desfase orbital inducido por gas con emisión de GW ($\delta\phi_{orb}^{(GW)}$) se mide en 0.12 rad sobre 600 ciclos.
  • El desfase sin emisión de GW ($\delta\phi_{orb}^{(NoGW)}$) es de 0.08 rad.
  • El desfase de GW resultante ($\delta\phi_{GW} = 2\delta\phi_{orb}^{(GW)}$) es de 0.24 rad.
  • Comparación Analítica: Un cálculo semi-analítico estándar (sumando linealmente las tasas de GW y gas) produce un desfase de $\sim 0.002$ rad. El resultado numérico es aproximadamente 60 veces mayor que esta expectativa analítica, resaltando la importancia del acoplamiento no lineal GW-gas.
• Detectabilidad: Para una fuente en $z=1$ con una relación señal-ruido (SNR) de $\sim 1300$, el desfase calculado de 0.24 rad excede el umbral de detección ($8/\text{SNR} \approx 0.006$ rad), lo que sugiere que LISA podría detectar esta huella ambiental.

Significado y Reivindicaciones
Los autores afirman que sus resultados tienen implicaciones directas para la astronomía de mensajeros múltiples. La detección de la modulación de la tasa de acreción por parte de estudios ópticos (ej. LSST/Roman) combinada con la detección del desfase por parte de LISA proporcionaría información crucial sobre el complejo entorno que rodea a las MBHB.

El artículo enfatiza que la gran discrepancia entre el desfase numérico y la predicción semi-analítica sugiere que simplemente sumar linealmente los torques de gas a la evolución de GW es insuficiente para un modelado preciso. El estudio sirve como un referente que demuestra la necesidad de incluir la emisión concurrente de GW en las simulaciones de hidrodinámica para predecir correctamente los observables de LISA. Los autores señalan que, aunque su enfoque asume una condición inicial de estado estacionario y viscosidad de Shakura-Sunyaev (ignorando la turbulencia magnetohidrodinámica), sigue siendo el método más factible para inferir predicciones para fuentes de LISA en separaciones donde son observables cientos de órbitas.