Interaction between disk and extended corona in a general relativistic framework

Technische Samenvatting: Interactie tussen Schijf en Uitgebreide Corona binnen een Algemeen Relativistisch Kader (Deel I)

Probleemstelling
Huidige modellen van twee-fasen accretiestromen in röntgen-dubbelsterren (XRBs) hebben vaak moeite om de observaties in de harde staat te verklaren, waarbij de optisch dikke schijf dicht bij het zwarte gat (BH) doorloopt, geïllumineerd door een hard röntgen-emitterende corona. Eerdere studies (bijv. Haardt & Maraschi 1991; Poutanen et al. 2018) stelden vast dat een statische plaatvormige (slab) corona boven een schijf een bovengrens heeft voor de Comptonized energie die het kan produceren in energetisch evenwicht. Als de corona te optisch dik of te heet wordt, treedt er een runaway-cascade op waarbij verhoogde illuminatie van de schijf leidt tot meer zaadfotonen, wat de corona-output verder stimuleert, waardoor een zelfconsistent evenwicht wordt verhinderd. Bovendien gaan bestaande modellen vaak uit van lokaal energiebalans, waarbij radiale variaties en algemene relativistische (GR) effecten worden genegeerd, of vertrouwen ze op codes die geen feedback van de corona naar de schijf bevatten. Dit werk adresseert de kloof in het bereiken van een globaal zelfconsistent evenwichtssolutie voor een schijf-corona systeem binnen een volledig GR-kader, waarbij specifiek de levensvatbaarheid van een statische slab-geometrie voor harde staten wordt getest.

Methodologie
De auteurs hebben de Monte Carlo radiatieve transfercode MONK (Zhang et al. 2019) geüpgraded om de wederzijdse interactie tussen de accretieschijf en een uitgebreide corona te incorporeren. De belangrijkste methodologische vooruitgang omvat:

1. Iteratieve Globale Energiebalans: De code berekent nu iteratief de Comptonisatie door de corona en de schijf-illuminatie totdat een globaal evenwicht is bereikt. De fundamentele beperking is dat de totale waargenomen luminositeit ($L_{tot}$) gelijk moet zijn aan de totale beschikbare accretie-energie ($L_{acc} = \eta \dot{M}c^2$).
2. Albedo en Feedbackmechanisme: Er is een albedo-parameter geïntroduceerd om het fractie van de invallende corona-flux te definiëren die door de schijf wordt gereflecteerd (gemodelleerd als een semi-oneindige elektronenatmosfeer met behulp van Chandrasekhar's formalisme). Het resterende deel wordt geabsorbeerd, wat de effectieve temperatuur en intrinsieke dissipatie van de schijf verhoogt.
3. Viskeuze Dissipatieparameter ($\alpha$): In tegenstelling tot eerdere werken die aannamen dat alle accretie-energie in de corona wordt gedissipeerd, past dit model het fractie $\alpha$ aan van de viskeuze energie die intrinsiek in de schijf wordt gedissipeerd. De code iterereert $\alpha$ zodat $L_{tot} = L_{acc}$. Als $L_{tot} > L_{acc}$, wordt $\alpha$ verlaagd om de zaadfotonflux te verlagen, en vice versa.
4. Algemeen Relativistisch Kader: Fotonpropagatie wordt gevolgd via nul-geodesieken in de Kerr-geometrie. De code houdt rekening met lichtbuiging, gravitationele roodverschuiving en de vier mogelijke fotonbanen: het bereiken van oneindigheid, het raken van de corona, het zelf-irradiëren van de schijf, of het vallen in het zwarte gat.
5. Spectrale Analyse: De fotonindex ($\Gamma$) van het Comptonized component wordt afgeleid door het emergente spectrum te fitten aan een machtswet ($N(E) \propto E^{-\Gamma}$) binnen een energiebereik dat wordt gedefinieerd door de overlap van meerdere verstrooiingsordes.

Belangrijkste Resultaten
Gebruikmakend van referentieparameters (BH massa $10 M_\odot$, spin $a=0$, corona temperatuur $T_e = 120$ keV, albedo $= 0.5$), levert de studie de volgende bevindingen op:

• Evenwichtssoluties: Er bestaan globale evenwichtssoluties voor zowel harde als intermediaire staten. Voor een statische slab corona die de schijf volledig bedekt, ligt de laagst haalbare fotonindex ongeveer op 1.7–1.8 onder typische condities.
• Optische Dikte Limiet: Er is een strikte bovengrens aan de corona-optische dikte ($\tau_{max} \approx 0.16$ voor de referentiecasus) waarboven geen fysieke evenwichtssolutie bestaat (wat een $\alpha < 0$ vereist). Naarmate $\tau$ toeneemt, moet de intrinsieke schijf-dissipatiefractie $\alpha$ afnemen om het energiebalans te handhaven, waarbij uiteindelijk bijna-nulwaarden worden bereikt waarbij de schijfemissie volledig wordt gereproduceerd door corona-illuminatie.
• Parameterafhankelijkheden:
  • BH Spin: Een hogere spin ($a=0.998$) leidt tot hardere spectra ($\Delta \Gamma \sim -0.15$) door toegenomen lichtbuiging en schijf-illuminatie, wat een lagere $\alpha$ afdwingt.
  • Corona Temperatuur: Een hogere $T_e$ resulteert in hardere spectra.
  • Albedo: Een hogere albedo hardt de spectra significant uit. Voor een albedo van 1.0 (volledige reflectie) kan het systeem veel hogere optische diktes handhaven ($\tau \approx 0.6$), waardoor de fotonindex kan dalen tot $\Gamma \sim 1.0$. Dit komt doordat een hoge albedo de gereproduceerde zachte zaadfotonen onderdrukt, waardoor de feedbackloop ontkoppeld wordt.
• Spectrale Vorm: Het Comptonized spectrum wijkt vaak af van een eenvoudige machtswet, waarbij een "anisotropie-breuk" bij hoge energieën zichtbaar is, een kenmerk dat wordt toegeschreven aan de eenzijdige injectie van zaadfotonen en GR-effecten.
• Radiale Energiestroom: Onder de aanname van een ruimtelijk uniforme, uitgebreide slab corona, is het vereist dat energie van de binnenste naar de buitenste ringen stroomt om lokaal evenwicht op elke radiale annulus te bereiken, naast het voldoen aan het globale evenwicht.

Betekenis en Conclusies
De auteurs concluderen dat het met de geüpgradede MONK-code mogelijk is om globaal energie-evenwicht te bereiken voor schijf-corona systemen, wat een robuust kader biedt voor het beperken van systeemparameters. Echter, de studie stelt vast dat een statische slab corona geen gunstige geometrie is voor het verklaren van de hardste geobserveerde spectra in XRBs, zelfs wanneer een globaal energiebalans wordt afgedwongen. Het onvermogen van de statische slab om fotonindices significant lager dan 1.7–1.8 te produceren (tenzij de albedo eenheid is), suggereert dat alternatieve geometrieën vereist zijn. De auteurs stellen dat toekomstig werk zich zal richten op geëxtraheerde (truncated) schijfgeometrieën en uitstromende (outflowing) corona's als potentiële oplossingen voor deze beperkingen. De aanpak biedt ook potentieel voor het beperken van de corona-geometrie met behulp van polarisatiegegevens naast flux.