Stochastic Variability of Binary Accretion

Resumo Técnico: Variabilidade Estocástica de Acreção Binária

Definição do Problema
Binários de buracos negros supermassivos (SMBH) são uma consequência prevista da formação hierárquica de galáxias, embora permaneçam observacionalmente elusivos. Embora a modulação de fluxo quase periódica próxima à frequência orbital do binário tenha sido perseguida como um sinalizador primário, tais buscas são complicadas pela natureza intrinsecamente estocástica das curvas de luz de Núcleos Ativos de Galáxias (AGN) e limitadas por bases de dados de levantamento que excluem a vasta maioria dos binários sub-parsec com períodos orbitais longos. Além disso, as taxas de recuperação de periodogramas padrão caem significamente para as formas de pulso não sinusoidais previstas por simulações hidrodinâmicas. Consequentemente, há uma necessidade de identificar sinais temporais na densidade espectral de potência (PSD) contínua das taxas de acreção que sejam distintos da variabilidade canônica de ruído vermelho de um AGN único, focando especificamente na variabilidade estocástica em escalas de tempo menores que o período orbital do binário.

Metodologia
Os autores empregam simulações de hidrodinâmica 2D de altíssima resolução utilizando o código Sailfish para modelar a acreção de gás em um binário circular de massas desiguais ($q = M_2/M_1 = 0.2$) cercado por um disco de gás circumbinário. As simulações utilizam uma equação de estado localmente isotérmica com um número de Mach orbital $M=10$ e uma prescrição de Shakura-Sunyaev $\alpha$ para viscosidade ($\alpha = 0.01$).

Um aspecto crítico da metodologia é a variação sistemática do parâmetro de tamanho do "sink" (sumidouro), $s$, que determina os raios de truncamento interno dos minidiscos componentes ($r_{sink} \propto s \cdot a$). O parâmetro $s$ é variado de 0,05 a 0,40. Os autores distinguem entre o fluxo de massa que entra no minidisco (o "driver" ou indutor) e a taxa de acreção de massa sobre o sumidouro do buraco negro (o "filter" ou filtro). As simulações são executadas através de uma fase de burn-in (2900 órbitas) seguida por uma fase de amostragem (100 órbitas) em alta resolução ($\Delta \ln r = 0.0025$), com séries temporais amostradas $\sim 10^3$ vezes por órbita para resolver frequências bem acima da frequência orbital.

Principais Resultados

1. Forma da PSD Contínua: No limite de tamanhos de sumidouro pequenos ($s \lesssim 0.1$), a densidade espectral de potência da taxa de acreção converge para uma forma universal de lei de potência quebrada. A PSD é plana (ruído branco) em baixas frequências, transiciona para uma inclinação de $-2$ em uma banda estreita e torna-se mais íngreme para uma inclinação de $-4$ em altas frequências.
2. Frequência de Quebra: A transição para a inclinação íngreme de $-4$ ocorre em uma frequência de quebra $f_d \approx 5$ vezes a frequência orbital do binário ($f_{bin}$). Esta frequência corresponde à frequência Kepleriana local na borda externa do minidisco do secundário ($r_{out} \approx 0.18\text{--}0.20 a$).
3. Dependência do Tamanho do Sumidouro: Para sumidouros grandes ($s \gtrsim 0.2$), a PSD não converge para a forma universal. Em vez disso, o minidisco atua como um buffer insuficiente, permitindo que o excesso de potência de alta frequência alcance o sumidouro, manifestando-se como picos na taxa de acreção. Os autores argumentam que esses picos são artefatos numéricos quando o raio do sumidouro excede a órbita circular mais interna estável (ISCO) física do buraco negro.
4. Mecanismo Driver-Filter: A variabilidade é explicada por um processo estocástico de dois estágios:
   • Driver: Os fluxos circumbinários entregam gás à borda externa do minidisco de uma maneira que se assemelha a um processo de caminhada aleatória amortecida (processo de Ornstein-Uhlenbeck) com um tempo de correlação comparável ao período orbital do binário.
   • Filter: O minidisco atua como um filtro passa-baixa, amortecendo flutuações mais rápido do que a escala de tempo dinâmica em sua borda externa. Este filtragem produz o adicional declive de $f^{-2}$, resultando na inclinação total de $f^{-4}$ em altas frequências.
5. Independência de Viscosidade: A forma da PSD, incluindo a frequência de quebra e a inclinação de alta frequência, é encontrada como insensível ao parâmetro de viscosidade $\alpha$ ao longo de uma ordem de magnitude, confirmando que a quebra é definida por escalas de tempo dinâmicas em vez de difusão viscosa.

Significância e Alegações
O artigo afirma que a PSD de lei de potência quebrada com uma inclinação de alta frequência de $-4$ serve como um potencial sinalizador contínuo para binários de SMBH sub-parsec. Esta morfologia difere dos modelos padrão de caminhada aleatória amortecida frequentemente aplicados a AGN únicos (que tipicamente preveem inclinações próximas de $-2$).

Os autores propõem que, se um AGN já é suspeito de hospedar um binário (por exemplo, via candidatos de periodicidade ou deslocamentos de velocidade de linhas largas), a detecção de uma quebra de contínuo próxima a $5 f_{bin}$ com uma inclinação íngreme de $-4$ apoiaria a hipótese do binário. Inversamente, a ausência de tal quebra pesaria contra ela.

No entanto, o artigo exerce cautela significativa em relação à conjectura mais ampla de que toda a variabilidade de AGN surge de binaridade generalizada. Embora as PSDs simuladas se assemelhem à variabilidade estocástica observada em muitos AGN, os autores observam que as frequências de quebra observadas em AGN (correspondendo a períodos orbitais de anos a décadas) colocariam esses sistemas na banda de frequência de ondas gravitacionais detectável por arrays de temporização de pulsares. Dadas as restrições atuais sobre o fundo estocástico de ondas gravitacionais, é improvável que uma grande fração dos AGN luminosos sejam binários nesses períodos específicos, pois isso superproduziria o fundo observado. Assim, embora o mecanismo seja um sinalizador viável para candidatos específicos, ele provavelmente não explica a variabilidade de toda a população de AGN.

Finalmente, o estudo estabelece um critério numérico para a metodologia de simulação: para capturar com precisão a variabilidade física da acreção binária, as simulações devem empregar tamanhos de sumidouro suficientemente pequenos ($s \lesssim 0.1$) para resolver o efeito de buffering do minidisco, evitando o ruído artificial de alta frequência gerado por sumidouros excessivamente grandes.