Interaction between disk and extended corona in a general relativistic framework

Technische Zusammenfassung: Interaktion zwischen Disk und erweiterter Korona in einem allgemeinen relativistischen Rahmen (Teil I)

Problemstellung
Aktuelle Modelle von zweiphasigen Akkretionsflüssen in Röntgentbinären (XRBs) haben oft Schwierigkeiten, die Beobachtungen im harten Zustand zu erklären, in dem die optisch dicke Akkretionsscheibe (Disk) nahe an das Schwarze Loch (BH) reicht und von einer hart röntgenemittierenden Korona beleuchtet wird. Vorherige Studien (z. B. Haardt & Maraschi 1991; Poutanen et al. 2018) etablierten, dass eine statische Slab-Korona über einer Disk eine Obergrenze für die Comptonisierte Leistung besitzt, die sie im Energiegleichgewicht erzeugen kann. Wenn die Korona zu optisch dick oder zu heiß wird, tritt eine Runaway-Kaskade auf, bei der eine verstärkte Bestrahlung der Disk die Anzahl der Seed-Photonen erhöht, was wiederum die koronale Leistung weiter steigert und ein konsistentes Gleichgewicht verhindert. Zudem nehmen bestehende Modelle oft ein lokales Energiegleichgewicht an, vernachlässigen radiale Variationen sowie allgemeinrelativistische Effekte oder nutzen Codes, die keine Rückkopplung von der Korona auf die Disk berücksichtigen. Diese Arbeit adressiert die Lücke bei der Erzielung einer global konsistenten Gleichgewichtslösung für ein Disk-Korona-System innerhalb eines vollen GR-Rahmens, indem sie spezifisch die Lebensfähigkeit einer statischen Slab-Geometrie für harte Zustände testet.

Methodik
Die Autoren haben den Monte-Carlo-Strahlungstransfer-Code MONK (Zhang et al. 2019) aufgerüstet, um die wechselseitige Interaktion zwischen der Akkretionsscheibe und einer erweiterten Korona zu implementieren. Die wesentlichen methodischen Fortschritte umfassen:

1. Iterativer globaler Energiebilanz: Der Code berechnet iterativ die Comptonisierung durch die Korona und die Bestrahlung der Disk, bis ein globales Gleichgewicht erreicht ist. Die fundamentale Randbedingung ist, dass die gesamte beobachtete Leuchtkraft ($L_{tot}$) der gesamten verfügbaren Akkretionsleistung ($L_{acc} = \eta \dot{M}c^2$) entsprechen muss.
2. Albedo- und Feedback-Mechanismus: Ein Albedo-Parameter wird eingeführt, um den Bruchteil des einfallenden koronalen Flusses zu definieren, der von der Disk reflektiert wird (modelliert als semi-infinites Elektronengas nach dem Chandrasekhar-Formalismus). Der verbleibende Anteil wird absorbiert, was die effektive Temperatur und die intrinsische Dissipation der Disk erhöht.
3. Viskoser Dissipationsparameter ($\alpha$): Im Gegensatz zu früheren Arbeiten, die annahmen, dass die gesamte Akkretionsenergie in der Korona dissipiert wird, passt dieses Modell den Bruchteil $\alpha$ der viskosen Leistung an, die intrinsisch in der Disk dissipiert wird. Der Code iteriert $\alpha$ so, dass $L_{tot} = L_{acc}$ gilt. Wenn $L_{tot} > L_{acc}$, wird $\alpha$ reduziert, um den Seed-Photonenfluss zu senken, und umgekehrt.
4. Allgemeinrelativistischer Rahmen: Die Photonenpropagation wird mittels Nullgeodäten in der Kerr-Geometrie verfolgt. Der Code berücksichtigt Lichtbeugung, gravitative Rotverschiebung und die vier möglichen Photonenbahnen: das Erreichen von Unendlich, das Auftreffen auf die Korona, die Selbstbestrahlung der Disk oder der Fall in das Schwarze Loch.
5. Spektralanalyse: Der Photonenindex ($\Gamma$) der Comptonisierten Komponente wird durch Anpassung des emittierten Spektrums an ein Potenzgesetz ($N(E) \propto E^{-\Gamma}$) innerhalb eines Energiebereichs abgeleitet, der durch den Überlapp mehrerer Streuordnungen definiert ist.

Wesentliche Ergebnisse
Unter Verwendung von Referenzparametern (BH-Masse $10 M_\odot$, Spin $a=0$, Koronatemperatur $T_e = 120$ keV, Albedo $= 0,5$) liefert die Studie die folgenden Ergebnisse:

• Gleichgewichtslösungen: Es existieren globale Gleichgewichtslösungen sowohl für harte als auch für intermediäre Zustände. Für eine statische Slab-Korona, die die Disk vollständig bedeckt, liegt der niedrigste erreichbare Photonenindex unter typischen Bedingungen bei etwa 1,7–1,8.
• Optische Dicke-Limit: Es gibt eine strikte Obergrenze für die koronale optische Dicke ($\tau_{max} \approx 0,16$ für den Referenzfall), jenseits derer keine physikalische Gleichgewichtslösung existiert (was ein $\alpha < 0$ erfordern würde). Mit zunehmendem $\tau$ muss der intrinsische Disk-Dissipationsanteil $\alpha$ sinken, um das Energiegleichgewicht aufrechtzuerhalten, bis er schließlich gegen Null geht, wenn die Disk-Emission vollständig durch die korale Bestrahlung ersetzt wird.
• Parameterabhängigkeiten:
  • BH-Spin: Ein höherer Spin ($a=0,998$) führt zu härteren Spektren ($\Delta \Gamma \sim -0,15$) aufgrund erhöhter Lichtbeugung und Disk-Bestrahlung, was ein niedrigeres $\alpha$ erzwingt.
  • Koronatemperatur: Eine höhere $T_e$ führt zu härteren Spektren.
  • Albedo: Eine höhere Albedo härtet die Spektren signifikant auf. Bei einer Albedo von 1,0 (vollständige Reflexion) kann das System wesentlich höhere optische Dicken ($\tau \approx 0,6$) aufrechterhalten, was den Photonenindex auf $\Gamma \sim 1,0$ absinken lässt. Dies geschieht, weil eine hohe Albedo die reprozessierten weichen Seed-Photonen unterdrückt und somit die Feedback-Schleife entkoppelt.
• Spektrale Form: Das Comptonisierte Spektrum weicht oft von einem einfachen Potenzgesetz ab und zeigt einen „Anisotropie-Bruch“ bei hohen Energien, ein Merkmal, das auf die einseitige Injektion von Seed-Photonen und GR-Effekte zurückzuführen ist.
• Radialer Energiefluss: Unter der Annahme einer räumlich uniformen, erweiterten Slab-Korona erfordert das Erreichen eines lokalen Gleichgewichts in jedem radialen Ring zusätzlich zum globalen Gleichgewicht einen Energietransport von den inneren zu den äußeren Ringen.

Bedeutung und Schlussfolgerungen
Die Autoren kommen zu dem Schluss, dass es mit dem aufgerüsteten MONK-Code nun möglich ist, ein globales Energiegleichgewicht für Disk-Korona-Systeme zu erreichen, was einen robusten Rahmen zur Bestimmung von Systemparametern bietet. Die Studie stellt jedoch fest, dass eine statische Slab-Korona keine günstige Geometrie ist, um die härtesten beobachteten Spektren in XRBs zu erklären, selbst wenn ein globales Energiegleichgewicht erzwungen wird. Die Unfähigkeit der statischen Slab, Photonenindizes signifikant unter 1,7–1,8 zu erzeugen (es sei denn, die Albedo beträgt 1,0), deutet darauf hin, dass alternative Geometrien erforderlich sind. Die Autoren geben an, dass zukünftige Arbeiten abgeschnittene Disk-Geometrien (truncated disks) und ausströmende Koronae (outflowing coronae) als potenzielle Lösungen untersuchen werden. Der Ansatz bietet zudem das Potenzial, die Koronageometrie zusätzlich zum Fluss auch mithilfe von Polarisationsdaten einzugrenzen.