Issue 
A&A
Volume 633, January 2020



Article Number  A35  
Number of page(s)  17  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201935033  
Published online  06 January 2020 
Warm and thick corona for a magnetically supported disk in galactic black hole binaries
Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, 00716 Warsaw, Poland
email: gronki@camk.edu.pl
Received:
8
January
2019
Accepted:
4
September
2019
Context. We selfconsistently model a magnetically supported accretion disk around a stellarmass black hole with a warm optically thick corona based on first principles. We consider the gas heating by magnetorotational instability dynamo.
Aims. Our goal is to show that the proper calculation of the gas heating by magnetic dynamo can build up the warm optically thick corona above the accretion disk around a black hole of stellar mass.
Methods. Using the vertical model of the disk supported and heated by the magnetic field together with radiative transfer in hydrostatic and radiative equilibrium, we developed a relaxation numerical scheme that allowed us to compute the transition form the disk to corona in a selfconsistent way.
Results. We demonstrate here that the warm (up to 5 keV) optically thick (up to 10 τ_{es}) Comptoncooled corona can form as a result of magnetic heating. A warm corona like this is stronger in the case of the higher accretion rate and the greater magnetic field strength. The radial extent of the warm corona is limited by local thermal instability, which purely depends on radiative processes. The obtained coronal parameters are in agreement with those constrained from Xray observations.
Conclusions. A warm magnetically supported corona tends to appear in the inner disk regions. It may be responsible for soft Xray excess seen in accreting sources. For lower accretion rates and weaker magnetic field parameters, thermal instability prevents a warm corona, giving rise to eventual clumpiness or ionized outflow.
Key words: accretion / accretion disks / Xrays: binaries / radiative transfer / instabilities / magnetic fields / methods: numerical
© ESO 2020
1. Introduction
There is growing evidence that a warm optically thick corona exists in accreting black holes when an accretion disk is present. It is observed in different types of accreting black holes across various masses: active galactic nuclei (AGN; Pounds et al. 1987; Magdziarz et al. 1998; Mehdipour et al. 2011; Done et al. 2012; Petrucci et al. 2013, 2018; Keek & Ballantyne 2016) including quasars (Madau 1988; Laor et al. 1994, 1997; Gierliński & Done 2004; Piconcelli et al. 2005), ultraluminous Xray sources (ULXs; Goad et al. 2006; Stobbart et al. 2006; Gladstone et al. 2009), and galactic black hole binaries (GBHBs; Gierliński et al. 1999; Zhang et al. 2000; Di Salvo et al. 2001). A warm corona is visible as a soft component (below 2 keV) in the Xray spectra of these objects. It is generally characterized by an excess with respect to the extrapolation of a hard Xray power law. This power law typically originates from a region called hot corona.
The observed spectral shape of the soft Xray excess provides us some insight into the radiative cooling mechanism that operates in the warm corona. Assuming that the thermal Comptonization is the dominant cooling process, we can model the observed spectra with socalled slab model, where soft photons from the disk enter the warm corona that is located above the disk and undergo Compton scattering with amplification factor y. During the fitting procedure, the temperature and optical depth of a warm corona can be determined. Observations of many sources indicate that a corona like this has an optical depth ranging from 4 up to 40 in different objects (Zhang et al. 2000; Jin et al. 2012; Petrucci et al. 2013, 2018, and references therein).
Furthermore, observations constrain the amount of energy that is dissipated in corona, f, in comparison to the total energy released in the whole disk or corona system. Observations of GBHBs and AGNs show that f ≈ 1 is required to explained the observed spectral index α = 0.9 in some sources (Haardt & Maraschi 1991; Życki et al. 2001; Petrucci et al. 2013). Petrucci et al. (2018) set a lower limit of f = 0.8. This behavior is not universal because it depends on the object and its spectral state, but it shows that in some circumstances, much energy is released outside of an accretion disk.
The more fundamental question is how this warm optically thick slab of gas can be created above an accretion disk that has a lower temperature. The additional process that heats the warm corona and maintains it in steady state with disk is also unclear.
Many attempts have been made to determine the universal mechanism of energy dissipation in the warm corona, but the problem is still not fully solved. A strong corona, which is responsible for most of the thermal energy that is released increases the disk stability (Svensson & Zdziarski 1994; Kusunose & Mineshige 1994; Begelman & Pringle 2007). However, detailed computations of radiative transfer in illuminated accretion disk atmospheres show that the outer warm or hot skin cannot be optically thick and stable at the same time (Ballantyne et al. 2001; Nayakshin & Kallman 2001; Różańska et al. 2002, 2015; Madej & Różańska 2004). When the irradiation increases, the ionized skin becomes unstable and the gas most probably flows out in the form of wind (Proga & Waters 2015). Computations show that the warm or hot skin with an optical depth higher than 3 cannot be thermally stable (in pressure equillibrium) with a cold accretion disk when the skin is heated only radiatively (Krolik & Kriss 1995; Nayakshin et al. 2000; Różańska et al. 2002, 2015, and references therein).
Additional heating of the warm corona layer by an accretion process has also been considered. However, Shakura & Sunyaev (1973) introduced the geometrically thindisk model in their seminal paper, where the kinetic energy of accreting gas is locally converted into thermal energy. According to this, within the standard disk model all energy is dissipated deep inside the disk at the equatorial plane. Vertical profiles of energy dissipation clearly show the exponential decrease toward the disk surface.
The development of numerical simulations has shown the role of the magnetic field in the structure of the accretion disk (Hawley & Balbus 1991; Balbus & Hawley 1991; Tout & Pringle 1992; Hawley 2001; Stone et al. 2008). Many magnetohydrodynamic (MHD) simulations are available, but none of them specifies the optical depth of the warm corona because most of them do not take radiative cooling into account (e.g., Bai & Stone 2013; Penna et al. 2013). The simulations that contain simplified radiative cooling (e.g., Turner 2004; Hirose et al. 2009; Ohsuga & Mineshige 2011; Noble et al. 2011) only discuss a hot corona without constrains on its observational properties (i.e., Jiang et al. 2014; Sądowski et al. 2017). Highly advanced radiative transfer calculations are usually made a posteriori after general relativistic MHD simulations are completed (Schnittman et al. 2016).
Różańska et al. (2015, hereafter RMB15) have developed an analytic model for an optically thick uniformly heated Comptoncooled corona in hydrostatic equilibrium. The heating mechanism was not specified in the model, but it allowed the authors to integrate vertically both disk and warm corona, coupled by mechanical heating and radiative cooling. This allowed them to determine the relations between the optical depth of the warm layer and the diskcorona energy budget. Furthermore, the authors have shown that when some part of the gas pressure is replaced by magnetic pressure, the optical depth of the stable corona increases.
We here assume a geometrically thin and optically thick accretion disk around a black hole (see Fig. 1). We also assume that the accretion disk is magnetized and the magnetorotational instability (MRI) is the primary source of viscosity and energy dissipation. Based on these assumptions, we develop the model of a slablike warm corona that covers the accretion disk and adopt realistic heating of the gas by magnetic field reconnection, which in contrast to the standard disk increases toward the disk surface (Hirose et al. 2006). We directly use the analytic formula derived by Begelman et al. (2015, hereafter BAR15), where the vertical profile of the magnetic heating of the accretion flow is determined. In addition to this major assumption, the vertical structure of the disk together with the radiative transfer equation in a gray atmosphere are fully solved with the relaxation method proposed by Henyey et al. (1964). We have developed a new numerical scheme that allows solving nonlinear differential equations in a relatively short time, which keeps the integration error low. Furthermore, with our new method we are able to pass through the thermally unstable regions that may arise when nonuniform heating and cooling mechanisms take place.
Fig. 1. Schematic illustration of a slice through the disk plane (only one side is shown due to symmetry). The black hole is shown as the circle on the left. The gray area represents the optically thick and geometrically thin disk where most of the accreting mass is located. The warm corona that covers the disk, in which the magnetic energy is released as radiation, is marked in orange. Magnetic and radiative energy flux are represented by blue and red arrows, respectively. The part of the corona that might collapse as a result of thermal instability is shown in magenta. The inner optically thin and geometrically thick hot flow that we did not consider is drawn in yellow for completion. 

Open with DEXTER 
As a result, we show for which range of parameters a warm corona can exist in case of a GBHB. We determine the optical depth of the warm layer, which is the main observable when Xray data are analyzed. Detailed calculations of the vertical structure allow us to estimate the amount of energy that is dissipated in the corona in comparison to the total energy that is released by accretion at a given radius. We show the radial structure of an optically thick corona and give tight constraints for conditions for which this warm layer can exist. All results are compared with observations.
A warm corona is not the only possible model with which the soft Xray excess can be explained. Other models assume relativistically smeared absorption (Gierliński & Done 2004) or reflection from an accretion disk that is illuminated by a hard Xray continuum (García et al. 2019) or bulk Comptonization in a central Compton cloud (Shaposhnikov & Titarchuk 2006; Seifina et al. 2016). All of these models provide good fits to Xray observations by incorporating an additional Comptonized component, but they do not explain the physical mechanism that feeds energy into the warm gas to compensate for the huge rate of the Compton cooling. Models that do provide the heating mechanism for a central compact corona exist and can successfully explain some observations without the need of involving the magnetic field (Seifina et al. 2018). However, we note that the role of the magnetic field in the accretion and jet formation process has been shown numerous times in the papers cited in this section, and we do not consider it a strong assumption. We show that an accretion disk with an operating magnetic dynamo is able to produce enough energy to heat the extended surface layer, whose physical parameters are consistent with those determined from spectral fitting of observations. Our results give a possible physical explanation for the phenomenon that was in most cases only considered from an observational perspective without analyzing the global energy budget. In order to verify whether our model can correctly reproduce the spectral features seen in Xray sources, spectral modeling involving the redistribution function of Compton scattering must be performed. This is beyond the scope of this paper, however.
The structure of the paper is as follows: Sect. 2 presents the formalism of a magnetically heated corona and the diskcorona radiative transfer equation. It ends with a full description of the differential equations and the relaxation method we used to solve them. Section 3 presents the results of our numerical computations. We display the vertical structure of the diskcorona system, but we also show the radial limitation within which an optically thick corona can exist. Discussion and conclusions are given in Sect. 4.
2. Model setup
2.1. Magnetically heated corona
Begelman & Pringle (2007) and later BAR15 have proposed a vertical model of magnetically supported disks (MSDs) that are driven by the MRI dynamo. The MRI dynamo operating near the equatorial plane in the presence of an external poloidal magnetic field generates the toroidal magnetic flux, which converts the mechanical energy of the accretion into the electromagnetic energy that is stored as tangled magnetic field lines. The total rate of this energy deposition per unit height is α_{B}ΩP_{tot}, where α_{B} is the toroidal magnetic field production parameter, Ω is the Keplerian angular velocity, and P_{tot} is the total pressure, which is the sum of gas, radiation, and magnetic pressure. Although the toroidal field production parameter α_{B} does resemble the effective viscosity parameter α = t_{rϕ}/P (known from the classical thindisk model) in the sense that it scales the energy release within the disk, its physical meaning is different, even though the expected values of both parameters are similar (Salvesen et al. 2016, hereafter SSAB16, Table 3).
The magnetic flux ropes that form near the midplane rise buoyantly toward the surface, and in the model by BAR15, the velocity is very coarsely approximated by
where 0 < η < 1 is a dimensionless parameter. A similar process of a buoyant emergence of magnetic fields toward the surface is observed on the Sun, but the dynamical properties and the origin of the magnetic field are vastly different.
During the upward motion of the ropes, the magnetic field decays by various processes, releasing its accumulated energy into heating the gas. When we assume a purely toroidal flux that rises vertically with speed v_{B}(z), the energy dissipation rate due to the induced current, ℋ_{mag}, can be derived from the elemental electrodynamics and yields
This expression is consistent with intuitive reasoning: the energy loss is proportional to the magnetic energy density gradient times velocity.
Another process that may heat the gas is magnetic reconnection due to cyclic reversals of the toroidal field polarity that is observed in simulations of the MRI dynamo. The energy output of this process is proportional to the magnetic pressure P_{mag} and can be estimated as
where ξ is a proportionality constant (Appendix A in BAR15). The rate of thermal energy release peaks in the proximity of the disk photosphere (τ = 1), but its exact distribution may change depending on the model parameters. Alongside ξ, we also define, after BAR15, the reconnection efficiency parameter ν, which describes the ratio of reconnection process to the total magnetic viscosity: ν = 2ξ/α_{B}, and it is very convenient to present our results (see Sect. 3.5).
One consequence of this picture is that the field outflow eventually reaches the point where heating becomes inefficient and some portion of the accretion energy may be carried away by the Poynting flux and is never released to heat the gas. This is contrary to nonmagnetic classical (neglecting advection) αviscosity models where all accretion power is released as radiation. The magnetic outflow may be associated with the matter outflow, but this association is not trivial: because the outflow speed v_{B} is lower than the freefall speed v_{ff} = Ωz, most of the matter will probably slide down the magnetic flux tubes, and only a small amount is carried away by the outflow. In contrast, during the mass ejections in the solar corona, the eruption timescale is much shorter than the dynamical timescale, which results in much more efficient mass transport.
Numerical modeling shows that the process described above can be very complex, but this analytic approximation enables us to build and study the atmospheres of MSDs, which is the purpose of this paper.
Following BAR15, the accretion energy is injected into the Poynting flux and consumed for gas heating, in varying proportion. When we express the Poynting flux as F_{mag} = 2P_{mag}v_{B}, the magnetic energy conservation equation can be written as
We substitute velocity and different energy rates by Eqs. (1)–(3) to obtain the final expression for the magnetic pressure gradient,
Finally, we can use this equation to obtain the expression for the heating rate ℋ [erg cm^{−3} s^{−1}] without using gradients,
This heating rate describes the amount of released magnetic energy that heats the gas and can then be converted into radiation. In this paper we neglect advection and assume that the system is fully cooled by radiation, as specified in the next section.
2.2. Radiative transfer in the diskcorona system
We consider a gray planeparallel optically thick atmosphere, which is additionally heated by dissipation of the magnetic field. We include freefree process and Compton scattering in our model and neglect synchrotron radiation, which is not important in this regime of magnetic pressure (see Sect. 3.1 for exact numbers). To determine the temperature structure, we consider three equations: a radiative equilibrium equation, a transfer equation, and an Eddington approximation.
The radiative equilibrium equation for a smallmatter volume heated by various dissipation processes at rate ℋ and by the local radiation field of mean intensity J_{ν} and cooled by the radiation emissivity j_{ν} reads
where j_{bol} ≡ ∫_{ν}j_{ν}dν, and ρ is gas density in g cm^{−3}. Here, we assumed that the entire dissipated thermal energy is locally converted into radiation (neglecting advection).
In an optically thick regime, when true absorption processes dominate scattering, the radiation field is fully thermalized and strictly local, therefore for matter at temperature T, we have J_{ν} = B_{ν}(T), where B_{ν}(T) is the Planck function for temperature T. However, because a corona is a strongly scattering medium, this is not the case (particularly near the surface), and the actual radiation field intensity may deviate from a Planck distribution toward a Wien distribution. The further deviation from the Planck law is caused by noncoherent inverse Compton scattering. We still assume that the spectrum of mean radiation intensity J_{ν} within the disk can be approximately described by the Planck function with temperature T_{rad}, that is, J_{ν} ≈ B_{ν}(T_{rad}). With these assumptions, we can define our frequencyaveraged quantities
Even when we assume that the radiation field has a Planck distribution, we do not equal gas and radiation temperatures because the corona is strongly dominated by scattering.
For the gray atmosphere we considered in this paper, we used both Rosseland and Planck averages, the latter denoted with a “P” superscript. The following values for electron scattering and freefree opacities were assumed:
where is a constant from Kramers opacity approximation in cgs units. Through this paper, we denote the total Rosseland mean opacity as κ = κ_{es} + κ_{ff}, while is the total Planck opacity.
Using the gray approximation, we can rewrite the radiative equilibrium Eq. (7) as
We refer to the function Λ_{rad}(ρ, T, T_{rad}) as the net radiative cooling rate.
To estimate the effect of inverse Compton scattering, we used the standard formula given by Rybicki & Lightman (2008), where by J_{C} we mean all incident photons that were scattered with thermalized electrons of the gas temperature T:
In this formula is the energy passed to the radiation from electrons by the inverse Compton scattering. The mean photon energy ⟨ε⟩ and mean squared energy ⟨ε^{2}⟩ depend on the photon energy distribution, and their ratio has the following values:
Although we assumed that the radiation has a Planck spectrum everywhere in the atmosphere, for a strongly scattering medium with a Comptoncooled corona, the radiation spectrum is shifted toward slightly higher energies, and we find the approximation ⟨ε^{2}⟩ = 4kT_{rad} ⋅ ⟨ε⟩ satisfactory for our model. When the temperatures are low, density is high, and freefree absorption dominates, this term is negligible and does not change the result. Therefore we replace J_{C} = κ_{es}J to obtain the Compton term for the emission function,
Finally, by taking into account all relevant opacities, we obtain the emission function j_{bol} of the following form:
We can now substitute the expressions for B, J, and j_{bol} into the cooling function Λ_{rad} to obtain its final form,
The net cooling rate can be split into two terms: Λ_{rad} = Λ_{B} + Λ_{C}, describing the contributions of bremsstrahlung and Compton scattering,
We solve here the frequencyintegrated radiative transfer equation (Eq. (4) in RMB15), where the zeroth moment of this equation reads
When we average this over the frequencies and replace the Eddington flux with a physical flux 4πH_{Edd} = F_{rad}, we obtain the expected result that the energy gained by the radiative flux is equal to the rate of radiative gas cooling,
Because we are interested in the optically thick regime, we assume the Eddington approximation, that is, J = 3 K, where K is the second moment of specific intensity. We therefore solve the first moment of the radiative transfer equation as
Taking into account Eq. (9) and the connection of the Eddington flux with the physical flux, this equation yields
Radiative transfer should be always coupled with the gas structure. We assume that the whole diskcorona system is in hydrostatic equilibrium, that is,
where the gas pressure typically is P_{gas} = k/μm_{H}Tρ, and the first part on the right side denotes the radiation pressure P_{rad}. This equation is consistently implemented into a full set of equations that is solved simultaneously from the midplane to the surface of the corona. This treatment proves that we do not treat the diskcorona system as a twoslab model, but solve the full gas structure together with proper heating processes and radiative cooling taken into account.
2.3. Equation set and boundary conditions
We solved the following set of five equations: (5), (12), (21), (23), and (24), using the expression for the heating rate given by Eq. (6) and that for the cooling function given by Eq. (17). The full numerical procedure of our numerical code is presented in Appendix A. We integrated them from the disk equatorial plane (z = 0) to the upper bound of the computational range (z = z_{max}). We estimated z_{max} by solving Eq. (38) in BAR15 for the height at which the magnetic pressure reaches p = 10^{−5} of the magnetic pressure at the midplane,
where H_{SS73} is the disk height derived from the αdisk model (Shakura & Sunyaev 1973) and
The requirement of hydrostatic equilibrium gives the balance between all pressures: magnetic, radiation, and gas pressure, and the gravity force in vertical direction. We determined the gas temperature by solving the balance equation between magnetic heating and Compton cooling in order to determine whether and in which circumstances the corona can form.
We show that the problem is better solved when other radiative processes, such as freefree emission, are taken into account. We analyzed two cases: case A, when Λ_{rad} ≡ Λ_{B} + Λ_{C} is given by the full Eq. (17) above, and case B, when Λ_{rad} ≡ Λ_{C} and thus the freefree term is set to zero, and the atmosphere is only cooled by Compton scattering.
In the case B, Eq. (17) is linear and we obtain exactly one temperature solution in each point of the vertical structure, as opposed to case A, for which multiple solution can exist. However, in case B, the cooling rate is decreased, therefore the gas temperature in the transition region between the disk and corona is overestimated. Despite this, because the density in the corona depends mostly on magnetic pressure gradient and becomes bounded to the gas pressure only in extremely weakly magnetized cases, case B allows us to obtain almost the exact density value compared to case A.
We adopted the following procedure: We solved the entire model using case B. Now to correct the temperature in the transition region where case B performs worst, we solved the Eq. (21) again, but kept the density fixed, which guarantees a unique solution. Because we corrected the temperature but did not update the density or gas pressure, we deviate slightly from hydrostatic equilibrium. However, the corona is dominated by the magnetic field pressure gradient, and this deviation is negligible and does not change the overall structure of the atmosphere. At this little cost, we are able to overcome the instability and investigate the influence of disk magnetization on its presence and strength.
We adopted boundary conditions appropriate for a cylindrical geometry of the accretion disk. The radiative flux carried through the midplane (z = 0) must vanish due to symmetry,
Following BAR15, we also assumed that the magnetic pressure gradient at the equatorial plane must be zero. Therefore, from Eq. (5) we obtain
If the radiation pressure is neglected, the relation derived by BAR15 between the parameters α_{B}, η, ν, and the magnetic parameter (β = P_{gas}/P_{mag}) at the equatorial plane β_{0} is
At the top of the atmosphere (z → ∞), we assumed that the sum of the flux carried away by radiation and of the magnetic field is equal to the flux obtained by Keplerian disk theory,
where G is the gravitational constant, M is the black hole mass, Ṁ is the accretion rate, and R is the radial distance from black hole. Typically, R_{Schw} denotes the Schwarzschild radius, which is given as 2GM/c^{2}. We also used the standard boundary condition for radiative transfer, J = 2H, which guarantees that the disk is not externally illuminated. Expressed in our convention, the boundary condition at z_{max} takes the form
This gives a total of four boundary conditions, which complete our set of five equations because Eq. (12) is not a differential equation, and Eqs. (5), (21), (23), and (24) are firstorder nonlinear ordinary differential equations.
Our numerical computations are parameterized by six input quantities: three magnetic, α_{B}, η and ν, and three disk parameters: black hole mass, radial distance, and the accretion rate. For better convenience, we used the accretion rate ṁ = Ṁκ_{es}c/48πGM, given in units of the Eddington rate with an accretion efficiency for a Newtonian potential.
3. Vertical structure of the disk with corona
3.1. Parameter regime
Throughout this paper we consider the case of GBHB with M = 10 M_{⊙}. The results are presented for a wide range of accretion rates from ṁ = 10^{−3} up to 1 and for radial distances from a marginally stable orbit up to R/R_{Schw} = 400. The value of the disk outer radius we adopted was taken from the fact that for larger distances corona does not exist at all or it is very weak.
The parameters α_{B}, η, and ν used in our computations have been determined by SSAB16 from numerical simulations for several runs of a varying net poloidal magnetic flux value, which resulted in different degrees of disk magnetization. Their results are plotted in Fig. 2 by points.
Fig. 2. Dependence of model parameters: toroidal field production efficiency, α_{B} (lower panel), reconnection efficiency, ν (middle panel), and the magnetic field gradient parameter, (upper panel) on the magnetic parameter in the disk midplane β_{0}. Points denote values obtained from numerical simulations in SSAB16 and lines are our results from the powerlaw relation between α_{B} and η (Eq. (33)) for four values of ξ. 

Open with DEXTER 
For strongly magnetized case, simulations show that α_{B} ≈ η ≃ 0.3, β_{0} ≃ 1 and ξ ≃ 0, which gives the reconnection efficiency parameter ν = 0. This means that the field production is efficient, the outflow is at a high rate, and dissipation by reconnection is low. On the other hand, for the weakly magnetized case, simulations result in α_{B} ≃ 0.02, η ≃ 0.04, β_{0} ≃ 50, and due to frequent dynamo reversals, ξ ≃ 0.5. For this case, the reconnection efficiency parameter ν ≃ 50. On the other hand, SSAB16 concluded that models with ν higher than a few tend to overestimate the magnetic field gradient for weakly magnetized models and that ν = 0 describes the magnetic pressure descent more accurately for these cases (see Fig. 13 from SSAB16). This means that it would be more consistent to assume the ν parameter to have relatively low values, on the order of a few, to be in agreement with the vertical shape of the magnetic pressure.
To reduce the number of magnetic parameters to consider, we assumed a constant value of ξ and made the simplification that α_{B} and η are tied by a power law relation. Using Eq. (29), we recomputed α_{B}, which is consistent with an assumed value of ξ and η^{sim} and taken from Tables 2 and 3 of the SSAB16 paper,
We then performed a powerlaw fit and found that for ξ = 0.05 the relation that best fits the six data points is
The dependence of the magnetic parameters on the use of these analytical formulae is plotted by lines in Fig. 2 for comparison with simulations. The figure shows that we can adopt the values of magnetic parameters, which are consistent with simulations.
The upper panel of Fig. 2 shows the magnetic pressure gradient, which is an important quantity influencing the magnetic energy release. The magnetic pressure gradient can be derived from Eq. (5) and yields
In the corona, when P_{gas} + P_{rad} ≪ P_{mag}, we can substitute P_{tot}/P_{mag} ≃ 1. The value of q is then constant, and when we express it in terms of model parameters only, we obtain a result that is compatible with Eq. (37) in BAR15:
With α_{B} remaining the only free magnetic parameter in our model, the value that we find most suitable for potentially interesting disks is α_{B} = 0.1. For this parameter, the radial distribution of the toroidal magnetic field at the equatorial plane B_{0} and at the base of the corona B_{cor} (where the gas temperature achieves minimum), depending on the accretion rate and the distance from black hole, is given in Fig. 3. The magnetic field in our computations never approached values of 10^{8} G.
Fig. 3. Toroidal magnetic field at the equatorial plane B_{0} and in the corona B_{cor} at optical depth τ_{cor} (see text for definition) in Gauss units. Logarithmic values of the magnetic field are plotted in the radius (R/R_{schw}) – accretion rate (ṁ) parameter plane. 

Open with DEXTER 
3.2. Characteristics of the vertical profile
The typical structure of the diskcorona system for ṁ = 0.05 and R/R_{Schw} = 10 is shown in Fig. 4. The results are presented with respect to the height above the disk midplane z/H, where half of the disk thickness is given by
Fig. 4. Vertical structure of the diskcorona system computed at the single radius R = 10 R_{Schw} for ṁ = 0.05. The horizontal axis is the height above the disk midplane given in scale height H defined by Eq. (36). Three significant locations are marked with vertical dotted lines: the photosphere (green), the temperature minimum (dark yellow), and the thermalization depth (red). The following panels display (a) radiation and gas temperature; (b) gas, radiation, magnetic, and total pressure; (c) gas density and pressure gradients; (d) radiative energy flux, magnetic (Poynting) flux, and gas heating rate all normalized to the total accretion flux F_{acc}; (e) net radiative cooling terms (Λ_{B} and Λ_{C}) and gas magnetic heating rates; and (f) electron scattering and total and effective optical depths. 

Open with DEXTER 
Furthermore, to better present the results, we defined the total optical depth of the system according to the standard formula dτ = −κ^{P}ρ dz, where the electron scattering optical depth is denoted as dτ_{es} = −κ_{es}ρ dz, and the density number per unit of volume is n_{H} = ρ/m_{H}.
In our model the transition in vertical structure between the disk and warm corona occurs naturally, therefore we show the location of the photosphere τ = 2/3 and thermalization depth τ^{⋆} = 1, where the seed photons originate (see Rybicki & Lightman 2008, Sect. 1.7), by the vertical lines in Fig. 4. The thermalization depth is equivalent to the effective optical depth, which in our case was computed using Planck opacities (see Sect. 2.2). We took the temperature minimum as the inner boundary of the corona, and we refer to electron scattering optical depth of the temperature minimum as τ_{cor}, which is also marked by a vertical line in all panels of Fig. 4.
Although the geometrical proportions and physical parameters can vary when the model input changes, we almost universally obtain a structure in which three distinct regions can be observed that directly correspond to the gas pressure structure: (i) the accretion disk, where most of the matter is located, characterized by high density near the midplane, the gas temperature that exactly follows the radiation temperature, and domination of the gas and radiation pressure over the magnetic pressure; panels a and b. (ii) The transition region, where bremsstrahlung is still an important cooling mechanism and the gas temperature is coupled to the radiation field temperature, but the magnetic pressure takes over in supporting the structure and slows the density decrease with height; panel c. (iii) The corona, where the density is low and therefore the temperature must increase with height so that the disk remains in energy equilibrium; matter is cooled mainly by the inverse Compton scattering, while the density structure is fully supported by the magnetic pressure; panels b and c. In the case presented in Fig. 4, the following regions extend within 0 ≤ z/H ≲ 3, 3 ≲ z/H ≲ 7, and 7 ≲ z/H, respectively. A noticeable increase in temperature appears in the transition region, when the warm corona starts to form due to the magnetic heating.
The overall magnetic pressure profile (see panel b) is in agreement with that obtained by BAR15, which confirms the correctness of our calculations. The density profile of the outermost coronal layers is fully shaped by pressure gradients, of which the magnetic pressure gradient is the largest (see panel c).
When we solve the vertical structure of the MSD, the magnetic energy flux F_{mag} needs to be considered in addition to the radiative energy flux F_{rad}, the latter being the only energy transport channel in a classical αdisk. Both fluxes normalized to F_{acc} are shown in Fig. 4d and are equal to zero in the disk midplane. In contrast to the αdisk, where the heating is mostly located near the disk midplane, we obtained a heating distribution that peaks at a few scale heights, although still below the temperature inversion point; see Fig. 4e. This behavior is mostly the effect of purely magnetic heating ℋ_{mag}, which depends on the magnetic pressure gradient. Heating by reconnection ℋ_{rec} is proportional to the magnetic pressure (rather than its gradient) and peaks around the midplane.
The total optical depth of the diskcorona system is about 10^{4}. The differences between various optical depth profiles are visible in panel f of Fig. 4. For the above example, the whole corona is dominated by electron scattering, that is, τ_{es} = τ. The corona starts at τ_{cor} ∼ 10, and the total optical depth τ ∼ 20 where full thermalization occurs, for instance, τ^{⋆} = 1. The corona is optically thick, which agrees with recent observations.
3.3. Thermal instability of the intermediate layer
Thermal instability in the medium can occur if under given conditions the cooling rate does not increase with temperature. The thermal balance equation then has three solutions, and if the matter is on the unstable branch, a thermal runaway or collapse occurs, following a small deviation from equilibrium. The coronal thermal instability in accretion disks has been investigated for many years (Krolik et al. 1981; Kusunose & Mineshige 1994; Nakamura & Osaki 1993; Różańska & Czerny 1996; Różańska et al. 1999), but it has never been treated in the presence of magnetic field.
Considering the static vertical structure of the accretion disk, it is reasonable to assume that any thermal collapse is an isobaric process. The condition for thermal instability reads
When we split the cooling into Λ_{B} and Λ_{C} contributions and substitute the opacities given by Eqs. (10) and (11), we obtain
Depending on the contributions from these two terms, we distinguish tree cases: i when Λ_{B} ≪ Λ_{C}, the solution is always thermally stable because Eq. (39) has positive sign, ii when Λ_{B} ≫ Λ_{C}, instability occurs for
and finally, iii when Λ_{B} ≈ Λ_{C} and T > T_{rad}, radiative heating by bremsstrahlung is negligible, and the criterion for instability reads
In the isobaric regime, when δP_{gas} = 0 or ρT = const., this defines the upper limit for the temperature for which thermal instability occurs. The transition between these regimes ii and iii occurs for the density
In order to show the transition between stable and unstable corona, the net cooling function in MSD for three slightly varying accretion rates 0.35, 0.50, and 0.72 is presented in Fig. 5 for the radiation temperature T_{rad} = 10^{6}K. The three models vary mainly by the shift in density; the highest density corresponds to the lowest accretion rate.
Fig. 5. Logarithm of the net cooling rate Λ_{rad} (Eq. (17)), plotted in the n_{H} – T parameter plane ranging from blue (low rate) to red (high rate). The points of constant pressure are connected with black dashdotted lines. The red dotted contours show the contribution of the Compton cooling to the total cooling. The solid white line shows where the cooling function is constant in isobaric perturbation, and separates the area of unstable solutions where the gradient of the cooling function ℒ_{rad} < 0. Dark blue lines of increasing width represent the vertical structure for three accretion rates 0.035, 0.050, and 0.072. 

Open with DEXTER 
Each model starts at the disk midplane where T ≈ T_{rad}, and only when the transition to the corona starts does the gas temperature diverge from the radiation temperature. When the density is higher than that given by Eq. (42), the model will inevitably enter the zone defined by Eq. (40). Otherwise, for hotter disks, the contribution of Compton cooling will yield at least half of the total cooling, and the model might graze the instability boundary defined by Eq. (41) but never exceed it. This means that all but the Comptondominated disks will be affected by thermal instability of the transition region. The only stable solution is possible for relatively rare disks for the low cooling rate, that is, the lefthand side of the white thick solid line.
The same three models are shown Fig. 6. In order to identify all thermal equilibrium solutions for each height above the disk midplane, the fraction of Compton cooling, the energy balance, and the cooling function gradient ℒ_{rad} are shown depending on the gas temperature T while maintaining constant gas pressure. For the hottest model (highest accretion rate), there is an unambiguous solution in each point. This changes when the temperature decreases: although the main continuous solution is still stable (compared to Fig. 5), two cooler solutions appear. The middle solution lies within the unstable zone (pink area in the third column), but the coldest solution is stable. This means that because this area is dominated by magnetic pressure, an abrupt density change would not greatly affect the hydrostatic equilibrium and the existence of twophase matter is possible. Finally, for the lowest accretion rate (upper panels of Fig. 6), no numerically consistent stable solution exists. We expect that matter in this region would remain cold, but above the unstable region, a corona might still form. This corona would be hot and optically thin, and in the example we discussed, the transition would occur at about z/H ∼ 20.
Fig. 6. Onset of thermal instability in a lowtemperature disk, ṁ = 0.035 shown in the top panels, and its stabilization with the increase of the temperature, ṁ = 0.05 and 0.072 shown in the middle and bottom panels. Other disk parameters are as given in Sect. 3.1. In each panel the vertical axis represents the possible gas temperature T assuming constant gas pressure (δP_{gas} = 0), and the horizontal axis is the height above the disk midplane given in scale height H (Eq. (36)). The possible solutions that satisfy the radiative energy balance are shown with thick black contours, while the solution obtained from our numerical computations is additionally marked with light gray dots. A thinner gray line is used to show the radiation temperature, T_{rad}. Color maps in the first panels column show the contribution from Compton scattering to the total cooling. In the second column, the overall heatingcooling balance is shown: the yelloworange areas are where the heating dominates cooling, and blue areas show where the cooling rate exceeds the heating. The third column presents the stability diagnostic, i.e., the cooling function gradient ℒ_{rad} (Eq. (37)), with green areas marking the stable regions and violet areas the unstable regions. Two limit solutions according to Eqs. (40) and (41) are plotted with orange lines (dashdotted and dashed, respectively) and show excellent agreement with the numerical solution. Additionally, the content of neighboring plots (from right panels to the middle, and from middle to the left) are superimposed as thin contours for the sake of clarity. 

Open with DEXTER 
Equations (40) and (41) define the lower and upper limits for the temperature that allows the occurrence of thermal instability. This is clearly shown in Fig. 5, as each isobaric contour slices the instability boundary in exactly two points. These solutions are also plotted in Fig. 6 with orange lines. They indeed suffice to describe the instability in most of the vertical structure of the disk.
3.4. Local properties of the corona
The diskcorona transition at the given radius depends on the value of the local temperature and density. We present here the evolution of the selected disk parameters with an accretion rate and magnetization state with the aim to demonstrate better how the warm corona can arise at radius R = 10 R_{Schw}. For this purpose we define here two additional quantities related to the radiation pressure because the radiation energy is the main observable of the accretion diskcorona systems. The first such quantity is the socalled radiation efficiency, that is, the ratio of the radiation flux to the total flux that is dissipated locally at a certain point of the vertical structure due to magnetic field reconnection F_{rad}/F_{acc}, while the second quantity is the ionization parameter Ξ (Krolik et al. 1981), which is defined as
In Figs. 7 and 8 we show the structure of the temperature, the magnetic field gradient, the radiation efficiency, the cooling function gradient, the contribution of Compton cooling to the total cooling, and the ionization parameter in panels a to f, respectively. For all cases of the accretion rate we achieve a hot outermost layer at the top of the inner accretion disk.
Fig. 7. Vertical structure of the diskcorona models computed at R = 10 R_{Schw} for different values of the accretion rate (ṁ), plotted against the total optical depth τ. The panels show (a) the gas temperature, (b) the gradient of the magnetic pressure (Eq. (34)), (c) the radiative energy flux normalized to the total accretion flux F_{acc}, (d) the cooling function gradient (Eq. (37)), (e) the contribution of Compton cooling to the total cooling, and (f) the ionization parameter Ξ (Eq. (43)). 

Open with DEXTER 
Fig. 8. Same as in Fig. 7, but for different values of the magnetic parameter α_{B}. 

Open with DEXTER 
The temperature of the corona is very stable with respect to the accretion rate, but the optical depth for which the diskcorona transition occurs increases with accretion rate and reaches the value of ≈10 for ṁ = 0.71. In contrast, the accretion rate determines the density of the disk, and consequently, the value of the magnetic pressure gradient that saturates in the disk at the highest value for the lowest accretion. At the point of diskcorona transition, the magnetic pressure slightly decreases, but the rate of the radiative energy release is roughly the same, see panel c. We conclude that the primary factors causing the optical depth of the corona to vary are the changes in disk temperature (and the transition region, consequently).
Panel e partly explains this behavior: for lower accretion rates, cooling in the disk is dominated by bremsstrahlung, which switches off at around τ ∼ 1, causing an abrupt temperature gradient in the corona boundary. At the same time, the cooling function gradient also becomes negative for these cases (panel d). We conclude here that this contribution from bremsstrahlung also causes the formation of an unstable region, as discussed in Sect. 3.3. For higher accretion rates, the contribution from bremsstrahlung is not that significant, and the transition from the disk to the corona is almost seamless. The ionization parameter presented in panel f does not indicate a thermally unstable region. It shows the inversion with optical depth, which is shallower for Comptondominated cases (high accretion rate). Finally, Ξ is larger for hotter and rarer disks.
The changes in the corona are much more dramatic when the dependence on the magnetic field is examined, as shown in Fig. 8. Panel a shows that corona optical depth and temperature are radically different for weak and strong magnetic fields. For the weak field, the disk is much hotter and the radiation flux almost saturates to its maximum value well inside the disk (panel c), blue line, in contrast to the strong field case, where the disk is much colder, and more than half of the thermal energy is released in the corona above τ = 1 (panel c), red line. For the corona the temperature behavior is exactly opposite.
When magnetic pressure play a dominant role in the corona, it changes as , the density changes as , and the heating rate changes as . High values of q mean a very steep decline in magnetic field strength and gas heating rate, which in turn means that most of the energy will be dissipated within the disk, and the corona will be weak. This behavior is nicely demonstrated in panel b of Fig. 8. A stronger magnetic field adds to the effective disk cooling and consequently results in the appearance of thermal instability in diskcorona transition (panel d). Despite large temperature differences, Ξ within the disk remains almost constant (panel e), while in the corona it decreases with magnetization even if the coronal temperature is higher.
3.5. Global properties of the corona
To analyze the global behavior of our model within the parameter space, we define below several diagnostic quantities as proxies for various properties of the diskcorona system. We intended to select these quantities so that they can be compared to observables. The maps in Figs. 9–11 present these diagnostics depending on distance from the black hole R/R_{schw}, accretion rate ṁ, and α_{B} parameter.
Fig. 9. Properties of the corona on the radius (R/R_{schw}) – accretion rate (ṁ) parameter plane. All models are calculated for M_{BH} = 10M_{⊙}, α_{B} = 0.1, and ξ = 0.05. The panels display (a) the average temperature of the corona (Eq. (44)) in keV (colors) and τ_{cor} (contours); (b) the ionization parameter Ξ (colors) and magnetic β parameter (contours), both at the disk midplane; (c) the number density at τ = 2/3 (colors) and the total column density Σ in g cm^{2} (contours); (d) the fraction of radiative energy produced by the corona f (Eq. (45), colors), the maximum value of the magnetic field gradient q (Eq. (35), contours); (e) the y parameter of the warm corona (Eq. (46)), computed up to the thermalization layer τ^{⋆} = 1 (colors) and to τ_{cor} (contours); the (f) photon destruction probability (colors) and electron scattering optical depth (τ_{es}) (contours) at the thermalization layer; (g) the minimum value of ℒ_{rad} (Eq. (37)) throughout the disk height (colors), here, green areas are stable disks, and magenta areas indicate thermal instability in the corona, the gray dotted line corresponds to ℒ_{rad} = 0, and the ionization parameter Ξ at the corona base (contours); and (h) the radiative energy flux normalized to the total accretion flux (colors), and the disk scale ratio H/R (contours). 

Open with DEXTER 
Similarly to RMB15, and opposite to most models, our corona does not have a constant temperature, but changes with the optical depth, reaching temperature minimum at τ_{cor}. From Xray spectral fitting we measure some averaged temperature of the warm medium cooled by Comptonization, therefore it is useful to define the average temperature of the corona as
with the aim to compare it with observations. Both quantities T_{avg} and τ_{cor} are shown in each panel a of three maps. A warm optically thick corona tends to form above the inner strongly magnetized disk of high accretion rate.
Analogously to Haardt & Maraschi (1991), we denote the fraction of thermal energy released in the corona (versus total thermal energy) as
f = 0 corresponds to the passive corona, and f = 1 corresponds to passive disk. This parameter is shown in panel d and is tightly correlated with the magnetic field gradient q (Eq. (35)), plotted with contours. Because f is dependent only on magnetic parameters, we suggest that a strong magnetic field is required to obtain the high values of f that are reported by observations. Although we never approach f ≈ 1 and therefore cannot reproduce the passivedisk scenario discussed by Petrucci et al. (2013, 2018) for the case of AGN, we show that up to half of the total accretion energy can be dissipated in the corona.
Another useful quantity is the Compton y parameter, which can also be determined from observations. To derive the y value for our corona, we used the following definition, which incorporates an additional term that accounts for optically thick medium:
The value of this parameter is shown in panels e of Figs. 9–11 and can range between 0 and 0.5. This is consistent with the values obtained from observations of Xray binaries: y ≈ 0.3 (Życki et al. 2001), y ≈ 0.4 (Goad et al. 2006), and up to y ∼ 1 (Zhang et al. 2000).
All these quantities (τ_{cor}, T_{avg}, f, and y_{avg}) positively correlate both the with magnetic field strength and with the released accretion power. This means that the most prominent corona would be observed in the case of a high accretion rate and a strong field. However, above some threshold of disk magnetization (around α_{B} = 0.2, which corresponds to q ≈ 1.3), the trend reverses, and magnetic flux escapes without releasing all of its energy into radiation, as is shown in panel i of Figs. 10 and 11.
Fig. 10. Same as in Fig. 9, but on the radius (R/R_{schw}) – toroidal field production parameter (α_{B}) plane. 

Open with DEXTER 
Fig. 11. Same as Fig. 9, but on toroidal field production parameter (α_{B}) – accretion rate (ṁ) plane. 

Open with DEXTER 
In the case of GBHBs, where the disk temperatures are generally high, the optical depth of the corona τ_{cor} is not usually limited by a thermalization layer τ^{⋆} = 1 because the medium is strongly scattering dominated, as presented in panel f of Fig. 9. The electronscattering optical depth of this layer is usually ten times higher than τ_{cor} for the parameters considered here. Photons become thermalized due to high freefree absorption, as we show in comparison to the total Planck opacity in the same panel.
Even if the corona might exist as a solution, thermal instability will occur in many cases, as discussed in Sects. 3.3 and 3.4. Panels g of Figs. 10 and 11 show the stability of the warm corona solutions. The dependence on accretion power is very abrupt and almost switchlike: if the matter becomes too cold or too dense and the contribution of bremsstrahlung cooling accordingly increases, thermal instability occurs, and the corona cannot exist as a stable uniform layer (magenta regions). The same effect is obtained when the magnetic field strength increases and the disk becomes colder (compare with Fig. 8), but the transition is less sensitive and more gradual. The instability appears mostly in cold disk regions, that is, for a low accretion rate and a high magnetic field, and will remove only the optically thick skin when bremsstrahlung contribution is significant, and will most likely not affect the optically thin hot corona, which will still be present, even if its optical depth may vary with radius and accretion rate. This result, particularly the existence of the corona only around the zone of the maximum energy dissipation rate (around R/R_{schw} = 10) is consistent with Kusunose & Mineshige (1994).
The disk thickness d = H/R, where H is defined by Eq. (36), is shown in panels i of the figure maps. It ranges between 10^{−3} and 10^{−1}, which means that we still remain in thindisk regime, consistent with the assumption. It is important to note that the disk thickness changes mostly with the accretion rate. The increase in magnetic field strength does not inflate the disk (Fig. 11i) when the radiative pressure is included, as the fraction of energy that is dissipated within the disk decreases (panel d), and a disk core that is supported by radiative pressure (as in a classical αdisk) is not present, as shown by the radiative to gas pressure ratio Ξ in panel b. Similarly, the total column density Σ, shown by contours in panels c, is decreased because the magnetic fields take over the role of supporting the disk structure.
4. Discussion and conclusions
A disk with an optically thick corona is a frequently used satisfactory representation of the accretion sources in the soft state, and we selfconsistently modeled this structure here, based on first principles. Using a semianalytical static model of the magnetically supported disk (MSD) with radiation, we developed a new relaxation numerical code to calculate the diskcorona transition. The relaxation method designed by us and described in Appendix A is general (Henyey et al. 1964) and solves any set of linear differential equations with assumed boundary conditions. For the case of an accretion disk with a corona, we solved a set of five equations, four of which are differential and one is algebraic, with four boundary conditions (Sect. 2.3).
Our code fully solves the vertical structure, taking into account that part of the magnetic energy generated in the MRI process is converted into radiation and can be dissipated in the warm Comptoncooled corona. The remaining part of the magnetic flux is channeled into a toroidal magnetic field and escapes from the system, as introduced by BAR15. The original model proposed by BAR15 included no energy balance checking or radiative transfer, however, which allowed reducing the model to dimensionless equations. In this paper, we combined the MSD of these authors with solving the radiation transfer through the gray atmosphere in radiative and hydrostatic equilibrium. Additional mean opacities due to freefree absorption are taken into account.
In case of GBHBs with a black hole of stellar mass, we were able to reproduce both an optically thick and thin corona, which is commonly detectable in observations of accreting sources. In contrast to our previous paper (RMB15), where coronal heating was artificially assumed and kept constant with height, here we obtained a selfconsistent transition through all three layers: (i) a deep accretion disk up to a few scale heights that is dense and fully thermalized, (ii) a warm corona, where magnetic pressure slows the density decrease and the temperature slowly starts to increase, and (iii) a hot Corona where the density is low. Compton cooling is fully taken into account in all layers, but it strives against freefree cooling depending on the local physical gas conditions (see panel e in Figs. 7 and 8).
The base of the corona at temperature minimum τ_{cor} may reach 10 electron scattering optical depths, while the thermalization zone can be even 60 times τ_{es} (panels a and f in Figs. 9–11). This agrees with some observations of a warm Comptonthick corona (Zhang et al. 2000; Życki et al. 2001). The temperature, optical depth, and energetic output of the corona are dependent on distance from the black hole, accretion rate, and toroidal magnetic field strength. A high accretion rate and a strong toroidal field are favorable conditions for the development of a warm ∼5 keV and optically thick τ_{cor} ∼ 11 corona, as shown in the upper left corner of panel a in Fig. 11.
The magnetic field strength, β_{0}, seems to have a leading role in regulating the formation of the corona because the energy is stored very efficiently in the poloidal magnetic field near the equatorial plane, and the Poynting flux, F_{mag}, quickly increases, reaching its maximum value at approximately z/H ∼ 2 for the case presented in Fig. 4, panel d. At this point, the magnetic energy stored in the outflowing field is gradually released by heating the gas. Farther out in the atmosphere, this heating rate exceeds the magnetic energy gain from the toroidal field production (MRI). Asymptotically, the Poynting flux will fully convert into the radiation flux, but practically, this convergence may be very slow in some cases (particularly for strongly magnetized disks, where the magnetic pressure declines slowly with height). Therefore, the magnetic pressure gradient q (Eq. (34)) is an important quantity that influences the release of magnetic energy and in consequence, the formation of the corona (see panel b) in Figs. 7 and 8. A stronger toroidal field, that is, a lower β_{0}, implies a milder gradient and more accretion energy that is transported from the disk to the atmosphere by Poynting flux. This is nicely demonstrated in panel b of Fig. 11, where the values of β_{0} are lowest, and where the warm corona has the highest optical thickness (panel a) of the same figure. On the other hand, the radiative flux F_{rad} will never saturate asymptotically to a constant value for q ≤ 1 (see panels d) and i in Fig. 11.
However, an optically thick corona may be subject to local thermal instability caused by contribution of freefree opacity in the transition layer. These instabilities have been problematic in radiative transfer computations of accretion disk atmospheres of AGN for decades (Krolik & Kriss 1995; Nayakshin et al. 2000; Różańska et al. 2002, 2015). We were expecting that in case of hotter disks around black holes of stellar mass those instabilities would be absent, and in addition, magnetic heating might remove them. Our results show the opposite, however: a stronger field expands in the corona vertically and at the same time prevents formation of a very hot disk that is supported by radiation pressure, which might be sensitive to thermal instabilities. When this happens, two outcomes are possible: the sharp transition between cold and hot phase will be smoothed out by thermal conduction and form a relatively plain and regular atmosphere (Różańska 1999), or prominencelike clumps of dense matter will form, and will coexist at the same level with hot medium, similar to lower solar corona. The clumps could condense or simply emerge from the disk, pulled by magnetic field buoyancy (Fig. 1 in Jiang et al. 2014). Given the violent and turbulent nature of an assumed engine that powers the disk (MRI), the clumpy twophase medium interpretation seems more likely, but also poses a huge challenge for spectrum modeling. Variations of this scenario that often consider not only a corona but an entire disk that is clumpy have been discussed in the literature (Schnittman & Krolik 2010; Yang et al. 2015; Wu et al. 2016). Regardless of the interpretation, we show that thermal instability occurs in very many cases, even though GBHB disks are hot and dense and not as prone to it as AGN disks.
Although thermal instability affects the very same transition region that is responsible for producing optically thick corona, our model places constraints on the parameters for which a warm or hot corona can exist. For this purpose, we drew the radial extension of the corona base determined by the usual temperature minimum when no thermal instability is present, or by the temperature of the upper stable branch where thermal instability starts to developed (the cross point in the upper panels of Fig. 6) with the black dotted line in Fig. 12. The gas density, average temperature, and the radial profiles of the Compton parameter at selected optical depths for three accretion rates are shown in Fig. 12.
Fig. 12. Radial profiles of selected quantities: gas density ρ [cm^{−3}], average temperature T_{avg} [K], and average Compton y_{avg} (Eq. (46)), measured at different depths in the disk atmosphere. Values at τ = 0.1, τ = 1, and τ = 10 are shown by dashed blue lines, at temperature minimum (τ_{cor}) by the green line, at thermalization depth by the red line, and at the corona base (determined by temperature minimum or thermal instability in case it is developed) with the black dotted line. 

Open with DEXTER 
When the accretion rate is low, thermal instability arises at each distance from the black hole and only a hot corona appears, which is four orders of magnitude less dense than the disk atmosphere (three upper panels of Fig. 12). A corona like this would be optically thin, and the transition between disk and corona occurs at about z/H ∼ 20, as demonstrated in the upper panels of Fig. 6. For ṁ = 0.05, a warm stable corona can form up to R ∼ 10 R_{Schw}, and this corona is dense and optically thick, reaching at the base almost τ_{es} ∼ 10 (three middle panels of Fig. 12). Nevertheless, farther out from the black hole, thermal instability arises and a stable corona is again very hot and rare. Such a hot and rare outer corona dissipates only about 1% of the local energy and cannot be associated with the usual hard Xray corona observed in accreting sources. Finally, for a high accretion rate ṁ = 0.5, the radial extent of the warm corona is the largest ∼60 R_{Schw} and optical depth τ_{es} = 10. Thermal instability and a hot stable corona appear farther away from the black hole. Our results strongly indicate that a warm optically thick corona tends to develop in the innermost disk regions, while the outer disk is covered only by a rare hot corona. Furthermore, a warm corona is more radially extended for disks of high accretion rates. This feature explains why we see soft Xray excess only in the most luminous accreting sources.
Our model predicts a corona formation that is adjusted selfconsistently by model parameters, therefore it places constraints on several quantities that can be directly compared with observations. The fist straightforward outcome is that the value of the temperature T_{avg} is in the range from 0.1–10 keV. The second constrained parameter is the density of the warm corona, which arises from the dense disk atmosphere. The predicted density from our model is about 10^{16 − 20} cm^{−3}. This density is very high, but it agrees with the density of the ionized outflow in accreting black holes that was determined from observations by Miller et al. (2014). Our model supports the scenario that eventually, a disk wind can arise from the upper dense layers of a magnetically driven accretion disk atmosphere.
The most important observable is the ratio of the thermal energy that is released in the corona f in comparison to the total energy that is generated in the system. This ratio is always higher than zero because we computed a Comptoncooled atmosphere, which is characterized by typical temperature inversion with increasing optical depth (Różańska et al. 2011). In the frame of our model, we showed that up to half of the total accretion energy can be dissipated in the corona in case of GBHBs (up to f = 0.75 for the highest value of α_{B}), therefore we are not able to reproduce the passivedisk scenario discussed by Petrucci et al. (2013, 2018) for the case of AGN. This is a natural consequence of our assumed heating mechanism, which can transport a large portion of energy to the corona, but requires very effective magnetic coupling with gas down to the equatorial plane. The value of the Compton parameter averaged over the optical depth of the warm corona we determined with our model is consistent with the values determined from Xray observations of several GBHBs: y ≈ 0.3 (Życki et al. 2001), y ≈ 0.4 (Goad et al. 2006), up to y ∼ 1 (Zhang et al. 2000).
One of the main postulates of the BAR15 model is that it can explain the spectral state transition of the Xray binary that is observed in many sources (e.g., Yamada et al. 2013) by a smooth transition from the soft to hard state by introducing an intermediate dead zone between the disk and the coronal accretion flow. They used the condition for MRI shutoff when the gas energy is not high enough to sustain the dynamo (Pessah & Psaltis 2005). In our model, we considered only one zone where MRI is always operating. We suggest here that the spectral state transition should occur depending on the strength and radial extension of the warm optically thick corona rather than the full MRI shutoff. In the frame of our model we did not consider the transition from optically thick disk to inner advectiondominated accretion flow.
As discussed in BAR15 and SSAB16, one of the key features of MSDs is that they require a strong global poloidal magnetic field for the magnetic dynamo to operate and generate the toroidal field that supports the disk. In case of our typical model α_{B} = 0.1, the toroidal magnetic field in the equatorial plane and at the base of corona typically is about 10^{8} G in the inner disk region and of 10^{5} G for radii above 70 R_{Schw}. These values are lower than those detected in pulsars, but they fully agree with the prediction for the poloidal magnetic filed value that is required for the jetlaunching mechanism.
Studies suggest that a poloidal magnetic fields is important in jetlaunching mechanisms (see also Beckwith et al. 2009; McKinney et al. 2012; Liska et al. 2018). Zhang et al. (2000) noted that in both objects where an optically thick corona was observed, jets were also present. Our study shows that the most prominent warm corona will occur in a state when the dynamo is very effective (high α_{B}) and the contribution of the magnetic field to the disk structure is high (low β_{0}). On the other hand, no jet was observed in Mrk 509 (Petrucci et al. 2013). Nevertheless, the case of AGN will be presented in a forthcoming paper.
Acknowledgments
We thank Bożena Czerny for all comments and helpful discussion. This research was supported by Polish National Science Center grants No. 2015/17/B/ST9/03422 and 2015/18/M/ST9/00541.
References
 Anderson, E., Bai, Z., Dongarra, J., et al. 1990, in Proceedings of the 1990 ACM/IEEE Conference on Supercomputing, Supercomputing ’90 (Los Alamitos, CA, USA: IEEE Computer Society Press), 2 [Google Scholar]
 Bai, X.N., & Stone, J. M. 2013, ApJ, 767, 30 [NASA ADS] [CrossRef] [Google Scholar]
 Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 [NASA ADS] [CrossRef] [Google Scholar]
 Ballantyne, D. R., Ross, R. R., & Fabian, A. C. 2001, MNRAS, 327, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Beckwith, K., Hawley, J. F., & Krolik, J. H. 2009, ApJ, 707, 428 [NASA ADS] [CrossRef] [Google Scholar]
 Begelman, M. C., & Pringle, J. E. 2007, MNRAS, 375, 1070 [NASA ADS] [CrossRef] [Google Scholar]
 Begelman, M. C., Armitage, P. J., & Reynolds, C. S. 2015, ApJ, 809, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Di Salvo, T., Done, C., Życki, P. T., Burderi, L., & Robba, N. R. 2001, ApJ, 547, 1024 [NASA ADS] [CrossRef] [Google Scholar]
 Done, C., Davis, S. W., Jin, C., Blaes, O., & Ward, M. 2012, MNRAS, 420, 1848 [NASA ADS] [CrossRef] [Google Scholar]
 García, J. A., Kara, E., Walton, D., et al. 2019, ApJ, 871, 88 [NASA ADS] [CrossRef] [Google Scholar]
 Gierliński, M., & Done, C. 2004, MNRAS, 349, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Gierliński, M., Zdziarski, A. A., Poutanen, J., et al. 1999, MNRAS, 309, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Gladstone, J. C., Roberts, T. P., & Done, C. 2009, MNRAS, 397, 1836 [NASA ADS] [CrossRef] [Google Scholar]
 Goad, M. R., Roberts, T. P., Reeves, J. N., & Uttley, P. 2006, MNRAS, 365, 191 [NASA ADS] [CrossRef] [Google Scholar]
 Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Hawley, J. F. 2001, ApJ, 554, 534 [NASA ADS] [CrossRef] [Google Scholar]
 Hawley, J. F., & Balbus, S. A. 1991, ApJ, 376, 223 [NASA ADS] [CrossRef] [Google Scholar]
 Henyey, L. G., Forbes, J. E., & Gould, N. L. 1964, ApJ, 139, 306 [NASA ADS] [CrossRef] [Google Scholar]
 Hirose, S., Krolik, J. H., & Stone, J. M. 2006, ApJ, 640, 901 [NASA ADS] [CrossRef] [Google Scholar]
 Hirose, S., Krolik, J. H., & Blaes, O. 2009, ApJ, 691, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Jiang, Y.F., Stone, J. M., & Davis, S. W. 2014, ApJ, 784, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Jin, C., Ward, M., Done, C., & Gelbord, J. 2012, MNRAS, 420, 1825 [NASA ADS] [CrossRef] [Google Scholar]
 Keek, L., & Ballantyne, D. R. 2016, MNRAS, 456, 2722 [NASA ADS] [CrossRef] [Google Scholar]
 Krolik, J. H., & Kriss, G. A. 1995, ApJ, 447, 512 [NASA ADS] [CrossRef] [Google Scholar]
 Krolik, J. H., McKee, C. F., & Tarter, C. B. 1981, ApJ, 249, 422 [NASA ADS] [CrossRef] [Google Scholar]
 Kusunose, M., & Mineshige, S. 1994, ApJ, 423, 600 [NASA ADS] [CrossRef] [Google Scholar]
 Laor, A., Fiore, F., Elvis, M., Wilkes, B. J., & McDowell, J. C. 1994, ApJ, 435, 611 [NASA ADS] [CrossRef] [Google Scholar]
 Laor, A., Fiore, F., Elvis, M., Wilkes, B. J., & McDowell, J. C. 1997, ApJ, 477, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Liska, M. T. P., Tchekhovskoy, A., & Quataert, E. 2018, ArXiv eprints [arXiv:1809.04608] [Google Scholar]
 Madau, P. 1988, ApJ, 327, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Madej, J., & Różańska, A. 2004, MNRAS, 347, 1266 [NASA ADS] [CrossRef] [Google Scholar]
 Magdziarz, P., Blaes, O. M., Zdziarski, A. A., Johnson, W. N., & Smith, D. A. 1998, MNRAS, 301, 179 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083 [NASA ADS] [CrossRef] [Google Scholar]
 Mehdipour, M., BranduardiRaymont, G., Kaastra, J. S., et al. 2011, A&A, 534, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meurer, A., Smith, C. P., Paprocki, M., et al. 2017, PeerJ Comput. Sci., 3, e103 [CrossRef] [Google Scholar]
 Miller, J. M., Raymond, J., Kallman, T. R., et al. 2014, ApJ, 788, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Nakamura, K., & Osaki, Y. 1993, PASJ, 45, 775 [NASA ADS] [Google Scholar]
 Nayakshin, S., & Kallman, T. R. 2001, ApJ, 546, 406 [NASA ADS] [CrossRef] [Google Scholar]
 Nayakshin, S., Kazanas, D., & Kallman, T. R. 2000, ApJ, 537, 833 [NASA ADS] [CrossRef] [Google Scholar]
 Noble, S. C., Krolik, J. H., Schnittman, J. D., & Hawley, J. F. 2011, ApJ, 743, 115 [NASA ADS] [CrossRef] [Google Scholar]
 Ohsuga, K., & Mineshige, S. 2011, ApJ, 736, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Penna, R. F., Kulkarni, A., & Narayan, R. 2013, A&A, 559, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pessah, M. E., & Psaltis, D. 2005, ApJ, 628, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Petrucci, P.O., Paltani, S., Malzac, J., et al. 2013, A&A, 549, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Petrucci, P. O., Ursini, F., De Rosa, A., et al. 2018, A&A, 611, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Piconcelli, E., JimenezBailón, E., Guainazzi, M., et al. 2005, A&A, 432, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pounds, K. A., Stanger, V. J., Turner, T. J., King, A. R., & Czerny, B. 1987, MNRAS, 224, 443 [NASA ADS] [CrossRef] [Google Scholar]
 Proga, D., & Waters, T. 2015, ApJ, 804, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Różańska, A. 1999, MNRAS, 308, 751 [NASA ADS] [CrossRef] [Google Scholar]
 Różańska, A., & Czerny, B. 1996, Acta Astron., 46, 233 [NASA ADS] [Google Scholar]
 Różańska, A., Czerny, B., Życki, P. T., & Pojmański, G. 1999, MNRAS, 305, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Różańska, A., Dumont, A.M., Czerny, B., & Collin, S. 2002, MNRAS, 332, 799 [NASA ADS] [CrossRef] [Google Scholar]
 Różańska, A., Madej, J., Konorski, P., & Sądowski, A. 2011, A&A, 527, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Różańska, A., Malzac, J., Belmont, R., Czerny, B., & Petrucci, P.O. 2015, A&A, 580, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rybicki, G., & Lightman, A. 2008, Radiative Processes in Astrophysics, Physics Textbook (New York: Wiley) [Google Scholar]
 Salvesen, G., Simon, J. B., Armitage, P. J., & Begelman, M. C. 2016, MNRAS, 457, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Sądowski, A., Wielgus, M., Narayan, R., et al. 2017, MNRAS, 466, 705 [NASA ADS] [CrossRef] [Google Scholar]
 Schnittman, J. D., & Krolik, J. H. 2010, ApJ, 712, 908 [NASA ADS] [CrossRef] [Google Scholar]
 Schnittman, J. D., Krolik, J. H., & Noble, S. C. 2016, ApJ, 819, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Seifina, E., Titarchuk, L., & Shaposhnikov, N. 2016, ApJ, 821, 23 [NASA ADS] [CrossRef] [Google Scholar]
 Seifina, E., Chekhtman, A., & Titarchuk, L. 2018, A&A, 613, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
 Shaposhnikov, N., & Titarchuk, L. 2006, ApJ, 643, 1098 [NASA ADS] [CrossRef] [Google Scholar]
 Stobbart, A.M., Roberts, T. P., & Wilms, J. 2006, MNRAS, 368, 397 [NASA ADS] [CrossRef] [Google Scholar]
 Stone, J. M., Gardiner, T. A., Teuben, P., Hawley, J. F., & Simon, J. B. 2008, ApJS, 178, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Svensson, R., & Zdziarski, A. A. 1994, ApJ, 436, 599 [NASA ADS] [CrossRef] [Google Scholar]
 Tange, O. 2011, The USENIX Magazine, 36, 42 [Google Scholar]
 Tout, C. A., & Pringle, J. E. 1992, MNRAS, 259, 604 [NASA ADS] [Google Scholar]
 Turner, N. J. 2004, ApJ, 605, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, M.C., Xie, F.G., Yuan, Y.F., & Gan, Z. 2016, MNRAS, 459, 1543 [NASA ADS] [CrossRef] [Google Scholar]
 Yamada, S., Torii, S., Mineshige, S., et al. 2013, ApJ, 767, L35 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, Q.X., Xie, F.G., Yuan, F., et al. 2015, MNRAS, 447, 1692 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, S. N., Cui, W., Chen, W., et al. 2000, Science, 287, 1239 [NASA ADS] [CrossRef] [Google Scholar]
 Życki, P. T., Done, C., & Smith, D. A. 2001, MNRAS, 326, 1367 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Numerical procedure
The solution of the set of five equations with four boundary conditions described in the previous sections allows determining the density ρ, gas temperature T_{gas}, radiation pressure P_{rad}, radiative flux F_{rad}, and magnetic pressure P_{mag} at each point of the vertical profile. The requirement of satisfying four boundary conditions, two on each side of the computational interval, requires multiple iterations when shooting methods such as a RungeKutta scheme are used. Moreover, the equations are highly nonlinear, and keeping the integration error low requires a high number of grid points (up to 10^{5}), even using a fourthorder RungeKutta scheme. Although the computation time using a shooting method is still less than a few minutes, it becomes a problem when the global behavior over the parameter space is investigated, which requires computing thousands of models for different disk and magnetic parameters. Therefore, to solve this issue, we developed a variant of the relaxation method that was originally described by Henyey et al. (1964) for stellar evolution problems.
The relaxation methods are iterative, when iterations are made for the functions defined in the whole argument space. In the first step, some approximate solution for the functions, which might be close to the exact solution, is needed. Then by applying the procedure described below, the correction to the entire profile is obtained and added to the original solution. Because we have many strongly nonlinear terms in our equations while the method is intrinsically linear, we applied a reduction factor for these corrections in the first few iterations to avoid overshooting. The new solution was then used as input to the next step and the iteration continued until the mean relative value of the obtained correction fell below 10^{−7}, which we consider a convergence criterion.
To apply the relaxation method to our problem, define a grid of N points, in each we have five physical quantities. We introduce the following notation: the quantity j = 1…5 in point i = 1…N is denoted . We can also represent it as a vector with five components:
Then we discretize the set of all five equations together with the boundary conditions. For nondifferential equations (balance equation and boundary conditions), we evaluate them at the respective grid points, that is, . For all four differential equations we replace quantities with midpoint values and their spatial derivatives with finite differences. As they are evaluated at midpoints, there are only N − 1 terms,
The core of the problem is to perturb the entire profile to obtain a new one that is closer to the solution. Using the chain rule, we can express the computed profile in terms of perturbations to actual physical quantities:
where i = 1…N − 1 and k = 1…5. Most of the terms in this sum are zero because equations evaluated at point i depend only on quantities in points i and i + 1, thus l ∈ {i, i + 1}. The goal of the relaxation method is to transform all equations where , which means that each of the equations k = 1…5 is satisfied in each point i = 1…N. However, having started from some initial guess, this is more likely not true, and the remainder .
Analogically, we do the same for boundary conditions. When we denote the remaining lower boundary conditions as L and upper boundary conditions as U, the expansions in terms of perturbations with respect to variables yield
At this point we have a mathematically welldefined problem: we must solve Eqs. (A.4)–(A.6) to obtain the profile corrections . We can express it as the perturbation vector δX, remainder vector R and Jacobian matrix M:
We arrange our five variables and five equations into a vector using the following layout:
Now the remaining issue is to fill the matrix accordingly. For linear equations, to evaluate all nonzero matrix terms for i = 1…N, j = 1…5, k = 5 and l = i, we directly use the formula
For differential equations, each point depends on two consecutive points, and we evaluate these terms for i = 1…N − 1, j = 1…5, k = 1…4 and l ∈ {i, i + 1} using the chain rule:
Finally, for boundary conditions, we evaluate the matrix coefficients at lower (L) and upper (U) boundary as
The relaxation method provides excellent precision even for very few computational points because the method adjusts the entire profile at once, allowing us to investigate even marginally stable cases. The method is also clearly expandable for higherorder derivatives.
Its major disadvantage is the requirement that the partial derivatives of the lefthand sides of all equations (and boundary conditions) need to be evaluated with respect to all variables and their spatial derivatives. This makes the code very difficult to modify and error prone. We solved this problem using a symbolic computation techniques within the sympy framework (Meurer et al. 2017), which also allowed us to generate readytocompile procedures to evaluate the matrix coefficients that we included in our code, which is written in Fortran 2008. We then solved the system of linear equations using the DGBSV procedure from LAPACK (Anderson et al. 1990) and obtained corrections to all five variables at each point. This method is very fast and allows us to obtain a precise solution in less than half a second on a typical desktop computer.
All Figures
Fig. 1. Schematic illustration of a slice through the disk plane (only one side is shown due to symmetry). The black hole is shown as the circle on the left. The gray area represents the optically thick and geometrically thin disk where most of the accreting mass is located. The warm corona that covers the disk, in which the magnetic energy is released as radiation, is marked in orange. Magnetic and radiative energy flux are represented by blue and red arrows, respectively. The part of the corona that might collapse as a result of thermal instability is shown in magenta. The inner optically thin and geometrically thick hot flow that we did not consider is drawn in yellow for completion. 

Open with DEXTER  
In the text 
Fig. 2. Dependence of model parameters: toroidal field production efficiency, α_{B} (lower panel), reconnection efficiency, ν (middle panel), and the magnetic field gradient parameter, (upper panel) on the magnetic parameter in the disk midplane β_{0}. Points denote values obtained from numerical simulations in SSAB16 and lines are our results from the powerlaw relation between α_{B} and η (Eq. (33)) for four values of ξ. 

Open with DEXTER  
In the text 
Fig. 3. Toroidal magnetic field at the equatorial plane B_{0} and in the corona B_{cor} at optical depth τ_{cor} (see text for definition) in Gauss units. Logarithmic values of the magnetic field are plotted in the radius (R/R_{schw}) – accretion rate (ṁ) parameter plane. 

Open with DEXTER  
In the text 
Fig. 4. Vertical structure of the diskcorona system computed at the single radius R = 10 R_{Schw} for ṁ = 0.05. The horizontal axis is the height above the disk midplane given in scale height H defined by Eq. (36). Three significant locations are marked with vertical dotted lines: the photosphere (green), the temperature minimum (dark yellow), and the thermalization depth (red). The following panels display (a) radiation and gas temperature; (b) gas, radiation, magnetic, and total pressure; (c) gas density and pressure gradients; (d) radiative energy flux, magnetic (Poynting) flux, and gas heating rate all normalized to the total accretion flux F_{acc}; (e) net radiative cooling terms (Λ_{B} and Λ_{C}) and gas magnetic heating rates; and (f) electron scattering and total and effective optical depths. 

Open with DEXTER  
In the text 
Fig. 5. Logarithm of the net cooling rate Λ_{rad} (Eq. (17)), plotted in the n_{H} – T parameter plane ranging from blue (low rate) to red (high rate). The points of constant pressure are connected with black dashdotted lines. The red dotted contours show the contribution of the Compton cooling to the total cooling. The solid white line shows where the cooling function is constant in isobaric perturbation, and separates the area of unstable solutions where the gradient of the cooling function ℒ_{rad} < 0. Dark blue lines of increasing width represent the vertical structure for three accretion rates 0.035, 0.050, and 0.072. 

Open with DEXTER  
In the text 
Fig. 6. Onset of thermal instability in a lowtemperature disk, ṁ = 0.035 shown in the top panels, and its stabilization with the increase of the temperature, ṁ = 0.05 and 0.072 shown in the middle and bottom panels. Other disk parameters are as given in Sect. 3.1. In each panel the vertical axis represents the possible gas temperature T assuming constant gas pressure (δP_{gas} = 0), and the horizontal axis is the height above the disk midplane given in scale height H (Eq. (36)). The possible solutions that satisfy the radiative energy balance are shown with thick black contours, while the solution obtained from our numerical computations is additionally marked with light gray dots. A thinner gray line is used to show the radiation temperature, T_{rad}. Color maps in the first panels column show the contribution from Compton scattering to the total cooling. In the second column, the overall heatingcooling balance is shown: the yelloworange areas are where the heating dominates cooling, and blue areas show where the cooling rate exceeds the heating. The third column presents the stability diagnostic, i.e., the cooling function gradient ℒ_{rad} (Eq. (37)), with green areas marking the stable regions and violet areas the unstable regions. Two limit solutions according to Eqs. (40) and (41) are plotted with orange lines (dashdotted and dashed, respectively) and show excellent agreement with the numerical solution. Additionally, the content of neighboring plots (from right panels to the middle, and from middle to the left) are superimposed as thin contours for the sake of clarity. 

Open with DEXTER  
In the text 
Fig. 7. Vertical structure of the diskcorona models computed at R = 10 R_{Schw} for different values of the accretion rate (ṁ), plotted against the total optical depth τ. The panels show (a) the gas temperature, (b) the gradient of the magnetic pressure (Eq. (34)), (c) the radiative energy flux normalized to the total accretion flux F_{acc}, (d) the cooling function gradient (Eq. (37)), (e) the contribution of Compton cooling to the total cooling, and (f) the ionization parameter Ξ (Eq. (43)). 

Open with DEXTER  
In the text 
Fig. 8. Same as in Fig. 7, but for different values of the magnetic parameter α_{B}. 

Open with DEXTER  
In the text 
Fig. 9. Properties of the corona on the radius (R/R_{schw}) – accretion rate (ṁ) parameter plane. All models are calculated for M_{BH} = 10M_{⊙}, α_{B} = 0.1, and ξ = 0.05. The panels display (a) the average temperature of the corona (Eq. (44)) in keV (colors) and τ_{cor} (contours); (b) the ionization parameter Ξ (colors) and magnetic β parameter (contours), both at the disk midplane; (c) the number density at τ = 2/3 (colors) and the total column density Σ in g cm^{2} (contours); (d) the fraction of radiative energy produced by the corona f (Eq. (45), colors), the maximum value of the magnetic field gradient q (Eq. (35), contours); (e) the y parameter of the warm corona (Eq. (46)), computed up to the thermalization layer τ^{⋆} = 1 (colors) and to τ_{cor} (contours); the (f) photon destruction probability (colors) and electron scattering optical depth (τ_{es}) (contours) at the thermalization layer; (g) the minimum value of ℒ_{rad} (Eq. (37)) throughout the disk height (colors), here, green areas are stable disks, and magenta areas indicate thermal instability in the corona, the gray dotted line corresponds to ℒ_{rad} = 0, and the ionization parameter Ξ at the corona base (contours); and (h) the radiative energy flux normalized to the total accretion flux (colors), and the disk scale ratio H/R (contours). 

Open with DEXTER  
In the text 
Fig. 10. Same as in Fig. 9, but on the radius (R/R_{schw}) – toroidal field production parameter (α_{B}) plane. 

Open with DEXTER  
In the text 
Fig. 11. Same as Fig. 9, but on toroidal field production parameter (α_{B}) – accretion rate (ṁ) plane. 

Open with DEXTER  
In the text 
Fig. 12. Radial profiles of selected quantities: gas density ρ [cm^{−3}], average temperature T_{avg} [K], and average Compton y_{avg} (Eq. (46)), measured at different depths in the disk atmosphere. Values at τ = 0.1, τ = 1, and τ = 10 are shown by dashed blue lines, at temperature minimum (τ_{cor}) by the green line, at thermalization depth by the red line, and at the corona base (determined by temperature minimum or thermal instability in case it is developed) with the black dotted line. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.