Stochastic Variability of Binary Accretion

Sintesi Tecnica: Variabilità Stocastica dell'Accrezione Binaria

Definizione del Problema
I binari di buchi neri supermassicci (SMBH) sono una conseguenza prevista della formazione gerarchica delle galassie, ma rimangono osservativamente elusivi. Sebbene la ricerca di una modulazione quasi-periodica del flusso vicino alla frequenza orbitale del binario sia stata perseguita come un segnale primario, tali ricerche sono complicate dalla natura intrinsecamente stocastica delle curve di luce degli AGN e limitate da basi temporali dei survey che escludono la stragrande maggioranza dei binari sub-parsec con lunghi periodi orbitali. Inoltre, i tassi di recupero dei periodogrammi standard diminuiscono significativamente per le forme d'impulso non sinusoidali previste dalle simulazioni idrodinamiche. Di conseguenza, è necessario identificare segnali temporali nella densità spettrale di potenza (PSD) dei tassi di accrezione che siano distinti dalla canonica variabilità a rumore rosso dei singoli AGN, concentrandosi specificamente sulla variabilità stocastica a scale temporali più brevi del periodo orbitale del binario.

Metodologia
Gli autori impiegano simulazioni idrodinamiche 2D ad altissima risoluzione utilizzando il codice Sailfish per modellare l'accrezione di gas in un binario circolare di massa diseguale ($q = M_2/M_1 = 0.2$) circondato da un disco di gas circumbinario. Le simulazioni utilizzano un'equazione di stato localmente isotermica e una prescrizione di viscosità di Shakura-Sunyaev ($\alpha = 0.01$) con un numero di Mach orbitale $M=10$.

Un aspetto critico della metodologia è la variazione sistematica del parametro della dimensione del "sink" (pozzo), che determina i raggi di troncamento interno dei minidischi componenti ($r_{sink} \propto s \cdot a$). Il parametro $s$ è variato da 0.05 a 0.40. Gli autori distinguono tra il flusso di massa che entra nel minidisco (il "driver") e il tasso di accrezione di massa sul sink del buco nero (il "filtro"). Le simulazioni vengono eseguite attraverso una fase di burn-in (2900 orbite) seguita da una fase di campionamento (100 orbite) ad alta risoluzione ($\Delta \ln r = 0.0025$), con serie temporali campionate $\sim 10^3$ volte per orbita per risolvere frequenze ben al di sopra della frequenza orbitale.

Risultati Chiave

1. Forma della PSD del Continuo: Nel limite di dimensioni del sink piccole ($s \lesssim 0.1$), la densità spettrale di potenza del tasso di accrezione converge a una forma universale a legge di potenza spezzata. La PSD è piatta (rumore bianco) alle basse frequenze, transisce verso una pendenza di $-2$ in una banda stretta e si inclina verso una pendenza di $-4$ alle alte frequenze.
2. Frequenza di Break: La transizione verso la ripida pendenza di $-4$ avviene a una frequenza di break $f_d \approx 5$ volte la frequenza orbitale del binario ($f_{bin}$). Questa frequenza corrisponde alla frequenza Kepleriana locale al bordo esterno del minidisco del secondario ($r_{out} \approx 0.18\text{--}0.20 a$).
3. Dipendenza dalla Dimensione del Sink: Per sink grandi ($s \gtrsim 0.2$), la PSD non converge alla forma universale. Inveve, il minidisco agisce come un buffer insufficiente, permettendo a un eccesso di potenza ad alta frequenza di raggiungere il sink, manifestandosi come picchi nel tasso di accrezione. Gli autori sostengono che questi picchi siano artefatti numerici quando il raggio del sink supera l'orbita circolare interna più stabile (ISCO) fisica del buco nero.
4. Meccanismo Driver-Filtro: La variabilità è spiegata come un processo stocastico a due stadi:
   • Driver: Gli stream circumbinari consegnano il gas al bordo esterno del minidisco in un modo che assomiglia a un processo di cammino casuale smorzato (processo di Ornstein-Uhlenbeck) con un tempo di correlazione comparabile al periodo orbitale del binario.
   • Filtro: Il minidisco agisce come un filtro passa-basso, smorzando le fluttuazioni più velocemente della scala temporale dinamica al suo bordo esterno. Questo filtraggio produce l'ulteriore inclinazione $f^{-2}$, risultando in una pendenza totale $f^{-4}$ alle alte frequenze.
5. Indipendenza dalla Viscosità: La forma della PSD, inclusa la frequenza di break e la pendenza ad alta frequenza, si trova essere insensibile al parametro di viscosità $\alpha$ su un ordine di grandezza, confermando che la break è determinata dalle scale temporali dinamiche piuttosto che dalla diffusione viscosa.

Significatività e Rivendicazioni
L'articolo sostiene che la PSD a legge di potenza spezzata con una pendenza ad alta frequenza di $-4$ funga da potenziale segnale del continuo per i binari SMB sub-parsec. Questa morfologia differisce dai modelli standard di cammino casuale smorzato spesso applicati ai singoli AGN (che tipicamente prevedono pendenze vicine a $-2$).

Gli autori propongono che se un AGN è già sospettato di ospitare un binario (ad esempio, tramite candidati di periodicità o spostamenti di velocità delle righe larghe), la rilevazione di una break nel continuo vicino a $5 f_{bin}$ con una ripida pendenza di $-4$ supporterebbe l'ipotesi del binario. Al contrario, l'assenza di tale break peserebbe contro di essa.

Tuttavia, il documento esercita una cautela significativa riguardo alla congettura più ampia secondo cui tutta la variabilità degli AGN derivi da una diffusa binarità. Sebbene le PSD simulate somiglino alla variabilità stocastica osservata in molti AGN, gli autori notano che le frequenze di break osservate negli AGN (corrispondenti a periodi orbitali di anni o decenni) porrebbero questi sistemi nella banda di frequenza delle onde gravitazionali rilevabili dagli array di timing dei pulsar. Date le attuali restrizioni sullo sfondo di onde gravitazionali stocastiche, è improbabile che una grande frazione di luminosi AGN siano binari a questi specifici periodi, poiché ciò produrrebbe un eccesso del fondo osservato. Pertanto, mentre il meccanismo è un segnale valido per candidati specifici, probabilmente non spiega l'intera popolazione di AGN nella sua variabilità.

Infine, lo studio stabilisce un criterio numerico per la metodologia di simulazione: per catturare accuratamente la variabilità fisica dell'accrezione binaria, le simulazioni devono impiegare dimensioni del sink sufficientemente piccole ($s \lesssim 0.1$) per risolvere l'effetto di buffering del minidisco, evitando il rumore artificiale ad alta frequenza generato da sink sovradimensionati.