Eclipsing the innermost accretion disc regions in AGN
Variable X–ray absorption has been observed in active galactic nuclei (AGN) on several time scales. Observations allow us to identify the absorber with clouds associated either with the clumpy torus (parsec scales, long timescales) or with the broad line region (BLR) (short timescales). In the latter, the cloud size has been estimated to be of the order of few gravitational radii from the observed absorption variability. Such small cloud sizes are comparable to the X–ray emitting regions so that a detailed modeling of occultation events in AGN has the potential of enabling us to infer accurately the geometry of the system. We have developed a relativistic X–ray spectral model for occultation events and we present here theoretical predictions on the different observables that can be inferred by studying X–ray eclipses in simulated XMM–Newton data. These include the size of the X–ray emitting regions as well as more fundamental parameters such as the black hole spin and the system inclination. We find that absorption varies as a function of the energy range and that its maximum takes place when the approaching part of the accretion disc is covered. Therefore we study the hard–to–soft (H / S) ratio light curves produced during an eclipse and use them to characterise the properties of the inner accretion disc in a new model–independent way.
Several spectral features shuch as Doppler boosting, gravitational redshift and light bending are accounted for in the X–ray emission reprocessed by the innermost regions of the accretion disc, which are in a strong gravity regime due to its proximity to the supermassive black hole (SMBH). In this context, both the black hole spin and the disc properties shape the K emission line profile (see e.g. Reynolds & Nowak 2003). However, there is some controversy on this interpretation: the partial covering scenario supported by Turner et al. (2007) claims that such spectral features can be explained by calling upon a variable wind composed by several ionised clouds that partially cover the X–ray continuum source. Within this paper we argue that the ultimate test to prove the detectability of general relativistic effects under extreme gravity conditions could be furnished by X–ray spectroscopy of AGN, which would definitely rule out the partial covering scenario.
On the side of special relativistic effects, Doppler boosting would have a clear and distinct effect on the reflected component: the emission from the parts of the disc approaching towards us would be enhanced while we would detect diminished flux from the receding parts. Due to the fact that the central engine is unresolved, the only way to distinguish the reflection from the various parts of the disc is by means of the obscuration of the different emitting regions by structures in our line of sight (LOS). Occultation of the X–ray source by optically thick matter was first proposed by McKernan & Yaqoob (1998). A method to investigate relativistic effects during X–ray eclipses in AGN spectra was developed by Risaliti et al. (2011a), consisting on measuring the emission from half of the disc while the other half is covered, and viceversa. Relativistic effects arise when the emission of each half are different, hinting an asymmetric emission from the disc due to Doppler boosting.
Occultation events are rather common in AGN. These objects show X–ray spectral variability on a variety of time–scales, independently of their luminosity or morphology, (see e.g. Marshall et al. 1981). Specifically, Narrow Line Seyfert 1 (NLS1) galaxies systematically show a higher amplitude variability (Leighly 1999). Also, variability can occur at time–scales as short as a few hours in all types of AGN, as shown by Elvis et al. (2004) in NGC 4388 (type 2), Bianchi et al. (2009) in NGC 7582 (type 2), Risaliti et al. (2009) and Maiolino et al. (2010) in NGC 1365 (Seyfert 1.8), Puccetti et al. (2007) in NGC 4151 (Seyfert 1.5), Risaliti et al. (2011b) in Mrk 766 (NLS1) and Sanfrutos et al. (2013) in SWIFT J2127.4+5654 (NLS1). The high variability characterised in the references above can be explained in the context of rapid column density changes, which are the result of the interception of material located in the BLR of the system to our LOS (see e.g. Miniutti et al. 2014). Hence, AGN are perfect laboratories to check general relativity (GR) effects by using eclipses to probe the innermost regions of the accretion disc.
2 The absorption model
The geometry that we assume for the system consists of a SMBH, characterised by its mass () and spin parameter (), around which an optically thick but physically thin ionised accretion disc extends from the innermost stable circular orbit (ISCO) to . About the geometry of the system formed by the disc and the X–ray source (or the so–called corona), several configurations have been provided in the literature: the hot and radiatively compact spherical corona with radius between and proposed by Fabian et al. (2015); the slab (Rozanska et al. 2015) or the patchy (Stern et al. 1995) geometries; and the jet base interpretation (Wilkins & Gallo 2015), just to mention some of them. Our goal in this work is only to perform a first exploration of the inner accretion disc obscuration, therefore we do not choose any particular geometry for the corona yet, leaving that point to a more extensive work (Sanfrutos et al. in preparation). However, the continuum X–ray source has to be located in the vicinity of the SMBH in order to properly illuminate the innermost regions of the accretion disc; i.e. in our approximation we assume that both the primary continuum and the reflection component arise together from every point on the disc.
The eclipsing cloud is assumed to be co–rotating with the disc at a velocity of 3000 km s, and therefore to be located much further away (), in the BLR of the system. Since these kind of obscuration events are more likely to be detected in unobscured sources, we disregard the torus of the Unification models within the framework of this paper.
The model we are using consists of a power law plus an X–ray reflection component arising from the accretion disc, which is represented by the XILLVER code in XSPEC (García et al. 2013). These two components are multiplied by the KYNCONV relativistic convolution model from Dovčiak et al. (2004), that produces accretion discs spectra in the strong gravity regime, allowing to obscure part of the emission with a circular cloud whose size and position can be determined. We model the accretion disc to cover at least one half of the sky as seen from the central engine, with a typical ionisation of erg cm s and solar abundances. The energy cutoff is frozen to 300 keV. The main parameters that can be tuned in the KYNCONV relativistic convolution model are: (i) the inclination of the system with respect to our LOS (: face–on, : edge–on); (ii) the black hole spin (: not rotating, : maximally rotating); (iii) the emissivity index , giving account of how the disc is illuminated (low: uniform, high: centrally concentrated); (iv) the cloud position set by the impact parameters and in units of with horizontal position being positive towards the approaching side of the accretion disc and negative towards the receding part, and with zero vertical position of the cloud (i.e. ); and (v) the cloud size in terms of its radius , also in units. The intrinsical physical properties of the obscuring cloud such as its column density and ionisation log are set by muliplying the components defined above by the ionised absorption code ZXIPCF (Reeves et al. 2008). We take into consideration the Galactic absorption too, by means of a neutral absorption component. We choose an intermediate arbitrary value of .
The procedure we perform is the following: for every position of the cloud we compute the flux in the keV energy band, then we compute the intrinsic flux in the same band by removing the cloud, and finally we derive the covering fraction (CF) at every position of the cloud all along the eclipse. Any asymmetry in these CF–light curves would be indicative of asymmetries in the emission from the disc, and therefore it would evidence relativistic effects.
2.1 The Compton–thick neutral absorber
Both low column densities and large ionisations are associated with more transparent clouds, and therefore with fainter CFs. Since our goal is to study the shape of the CF–profiles in order to look for asymmetries pointing to relativistic effects, we decided to explore the parameters space of an absorber able to produce noticeable coverings of the flux emitted by the accretion disc. Hence, the parameters of the absorber that we consider are those of an unphysical opaque neutral cloud, i.e. cm and log . The study of realistic clouds (Compton–thin and Compton–thick up to cm, and ionisations of up to log ) is beyond the purpose of this paper and will be detailed elsewhere (Sanfrutos et al. in preparation).
The illumination profiles of the disc play an important role too: the steeper the emissivity index is, the more centrally concentrated is the corona and the flux is reflected by material within regions closer to the central engine. Therefore, the higher is the emissivity index parameter the deeper is the eclipse and the larger the CF. Below we present results from some CF–profiles computed with an intermediate emissivity index of . The CF–profiles are compiled in Fig. 1.
The first parameter about which we can establish some conclusions is the inclination of the system, in terms of the angle between the normal to the disc and our LOS. When we observe the system face–on (i.e. ), Doppler boosting is undetectable from our point of view: there are no approaching nor receding regions in the disc, since its orbital plane is perpendicular to our LOS. Therefore the CF–profile observed would be symmetric with respect to the position of the black hole. The larger the inclination is, the greater is the Doppler boosting effect, so that the more noticeable is the asymmetry of the CF–profiles, as can be seen in every panel of Fig. 1, independently of every other parameter. The CF–profiles during an eclipse in a system at an inclination of (near to pole–on) are shown in dashed red lines; at a greater inclination of in solid green lines, showing a greater asymmetry; and finally at an inclination of in dash–dotted blue lines, when the asymmetry is the largest. We do not show inclinations lower than because one would not expect to find BLR clouds at such high latitudes. We do not show inclinations larger than either because the torus would make impossible any observation.
The space of parameters was explored for three characteristic values of the spin. On the left column of Fig. 1 we show the CF–profiles of a Schwarzschild black hole (ISCO = 6 ). On the column in the middle we show the same for an intermediate Kerr black hole (ISCO = 4.23 ). On the right column we show the CFs with a maximally rotating Kerr black hole (ISCO = 1.25 ). In the cases where the cloud size is comparable to the ISCO (upper row in Fig. 1, with a cloud radius of 1 , and also in the left and even the middle panels in the second row, with a radius of 6 ), we could measure the ISCO with great precision in a completely new way when the size of the cloud is known, just by measuring the distance between the two maxima in any of those panels. When larger clouds are involved in the eclipse, general relativistic effects are still detectable, as highlights the asymmetry of the characteristic top hat–shaped CF–profiles shown on any of the panels in the third () and fourth () rows in Fig. 1.
It is important to mention that some of the maxima in the CF–profiles shown in Fig. 1 are beyond the limits of detectability of available instruments (mainly the upper row cases, when the cloud radius is too small). However, most of them still represent interesting cases that we could discern by using current technology.
3 Spectral analysis
We focus now on the case of an intermediate Kerr black hole seen at an inclination of and obscured by two different clouds: one () slightly larger than the ISCO and one slightly smaller (). We choose a relatively steep emissivity index in order to get a more centrally concentrated emission from the accretion disc, which will allow the occultation event to reach noticeable CFs with small–sized clouds. Such steep emissivities or even larger were reported in some of the AGN, e.g. 1H 0707-495 (Fabian et al. 2009) or IRAS13224-3809 (Ponti et al. 2010).
In order to check if the existence of such a cloud is physically possible attending to the physical constraints of the X–ray source and absorber sizes and the distance between them, we assume that the central SMBH mass is , hence its gravitational radius is cm. With a cloud of column density cm, this implies densities of cm and cm respectively, which are typical values for the clouds in the BLR. Assuming an orbital velocity of km s, as already mentioned in Section 2, the cloud would be located at a distance of cm from the SMBH, and would travel the region of our study () in ks, at a velocity of ks.
3.1 Simulated spectra
Next we generate the corresponding CF–profiles for the situation under study, which are shown in the upper panel of Fig. 2. Beyond the covering fraction profiles, it is more interesting to determine whether a cloud involved in this kind of eclipses could imprint any effect on the spectra taken all along one obscuration event, and if so, discern if it would be detectable by current or future observatories. To this extent, we perform several 10–ks simulations of spectra as they would be observed with the EPIC-pn instrument on board XMM–Newton during different stages of an eclipse caused by a BLR cloud of and . We show the most characteristic ones in the middle and bottom panels of Fig. 2.
In both cases, the spectra out of the eclipse (at ) are shown in black X–shaped markers. The maximum absorption occurs when the cloud covers the approaching part of the inner accretion disc, at in the case of the smallest cloud and at in the case of the biggest, shown in blue circles. The hard band ( keV) remains unabsorbed during all the event, while the soft band ( keV) shows variability. Since the radius of the small cloud is lower than the ISCO, when the cloud is around the position in the reference grid with respect to the black hole, it does not cover any flux, hence the spectrum is again unabsorbed (not shown for clarity, but consistent with the black spectrum in the middle panel of Fig. 2). However, when the cloud is larger than the ISCO, it intercepts flux from both the approaching and the receding parts of the disc, and absorption is still notable (though lower than when the cloud was in the position), as shown in green squares in the bottom panel of Fig. 2. When the cloud is over the receding part of the disc (), absorption enhances as expected in the case of the small cloud, as shown in red triangles in the middle panel of Fig. 2, although not as much as when the approaching part was eclipsed. In the case of the biggest cloud, when it occupies the same position absorption does not enhance since the cloud is larger than the ISCO and covers part of the inner accretion disc all the time. Finally, as the cloud goes away the spectra recover their initial flux values at all energy ranges (not shown for clarity reasons, but consistent with the black spectra). From this analysis, specifically of the small cloud transit, we conclude that even though the CF is low, both eclipses are deep enough to be detected spectroscopically.
3.2 Model–independent H / S light curves
It is even more interesting to compare the different behaviour between the soft and the hard bands: while the clouds are transparent to hard photons, part of the soft ones are absorbed. Therefore these kind of obscuration events drive absorption variability in the soft X–ray band. With the aim of better understanding how this effect takes place, we simulate the H / S ratios too. The process is as follows: from the spectra generated in Section 3.1 we compute the count rates within the hard band in the keV range, and also in the soft band in the keV range, and then we compute the H / S ratio light curve of each observation simply by computing the fraction of the hard band over the soft one. The light curves are shown in Fig. 3, in terms of the net count rates computed from the simulated spectra (red small circles for the –sized cloud, green large circles for the –sized cloud), together with the model–predicted count rates (dotted and dashed lines respectively). The advantage of this technique is that the eclipse can be identified and characterised just computing the total count rates in different bands of each one of the spectra, in a completely model–independent way. By comparing Fig. 2 and Fig. 3, we conclude that an obscuration event is associated with the existence of high H / S ratio regions in the light curves. In addition, the H / S light curves and the CF profiles share the same shape, so that the eclipse parameters could be measured in a model–independent way studying the H / S ratio of several short spectra all along one observation. As an example, from the CF–profile in Fig. 2, we find the distance between the two maxima to be , which corresponds with a time interval of ks for a cloud moving at 3000 km s. From the H / S light curve we find this interval to be ks for the same cloud, which is consistent with the value computed from the model–dependent CFs. Designing a finer–scale sampling is not possible for the chosen set of parameters, but it would not be an issue to do it for a system with a slower cloud or a larger black hole mass. We let the study of such systems for a more exhaustive future paper (Sanfrutos et al. in preparation).
We have shown some promising preliminar results. In Section 2 we defined a model for the continuum X–ray emission plus the reflection component convolved with the relativistic kernel KYNCONV (Dovčiak et al. 2004), giving account of GR effects due to the proximity of the material in the inner accretion disc to the SMBH. The eclipse of these regions by a Compton–thick neutral cloud produces peculiar CF–profiles depending on the inclination of the system, on the black hole spin and emissivity index, and on the size of the absorber. The study of such profiles allows us to discern how much flux is emitted from every region within the disc. This way we find that the emission measured from our point of view is anisotropic: it is much larger from the parts of the disc that are running towards us in their orbital motion around the SMBH than from its receding parts. This is an effect of Doppler boosting, and its detection would definitely prove that GR effects imprint the AGN X–ray spectra. We may have already detected these disc emission anisotropies during some of these BLR clouds eclipses, since the shapes of the CF–profiles shown in Fig. 1 are remarkably similar to those shown in Fig. 6 in Risaliti et al. (2009) and Fig. 6 in Sanfrutos et al. (2013). This encourages us to deepen the study of these relativistic effects by fitting simulated and real spectra, with the caveat of the detectability limit.
In Section 3.1 we explored the detectability of GR effects in a set of simulated X–ray spectra along one occultation event, finding that (i) a Compton-thick ( cm) neutral cloud is transparent to hard X–rays ( keV), but produces absorption in the soft band ( keV); and that (ii) the maximum absorption occurs when the cloud intercepts our LOS to the approaching part of the accretion disc, as expected due to Doppler boosting. The difference in absorption between the hard and the soft bands led us to study the H / S ratio light curves produced during an eclipse. In Section 3.2 we showed that the behaviour of the H / S light curves follows that of the CF–profiles, hence we concluded that it is possible to characterise the properties of the innermost regions of the accretion disc in a new model–independent way just by computing the count rates of the hard and soft bands of spectra taken all along an eclipse.
Only XMM–Newton simulations were studied at the moment. In a forthcoming work (Sanfrutos et al. in preparation) we will fit real XMM–Newton data and will also simulate spectra from other instruments. In this context we will study the effect of eclipses produced by thicker and / or ionised clouds in harder bands by means of NuSTAR. We will also perform simulations of spectroscopic data from Athena, considering more complex models giving account of resolved spectral absorption lines (PHASE code by Krongold et al. 2003) in order to take advantage of the outstanding spectral resolution of the future X–ray observatory.
Finally, among the technical issues to take into account, it will be important to refine any possible degeneracies among the parameters, to study spectra in several bands of interest in the XMM–Newton archive, and to study variations in the soft excess band, where the sensitivity of our instruments gets its highest.
Financial support for this work was provided by the European Union through the Seventh Framework Programme (FP7/2007–2013) under grant n. 312789. MS thanks CSIC for support through a JAE–Predoc grant, and rejects cuts to public science budget.
- Bianchi et al. (2009) Bianchi, S., Piconcelli, E., Chiaberge, M., et al. 2009, ApJ, 695, 781
- Dovčiak et al. (2004) Dovčiak, M., Karas, V., & Yaqoob, T. 2004, ApJS, 153, 205
- Elvis et al. (2004) Elvis, M., Risaliti, G., Nicastro, F., et al. 2004, ApJ Letters, 615, L25
- Fabian et al. (2015) Fabian, A. C., Lohfink, A., Kara, E., et al. 2015, MNRAS, 451, 4375
- Fabian et al. (2009) Fabian, A. C., Zoghbi, A., Ross, R. R., et al. 2009, Nature, 459, 540
- García et al. (2013) García, J., Dauser, T., Reynolds, C. S., et al. 2013, ApJ, 768, 146
- Krongold et al. (2003) Krongold, Y., Nicastro, F., Brickhouse, N. S., et al. 2003, ApJ, 597, 832
- Leighly (1999) Leighly, K. M. 1999, ApJS, 125, 317
- Maiolino et al. (2010) Maiolino, R., Risaliti, G., Salvati, M., et al. 2010, A&A, 517, A47
- Marshall et al. (1981) Marshall, N., Warwick, R. S., & Pounds, K. A. 1981, MNRAS, 194, 987
- McKernan & Yaqoob (1998) McKernan, B. & Yaqoob, T. 1998, ApJ Letters, 501, L29
- Miniutti et al. (2014) Miniutti, G., Sanfrutos, M., Beuchert, T., et al. 2014, MNRAS, 437, 1776
- Ponti et al. (2010) Ponti, G., Gallo, L. C., Fabian, A. C., et al. 2010, MNRAS, 406, 2591
- Puccetti et al. (2007) Puccetti, S., Fiore, F., Risaliti, G., et al. 2007, MNRAS, 377, 607
- Reeves et al. (2008) Reeves, J., Done, C., Pounds, K., et al. 2008, MNRAS, 385, 108
- Reynolds & Nowak (2003) Reynolds, C. S. & Nowak, M. A. 2003, PhR, 377, 389
- Risaliti et al. (2009) Risaliti, G., Miniutti, G., Elvis, M., et al. 2009, ApJ, 696, 160
- Risaliti et al. (2011a) Risaliti, G., Nardini, E., Elvis, M., Brenneman, L., & Salvati, M. 2011a, Monthly Notices of the Royal Astronomical Society, Volume 417, Issue 1, pp. 178-183., 417, 178
- Risaliti et al. (2011b) Risaliti, G., Nardini, E., Salvati, M., et al. 2011b, MNRAS, 410, 1027
- Rozanska et al. (2015) Rozanska, A., Malzac, J., Belmont, R., Czerny, B., & Petrucci, P. O. 2015, ArXiv e-prints
- Sanfrutos et al. (2013) Sanfrutos, M., Miniutti, G., Agís-González, B., et al. 2013, MNRAS, 436, 1588
- Stern et al. (1995) Stern, B. E., Poutanen, J., Svensson, R., Sikora, M., & Begelman, M. C. 1995, ApJ Letters, 449, L13
- Turner et al. (2007) Turner, T. J., Miller, L., Reeves, J. N., & Kraemer, S. B. 2007, A&A, 475, 121
- Wilkins & Gallo (2015) Wilkins, D. R. & Gallo, L. C. 2015, MNRAS, 449, 129