Acta Polytechnica doi:10.14311/AP.2014.54.0301 Acta Polytechnica 54(4):301–304, 2014 © Czech Technical University in Prague, 2014 available online at http://ojs.cvut.cz/ojs/index.php/ap BEYOND THE STANDARD MODEL OF THE DISC–LINE SPECTRAL PROFILES FROM BLACK HOLE ACCRETION DISCS Vjačeslav Sochora∗, Vladimír Karas, Jiří Svoboda, Michal Dovčiak Astronomical Institute, Academy of Sciences of the Czech Republic, Boční II 1401, Prague, Czech Republic ∗ corresponding author: sochora@astro.cas.cz Abstract. The strong gravitational field of a black hole has distinct effects on the observed profile of a spectral line from an accretion disc near a black hole. The observed profile of the spectral line is broadened and skewed by a fast orbital motion and redshifted by a gravitational field. These effects can help us to constrain the parameters of a system with a black hole, both in active galactic nuclei and in a stellar-mass black hole. Here we explore the fact that an accretion disc emission can be mathematically imagined as a superposition of radiating accretion rings that extend from the inner edge to the outer rim of the disc, with some radially varying emissivity. In our work, we show that a characteristic double-horn profile of several radially confined (relatively narrow) accretion rings or belts could be recognized by the planned instruments onboard future satellites (such as the proposed ESA Large Observatory for X-ray Timing). Keywords: black hole physics, accretion discs, galactic nuclei. 1. Introduction An observation of spectral lines from the inner regions of an accretion disc around a black hole, both in active galactic nuclei (AGN) [1, 2] and in Galactic black holes [3, 4], gives us information about matter in extreme conditions. These spectral lines are broadened and skewed by a fast orbital motion and redshifted by a strong gravitational field. According to the standard scenario [5, 6], line emis- sivity is assumed to be a simple power-law of the radius. With a typically moderate inclination angle of the source, a broad profile is formed with an ex- tended red wing and a dominant well-defined blue peak. However, the radial emissivity of an astrophysi- cally realistic accretion disc cannot be a simple smooth function of the radius. Instead, it is expected to have peaks of enhanced emissivity occurring at particular radii, e.g. due to localized irradiation by magnetic flares [7, 8]. We address the question whether the emission ex- cesses on top of the standard emission profile can be resolved in observed spectra and used to further constrain the black hole spin to better precision. We discussed in [9] whether the proposed Large Obser- vatory for X-ray Timing (LOFT), [10, 11], will have the necessary capability to reconstruct the parameters from a model spectrum. We produced artificial data with appropriate properties and then we analyzed them by using a preliminary response file for LOFT. Two scientific instruments form the payload of the satellite: LAD (Large Area Detector) with a large effective area (designed to reach ' 12 m2) and the energy resolution should be about 200–300 eV; and WFM (Wide Field Monitor), which will observe about 50 % of the sky available to the LAD in the same energy band at any time. In this paper, we compare the results obtained in [9] and discuss how variability of the background can affect the spectrum. Correct modeling of the back- ground is crucial for the success of the measurements, because LOFT does not contain any telescope that could measure the background from the neighborhood of the observed object. 2. Model spectrum and analysis by LOFT 2.1. Test case We took our fiducial model (Figure 1, left panel) from [9]: it was phabs * (powerlaw + 4 * kyrline), i.e., a photo-absorbed power-law continuum and four line components blurred by relativistic effects (we used XSPEC v. 12.6.0). One of the kyrline components originates over the entire disc surface, and it has been fixed to its default parameters (rISCO ≤ r ≤ 400, radial emissivity index α = 3). We set the model parameters to: a = 0.93 (rapidly spinning Kerr black hole in prograde rotation), i = 30 deg (moderate inclination typical for the Seyfert 1 nucleus) and three rings of width 0.5rg at the positions r1 = 3rg, r2 = 4rg and r3 = 6rg. We produced the simulated spectrum (Figure 1, right panel) by assuming a source flux of approximately 1.3 mCrab (' 3 × 10−11 erg/cm2 in the energy range 2–10 eV), a photo-absorbed power-law continuum (photon index Γ = 1.9, hydrogen density nH = 4 × 1021 cm−2) and the rest energy Erest = 6.4 keV. The exposure time was set to 20 ksec. 301 http://dx.doi.org/10.14311/AP.2014.54.0301 http://ojs.cvut.cz/ojs/index.php/ap V. Sochora, V. Karas, J. Svoboda, M. Dovčiak Acta Polytechnica Figure 1. Left panel: The complete theoretical model and the model components: a power-law continuum and the individual line profiles from which the energy shifts of the components are derived. Right panel: Simulated data and the ratio to the baseline model consisting of the power-law and the disc-line components. Residuals related to the three additional narrow rings are clearly visible. Taken from [9]. Ring gmin gmax rin rout a = 0.76 a = 1.00 a = 0.76 a = 1.00 1 0.36 0.81 3.1 2.8 3.7 3.4 2 0.48 0.91 4.1 3.9 4.9 4.7 3 0.59 0.98 5.8 5.6 7.1 6.9 Table 1. Parameters of the model inferred from the energy positions of the spectral peaks in the test spectrum from Figure 1. We identified the visible features with the horns of the line components. We imposed the same inclination i = 30 deg for all three rings and required the inferred spin values to be consistent with each other. The spin turns out to be constrained only partially, with the values from 0.77 up to 1.00 being consistent with the positions of peaks in the model spectrum when the radius is set appropriately. The fiducial test spectrum was generated for rings positioned at radii rin = 3 rg, 4 rg, and 6 rg, respectively. The tabulated values demonstrate the accuracy of the fitting procedure. See the text for details. Taken from [9]. 2.2. Determination of shifts from the spectral profile To determine the relativistic energy shifts of pho- tons, we adopt the method described in our recent paper [12], where we considered the propagation of photons from the source in the limit of geometrical optics in the Kerr metric [13]. There is a partial degeneracy of the parameter values. In our case this exhibits itself by the fact that, in order to ob- tain the red peaks of the line in the right position, the spin has to be greater than the lower limit of a = 0.76. However, the upper bound remains unde- termined. For 0.76 ≤ a ≤ 1, i.e., up to the maximum spin of the Kerr black hole, we can reproduce the peaks by rearranging the ring radii. This is shown in Table 1 by giving two possible values of rin and rout that are consistent simultaneously with the men- tioned minimum and maximum spin values. One can see that the uncertainty in the inferred radii is be- low 10%, while for spin the relative error represents about 25%. Figure 2. The LAD background and its various contributions. The results of a study in which the behavior of each background component was modeled along the satellite orbit. The spectrum of a 10 mCrab source is also shown. [14] 302 vol. 54 no. 4/2014 Beyond the Standard Model of the Disc–Line Spectral Profiles 2 5 0 .9 5 1 1 .0 5 1 .1 1 .1 5 ra ti o Energy (keV) data/model 2 5 0 .9 5 1 1 .0 5 1 .1 ra ti o Energy (keV) data/model 2 5 0 .9 5 1 1 .0 5 1 .1 ra ti o Energy (keV) data/model Figure 3. The ratio between the simulated data and the baseline model with a background level of 5% (top), 7% (middle) and 10% (bottom). 0.85 0.9 0.95 1 2 .8 3 3 .2 ri n ( G M /c 2 ) a/M (GM/c) + 0.85 0.9 0.95 1 2 .9 3 3 .1 3 .2 ri n ( G M /c 2 ) a/M (GM/c) + 0.85 0.9 0.95 1 2 .8 2 .9 3 3 .1 ri n ( G M /c 2 ) a/M (GM/c) + Figure 4. Constraints on the best-fit model with a background level of 5% (top), 7% (middle) and 10% (bottom). 303 V. Sochora, V. Karas, J. Svoboda, M. Dovčiak Acta Polytechnica 3. Model spectrum considering background 3.1. LAD background The LAD background has been analyzed and com- puted by [14] using Monte Carlo simulations of a mass model of the whole LOFT spacecraft and all known ra- diation sources in the LOFT orbit, see Figure 2. From the simulations, it is obvious that the background is dominated (> 70 %) by high energy photons of cos- mic X-ray background and Earth albedo leaking and scattering through the collimator structure, which be- comes less efficient at high energies. These two sources are stable and predictable, although there exist small modulations of these components due to the orbital motion of the satellite around the Earth. One of the varying sources is the particle induced background (< 6 % of the overall background). The largest modu- lation of the total LAD background is estimated as < 20 %, and can be effectively described by a geomet- rical model that should predict the background at the level of 1 % or better (1–20 keV). 3.2. Data to model ratio We considered the model presented in the previous section and tested the expected impact of background contamination by applying corrnorm to the back- ground file in XSPEC. Figure 3 shows the data to model ratio with the randomized background at the 5, 7 and 10 % levels. The blue peaks can be recognized in the first two graphs, while the red peaks are lost in the signal-to-noise. The 10 % inaccuracy of the back- ground degrades the visibility of the peaks. Figure 4 demonstrates the expected accuracy with which the model parameters are constrained. The results can seem to be counter-intuitive (the 10 % level is better than 5 %). This is caused by the method of fitting. In our case, we know the parameters of the system and the constraints would look similar. The study of the background deserves more investigation. The con- straints on the best-fit model parameters are derived from the simulation data. Confidence contours are shown (1, 2, and 3σ) of the inner ring radius rin vs. dimension-less spin a. 4. Conclusions In this paper, we have applied the background re- sponse file of the proposed satellite LOFT to the test model from [9], and we have studied the variability of the model spectrum. Our result is that 10 % in- accuracy is a limit value of the background, and it should not be exceeded. Accuracy of the background at the 5% level is sufficient to recognize and fit the structures in the model spectrum, and the peaks of the energy shifts, and also to determine the parameters, as described in [12]. Acknowledgements We would like to acknowledge support from GAČR 13- 00070J. References [1] A. C. Fabian, K. Iwasawa, R. C. S., A. J. Young. Broad Iron Lines in Active Galactic Nuclei. PASP 112:1145–1161, 2000. doi:10.1086/316610. [2] C. S. Reynolds, M. A. Nowak. Fluorescent iron lines as a probe of astrophysical black hole systems. Physics Reports 377:389–466, 2003. doi:10.1016/S0370-1573(02)00584-7. [3] J. M. Miller, A. C. Fabian, R. Wijnands, et al. Resolving the Composite Fe Kα Emission Line in the Galactic Black Hole Cygnus X-1 with Chandra. ApJ 578:348–356, 2002. doi:10.1086/342466. [4] J. E. McClintock, R. A. Remillard. Compact Stellar X-ray Sources. Cambridge University Press, Cambridge, 2006. [5] I. D. Novikov, K. S. Thorne. Black Holes. Gordon and Breach Publishers, New York, 1973. [6] D. N. Page, K. S. Thorne. Disk-Accretion onto a Black Hole. Time-Averaged Structure of Accretion Disk. ApJ 191:499–506, 1974. doi:10.1086/152990. [7] B. Czerny, A. Różańska, M. Dovčiak, et al. The structure and radiation spectra of illuminated accretion disks in AGN. II. Flare/spot model of X-ray variability. A&A 420:1–16, 2004. doi:10.1051/0004-6361:20035741. [8] D. A. Uzdensky, J. Goodman. Statistical Description of a Magnetized Corona above a Turbulent Accretion Disk. ApJ 682:608–629, 2008. doi:10.1086/588812. [9] V. Sochora, V. Karas, J. Svoboda, M. Dovčiak. Black hole accretion rings revealed by future X-ray spectroscopy. MNRAS 418:276–283, 2011. doi:10.1111/j.1365-2966.2011.19483.x. [10] M. Feroci, L. Stella, A. Vacchi, et al. LOFT: a large observatory for x-ray timing. In Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 7732. 2010. arXiv:1008.1009 doi:10.1117/12.857337. [11] M. Feroci, L. Stella, M. van der Klis, et al. The Large Observatory for X-ray Timing (LOFT). Experimental Astronomy 34:415–444, 2012. arXiv:1107.0436 doi:10.1007/s10686-011-9237-2. [12] V. Karas, V. Sochora. Extremal Energy Shifts of Radiation from a Ring Near a Rotating Black Hole. ApJ 725:1507–1515, 2010. doi:10.1088/0004-637X/725/2/1507. [13] S. Kato, J. Fujue, S. Mineshige. Black-Hole Accretion Disks. Kyoto University Press, Kyoto, 1998. [14] R. Campana, M. Feroci, E. Del Monte, et al. The LOFT (Large Observatory for X-ray Timing) background simulations. In Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 8443. 2012. arXiv:1209.1661 doi:10.1117/12.925999. 304 http://dx.doi.org/10.1086/316610 http://dx.doi.org/10.1016/S0370-1573(02)00584-7 http://dx.doi.org/10.1086/342466 http://dx.doi.org/10.1086/152990 http://dx.doi.org/10.1051/0004-6361:20035741 http://dx.doi.org/10.1086/588812 http://dx.doi.org/10.1111/j.1365-2966.2011.19483.x arXiv:1008.1009 http://dx.doi.org/10.1117/12.857337 arXiv:1107.0436 http://dx.doi.org/10.1007/s10686-011-9237-2 http://dx.doi.org/10.1088/0004-637X/725/2/1507 arXiv:1209.1661 http://dx.doi.org/10.1117/12.925999 Acta Polytechnica 54(4):301–304, 2014 1 Introduction 2 Model spectrum and analysis by LOFT 2.1 Test case 2.2 Determination of shifts from the spectral profile 3 Model spectrum considering background 3.1 LAD background 3.2 Data to model ratio 4 Conclusions Acknowledgements References