Stochastic Variability of Binary Accretion

Resumen Técnico: Variabilidad Estocástica de la Acreción Binaria

Planteamiento del Problema
Los binarios de agujeros negros supermasivos (SMBH, por sus siglas en inglés) son una consecuencia predicha de la formación jerárquica de galaxias, aunque siguen siendo observacionalmente elusivos. Si bien se ha perseguido la modulación de flujo cuasi-periódica cerca de la frecuencia orbital del binario como un signo distintivo primario, tales búsquedas se ven complicadas por la naturaleza intrínsecamente estocástica de las curvas de luz de los Núcleos Galácticos Activos (AGN) y limitadas por líneas de base de observación que excluyen a la gran mayoría de los binarios sub-parsec con periodos orbitales largos. Además, las tasas de recuperación de periodogramas estándar disminuyen significamente para las formas de pulso no sinusoidales predichas por simulaciones hidrodinámicas. En consecuencia, existe la necesidad de identificar signos temporales en la densidad espectral de potencia (PSD) de las tasas de acreción que sean distintos de la variabilidad canónica de ruido rojo de un AGN simple, centrándose específicamente en la variabilidad estocástica en escalas de tiempo más cortas que el periodo orbital del binario.

Metodología
Los autores emplean simulaciones de hidrodinámica 2D de muy alta resolución utilizando el código Sailfish para modelar la acreción de gas en un binario circular de masas desiguales ($q = M_2/M_1 = 0.2$) rodeado por un disco de gas circumbinario. Las simulaciones utilizan una ecuación de estado localmente isotérmica con un número de Mach orbital $M=10$ y una prescripción de Shakura-Sunyaev $\alpha$ para la viscosidad ($\alpha = 0.01$).

Un aspecto crítico de la metodología es la variación sistemática del parámetro de tamaño del "sumidero" (sink), $s$, que determina los radios de truncamiento interno de los minidiscos componentes ($r_{sink} \propto s \cdot a$). El parámetro $s$ se varía de 0.05 a 0.40. Los autores distinguen entre el flujo de masa que entra en el minidisco (el "impulsor" o driver) y la tasa de acreción de masa sobre el sumidero del agujero negro (el "filtro"). Las simulaciones se ejecutan a través de una fase de calentamiento (burn-in) de 2900 órbitas seguida de una fase de muestreo de 100 órbitas a alta resolución ($\Delta \ln r = 0.0025$), con series temporales muestreadas $\sim 10^3$ veces por órbita para resolver frecuencias muy por encima de la frecuencia orbital.

Resultos Clave

1. Forma de la PSD del Continuo: En el límite de tamaños de sumidero pequeños ($s \lesssim 0.1$), la densidad espectral de potencia de la tasa de acreción converge a una forma universal de ley de potencia rota. La PSD es plana (ruido blanco) en bajas frecuencias, transiciona a una pendiente de $-2$ en una banda estrecha y se vuelve más empinada hacia una pendiente de $-4$ en altas frecuencias.
2. Frecuencia de Ruptura: La transición a la pendiente pronunciada de $-4$ ocurre en una frecuencia de ruptura $f_d \approx 5$ veces la frecuencia orbital del binario ($f_{bin}$). Esta frecuencia corresponde a la frecuencia Kepleriana local en el borde exterior del minidisco del secundario ($r_{out} \approx 0.18\text{--}0.20 a$).
3. Dependencia del Tamaño del Sumidero: Para sumideros grandes ($s \gtrsim 0.2$), la PSD no converge a la forma universal. En su lugar, el minidisco actúa como un amortiguador insuficiente, permitiendo que el exceso de potencia de alta frecuencia alcance el sumidero, lo que se manifiesta como picos en la tasa de acreción. Los autores argumentan que estos picos son artefactos numéricos cuando el radio del sumidero excede la órbita circular más interna estable (ISCO) física del agujero negro.
4. Mecanismo Impulsor-Filtro: La variabilidad se explica mediante un proceso estocástico de dos etapas:
   • Impulsor (Driver): Los flujos circumbinarios entregan gas al borde exterior del minidisco de una manera que se asemeja a un proceso de paseo aleatorio amortiguado (proceso de Ornstein-Uhlenbeck) con un tiempo de correlación comparable al periodo orbital del binario.
   • Filtro (Filter): El minidisco actúa como un filtro de paso bajo, amortiguando las fluctuaciones más rápido que la escala de tiempo dinámica en su borde exterior. Este filtrado produce el endurecimiento adicional de $f^{-2}$, resultando en la pendiente total de $f^{-4}$ en altas frecuencias.
5. Independencia de la Viscosidad: Se encontró que la forma de la PSD, incluyendo la frecuencia de ruptura y la pendiente de alta frecuencia, es insensible al parámetro de viscosidad $\alpha$ a lo largo de un orden de magnitud, confirmando que la ruptura está determinada por escalas de tiempo dinámicas en lugar de la difusión viscosa.

Significancia y Reclamaciones
El artículo afirma que la PSD de ley de potencia rota con una pendiente de alta frecuencia de $-4$ sirve como un potencial signo distintivo del continuo para binarios de SMBH sub-parsec. Esta morfología difiere de los modelos estándar de paseo aleatorio amortiguado aplicados comúnmente a AGN simples (que típicamente predicen pendientes cercanas a $-2$).

Los autores proponen que si ya se sospecha que un AGN alberga un binario (por ejemplo, mediante candidatos de periodicidad o desplazamientos de velocidad de líneas anchas), la detección de una ruptura de continuo cerca de $5 f_{bin}$ con una pendiente pronunciada de $-4$ apoyaría la hipótesis del binario. Por el contrario, la ausencia de tal ruptura pesaría en contra de ella.

Sin embargo, el artículo ejerce una cautela significativa respecto a la conjetura más amplia de que toda la variabilidad de los AGN surge de una binaridad generalizada. Si bien las PSD simuladas se asemejan a la variabilidad estocástica observada en muchos AGN, los autores señen que las frecuencias de ruptura observadas en los AGN (que corresponden a periodos orbitales de años a décadas) situarían a estos sistemas en la banda de frecuencia de ondas gravitacionales detectable por arreglos de temporización de pulsares. Dadas las restricciones actuales sobre el fondo estocástico de ondas gravitacionales, es poco probable que una gran fracción de los AGN luminosos sean binarios en estos periodos específicos, ya que esto sobreproduciría el fondo observado. Por lo tanto, aunque el mecanismo es un signo distintivo viable para candidatos específicos, es probable que no explique la variabilidad de toda la población de AGN.

Finalmente, el estudio establece un criterio numérico para la metodología de simulación: para capturar con precisión la variabilidad física de la acreción binaria, las simulaciones deben emplear tamaños de sumidero lo suficientemente pequeños ($s \lesssim 0.1$) para resolver el efecto de amortiguación del minidisco, evitando el ruido artificial de alta frecuencia generado por sumideros excesivamente grandes.