Interaction between disk and extended corona in a general relativistic framework

Résumé Technique : Interaction entre le disque et la couronne étendue dans un cadre de relativité générale (Partie I)

Énoncé du problème
Les modèles actuels de flux d'accrétion à deux phases dans les binaires X (XRB) peinent souvent à expliquer les observations en état dur (hard state) où le disque optiquement épais s'étend près du trou noir (BH), illuminé par une couronne émettant des rayons X durs. Des études antérieures (par exemple, Haardt & Maraschi 1991 ; Poutanen et al. 2018) ont établi qu'une couronne de type nappe statique au-dessus d'un disque possède une limite supérieure sur la puissance Comptonisée qu'elle peut produire en équilibre énergétique. Si la couronne devient trop épaisse ou trop chaude, une cascade incontrôlée se produit : l'augmentation de l'illumination du disque accroît le flux de photons de semence, ce qui booste davantage la production de la couronne, empêchant ainsi un équilibre auto-cohérent. De plus, les modèles existants supposent souvent un équilibre énergétique local, négligeant les variations radiales et les effets de relativité générale (RG), ou reposent sur des codes qui n'intègrent pas le feedback de la couronne vers le disque. Ce travail comble la lacune consistant à atteindre une solution d'équilibre globalement auto-cohérente pour un système disque-couronne dans un cadre complet de RG, en testant spécifiquement la viabilité d'une géométrie de nappe statique pour les états durs.

Méthodologie
Les auteurs ont mis à niveau le code de transfert radiatif de Monte Carlo MONK (Zhang et al. 2019) pour incorporer l'interaction mutuelle entre le disque d'accrétion et une couronne étendue. Les avancées méthodologiques clés incluent :

1. Équilibre Énergétique Global Itératif : Le code calcule désormais de manière itérative la Comptonisation par la couronne et l'illumination du disque jusqu'à ce qu'un équilibre global soit atteint. La contrainte fondamentale est que la luminosité totale observée ($L_{tot}$) doit être égale à la puissance d'accrétion totale disponible ($L_{acc} = \eta \dot{M}c^2$).
2. Albédo et Mécanisme de Feedback : Un paramètre d'albédo est introduit pour définir la fraction du flux coronal incident réfléchie par le disque (modélisé comme une atmosphère électronique semi-infinie selon le formalisme de Chandrasekhar). La fraction restante est absorbée, augmentant la température effective et la dissipation intrinsèque du disque.
3. Paramètre de Dissipation Visqueuse ($\alpha$) : Contra-irement aux travaux précédents qui supposaient que toute l'énergie d'accrétion était dissipée dans la couronne, ce modèle ajuste la fraction $\alpha$ de la puissance visqueuse dissipée intrinsèquement dans le disque. Le code itère sur $\alpha$ de telle sorte que $L_{tot} = L_{acc}$. Si $L_{tot} > L_{acc}$, $\alpha$ est réduit pour abaisser le flux de photons de semence, et inversement.
4. Cadre de Relativité Générale : La propagation des photons est suivie via des géodésiques nulles dans une géométrie de Kerr. Le code prend en compte la déflexion de la lumière, le décalage vers le rouge gravitationnel et les quatre trajectoires possibles des photons : atteindre l'infini, frapper la couronne, auto-irradier le disque, ou tomber dans le trou noir.
5. Analyse Spectrale : L'indice de photon ($\Gamma$) de la composante Comptonisée est dérivé en ajustant le spectre émergent à une loi de puissance ($N(E) \propto E^{-\Gamma}$) dans une plage d'énergie définie par le chevauchement de multiples ordres de diffusion.

Résultats Clés
En utilisant des paramètres de référence (masse du BH $10 M_\odot$, spin $a=0$, température coronale $T_e = 120$ keV, albédo $= 0,5$), l'étude produit les résultats suivants :

• Solutions d'Équilibre : Des solutions d'équilibre global existent pour les états durs et intermédiaires. Pour une couronne de nappe statique couvrant entièrement le disque, l'indice de photon le plus bas atteignable est d'environ 1,7–1,8 dans des conditions typiques.
• Limite d'Épaisseur Optique : Il existe une limite supérieure stricte sur l'épaisseur optique coronale ($\tau_{max} \approx 0,16$ pour le cas de référence) au-delà de laquelle aucune solution d'équilibre physique n'existe (nécessitant un $\alpha < 0$). À mesure que $\tau$ augmente, la fraction de dissipation intrinsèque du disque $\alpha$ doit diminuer pour maintenir l'équilibre énergétique, atteignant des valeurs proches de zéro où l'émission du disque est entièrement issue de la réémission de l'illumination coronale.
• Dépendances de Paramètres :
  • Spin du BH : Un spin plus élevé ($a=0,998$) conduit à des spectres plus durs ($\Delta \Gamma \sim -0,15$) en raison de l'augmentation de la déflexion de la lumière et de l'illumination du disque, ce qui force un $\alpha$ plus faible.
  • Température Coronale : Une $T_e$ plus élevée produit des spectres plus durs.
  • Albédo : Un albédo élevé durcit considérablement les spectres. Pour un albédo de 1,0 (réflexion complète), le système peut soutenir des épaisseurs optiques beaucoup plus élevées ($\tau \approx 0,6$), permettant à l'indice de photon de descendre à $\Gamma \sim 1,0$. Cela se produit car un albédo élevé supprime les photons de semence mous reprocessés, découplant ainsi la boucle de feedback.
• Forme du Spectre : Le spectre Comptonisé dévie souvent d'une simple loi de puissance, présentant une « rupture d'anisotropie » aux hautes énergies, une caractéristique attribuée à l'injection unidirectionnelle des photons de semence et aux effets de la RG.
• Flux d'Énergie Radial : Sous l'hypothologie d'une couronne de nappe étendue et spatialement uniforme, atteindre un équilibre local à chaque anneau radial nécessite un flux d'énergie des anneaux internes vers les anneaux externes, en plus de satisfaire l'équilibre global.

Signification et Conclusions
Les auteurs concluent qu'avec la version améliorée du code MONK, il est désormais possible d'atteindre un équilibre énergétique global pour les systèmes disque-couronne, fournissant un cadre robuste pour contraindre les paramètres du système. Cependant, l'étude démontre qu'une couronne de nappe statique n'est pas une géométrie favorable pour expliquer les spectres les plus durs observés dans les XRB, même lorsque l'équilibre énergétique global est imposé. L'incapacité de la nappe statique à produire des indices de photon significativement inférieurs à 1,7–1,8 (à moins que l'albédo ne soit de l'ordre de l'unité) suggère que des géométries alternatives sont nécessaires. Les auteurs déclarent que les travaux futurs exploreront les géométries de disque tronqué et les couronnes sortantes comme solutions potentielles à ces limitations. L'approche permet également de contraindre la géométrie de la couronne en utilisant les données de polarisation en plus du flux.