Interaction between disk and extended corona in a general relativistic framework

Sintesi Tecnica: Interazione tra Disco e Corona Estesa in un Framework Relativistico Generale (Parte I)

Enunciato del Problema
Gli attuali modelli di flussi di accrezione a due fasi nelle binarie a raggi X (XRB) spesso faticano a spiegare le osservazioni dello stato "hard" (duro), in cui il disco otticamente spesso si estende vicino al buco nero (BH), illuminato da una corona che emette raggi X duri. Studi precedenti (ad esempio, Haardt & Maraschi 1991; Poutanen et al. 2018) hanno stabilito che una corona a slab statica sopra un disco ha un limite superiore sulla potenza Comptonizzata che può produrre in equilibrio energetico. Se la corona diventa troppo otticamente spessa o calda, si verifica un cascata incontrollata in cui l'incremento dell'illuminazione del disco aumenta i fotoni di seed, potenziando ulteriormente l'output coronale, impedendo un equilibrio autosufficiente. Inoltre, i modelli esistenti spesso assumono un equilibrio energetico locale, trascurando le variazioni radiali e gli effetti di relatività generale (GR), o si affidano a codici che non incorporano il feedback dalla corona al disco. Questo lavoro affronta la lacuna nel raggiungere una soluzione di equilibrio globalmente autosufficiente per un sistema disco-corona all'interno di un pieno framework GR, testando specificamente la viabilità di una geometria a slab statica per gli stati hard.

Metodologia
Gli autori hanno aggiornato il codice di trasferimento radiativo Monte Carlo MONK (Zhang et al. 2019) per incorporare l'interazione reciproca tra il disco di accrezione e una corona estesa. I principali avanzamenti metodologici includono:

1. Equilibrio Energetico Globale Iterativo: Il codice calcola ora iterativamente la Comptonizzazione della corona e l'illuminazione del disco fino al raggiungimento di un equilibrio globale. Il vincolo fondamentale è che la luminosità totale osservata ($L_{tot}$) deve essere uguale alla potenza di accrezione totale disponibile ($L_{acc} = \eta \dot{M}c^2$).
2. Meccanismo di Albedo e Feedback: Viene introdotto un parametro di albedo per definire la frazione di flusso coronale incidente riflessa dal disco (modellato come un'atmosfera elettronica semi-infinita utilizzando il formalismo di Chandrasekhar). La frazione rimanente viene assorbita, aumentando la temperatura effettiva e la dissipazione intrinseca del disco.
3. Parametro di Dissipazione Viscosa ($\alpha$): A differenza dei lavori precedenti che assumevano che tutta l'energia di accrezione fosse dissipata nella corona, questo modello regola la frazione $\alpha$ della potenza viscosa dissipata intrinsecamente nel disco. Il codice itera $\alpha$ affinché $L_{tot} = L_{acc}$. Se $L_{tot} > L_{acc}$, $\alpha$ viene ridotto per abbassare il flusso di fotoni di seed, e viceversa.
4. Framework di Relatività Generale: La propagazione dei fotoni è tracciata tramite geodetiche nulle in geometria di Kerr. Il codice tiene conto della curvatura della luce (light bending), del redshift gravitazionale e delle quattro possibili traiettorie dei fotoni: raggiungere l'infinito, colpire la corona, auto-irradiare il disco o cadere nel BH.
5. Analisi Spettrale: L'indice del fotone ($\Gamma$) della componente Comptonizzata è derivato adattando lo spettro emergente a una legge di potenza ($N(E) \propto E^{-\Gamma}$) entro un intervallo di energia definito dalla sovrapposizione di molteplici ordini di scattering.

Risultati Chiave
Utilizzando parametri di riferimento (massa del BH $10 M_\odot$, spin $a=0$, temperatura coronale $T_e = 120$ keV, albedo $= 0.5$), lo studio produce i seguenti risultati:

• Soluzioni di Equilibrio: Esistono soluzioni di equilibrio globale sia per gli stati hard che per quelli intermedi. Per una corona a slab statica che copre interamente il disco, l'indice del fotone minimo raggiungibile è approssimativamente 1.7–1.8 in condizioni tipiche.
• Limite di Spessore Ottico: Esiste un limite superiore rigoroso sulla profondità ottica coronale ($\tau_{max} \approx 0.16$ per il caso di riferimento) oltre il quale non esiste alcuna soluzione di equilibrio fisica (richiedendo $\alpha < 0$). All'aumentare di $\tau$, la frazione di dissipazione intrinseca del disco $\alpha$ deve diminuire, raggiungendo infine valori prossimi allo zero dove l'emissione del disco è interamente dovuta alla rielaborazione dell'illuminazione coronale.
• Dipendenze dai Parametri:
  • Spin del BH: Uno spin più elevato ($a=0.998$) porta a spettri più duri ($\Delta \Gamma \sim -0.15$) a causa dell'aumento della curvatura della luce e dell'illuminazione del disco, che forza un $\alpha$ più basso.
  • Temperatura Coronale: Una $T_e$ più alta produce spettri più duri.
  • Albedo: Un'albedo più alta indurisce significativamente gli spettri. Per un'albedo di 1.0 (riflessione completa), il sistema può sostenere profondità ottiche molto più elevate ($\tau \approx 0.6$), permettendo all'indice del fotone di scendere a $\Gamma \sim 1.0$. Ciò avviene perché un'alta albedo sopprime i fotoni di seed morbidi rielaborati, disaccoppiando il ciclo di feedback.
• Forma Spettrale: Lo spettro Comptonizzato spesso devia da una semplice legge di potenza, mostrando un "break di anisotropia" ad alte energie, una caratteristica attribuita all'iniezione unidirezionale di fotoni di seed e agli effetti di GR.
• Flusso di Energia Radiale: Sotto l'ipotesi di una corona a slab estesa e spazialmente uniforme, raggiungere l'equilibrio locale in ogni anello radiale richiede che l'energia fluisca dagli anelli interni a quelli esterni, oltre a soddisfare l'equilibrio globale.

Significatività e Conclusioni
Gli autori concludono che, con l'aggiornamento del codice MONK, è ora possibile raggiungere l'equilibrio energetico globale per i sistemi disco-corona, fornendo un framework robusto per vincolare i parametri del sistema. Tuttavia, lo studio rileva che una corona a slab statica non è una geometria favorevole per spiegare gli spettri più duri osservati nelle XRB, anche quando viene imposto l'equilibrio energetico globale. L'incapacità della slab statica di produrre indici di fotone significativamente inferiori a 1.7–1.8 (a meno di un'albedo unitaria) suggerisce che siano necessarie geometrie alternative. Gli autori affermano che i lavori futuri esploreranno geometrie di disco troncato e corone in outflow come potenziali soluzioni a queste limitazioni. L'approccio ha anche il potenziale per vincolare la geometria coronale utilizzando i dati di polarizzazione oltre al flusso.