Dynamics and detectability of long-lived non-accretion phases for massive black hole binaries in cold, thermally regulating disks

Resumen Técnico: Dinámica y Detectabilidad de las Fases de No Acreción de Larga Duración para Binarias de Agujeros Negros Masivos

Planteamiento del Problema
Las simulaciones hidrodinámicas previas de binarias de agujeros negros supermasivos (SMBHB) interactuando con discos circumbinarios (CBD) han dependido en gran medida de supuestos localmente isotérmicos o de discos con razones de aspecto relativamente altas ($h/r \sim 0.1$). Estos estudios sugirieron que, si bien el potencial de la binaria trunca el disco interno, el gas sigue acreciendo mediante corrientes de marea. Sin embargo, se espera que los discos de las SMBHB realistas sean significativamente más delgados ($h/r \sim 10^{-2} - 10^{-3}$), donde la reducción del soporte de presión puede, funcionalmente, cesar la acreción. Una cuestión crítica abierta es si las fases de "no acreción desbocada" observadas en simulaciones de discos isotérmicos y delgados sobreviven cuando la termodinámica del disco se evoluciona de manera autoconsistente. Específicamente, no está claro si el calentamiento viscoso e hidrodinámico en el disco interno puede calentar el gas lo suficiente como para restaurar el soporte de presión y permitir que la acreción se reanude, o si el sistema permanece en un estado de acreción suprimida de larga duración.

Metodología
Los autores presentan simulaciones hidrodinámicas de alta resolución basadas en rejilla utilizando el código de volumen finito Sailfish. El estudio resuelve las ecuaciones de Navier-Stokes compresibles en 2D, promediadas verticalmente, para un disco de gas alrededor de una binaria de masas iguales. Las características metodológicas clave incluyen:

• Ecuación de Energía: A diferencia de los estudios isotérmicos previos, los autores evolucionan una ecuación de energía que rastrea de forma autoconsistente el calentamiento viscoso e hidrodinámico acoplado con el enfriamiento de cuerpo negro radiativo. Se asume que el disco es ópticamente grueso con una opacidad dominada por la dispersión de electrones.
• Configuración del Sistema: Las simulaciones modelan una binaria con una masa total $M = 10^6 M_\odot$ y una separación $a = 0.001$ pc (periodo orbital $\sim 3$ años) dentro de un disco $\alpha$ dominado por la presión de gas ($\alpha=0.1$).
• Espacio de Parámetros: Los autores exploran tres números de Mach ($M_a = \{12, 50, 80\}$), que corresponden a fracciones de Eddington efectivas que van desde regímenes super-Eddington ($\sim 10^5$) hasta casi-Eddington ($\sim 10$).
• Prescripción de Sumidero (Sink): Los autores implementan una prescripción de sumidero "adiabática libre de torque" modificada. Esta actualización corrige las pérdidas espurias de energía interna asociadas con la eliminación de masa libre de torque en simulaciones previas de evolución de energía, asegurando soluciones de disco de un solo agujero negro estables fuera del radio del sumidero.
• Pruebas de Convergencia: Las simulaciones se ejecutaron en múltiples resoluciones ($\Delta x = 0.015a, 0.0075a, 0.005a$) para verificar la convergencia del parámetro de torque y las tasas de acreción.

Resultados Clave

1. Persistencia de las Fases de No Acreción: El estudio confirma que las fases de no acreción de larga duración persisten incluso con una ecuación de energía autoconsistente. Mientras que el gas se acumula y se calienta en el borde lejano (apocentro) de la cavidad circumbinaria, el material que lanza las corrientes de acreción cerca del pericentro de la cavidad permanece comparativamente frío. En consecuencia, estas corrientes carecen del soporte de presión interna suficiente para resistir la compactación de marea por el potencial de la binaria, lo que conduce a una supresión sostenida de la tasa de acreción de la binaria ($\dot{M}_b \ll \dot{M}_\infty$).
2. Moderación de la No Acreción Desbocada: La inclusión de la ecuación de energía debilita, pero no elimina, el problema de la "no acreción desbocada" encontrado en las soluciones isotérmicas. El parámetro de torque ($\ell_0$), que gobierna la supresión de la acreción, todavía crece exponencialmente con el número de Mach, aunque a un ritmo más lento ($\ln \ell_0 \propto 0.085 M_a$) en comparación con las soluciones isotérmicas ($\ln \ell_0 \propto 0.21 M_a$).
3. Estructura Térmica y Formación de Corrientes: La cavidad exhibe una asimetría térmica distintiva. El apocentro es calentado por acumulaciones de densidad y la generación de corrientes, mientras que el pericentro se enfría eficientemente mediante expansión adiabática a medida que los elementos de fluido regresan en órbitas excéntricas. Este mecanismo de enfriamiento es dominante sobre las pérdidas radiativas en la región de lanzamiento de corrientes.
4. Propiedades Electromagnéticas: A pesar de las bajas tasas de acreción, estos sistemas no son electromagnéticamente oscuros. El borde interno calentado del CBD permanece luminoso en las frecuencias de óptico e infrarrojo cercano (NIR). Sin embargo, los sistemas son intrínsecamente débiles en rayos X debido a la truncación del disco interno y la falta de un continuo de alta energía y fotoionizante.
5. Escalas de Tiempo de Fusión vs. Activación: Para números de Mach altos ($M_a \gtrsim 90$), el tiempo requerido para que la binaria se "active" (alcance un estado de acreción estacionario) puede exceder la vida de fusión de la binaria impulsada por ondas gravitacionales. Esto sugiere que las binarias masivas en discos dominados por la presión de gas pueden fusionarse mientras aún se encuentran en una fase de no acreción.

Significancia y Reivindicaciones
El artículo afirma que el problema de la "no acreción desbocada" es una característica robusta de los discos delgados dominados por la presión de gas, sobreviviendo incluso cuando se incluye la retroalimentación térmica. Esto tiene profundas implicaciones para la identificación electromagnética de las SMBHB:

• Detectabilidad: Las binarias en fase de no acreción pueden ser detectables en los próximos sondeos de campo amplio (por ejemplo, LSST, Telescopio Roman) como fuentes periódicamente variables en las bandas óptica e infrarrojo cercano (NIR), a pesar de su baja luminosidad bolométrica.
• Firmas Espectrales: Estos sistemas pueden aparecer como núcleos galácticos activos (AGN) intrínsecamente débiles en rayos X, con líneas de emisión débiles o ausentes, lo que potencialmente los vincula con la población de cuásares de líneas débiles.
• Fenómenos de "Cambio de Aspecto" (Changing-Look): La alimentación intermitente de minidiscos, impulsada por la compleja dinámica de la cavidad truncada, ofrece un mecanismo plausible para el comportamiento de los AGN de "cambio de aspecto", donde los flujos ionizantes varían significamente con el tiempo.

Los autores concluyen que, si bien sus resultados se basan en modelos de presión de gas, promediados verticalmente y no magnetizados con enfriamiento simplificado, proporcionan una línea base crítica para comprender cómo la autoconsistencia termodinámica afecta la dinámica de acreción de las binarias. Sugieren que el trabajo futuro debe incorporar campos magnéticos, presión de radiación y geometría 3D para caracterizar completamente la transición al estado de acreción estacionario y la fase final de fusión.