Issue |
A&A
Volume 522, November 2010
|
|
---|---|---|
Article Number | A64 | |
Number of page(s) | 8 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/201015207 | |
Published online | 03 November 2010 |
The broad-band X-ray spectrum of the Seyfert 1 galaxy, MCG+8-11-11
1
Dipartimento di FisicaUniversità degli Studi Roma Tre,
via della Vasca Navale 84,
00146
Roma,
Italy
e-mail: bianchi@fis.uniroma3.it
2
INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica, via U. La
Malfa 153, 90146
Palermo,
Italy
3
INAF/IASF-Roma, via del Fosso del Cavaliere, 00133
Roma,
Italy
4
INAF-IASF Bologna, via Gobetti 101, 40129
Bologna,
Italy
5
Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal
70-264, 04510
Mexico DF,
Mexico
6
Osservatorio Astronomico di Roma (INAF), via Frascati 33,
00040
Monte Porzio Catone,
Italy
Received: 14 June 2010
Accepted: 9 July 2010
Aims. Evidence is mounting that the main ingredients of the unification models of active galactic nuclei may behave differently from expectations, and be intimately related to fundamental physical parameters. The availability of high signal-to-noise broad-band X-ray spectra provides us with the opportunity to study in detail all the contributions from the materials invoked in these models, and infer their general properties, including whether their presence or absence is related to other quantities.
Methods. We present a long (100 ks) Suzaku observation of one of the X-ray brightest AGN, MCG+8-11-11. These data are complemented with the 54-month Swift BAT spectrum, allowing us to perform a broad-band fit to the 0.6–150 keV range.
Results. The fits performed in the 0.6–10 keV band provide results consistent with those of a a previous XMM-Newton observation, i.e. a lack of a soft excess, warm absorption along the line of sight, a large Compton reflection component (R ≃ 1), and an absence of a relativistic component of the neutral iron Kα emission line. However, when the PIN and Swift BAT data are included, the reflection amount drops significantly (R ≃ 0.2 − 0.3), and a relativistic iron line is required, the latter being confirmed by a phenomenological analysis in a restricted energy band (3–10 keV). When a self-consistent model is applied to the whole broad-band data, the observed reflection component appears to be entirely associated with the relativistic component of the iron Kα line.
Conclusions. The implied scenario, though strongly model-dependent, requires that all the reprocessing spectral components from Compton-thick material be associated with the accretion disc, and no evidence of a classical pc-scale torus is found. The narrow core of the neutral iron Kα line is therefore produced in Compton-thin material, such as the BLR, similarly to what is found in another Seyfert galaxy, NGC 7213, but with the notable difference that MCG+8-11-11 presents spectral signatures from an accretion disc. The very low accretion rate of NGC 7213 could explain the lack of relativistic signatures in its spectrum, but the absence of the torus in both sources is more difficult to explain, since their luminosities are comparable, and their accretion rates are completely different.
Key words: galaxies: active / galaxies: Seyfert / X-rays: general
© ESO, 2010
1. Introduction
The X-ray spectrum of unobscured AGN is characterised by ubiquitous and, generally, non-varying features, such as the neutral iron Kα narrow core and the Compton reflection component (e.g., Perola et al. 2002; Bianchi et al. 2007), whose origin appears quite clearly established within the framework of standard unification models as the pc-scale “torus” (Antonucci 1993). On the other hand, there are other spectral components, such as the warm absorber, the soft excess, and the relativistic component of the iron Kα line, which are present in only a fraction of the sources, and may vary in different observations of the same object (e.g. Piconcelli et al. 2005; Nandra et al. 2007). Their origin is less clear, and, as in the case of the relativistic iron line, may pose some crucial problems to the standard views.
X-ray emission from the accretion disc, as a result of the reprocessing of the X-ray photons produced in the corona, is a fundamental component of any AGN model, and should be characterised by a strong iron Kα line, whose profile can be clearly distinguished from the narrow core, because of well-known relativistic effects occurring near the black hole (see e.g. Fabian et al. 2000). However, this relativistic profile, when large samples of AGN are systematically analysed, is not observed in all the sources (see e.g. Nandra et al. 2007; Longinotti et al. 2008). Although at least a part of this problem must be related to the signal-to-noise ratio of the available X-ray spectra (indeed the fraction of sources with relativistic signatures appears to be larger when only sufficiently exposed X-ray spectra are analysed: see, e.g., de la Calle Pérez et al. 2010), the relativistic component of the iron Kα line still appears not to be ubiquitous, and, even when present, its equivalent width (EW) is often lower than expected.
These results imply that a revision of the standard models of accretion discs is required. Truncation of the disc, a highly ionization state, and complex illumination may all be invoked to explain the observed deviations from the current picture, but it is most fundamental to understand the ultimate drivers of this behaviour, such as the accretion rate, the luminosity, and the black hole (BH) mass. In this respect, evidence is mounting that other ingredients of the unification models may behave differently from expected, and be closely related to the above-mentioned physical parameters. For example, the broad line region (BLR) was found to be most likely absent in some sources (e.g. Bianchi et al. 2008a; Panessa et al. 2009), to act as a Compton-thick X-ray absorber in others (e.g. Elvis et al. 2004; Risaliti et al. 2005; Puccetti et al. 2007; Bianchi et al. 2009b), or even to produce most of the observed neutral iron Kα line, when no other evidence of the torus could be found in the data (Bianchi et al. 2008b).
MCG+8-11-11 (z = 0.0205) is one of the brightest AGN in the X-ray band, detected by INTEGRAL with a F20 − 100 keV ≃ 1.2 × 10-10 cgs (Beckmann et al. 2009). ASCA (Grandi et al. 1998) and BeppoSAX (Perola et al. 2000, 2002) found that the spectrum is well fitted by a pretty standard model composed of a power law, a warm absorber, a Compton reflection component, and an iron Kα line. The higher quality of the XMM-Newton spectrum (Matt et al. 2006) revealed some interesting features: the lack of a soft excess, a large reflection component, and a narrow iron line with a low EW. The presence of a large, relativistically broadened line arising from the accretion disk could also be excluded in these data.
In this paper, we present a long (100 ks) Suzaku observation of MCG+8-11-11. In addition to its 54-month Swift BAT spectrum, we explore different energy bands, in order to understand the origin of the X-ray spectral components of this bright Seyfert 1. The availability of a high signal-to-noise spectrum in such a broad band (0.6–150 keV) provides us with the opportunity to study in detail all the contributions from the materials invoked in unification models, and infer their presence/absence and general properties.
2. Data reduction
MCG+8-11-11 was observed by Suzaku on 2007, September 17, for 100 ks (obsid 702112010). X-ray Imaging Spectrometer (XIS) and Hard X-ray Detector (HXD) event files were reprocessed with the latest calibration files available (2010-03-23 release), using ftools 6.9 and Suzaku software Version 16, adopting standard filtering procedures. Source and background spectra for all the three XIS detectors were extracted from circular regions with radius of 250 pixels ( ≃ 260 arcsec), avoiding the calibration sources. Response matrices and ancillary response files were generated using xisrmfgen and xissimarfgen. We downloaded the tuned non-X-ray background (NXB) for our HXD/PIN data provided by the HXD team and extracted source and background spectra using the same good time intervals. The PIN spectrum was then corrected for dead time, and the exposure time of the background spectrum was increased by a factor of 10, as required. Finally, the contribution from the cosmic X-ray background (CXB) was subtracted from the source spectrum, simulating it as suggested by the HXD team.
We extracted the light curves of the three XIS in the soft and the hard energy band (see Fig. 1 for the XIS0). The source exhibits some variability in both bands, but the hardness ratio is practically constant within 10%. On the other hand, the HXD pin light curve does not show any significant evidence of variability. Therefore, we decided to analyse the total spectra, without any additional timing analysis.
![]() |
Fig. 1 Soft (0.3–2 keV) and hard (2–10 keV) light-curves of the Suzaku 100 ks observation of MCG+8-11-11 (only the XIS0 is plotted for clarity). The bottom panel shows the ratio between the two curves. |
The front-illuminated XIS0 and XIS3 were fitted between 0.6 and 10 keV, while the back-illuminated XIS1 was used up to 8 keV, because the background above that energy becomes significant with respect to the source spectrum for that instrument. The band 1.6–2.1 keV was excluded from all the fits, because of known inter-calibration issues at those energies due to the Si K edge. The normalizations of XIS1 and XIS3 with respect to XIS0 were left free in the fits, and always resulted in agreement to within 2%. A constant factor of 1.18 was used instead between the PIN and the XIS0, as recommended for observations taken in the HXD nominal position. In the following, all the PIN fluxes are given with respect to the XIS0 flux scale, which is 1.18 times lower than the HXD absolute flux scale.
In the following, errors correspond to the 90% confidence level for one interesting parameter (Δχ2 = 2.71), unless not otherwise stated. The adopted cosmological parameters are H0 = 70 km s-1 Mpc-1, ΩΛ = 0.73, and Ωm = 0.27 (i.e., the default ones in xspec 12.6.0: Arnaud 1996). We use the Anders & Grevesse (1989) abundances and the photoelectric absorption cross-sections by Balucinska-Church & McCammon (1992).
3. Data analysis
3.1. The iron line complex
As a first step, we fitted the 3–10 keV spectra with a simple power law. The photon index is quite flat (1.62 ± 0.02), and clear residuals remain between 6 and 7 keV (see Fig. 2). Three Gaussian lines are required by the data, at 6.401 ± 0.008, 6.66 ± 0.06, and 6.96 ± 0.02 keV, that are consistent with emission from neutral, He-like, and H-like Fe Kα. An absorption edge from neutral iron (its energy fixed to 7.11 keV: Bearden 1967) is also marginally detected, with an optical depth τ = 0.04 ± 0.02 (model I). On the other hand, a Fe Kβ emission line is not required by the data, with an upper limit of ≃ 7% with respect to the flux of the Fe Kα. In an analogous way, the Kα emission from nickel is not detected. While the two ionised iron lines are consistent with being unresolved (σ < 80 eV), the neutral Fe Kα is resolved, with eV. This width is larger than the calibration uncertainties in the XIS redistribution matrix1, and corresponds to a full width at half maximum FWHM = 6500 ± 1900 km s-1, if caused by Doppler broadening.
![]() |
Fig. 2 Left: residuals in the (rest-frame) iron line band, after a 3–10 keV fit on the three Suzaku XIS spectra, with a simple power law. Right: the same as above, but for model III, after the relativistic component of the iron line is removed (see Table 1). Data have been rebinned for clarity. |
The neutral iron Kα width can be explained by the contribution of a Compton shoulder (CS), as suggested by some residuals redwards of the line. We therefore added another Gaussian line with centroid at 6.3 keV and σ = 40 eV (both fixed: see Matt 2002). The improvement of the fit is marginal (model II: Δχ2 = 4), but the Fe Kα line is now unresolved (σ < 45 eV). However, its centroid is shifted to keV, and the flux ratio of the CS to the core is
%, which is larger than the theoretical expectations for Compton-thick matter (Matt 2002).
The width of the iron line and the redwards residuals may alternatively be signatures of a mildly relativistic line. We tested this scenario by adding to the unresolved core a kyrline model component (Dovčiak et al. 2004), adopting a fixed centroid energy at 6.4 keV, an outer radius of 400 rg, and a power-law index β = 3 for the radial dependence of the emissivity. This fit (model III) is statistically preferred to the others (Δχ2 = 15 and 11 with respect to models I and II, and a null hypothesis probability of 0.43, which is significantly better than 0.22 and 0.26, respectively). The inner radius (< 20rg) and the inclination angle (θ = 30 ± 3 deg) can be constrained in the fit, while the BH spin is completely unconstrained, and therefore fixed to be zero. The Fexxv line is no longer required by the data: if the centroid energy is allowed to vary between or be fixed to the theoretical values for each component of this multiplet (see e.g. Bianchi et al. 2005), only upper limits are recovered (of the order of 1 × 10-5 ph cm-2 s-1). The iron K edge is also an upper limit in this model (τ < 0.03). The final fit is very good (χ2 = 301 / 298 d.o.f.), and is applied in the following broad-band fits, where the presence of the relativistic component is investigated again once a more sophisticated model for the continuum is adopted. Table 1 summarizes the best-fit parameters for the models described in this section.
Best-fit model parameters for the emission features in the Suzaku XIS spectra (3–10 keV).
3.2. The 0.6–10 keV band
The best fit adopted for the 3–10 keV band (model III in Table 1) can be applied to the whole XIS band, once Galactic absorption along the line of sight is taken into account. We recover a neutral column density2 of ≃ 2.1 − 2.4 × 1021 cm-2, depending on the adopted complete model (see below and Table 2). This measure is in good agreement with the estimate of Dickey & Lockman (1990), 2.09 × 1021 cm-2, but somewhat larger than the estimate of Kalberla et al. (2005), 1.76 × 1021 cm-2. Fixing the Galactic column density to the latter value, a comparable χ2 is recovered by adding an additional neutral column density at the redshift of the source ( cm-2), without affecting any other parameter of the fit.
However, there is clearly some curvature left in the residuals at low energies, resulting in a χ2 = 681 / 515. Since additional neutral absorption at the redshift of the source is not required by the data, we added a warm absorber component (at first we adopted absori, with the temperature of the material being kept fixed to 106 K, which can be directly compared to the results reported by Matt et al. 2006), improving our fit down to χ2 = 641 / 513, by using best-fit model parameters for the warm absorber in good agreement with those found by Matt et al. (2006). As already noted by these authors, a single-zone absorber cannot properly describe the data. Some residuals remain around 0.8 keV: they can be fitted by an edge, found at keV, with τ = 0.050 ± 0.017. The energy of the edge is intermediate between K absorption from O vii (0.7393 keV) and O viii (0.8714 keV), but inconsistent with both. We note here that the energy of this feature is not even close to the absorption features observed in the XMM-Newton RGS spectra (Matt et al. 2006).
![]() |
Fig. 3 Best-fit model and Δχ2 for the three XIS spectra, with the same model as in the 3–10 keV fit, but including Galactic absorption, a warm absorber and an unidentified absorption edge (see text for details). Note the residuals above 8 keV. |
The resulting χ2 = 616 / 511 would be acceptable, since most of the residuals appears to be caused mainly by poor background subtraction (such as between 7 and 8 keV for the XIS0) and some inter-calibration issues between the instruments. However, above 8 keV the data-points appears to be systematically higher than the model (see Fig. 3). We therefore added a Compton reflection component (model pexrav, with inclination angle fixed to 30° and cutoff energy to 150 keV, see next section), which indeed improves the fit, for a final value of χ2 = 586 / 509. The photon index is now steeper (Γ = 1.772 ± 0.017) than the one found in the previous section for the limited band, and the amount of Compton reflection is R = 0.9 ± 0.3, with an iron abundance of . The need for a relativistic component of the neutral iron line is far less significant, given the large inner radius (
rg) and the small EW (30 ± 12 eV). To understand whether the disappearance of the relativistic line and the other fitting parameters depend strongly on our choice of a simple model for warm absorption, we tried a more sophisticated warm absorber model with respect to absori, which only takes into account absorption edges.
![]() |
Fig. 4 Left: the best-fit model for the 0.6–10 keV band, when extrapolated to the band covered by the PIN spectrum. Right: the Suzaku+Swift BAT best-fit model in the whole 0.6–150 keV band. |
Best-fits for the Suzaku XIS spectra (0.6–10 keV, first two columns) and the Suzaku+Swift BAT spectra (0.6–150 keV, last two columns).
We prepared an ad-hoc table produced with the photoionization code cloudy C08.00 (Ferland et al. 1998). The model ingredients are: (1) a multi-component continuum mimicking a typical AGN spectral energy distribution (SED), with αox = 1.26 and Γ2 − 10keV = 1.8, as extracted from the results of the simultaneous UV and X-ray data of the XMM-Newton observation of MCG+8-11-11 (Matt et al. 2006); (2) constant electron density3ne = 104 cm-3; (3) ionization parameter4 in the range log U = − 2.0:2.5; (4) intervening column density in the range log NNH = 20.0:22.0. The resulting fit is marginally better than the one found with absori (Δχ2 = 5) and we note that the unidentified absorption edge is still required by the data, with parameters consistent with those found in the other fit. However, the properties of the warm absorber gas are quite different, both the column density and the ionization parameter being significantly lower in the cloudy best fit, confirming that the signatures of warm absorption are weak in this source.
All the other best-fit parameters are indeed consistent with those found with absori, including the narrowish and weak relativistic component of the Fe Kα emission line. The fit is improved neither when a velocity shift is allowed for the warm absorbing gas, nor with the addition of a second warm absorber component. The first two columns of Table 2 summarize and compare the results of these two fits.
3.3. The broad-band fit
When the best-fit model found for the 0.6–10 keV XIS spectra is extrapolated to higher energies, the PIN data points all fall well below the expectations (see Fig. 4). Allowing the model to be adjusted for the high-energy data, the best-fit model is acceptable (χ2 = 625 / 524), the main difference with the limited band fit being the very low value of the Compton reflection component . The relativistic component of the iron Kα line is now very significant, recovering the same properties found in model III of the 3–10 keV fit (see Sect. 3.1 and Table 1). The cutoff energy cannot be measured, with a 99% lower limit of ≃ 130 keV (for two interesting parameters, when we assume the reflection amount to be the second interesting parameter).
To constrain the cutoff energy, we then added the Swift BAT spectrum (in the range 15–150 keV), from the 54-month source catalog (Cusumano et al., A&A submitted). The normalization factor between the XIS0 and the BAT spectrum was allowed to vary, to take into account source variability during the 54 months of integration. The best-fit value of this factor (0.76 ± 0.03) indeed clearly suggests that the Suzaku observation caught the source in a flux state significantly higher than the average of the BAT campaign. This is confirmed by the BAT lightcurve, which show significant variability in the 54 months, and a flux comparable to the one measured by the PIN at the date of the Suzaku observation (see Fig. 5). After assuming that this flux variability does not imply a spectral variability, we refitted the best-fit model including the BAT spectrum, and found that keV. This value, which remains consistent with the lower limit found with the Suzaku data only, is also in agreement with that measured by BeppoSAX (
keV: Perola et al. 2000) and only marginally inconsistent with OSSE (
keV: Grandi et al. 1998).
![]() |
Fig. 5 The 15–50 keV Swift BAT 54-month light curve of MCG+8-11-11 (a bin is 15 days). The countrate was converted to a 15–50 keV flux, assuming no spectral variation. The Suzaku PIN flux is shown in comparison, with its uncertainty, with the dotted lines and the red symbol. The agreement at the date of the Suzaku observation is good. |
In contrast to the 0.6–10 keV fit discussed in the previous section, the best-fit model for the broad-band 0.6–150 keV spectra is achieved with the warm absorber parametrised by absori: a significantly poorer fit is found when the cloudy table is used (Δχ2 = 18: see Table 2). In any case, all the other parameters appears to be independent of the choice of the warm absorber. In particular, both fits recover a low reflection component (R ≃ 0.2 − 0.3) and an iron abundance significantly lower than the solar one. Their inter-correlation and dependence upon other relevant parameters are shown in Fig. 6. The relativistic component of the iron line is now interestingly very significant, and its best-fit model parameters are similar to those found in our first phenomenological fits of the iron line complex (see Sect. 3.1 and Table 1).
The broad-band best fit is not as statistically good as the one constrained by the XIS band-pass, being χ2 = 640 / 531. Most of the residuals are due to the XIS data-points above 8 keV, which are on average higher than the model: as we have seen in the previous section, they were most likely the main reasons for requiring an high Compton reflection in the 0.6–10 keV band, but they are apparently inconsistent with the PIN flux. We alternatively tried to keep the iron abundance fixed to the solar value, leaving the inclination angle of the reflecting material free in the fit. The resulting χ2 is comparable to the best fit (χ2 = 646 / 531) with a very high inclination angle (cosi < 0.30). The reflection fraction is now much larger, but loosely constrained ().
4. Discussion
4.1. Comparison with previous observations
In the past 15 years, MCG+8-11-11 has been observed in the X-rays by ASCA (simultaneously with OSSE, Grandi et al. 1998), BeppoSAX (Perola et al. 2000), and XMM-Newton (Matt et al. 2006). Despite significant flux variability (a factor ≃ 3 between ASCA and the present Suzaku observation), there is no clear evidence of variations of the spectral shape of the source. The photon index of the primary continuum was also measured in the range 1.7–1.8, a difference of 0.05 arising in this Suzaku observation whether the PIN data are considered or not in the fit. The warm absorber, if modelled in the same way (i.e. absori) also gives consistent results between XMM-Newton and Suzaku. The notable absence of a soft excess is confirmed in all the observations, even if the large Galactic column density may help to conceal it in the data.
The iron Kα neutral line has a complex profile, and is probably composed of different components. Both ASCA and BeppoSAX data suggested a broadish profile (σ ≃ 0.2 − 0.3 keV, but marginally resolved) with EWs consistent with a constant line flux (Grandi et al. 1998; Perola et al. 2000). The same applies to the XMM-Newton and the Suzaku iron Kα flux, once all the components required by the data are taken into account. However, while these results suggest that the iron Kα flux did not vary significantly, the errors and the differences in modelling its profile make it difficult to extract any other useful information from this apparent lack of response to the primary flux variability.
![]() |
Fig. 6 Suzaku+Swift BAT best fit contour plots (black, red, green lines correspond to 68, 90, and 99% confidence level for two interesting parameters). Left: Compton reflection amount R vs. cutoff energy Ec; middle: Compton reflection amount R vs. iron abundance AFe; right: Compton reflection amount R vs. spectral photon index Γ. |
As for the Compton reflection component, all previous observations are consistent with values around R ≃ 1 − 1.5, but with a good (statistical) constraint only on the XMM-Newton data (Matt et al. 2006). The Suzaku best-fit model limited to the 0.6–10 keV band is consistent with these values, but still in agreement, within the errors, with a scenario where the reflection component flux is constant, and the value of R anticorrelates with the flux of the primary continuum, as expected if the reprocessing originates in a distant material, such as the torus. However, when the PIN data are taken into account, the Compton reflection drops to a significantly lower value, around 0.2–0.3 (see Table 2). This is also supported by the preference for very low iron abundance and/or large inclination angles with respect to the line-of-sight, probably caused by the iron edge being not particularly deep, as confirmed by the phenomenological fit we performed in the 3–10 keV band (see Sect. 3.1), where the edge is significant only when the relativistic component of the iron Kα line is not included in the model. We note that an iron abundance lower than solar is also preferred in the 0.6–10 keV fit, as found by Matt et al. (2006) in the XMM-Newton data. In these cases, though, the relativistic component of the iron Kα line is completely negligible.
4.2. A global scenario
To test the data with a self-consistent global scenario, we replaced in the model described in Sect. 3.3 (see also Table 2) the pexrav and the narrow core of the iron Kα line with the pexmon model (Nandra et al. 2007), which includes at the same time the incident primary power law, and all the reprocessing components produced by Compton-thick material: the emission features (Fe Kα, Kβ and Compton shoulder, and the Ni Kα emission line) and the Compton reflection hump, with all the relevant inter-dependencies (inclination, index, and abundance). The resulting fit is worse than the best fit, but still acceptable (χ2 = 650 / 532). The iron abundance is now much higher (), and the amount of reflection is somewhat higher (
). This is due to the need to model the narrow component of the iron neutral Kα line self-consistenly with the Compton reflection hump. On the other hand, all the other parameters are consistent with those found in our best fit presented in the previous section, within the errors, including the warm absorber (absori) and the absorption edge.
We then tested whether an additional component of the neutral iron line narrow core is needed by the data. A possibility is a contribution from the BLR. In the case of MCG+8-11-11, the FWHM of the Hα is 2920 km s-1, and that of the Hβ 3630 km s-1 (Osterbrock & Shuder 1982). We therefore added a neutral Fe Kα line with σ fixed at 30 eV (consistent with a production in the same gas in which Hα and Hβ originate), and fixed the iron abundance at the Solar value. The fit only slightly improves (χ2 = 647 / 532), due to the marginal significance of the iron Kα line supposedly produced in the BLR (EW = 14 ± 12 eV). This test appears to suggest that the observed iron Kα line can be satisfactorily accompanied by the measured low reflection component, in a self-consistent way. But this is due to only the narrow core needing to be modelled, with EW = 60 − 70 eV (see Table 2).
Indeed, in all these fits, a significant contribution from a relativistic component to the iron line is still required. This component should be accompanied by its own Compton reflection, produced in the accretion disc. To test this scenario, we therefore added another pexmon component (only the reprocessing components), convolved with a relativistic kernel (kdblur) to model the relativistic effects occurring in the inner accretion disc, close to the BH. The final best-fit model improves, and it is almost statistically equivalent to the best fit found in Sect. 3.3 (χ2 = 644 / 532). All the reflection component originates in the disc, with , while the contribution from the torus is negligible (R < 0.1). Therefore, all the narrow component of the iron Kα line is now modelled by emission from the BLR (σ = 30 eV as above), with an EW ≃ 70 eV. The inner radius of the accretion disc is just beyond the innermost stable orbit for a non-spinning BH (
).
The resulting scenario is very interesting. If the reflection component originates entirely in the accretion disc, consistently with the observed relativistic component of the iron Kα line, it is unavoidable that the narrow core of the emission line is produced in Compton-thin material, such as the BLR, whose contribution to the reflection amount would be negligible. This, in turns, implies that there would be no evidence of the classical Compton-thick torus in this source. Given the high-flux state of the source during the Suzaku observation (see Sect. 3.3 and Fig. 5), it is possible that the reprocessing components from the torus are diluted by the stronger primary continuum. Although this would justify the low (or absent) Compton reflection from the torus, most of the narrow iron Kα line would still be produced in Compton-thin material, because the component originating from the torus would be diluted.
This situation is similar to the one observed in another Seyfert galaxy, NGC 7213, the only Seyfert 1 galaxy with a negligible amount of Compton reflection (Bianchi et al. 2003, 2008b; Lobban et al. 2010). The observed iron neutral line in this source (EW ≃ 100 eV) was explained by Bianchi et al. (2008b) as being produced mainly in the BLR, as confirmed by the consistency between the FWHM of the broad component of the Hα line and the neutral iron Kα line, in quasi-simultaneous optical/Chandra observations. The same scenario would apply to MCG+8-11-11, but with the notable difference that this object exhibits spectral signatures of an accretion disc.
We found two BH mass estimates for MCG+8-11-11 in the literature, both derived by means of the BLR radius vs. 5100 Å relationship, and the Hβ, but using different data, i.e. 1.5 × 107 M⊙ (Bian & Zhao 2003) and 1.2 × 108 M⊙ (Winter et al. 2010). Applying a bolometric correction of 30 to the observed 2–10 keV luminosity (as appropriate in this luminosity range, according to Marconi et al. 2004), we recover Lbol/Ledd ≃ 0.1 − 1. Although a better BH mass estimate is needed to assess the reliability of this accretion rate, the difference with respect to NGC 7213 is clear, being it characterised by a low accretion rate (≃ 3 × 10-3: Bianchi et al. 2009a). This could explain the lack of relativistic signatures in the latter, where the accretion disc may be truncated and its inner part replaced by some sort of inefficient disc (see e.g. Lobban et al. 2010). The absence of the torus in both sources is instead more puzzling. The formation of the torus may be suppressed at low accretion rates and/or luminosities, as predicted for the BLR (e.g. Nicastro 2000; Elitzur & Shlosman 2006). However, as we have already seen, the accretion rates of these two sources have nothing in common, and their luminosities are by no means low.
5. Conclusions
We have presented a long Suzaku observation of one of the X-ray brightest AGN, MCG+8-11-11. The spectrum is characterised by a standard power-law photon index of 1.7–1.8, and a warm absorber which cannot be satisfactorily modelled with a single gas component. The notable absence of a soft excess, which is nearly ubiquitous in unobscured Seyfert galaxies, may be partly due to the high Galactic column density along the line-of sight.
The fits performed in the 0.6–10 keV band provide results consistent with a previous XMM-Newton observation, i.e. a large Compton reflection component (R ≃ 1) and the absence of a relativistic component of the neutral iron Kα emission line. However, when the PIN data are included, the reflection amount drops significantly (R ≃ 0.2 − 0.3), and a relativistic iron line is required, the latter confirmed by a phenomenological analysis in a restricted energy band (3–10 keV). The addition of the Swift BAT 54-month spectrum up to 150 keV allowed us to measure an high-energy cutoff at 150 keV, in agreement with previous results.
When a self-consistent model is applied to the whole data, the observed reflection component appears to be associated entirely with the relativistic component of the iron Kα line. Therefore, all reprocessing from Compton-thick material must be associated with the accretion disc, and no evidence for the classical pc-scale torus is found. It is unavoidable that the narrow core of the neutral iron Kα line is produced in Compton-thin material, such as the BLR, whose contribution to the reflection amount would be negligible. Although more model-dependent, this situation is similar to the one observed in another Seyfert galaxy, NGC 7213, the only Seyfert 1 galaxy with a negligible amount of Compton reflection, but with the notable difference that MCG+8-11-11 exhibits spectral signatures of an accretion disc. The very low accretion rate of NGC 7213 could explain the lack of relativistic signatures in its spectrum, but the absence of the torus in both sources is more difficult to explain, since their luminosities are comparable, and their accretion rates are completely different.
We extracted the calibration spectra for this observation, and fitted them separately for each XIS with a model consisting of three emission lines: the Mn Kα doublet (the separation of the lines fixed to be 0.0111 keV: Bearden 1967) and the Kβ. The width of the three lines were kept linked. We recovered upper limits in all the spectra: 23 eV (XIS0), 25 eV (XIS1), 14 eV (XIS3).
This value of column density was found by adopting the tbabs model, and the abundances and cross-sections minimizing the χ2 (see Sect. 2 for the final choice).
Models with higher electron densities up to ne = 109 cm-3 have also been produced, but all the following fits were found to be insensitive to this parameter, as expected (see e.g. Nicastro et al. 1999).
We use the standard definition where c is the speed of light, r the distance of the gas from the illuminating source, ne the electron density, and νR the frequency corresponding to 1 Rydberg (Osterbrock & Ferland 2006).
Acknowledgments
This research has made use of data obtained from the Suzaku satellite, a collaborative mission between the space agencies of Japan (JAXA) and the USA (NASA). We thank the anonymous referee for useful suggestions which improved the clarity of the paper. S.B., E.P. and G.M. acknowledge financial support from ASI (grant I/088/06/0). We would like to thank K. Nandra, G. Miniutti, K. Arnaud and C. Gordon for the xspec model pexmon, and K. Hamaguchi for his support on Suzaku data reduction and analysis.
References
- Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 [Google Scholar]
- Antonucci, R. 1993, ARA&A, 31, 473 [Google Scholar]
- Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, ASP Conf. Ser., 101, 17 [Google Scholar]
- Balucinska-Church, M., & McCammon, D. 1992, ApJ, 400, 699 [NASA ADS] [CrossRef] [Google Scholar]
- Bearden, J. A. 1967, Rev. Mod. Phys., 39, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Beckmann, V., Soldi, S., Ricci, C., et al. 2009, A&A, 505, 417 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bian, W., & Zhao, Y. 2003, MNRAS, 343, 164 [NASA ADS] [CrossRef] [Google Scholar]
- Bianchi, S., Matt, G., Balestra, I., & Perola, G. C. 2003, A&A, 407, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bianchi, S., Matt, G., Nicastro, F., Porquet, D., & Dubau, J. 2005, MNRAS, 357, 599 [NASA ADS] [CrossRef] [Google Scholar]
- Bianchi, S., Guainazzi, M., Matt, G., & Fonseca Bonilla, N. 2007, A&A, 467, L19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bianchi, S., Corral, A., Panessa, F., et al. 2008a, MNRAS, 385, 195 [NASA ADS] [CrossRef] [Google Scholar]
- Bianchi, S., La Franca, F., Matt, G., et al. 2008b, MNRAS, 389, L52 [NASA ADS] [Google Scholar]
- Bianchi, S., Guainazzi, M., Matt, G., Fonseca Bonilla, N., & Ponti, G. 2009a, A&A, 495, 421 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bianchi, S., Piconcelli, E., Chiaberge, M., et al. 2009b, ApJ, 695, 781 [NASA ADS] [CrossRef] [Google Scholar]
- de la Calle Pérez, I., Longinotti, A. L., Guainazzi, M., et al. 2010, A&A, in press [Google Scholar]
- Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Dovčiak, M., Karas, V., & Yaqoob, T. 2004, ApJS, 153, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Elitzur, M., & Shlosman, I. 2006, ApJ, 648, L101 [NASA ADS] [CrossRef] [Google Scholar]
- Elvis, M., Risaliti, G., Nicastro, F., et al. 2004, ApJ, 615, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Fabian, A. C., Iwasawa, K., Reynolds, C. S., & Young, A. J. 2000, PASP, 112, 1145 [NASA ADS] [CrossRef] [Google Scholar]
- Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761 [Google Scholar]
- Grandi, P., Haardt, F., Ghisellini, G., et al. 1998, ApJ, 498, 220 [NASA ADS] [CrossRef] [Google Scholar]
- Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lobban, A. P., Reeves, J. N., Porquet, D., et al. 2010, MNRAS, 408, 551 [NASA ADS] [CrossRef] [Google Scholar]
- Longinotti, A. L., de La Calle, I., Bianchi, S., Guainazzi, M., & Dovčiak, M. 2008, Mem. Soc. Astron. Ital., 79, 259 [Google Scholar]
- Marconi, A., Risaliti, G., Gilli, R., et al. 2004, MNRAS, 351, 169 [NASA ADS] [CrossRef] [Google Scholar]
- Matt, G. 2002, MNRAS, 337, 147 [NASA ADS] [CrossRef] [Google Scholar]
- Matt, G., Bianchi, S., de Rosa, A., Grandi, P., & Perola, G. C. 2006, A&A, 445, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nandra, K., O’Neill, P. M., George, I. M., & Reeves, J. N. 2007, MNRAS, 382, 194 [NASA ADS] [CrossRef] [Google Scholar]
- Nicastro, F. 2000, ApJ, 530, L65 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Nicastro, F., Fiore, F., Perola, G. C., & Elvis, M. 1999, ApJ, 512, 184 [NASA ADS] [CrossRef] [Google Scholar]
- Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of gaseous nebulae and active galactic nuclei, ed. G. J. Osterbrock, & D. E. Ferland [Google Scholar]
- Osterbrock, D. E., & Shuder, J. M. 1982, ApJS, 49, 149 [NASA ADS] [CrossRef] [Google Scholar]
- Panessa, F., Carrera, F. J., Bianchi, S., et al. 2009, MNRAS, 398, 1951 [NASA ADS] [CrossRef] [Google Scholar]
- Perola, G. C., Matt, G., Fiore, F., et al. 2000, A&A, 358, 117 [NASA ADS] [Google Scholar]
- Perola, G. C., Matt, G., Cappi, M., et al. 2002, A&A, 389, 802 [Google Scholar]
- Piconcelli, E., Jimenez-Bailón, E., Guainazzi, M., et al. 2005, A&A, 432, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Puccetti, S., Fiore, F., Risaliti, G., et al. 2007, MNRAS, 377, 607 [NASA ADS] [CrossRef] [Google Scholar]
- Risaliti, G., Elvis, M., Fabbiano, G., Baldi, A., & Zezas, A. 2005, ApJ, 623, L93 [NASA ADS] [CrossRef] [Google Scholar]
- Winter, L. M., Lewis, K. T., Koss, M., et al. 2010, ApJ, 710, 503 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Best-fit model parameters for the emission features in the Suzaku XIS spectra (3–10 keV).
Best-fits for the Suzaku XIS spectra (0.6–10 keV, first two columns) and the Suzaku+Swift BAT spectra (0.6–150 keV, last two columns).
All Figures
![]() |
Fig. 1 Soft (0.3–2 keV) and hard (2–10 keV) light-curves of the Suzaku 100 ks observation of MCG+8-11-11 (only the XIS0 is plotted for clarity). The bottom panel shows the ratio between the two curves. |
In the text |
![]() |
Fig. 2 Left: residuals in the (rest-frame) iron line band, after a 3–10 keV fit on the three Suzaku XIS spectra, with a simple power law. Right: the same as above, but for model III, after the relativistic component of the iron line is removed (see Table 1). Data have been rebinned for clarity. |
In the text |
![]() |
Fig. 3 Best-fit model and Δχ2 for the three XIS spectra, with the same model as in the 3–10 keV fit, but including Galactic absorption, a warm absorber and an unidentified absorption edge (see text for details). Note the residuals above 8 keV. |
In the text |
![]() |
Fig. 4 Left: the best-fit model for the 0.6–10 keV band, when extrapolated to the band covered by the PIN spectrum. Right: the Suzaku+Swift BAT best-fit model in the whole 0.6–150 keV band. |
In the text |
![]() |
Fig. 5 The 15–50 keV Swift BAT 54-month light curve of MCG+8-11-11 (a bin is 15 days). The countrate was converted to a 15–50 keV flux, assuming no spectral variation. The Suzaku PIN flux is shown in comparison, with its uncertainty, with the dotted lines and the red symbol. The agreement at the date of the Suzaku observation is good. |
In the text |
![]() |
Fig. 6 Suzaku+Swift BAT best fit contour plots (black, red, green lines correspond to 68, 90, and 99% confidence level for two interesting parameters). Left: Compton reflection amount R vs. cutoff energy Ec; middle: Compton reflection amount R vs. iron abundance AFe; right: Compton reflection amount R vs. spectral photon index Γ. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.