Volume 540, April 2012
|Number of page(s)||4|
|Published online||16 March 2012|
The physical origin of the X-ray power spectral density break timescale in accreting black holes
1 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK
2 ISDC Data Centre for Astrophysics, ch. d’Ecogia 16, 1290 Versoix, Switzerland
3 Geneva Observatory, Geneva University, ch. des Maillettes 51, 1290 Sauverny, Switzerland
Received: 26 January 2012
Accepted: 1 March 2012
X-ray variability of active galactic nuclei (AGN) and black hole binaries can be analysed by means of the power spectral density (PSD). The break observed in the power spectrum defines a characteristic variability timescale of the accreting system. The empirical variability scaling that relates characteristic timescale, black hole mass, and accretion rate () extends from supermassive black holes in AGN down to stellar-mass black holes in binary systems. We suggest that the PSD break timescale is associated with the cooling timescale of electrons in the Comptonisation process at the origin of the observed hard X-ray emission. We find that the Compton cooling timescale directly leads to the observational scaling and naturally reproduces the functional dependence on black hole mass and accretion rate (). This result simply arises from general properties of the emission mechanism and is independent of the details of any specific accretion model.
Key words: black hole physics / radiation mechanisms: general / galaxies: active / X-rays: binaries / X-rays: galaxies
© ESO, 2012
X-ray variability is a characteristic property of all accreting black hole systems, including supermassive black holes in active galactic nuclei (AGN) and stellar-mass black holes in binary systems (black hole binaries, BHB). The observation of rapid and aperiodic variations indicates that X-ray emission is produced in the innermost regions of the accretion flow, close to the central black hole. The study of X-ray variability is thus thought to provide information on the physical emission mechanisms and possibly on the nature of the underlying accretion flows.
Quantitative analysis of X-ray variability is often performed by means of the power spectral density (PSD). The observed variability power spectrum can be roughly modelled by power laws of the form Pν ∝ να, where ν = 1/T is the temporal frequency and α the power law slope. A typical AGN power spectrum is characterised by a slope of α ~ −2 at high frequencies (corresponding to short timescales), flattening to a slope of α ~ −1 at lower frequencies (or equivalently longer timescales). A similar shape is observed in the power spectra of galactic black holes in the high/soft state. In addition, a second break to a slope of α ~ 0, occurring at still lower frequencies, is seen in the low/hard state of black hole binaries. In the following we consider the high-frequency PSD break, νB (from α ~ −1 to α ~ −2), in both AGN and black hole binaries. The corresponding break timescale, TB = 1/νB, is considered to be a characteristic X-ray variability timescale of the accreting system. Typical break timescales measured in AGN are of ~day-timescales, while the corresponding break timescales seen in black hole binaries are of the order of fractions of seconds (McHardy 2010; Gilfanov 2010). The difference in break timescales observed between supermassive and stellar-mass black holes is broadly consistent with a simple scaling with black hole mass, and linear mass-timescale relations have been proposed (Markowitz et al. 2003; Papadakis 2004).
A more quantitative X-ray variability scaling has been obtained by McHardy et al. (2006), who considered the explicit dependence of the characteristic timescale on both black hole mass and accretion rate. A major result is that this observational scaling relation extends from supermassive black holes in AGN down to stellar-mass black holes in binary systems, covering many orders of magnitude in black hole mass. The fact that both AGN and black hole binaries seem to follow the same empirical scaling has been claimed to be a strong support for the scale invariance of accreting black hole systems, and led to the suggestion that AGN are simply scaled-up galactic black holes (McHardy et al. 2006). This in turn led to the proposal of the so-called “variability fundamental plane”, which relates characteristic timescale, black hole mass, and accretion rate.
However, the physical origin of the PSD break is not definitively established yet. Various interpretations of the PSD break timescale have been proposed in the literature, trying to link the observed break timescales with different accretion model timescales. At first, the measured break timescales have been compared with characteristic physical timescales of the accretion disc, including orbital, thermal, and viscous timescales at different radii (Markowitz et al. 2003; Papadakis 2004). More recently, the PSD break has been discussed in the framework of the inner propagating accretion flow fluctuations model, in which fluctuations originating at different outer radii propagate inwards, leading to modulations of the accretion rate at inner radii and hence to variations in the observed X-ray flux (Lyubarskii 1997; Churazov et al. 2001). The optically thick accretion disc is assumed to be truncated at some distance from the centre and the inner region is replaced by a hot, optically thin flow component. In this picture, the break timescale is associated with a characteristic timescale, usually the viscous time, at the truncation radius of the standard accretion disc. In an alternative approach, we have discussed the characteristic X-ray variability timescale within the framework of a clumpy accretion model in AGN (Ishibashi & Courvoisier 2009b). In this scenario, accretion onto the central black hole occurs through a sequence or cascade of shocks, which are at the origin of the observed radiation. Accreting elements first interact in the form of optically thick shocks, giving rise to optical/UV emission, later followed by optically thin shocks giving rise to X-rays through inverse Compton emission. Optical/UV photons emitted in the optically thick shocks form the seed photons for Compton-cooling of the electrons heated in the optically thin shocks. The characteristic heating/cooling timescale of the electrons in the optically thin shocks has been associated with the PSD break timescale.
Although the exact geometry of the inner coronal flow is still unknown, hard X-ray emission is widely attributed to Compton processes. In the inverse Compton process, soft seed photons (possibly emitted as thermal emission from the accretion disc) are upscattered to higher energies by energetic electrons in a hot, optically thin corona. The observed hard X-ray power law component of accreting black holes is indeed well represented by Comptonisation models (Titarchuk 1994, and references therein). As hard X-rays are produced by the Comptonisation mechanism, it seems natural to consider physical timescales associated with the emission process itself as possible candidates for the identification with the observed X-ray PSD break timescale. In the following we consider the characteristic timescale associated with the Compton cooling of the hot electrons.
The power emitted by a single electron in the inverse Compton scattering is given by (Rybicki & Lightman 1979) (1)where σT is the Thomson cross section and uγ the photon energy density.
The Compton cooling timescale, which sets the characteristic timescale for the cooling of electrons by inverse Compton emission, is defined as (2)where Ee is the average electron energy. In the non-relativistic limit (γ ~ 1, β2 ∝ Ee), the dependence on the electron energy cancels, and Eq. (2) leads to (3)The radiation density at distance R is given by (4)where Ls is the seed photon luminosity for inverse Compton scattering. Assuming that the seed photons are generated in the radiatively efficient accretion process, the seed photon luminosity is directly proportional to the accretion rate, Ls = ηṀc2, where η ~ 0.1 is the canonical radiative efficiency. The radial distance R can be expressed in units of the Schwarzschild radius as R = ζRS, where .
The Compton cooling timescale (Eq. (3)) is now written as (5)We see that the Compton cooling time is mainly determined by two of the black hole fundamental parameters, i.e. mass and accretion rate. As the Compton cooling time is directly linked to the physical process of X-ray emission, it can be considered to be a characteristic X-ray timescale of the accreting system.
For a typical AGN of mass MBH ~ 107 M⊙ and accreting in the radiatively efficient regime (with typical Eddington ratios of L/LE ~ 0.1), the numerical value of the characteristic timescale (Eq. (5)) is of the order of (6)where MBH = M7 × 107 M⊙ and Ṁ = Ṁ-2 × 0.01 M⊙/yr. We see that the characteristic Compton cooling timescale in AGN is of ~day-timescales. We can also estimate the Compton cooling timescale in black hole binaries by introducing parameters appropriate for galactic black holes in Eq. (5). Considering an object of about MBH ~ 10 M⊙ accreting at a few percent of the Eddington limit, and taking into account the fact that the seed photon luminosity is only a fraction of the bolometric luminosity, the characteristic timescale is of the order of tBHB ~ 0.07 s.
We observe that the Compton cooling timescale is shorter for small black hole mass and/or high accretion rate, scaling as (7)The Compton cooling of the hot electrons leads to X-ray emission. The average electron energy can be estimated by assuming an equilibrium situation between Compton cooling and Coulomb heating rates: LC = LH (Ishibashi & Courvoisier 2009a). The resulting electron energy scales as , where ne is the electron number density. The average X-ray luminosity is given by the Compton luminosity of a single electron (Eq. (1)) multiplied by the average number of electrons: (8)and we see that the dependence on the electron number density is stronger than that on the photon energy density. Thus X-ray luminosity variations are mainly driven by electron density modulations associated with the inhomogeneities of the upscattering medium. In the clumpy accretion scenario, the origin of the density modulations was associated with the inhomogeneous structure of the accretion flow resulting from the succession of shocks (Ishibashi & Courvoisier 2009b). Here we note that the physical origin of the modulations can be diverse and does not need to be specified. Density fluctuations on timescales shorter than the electron cooling time (tf < tC), which at equilibrium is equal to the electron heating time, cannot be propagated and are smoothed out. The suppression of fluctuations on timescales shorter than a critical value, given by the Compton cooling time, induces a break in the observed variability power spectrum, which may be associated with the PSD break.
Observational measurements of the high-frequency PSD break (from α ~ −1 to α ~ −2) have been performed for a number of AGN and galactic black holes. The observed break timescales in Seyfert galaxies are usually in the range ~0.01 to a few days (Markowitz et al. 2003; McHardy et al. 2004; Uttley & McHardy 2005), while the corresponding break timescales measured in black hole binaries are of the order of seconds or less (Gilfanov 2010, and references therein). The Compton cooling timescales estimated in the previous section (~day-timescales in AGN and fractions of seconds in BHB) thus correspond to the typically measured values of the PSD break in both AGN and galactic black holes. The characteristic timescale therefore seems to scale roughly linearly with black hole mass, from stellar-mass black holes up to supermassive black holes. However, an additional dependence on the accretion rate is required by more detailed observations of AGN: for a given black hole mass, higher accretion rate objects have shorter break timescales than lower accretion rate counterparts (Uttley & McHardy 2005). The empirical relationship between variability timescale, black hole mass, and accretion rate has been obtained by McHardy et al. (2006) based on a main sample of AGN including two galactic black holes in radio-quiet states. The resulting empirical scaling can be written as (9)where ṁE = Lbol/LE is the ratio of the bolometric luminosity to Eddington luminosity. This observational scaling is shown to extend from AGN down to stellar-mass black holes in binary systems. As the Eddington luminosity is given by the black hole mass (LE ∝ MBH) and assuming that the bolometric luminosity is proportional to the accretion rate (Lbol ∝ Ṁ), the empirical scaling (Eq. (9)) can be re-expressed directly in terms of the physical parameters as (10)Comparison of the Compton cooling timescale, (derived in Eq. (5)), with the observational scaling (Eq. (10)) shows a very close match. We observe that the characteristic timescale associated with the Compton cooling process directly reproduces the observed dependence on black hole mass and accretion rate. We note that a similar result has been previously derived but within the specific framework of the clumpy accretion model in AGN (Ishibashi & Courvoisier 2009b).
X-ray variability properties are also closely coupled with X-ray spectral properties. Observations of “spectral-timing” correlations have been reported in both AGN and black hole binaries, whereby an increase in the characteristic frequency is correlated with a steepening of the X-ray spectrum. The characteristic noise frequency is observed to increase with the steepening of the energy spectrum in Cygnus X-1 (Gilfanov et al. 1999), and similar patterns are also observed in other X-ray binary systems, e.g. GX 339-4 (Revnivtsev et al. 2001). In the case of AGN, the average X-ray spectral slope is found to be positively correlated with the normalized characteristic frequency defined as νnorm = νB × MBH, where νB is the break frequency observed in the X-ray PSD (Papadakis et al. 2009).
It is known that Comptonisation processes lead to power law spectra, and the resulting X-ray spectra observed in the hard energy band are usually described by power laws, characterised by the photon index Γ. A significant correlation between photon index and Eddington ratio has been reported in numerous AGN samples (Shemmer et al. 2008; Kelly et al. 2008; Risaliti et al. 2009), suggesting that the X-ray spectral slope is correlated with some form of accretion rate. Within the Comptonisation scenario, the physical interpretation of the observed Γ − L/LE correlation is that an increase in the accretion rate leads to an enhancement in the seed photon emission, which in turn leads to efficient Compton cooling of the hot corona, which finally leads to a steepening of the hard X-ray spectrum. On the other hand, we have seen that the Compton cooling time scales inversely with the seed photon energy density (Eq. (3)), thus an increase in the seed photon emission also leads to a decrease in the characteristic timescale. Combining the two results we may account for the observed correlation between characteristic frequency and spectral steepening. Therefore in the context of Comptonisation, both the X-ray characteristic timescale and the spectral slope seem to be determined by the supply of seed photons, which in turn depends on the accretion rate. Hence, for a given black hole mass, higher accretion rate systems tend to have both steep spectra and short characteristic timescales, in agreement with the observed “spectral-timing correlation” (Papadakis et al. 2009). This suggests that both timing and spectral properties of accreting black hole systems are mainly determined by the accretion rate parameter.
The break observed in the X-ray PSD is thought to provide insight into the underlying variability mechanism and various physical interpretations of the break timescale have been proposed. As a first approach, the observed PSD break timescales have been compared with physical timescales of the accretion disc in the innermost regions of the accretion flow where X-ray emission is thought to be produced. The shortest timescale is given by the dynamical or orbital timescale, tφ. This is of the order of ~hours for typical AGN parameters and thus too short compared to the break timescales of the order of ~days observed in Seyfert galaxies. On the other hand, for a standard geometrically thin disc (with aspect ratio H/R ≪ 1) the viscous timescales (of the order of ≳ years) are found to be a few orders of magnitude longer compared to the observed values (Markowitz et al. 2003). The closest values to the measured break timescales are perhaps given by thermal timescales, and thermal instabilities have also been considered (Papadakis 2004). We note that the physical timescales of the accretion disc are all proportional to the black hole mass and thus the approximate scaling of the characteristic timescale with black hole mass is easily obtained; however, the dependence on the accretion rate is less obvious and cannot be simply accounted for.
The currently favoured interpretation of the PSD break is based on the inner propagating fluctuations model (Lyubarskii 1997; Churazov et al. 2001). The variability power is found to decrease on short timescales, with a rapid fall-off below the break, suggesting that the break timescale might be associated with some cut-off in the accreting system. Accordingly, the break timescale has been associated with the viscous timescale (tν ~ tφ/αss (H/R)2) at the inner edge of the optically thick accretion disc. The accretion disc is assumed to be truncated at some distance from the centre, the location of the truncation radius being connected with the accretion rate: at high accretion rates the optically thick disc extends down to the last stable orbit, while at low accretion rates the disc is truncated at large distances (Zdziarski & Gierliński 2004, and references therein). In this picture, the value of the characteristic timescale is related to the position of the disc truncation radius. If the disc is truncated far away from the centre, the associated viscous timescale is long; while with increasing accretion rate, the truncation radius becomes smaller and the associated characteristic timescale correspondingly shorter. Once the accretion disc has reached the last stable orbit, a further decrease of the viscous timescale can only be obtained by varying the disc thickness parameter (H/R), assuming that the H/R ratio increases with increasing accretion rate (McHardy 2010). Therefore one may qualitatively account for the observed trend with accretion rate by combining variations of the disc truncation radius with variations of the disc thickness parameter.
Although the PSD break timescale has usually been associated with some characteristic timescale of the accretion disc, it is well known that the accretion disc itself cannot account for the hard X-ray emission component of accreting black holes. Consequently, the identification of a characteristic disc timescale at a particular radius with the characteristic X-ray variability timescale is not straightforward, and other interpretations may be considered. On the other hand, it is widely agreed that hard X-ray emission is produced by Comptonisation processes. It seems therefore interesting to see whether there is any characteristic timescale associated with the emission process itself that might be compared with the observed PSD break timescale.
We have previously associated the PSD break timescale with the heating/cooling time of the electrons in the optically thin shocks within the clumpy accretion scenario (Ishibashi & Courvoisier 2009b). Here we suggest that the X-ray PSD break timescale could be associated with the cooling time of electrons in the Comptonisation process, independently of the details of any specific accretion model. The characteristic X-ray variability timescale is simply given by the Compton cooling time; the observed scaling with black hole mass and accretion rate is quantitatively reproduced, without additional assumptions. It is quite remarkable that the observational scaling simply derives from the definition of the Compton cooling timescale with a basic assumption of a radiatively efficient accretion and a quasi-spherical distribution of the seed photon population.
The Comptonisation process is thought to operate in a similar way in supermassive and stellar-mass black holes (Gliozzi et al. 2011). We have seen that numerical estimates of the Compton cooling time give correct order-of-magnitude values of the PSD break timescales in AGN and binary systems. In our picture, the characteristic timescale is set by the Compton cooling time in AGN and galactic black holes, suggesting that the underlying emission mechanism is scale-invariant. In addition, if the characteristic X-ray variability timescale is given by the
Compton cooling timescale, this may additionally support the hypothesis that the main mechanism of hard X-ray emission is indeed Comptonisation in both classes of accreting black holes.
Various physical interpretations of the PSD break have been proposed in the literature. The observed break timescale is usually associated with some characteristic timescale of the accretion disc, such as the viscous timescale at the truncation radius. Here we suggest another possibility and associate the break timescale with a physical timescale directly linked to the X-ray emission mechanism, i.e. the cooling time of the electrons in the Comptonisation process. The association seems natural since hard X-rays are generated in the Comptonisation process and it also provides a possible physical interpretation of the PSD break. We obtain that the Compton cooling time directly reproduces the correct functional dependence on black hole mass and accretion rate, defining the empirical variability scaling of McHardy et al. (2006). We show that this result is model-independent while the detailed physical realisations of the accretion flows might take different forms in different accreting systems.
W.I. acknowledges support from the Swiss National Science Foundation.
- Churazov, E., Gilfanov, M., & Revnivtsev, M. 2001, MNRAS, 321, 759 [NASA ADS] [CrossRef] [Google Scholar]
- Gilfanov, M. 2010, in Lecture Notes in Physics 794, ed. T. Belloni (Berlin: Springer Verlag), 17 [Google Scholar]
- Gilfanov, M., Churazov, E., & Revnivtsev, M. 1999, A&A, 352, 182 [NASA ADS] [Google Scholar]
- Gliozzi, M., Titarchuk, L., Satyapal, S., Price, D., & Jang, I. 2011, ApJ, 735, 16 [NASA ADS] [CrossRef] [Google Scholar]
- Ishibashi, W., & Courvoisier, T. J.-L. 2009a, A&A, 495, 113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ishibashi, W., & Courvoisier, T. J.-L. 2009b, A&A, 504, 61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kelly, B. C., Bechtold, J., Trump, J. R., Vestergaard, M., & Siemiginowska, A. 2008, ApJS, 176, 355 [NASA ADS] [CrossRef] [Google Scholar]
- Lyubarskii, Y. E. 1997, MNRAS, 292, 679 [NASA ADS] [CrossRef] [Google Scholar]
- Markowitz, A., Edelson, R., Vaughan, S., et al. 2003, ApJ, 593, 96 [NASA ADS] [CrossRef] [Google Scholar]
- McHardy, I. 2010, in Lecture Notes in Physics 794, ed. T. Belloni (Berlin: Springer Verlag), 203 [Google Scholar]
- McHardy, I. M., Papadakis, I. E., Uttley, P., Page, M. J., & Mason, K. O. 2004, MNRAS, 348, 783 [NASA ADS] [CrossRef] [Google Scholar]
- McHardy, I. M., Koerding, E., Knigge, C., Uttley, P., & Fender, R. P. 2006, Nature, 444, 730 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Papadakis, I. E. 2004, MNRAS, 348, 207 [NASA ADS] [CrossRef] [Google Scholar]
- Papadakis, I. E., Sobolewska, M., Arevalo, P., et al. 2009, A&A, 494, 905 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Revnivtsev, M., Gilfanov, M., & Churazov, E. 2001, A&A, 380, 520 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Risaliti, G., Young, M., & Elvis, M. 2009, ApJ, 700, L6 [NASA ADS] [CrossRef] [Google Scholar]
- Rybicki, G. B., & Lightman, A. P. 1979, Radiative processes in astrophysics, ed. G. B. Rybicki, & A. P. Lightman [Google Scholar]
- Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., & Kaspi, S. 2008, ApJ, 682, 81 [NASA ADS] [CrossRef] [Google Scholar]
- Titarchuk, L. 1994, ApJ, 434, 570 [NASA ADS] [CrossRef] [Google Scholar]
- Uttley, P., & McHardy, I. M. 2005, MNRAS, 363, 586 [NASA ADS] [CrossRef] [Google Scholar]
- Zdziarski, A. A., & Gierliński, M. 2004, Prog. Theor. Phys. Suppl., 155, 99 [NASA ADS] [CrossRef] [Google Scholar]
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.