Testing wind as an explanation for the spin problem in the continuum-fitting method
The continuum-fitting method is one of the two most advanced methods of determining the black hole spin in accreting X-ray binary systems. There are, however, still some unresolved issues with the underlying disk models. One of them manifests as an apparent decrease in spin for increasing source luminosity. Here, we perform a few simple tests to establish whether outflows from the disk close to the inner radius can address this problem. We employ four different parametric models to describe the wind and compare these to the apparent decrease in spin with luminosity measured in the sources LMC X-3 and GRS 1915+105. Wind models in which parameters do not explicitly depend on the accretion rate cannot reproduce the spin measurements. Models with mass accretion rate dependent outflows, however, have spectra that emulate the observed ones. The assumption of a wind thus effectively removes the artifact of spin decrease. This solution is not unique; the same conclusion can be obtained with a truncated inner disk model. To distinguish among valid models, high resolution X-ray data and a realistic description of the Comptonization in the wind will be needed.
Subject headings:accretion, accretion disks, wind, black hole physics
The determination of the black hole spin is one of the most intriguing issues in astrophysics. Knowledge of the spin sheds light on the evolution of X-ray binary systems (XRBs) and active galactic nuclei (AGN) (King et al., 2008; Volonteri et al., 2008; Wang et al., 2009; Cao, 2010), and it allows to assess the role of the black hole angular momentum in the formation of relativistic jets (Blandford & Znajek, 1977; Wilms et al., 2001; Miller et al., 2002; Wu et al., 2011; Li & Cao, 2012; Narayan & McClintock, 2012; Gu et al., 2015). There are several ways to determine the spin of a black hole. The two most reliable methods to date are fitting the broad Fe K line from reflected emission (Brenneman & Reynolds, 2006; Miller et al., 2008; Liu et al., 2015) and fitting the continuum emission (Zhang et al., 1997; Gou et al., 2014; Steiner et al., 2014b). Further ways to estimate the black hole spin include X-ray timing, i.e. the modeling of high-frequency quasi-periodic oscillations (HFQPOs) (Stella et al., 1999; Motta et al., 2014; Pasham et al., 2014) and X-ray polarimetry of the inner disk region (Connors et al., 1980; Dovčiak et al., 2008; Li et al., 2009). A detailed discussion of these four methods can be found in the review by McClintock & Remillard (2009). In this paper we concentrate on the problems connected to the continuum-fitting method.
Spin measurements for a number of (mostly) Galactic XRB sources based on the continuum-fitting method are available in the literature, for instance for GRS 1915+105 (e.g. McClintock et al., 2006; Middleton et al., 2006) and LMC X-3 (e.g. Davis et al., 2006; Steiner et al., 2010), two sources that attain very high luminosities. The method was also successfully used in the case of the Weak Line Quasar SDSS J094533.99+100950.1 (Czerny et al., 2011), the radio-loud narrow-line Seyfert 1 galaxy RX J1633+4718 (Yuan et al., 2010), and several other intermediate redshift quasars (Capellupo et al., 2015).
Despite the overall success of the continuum-fitting method, one fundamental problem remains unsolved, namely the apparent dependence of the spin on the source luminosity. During an outburst, XRBs go through the high/soft state where their spectra are dominated by emission that peaks at soft X-ray energies. This peak can be modeled by a multi-color blackbody disk. Thermal spectra with a disk flux contribution above to the total flux are suitable for the analysis (Remillard & McClintock, 2006). The continuum-fitting method compares the observed shape of the thermal continuum spectrum to an accretion disk model spectrum. A sequence of such thermal state spectra have been modeled, e.g. for GRS 1915+105 (McClintock et al., 2006). There, the fits resulted in a black hole spin that was significantly decreasing with the modeled disk luminosity above , where erg/s.
A plausible explanation is that the assumptions underlying the disk model in McClintock et al. (2006) break down at higher Eddington ratios. The standard accretion model describes an optically thick and geometrically razor thin () gas disk that reaches all the way down to the innermost stable circular orbit (ISCO) and efficiently emits at each annulus (Shakura & Sunyaev, 1973; Novikov & Thorne, 1973). With increasing mass accretion rate the standard disk is unable to cool efficiently. As the disk overheats, it expands vertically to scale heights that are not consistent with the initial assumptions anymore and a significant fraction of energy is advected inward. The disk structure is then better described by the slim disk model (Abramowicz et al., 1988; Sa̧dowski et al., 2011), which has a radiative efficiency that decreases with increasing mass accretion rate. The corresponding analysis has been done, e.g. for LMC X-3 (Straub et al., 2011). The estimated black hole spin, however, still decreases with luminosity also for the advective slim disk.
There are two reasons which could explain this apparent decrease in the spin. The first one is the description of the disk atmosphere. Radiative transfer in the disk atmosphere is usually incorporated in accretion disk spectral models via the TLUSTY code (Hubeny & Lanz, 1995) which locally computes absorption opacities of various chemical elements in plane-parallel layers. In the X-ray spectral fitting software XSPEC (Arnaud, 1996) a well known and widely used thin disk model with full radiative transfer is bhspec (Davis & Hubeny, 2006), while slim disks with full radiative transfer are available as slimbh111http://astro.cas.cz/slimbh (Sa̧dowski et al., 2011; Straub et al., 2011). However, the physics of the disk atmosphere is very complicated and largely unknown, with magnetic fields, turbulent stress, and dissipation taking place (see e.g Blaes et al., 2011; Tao & Blaes, 2013). This manifests in the fact that both bhspec and slimbh seem to overestimate the amount of photons at the spectral peak (see Fig. 4 in Straub et al., 2013). The second reason is mass-loss due to winds/outflows from the disk surface. In general, outflows are thought to become significant when the mass accretion rate through the disk increases as discussed in a number of papers (King & Pounds, 2003; Slone & Netzer, 2012; Yuan et al., 2012; Cao, 2014; Laor & Davis, 2014; Yuan et al., 2015; Cao, 2016). Outflowing winds have recently been detected in several XRBs and AGNs via blueshifted X-ray absorption lines which are observed in high resolution spectra, for instance in GRS 1915+105 (e.g. Neilsen & Lee, 2009), GRO J165540 (e.g. Miller et al., 2006a; Kallman et al., 2009), IGR J17091-3624 (e.g. King et al., 2012), NGC 3783 (e.g. Kaspi et al., 2002), and PG 1211+143 (e.g. Pounds et al., 2003).
In the present paper we perform simple estimates to check whether a wind scenario is likely to be responsible for the observed apparent spin decrease over luminosity. Given the fact that spin measurements from both thin and slim disks show the same type of deviation we use the mathematically simpler standard Novikov & Thorne (1973) model to describe the underlying accretion disk. We then modify the mass accretion rate to incorporate mass loss due to an unspecified wind and calculate the energy spectra of such flows. In Sec. 2 and 3 we present an estimate of the energy that is driven away by wind and possible disk models with wind. Sec. 4 shows Comptonization of the disk spectrum by the wind. In the extended discussion in Sec. 6 we test our models against observational data of the two most explored sources LMC X-3 and GRS 1915+105. We give our conclusions in Sec. 7.
2. Analytical estimates
The apparent decrease of the black hole spin with increasing source luminosity reported in GRS 1915+105 and LMC X-3 is significant in both cases (McClintock et al., 2006; Steiner et al., 2010; Straub et al., 2011). For this reason the disk is usually assumed to be reliably optically thick and geometrically thin only for modeled disk luminosities in the range . In the relativistic standard model (Novikov & Thorne, 1973), the apparent decline of spin with luminosity can be directly translated into a decreasing accretion efficiency, . We now postulate that (i) this decrement in spin roots in the existence of wind and that (ii) the amount of energy driven away by this wind corresponds to the drop in accretion efficiency. The efficiency to convert mass into radiation is given by
where is the co-variant time component of the four-velocity of the gas. For convenience, we show the theoretical radiation efficiency for a relativistic thin disk in Fig. 1. A more realistic modeling of how efficiency depends on spin one can be found in Li et al. (2005). Wind then has to carry away the fraction of the total accretion energy
where is the asymptotic flow efficiency at low luminosity.
Recent optical data lead to an improvement of the dynamical parameters of LMC X-3 (Orosz et al., 2014). The black hole mass, and the inclination angle, , result in a black hole spin, (Steiner et al., 2014b). In the framework of our forthcoming paper II (Straub et al., in preparation) we estimate the black hole spin in LMC X-3 based on the improved parameters using the slim disk model, following the methodology of Straub et al. (2011). We find for a mean spin value that trails off to very quickly beyond which is consistent with previous work. Given the flat theoretical profile of at small spins, the new values imply .
Recent measurements of the GRS 1915+105 parameters yield a black hole mass of , an inclination of , and a source distance of about 8.6 kpc (Reid et al., 2014). These values entail a spin of (Fragos & McClintock, 2015) which falls off to (assuming , Straub et al., in preparation). This translates to an asymptotic efficiency that decreases to resulting in .
We adopt throughout the paper the above quoted dynamical parameters for LMC X-3 and GRS 1915+105 and the respective mean spin values which we found using the slim disk model.
3. Simple wind models
In this section we introduce a simple parametric description of the wind. First, we concentrate on the effect the wind has on the intrinsic shape of the thermal emission from an accretion disk. We describe the wind as a radially dependent quantity that decreases the accretion rate towards the black hole. In the following calculations we adopt the code by Czerny et al. (2011) which is based on Novikov & Thorne (1973) and neglects self-irradiation of the disk as well as Comptonization in the disk atmosphere. Moreover, we assume that the disk thickness is small so that the photons are emitted from the equatorial plane (Sa̧dowski, 2009). Therefore, the local disk emissivity is independent of any assumptions about the disk’s vertical structure, including viscosity. We assume zero torque at the ISCO.
The wind from the accretion disk surface that is required to cause the apparent change in the measured spin must have the appropriate radial distribution. Since the observed X-ray emission comes predominantly from the innermost tens of gravitational radii, the type of wind that can solve the spin drop problem must be localized there.
There are a number of ways to describe the mass loss due to wind. We define and as the accretion rates at the inner and outer disk radius, and , respectively. and are the inflow at a given radius, , and total mass loss beyond this radius, respectively. We consider several wind models and calculate their spectra for a range of accretion rates and fixed spin parameter . Then, we fit the emergent spectra using the relativistic standard model (Novikov & Thorne, 1973) without any wind, leaving accretion rate and spin as free parameters. This procedure allows to verify whether the measured spin is lower than the input spin in the disk-wind model as it is in the fits to observational data, i.e. if .
We define the Eddington limit on the mass accretion rate as g/s, where () for GRS 1915+105 and g/s, where () for LMC X-3. With the above definitions of , we have assuming the efficiency is constant. In the following four wind basic models we assume and , where . The corresponding inflow and mass loss profiles as a function of are shown in Fig. 2 for , while Fig. 3 illustrates the respective spectra.
In the first wind model, we simply assume that the radial distribution of the accretion rate has the form of an exponential function,
where and are two free variables satisfying and . If there is no wind, which implies . Figure 2 shows the profile for , , , and total mass loss rate .
In the third model we adopt the description of the wind suggested by Blandford & Payne (1982),
The continuity equation reads
where is a free parameter (see Fig. 2 for , , and total mass loss rate ).
Finally, in the fourth model we release the assumption that the disk extends to the ISCO. Instead, we assume that the disk is truncated at due to the wind and that the accretion rate is constant at all disk radii. For simplicity we assume that depends linearly on accretion rate
where is in units of . and are two arbitrary free variables (see Fig. 2 for inner truncation at , and the total mass loss rate ). This model may seem particularly unphysical, but it is useful as a parametrization. Its plausible physical interpretation is later discussed in Sec. 6.
We now try to recreate the measured accretion rate dependency of the black hole spin by fitting disk spectra based on the four wind models above with a standard wind-less thin disk model (Novikov & Thorne, 1973). This part of analysis does not involve any real spectral data but is based on comparison of models with and without wind. We use the Fortran ray-tracing code described in Czerny et al. (2011) which includes all general relativistic effects. The local emissivity follows the Novikov-Thorne prescription but the local accretion rate decreases inwards with radius in the presence of winds. The disk spectra are parametrized by input mass accretion rate and spin (, ) for an assumed viewing angle, , plus the parameters of the individual wind model. We calculate the resulting spectra including light bending effects. In the next step we treat these synthetic spectra as “observed spectra” and apply a maximum likelihood method to fit them (still using Fortran) with the original wind-less Novikov-Thorne prescription for the local flux. We compare the two models in the range 0.1-100.0 keV. The accretion rate and the spin are now the fit parameters. These are plotted in Fig. 4 and Fig. 5 for GRS 1915+105 and LMC X-3, respectively. In both figures the red dashed line indicates the best linear fit to the bulk of spin measurements obtained from fitting the slim disk model to the actual source spectra. The question is whether the derived sets of (, ) from fitting the synthetic spectra follow these reference lines from fits to real spectra. There are no observational errors or instrument sensitivity effects, so the absolute values of are meaningless but the procedure allows to find the best match between the models with and without wind.
The radial mass loss profiles presented in Fig. 2 imply that Model (2) shows the steepest dependence on the radius. However, the inflow rate only significantly decreases in a small portion of the disk close to ISCO and has little influence on the shape of the model spectrum compared to other models (see Fig. 3). The radial change of the accretion rate in Model (1) is in general steeper than that in Model (3) at the inner radius. This means that Model (1) will affect the disk spectrum, especially the contribution from the inner disk radius, more significantly which makes it more interesting for our case. However, extending the analysis to spectra from to we notice that Model (1) overproduces the mass loss rate which results in much lower black hole spins (green crosses in Fig. 4) compared to the measured ones, since the character of the current radial change of the accretion rate is not coupled strongly enough with the accretion rate. The - dependent Model (4) is the simplest possible model that is able to reproduce the observed apparent spin decrease. In the parameter space (D, E) when simulating the decline of spin with mass accretion rate, and give the best fits (the blue circles in Fig. 4) to the red dashed line. The uncertainties are quoted at the 90% confidence level by assuming an error of in the spin of the red dashed line. Interestingly, we can achieve an equally good match with Model (1), by assuming that parameter depends linearly on the mass accretion rate, i.e. Model (1*), . On the one hand, this means that when no winds are present, the exponential function in Eq. 3 must equal 1 and thus . On the other hand, an increasing reflects an increasing contribution of winds that reduces the local mass accretion rate and modifies the spectrum. This works because a smaller at low mass accretion rates gives a shallower shape to the radial distribution of the accretion rate. In this manner the - dependent Model () reduces the influence of the wind on the Novikov-Thorne spectrum at low accretion rates.
When simulating the behavior of versus for GRS 1915+105 in the range with Model () we explore the parameter space (). Parameter B is arbitrarily fixed at B = 0.15 since the radial distribution of the accretion rate is insensitive to it. We find that and give the best fits to the observed spin trend (see the green stars on the red dashed line in Fig. 4). We repeat the same analysis for LMC X-3 in the range of . We find again that Model () can describe very well the observed spin trend with . The reconstructed values of (, ) are plotted in Fig. 5 as green stars. The use of fixed values A1 and A2 works only in a limited range of accretion rates but this simple exercise shows that the requested wind is highly concentrated towards the ISCO: in GRS 1915+105 more than 50 per cent of wind mass loss takes place below 5 and in LMC X-3 most of the wind occurs inside 14 . Thus, in both cases the radial range for the apparent wind is only about three times the size of the ISCO.
From a cross-correlation analysis of optical/infrared (OIR) data and contemporaneous X-ray data, Steiner et al. (2014a) found that there is a non-linear relation between the OIR flux and the time-lagged X-ray emission: with . This may indicate that the rates of inflow at the inner and outer radii are not matched and that the ratio of the inflow rate scales inversely with the OIR flux, i.e. . Disk winds are suggested to carry the accreted gas away from the disk. Given the above - dependence on the parameter and Eq. 3, our Model () predicts that the ratio of , which is fully consistent with the anti-correlation between the inflow rate and found by Steiner et al. (2014a).
Our simple parametric discussion above shows that Model (), Model (4), or possibly more sophisticated models not considered here can describe such a mass loss and successfully reproduce the observed decrease in spin over a range of mass accretion rates. The removed material will be present between the disk and the observer and may contribute to additional modifications to the observed spectrum.
4. Physics of wind and Comptonization of the disk spectrum by the surrounding material
Our simple parametric description of the thermal emission from an accretion disk with mass losses used in Sec. 3 is too simplistic to predict the resulting effects of the removed material on the observed spectra.
The most general description of a 1-D disk structure with the wind should include both the mass loss from the disk, implying the dependence of the accretion rate on the radius, , as well as the energy extraction by the wind. Thus, in general, the local radiation flux, , in the disk is given by
if the particles in the wind carry more energy than the remaining particles in the disk. Here (r) describes the standard term for relativistic effects in the Novikov-Thorne disk, and the energy fraction accounts for additional removal of the energy by the wind. We need to specify two arbitrary radial functions, and , to fully account for the impact of the wind on the disk structure. The relative importance of the two effects likely depends on the wind driving mechanism.
There are a number of wind mechanisms proposed which may act on a cold accretion disk surface: (i) thermally driven winds (Parker, 1958; Balsara & Krolik, 1993; Woods et al., 1996); (ii) magnetically driven winds (Blandford & Payne, 1982; Konigl & Kartje, 1994; Contopoulos & Lovelace, 1994); (iii) radiation pressure driven winds, which in turn include Compton-driven and line-driven winds (Shlosman et al., 1985; Proga et al., 2000; Higginbottom et al., 2014). Thermally driven winds and radiation pressure driven winds do not require considerable energy segregation at the launch radius and in this case the role of the factor can be possibly neglected although outflowing plasma needs extra energy to achieve escape velocity. In a Newtonian Keplerian motion, it requires doubling the test particle energy to unbind it. Thus the effect is small as long as the ratio of the total wind mass flux and the accreted mass flux is not too high. On the other hand, most energy and angular momentum can be efficiently extracted from the accretion flow by a small fraction of the magnetically driven plasma (winds or jets).
In our description of the disk in Sec. 3 we assumed for simplicity. If the wind is actually a thermally driven wind, our expression for the mass loss indeed implies the true mass loss from the system. If the magnetic mechanism dominates, our simplistic method gives an upper limit on the true mass loss from observational constraints.
In any of these cases the material may, or may not, be present along the line of sight since this will depend of the wind geometry for a given inclination of an observer. If the wind is mostly collimated, the line of sight may miss most of the removed material, but if the wind is roughly spherical (or, accidentally, mostly towards an observer) we should see signatures of this wind in the form of Comptonization effect and/or absorption lines.
An example of a continuum radiation driven wind has been discussed in detail by King & Pounds (2003) *[also see]2016MNRAS.455.1211K. Such an outflow is expected in sources radiating close or above the Eddington limit. If most of the material does not reach the ISCO but flows out, it may form a roughly spherical outflow or a collimated wind in a double cone. The analytic solutions for both geometries showed that the outflow may be Compton-thick and is a promising candidate for the soft excess observed in many AGNs and ultraluminous X-ray sources, given that much of the accretion energy must emerge as blackbody emission due to its large optical depth. The parameters of such an outflow were discussed by King & Pounds (2003). If the wind is spherical, its optical depth can be described as
where is the total mass outflow rate, and is the wind velocity which is higher than the escape velocity from a given radius. For sources that accrete close to the Eddington limit this optical depth is relatively large. However, the temperature, , of such an outflow is not very high since there is not enough additional energy to heat the plasma. So, the temperature of the outflow stays close to the disk temperature. Another possibility which is also discussed in King & Pounds (2003) is that the outflow is highly collimated and of a jet-like form.
In Fig. 6 we show the predicted disk spectrum seen through such a cold optically thick emerging wind for GRS 1915+105 predicted by Model (, red long-dashed line), and compare it to the initial Model () spectrum (black solid line). A cold wind softens the low-energy part of the spectrum and hardens the spectral peak and the high-energy part, and does not produce a hard X-ray power law tail. Since the wind density is relatively high, we have to estimate the possible role of the bremsstrahlung radiation. If the outflow consists of fully ionized pure hydrogen gas, the luminosity of the bremsstrahlung radiation is
where the density of the wind is assumed to be , and then the optical depth . For keV, , , and (at ) one gets erg/s, while the luminosity of the outflow spectra (black solid line) is erg/s. Therefore, the emergent spectra from the outflow are dominated by the disk emission and the Comptonization, and the model neglecting the bremsstrahlung is self-consistent. Thus the bremsstrahlung radiation will not be taken into account in further considerations.
On the other hand, the emerging wind may be magnetic in nature, or a weak static magnetic corona can be present above the accretion disk, as frequently postulated for soft states (Done et al., 2007; Cao, 2009; You et al., 2012). We then have no particular predictions of the properties of the Comptonizing medium from our model. Therefore, we arbitrarily assume the values of keV, , which implies a small fraction of energy in this component (Comptonization parameter ). We show this example for the same disk model in Fig. 6 (blue short-dashed line). A hot Comptonizing medium marginally hardens the spectral peak and produces a hard high-energy tail to the thermal component which is in general rather steep due to the efficient supply of the soft photon from the disk and the negligence of the outflow velocity.
In order to show the magnitude of the effect of Comptonization inside the wind more clearly, we plot the renormalized and shifted spectra in a linear scale, with peaks at the same energy, in Fig. 7. We can now compare them at the energy where the initial Model () disk flux is a factor 2 below its peak flux. At this location the differences between the models should show more clearly in the data, since at higher energies the data usually requires additional hard X-ray component. We see that the spectrum where photons go through a cold optically thick wind is visibly steeper with fewer photons at the given energy, while the spectrum where photons go through a hot optically thin wind is less steep with slightly elevated flux. For comparison, we also added the Novikov-Thorne model without any form of outflow (magenta dotted line). This spectrum is also slightly harder than Model () and the cold wind case. The difference between the two disk models with/without outflows is at the level of , again at the energy where the flux of the disk before Comptonization is a factor 2 below the peak (see Fig. 7). However, the difference between the disk without wind and the disk with wind and optically thin Comptonization (i.e. between magenta dotted line and the blue short-dashed line) is very small, only , almost invisible in the figure scale. Observational data of good quality could thus allow us to discriminate (within the observational error) between the different spectral shapes shown in Fig. 7 and enable us to assess not only whether wind resolves the spin problem in accretion disks, but also which of the outflow/wind scenarios presented above is the preferred one.
Simple models discussed above are order of magnitude estimates and do not rely on detailed physical studies; they consider only the simplest geometry which almost certainly differs from real accreting systems. Recent modeling of hard X-ray emission in soft state XRBs and AGN with Nuclear Spectroscopic Telescope Array (NuSTAR) indicates that the emission comes from a very compact region, 3 to 10 in size, that is well approximated by a lamp-post located on the symmetry axis (e.g. Fabian et al., 2015). This region can be identified as the jet base. Although the jet seems to be suppressed in the high/soft state in XRBs (Fender et al., 1999, 2004), the jet base may remain present and be related to uncollimated winds and significant energy dissipation near the black hole. The material presumably comes from an accreted hot corona as it has been discussed in a number of papers over the past decade (e.g. Chakrabarti & Titarchuk, 1995; Zycki et al., 1995; Liu et al., 2015; Wilkins & Gallo, 2015). Since winds due to radiative or magnetic acceleration from the inner disk region must be Compton-thin as demonstrated by Reynolds (2012), the general scenario of a cool, optically thick disk wrapped in a hot, optically thin corona, where a fraction of the coronal material is removed in form of winds while the rest is accreted, applies. We sketch this in Fig. 8 (see the discussion in the following section). A detailed computation of the emission from such a complex model is beyond the scope of the present paper as this would require knowledge of the dissipation in the coronal hot flow. However, we can verify if the simple wind model discussed in Sec. 3 is self-consistent and determine which part of the flow is actually seen. We present such observational tests in the following Sec. 5 for LMC X-3 and GRS 1915+105 using Model ().
5. Testing the wind scenario in LMC X-3 and GRS 1915+105
In the previous section it was shown that in principle one can test the wind scenario if the observational data determine not just the position of the spectral peak and the normalization but if they are also sensitive to the precise shape of the thermal component. In this section we fit Model () against observational data of LMC X-3 and GRS 1915+105. The spectral analysis on all data sets is performed with the X-ray spectral fitting software package, XSPEC v.12.8.2. In order to implement our disk-wind model into the framework of XSPEC we construct a table model diskw for which each individual spectrum is calculated with a modified version of the code by Czerny et al. (2011) based on Novikov & Thorne (1973), which accounts for the modified accretion rate given in Model (). The description of the hardening factor adopted in the model spectrum is done through the parametrization of bhspec results for the viscosity parameter
where and K. It partially comes from Done et al. (2012) and Eq. (A13) in Davis et al. (2006), but includes the explicit term in accretion rate which is absent in those papers. The hardening factor increases monotonically with temperature. Therefore, in our table model, there are only two free parameters, parameter (see Eq. 3) and the initial mass accretion rate at the outer disk edge, i.e. before the disk is affected by mass loss due to the wind. The mass accretion rate, , is scaled by the Eddington limit of a given source, (see Sec. 3 for the definition of ). Unlike in the initial Model () in Sec. 3, the free input parameters and are not coupled in the fitting here. Parameter is fixed at , since the radial distribution of the accretion rate is insensitive to it. For each source we use the respective up-to-date binary parameters given in Sec. 2 and the resulting mean spin values derived from the spectra obtained in low disk luminosities, , in fits with the slim disk model (Straub et al., in preparation). Our base model to fit the soft state spectra of both sources reads
tbabs * (diskw + nthcomp).
The component tbabs (Wilms et al., 2000) accounts for photon absorption by neutral hydrogen in the direction of the source. The column density towards LMC X-3 is cm (Page et al., 2003) and towards GRS 1915+105 it is cm (Lee et al., 2002). The disk component is normalized by a constant factor 0.0575 which converts the units between the table model and the energy bins. Both sources require an additional component to account for the high energy photons. We use the thermal Comptonization model nthcomp (Zdziarski et al., 1996; Życki et al., 1999). The photon counts at high spectral energies have fairly large uncertainties, so in order to avoid degeneracies with other parameters the photon index which defines the powerlaw slope is allowed to vary only in the range . The electron temperature, and the normalization are left free, whereas the seed photon temperature, , is fixed. We get lower limit on by first running fits with the non-relativistic disk model diskbb (Mitsuda et al., 1984) which are found to be reasonable values as for the temperatures of the disk with wind, and those inner disk temperatures are used as fixed seed photon temperatures in nthcomp.
In the spectral fits with diskw the vertical structure of the disk is not taken into account. In particular the disc thickness is assumed to be small so that the photons are emitted from the equatorial plane, which is most likely inappropriate for high mass accretion rates.
|source||Obs 1||Obs 2||Obs 3||Obs 4|
5.1. Observational tests with LMC X-3
We look at the sample of 712 X-ray spectra that have been obtained using the large-area Proportional Counter Array (PCA) aboard the Rossi X-ray Timing Explorer (RXTE). These data were extracted from the Proportional Counter Unit 2 (PCU-2), background subtracted and corrected for a detector dead time of , and classified as the soft states (a detailed description is given in Steiner et al., 2010, 2014b). We fit these soft state spectra with the above given base model tbabs * (diskw + nthcomp) over the energy range 2.51-25.0 keV. A good fit satisfies (i) and (ii) the unabsorbed disks flux must be at least of the total unabsorbed flux. Out of the total number of soft state observation 381 spectra have mass accretion rates in our range of interest , where winds are expected to play a role. We obtain 336/381 good fits. Parameter decreases monotonically with mass accretion rate from about 4.5-0.03 while decreases from 1.5-0.15 (see Fig. 11). The behavior of parameter suggests that the wind in this source decreases markedly, but never entirely dies down. The estimated seed photon temperatures take values keV and the electron temperatures lie between keV. From the goodness of fit and the residuals we see that a disk-wind model based on a fixed spin, , represents the LMC X-3 data well. We choose four representative spectra in a wide range of mass accretion rates and show their data, the models, and the residuals in Fig. 9. The results are summarized in Table 2.
|model||Obs 1||Obs 2||Obs 3||Obs 4|
|/dof ()||50.24/45 (1.12)||45.03/44 (1.02)||42.68/49 (0.87)||21.32/44 (0.49)|
Note. – The best fitting spectral parameters. All errors are quoted at the confidence level (). The asterisk, , indicates an unconstrained upper or lower limit. The blackbody temperature, , has been estimated with diskbb and was subsequently kept frozen during the fit.
The best fits for Observations 2 and 3 give strong constraints on the Comptonizing medium since the electron temperature in nthcomp is well measured. As an example, we analyze in detail the solution obtained for Observation 2. nthcomp uses the slope, , as the second parameter but for a given electron temperature, the optical depth is uniquely defined for a fixed geometry. Using the formulae given in Czerny & Zbyszewska (1991) for a point-like source of photons and a spherical geometry of the Comptonizing medium (the same as adopted in nthcomp), we obtain the optical depth and the Compton parameter . The medium is therefore optically thick, but the Comptonization is still unsaturated due to the low value of the electron temperature. In the data, however, the normalization of the Compton component is much lower (below 10 %) than 63 % expected from . This reveals that the material responsible for the fitted Comptonization does not cover the whole inner disk. In addition, the wind in this solution is strong; 84 % of the material is removed from the disk at the ISCO.
On the one hand, the maximum of the local disk temperature is at 10 where the mass fraction removed from the disk is 40%. On the other hand, half of the total flux is emitted at 40 where the removed mass fraction is 22%. If the hot, coronal material is predominantly inflowing (as a sub-Keplerian flow), then the second value is more representative for the overall Comptonization of the disk spectrum. From this we can estimate the vertical optical depth of the inflow at a given radius using
where we only have to assume an inflow velocity, , as both the corona density, , and the corona thickness, , cancel out. Assuming a velocity of 0.1 c, one then obtains at 40 an optical depth of only . This type of material would cover the whole inner disk. Given the low normalization of the Compton component in the data, however, it cannot be the material responsible for the fitted Comptonization. The optically thick Comptonizing medium seen in the data should rather be identified with a central hard X-ray source (see Fig. 8). In Observation 3 the optical depth, , of the medium is somewhat lower than in Observation 2, but otherwise the solution is very similar.
5.2. Observational tests with GRS 1915+105
RXTE made pointed observations of GRS 1915+105 from April 1996 to April 2009. We use the steady-soft observations of Peris et al. (2015). These energy spectra were extracted from the data of the top layer of PCU-2 and the background was subtracted using the model applicable for bright sources. A systematic error of 1% has been added for each energy channel. From the total set of over 2000 continuous exposure segments (mean exposure time = 2.1 ks) 1257 spectra are steady and of these 264 are soft. The majority of these show absorption and/or emission as well as disk reflection features. In order to fit the soft state spectra of GRS 1915+105 we therefore include a broad Fe K emission line fixed at keV and modeled with laor (Laor, 1991) as well as a Fe absorption component (gabs) with a central energy fitted between 6.5-7.5 keV and fixed width keV. We use smedge to describe a smeared Fe edge where we fit the edge energy between 7 - 9 keV, leave the maximum absorption factor free, and fix the smearing width at keV. We fit the full model tbabs * (diskw + nthcomp + laor) * gabs * smedge in the energy range 2.51-45.0 keV. 233 spectra have luminosities in our range of interest () and we obtain 178/233 good fits. The modeled mass accretion rates lie between over the range of which parameter increases from near zero to 2.5 (see Fig. 11). We note in particular that becomes infinitesimal already at around of the Eddington limit. Below, the disk is effectively a standard thin disk without winds. The electron temperature slightly increases with mass accretion rate from keV together with the seed photon temperature keV. The latter stays roughly constant for and is substantially higher than in other Galactic XRBs but typical for the source. The absorption line energy which lies at about keV is practically independent of the accretion rate. Its strength, however, decreases with increasing mass accretion rate. We remind the reader that in gabs the line strength, , refers to the line depth and is related to the optical depth at the line center via . The anti-correlation between absorption line strength and mass accretion rate implies that an increasing amount of wind can effectively remove absorption line features. The laor normalization which parametrizes the number of photons per area per time increases with mass accretion rate. The Fe K line is present in all spectra but more pronounced when wind is present. The central edge energy shows an increasing trend with mass accretion rate, keV. Although the goodness of fit is already good with without Fe line component, there are still small residuals which seem barely visible in some of the spectra. Adding the laor model component to account for the Fe line is in fact improving the fits. The disk-wind model based on a fixed spin, , can fit the GRS 1915+105 data very well. We choose again four representative spectra at widely different mass accretion rates and show the data, models, and model residuals in Fig. 10. The results are summarized in Table 3.
The best fit found for Observation 1 is consistent with recent results based on data recorded with NuSTAR which operates in a wider bandpass (3 - 79 keV) than RXTE: Miller et al. (2013) measure a slope and an electron temperature 16 keV (see their Fig. 2) while we have only a lower limit for the electron temperature and a very simplified description of reflection due to the lack of high energy data in RXTE, which imply an optical depth of for the geometry adopted in nthcomp ( for Observation 1 in this work); and a similar value of the optical depth () is obtained by Miller et al. (2013) using the eqpair model (Coppi, 1999), where the optical depth is a free parameter. The derived Compton parameter for the NuSTAR data, ( for Observation 1 in this work), is relatively high but, like in the case of LMC X-3, the Compton component contributes less than 10% to the total spectrum in the data; most of the emission stems from the disk component. This means that in Observation 1 the Comptonizing medium covers only a small part of the disk. Based on our model, the wind contribution during Observation 1 is extreme with over 90 % of the material being removed from the inner disk. On the one hand, half of the material has left the disk by the time the disk flow reaches 6.4 . On the other hand, half of the flux is emitted at 14.3 , where the disk has lost only 26 % of its matter to the wind. The optical depth at the latter radius is expected to be for an adopted corona flow speed , where the high corona speed reflects the fact that the inner disk region in a high spin source like GRS 1915+105 is closer to the black hole than in a low spin source like LMC X-3. A spherical wind would have to have a very high optical depth (see Eq. 11), therefore, the self-consistent picture is the following: The Comptonized emission seen in Observation 1 originates in fact from a central hard X-ray source and the wind material removed from the inner disk forms an inflowing coronal layer that is accreted onto the black hole. The emission from the optically thin corona is in this scenario not directly visible in the spectral data.
The situation during Observation 4 is a quite different one. Firstly, the mass accretion rate is small and there is no significant presence of wind. Secondly, we obtain good limits on the temperature of the Comptonizing plasma since it is sufficiently low for the data to be properly recorded within the energy bandpass of RXTE. We derive the optical depth from the measured electron temperature, keV, and the lower limit of the photon index, , for the geometry adopted in nthcomp. The corresponding Compton parameter is then only , much lower than in Observation 1, and indicates only a week Comptonization. The normalization of the Compton component is much higher, though, suggesting that a significant part of the inner disk is covered. The upper limit on the wind parameter implies an upper limit on the total wind mass loss of only 20 % while merely 4 % of the disk mass is liberated at the location () where half of the flux is dissipated. This small amount of material is not enough to form an optically thick Comptonizing zone, in particular if the material flows in as fast, as we assumed before. The coronal optical depth estimated at this radius is which is much less than what is measured for the spherical Comptonizing medium. It may imply, again, the presence of a central hard X-ray source. Given the high normalization during Observation 4, however, the result rather suggests that a significant part of the Comptonization may take place in a plane-parallel optically thin corona layer. To obtain such a geometry, the hot, coronal inflow must be more than order of magnitude slower than the free fall velocity.
|model||Obs 1||Obs 2||Obs 3||Obs 4|
|/dof ()||43.23/70 (0.62)||48.97/70 (0.7)||44.53/69 (0.65)||63.47/62 (1.02)|
Note. – The best fitting spectral parameters. All errors are quoted at the confidence level. The asterisk, , indicates an unconstrained upper or lower limit. The blackbody temperature, , in nthcomp has been estimated with diskbb and was subsequently kept frozen during the fit.
We fitted the bulk of our thermal spectra of LMC X-3 and GRS 1915+105 in a broad luminosity range with a simple disk-wind model that includes an inward radially decreasing accretion rate but neglects the problem of additional energy extraction from the disk by the wind. In data fitting we do not use any coupling between the amount of the removed material and the properties of the Comptonizing medium (i.e. the possible geometry, the temperature and optical depth). The model fits the data well and there is no effect of decreasing spin with luminosity since the spin is fixed by definition. The required wind is negligible at low mass accretion rates, particularly in GRS 1915+105. For increasing (outer) mass accretion rates as much as 80% of the material can be lost in the wind (see Tab. 2, Tab. 3, and Eq. 3).
Both sources show significant Comptonization effects in the spectra. The Comptonization is consistent with the upscattering of the soft disk photons by the hot medium since the soft photon input temperature is comparable to the disk temperature. The temperature of the Comptonizing medium is higher than the maximum disk temperature by a factor 2 to 10, and in general becomes higher for higher luminosities. The power law slopes measured for GRS 1915+105 (see Tab. 3) range from = 1.5 to 3.5 and are often not determined accurately due to the limited energy band of RXTE. The well determined slopes combined with the corresponding plasma temperatures yield the optical depths of the Comptonizing plasma of the order of = 2 - 3. Such a plasma phase is typical for many soft state sources (Walter & Fink, 1993; Magdziarz et al., 1998; Zhang et al., 2000; Kubota & Done, 2004; Done et al., 2012; Petrucci et al., 2013). A similar result was obtained by Miller et al. (2013) with NuSTAR data, where the slope of the Comptonized component is determined more reliably than with our data. For LMC X-3 we arrive at the same conclusion: whenever is well constrained, the optical depth of the Comptonizing medium is of order of 3 - 4. The exact values depend on the geometry (plane-parallel or spherical) and the assumed photon distribution. We are using here a spherical symmetry, with photons at the center like in nthcomp, for consistency.
The problem is that the optical depth of the Comptonizing medium implied by the fits does not increase with the modeled disk luminosity but remains practically constant despite the sharp rise in the wind from the disk itself. Thus, the Comptonizing medium is present at low Eddington ratios, when mass loss from the disk is not required as well as at high Eddington ratios, when the implied disk mass loss is strong. Moreover, the medium does not become optically thicker as the mass accretion rate rises and winds become more prominent. This seems to suggest that the wind material is not located along our line of sight. The behavior of the normalization of the Compton component points towards an absence of Comptonizing material along the line of sight, too. Given the Compton parameter (LMC X-3), for instance, the Comptonization is strong and the material should be optically thick. If this kind of material would cover large part of the disk, the disk itself would be barely visible as almost all its energy would end up in the Compton component. This contradicts the nature of the data where the Compton component makes up only a very small fraction of the total flux and the disk is seen well, including the reflected component. Therefore, the material that has been removed from the disk by the wind cannot be fully identical to the Comptonizing medium. In general, this medium likely forms a kind of magnetized dissipative skin (e.g. Różańska et al., 2015), a small fraction of it may be outflowing, but most of it would be inflowing. We sketch this geometry in Fig. 8. If the hard emission comes from a lamp-post style hard X-ray source, the question remains whether we have any method to prove or disprove the mass outflow from the disk. Our consistency check made for LMC X-3 and GRS 1915+105 implies that the actual wind material may be unnoticed since the data do not require a second Comptonizing medium and the model does not overpredict the optical depth of the corona skin at the radius where half of the disk flux is emitted.
In the high spin source GRS 1915+105 we detect an absorption feature along the line of sight but we do not see it in the low spin source LMC X-3 despite their similar inclination angles, . If this absorption line at 7.2 keV is caused by a hot plasma it must be blueshifted from its original value of 6.7 keV. This implies that the absorbing (wind) material has a speed of 0.07 c. Such a wind may originate from the innermost part of the disk but does not have to. Hydrogen-like and helium-like iron lines are resonant lines and the line-locking mechanism can accelerate the flow very efficiently independent of the launching radius. In comparison, broad absorption line (BAL) features observed in AGN have blueshifts that correspond to velocities of the order of 0.1 - 0.3 c, but the launching radius is likely further than 0.1 pc from the black hole. Recent studies of a large sample of AGN suggest typical launching radii of 10 to 1000 light days and some features can come from the distances as far as 100 - 3000 pc (Filiz Ak et al., 2013). The lack of line absorption in LMC X-3 might be due to its low spin. Steiner et al. (2014b) tentatively suggested that there is a possible link between black hole spin and spectral complexity. They argued that some sources with known low spins have remarkably simple spectra with few spectral features whereas high spin sources exhibit strong Comptonisation/reflection and a rather large variety of spectral features.
Previous models of the outflows in GRS 1915+105 assumed magnetically driven outflows to reproduce the amplitude of the heartbeat regular bursts within the scenario of the radiation pressure instability (Nayakshin et al., 2000; Janiuk et al., 2002). However, models with a static corona and a time-dependent mass exchange between the disk and the corona (e.g. Janiuk & Czerny, 2005) or models with different viscosity prescription (e.g. Merloni & Nayakshin, 2006) worked equally well.
Model (4) which we used in Sec. 3 as one of the simple parametric prescription implicitly refers to a magnetic outflow. In this model, formally, all the matter is outflowing at the disk truncation radius. From an observational point of view, by fitting the broad asymmetrical Fe K line, Cowperthwaite & Reynolds (2012) find that the accretion disk in 3C 120 might be truncated at independent of the spin configuration and at a moderately high luminosity . The authors suggest that the material in the inner disk, instead of transforming into an advection dominated accretion flow (ADAF), could be ejected in form of a jet. This could explain the observed periodic dips in X-ray luminosity that are accompanied by large radio bursts (King et al., 2011); From a theoretical point of view, models with an inner truncation radius are frequently considered, but in the context of a two-zone geometry where an inner hot optically thin and geometrically thick ADAF is enclosed by a truncated outer standard disk (Esin et al., 1997; Yuan et al., 2005; Liu & Taam, 2009, 2013; Yuan & Narayan, 2014). In this case, the ADAF fully describes the low luminosity/hard emission state where the accretion rate drops below a few percent of the Eddington limit (Zdziarski et al., 2004). However, this trend goes in the opposite direction compared to the trend required to solve the spin paradox in our paper: the transition radius in the case of standard disk/ADAF hybrid solution decreases with the increase of the Eddington ratio while we need an opposite trend. Therefore, the mechanism that acts in high Eddington ratio sources must be a different one. Observationally, the presence of the relativistically broadened reflection provides the support for the disk close to the black hole. We see it in some of the RXTE spectra of GRS 1915+105, and soft state reflection has been detected in some other galactic sources in the soft state, e.g. Cyg X-1 (Tomsick et al., 2014); XTE J1908+094 (Tao et al., 2015); LMC X-1 (Koyama et al., 2015); 4U 1543-47 (Morningstar & Miller, 2014), as well as in AGNs, e.g. Mkn 335 (Wilkins et al., 2015). In this case only a small fraction of material should be outflowing, but this is not fully self-consistently treated in the Model (4) as the temperature of the disk with magnetically driven wind should be much lower than the one for the absence of the term in Model (4), which affects the color-correction term. Therefore, we do not use this model when fitting the observational data in Sec. 5.
6.1. Comparison between disks with and without wind
There is only one significant difference between disks with and without wind, namely the predicted mass accretion rate. For a given observation, a disk with initial mass accretion rate modified by wind can exhibit accretion rates significantly larger than a disk without wind. This is a direct consequence of being measured at the outer disk edge and not being constant with radius due to wind mass loss. Accretion disks that self-consistently include winds are thus much less efficient at converting mass into radiation than a standard thin disk.
The spectral difference in the shape between the models with and without wind is very small, when the model parameters are appropriately adjusted. In a given observation, the spin, accretion rate, and wind prescription are strongly degenerate. They are additionally masked by the Comptonization at higher energies. A sequence of disk-dominated spectra is necessary, since then the spin should remain the same, independently from the source luminosity. Models without wind and with hardening factor taken from bhspec do not provide the correct solution for such a sequence for GRS 1915+105 and LMC X-3, while the wind model presented in this paper is satisfactory.
The comparison of the disk wind model with the data is a first step towards the assessment of the flow geometry. However, next steps are very difficult, since the material removed from the disk flows out as a wind but also (mostly) flows in as a corona flow. Comptonization computations should be done in 3-D, for the adopted flow pattern. The basic difficulty is that we then need to know the temperature distribution inside the hot flow which depends on the plausible local dissipation, and the global magnetic field will affect the global flow.
The success of the wind model in reproducing the right trends over a range of mass accretion rates does not mean that the wind is the only solution to the spin problem. Straub et al. (2011) have shown that if a constant hardening factor is used instead of a variable hardening factor predicted by bhspec and slimbh, the spin problem is significantly alleviated. We find that fits with slimbh are just as good as those with diskw and, consequently, that the RXTE data are not sensitive to the shape of the disk model spectrum. Suzaku (Mitsuda et al., 2007) is a satellite that has a moderate spectral resolution CCD for the 0.7-10 keV band and a hard X-ray detector for the 12-30 keV band. Kubota et al. (2010) used Suzaku observations of LMC X-3 to constrain the spectral shape of the intrinsic disk emission in various models. They found that bhspec which self-consistently includes radiative transfer through the vertical structure of the disk gives an excess at the absorption edge energy below 1 keV compared to the data. This may suggest the need to include more physical processes in the disk atmosphere, i.e. atomic physics and better assumptions for the disk density and emissivity profile with respect to the disk height.
One way to prove the physical existence of the wind close to the black hole horizon would be to see it directly in the spectral shape of the pure disk component, due to general relativity effects. However, this would require a data set with very low systematic errors, below 0.2%. In addition, the precision in the description of the Comptonization process would also matter in this case, including the 3-D distribution of the static Comptonizing medium and taking into account the material motion would pose yet another challenge (Beloborodov, 1999).
6.2. Observational evidence of winds
Outflows are expected and actually observed in accreting systems. Their occurrence in the form of collimated jets is characteristic of low Eddington rate sources. While jets seem to be suppressed in the high/soft state of black hole XRBs (Qiao & Liu, 2015), uncollimated outflows/winds are still found in high luminosity sources (Miller et al., 2006b; King et al., 2014). However, the presence of the specific outflow we request here to solve the problem of the apparent spin decrease is difficult to test. The outflow temperature must be at least K or higher and the outflow may, or may not, be located along the line of sight, depending on the collimation which is unspecified. The material is thus highly, or fully ionized leaving only a possibility of emission/absorption from highly ionized iron.
There are ample observational evidences for outflows from accretion disks. Outflowing winds have recently been detected in several XRBs via blueshifted X-ray absorption lines which are observed in high resolution spectra (see Miller et al., 2015, and references therein). With modeling the observed narrow absorption features recorded in GRS 1915+105 by Chandra HETGS (High Energy Transmission Grating spectra), the total mass outflow rate is estimated to be comparable to the mass accretion rate in the inner part of the disk (Ueda et al., 2009). A re-analysis of the Chandra HETGS spectroscopy data of GRS 1915+105 by Miller et al. (2015) with improved multi-zone photoionization models and with different ionization parameters and velocities reveals that when the accreted gas rate , where the Eddington accretion rate in that work is defined as: .
These detected winds, however, clearly do not correspond to the winds requested in our paper to solve the problem of the apparent spin evolution. Our winds should happen very close to the black hole, at a few gravitational radii, and with the terminal speed of the order of the local Keplerian velocity at these radii, i.e. at a large fraction of the light speed. All winds discussed above show narrow absorption features with shifts corresponding to a few hundreds to a few thousands km s, e.g., 4U 1630472: the velocity and the launching radius ; GRO J165540: and ; H 1743322: and ; GRS 1915+105: and (see Table 4-6 and 12 of Miller et al., 2015). Thus, observations usually give constraints on the partially ionized winds from the outer disk while here we deal with the highly ionized winds from the inner disk.
Magnetically driven outflows have been postulated and studied from a theoretical point of view (e.g. Blandford & Payne, 1982; Konigl & Kartje, 1994; Contopoulos & Lovelace, 1994), and efficient outflow close to the inner disk radius due to the bending of the poloidal field is expected (Campbell, 2010). Observationally, magnetic assistance has been advocated in some outflows, like ultra-fast outflows in active galaxies. Ultra-fast outflows have been detected in a number of AGN through the analysis of absorptions features. For example, in PG 1211+143 Pounds et al. (2003) identified highly ionized outflows which had a velocity of and mass flux comparable to the mass accretion rate. A more recent spectral analysis of XMM-Newton/EPIC data of PG 1211+143 revealed an outflow velocity of , with most of the wind launched at 200 and an inner wind truncation radius of (Fukumura et al., 2015), which favors magnetic driving. The location of this wind, however, is still too far from the black hole for our purpose.
In principle, the motion of the plasma is encoded in the way how a fully ionized plasma up-scatters the photons. Since the terminal outflow velocities are large while the initial denser part of the wind has a lower velocity, the expected effects of the bulk Comptonization are difficult to model, without assuming a specific velocity profile. The issues of bulk Comptonization were occasionally discussed both in the case of an inflow (e.g. Chakrabarty & Titarchuk 1995; Liu et al. 2015) and an outflow (e.g. Beloborodov, 1999). Computations of the Comptonization of the disk photons in the postulated winds with density/velocity/temperature gradients is beyond the scope of the present simple work.
The determination of the black hole spin in XRBs using advanced fully relativistic disk models including full radiative transfer to account for Comptonization in the disk atmosphere leads to a paradox: the spin decreases with increasing source luminosity (McClintock et al., 2006; Steiner et al., 2010; Straub et al., 2011). Here, we analyze whether this effect might be explained by the presence of the winds from the innermost part of the accretion disk. We find that relatively smooth radial winds without explicit dependence on mass accretion rate (Models (1)-(3) in Sec. 3) do not provide the solution for this effect. However, if we assume in Model that parameter depends on (the mass accretion rate at the outer disk radius) the apparent spin decrease can be reproduced (green stars in Fig. 4). We can recover the observed trend also in Model (4) where the truncated inner disk radius depends linearly on . The physical interpretation of such a model is that either most of the material, or at least most of the energy and angular momentum, has to be removed from the disk in a form of a wind or an uncollimated jet, leaving behind a cold non-radiating flow close to the black hole horizon.
In all these scenarios the shape of the thermal disk component is modified, leaving a trace of the process in the spectral shape. We show that the spectra for the two cases of outflow discussed in Sec. 4, namely the hot optically thin and the cold optically thick wind, are in comparison to a spectrum without Comptonization either visibly harder or softer, respectively. However, if the Comptonizing medium is chosen independently from the predicted outflow, the difference between the model with, and without outflow is very small, at the energy where the disk flux before Comptonization drops by a factor of two below the peak flux (see Fig. 7).
We test in more detail the disk-wind scenario against the data for two bright XRBs, LMC X-3 and GRS 1915+105, using the Model and assuming a Comptonizing medium surrounding the disk. For each source we fit several hundreds of high/soft spectra assuming an intrinsically constant black hole spin and the presence of the wind that correlates with mass accretion rate. We obtain good fit statistics as demonstrated for four representative spectra for each source (Figures 9, 10 and Tables 2, 3). The apparent decrease of the black hole spin over a range of mass accretion rates in LMC X-3 and GRS 1915+105 could originate from the presence of wind or be a consequence of an incomplete description of the disk atmosphere. The latter issue will be addressed in a forthcoming paper. Here, we show that the wind model with a fixed spin can fit the data.
However, the implied properties of the Comptonizing medium do not correlate with the required wind, as the wind intensity rises sharply with the luminosity while the Comptonizing medium properties remain almost unchanged. This implies the presence of a separate Comptonizing medium, likely in the form of an X-ray source located on the symmetry axis, while the material removed from the disk, albeit in large amount, can escape detection if it flows predominantly inward as a corona flow. Thus we have no direct observational evidence of the requested wind from the inner few gravitational radii.
Our warmest and most heartfelt thanks go to R. Remillard and J. F. Steiner who shared their data for this study. We thank F. Yuan for his helpful comments. Furthermore, we are greatly indebted to the anonymous referee whose to the point comments and suggestions helped to significantly improve the manuscript. The project has received funding from the European Union Seventh Framework Program (FP7/2007-2013) under grant agreement No. 312789. This work was also supported by the Ministry of Science and Higher Education grant W30/7.PR/2013 and the Polish National Science Center grants No. 2011/03/B/ST9/03281 and No. 2013/08/A/ST9/00795.
- Abramowicz et al. (1988) Abramowicz, M. A., Czerny, B., Lasota, J. P., et al. 1988, ApJ, 332, 646
- Arnaud (1996) Arnaud, K. A. 1996, Astronomical Data Analysis Software and Systems V, 101, 17
- Balsara & Krolik (1993) Balsara, D. S., & Krolik, J. H. 1993, ApJ, 402, 109
- Beloborodov (1999) Beloborodov, A. M. 1999, ApJ, 510, L123
- Blaes et al. (2011) Blaes, O., Krolik, J. H., Hirose, S., & Shabaltas, N. 2011, ApJ, 733, 110
- Blandford & Znajek (1977) Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433
- Blandford & Payne (1982) Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883
- Brenneman & Reynolds (2006) Brenneman, L. W., & Reynolds, C. S. 2006, ApJ, 652, 1028
- Campbell (2010) Campbell, C. G. 2010, MNRAS, 401, 177
- Cao (2009) Cao, X. 2009, MNRAS, 394, 207
- Cao (2010) Cao, X. 2010, ApJ, 725, 388
- Cao (2014) Cao, X. 2014, ApJ, 783, 51
- Cao (2016) Cao, X. 2016, ApJ, 817, 71
- Capellupo et al. (2015) Capellupo, D. M., Netzer, H., Lira, P., et al. 2015, MNRAS, 446, 3427
- Chakrabarti & Titarchuk (1995) Chakrabarti, S., & Titarchuk, L. G. 1995, ApJ, 455, 623
- Connors et al. (1980) Connors, P. A., Stark, R. F., & Piran, T. 1980, ApJ, 235, 224
- Contopoulos & Lovelace (1994) Contopoulos, J., & Lovelace, R. V. E. 1994, ApJ, 429, 139
- Coppi (1999) Coppi, P. S. 1999, High Energy Processes in Accreting Black Holes, 161, 375
- Cowperthwaite & Reynolds (2012) Cowperthwaite, P. S., & Reynolds, C. S. 2012, ApJ, 752, L21
- Czerny & Zbyszewska (1991) Czerny, B., & Zbyszewska, M. 1991, MNRAS, 249, 634
- Czerny et al. (2011) Czerny, B., Hryniewicz, K., Nikołajuk, M., et al. 2011, MNRAS, 415, 2942
- Davis et al. (2006) Davis, S. W., Done, C., & Blaes, O. M. 2006, ApJ, 647, 525
- Davis & Hubeny (2006) Davis, S. W., & Hubeny, I. 2006, ApJS, 164, 530
- Done et al. (2007) Done, C., Gierliński, M., & Kubota, A. 2007, A&A Rev., 15, 1
- Done et al. (2012) Done, C., Davis, S. W., Jin, C., et al. 2012, MNRAS, 420, 1848
- Dovčiak et al. (2008) Dovčiak, M., Muleri, F., Goosmann, et al. 2008, MNRAS, 391, 32
- Esin et al. (1997) Esin, A. A., McClintock, J. E., & Narayan, R. 1997, ApJ, 489, 865
- Fabian et al. (2015) Fabian, A. C., Lohfink, A., Kara, E., et al. 2015, MNRAS, 451, 4375
- Fender et al. (1999) Fender, R., Corbel, S., Tzioumis, T., et al. 1999, ApJ, 519, L165
- Fender et al. (2004) Fender, R. P., Belloni, T. M., & Gallo, E. 2004, MNRAS, 355, 1105
- Filiz Ak et al. (2013) Filiz Ak, N., Brandt, W. N., Hall, P. B., et al. 2013, ApJ, 777, 168
- Fragos & McClintock (2015) Fragos, T., McClintock, J. E. 2015, ApJ, 800, 17
- Fukumura et al. (2015) Fukumura, K., Tombesi, F., Kazanas, D., et al. 2015, ApJ, 805, 17
- Gou et al. (2014) Gou, L., McClintock, J. E., Remillard, R. A., et al. 2014, ApJ, 790, 29
- Gu et al. (2015) Gu, M., Chen, Y., Komossa, S., et al. 2015, ApJS, 221, 3
- Higginbottom et al. (2014) Higginbottom, N., Proga, D., Knigge, C., et al. 2014, ApJ, 789, 19
- Hubeny & Lanz (1995) Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875
- Janiuk et al. (2002) Janiuk, A., Czerny, B., & Siemiginowska, A. 2002, ApJ, 576, 908
- Janiuk & Czerny (2005) Janiuk, A., & Czerny, B. 2005, MNRAS, 356, 205
- Kallman et al. (2009) Kallman, T. R., Bautista, M. A., Goriely, S., et al. 2009, ApJ, 701, 865
- Kaspi et al. (2002) Kaspi, S., Brandt, W. N., George, I. M., et al. 2002, ApJ, 574, 643
- King & Pounds (2003) King, A. R., & Pounds, K. A. 2003, MNRAS, 345, 657
- King et al. (2008) King, A. R., Pringle, J. E., & Hofmann, J. A. 2008, MNRAS, 385, 1621
- King et al. (2011) King, A. L., Miller, J. M., Cackett, E. M., et al. 2011, ApJ, 729, 19
- King et al. (2012) King, A. L., Miller, J. M., Raymond, J., et al. 2012, ApJ, 746, L20
- King et al. (2014) King, A. L., Walton, D. J., Miller, J. M., et al. 2014, ApJ, 784, L2
- King & Muldrew (2016) King, A., & Muldrew, S. I. 2016, MNRAS, 455, 1211
- Konigl & Kartje (1994) Konigl, A., & Kartje, J. F. 1994, ApJ, 434, 446
- Koyama et al. (2015) Koyama, S., Yamada, S., Kubota, A., et al. 2015, PASJ, 67, 46
- Kubota & Done (2004) Kubota, A., & Done, C. 2004, MNRAS, 353, 980
- Kubota et al. (2010) Kubota, A., Done, C., Davis, S. W., et al. 2010, ApJ, 714, 860
- Laor (1991) Laor, A. 1991, ApJ, 376, 90
- Laor & Davis (2014) Laor, A., & Davis, S. W. 2014, MNRAS, 438, 3024
- Lee et al. (2002) Lee, J. C., Reynolds, C. S., Remillard, R., et al. 2002, ApJ, 567, 1102
- Li et al. (2005) Li, L.-X., Zimmerman, E. R., Narayan, R., et al. 2005, ApJS, 157, 335
- Li et al. (2009) Li, L.-X., Narayan, R., & McClintock, J. E. 2009, ApJ, 691, 847
- Li & Cao (2012) Li, S.-L., & Cao, X. 2012, ApJ, 753, 24
- Liu & Taam (2009) Liu, B. F., & Taam, R. E. 2009, ApJ, 707, 233
- Liu & Taam (2013) Liu, B. F., & Taam, R. E. 2013, ApJS, 207, 17
- Liu et al. (2015) Liu, B. F., Taam, R. E., Qiao, E., & Yuan, W. 2015, ApJ, 806, 223
- Liu et al. (2015) Liu, Z., Yuan, W., Lu, Y., & Zhou, X. 2015, MNRAS, 447, 517
- Magdziarz et al. (1998) Magdziarz, P., Blaes, O. M., Zdziarski, A. A., et al. 1998, MNRAS, 301, 179
- McClintock et al. (2006) McClintock, J. E., Shafee, R., Narayan, R., et al. 2006, ApJ, 652, 518
- McClintock & Remillard (2009) McClintock, J. E., & Remillard, R. A. 2009, astro2010: AADS, 2010, 197
- Merloni & Nayakshin (2006) Merloni, A., & Nayakshin, S. 2006, MNRAS, 372, 728
- Middleton et al. (2006) Middleton, M., Done, C., Gierliński, M., & Davis, S. W. 2006, MNRAS, 373, 1004
- Miller et al. (2002) Miller, J. M., Fabian, A. C., Wijnands, R., et al. 2002, ApJ, 570, L69
- Miller et al. (2006a) Miller, J. M., Raymond, J., Fabian, A., et al. 2006a, Nature, 441, 953
- Miller et al. (2006b) Miller, J. M., Raymond, J., Homan, J., et al. 2006b, ApJ, 646, 394
- Miller et al. (2008) Miller, J. M., Reynolds, C. S., Fabian, A. C., et al. 2008, ApJ, 679, L113
- Miller et al. (2013) Miller, J. M., Parker, M. L., Fuerst, F., et al. 2013, ApJ, 775, L45
- Miller et al. (2015) Miller, J. M., Fabian, A. C., Kaastra, J., et al. 2015, ApJ, 814, 87
- Mitsuda et al. (1984) Mitsuda, K., Inoue, H., Koyama, K., & Makishima, K. 1984, PASJ, 36, 741
- Mitsuda et al. (2007) Mitsuda, K., Bautz, M., Inoue, H., et al. 2007, PASJ, 59, 1
- Morningstar & Miller (2014) Morningstar, W. R., & Miller, J. M. 2014, ApJ, 793, L33
- Motta et al. (2014) Motta, S. E., Belloni, T. M., Stella, L., et al. 2014, MNRAS, 437, 2554
- Narayan & McClintock (2012) Narayan, R., & McClintock, J. E. 2012, MNRAS, 419, L69
- Nayakshin et al. (2000) Nayakshin, S., Rappaport, S., & Melia, F. 2000, ApJ, 535, 798
- Neilsen & Lee (2009) Neilsen, J., & Lee, J. C. 2009, Nature, 458, 481
- Novikov & Thorne (1973) Novikov, I. D., & Thorne, K. S. 1973, Black Holes (Les Astres Occlus), 343
- Orosz et al. (2014) Orosz, J. A., Steiner, J. F., McClintock, J. E., et al. 2014, ApJ, 794, 154
- Page et al. (2003) Page, M. J., Soria, R., Wu, K., et al. 2003, MNRAS, 345, 639
- Parker (1958) Parker, E. N. 1958, ApJ, 128, 664
- Pasham et al. (2014) Pasham, D. R., Strohmayer, T. E., & Mushotzky, R. F. 2014, Nature, 513, 74
- Peris et al. (2015) Peris, C. S., Remillard, R. A., Steiner, J. F., et al. 2015, arXiv:1509.08941
- Petrucci et al. (2013) Petrucci, P.-O., Paltani, S., Malzac, J., et al. 2013, A&A, 549, A73
- Pounds et al. (2003) Pounds, K. A., Reeves, J. N., King, A. R., et al. 2003, MNRAS, 345, 705
- Proga et al. (2000) Proga, D., Stone, J. M., & Kallman, T. R. 2000, ApJ, 543, 686
- Qiao & Liu (2015) Qiao, E., & Liu, B. F. 2015, MNRAS, 448, 1099
- Reid et al. (2014) Reid, M. J., McClintock, J. E., Steiner, J. F., et al. 2014, ApJ, 796, 2
- Remillard & McClintock (2006) Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 49
- Reynolds (2012) Reynolds, C. S. 2012, ApJ, 759, L15
- Różańska et al. (2015) Różańska, A., Malzac, J., Belmont, R., et al. 2015, A&A, 580, A77
- Sa̧dowski (2009) Sa̧dowski, A. 2009, ApJS, 183, 171
- Sa̧dowski et al. (2011) Sa̧dowski, A., Abramowicz, M., Bursa, M., et al. 2011, A&A, 527, AA17
- Shakura & Sunyaev (1973) Shakura N. I., & Sunyaev R. A. 1973, A&A, 24, 337
- Shlosman et al. (1985) Shlosman, I., Vitello, P. A., & Shaviv, G. 1985, ApJ, 294, 96
- Slone & Netzer (2012) Slone, O., & Netzer, H. 2012, MNRAS, 426, 656
- Steiner et al. (2010) Steiner, J. F., McClintock, J. E., Remillard, R. A., et al. 2010, ApJ, 718, L117
- Steiner et al. (2014a) Steiner, J. F., McClintock, J. E., Orosz, J. A., et al. 2014a, ApJ, 783, 101
- Steiner et al. (2014b) Steiner, J. F., McClintock, J. E., Orosz, J. A., et al. 2014b, ApJ, 793, LL29
- Stella et al. (1999) Stella, L., Vietri, M., & Morsink, S. M. 1999, ApJ, 524, L63
- Straub et al. (2011) Straub, O., Bursa, M., Sa̧dowski, A., et al. 2011, A&A, 533, AA67
- Straub et al. (2013) Straub, O., Middleton, M., & Done, C. 2013, A&A, 553, A61
- Tao et al. (2015) Tao, L., Tomsick, J. A., Walton, D. J., et al. 2015, ApJ, 811, 51
- Tao & Blaes (2013) Tao, T., & Blaes, O. 2013, ApJ, 770, 55
- Tomsick et al. (2014) Tomsick, J. A., Nowak, M. A., Parker, M., et al. 2014, ApJ, 780, 78
- Ueda et al. (2009) Ueda, Y., Yamaoka, K., & Remillard, R. 2009, ApJ, 695, 888
- Volonteri et al. (2008) Volonteri, M., Lodato, G., & Natarajan, P. 2008, MNRAS, 383, 1079
- Walter & Fink (1993) Walter, R., & Fink, H. H. 1993, A&A, 274, 105
- Wang et al. (2009) Wang, J.-M., Hu, C., Li, Y.-R., et al. 2009, ApJ, 697, L141
- Wilkins & Gallo (2015) Wilkins, D. R., & Gallo, L. C. 2015, MNRAS, 449, 129
- Wilkins et al. (2015) Wilkins, D. R., Gallo, L. C., Grupe, D., et al. 2015, MNRAS, 454, 4440
- Wilms et al. (2000) Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914
- Wilms et al. (2001) Wilms, J., Reynolds, C. S., Begelman, M. C., et al. 2001, MNRAS, 328, L27
- Woods et al. (1996) Woods, D. T., Klein, R. I., Castor, J. I., McKee, C. F., & Bell, J. B. 1996, ApJ, 461, 767
- Wu et al. (2011) Wu, Q., Cao, X., & Wang, D.-X. 2011, ApJ, 735, 50
- You et al. (2012) You, B., Cao, X., & Yuan, Y.-F. 2012, ApJ, 761, 109
- Yuan et al. (2005) Yuan, F., Cui, W., & Narayan, R. 2005, ApJ, 620, 905
- Yuan et al. (2012) Yuan, F., Bu, D., & Wu, M. 2012, ApJ, 761, 130
- Yuan et al. (2015) Yuan, F., Gan, Z., Narayan, R., et al. 2015, ApJ, 804, 101
- Yuan & Narayan (2014) Yuan, F., & Narayan, R. 2014, ARA&A, 52, 529
- Yuan et al. (2010) Yuan, W., Liu, B. F., Zhou, H., & Wang, T. G. 2010, ApJ, 723, 508
- Zdziarski et al. (1996) Zdziarski, A. A, Johnson, W. N., & Magdziarz, P. 1996, MNRAS, 283, 193
- Zdziarski et al. (2004) Zdziarski, A. A., Gierliński, M., Mikołajewska, J., et al. 2004, MNRAS, 351, 791
- Zhang et al. (1997) Zhang, S. N., Cui, W., & Chen, W. 1997, ApJ, 482, L155
- Zhang et al. (2000) Zhang, S. N., Cui, W., Chen, W., et al. 2000, Science, 287, 1239
- Zycki et al. (1995) Zycki, P. T., Collin-Souffrin, S., & Czerny, B. 1995, MNRAS, 277, 70
- Życki et al. (1999) Życki P. T., Done, C., & Smith, D. A., 1999, MNRAS, 309, 561