Free Access
Issue
A&A
Volume 614, June 2018
Article Number A79
Number of page(s) 12
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201732345
Published online 18 June 2018

© ESO, 2018

1 Introduction

Galactic black hole (BH) X-ray binaries (XRBs) are important cosmic laboratories which allow us to study many aspects of accretion. Thanks to their relative proximity, their X-ray brightness is often high, and thus their spectra and timing can be studied in detail. Given their size, their variability can be observed on a wide range of timescales, from the light travel time across the BH horizon to years. This allows us to observe their outbursts and state transitions between the two main spectral states, hard and soft.

Most of the studies agree on the geometry in the soft state. The accretion flow consists of an optically thick and geometrically thin accretion disc (Shakura & Sunyaev 1973; Novikov & Thorne 1973) extending down to the innermost stable circular orbit (ISCO). The blackbody-like spectrum of the disc is followed by a hard X-ray tail, which can have varying amplitude. It is probably due to coronal emission by hybrid, i.e. both thermal and non-thermal, electrons, which Compton scatter disc blackbody photons (e.g. Poutanen & Coppi 1998; Gierliński et al. 1999; Poutanen & Vurm 2009).

The hard-state data are well described by Comptonization proceeding in a hot medium. Location of this medium, above the accretion disc or within its truncation radius, has been a matter of debates for decades. Related to this, there is an intense controversy regarding the value of the truncation radius. In transients, there is no inner disc during quiescence (Lasota et al. 1996; Dubus et al. 2001), so there has to be a truncated disc during initial parts of the outbursts. For many years, the dominant thinking was that the disc reaches the ISCO during the transition to the soft state. Lately, a number of publications have claimed (see the references below) that the disc is already close to the ISCO in the hard state, forming a sandwich-corona geometry. If this is the case, there are two questions that remain unanswered: at what Eddington ratio (between the quiescence and the luminous hard state) does the disc reach the ISCO, and what are observational signatures of this phenomenon?

Sandwich-corona models (see Fig. 1a) were proposed in Haardt & Maraschi (1991, 1993) in application to active galactic nuclei. These papers and many subsequent ones studied the disc-corona coupling in the limiting case of a passive disc with all the dissipation taking place in a corona above it (e.g. Haardt & Maraschi 1991, 1993; Stern et al. 1995; Poutanen & Svensson 1996; Zdziarski et al. 1998), and some treated the fraction of the dissipation in the corona above the disc as a free parameter (e.g. Svensson & Zdziarski 1994). Models of self-consistent coronal emission and reflection were considered in Stern et al. (1995), Poutanen & Svensson (1996), Svensson (1996), Malzac et al. (2001), and Malzac et al. (2005). These models take into account up-scattering of reflected photons, but the intrinsic (viscous) dissipation in the disc was assumed to be zero.

A lot of attention was paid to the energy balance between the disc and the corona. The emitted spectra were calculated self-consistently with the electron temperature in the corona, accounting for the feedback from the disc. The radiation intercepted by the disc is partly reflected by Compton scattering in the upper layers, and is partly reprocessed and emitted at lower energies. The calculation of reflection in most of these models was done assuming neutral gas and neglecting line emission (the approach is described in Magdziarz & Zdziarski 1995 and Poutanen et al. 1996), and the reprocessed emission was approximated with a diluted blackbody. These works showed that the models of the corona sandwiching the disc (slab-corona models) produce soft spectra, with photon index Γ ≳ 2.1 (see Fig. 2 in Malzac et al. 2001). The main reason for this is the large flux of the reprocessed photons which cool electrons in the corona. Specifically, since the corona emits about a half of its power towards the disc, there is an approximate equipartition between the power of the seed photons irradiating the corona and the power dissipated in the corona, which then leads to spectra roughly flat in νFν, i.e. with Γ ≃ 2. Adding intrinsic dissipation then leads to even softer spectra (e.g. Svensson 1996). This is in contradiction to the observed hard-state spectra, with Γ ~ 1.7. Specifically, in the case of a hard-state spectrum of GX 339–4, Zdziarski et al. (1998) found that the covering factor of the cold medium as seen from the hot plasma has to be ≲0.3 (see their Fig. 6), which rules out the geometry of a corona above the disc.

A more advanced self-consistent slab-corona model was considered in Malzac et al. (2005), who computed reprocessing in the disc in thermal and ionization equilibrium accounting for radiative transfer in lines using a two-stream approximation. They considered two models for the disc, one at a constant density and one at a constant pressure of the cold phase, showing that in the first case the ionization parameter log ξ has to be at least 3.5 for the Comptonization spectrum to be consistent with the spectra observed in the hard state. On the other hand, they also showed that the Comptonization spectra for the second more physically realistic model of the constant pressure disc depend very weakly on the ionization parameter and cannot be reconciled with the data even for very high ionization. Regardless of the specific set-up, all the proposed models in this category take into account the re-emission of the coronal radiation irradiating the disc as well as scattering of the reprocessed radiation in the corona.

Harder spectra can be obtained in the patchy corona models (Fig. 1b), which can represent the magnetic flares above the accretion disc (Stern et al. 1995; Svensson 1996). In this case, most of the reprocessed photons directly escape towards the observer and are not intercepted by the coronal plasma, leading to the higher electron temperature and, subsequently, hard spectra. However, the luminosity of the soft reprocessed photons is similar to the luminosity escaping in the form of hard X-rays, again inconsistent with the hard-state data. This conclusion is drawn assuming zero flux from the viscous dissipation in the cold disc, which should be added to the soft reprocessed flux, increasing the soft flux even further. Another possibility is to assume that the corona is outflowing with mildly relativistic velocity (Figs. 1c and d; Beloborodov 1999; Malzac et al. 2001). Because of the relativistic aberration, the coronal flux is beamed away from the disc, which reduces the reprocessed and the reflected components, making this model a feasible alternative for the hard-state spectra.

An alternative scenario considers the hot medium in the form of the hot flow surrounded by the cold disc truncated at large radius (rin ≳ 10Rg, where RgGMc2 is the gravitational radius, see Fig. 1e, Esin et al. 1997, 1998; Poutanen et al. 1997; Dove et al. 1997). However, unlike the sandwich corona, this geometry has the problem of seed photon starving, as only a small fraction of soft disc photons is intercepted by the hot flow. For example, if the disc truncation radius is rin = 100Rg, only about 0.1% of the disc photons enter the zone of major energy release (around 10Rg) as the result of reduction of the solid angle and Lambert law. Three solutions for this problem were proposed. In the first the disc overlaps with the hot medium (Fig. 1f; Poutanen et al. 1997). The second option proposes additional soft photons from the cold medium in the form of clumps within the hot flow (Celotti et al. 1992; Zdziarski et al. 1998) or a residual accretion disc (Liu et al. 2007), see Fig. 1g. The third possibility is to consider cyclo-synchrotron radiation from the hot flow itself as the additional or even dominant source of seed photons (Fig. 1h, Esin et al. 1997; Wardziński & Zdziarski 2001; Poutanen & Vurm 2009; Malzac & Belmont 2009; Veledina et al. 2013).

Recently, the truncated disc/hot flow geometry has been questioned in a number of works. Papers claiming the accretion disc in the hard state of BH binaries extends to the ISCO, or almost to it, fall into two categories. In one category, the hard X-rays are supposed to be emitted by a very compact source located on the rotation axis of the BH (e.g. Parker et al. 2015; Fürst et al. 2015; Miller et al. 2015; Walton et al. 2017). The way the compact source, known as the lamppost, emits, is not specified. If its emission is due to Compton up-scattering of blackbody photons emitted by the disc, the size of the lamppost has to be relatively large in order to intercept enough photons (Dovčiak & Done 2016), which is in conflict with most of the fitting results of this model. Also, the proximity of the lamppost to the horizon, found in many fitted models, gives rise to the problem of photon trapping of most of the emitted photons by the BH, and runaway e± pair production within the lamppost (Niedźwiecki et al. 2016), which reduces the electron temperature to low values (~20 keV) inconsistent with the data.

A different argument is put forward by Steiner et al. (2017), who present a model for the accretion flow in the hard state attempting to self-consistently account for mutual interactions of the disc and the corona. In particular, they quantitatively take into account and stress the importance of Compton scattering by the corona of photons reflected from the disc. In order to compute the reprocessed emission that comes in the form of lines and recombination continua, they considered a constant density model of the disc. They also linked the emission due to viscous dissipation within a disc truncated at an inner radius, rin, to the coronal emission by taking into account conservation of photons. They concluded that for the hard state a large truncation radius, rinRg ≫ 1, would require very high values, even super-Eddington values, of the accretion rate. This would argue strongly against the correctness of hard-state fit results obtaining a large rin in luminous hard states, in particular those for the low-mass X-ray binary GX 339–4 (Basak & Zdziarski 2016). We note that the Steiner et al. (2017) approach relies on the assumption that most of the seed photons come from the viscous dissipation in the accretion disc, while the amount of soft photons coming from reprocessing was heavily underestimated, and other sources of seed photons in the hard state, for example from the synchrotron emission, were ignored (see Sect. 3.2 for details). A major shortcoming of this work is that the cooling of the electrons in the corona by the reprocessed emission coming from the disc and the resulting coronal spectrum were not computed self-consistently.

In this paper, we reconsider the issue of energy balance between the disc and the corona. We calculate the self-consistent Comptonization spectra of the hot medium in various geometries and convolve them with the most recent reflection model xillver (García & Kallman 2010; García et al. 2013). We show that the sandwich-corona geometry is not a viable alternative for the hard state, even accounting for the Comptonization of the ionized reflection emission, the main reason being the substantial soft photon flux from the reprocessed emission, capable of cooling the coronal medium and leading to too soft spectra. We compare our results to previous works involving the geometry of sandwiching corona and revisit the arguments against the truncated disc geometry. We discuss different accretion flow geometries in the context of recent multiwavelength discoveries and conclude that the truncated disc/hot flow geometry with cyclo-synchrotron seed photons is the only scenario capable of explaining the whole complex of the hard-state data.

thumbnail Fig. 1

Geometries of the central parts of the accretion flow proposed over time. (a) Static sandwiching corona, (b) static patchy corona, (c) outflowing sandwiching corona, (d) outflowing patchy corona, (e) cold accretion disc detached from the hot flow and (f) intersecting with it, (g) hot flow with cold clumps, (h) truncated disc and hot flow with substantial cyclo-synchrotron radiation.

2 Self-consistent modelling of the hot medium in various geometries

2.1 Sandwich corona

We discuss here the simplest model of the slab corona above a cold disc. The exact solution for this problem was first obtained by Stern et al. (1995), Poutanen & Svensson (1996), and Malzac et al. (2001) assuming reflection from neutral material (Magdziarz & Zdziarski 1995; Poutanen et al. 1996). Here we use the newer reflection model, xillverCp1 (García & Kallman 2010; García et al. 2013; Dauser et al. 2016) in the form of precomputed tables.2 In this model, the reflection spectra from a constant density slab are computed for the intrinsic Comptonization continuum described by the nthcomp model (Zdziarski et al. 1996) incorporated to the xspec package (Arnaud 1996). The main parameters of nthcomp are the photon spectral index, Γ, and the electron temperature, kTe. The blackbody seed photons are assumed to have a temperature kTseed of 50 eV. The reflection spectrum depends on the ionization state of the material, described by the parameter ξ = 4πFXn, where FX is the integrated flux over the incident spectrum in the range 1–103 Ry and n (cm−3) is the gas density (particle concentration). It also depends on the metal abundance, which we assume here to be solar. The normalization of the reflection spectrum is given in the Appendix of Dauser et al. (2016). The reflection spectrum contains the continuum reflected by coherent electron scattering as well as lines and recombination continua. Because the Klein–Nishina effect is ignored in the calculation of the reflection, the shape of the cut-off of the reflected spectrum energies above 50 keV is not correct (compare spectra from Fig. 2 in García et al. 2013 with the spectra presented in Magdziarz & Zdziarski 1995; Poutanen et al. 1996; see also Fig. 1 in Niedźwiecki et al. 2016). This may introduce an error of the order of 10–20% to the reflection albedo, a, and the reprocessing efficiency. The reflection spectra in the xillverCp model extend down to only 70 eV. This sometimes results in a significant difference between the angle-integrated reflected and illuminating fluxes which is inconsistent with the energy conservation. We solved this problem by adding a soft component with the blackbody-like spectrum of temperature 50 eV to the reflected spectrum. Thus, the emission from the cold disc consists of two components: the reflection of the coronal X-rays and an additional soft component (see red dashed curve in Fig. 2).

Electrons in the corona Comptonize disc emission into a power-law-like spectrum, which turns over at energies slightly above kTe (Stern et al. 1995). The main parameter defining the Comptonization spectrum is the Thomson optical depth, τT, of the corona. The coronal electron temperature is not a model parameter, but is computed self-consistently by solving for the energy balance (Poutanen & Svensson 1996). The angle-dependent Comptonization spectrum in a slab geometry is computed using the code described in Poutanen& Svensson (1996), which utilizes an exact Compton scattering redistribution function. The code is similar to the compps model used in the xspec package, the only difference being that reflected photons are now scattered in the corona. To evaluate the reflection spectrum, we first compute the angle-averaged Comptonized spectrum of radiation going from the corona down towards the disc. This spectrum is then fit in the 2–10 keV range by a power law to get the photon index Γ, which together with the current value of kTe is used to obtain the angle-dependent reflection spectrum by interpolation from the precomputed tables. We note that the angle-averaged Comptonization spectra computed using nthcomp and compps are very close for kTe ≲ 100 keV (Zdziarski et al. 1996, 2003). This justifies using the reflection tables for nthcomp. We then compute the reflected flux in the 0.07–103 keV range (where the xillverCp model is computed). The remaining fraction of the bolometric Comptonized downward-directed flux is assumed to be thermalized in the disc and emitted as a blackbody of the temperature of 50 eV. This guarantees the correct energy balance in the disc. The self-consistent solution for the corona, which requires that the flux emitted by the disc is equal to the Comptonizedflux directed towards the disc, is obtained by iterating the coronal temperature. Thus, we obtain the equilibrium kTe and the escaping spectrum.

Examples of the resulting spectra from the slab corona as observed at some typical inclination of i = 60° are shown in Fig. 2 for two values of log ξ = 1 and 3 and a fixed Thomson optical depth of τT = 0.4. This value of τT gives the equilibrium coronal temperature of kTe = 80–100 keV (see blue open squares in Fig. 3), which is typical for the BH XRBs in the hard state (Wardziński et al. 2002; Zdziarski & Gierliński 2004; Ibragimov et al. 2005; Torii et al. 2011; Burke et al. 2017). We see that in a weakly ionized case with log ξ < 2, the Comptonized spectrum is rather soft (see the red squares in Fig. 3) with Γ ≈ 2.2 (slightly depending on the inclination i), which is consistent with previous results obtained for neutral reflection (see e.g. the high compactness range of the slab model in Fig. 6 in Svensson 1996 or the hr = 0 case in Fig. 2 in Malzac et al. 2001). The reflection amplitude in the total spectrum is rather small here because only a fraction of the reflected flux, exp(−τT∕ cos i), is transmitted through the corona without scattering. At the intermediate ionization of log ξ = 3, the coronal spectrum becomes harder with Γ ≈ 2.04. In both casesthe spectrum is much softer than those observed in the hard state of BH XRBs with Γ ≈ 1.6–1.8. We also note that the Comptonized spectra computed using the xillverCp reflection model are indistinguishable from those computed using the simpler reflection model pexriv for log ξ < 2; only at high ionization will the presence of lines and recombination continua in the xillverCp model start to make a difference (see Sect. 3.1). However, the resulting coronal spectrum is soft with Γ > 2.

The situation changes at high ionizations with log ξ ≥ 3.5 (see Fig. 3). Now the disc is strongly ionized and reprocesses very little radiation; most of the X-rays are scattered back. For log ξ = 3.5, the corona becomes photon-starved resulting in a harder spectrum with Γ ≈ 1.75, with no visible spectral features seen in the xillverCp spectrum. At even higher ionization the disc acts like a mirror resulting in the very hard coronal spectrum. Thus, we find that the slab corona can produce spectra similar to the hard states of BH XRBs only for very high ionizations log ξ ≳ 3.5 (see also Fig. 5 in Malzac et al. 2005), the lower boundary coinciding with the upper limit of the observed ionization levels for sandwich corona fits 2.0 ≲ log ξ ≲ 3.4 (Parker et al. 2015; Basak & Zdziarski 2016; Steiner et al. 2017; Basak et al. 2017).

Higher coronal optical depth, τT = 0.8, results in even softer spectra (see red triangles in Fig. 3) and predicts very low coronal temperatures of below 50 keV (blue open triangles), not consistent with the data. The main reason for spectral softening is that at higher τT a larger fraction of disc radiation would be scattered back to the disc. For the same spectral slope this would mean that more radiation is striking the disc than is reflected, which breaks the energy balance. Thus, the corona has to adjust by cooling and producing a softer spectrum. For the same reason a lower optical depth, τT = 0.2, gives slightly harder spectra with Γ = 2.03 at low ξ. Harder spectra with Γ < 1.8 can be achieved at log ξ > 3. However, the coronal temperature in this case is in excess of 150 keV (blue open circles), which is well above the observed temperature values in most of the observations (Wardziński et al. 2002; Zdziarski & Gierliński 2004; Ibragimov et al. 2005; Torii et al. 2011; Burke et al. 2017).

The observed range of Γ (1.6–1.8) and kTe (50–120 keV) is shown as a shaded area in Fig. 3. We see that most of the models demonstrate either too soft spectra or too high electron temperatures. Hence, the slab-corona model with no intrinsic disc dissipation requires much fine-tuning in both optical depth, τT = 0.4 ± 0.1, and ionization parameter, log ξ > 3.5, to be consistent with the observed slope of the Comptonized continuum. Furthermore, in the hard state the best-fitting ionization parameter is usually much smaller (e.g. Plant et al. 2014; Parker et al. 2015; Basak & Zdziarski 2016; Steiner et al. 2017; Basak et al. 2017).

We note here that the Compton reflection is mirror-like only in the Thomson limit, i.e. for spectral cut-off below E ≈ 50 keV. In reality, the Klein–Nishina effect results in the absorption of a fraction of the energy of incident photons at E ≳ 50 keV (see Fig. 11 in Zdziarski et al. 2003). This effect, not treated in xillverCp, leads to additional cooling and softening of the observed spectrum. For example, in the case of kTe = 120 keV and Γ≃ 1.5 considered inZdziarski et al. (2003), the reflection albedo is still as low as 0.48, which means that 52% of the irradiating luminosity is re-emitted at low energies and serve as seed photons, cooling the plasma. Thus, the spectral indices computed using xillverCp reflection are the lower limits and the actual spectra should be softer. Furthermore, the approximation of the constant-density slab used by xillverCp predicts too few soft photons compared to a more physically realistic case of reflection from a constant pressure slab (Fig. 5 in Malzac et al. 2005) or a slab in hydrostatic equilibrium (Fig. 8 in Nayakshin et al. 2000). The latter models would thus produce softer Comptonization coronal spectra and would rule out the slab-corona model for the hard state completely (see also Nayakshin & Dove 2001).

thumbnail Fig. 2

Self-consistent spectra from a homogeneous slab corona of τT = 0.4 above a disc as observed at inclination i = 60°. The red and blue curves correspond to log ξ = 1 and 3, respectively. The dotted lines show the Comptonization continuum, the dashed lines show the reflection xillverCp spectrum together with reprocessed blackbody emission, and the solid lines give the total spectrum. The reflection spectra have been re-binned at a resolution Δ log E = 0.025. The bolometric luminosity of the reflection spectrum here is smaller than that of the Comptonization continuum because part of the reflected/reprocessed emission is Compton scattered in the corona.

thumbnail Fig. 3

Dependence of the 2–10 keV photon spectral index Γ of the Comptonized spectra from the sandwich corona (red filled symbols and left axis) and of the coronal temperature kTe (blue open symbols and right axis) on the ionization parameter for various coronal optical depths. Circles, squares, and triangles correspond to τT = 0.2, 0.4, and 0.8, respectively. The shaded area represents the typical observed range for Γ (1.6–1.8) and kTe (50–120 keV) (Wardziński et al. 2002; Zdziarski & Gierliński 2004; Ibragimov et al. 2005; Torii et al. 2011; Burke et al. 2017).

2.2 Hot flow with cold clumps or overlapping with the disc

To make the spectra harder than typically predicted by the slab corona, we need to somehow reduce the amount of soft photons coming from the disc. It is straightforward to assume that the disc does not reside under the hot medium all the way down to the black hole, but that it is truncated at some radius. Self-consistent models for the central hot corona–approximated as either a sphere of radius rc irradiating an outer disc with some viscous dissipation (or without) and truncated at some radius rin (allowed to overlap with the corona) or as a hot flow with embedded clouds–have been considered by Zdziarski et al. (1998, 1999), and see also Poutanen et al. (1997). In this model, the corona Compton-scatters blackbody and reflected photons emitted by either the outer disc or cold clouds, reaching the luminosity Lc=fscA(Γ,Tbb,Te)Lbb,\begin{equation*} L_{\rm{c}}=f_{\rm{sc}} A(\mathrm{\Gamma}, T_{\rm{bb}}, T_{\rm{e}}) L_{\rm{bb}},\end{equation*}(1)

where fsc is the fraction of the blackbody photons that irradiate the corona (which depends on the geometry) and A(Γ, Tbb, Te) is the Comptonization amplification factor (which depends on the spectral index, Γ; the blackbody temperature, Tbb; and the electron temperature, Te). The blackbody luminosity in turn originates from reprocessing of the coronal Compton-scattered photons irradiating the cold medium and the intrinsic viscous dissipation in the disc: Lbb=Lintr+Lirr,Lirr=RF(1a)Lc.\begin{equation*}L_{\rm{bb}}=L_{\rm{intr}}&#x002B;L_{\rm{irr}}, \quad L_{\rm{irr}} = R_{\rm{F}}(1-a)L_{\rm{c}}. \end{equation*}(2)

The reflection fraction RF follows from the assumed geometry. The albedo, a, depends on the incident spectrum and the ionization properties of the medium defined by the ionization parameter ξ. For a given geometry, the model determines the X-ray spectral index, Γ, and the strength of reflection, which are found to be correlated, as observed (see e.g. Zdziarski et al. 2003; Gilfanov 2010).

To illustrate the model, we consider a hot electron slab with the central plane filled by cold clouds (i.e. fsc =1 in this case) with the covering factor fcl, with fcl = 1 corresponding to the slab-corona case described in Sect. 2.1. As in the case of the sandwich-corona geometry, we inject blackbody and reflected photons from the cloud surface.

They are then Comptonized by the hot medium described by τT (now measured from the central plane to the outer surface) and the electron temperature, kTe. We compute the Comptonization angle-dependent spectrum and the reflection exactly as in Sect. 2.1 with compps and xillverCp models, respectively. Because fcl can now be below unity, only a fraction of the Comptonizing radiation crossing the central plane produces blackbody and reflected photons; the rest pass to the other side of the hot slab (which we simulate by mirror reflection). The value of fcl close to zero means that there is very little reprocessing and very few seed soft photons from the clouds cooling the hot flow. Similarly to the slab corona, we search forthe equilibrium coronal electron temperature by iterations requiring that the Comptonized coronal flux through the central plane times the factor fcl equals the flux emitted by the cold clouds.

Figure 4 shows example spectra in the case of hot flow with embedded cold clouds. The decreasing covering factor leads to the hardening of the spectrum because the same power has to be emitted by the flow with far fewer soft seed photons. This is possible only if the coronal temperature is higher (i.e. the Comptonization spectrum is harder) as demonstratedin Fig. 5. We also note that the reflection spectrum and, in particular, the line strength vary with fcl because the softer spectra correspond to effectively greater ionization (see Sect. 3.3 in García et al. 2013). In order to reproduce the hard-state spectra of GX 339–4 with Γ ≈ 1.8, the factor fcl ~ 0.4 is required, which is consistent with the results of Zdziarski et al. (1998).

Figure 5 shows the dependence of the spectral index and the coronal temperature on the ionization parameter for different values of fcl. The general trend is the following: smaller fcl leads to harder spectrum: it evolves from Γ ≈ 2.2 for fcl = 1 to Γ ≈ 1.6 for fcl = 0.2. The coronal temperature correspondingly increases from 80 keV to 140 keV. As in the case of the slab corona, increasing ionization leads to a higher reflection albedo and a smaller soft photon luminosity, which result in hardening of the spectra. Thus, as log ξ changes from 1 to 3, Γ decreases by 0.2. We see in the context of this model that in order to reproduce the observed spectral indices Γ ≈ 1.7 ± 0.1 in the hard state, we either need to have rather small fcl < 0.4 or very high ionization with log ξ > 3. As in the slab-corona model, high ξ is unlikely to be a solution because no high-ionization spectral features are seen in the hard-state data (Plant et al. 2014; Parker et al. 2015; Basak & Zdziarski 2016; Steiner et al. 2017; Basak et al. 2017). Thus, the only physically realistic solution is to assume small fcl. However, the reflection amplitude in this case, ∝ fcl exp(−τT∕ cos i), is much smaller than the typically observed hard-state reflection fraction of ≳ 0.2 (see e.g. Parker et al. 2015; Fürst et al. 2015; Basak & Zdziarski 2016; Basak et al. 2017 for recent studies) requiring an additional source of reflection, likely from the truncated disc.

The case of a geometrically thick hot flow overlapping with an outer disc is similar, except that there is a significant fraction of blackbody photons that are not scattered in the hot flow, but escape to the observer directly. Thus, the relative amplitude of the soft component is higher.

thumbnail Fig. 4

Self-consistent spectra from hot thermal plasmas mutually interacting with cold media without internal dissipation. The Thomson optical depth of the hot flow down to the central plane is τT = 0.4 and the spectra are observed at inclination of i = 60°. Red solid, blue dotted, green dashed, brown dot-dashed and pink triple-dot-dashed lines correspond to the covering fraction of the cold clumps of fcl = 1.0, 0.8, 0.6, 0.4 and 0.2, respectively. The ionization parameter here is log ξ = 1. All spectra correspond to the same bolometric luminosity.

thumbnail Fig. 5

Same as Fig. 3, but for the hot flow with cold clumps. Here τT = 0.4 and squares, circles, triangles, diamonds and stars correspond to cold clump covering factors fcl = 1.0, 0.8, 0.6, 0.4 and 0.2, respectively.

2.3 Patchy and outflowing corona models

As we showed above, a reduction in the reprocessing efficiency of geometries where the hot flow is surrounded by the cold disc and in the hot flow with cold clumps produces hard spectra. A similar reduction can be achieved in patchy corona models where the hot phase is situated in localized regions (flares) above the cold accretion disc (Galeev et al. 1979; Haardt et al. 1994; Stern et al. 1995; Svensson 1996). In this case, the higher the flaring region above the disc, the smaller the fraction of reprocessed photons that come back to the flare. This reduces the cooling and produces harder spectra. This model, however, predicts an anticorrelation between the amplitude of reflection R and spectral index Γ in contradiction with the observations (Zdziarski et al. 2003). Furthermore, as the flares emit nearly isotropically, the reprocessed blackbody component from the accretion disc is expected to be of a comparable luminosity to the observed Comptonized emission from the flare, also contradicting the data, which show that the soft component is an order of magnitude less luminous (e.g. Frontera et al. 2001; Kajava et al. 2016).

Alternatively, the dynamical outflowing corona (Beloborodov 1999; Malzac et al. 2001) can satisfy the observational spectral constraints. In this model, the higher the velocity of the flaring material away from the disc, the more the radiation is beamed away from the disc and the lower the reprocessing efficiency leading to harder spectra. This model can reproduce the observed R–Γ correlation (Zdziarski et al. 1999, 2003; Gilfanov et al. 1999; Revnivtsev et al. 2001; Gilfanov 2010) and is consistent with the observed range of spectral indices and electron temperatures. However, it requires the velocity to decrease on transition from the hard to the soft state; the reasons for this are unknown.

Furthermore, this model heavily relies on the assumption that the energy dissipation occurs exclusively in the corona with no intrinsic disc dissipation. The intrinsic disc radiation in this model will significantly influence the resulting spectra, because the energy density of the disc radiation as observed in the comoving frame of the plasma will actually be higher than that measured in the disc frame relevant for the static corona models (see Sect. 2.4 below). Thus, the model also requires that the location of the energy release should change from the cold disc in the soft state to the corona in the hard state for reasons that are not understood or discussed.

2.4 Intrinsic disc dissipation

The models discussed above assumed that all the energy is dissipated in the hot phase (corona or hot flow) and there is no intrinsic viscous dissipation happening within the disc or cold clouds, which is not realistic. An additional dissipation of energy in the dense matter leads to an increase in the flux of soft photons, in particular those irradiating the hot plasma. This will then lead to a decrease in the electron temperature and the associated softening of the scattered spectra (at a given τT).

In the framework of the model of the hot flow overlapping with an outer disc (see e.g. Zdziarski et al. 1999), an intrinsic dissipation inthe disc was considered by Ibragimov et al. (2005), who also took into account the decrease in the strength of the observed Compton-reflection feature due to scattering in the corona. For completeness, here we also calculate the Comptonization spectra from the slab-corona model and the hot flow with cold clumps model taking into account the intrinsic disc dissipation.

The main new parameter in these models is the ratio of the luminosity associated with the intrinsic disc dissipation Lintr to that of the total luminosity. The intrinsic disc spectrum is assumed to have the blackbody shape of temperature 50 eV. In Fig. 6, we show the spectral slopes and the equilibrium electron temperatures as functions of this ratio for three different pairs of ionization parameter and cold clump covering factor. For large ionization parameter of log ξ = 4 in a slab-corona (sandwich) model (fcl = 1.0), the intrinsic disc dissipation rate cannot exceed 10% of the total luminosity, otherwise the Comptonization spectrum becomes too soft. In the case of the hot flow with cold clumps of covering factor fcl = 0.1, we get a similar constraint on the intrinsic disc luminosity nearly independently of the ionization parameter (compare squares and triangles in Fig. 6). Thus, it is clear that in order to be consistent with the observed range of spectral slopes and electron temperatures, for the considered models where all seed soft photons propagate through the hot phase, the intrinsic disc dissipation cannot exceed 10% of the total luminosity. Additionally, as shown before, for the sandwich-corona model the condition of log ξ > 3.5 has to be satisfied, while for the hot flow with cold clumps model either log ξ > 3.5 or fcl < 0.4 are required.

The intrinsic disc dissipation also affects the spectra from the outflowing corona discussed in Sect. 2.3. Relativistic aberration effects increase the solid angle of the disc seen from the outflowing medium (up to a factor of 2 for an infinitely flat disc and an ultrarelativistic outflow velocity). In addition, the photon spectrum and total energy are modified. The total spectrum thus softens in the case of no dissipation, similar to the case of zero outflow velocity (sandwich corona). The strong dependence of the resulting spectral slope on the dissipation fraction is complemented by the dependence on the outflow velocity (see e.g. Fig. 8 in Malzac et al. 2001), requiring some fine-tuning of parameters to place the modelled spectral slopes in the observed range.

thumbnail Fig. 6

Same as Fig. 3, but the photon index and the equilibrium electron temperature of the hot phase are shown as a function of the ratio of intrinsic disc dissipation rate to the total luminosity. Squares show the case of log ξ = 1, fcl = 0.1; circles log ξ = 3, fcl = 0.1; and triangles log ξ = 4, fcl = 1.0. In all cases τT = 0.4.

2.5 Hot flow with cyclo-synchrotron radiation

The above discussion concerns models in which the disc blackbody is the sole source of seed photons for Comptonization. On the other hand, an additional source is provided by the cyclo-synchrotron process, which is likely to be copious given that magnetic fields are required in accretion flows for an outward transfer of the angular momentum. Various types of radiatively inefficient accretion models rely exclusively on partially self-absorbed cyclo-synchrotron photons as seeds for Comptonization, in particular the advection-dominated accretion flow (ADAF) models, starting with Narayan & Yi (1994), and followed by a very large number of works (see Yuan & Narayan 2014 for a review).

These models assume that the electron distribution in the hot flow is Maxwellian. This results in a relatively low rate of the synchrotron emission at low temperatures,kTe ≲ 100 keV (Wardziński& Zdziarski 2000). Low synchrotron luminosity in turn results in inefficient electron cooling (in the absence of other sources, such as cold accretion disc), and the electron temperature increases to the point where the radiative losses become comparable to characteristic heating rates, which for ADAF models is typically achieved at kTe ≳ 200 keV (e.g. Yuan & Zdziarski 2004; Yuan et al. 2007). However, the electron distribution is very likely to have a non-thermal tail (i.e. to be hybrid; Poutanen& Coppi 1998; Coppi 1999; Gierliński et al. 1999) due to the energy release proceeding in the form of magnetic reconnection (e.g. Sironi & Spitkovsky 2014; Beloborodov 2017). This is supported by the observations of MeV emission in both hard and soft states of black holes (McConnell et al. 1994, 2002; Ling et al. 1997; Grove et al. 1998; Wardziński et al. 2002; Droulans et al. 2010; Zdziarski et al. 2017). The presence of the high-energy electron tail then significantly enhances the synchrotron emission (Wardziński& Zdziarski 2001).

The presence of an additional source of seed photons not directly connected to the reflecting surface of the disc/clumps allows a certain freedom in the spectral modelling. It easily resolves the photon starving problem. Calculations of emitted spectra from hybrid electrons are, however, more complex compared to the case of fully thermal plasma, as not only the electron temperature, but also the whole distribution function for electrons has to be found self-consistently with the photon field they produce. Hybrid models proved to give a good quantitative description for the X-ray/γ-ray spectra of both black hole binaries and active galactic nuclei (Poutanen & Vurm 2009; Malzac & Belmont 2009; Veledina et al. 2011b), and at longer wavelengths (Veledina et al. 2013; Poutanen & Veledina 2014).

We simulate the spectrum of the inner parts of hot accretion flow, which is approximated by a sphere with radius R = 3 × 107 cm, corresponding to 10 Schwarzschild radii for a 10 M black hole, of the Thomson optical depth of τT = 1. Particle and photon distributions are assumed to be homogeneous and isotropic. The magnetic field is assumed to be tangled and characterized by its strength of B = 5 × 105 G. We assume that the energization proceeds through the injection of electrons with power-law energy distribution, e (γ) ∝ γ−Γinj; we take Γinj = 2.5. The injected electron distribution evolves because of the cooling and heating by cyclo-synchrotron, Compton scattering, bremsstrahlung and pair production, and by electron-electron Coulomb collisions. The same number of electrons is removed from the evolved distribution to keep the Thomson optical depth constant. The injected energy is regulated through the total amount of injected particles related to the total emitted luminosity of L = 1037 erg s−1. The spectrum is calculated by solving a system of coupled relativistic kinetic equations for photons, electrons, and positrons. Further details of modelling can be found in Vurm & Poutanen (2009).

The resulting spectrum is shown in Fig. 7. The 2–10 keV spectral index is Γ = 1.76 and the equilibrium temperature (of the Maxwellian part of electron distribution) is kTe = 53 keV. The dotted curve shows the continuum that is produced by Comptonization of the synchrotron radiation from the hot flow. An additional component related to the reflection and reprocessing of the Comptonized radiation in the outer cold disc is also added (dashed curve). In this case, we do not consider mutual interaction between reflected/reprocessed radiation and the hot flow because these components are likely physically separated. As was shown in Poutanen & Vurm (2009) and Veledina et al. (2013), the hybrid thermal/non-thermal model predicts rather stable slopes of the Comptonized continuum if the total optical depth of the hot flow is of the order unity. The model is thus consistent with observations.

thumbnail Fig. 7

Example of a synchrotron self-Compton spectrum from hot plasma with a hybrid electron distribution. The components due to reflection/reprocessing with the reflection fraction of 0.3 and the ionization parameter of log ξ = 1 and the blackbody of the temperature of 50 eV due to the reprocesses emission are shown together by the blue dashed curve. The escaping synchrotron spectrum is shown by the green dot-dashed line, synchrotron Comptonization spectrum is shown by the red dotted curve, and the total spectrum is given by the red solid curve.

3 Discussion

3.1 Comparison of reflection models

Many previous models of the X-ray emission from the accreting black holes used the pexrav/pexriv model (Magdziarz & Zdziarski 1995). This is the simplest reflection model that correctly represents the angular-distribution of the reflected radiation using exact Klein–Nishina cross-section. It provides analytical formulae for the Green function for Compton reflection based on a Monte Carlo simulation of photon propagation in a homogeneous slab with an isotropic source above it. The reflected spectrum is computed as a convolution of the Green function with any incident spectrum (not only cut-off power law). The pexrav assumes neutral matter, while the pexriv model is for the ionized material. The latter model ignores radiative transfer effects for the computation of the ionization/recombination balance (Done et al. 1992).

Radiative transfer effects have been accounted for by various authors (e.g. Ross & Fabian 1993; García & Kallman 2010; García et al. 2013). For our analysis here we use the xillverCp model, which is provided as a tabulated set of reflection spectra for a family of the incident Comptonization continua described by the approximate nthcomp model. To test the role of the reflection model and the effect of the Comptonization model we ran two test cases with log ξ = 1 and 3 and computed self-consistent sandwich-corona spectra using pexriv and xillverCp reflection models.

In the first case (log ξ = 1), the greatest difference is observed in the reflection spectrum above 30 keV (see Fig. 8a). The most likely reasons being that the Klein–Nishina effects are not accounted for in xillverCp and that the incident photons are injected at 45° to the slab normal, whereas pexriv assumes an isotropic source. The difference at energies below 10 keV is likely caused by the choice of photoelectric cross sections (see also Fig. 1 in Niedźwiecki et al. 2016). However, the equilibrium coronal electron temperature differs by less than 1% and the total spectra are nearly identical up to 20 keV. Thus the self-consistent spectrum does not depend significantly on the choice of the reflection or Comptonization model. At higher ξ (see Fig. 8b), the xillverCp reflected spectrum has substantially more photons in the 1–10 keV range and fewer photons below 1 keV than the pexriv model. However, this difference in the reflection spectrum only results in a 5% higher electron temperature and somewhat harder self-consistent Comptonization spectra by ΔΓ = 0.18.

One of the major drawbacks of the considered reflection models is the assumption of the constant density matter. A better approximation would be a constant-pressure material (Malzac et al. 2005). In reality, of course, the disc atmosphere is in hydrostatic equilibrium with both density and pressure decreasing towards the surface. The X-ray illumination also changes the atmosphere structure causing the top layers to expand, further lowering their density. Thermal ionization instability can then lead to a rapid change in theionization state of the material. The top layers are heated to nearly Compton temperature determined by the balance between Compton scattering heating and cooling and bremsstrahlung emission. These layers are almost completely ionized, while the layers below are much denser and nearly neutral. The reflection models that take into account all these effects (Nayakshin et al. 2000; Nayakshin & Kallman 2001; Done & Nayakshin 2001) are much more time consuming to compute. This makes them less flexible and limits their application to the data (Barrio et al. 2003; Done & Gierliński 2006). The hydrostatic-equilibrium models, however, are not much different from the constant-density model at low ionization log ξ < 2 because in that case the gas is nearly neutral (Nayakshin et al. 2000). At higher ionization, the approximation of the constant-density slab used by xillverCp predicts too few soft photons compared to a more physically realistic case of reflection from a constant-pressure slab (see Fig. 5 in Malzac et al. 2005) or a slab in hydrostatic equilibrium (see Fig. 8 in Nayakshin et al. 2000). Thus, we conclude here that the results obtained using the xillverCp reflection model are robust for low-ionization discs, but at high ionizations they predict spectra that are too hard.

thumbnail Fig. 8

Comparison of the self-consistent spectra from the sandwich corona for two reflection models, xillverCp (solid blue curves) and pexriv (dashed red curves), for two different ionization parameters (a) log ξ = 1 and (b) log ξ = 3. The reflected and reprocessed blackbody spectra are shown by thin lines and the total spectrum by the thick lines. The upper and lower curves correspond to two inclinations, i = 18° and 60°, respectively.

3.2 Comparison to previous sandwich-corona models

Energy balance in the slab-corona geometry was first studied by Haardt & Maraschi (1991, 1993), but their finding was that the resulting spectra are softer than those of the hard state of Galactic BH XRBs. As we discussed above, this conclusion was confirmed in a number of subsequent papers, as it is here.

On the other hand, it was proposed that the hard-state X-ray spectra in this geometry can be recovered by accounting for the reflected emission of the cold disc, which is then Compton up-scattered again as it passes through the hot medium (e.g. Young et al. 2001 and more recently Steiner et al. 2017). In the latter paper, it was argued that the truncated disc/hot flow geometry is not capable of reproducing the luminosities observed in the bright hard state of GX 339–4 for large truncation radii as this would require super-Eddington mass accretion rates. We analyse these arguments and compare their findings to our results.

The following set-up is adopted in Steiner et al. (2017). The optically thick accretion disc in the hard-state accretion flow is described by the model of Shakura & Sunyaev (1973), but with the truncation radius being a free parameter, constrained to rinrISCO, and with the zero-stress boundary condition imposed at rin (as implemented in the spectral-fitting routine of Zimmerman et al. 2005). Then, the disc luminosity from the viscous dissipation from both sides of the disc is (in the Newtonian approximation) Lintr=M˙dc22rin/Rg,\begin{equation*} L_{\rm{intr}}=\frac{\dot M_{\rm{d}} c^2}{2 r_{\rm{in}}/R_{\rm{g}}},\end{equation*}(3)

implying that the accretion rate required to account for a given disc luminosity is drin, where d is the mass flow rate through the disc.

The authors then consider Compton scattering of a fraction of the disc blackbody photons in the corona requiring conservation of photons, and using either a simple description of Comptonization as an e-folded power law or the thermal Comptonization model of Zdziarski et al. (1996). The Comptonized spectrum is emitted partly outside and partly towards the disc, which then gives rise to the reflection component, using the model of García & Kallman (2010). We denote the total Comptonization luminosity emitted both outside and towards the disc as Lc, and the fraction of photons going towards the disc as RF. The reflection component is smeared relativistically (Dauser et al. 2010), and then Compton scattered in the corona, which is a major feature of the model. As we noted above, many theoretical self-consistent models of hot plasma interacting with cold gas (see Stern et al. 1995; Poutanen & Svensson 1996; Malzac et al. 2001) have accounted for scattering of the reflected radiation. However, this is not accounted for in the xspec models like compps. The authors argue that the inclusion of Compton scattering of the reflected photons allows us to increase the relative strength of the fitted reflected component to values expected in the geometry of the corona lying above the disc. The value of reflection fraction, RF, fitted to the GX 339–4 data used in this study is ≃0.8.

In the approach by Steiner et al. (2017) the seed photons for Compton scattering are dominated by the disc blackbody photons generated by viscous dissipation within the disc, which implies that LcLintr. At the same time a large truncation radius implies a large accretion rate in order to generate enough blackbody photons far away from the black hole, with the relationship of drin. This argument is very similar to the well-known photon starvation problem for the truncated disc/hot flow geometry. However, we point out that this problem appears only if we assume that the seed blackbody luminosity originates exclusively from viscous dissipation. In the considered slab-corona geometry, the disc is irradiated by the coronal emission, which in the hard state of BH XRBs has a much higher luminosity, typically one order of magnitude, than that of the disc. Steiner et al. (2017) take into account the reflection from the disc, with the reflection fraction of ≃0.8, i.e. with similar fluxes going to the disc and outside in their best-fit model for GX 339–4. The typical albedo for the reflection in medium/hard X-rays is a ≃ 0.2–0.3, and thus most of the irradiating flux is actually absorbed and then re-emitted as soft photons, with only a relatively small fraction of it Compton-scattered within the disc and emitted outside, resulting in LirrLintr. Therefore, the disc will emit not only the viscously dissipated power (given by d), but also the power from the irradiation by the corona, which is linked to the reflection fraction (see e.g. Gierliński et al. 2008). Therefore, d does not need to be high for a sufficient emission of blackbody photons.

Observationally, it has been found that the disc component in the spectra of BH transients in the hard state does not follow the disc-blackbody dependence of the luminosity on the inner radius (see e.g. Gierliński et al. 2008; Basak & Zdziarski 2016), which confirms the dominant role of the disc irradiation. Irradiation is important in the inner parts of the disc and also in the outer parts. It has also been a general feature of accretion disc models since Shakura & Sunyaev (1973) and it is crucial to explain the light curves, optical fluxes, and recurrence of BH transient systems (e.g. van Paradijs 1996; King et al. 1996; Gierliński et al. 2009; Coriat et al. 2012; Poutanen et al. 2014). In particular, irradiation reduces the critical luminosity below which BH or neutron-star LMXBs are transient by an order of magnitude.

Our results (Fig. 2) show that for a sandwich-corona geometry with zero viscous dissipation the spectra are soft. This should be contrasted with the model of the hard state of GX 339–4 presented in Fig. 4 and Table 1 in Steiner et al. (2017). There both high reflection and high scattering fractions imply a geometry close to the slab corona, but the resulting spectrum is hard. The reason for this discrepancy is our self-consistent calculation of the electron temperature of the corona and the resulting Comptonization spectra, omitted by Steiner et al. (2017). Their treatment of reflection also has multiple problems. First, the photon energy grid in the xillver and relxill models only extends down to 70 eV, as discussed in Sect. 2.1, but the number of photons in the reflected spectrum below that energy exceeds that above 70 eV by orders of magnitude (see Fig. 3 in García et al. 2013). This then implies that the normalization of the Comptonization continuum relative to the reflection spectrum is wrong by orders of magnitude and that Comptonization does not conserve photons. Furthermore, the limited photon energy range also leads to the angle-integrated reflected luminosity being different from the illuminating luminosity, i.e. the energy in these models is not conserved either. And finally, the angular distribution of reflected and Comptonized photons is different and therefore it is meaningless to demand the photon conservation for the spectra observed at one specific inclination.

We note here that spectra obtained for the geometries where the hot flow contains cold clumps or is overlapping with the disc (Fig. 4) can be made hard, similar to the observed ones. The problem of requiring a super-Eddington accretion rate for the truncated disc geometry disappears if we take into account soft photons from other sources, such as reprocessing and internal synchrotron emission.

3.3 Preferred geometry

As we argued above, sandwich, patchy, and outflowing corona models need a lot of fine-tuning to reproduce commonly observed hard spectra and electron temperatures. Moreover, the models assuming the disc is going down to the ISCO are inconsistent with the spectral modelling of reflection and the Kα line (Gilfanov et al. 2000; Plant et al. 2015; Basak & Zdziarski 2016), which suggest rin ~ 100Rg in the hard state. The model where the disc does not overlap with the hot flow is yet another extreme, expected to have too hard spectra. Three other models (hot flow with cold clumps, hot flow overlapping with the cold disc, and hot flow with cyclo-synchrotron photons) are degenerate in the X-ray spectral data alone. However, they are expected to make different predictions for the X-ray short- and long-term variability, optical/infrared (OIR) spectra, and interconnection of emission at different wavelengths. Many discoveries in terms of multiwavelength behaviour were made after those geometries were proposed, and properties for some of them were not discussed in application to the new data. We now summarize the timing and multiwavelength properties critical to discriminate between the geometries.

3.3.1 X-ray timing properties

Significant hard time lags between different X-ray energy bands were detected, with a peculiar power-law dependence on frequency (Miyamoto et al. 1988; Nowak et al. 1999) and a logarithmic dependence on the energy band (Kotov et al. 2001). This behaviour was interpreted both in terms of flares with pivoting power law (Poutanen & Fabian 1999), corresponding to a patchy corona geometry, and in terms of radial dependence of the emitted spectrum (Kotov et al. 2001), which is expected in the truncated disc/hot flow geometry. The observed power-law dependence of the X-ray variability power on Fourier frequency (e.g. Nowak et al. 1999) and the linear dependence of root mean square variability amplitude on flux (Uttley & McHardy 2001) argue against the model where X-rays are produced in the flares, and support the picture of accretion rate fluctuations propagating in the radially stratified hot medium (Lyubarskii 1997; Kotov et al. 2001). The inferred radial dependence of the X-ray spectral hardness, with harder spectra coming from the innermost parts, naturally arises in the geometries with hot flow overlapping with the disc and in the hot flow with non-thermal particles. On the other hand, this dependence is not expected to be automatically satisfied in the case of clumpy hot flow.

Another powerful technique used to probe the accretion flow geometry is X-ray reverberation (Reynolds et al. 1999; Poutanen 2002; Uttley et al. 2014). Measurements of the time delays between X-ray continuum and the signatures of the cold accretion disc, namely a soft thermal component and iron Kα line (De Marco et al. 2015, 2017; De Marco & Ponti 2016), suggest the truncated disc geometry, where the truncation radius decreases with rising luminosity, i.e. towards the soft state.

An independent line of arguments for a truncated disc geometry comes from the studies of the equivalent width of the iron Kα line in a range of Fourier frequencies (Gilfanov et al. 2000). A substantial suppression of the iron line variability at high Fourier frequencies corresponding to the innermost parts of the accretion flow was found, which implies that the cold disc truncation radius is around ~ 100Rg in the hard state.

These findings disfavour clumpy accretion flow geometry, unless these clumps are fully ionized and are acting as mirrors from the spectral point of view. On the other hand, such high values for the truncation radius results in small number of seed photons entering the zone of major energy release (within ~10Rg), leading to the standard problem with photon starvation. Thus, high truncation radii in the hard state disfavour the models where seed photons come solely from the disc.

3.3.2 Multiwavelength properties

We now discuss various possibilities involving truncation of the disc. The major question here concerns the source of seed photons, either thermal (viscous or reprocessed) photons from clumps/accretion disc or synchrotron photons. At the outburst peak, close to the soft state, the accretion disc is likely to be the major source of seed photons. As the outburst declines, the spectra are found to harden, and the trend is reversed at very low luminosities (below ~ 10−3LEdd, Sobolewska et al. 2011), i.e. the spectra are softer for lower luminosities. Spectral hardening is expected from the outward movement of the disc inner radius; however, the softening is not, and it instead requires an additional source of synchrotron seed photons (Veledina et al. 2011b). The softening can be explained in the clumpy hot flow scenario if ionization drops as the luminosity decreases. It is difficult to explain the constant hardness ratio during the rising phase (when the observed luminosity is in the range ~ 10−3 − 10−2LEdd) using the geometry where seed photons come solely from the disc overlapping with the hot flow or from the cold clumps. The synchrotron self-Compton mechanism is preferred here. The switch of the source of seed photons likely occurs close to the point where the spectra start to soften on the rising phase, or achieve the hardest spectra on the decline stage. The critical luminosity (of about 2% LEdd) where these two seed photon populations switch has recently been identified at the decline stage of SWIFT J1753.5–0127 (Kajava et al. 2016).

Evolution of the X-ray spectrum towards and from the soft states proceeds through the intermediate states when the spectra are flat (Γ ~ 2). These spectraare not well fitted with one Comptonization continuum, and two continua are required instead (Ibragimov et al. 2005; Shidatsu et al. 2011). It is then natural to interpret these continua as being produced in two physically separated regions with different sources of seed photons: discand synchrotron. These two continua are expected to change in response to fluctuations in mass accretion rate, one delayed with respect to another, producing an interference picture in the power spectrum, which we see in the data (Veledina 2016).

In the hard state, the scenario with the synchrotron seed photons predicts that OIR radiation is partially produced by the synchrotron component from the hot accretion flow, thus properties of emission at longer wavelengths can serve as a diagnostics of accretion geometry. The OIR excesses above the possible accretion disc fluxes were reported in a number of sources (Chaty et al. 2003; Chiang et al. 2010; Gandhi et al. 2011). An excess flux at these wavelengths can be produced by the jet (Blandford & Königl 1979); however, in some sourceseither the extrapolation of the power law from radio significantly underpredicts OIR fluxes or the extrapolated slope is different (Chaty et al. 2003; Durant et al. 2009). The jet scenario also has major difficulties in explaining the optical/X-ray cross-correlation function with broad anti-correlation observed in SWIFT J1753.5–0127 (Durant et al. 2008, 2011), which perhaps possesses the faintest radio jet. On the other hand, this shape of the cross-correlation function is expected to be seen if the hot accretion flow contributes to both optical and X-ray energy bands (Veledina et al. 2011a). The observed evolution of the shape of the cross-correlation function, from demonstrating a positive correlation at the outburst peak to the anti-correlation feature in the outburst tail, is consistent with the scenario where the cold disc recedes and is replaced by the hot flow (Veledina et al. 2017). Neither high OIR fluxes, nor short-term variability in the hard state can be understood in terms of the scenario of clumps/intersecting disc.

The long-term multiwavelength behaviour can also be used to constrain the accretion geometry. The outburst X-ray light-curve resembles a fast rise exponential decay profile, and the light-curves in OIR energy bands often demonstrate additional prominent flares lasting a substantial fraction of the outburst duration. The origin of these flares has been discussed in a context of hot flow and jet synchrotron emission (Buxton et al. 2012; Poutanen et al. 2014). The jet hypothesis is inconsistent with the spectrum of the flare at least in XTE J1550–564, which was shown to be very hard when the flare starts (Poutanen et al. 2014). On the other hand, the evolution of the flare spectrum is in good agreement with the hot flow scenario. Neither a cold accretion disc nor emission of cold clumps is able to account for the appearance of these flares, which favours the scenario of hot flow with cyclo-synchrotron radiation in the hard state.

4 Summary

We have modelled Comptonization and photon reprocessing processes in various geometries relevant to the hard state of BH XRBs. We have used a modern (constant-density) reflection code, xillverCp, to account for Compton reflection and reprocessing. Our main results are as follows.

Our self-consistent treatment confirms earlier studies, in particular that a static disc-corona is not a viable geometry for the hard state of BH X-ray XRBs, even accounting for scattering of reflection radiation in the corona. In the absence of intrinsic disc dissipation (energy is dissipated entirely in the corona), the sandwich-corona models require a high ionization parameter (log ξ ≳ 3.5) to agree with the observed spectral slopes. More physically realistic reflection models that account for density gradients and allow for thermal instability to operate in the disc, predict spectra softer than the observed ones even at such high ionizations. Alternative models, such as patchy or outflowing corona and clumpy hot flow can reproduce the observed spectra and electron temperatures after fine-tuning the parameters. We show that the spectral slopes heavily depend on the fraction of energy dissipated intrinsically in the disc. Assuming only a 10% energy release proceeding through the viscous heating in the cold disc, we obtain spectra that are too soft for the above-mentioned models. Given that the energy liberation almost entirely proceeds through the viscous heating in the soft state, it is not clear why such a tight upper boundary should apply to the hard state.

We conclude that the disc has to be truncated before the ISCO, and that it is replaced by some kind of a hot flow. This assures an agreement with the X-ray observations of BH XRBs in the hard state, namely their spectral slopes, the range of electron temperatures, and the ionization levels of the cold matter. Some mixing of the cold and hot media, either the disc overlapping with the hot flow or the hot flow containing some cold medium, can reproduce the observed ranges of parameters but does not reproduce the observed evolution of these parameters over the course of the outburst (nearly constant hardness ratio at moderate LLEdd and spectral softening towards low luminosities). We also reiterate the importance of the cyclo-synchrotron emission by hybrid electrons (mostly thermal with a small non-thermal tail) as the process providing both seed photons for Comptonization and low-energy (OIR) photons. The latter mechanism is favoured by the multiwavelength spectra and variability properties.

Acknowledgements

We thank A. Beloborodov and J. Steiner for valuable discussions, and J. García for help with the xillverCp model. This research has been supported in part by the Polish National Science Centre grants 2013/10/M/ST9/00729, 2015/18/A/ST9/00746, and 2016/21/P/ST9/04025; the Academy of Finland grant 309308; and grant 14.W03.31.0021 of the Ministry of Education and Science of the Russian Federation.

References

  1. Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, eds. G. H. Jacoby, & J. Barnes, ASP Conf. Ser., 101, 17 [NASA ADS] [Google Scholar]
  2. Barrio, F. E., Done, C., & Nayakshin, S. 2003, MNRAS, 342, 557 [NASA ADS] [CrossRef] [Google Scholar]
  3. Basak, R., & Zdziarski, A. A. 2016, MNRAS, 458, 2199 [NASA ADS] [CrossRef] [Google Scholar]
  4. Basak, R., Zdziarski, A. A., Parker, M., & Islam, N. 2017, MNRAS, 472, 4220 [NASA ADS] [CrossRef] [Google Scholar]
  5. Beloborodov, A. M. 1999, ApJ, 510, L123 [NASA ADS] [CrossRef] [Google Scholar]
  6. Beloborodov, A. M. 2017, ApJ, 850, 141 [NASA ADS] [CrossRef] [Google Scholar]
  7. Blandford, R. D., & Königl, A. 1979, ApJ, 232, 34 [NASA ADS] [CrossRef] [Google Scholar]
  8. Burke, M. J., Gilfanov, M., & Sunyaev, R. 2017, MNRAS, 466, 194 [NASA ADS] [CrossRef] [Google Scholar]
  9. Buxton, M. M., Bailyn, C. D., Capelo, H. L., et al. 2012, AJ, 143, 130 [NASA ADS] [CrossRef] [Google Scholar]
  10. Celotti, A., Fabian, A. C., & Rees, M. J. 1992, MNRAS, 255, 419 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chaty, S., Haswell, C. A., Malzac, J., et al. 2003, MNRAS, 346, 689 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Chiang, C. Y., Done, C., Still, M., & Godet, O. 2010, MNRAS, 403, 1102 [NASA ADS] [CrossRef] [Google Scholar]
  13. Coppi, P. S. 1999, in High Energy Processes in Accreting Black Holes, eds. J. Poutanen, & R. Svensson, ASP Conf. Ser., 161, 375 [NASA ADS] [Google Scholar]
  14. Coriat, M., Fender, R. P., & Dubus, G. 2012, MNRAS, 424, 1991 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dauser, T., Wilms, J., Reynolds, C. S., & Brenneman, L. W. 2010, MNRAS, 409, 1534 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dauser, T., García, J., Walton, D. J., et al. 2016, A&A, 590, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. De Marco, B., & Ponti, G. 2016, ApJ, 826, 70 [NASA ADS] [CrossRef] [Google Scholar]
  18. De Marco, B., Ponti, G., Muñoz-Darias, T., & Nandra, K. 2015, ApJ, 814, 50 [NASA ADS] [CrossRef] [Google Scholar]
  19. De Marco, B., Ponti, G., Petrucci, P. O., et al. 2017, MNRAS, 471, 1475 [NASA ADS] [CrossRef] [Google Scholar]
  20. Done, C., & Gierliński, M. 2006, MNRAS, 367, 659 [NASA ADS] [CrossRef] [Google Scholar]
  21. Done, C., Mulchaey, J. S., Mushotzky, R. F., & Arnaud, K. A. 1992, ApJ, 395, 275 [NASA ADS] [CrossRef] [Google Scholar]
  22. Done, C., & Nayakshin, S. 2001, ApJ, 546, 419 [NASA ADS] [CrossRef] [Google Scholar]
  23. Dove, J. B., Wilms, J., Maisack, M., & Begelman, M. C. 1997, ApJ, 487, 759 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dovčiak, M., & Done, C. 2016, Astron. Nachr., 337, 441 [NASA ADS] [CrossRef] [Google Scholar]
  25. Droulans, R., Belmont, R., Malzac, J., & Jourdain, E. 2010, ApJ, 717, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  26. Dubus, G., Hameury, J.-M., & Lasota, J.-P. 2001, A&A, 373, 251 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Durant, M., Gandhi, P., Shahbaz, T., et al. 2008, ApJ, 682, L45 [NASA ADS] [CrossRef] [Google Scholar]
  28. Durant, M., Gandhi, P., Shahbaz, T., Peralta, H. H., & Dhillon, V. S. 2009, MNRAS, 392, 309 [NASA ADS] [CrossRef] [Google Scholar]
  29. Durant, M., Shahbaz, T., Gandhi, P., et al. 2011, MNRAS, 410, 2329 [NASA ADS] [CrossRef] [Google Scholar]
  30. Esin, A. A., McClintock, J. E., & Narayan, R. 1997, ApJ, 489, 865 [NASA ADS] [CrossRef] [Google Scholar]
  31. Esin, A. A., Narayan, R., Cui, W., Grove, J. E., & Zhang, S.-N. 1998, ApJ, 505, 854 [NASA ADS] [CrossRef] [Google Scholar]
  32. Frontera, F., Zdziarski, A. A., Amati, L., et al. 2001, ApJ, 561, 1006 [NASA ADS] [CrossRef] [Google Scholar]
  33. Fürst, F., Nowak, M. A., Tomsick, J. A., et al. 2015, ApJ, 808, 122 [NASA ADS] [CrossRef] [Google Scholar]
  34. Galeev, A. A., Rosner, R., & Vaiana, G. S. 1979, ApJ, 229, 318 [NASA ADS] [CrossRef] [Google Scholar]
  35. Gandhi, P., Blain, A. W., Russell, D. M., et al. 2011, ApJ, 740, L13 [NASA ADS] [CrossRef] [Google Scholar]
  36. García, J., & Kallman, T. R. 2010, ApJ, 718, 695 [NASA ADS] [CrossRef] [Google Scholar]
  37. García, J., Dauser, T., Reynolds, C. S., et al. 2013, ApJ, 768, 146 [NASA ADS] [CrossRef] [Google Scholar]
  38. Gierliński, M., Zdziarski, A. A., Poutanen, J., et al. 1999, MNRAS, 309, 496 [NASA ADS] [CrossRef] [Google Scholar]
  39. Gierliński, M., Done, C., & Page, K. 2008, MNRAS, 388, 753 [NASA ADS] [CrossRef] [Google Scholar]
  40. Gierliński, M., Done, C., & Page, K. 2009, MNRAS, 392, 1106 [NASA ADS] [CrossRef] [Google Scholar]
  41. Gilfanov, M. 2010, in The Jet Paradigm, ed. T. Belloni (Berlin: Springer Verlag), Lect. Notes Phys., 794, 17 [Google Scholar]
  42. Gilfanov, M., Churazov, E., & Revnivtsev, M. 1999, A&A, 352, 182 [NASA ADS] [Google Scholar]
  43. Gilfanov, M., Churazov, E., & Revnivtsev, M. 2000, MNRAS, 316, 923 [NASA ADS] [CrossRef] [Google Scholar]
  44. Grove, J. E., Johnson, W. N., Kroeger, R. A., et al. 1998, ApJ, 500, 899 [NASA ADS] [CrossRef] [Google Scholar]
  45. Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 [NASA ADS] [CrossRef] [Google Scholar]
  46. Haardt, F., & Maraschi, L. 1993, ApJ, 413, 507 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  47. Haardt, F., Maraschi, L., & Ghisellini, G. 1994, ApJ, 432, L95 [NASA ADS] [CrossRef] [Google Scholar]
  48. Ibragimov, A., Poutanen, J., Gilfanov, M., Zdziarski, A. A., & Shrader, C. R. 2005, MNRAS, 362, 1435 [NASA ADS] [CrossRef] [Google Scholar]
  49. Kajava, J. J. E., Veledina, A., Tsygankov, S., & Neustroev, V. 2016, A&A, 591, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. King, A. R., Kolb, U., & Burderi, L. 1996, ApJ, 464, L127 [NASA ADS] [CrossRef] [Google Scholar]
  51. Kotov, O., Churazov, E., & Gilfanov, M. 2001, MNRAS, 327, 799 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lasota, J.-P., Narayan, R., & Yi, I. 1996, A&A, 314, 813 [NASA ADS] [Google Scholar]
  53. Ling, J. C., Wheaton, W. A., Wallyn, P., et al. 1997, ApJ, 484, 375 [NASA ADS] [CrossRef] [Google Scholar]
  54. Liu, B. F., Taam, R. E., Meyer-Hofmeister, E., & Meyer, F. 2007, ApJ, 671, 695 [NASA ADS] [CrossRef] [Google Scholar]
  55. Lyubarskii, Y. E. 1997, MNRAS, 292, 679 [NASA ADS] [CrossRef] [Google Scholar]
  56. Magdziarz, P., & Zdziarski, A. A. 1995, MNRAS, 273, 837 [NASA ADS] [CrossRef] [Google Scholar]
  57. Malzac, J., & Belmont, R. 2009, MNRAS, 392, 570 [NASA ADS] [CrossRef] [Google Scholar]
  58. Malzac, J., Beloborodov, A. M., & Poutanen, J. 2001, MNRAS, 326, 417 [NASA ADS] [CrossRef] [Google Scholar]
  59. Malzac, J., Dumont, A. M., & Mouchet, M. 2005, A&A, 430, 761 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. McConnell, M., Forrest, D., Ryan, J., et al. 1994, ApJ, 424, 933 [NASA ADS] [CrossRef] [Google Scholar]
  61. McConnell, M. L., Zdziarski, A. A., Bennett, K., et al. 2002, ApJ, 572, 984 [NASA ADS] [CrossRef] [Google Scholar]
  62. Miller, J. M., Tomsick, J. A., Bachetti, M., et al. 2015, ApJ, 799, L6 [NASA ADS] [CrossRef] [Google Scholar]
  63. Miyamoto, S., Kitamoto, S., Mitsuda, K., & Dotani, T. 1988, Nature, 336, 450 [NASA ADS] [CrossRef] [Google Scholar]
  64. Narayan, R., & Yi, I. 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
  65. Nayakshin, S., & Dove, J. B. 2001, ApJ, 560, 885 [NASA ADS] [CrossRef] [Google Scholar]
  66. Nayakshin, S., & Kallman, T. R. 2001, ApJ, 546, 406 [NASA ADS] [CrossRef] [Google Scholar]
  67. Nayakshin, S., Kazanas, D., & Kallman, T. R. 2000, ApJ, 537, 833 [NASA ADS] [CrossRef] [Google Scholar]
  68. Niedźwiecki, A., Zdziarski, A. A., & Szanecki, M. 2016, ApJ, 821, L1 [NASA ADS] [CrossRef] [Google Scholar]
  69. Novikov, I. D., & Thorne, K. S. 1973, in Black Holes (Les Astres Occlus), ed. C. Dewitt & B. S. Dewitt, 343 [Google Scholar]
  70. Nowak, M. A., Vaughan, B. A., Wilms, J., Dove, J. B., & Begelman, M. C. 1999, ApJ, 510, 874 [NASA ADS] [CrossRef] [Google Scholar]
  71. Parker, M. L., Tomsick, J. A., Miller, J. M., et al. 2015, ApJ, 808, 9 [NASA ADS] [CrossRef] [Google Scholar]
  72. Plant, D. S., Fender, R. P., Ponti, G., Muñoz-Darias, T., & Coriat, M. 2014, MNRAS, 442, 1767 [NASA ADS] [CrossRef] [Google Scholar]
  73. Plant, D. S., Fender, R. P., Ponti, G., Muñoz-Darias, T., & Coriat, M. 2015, A&A, 573, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Poutanen, J. 2002, MNRAS, 332, 257 [NASA ADS] [CrossRef] [Google Scholar]
  75. Poutanen, J., & Coppi, P. S. 1998, Phys. Scr., T77, 57 [NASA ADS] [CrossRef] [Google Scholar]
  76. Poutanen, J., & Fabian, A. C. 1999, MNRAS, 306, L31 [NASA ADS] [CrossRef] [Google Scholar]
  77. Poutanen, J., & Svensson, R. 1996, ApJ, 470, 249 [NASA ADS] [CrossRef] [Google Scholar]
  78. Poutanen, J., & Veledina, A. 2014, Space Sci. Rev., 183, 61 [NASA ADS] [CrossRef] [Google Scholar]
  79. Poutanen, J., & Vurm, I. 2009, ApJ, 690, L97 [NASA ADS] [CrossRef] [Google Scholar]
  80. Poutanen, J., Nagendra, K. N., & Svensson, R. 1996, MNRAS, 283, 892 [NASA ADS] [CrossRef] [Google Scholar]
  81. Poutanen, J., Krolik, J. H., & Ryde, F. 1997, MNRAS, 292, L21 [NASA ADS] [Google Scholar]
  82. Poutanen, J., Veledina, A., & Revnivtsev, M. G. 2014, MNRAS, 445, 3987 [NASA ADS] [CrossRef] [Google Scholar]
  83. Revnivtsev, M., Gilfanov, M., & Churazov, E. 2001, A&A, 380, 520 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Reynolds, C. S., Young, A. J., Begelman, M. C., & Fabian, A. C. 1999, ApJ, 514, 164 [NASA ADS] [CrossRef] [Google Scholar]
  85. Ross, R. R., & Fabian, A. C. 1993, MNRAS, 261, 74 [NASA ADS] [CrossRef] [Google Scholar]
  86. Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
  87. Shidatsu, M., Ueda, Y., Tazaki, F., et al. 2011, PASJ, 63, 785 [NASA ADS] [Google Scholar]
  88. Sironi, L., & Spitkovsky, A. 2014, ApJ, 783, L21 [NASA ADS] [CrossRef] [Google Scholar]
  89. Sobolewska, M. A., Papadakis, I. E., Done, C., & Malzac, J. 2011, MNRAS, 417, 280 [NASA ADS] [CrossRef] [Google Scholar]
  90. Steiner, J. F., García, J. A., Eikmann, W., et al. 2017, ApJ, 836, 119 [NASA ADS] [CrossRef] [Google Scholar]
  91. Stern, B. E., Begelman, M. C., Sikora, M., & Svensson, R. 1995, MNRAS, 272, 291 [NASA ADS] [CrossRef] [Google Scholar]
  92. Svensson, R. 1996, A&AS, 120, C475 [NASA ADS] [Google Scholar]
  93. Svensson, R., & Zdziarski, A. A. 1994, ApJ, 436, 599 [NASA ADS] [CrossRef] [Google Scholar]
  94. Torii, S., Yamada, S., Makishima, K., et al. 2011, PASJ, 63, S771 [NASA ADS] [CrossRef] [Google Scholar]
  95. Uttley, P. & McHardy, I. M. 2001, MNRAS, 323, L26 [NASA ADS] [CrossRef] [Google Scholar]
  96. Uttley, P., Cackett, E. M., Fabian, A. C., Kara, E., & Wilkins, D. R. 2014, A&ARv, 22, 72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  97. van Paradijs J. 1996, ApJ, 464, L139 [NASA ADS] [CrossRef] [Google Scholar]
  98. Veledina, A. 2016, ApJ, 832, 181 [NASA ADS] [CrossRef] [Google Scholar]
  99. Veledina, A., Poutanen, J., & Vurm, I. 2011a, ApJ, 737, L17 [NASA ADS] [CrossRef] [Google Scholar]
  100. Veledina, A., Vurm, I., & Poutanen, J. 2011b, MNRAS, 414, 3330 [NASA ADS] [CrossRef] [Google Scholar]
  101. Veledina, A., Poutanen, J., & Vurm, I. 2013, MNRAS, 430, 3196 [NASA ADS] [CrossRef] [Google Scholar]
  102. Veledina, A., Gandhi, P., Hynes, R., et al. 2017, MNRAS, 470, 48 [NASA ADS] [CrossRef] [Google Scholar]
  103. Vurm, I., & Poutanen, J. 2009, ApJ, 698, 293 [NASA ADS] [CrossRef] [Google Scholar]
  104. Walton, D. J., Mooley, K., King, A. L., et al. 2017, ApJ, 839, 110 [NASA ADS] [CrossRef] [Google Scholar]
  105. Wardziński, G. & Zdziarski, A. A. 2000, MNRAS, 314, 183 [NASA ADS] [CrossRef] [Google Scholar]
  106. Wardziński, G. & Zdziarski, A. A. 2001, MNRAS, 325, 963 [NASA ADS] [CrossRef] [Google Scholar]
  107. Wardziński, G., Zdziarski, A. A., Gierliński, M., et al. 2002, MNRAS, 337, 829 [NASA ADS] [CrossRef] [Google Scholar]
  108. Young, A. J., Fabian, A. C., Ross, R. R., & Tanaka, Y. 2001, MNRAS, 325, 1045 [NASA ADS] [CrossRef] [Google Scholar]
  109. Yuan, F., & Narayan, R. 2014, ARA&A, 52, 529 [NASA ADS] [CrossRef] [Google Scholar]
  110. Yuan, F., & Zdziarski, A. A. 2004, MNRAS, 354, 953 [NASA ADS] [CrossRef] [Google Scholar]
  111. Yuan, F., Zdziarski, A. A., Xue, Y., & Wu, X.-B. 2007, ApJ, 659, 541 [NASA ADS] [CrossRef] [Google Scholar]
  112. Zdziarski, A. A., & Gierliński, M. 2004, Prog. Theor. Phys. Suppl., 155, 99 [Google Scholar]
  113. Zdziarski, A. A., Johnson, W. N., & Magdziarz, P. 1996, MNRAS, 283, 193 [NASA ADS] [CrossRef] [Google Scholar]
  114. Zdziarski, A. A., Poutanen, J., Mikolajewska, J., et al. 1998, MNRAS, 301, 435 [NASA ADS] [CrossRef] [Google Scholar]
  115. Zdziarski, A. A., Lubiński, P., & Smith, D. A. 1999, MNRAS, 303, L11 [NASA ADS] [CrossRef] [Google Scholar]
  116. Zdziarski, A. A., Lubiński, P., Gilfanov, M., & Revnivtsev, M. 2003, MNRAS, 342, 355 [NASA ADS] [CrossRef] [Google Scholar]
  117. Zdziarski, A. A., Malyshev, D., Chernyakova, M., & Pooley, G. G. 2017, MNRAS, 471, 3657 [NASA ADS] [CrossRef] [Google Scholar]
  118. Zimmerman, E. R., Narayan, R., McClintock, J. E., & Miller, J. M. 2005, ApJ, 618, 832 [NASA ADS] [CrossRef] [Google Scholar]

2

We note that xillverCp table model contain errors, in particular, all spectra at one particular inclination of cos i = 0.65 are in error and we interpolated the spectra to this angle using the results for the neighbouring angles. This error may cause the fits with xillverCp or its relativistically smeared analog relxillCp to be stacked close to this specific inclination leading to unrealistically small errors on the inclination.

All Figures

thumbnail Fig. 1

Geometries of the central parts of the accretion flow proposed over time. (a) Static sandwiching corona, (b) static patchy corona, (c) outflowing sandwiching corona, (d) outflowing patchy corona, (e) cold accretion disc detached from the hot flow and (f) intersecting with it, (g) hot flow with cold clumps, (h) truncated disc and hot flow with substantial cyclo-synchrotron radiation.

In the text
thumbnail Fig. 2

Self-consistent spectra from a homogeneous slab corona of τT = 0.4 above a disc as observed at inclination i = 60°. The red and blue curves correspond to log ξ = 1 and 3, respectively. The dotted lines show the Comptonization continuum, the dashed lines show the reflection xillverCp spectrum together with reprocessed blackbody emission, and the solid lines give the total spectrum. The reflection spectra have been re-binned at a resolution Δ log E = 0.025. The bolometric luminosity of the reflection spectrum here is smaller than that of the Comptonization continuum because part of the reflected/reprocessed emission is Compton scattered in the corona.

In the text
thumbnail Fig. 3

Dependence of the 2–10 keV photon spectral index Γ of the Comptonized spectra from the sandwich corona (red filled symbols and left axis) and of the coronal temperature kTe (blue open symbols and right axis) on the ionization parameter for various coronal optical depths. Circles, squares, and triangles correspond to τT = 0.2, 0.4, and 0.8, respectively. The shaded area represents the typical observed range for Γ (1.6–1.8) and kTe (50–120 keV) (Wardziński et al. 2002; Zdziarski & Gierliński 2004; Ibragimov et al. 2005; Torii et al. 2011; Burke et al. 2017).

In the text
thumbnail Fig. 4

Self-consistent spectra from hot thermal plasmas mutually interacting with cold media without internal dissipation. The Thomson optical depth of the hot flow down to the central plane is τT = 0.4 and the spectra are observed at inclination of i = 60°. Red solid, blue dotted, green dashed, brown dot-dashed and pink triple-dot-dashed lines correspond to the covering fraction of the cold clumps of fcl = 1.0, 0.8, 0.6, 0.4 and 0.2, respectively. The ionization parameter here is log ξ = 1. All spectra correspond to the same bolometric luminosity.

In the text
thumbnail Fig. 5

Same as Fig. 3, but for the hot flow with cold clumps. Here τT = 0.4 and squares, circles, triangles, diamonds and stars correspond to cold clump covering factors fcl = 1.0, 0.8, 0.6, 0.4 and 0.2, respectively.

In the text
thumbnail Fig. 6

Same as Fig. 3, but the photon index and the equilibrium electron temperature of the hot phase are shown as a function of the ratio of intrinsic disc dissipation rate to the total luminosity. Squares show the case of log ξ = 1, fcl = 0.1; circles log ξ = 3, fcl = 0.1; and triangles log ξ = 4, fcl = 1.0. In all cases τT = 0.4.

In the text
thumbnail Fig. 7

Example of a synchrotron self-Compton spectrum from hot plasma with a hybrid electron distribution. The components due to reflection/reprocessing with the reflection fraction of 0.3 and the ionization parameter of log ξ = 1 and the blackbody of the temperature of 50 eV due to the reprocesses emission are shown together by the blue dashed curve. The escaping synchrotron spectrum is shown by the green dot-dashed line, synchrotron Comptonization spectrum is shown by the red dotted curve, and the total spectrum is given by the red solid curve.

In the text
thumbnail Fig. 8

Comparison of the self-consistent spectra from the sandwich corona for two reflection models, xillverCp (solid blue curves) and pexriv (dashed red curves), for two different ionization parameters (a) log ξ = 1 and (b) log ξ = 3. The reflected and reprocessed blackbody spectra are shown by thin lines and the total spectrum by the thick lines. The upper and lower curves correspond to two inclinations, i = 18° and 60°, respectively.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.