Abstract

We methodically model the broad-band Suzaku spectra of a small sample of six ‘bare’ Seyfert galaxies: Ark 120, Fairall 9, MCG −02-14-009, Mrk 335, NGC 7469 and SWIFT J2127.4+5654. The analysis of bare Seyferts allows a consistent and physical modelling of active galactic nuclei due to a weak amount of any intrinsic warm absorption, removing the degeneracy between the spectral curvature due to warm absorption and the red wing of the Fe K region. Through effective modelling of the broad-band spectrum and investigating the presence of narrow neutral or ionized emission lines and reflection from distant material, we obtain an accurate and detailed description of the Fe K line region using models such as laor, kerrdisk and kerrconv.

Results suggest that ionized emission lines at 6.7 and 6.97 keV (particularly Fe xxvi) are relatively common, and the inclusion of these lines can greatly affect the parameters obtained with relativistic models, i.e. spin, emissivity, inner radius of emission and inclination. Moderately broad components are found in all objects, but typically the emission originates from tens of rg, rather than within <6 rg of the black hole. Results obtained with kerrdisk line profiles suggest an average emissivity of q ∼ 2.3 at intermediate spin values with all objects ruling out the presence of a maximally spinning black hole at the 90 per cent confidence level. We also present new spin constraints for Mrk 335 and NGC 7469 with intermediate values of a= 0.70+0.12−0.01 and 0.69+0.09−0.09, respectively.

1 INTRODUCTION

Analysis of the Fe K line profile can help us determine some of the intrinsic properties of the central black hole (BH) in active galactic nuclei (AGN) and its accretion disc (e.g. Fabian et al. 1989; Laor 1991). Foremost in the current climate is the determination of BH spin; for example, recent observations by Miniutti et al. (2007, 2009) and Schmoll et al. (2009) have analysed the broad Fe Kα region in order to constrain the spin of the central BH. A BH can be characterized simply by its mass and spin. Many objects have now been classified in terms of their mass into three categories: Galactic BH, intermediate-mass BH and supermassive BH (SMBH). The spin of the BH determines the nature of the space–time metric in the regions close to it. The spin parameter a=cJ/GM2 (where J = angular momentum and 0 < a < 0.998) is used to describe the spin of the BH, where a is limited to a maximum value of a = 0.998 at the Thorne limit. This is due to photon capture in which photons travelling on ‘negative’ angular momentum orbits are preferentially captured by the BH, therefore producing an upper bound to the spin parameter a and hence limiting the innermost stable circular orbit (ISCO) to a minimum of rISCO = 1.235 rg (Thorne 1974). Constraining the spin of SMBHs in AGN and studying the distribution of BH spin can aid our understanding of the evolution of AGN and the BHs themselves, e.g. mergers, relativistic jets and variability (Blandford & Znajek 1977; Volonteri, Sikora & Lasota 2007; King, Pringle & Hofmann 2008).

Line emission from the inner regions of the accretion disc can become broadened due to relativistic effects and Doppler motions (Fabian et al. 1989) resulting in an asymmetric profile. Evidence for such broadening was found using the X-ray CCD detectors onboard ASCA (Tanaka, Inoue & Holt 1994) by Mushotzky et al. (1995), Tanaka et al. (1995) and Nandra et al. (1997), typically over the 0.5–10.0 keV range. Now with X-ray spectra of increasing quality and over wider energy ranges such as that obtained with the XIS and high-energy HXD detectors onboard Suzaku (0.5–60.0 keV; Mitsuda et al. 2007), the Fe K region of AGN can be examined in detail. High-energy X-ray data are important since they allow the reflection component and its strength to be properly fitted, assessing its contribution to the continuum and Fe K region (e.g. Reeves et al. 2007). Fitting features such as the Compton hump at ∼30 keV allows, for example, the ionization state of the reflecting material to be determined (Ross & Fabian 2005). With the aim of measuring the properties of the accretion disc and the central BH itself, broad-band data allow us to start making constraints on parameters in these regions based upon the shape of the Fe K line profile.

Previous studies of iron lines have been made using data from XMM–Newton over the 2.5–10.0 keV (Nandra et al. 2007) and 0.6–10.0 keV energy ranges (Brenneman & Reynolds 2009), finding complex emission in the Fe K band in the majority of Type 1 Seyfert AGN over and above narrow-line components originating from distant material. In a sample of 26 objects, Nandra et al. (2007) found that narrow 6.4-keV emission is ubiquitous amongst AGN, and the broad Fe K lines feature in ∼40–50 per cent of AGN and ionized emission due to Fe xxvi and Fe xxv is relatively rare amongst AGN. Brenneman & Reynolds (2009) found that 4/8 AGN were best fitted by a model consisting of relativistically blurred reflection from the inner regions of the accretion disc with 2/8 objects suggesting non-zero spin; however they note that the 10.0 keV cut-off with the EPIC-pn camera limits their findings of the reflection continuum.

This small sample features predominantly bare Seyferts, i.e. with a very weak warm absorber. This is in an attempt to simplify the modelling of the broad-band spectrum and therefore provide a basic understanding of the properties of the Fe K region and any accompanying broad red wing without the need for debate over differing interpretations of the origin of various absorption components in the spectrum (Turner & Miller 2009). Using data from Suzaku’s XIS (Koyama et al. 2007) and HXD (Takahashi et al. 2007) detectors spanning 0.5–60.0 keV and Burst Alert Telescope (BAT) data (from the Swift 22-month all-sky survey; Tueller et al. 2010) over 20.0–100.0 keV provides the broad-band spectra necessary for detailed modelling and measurement of the Fe K region and the associated Compton reflection hump. This approach ensures that a robust model of the Fe K line region (and the broad-band spectrum) can be made with the aim of establishing the degree to which narrow ionized emission lines and relativistically broadened components are required.

2 OBSERVATIONS AND DATA REDUCTION

2.1 Observations

The objects in this sample are described in Table 1 and are all Seyfert 1, radio-quiet AGN with low intrinsic absorption available in the public Suzaku data archive with spectra from the XIS and HXD detectors. All objects are detected above 15 keV and are relatively nearby with redshift z < 0.05. Hard X-ray data from the BAT instrument onboard Swift are used in addition to the HXD data in all observations other than MCG −02-14-009, for which no BAT spectrum was publically available (it was, however, detected by BAT in the 39-month survey with a 14–150 keV flux of 1.5 × 10−11 erg cm−2 s−1; Cusumano et al. 2010). In addition to the Suzaku data, observations with XMM–Newton are also used for comparison, except for SWIFT J2127.4+5654 which was not publically available. The details of the observations used in this analysis are given in Table 2.

Table 1

The Suzaku Seyfert sample.

ObjectRA (J2000)Dec. (J2000)RedshiftNH (Gal) (1022 cm−2)
Ark 12005 16 11.4−00 09 000.0330.0978
Fairall 901 23 45.8−58 48 210.0470.0316
MCG −02-14-00905 16 21.2−10 33 410.0280.0924
Mrk 33500 06 19.5+20 12 100.0260.0356
NGC 746923 00 44.4+08 36 170.0160.0445
SWIFT J2127.4+565421 27 45.0+56 56 400.0140.7650
ObjectRA (J2000)Dec. (J2000)RedshiftNH (Gal) (1022 cm−2)
Ark 12005 16 11.4−00 09 000.0330.0978
Fairall 901 23 45.8−58 48 210.0470.0316
MCG −02-14-00905 16 21.2−10 33 410.0280.0924
Mrk 33500 06 19.5+20 12 100.0260.0356
NGC 746923 00 44.4+08 36 170.0160.0445
SWIFT J2127.4+565421 27 45.0+56 56 400.0140.7650
Table 1

The Suzaku Seyfert sample.

ObjectRA (J2000)Dec. (J2000)RedshiftNH (Gal) (1022 cm−2)
Ark 12005 16 11.4−00 09 000.0330.0978
Fairall 901 23 45.8−58 48 210.0470.0316
MCG −02-14-00905 16 21.2−10 33 410.0280.0924
Mrk 33500 06 19.5+20 12 100.0260.0356
NGC 746923 00 44.4+08 36 170.0160.0445
SWIFT J2127.4+565421 27 45.0+56 56 400.0140.7650
ObjectRA (J2000)Dec. (J2000)RedshiftNH (Gal) (1022 cm−2)
Ark 12005 16 11.4−00 09 000.0330.0978
Fairall 901 23 45.8−58 48 210.0470.0316
MCG −02-14-00905 16 21.2−10 33 410.0280.0924
Mrk 33500 06 19.5+20 12 100.0260.0356
NGC 746923 00 44.4+08 36 170.0160.0445
SWIFT J2127.4+565421 27 45.0+56 56 400.0140.7650
Table 2

Summary of observations for the objects in the sample.1 The observed 2–10 keV flux for XIS and EPIC-pn instruments, 15–50 keV flux for HXD and 20–100 keV flux for BAT, in units of 10−11 erg cm−2 s−1 from model A.

ObjectMissionInstrumentDateExposure (s)Count rateFlux1Obs. ID
Ark 120XIS100 8641.896 ± 0.0033.06
Suzaku2007/04/01702014010
HXD89 4700.114 ± 0.0033.46
XMMEPIC-pn2003/08/2478 17025.04 ± 0.0183.87014719010
SwiftBAT2453 000(6.9 ± 0.4) × 10−44.88
Fairall 9XIS167 8141.718 ± 0.0022.53
Suzaku2007/06/07702043010
HXD127 3100.089 ± 0.0022.97
XMMEPIC-pn2000/07/0525 8305.825 ± 0.0151.170101040201
SwiftBAT3280 000(4.9 ± 0.4) × 10−44.35
MCG-02-14-009XIS142 1520.216 ± 0.0010.42
Suzaku2008/08/28703060010
HXD120 0280.017 ± 0.0020.61
XMMEPIC-pn2009/02/2782 7901.946 ± 0.0050.430550640101
Mrk 335XIS151 2961.324 ± 0.0021.49
Suzaku2006/08/21701031010
HXD131 7440.012 ± 0.0011.31
XMMEPIC-pn2006/01/0380 36016.43 ± 0.0141.820306870101
SwiftBAT3273 000(2.5 ± 0.3) × 10−41.74
NGC 7469XIS112 1131.091 ± 0.0022.11
Suzaku2008/06/24703028010
HXD85 3150.068 ± 0.0023.22
EPIC-pn2004/11/3053 01424.09 ± 0.0212.740207090101
XMM
EPIC-pn2004/12/0354 97716.32 ± 0.0172.890207090201
SwiftBAT3286 000(6.6 ± 0.3) × 10−44.82
SWIFT J2127.4+5654XIS91 7301.373 ± 0.0043.77
Suzaku2007/12/09702122010
HXD83 3210.074 ± 0.0023.33
SwiftBAT3999 000(4.2 ± 0.3) × 10−44.16
ObjectMissionInstrumentDateExposure (s)Count rateFlux1Obs. ID
Ark 120XIS100 8641.896 ± 0.0033.06
Suzaku2007/04/01702014010
HXD89 4700.114 ± 0.0033.46
XMMEPIC-pn2003/08/2478 17025.04 ± 0.0183.87014719010
SwiftBAT2453 000(6.9 ± 0.4) × 10−44.88
Fairall 9XIS167 8141.718 ± 0.0022.53
Suzaku2007/06/07702043010
HXD127 3100.089 ± 0.0022.97
XMMEPIC-pn2000/07/0525 8305.825 ± 0.0151.170101040201
SwiftBAT3280 000(4.9 ± 0.4) × 10−44.35
MCG-02-14-009XIS142 1520.216 ± 0.0010.42
Suzaku2008/08/28703060010
HXD120 0280.017 ± 0.0020.61
XMMEPIC-pn2009/02/2782 7901.946 ± 0.0050.430550640101
Mrk 335XIS151 2961.324 ± 0.0021.49
Suzaku2006/08/21701031010
HXD131 7440.012 ± 0.0011.31
XMMEPIC-pn2006/01/0380 36016.43 ± 0.0141.820306870101
SwiftBAT3273 000(2.5 ± 0.3) × 10−41.74
NGC 7469XIS112 1131.091 ± 0.0022.11
Suzaku2008/06/24703028010
HXD85 3150.068 ± 0.0023.22
EPIC-pn2004/11/3053 01424.09 ± 0.0212.740207090101
XMM
EPIC-pn2004/12/0354 97716.32 ± 0.0172.890207090201
SwiftBAT3286 000(6.6 ± 0.3) × 10−44.82
SWIFT J2127.4+5654XIS91 7301.373 ± 0.0043.77
Suzaku2007/12/09702122010
HXD83 3210.074 ± 0.0023.33
SwiftBAT3999 000(4.2 ± 0.3) × 10−44.16
Table 2

Summary of observations for the objects in the sample.1 The observed 2–10 keV flux for XIS and EPIC-pn instruments, 15–50 keV flux for HXD and 20–100 keV flux for BAT, in units of 10−11 erg cm−2 s−1 from model A.

ObjectMissionInstrumentDateExposure (s)Count rateFlux1Obs. ID
Ark 120XIS100 8641.896 ± 0.0033.06
Suzaku2007/04/01702014010
HXD89 4700.114 ± 0.0033.46
XMMEPIC-pn2003/08/2478 17025.04 ± 0.0183.87014719010
SwiftBAT2453 000(6.9 ± 0.4) × 10−44.88
Fairall 9XIS167 8141.718 ± 0.0022.53
Suzaku2007/06/07702043010
HXD127 3100.089 ± 0.0022.97
XMMEPIC-pn2000/07/0525 8305.825 ± 0.0151.170101040201
SwiftBAT3280 000(4.9 ± 0.4) × 10−44.35
MCG-02-14-009XIS142 1520.216 ± 0.0010.42
Suzaku2008/08/28703060010
HXD120 0280.017 ± 0.0020.61
XMMEPIC-pn2009/02/2782 7901.946 ± 0.0050.430550640101
Mrk 335XIS151 2961.324 ± 0.0021.49
Suzaku2006/08/21701031010
HXD131 7440.012 ± 0.0011.31
XMMEPIC-pn2006/01/0380 36016.43 ± 0.0141.820306870101
SwiftBAT3273 000(2.5 ± 0.3) × 10−41.74
NGC 7469XIS112 1131.091 ± 0.0022.11
Suzaku2008/06/24703028010
HXD85 3150.068 ± 0.0023.22
EPIC-pn2004/11/3053 01424.09 ± 0.0212.740207090101
XMM
EPIC-pn2004/12/0354 97716.32 ± 0.0172.890207090201
SwiftBAT3286 000(6.6 ± 0.3) × 10−44.82
SWIFT J2127.4+5654XIS91 7301.373 ± 0.0043.77
Suzaku2007/12/09702122010
HXD83 3210.074 ± 0.0023.33
SwiftBAT3999 000(4.2 ± 0.3) × 10−44.16
ObjectMissionInstrumentDateExposure (s)Count rateFlux1Obs. ID
Ark 120XIS100 8641.896 ± 0.0033.06
Suzaku2007/04/01702014010
HXD89 4700.114 ± 0.0033.46
XMMEPIC-pn2003/08/2478 17025.04 ± 0.0183.87014719010
SwiftBAT2453 000(6.9 ± 0.4) × 10−44.88
Fairall 9XIS167 8141.718 ± 0.0022.53
Suzaku2007/06/07702043010
HXD127 3100.089 ± 0.0022.97
XMMEPIC-pn2000/07/0525 8305.825 ± 0.0151.170101040201
SwiftBAT3280 000(4.9 ± 0.4) × 10−44.35
MCG-02-14-009XIS142 1520.216 ± 0.0010.42
Suzaku2008/08/28703060010
HXD120 0280.017 ± 0.0020.61
XMMEPIC-pn2009/02/2782 7901.946 ± 0.0050.430550640101
Mrk 335XIS151 2961.324 ± 0.0021.49
Suzaku2006/08/21701031010
HXD131 7440.012 ± 0.0011.31
XMMEPIC-pn2006/01/0380 36016.43 ± 0.0141.820306870101
SwiftBAT3273 000(2.5 ± 0.3) × 10−41.74
NGC 7469XIS112 1131.091 ± 0.0022.11
Suzaku2008/06/24703028010
HXD85 3150.068 ± 0.0023.22
EPIC-pn2004/11/3053 01424.09 ± 0.0212.740207090101
XMM
EPIC-pn2004/12/0354 97716.32 ± 0.0172.890207090201
SwiftBAT3286 000(6.6 ± 0.3) × 10−44.82
SWIFT J2127.4+5654XIS91 7301.373 ± 0.0043.77
Suzaku2007/12/09702122010
HXD83 3210.074 ± 0.0023.33
SwiftBAT3999 000(4.2 ± 0.3) × 10−44.16

The objects featuring in this sample, however, are all objects with a very low degree of warm absorption. Three objects in this analysis have been noted as having some additional absorption in previous observations: MCG −02-14-009, Mrk 335 and NGC 7469. Gallo et al. (2006) found a small neutral absorber (NH < 1021 cm−2) and an O vii edge in a re-analysis of a short (∼5 ks net exposure) observation of MCG −02-14-009 with XMM–Newton. However, including these components in the Suzaku data used here makes no improvement (NH < 3 × 1020 cm−2) and the optical depth of the O vii edge can only be constrained to τ < 0.04. Therefore, MCG −02-14-009 can be considered as ‘bare’ for the purposes of this paper.

Mrk 335 has also shown evidence for a warm emitter whilst being in a low state when observed with XMM–Newton (Grupe et al. 2008); however in the Suzaku data used here, Mrk 335 has a 0.5–2.0 keV flux 14 × higher than in the low state observation and as such any warm emission features are entirely dominated by the continuum. As a result of this, Mrk 335 is also suitable for inclusion within this sample. Previous observations of NGC 7469 have noted some degree of X-ray and UV absorption (NH∼ 1020–1021 cm−2; Scott et al. 2005; Blustin et al. 2007). The Suzaku data used here are consistent with the previous work using data from XMM–Newton, requiring an O vii edge depth of τ < 0.1; however, there is no effect upon the Fe K parameters and NGC 7469 has also been included within this sample. No significant additional absorption features below 2 keV can be seen in the Suzaku data.

2.2 Data reduction

All the Suzaku data in this paper were reduced using the heasoft reduction and analysis package (version 6.8). XIS source spectra were extracted from circular regions of 3.0 arcmin within xselect centred upon the source at the on-axis pointing position. Similarly, background spectra were extracted from 3.0 arcmin circular regions, taking care not to include the source or the Fe 55 calibration sources in the corners of the CCD’s field of view. Only data from the front-illuminated XIS cameras were used, i.e. the XIS 0 and XIS 3 cameras, due to their greater sensitivity at Fe K energies; however, the observation of Mrk 335 (Obs. ID 701031010) also includes the now non-operational front-illuminated XIS 2. It should be noted that only data from the XIS 3 camera were available for the Suzaku observation of SWIFT J2127.4+5654 (Obs. ID 702122010) since the XIS 0 camera was not operating properly.

XIS redistribution matrix files were created using the heasoft tool xisrmfgen and the ancillary response files using xissimarfgen. Data from the 3×3 and 5×5 modes were grouped together and the data from each of the XIS front-illuminated CCDs were co-added in order to increase the signal-to-noise ratio (S/N) using mathpha, addrmf and addarf. We ignore all XIS data below 0.5 keV, above 10.0 keV and between 1.7 and 1.95 keV due to uncertainties in the calibration of the detectors around the Si K edge.

The HXD/PIN spectrum was extracted from the cleaned HXD/PIN event files and corrected for dead time using the tool hxddtcor. The tuned HXD/PIN background events were used for background subtraction (Fukazawa et al. 2009) using identical good time intervals as per the source events, generated with 10 × the actual background count rate to minimize photon noise. A simulated cosmic X-ray background was also produced using xspec v 12.5.1n with a spectral form identical to Gruber et al. (1999) and then added to the corrected non-X-ray background file to create a single background file. In the analysis of the HXD/PIN spectra, we typically consider data between 15.0 and 60.0 keV; however, in some cases such as Mrk 335, the data were also ignored below 20.0 keV due to thermal noise. Hard X-ray data also obtained from the Swift 22-month BAT catalogue were included for all objects (other than MCG −02-14-009) over the energy range of 20.0–100.0 keV.

3 ANALYSIS AND RESULTS

Spectral analysis and model fitting are performed with xspec v 12.5.1n (Arnaud 1996); all models were modified by Galactic absorption via the wabs multiplicative model (Morrison & McCammon 1983) using a Galactic column density obtained using the nh ftool appropriate for each source giving the weighted average NH value of the (LAB) Survey of Galactic H i (Kalberla et al. 2005), using abundances from Anders & Grevesse (1989). In all fits, a constant factor was introduced to account for the cross-normalization between the XIS, HXD/PIN and BAT detectors, fixed at 1.16 or 1.18 between XIS and HXD/PIN according to the nominal pointing position1 and allowed to vary between XIS and BAT detectors (typically ∼1.0, indicating little variability between the Suzaku and BAT data sets). Data are fitted over the full 0.5–100.0 keV range, excluding those regions mentioned above. The χ2 minimization technique is used throughout, quoting 90 per cent errors for one interesting parameter (▵χ2 = 2.71) unless otherwise stated. Where the significance of components is stated according to ▵χ2, the component has been subtracted from the final model and then refitted to ensure that the order in which components are added does not affect the quoted statistical significance. Table 3 summarizes the subsequent models and the components used to model the continuum, soft excess and Fe K region.

Table 3

Summary of models and components used for each feature.

ModelContinuumDistant reflectionSoft excess6.4 keV Core6.4 keV broad
APowerlawpexravcompttGaussianGaussian
BPowerlawreflionxcompttWithin reflionxNone
CPowerlawreflionxcompttWithin reflionxlaor
DPowerlawreflionxcompttWithin reflionxkerrdisk
EPowerlawreflionxBlurred reflionxWithin reflionxBlurred reflionx
FPowerlawcompttGaussianBlurred reflionx
ModelContinuumDistant reflectionSoft excess6.4 keV Core6.4 keV broad
APowerlawpexravcompttGaussianGaussian
BPowerlawreflionxcompttWithin reflionxNone
CPowerlawreflionxcompttWithin reflionxlaor
DPowerlawreflionxcompttWithin reflionxkerrdisk
EPowerlawreflionxBlurred reflionxWithin reflionxBlurred reflionx
FPowerlawcompttGaussianBlurred reflionx
Table 3

Summary of models and components used for each feature.

ModelContinuumDistant reflectionSoft excess6.4 keV Core6.4 keV broad
APowerlawpexravcompttGaussianGaussian
BPowerlawreflionxcompttWithin reflionxNone
CPowerlawreflionxcompttWithin reflionxlaor
DPowerlawreflionxcompttWithin reflionxkerrdisk
EPowerlawreflionxBlurred reflionxWithin reflionxBlurred reflionx
FPowerlawcompttGaussianBlurred reflionx
ModelContinuumDistant reflectionSoft excess6.4 keV Core6.4 keV broad
APowerlawpexravcompttGaussianGaussian
BPowerlawreflionxcompttWithin reflionxNone
CPowerlawreflionxcompttWithin reflionxlaor
DPowerlawreflionxcompttWithin reflionxkerrdisk
EPowerlawreflionxBlurred reflionxWithin reflionxBlurred reflionx
FPowerlawcompttGaussianBlurred reflionx

3.1 Model A – Parametrization of spectra

Model A is intended to provide a simple parametrization of the spectra and give an insight into the presence of basic components of the spectra and the extent to which a possibly relativistically broadened component is required to model the Fe K line region. None of the spectra required any significant warm absorber (see Section 2.1), thereby simplifying any present broadening in the region. Additionally, the comptt model representing Comptonization of soft photons in a hot plasma above the disc (Titarchuk 1994) with a soft photon input temperature of 0.02 keV is employed to account for the soft excess if present in the spectra (see Fig. 1). Porquet et al. (2004) also found that in a sample of PG quasars the soft excess was better modelled in this way, rather than by thermal emission from the accretion disc. A second soft powerlaw component instead of comptt also gives a similar parametrization of the soft excess. The narrow 6.4 keV core due to reflection from distant material is present in all six objects and has been modelled with a narrow Gaussian with width σ free to vary. The narrow component is not resolved in any of the spectra and as such the width is fixed at σ = 0.01 keV. Emission resulting from Fe Kβ is also accounted for with the line energy fixed at 7.056 keV, width fixed to that of the narrow Kα and flux tied to 13 per cent of the Kα component.

The 0.5–70.0 keV spectrum and residuals after modelling of the continuum with a simple powerlaw and wabs to account for Galactic absorption; the entire Fe K region and any soft excess is left unmodelled. XIS data are in black, HXD in red and BAT data are represented by blue circles.
Figure 1

The 0.5–70.0 keV spectrum and residuals after modelling of the continuum with a simple powerlaw and wabs to account for Galactic absorption; the entire Fe K region and any soft excess is left unmodelled. XIS data are in black, HXD in red and BAT data are represented by blue circles.

Consistent with this, neutral distant reflection is accounted for using the pexrav model (Magdziarz & Zdziarski 1995) applied to the broad-band 0.5–100 keV spectra. This model requires the input of a photon index Γ which is tied to the continuum power law, the normalization of the pexrav component is also tied to that of the power law, abundances are assumed to be solar (Anders & Grevesse 1989) and the disc inclination to the observer is fixed at cos i = 0.87 throughout. The reflection fraction R=Ω/2π is left as a free parameter (where R = 1 denotes reflection from material subtending 2π sr). The cut-off energy for the pexrav component is fixed at 1000 keV, except for NGC 7469 and SWIFT J2127.4+5654 which show indications of a roll-over at high energies, occurring at Ec = 119+65−31 and 49+43−14 keV, respectively. Such a roll-over was found in SWIFT J2127.4+5654 previously by Malizia et al. (2008) at Ec = 33+19−10 keV.

Examining the residuals after the application of these components, some objects have more complex features such as excess emission around energies of 6.7 and 6.97 keV relating to narrow ionized emission from Fe xxv and Fe xxvi, respectively, again due to distant photoionized gas (see Figs 2 and 3). Similarly to the modelling of the narrow 6.4-keV core, the energies of these lines are free to vary as is the normalization; however, the width remains fixed to σ = 0.01 keV. A broad Gaussian is also added to the model to account for a red wing in the spectra with energy, normalization and σBroad as free parameters. In this case, the energy is allowed to drop below 6.4 keV in the rest frame as this model is only intended as a simple parametrization of the Fe K region and a test of the significance of the components used to model the region.

The 0.5–70.0 keV spectrum and residuals after modelling of the continuum with a powerlaw, neutral reflection from the pexrav model, comptt to model the soft excess and wabs to account for Galactic absorption; the entire Fe K region is left unmodelled. XIS data are in black, HXD in red and BAT data are represented by blue circles.
Figure 2

The 0.5–70.0 keV spectrum and residuals after modelling of the continuum with a powerlaw, neutral reflection from the pexrav model, comptt to model the soft excess and wabs to account for Galactic absorption; the entire Fe K region is left unmodelled. XIS data are in black, HXD in red and BAT data are represented by blue circles.

The 4–9 keV residuals after modelling of the continuum with a powerlaw, comptt to model the soft excess and Galactic photoelectric absorption. The dashed vertical line represents 6.4 keV in the rest frame.
Figure 3

The 4–9 keV residuals after modelling of the continuum with a powerlaw, comptt to model the soft excess and Galactic photoelectric absorption. The dashed vertical line represents 6.4 keV in the rest frame.

A good fit is obtained by model A in all objects, particularly Fairall 9, NGC 7469 and SWIFT J2127.4+5654. No significant soft excess is found in MCG −02-14-009 and SWIFT J2127.4+5654; the latter (if present) is likely to be mostly absorbed due to the relatively high amount of Galactic absorption (NH≈ 7.65 × 1021 cm−2). In SWIFT J2127.4+5654, there is some indication of the presence of an intrinsic neutral absorber, albeit of relatively low column density (as found in Miniutti et al. 2009). Instead, modelling the spectrum with an absorbed powerlaw using the zphabs model at the redshift of the source improves the fit by ▵χ2∼ 59 for one additional free parameter with intrinsic column density NH≈ 8.2+0.6−0.6× 1020 cm−2 and photon index Γ = 2.19+0.05−0.05.

Statistically significant narrow ionized emission is found in most objects, with SWIFT J2127.4+5654 and MCG −02-14-009 having very little requirement for these features. These lines occur at energies likely originating from Fe xxv and Fe xxvi with emission at ∼6.97 keV being particularly evident in the spectra, proving strong in Fairall 9 for which the fit improves by ▵χ2∼ 24 with the introduction of a narrow Gaussian (and therefore two additional free parameters) at 6.98+0.02−0.02 keV and to a lesser extent in Ark 120 (▵χ2∼ 12 at 6.96+0.04−0.04 keV). The only objects not showing residuals at ∼6.7 or ∼6.97 keV are NGC 7469 and SWIFT J2127.4+5654, also in agreement with Miniutti et al. (2009). Residuals at ∼6.7 keV are only found to be significant in two objects: Fairall 9 and Mrk 335 with improvements of ▵χ2∼ 7 and ∼14, respectively.

With the aim of determining the extent to which emission from further in to the BH is required, a broad component significantly improves the quality of the fit for most objects, typically ▵χ2 > 20 for three additional free parameters and with linewidths of the order σBroad≳ 0.3 keV. Only in MCG −02-14-009 does the addition of a broad Gaussian not make a particularly significant improvement (▵χ2∼ 3). This is surprising given evidence to the contrary by Porquet (2006) in which a prominent broad and statistically significant iron line was found in a short XMM–Newton observation with an equivalent width of EW ∼ 527+277−248 eV whereas only EW ∼ 92+59−56 eV is found here. However, the XMM–Newton observation was only 5-ks net exposure and as a result the Fe K line parameters are poorly constrained, while the Suzaku data also allow the broad-band continuum to be better constrained (however, we cannot rule out some variability between these two observations). In SWIFT J2127.4+5654, this may be due to the relatively high best-fitting value of the reflection component with R = 2.5+1.1−0.8 which may reduce the significance of a broad component. Mrk 335 features a relatively broad Gaussian with σBroad = 0.50+0.13−0.11 keV and EW = 134+42−38 eV; however, this feature is not as strong as that found by Larsson et al. (2008) with σ = 0.45+0.10−0.06 keV and EW = 250+40−39 eV. Larsson et al. (2008) also found that the Fe xxvi emission line does not improve the fit; however, an improvement of ▵χ2∼ 8 for three additional free parameters is found here with the introduction of a line at 6.98+0.06−0.14 keV.

In general, the broad emission found in these objects is typically too strong to be modelled as purely a Compton shoulder despite the similarity of the line energy of the broad Gaussian with the first-order Compton shoulder peak energy of ∼6.24 keV (see Table 4). However, we cannot rule out the possibility of a contribution from a Compton shoulder. Given the relatively high EW of the broad-line component (mean EW∼ 86 eV) compared to the narrow core, the majority of the emission is likely to arise from a broadened component. Indeed, the Compton shoulder is unlikely to contribute more than ∼20 per cent of the 6.4 keV core flux (e.g. George & Fabian 1991; Matt 2002). Note that the Compton shoulder is included in subsequent models through the use of a reflionx reflection component.

Table 4

Model A components for Suzaku XIS, HXD and BAT data from Swift.

Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ1.99+0.04−0.041.87+0.04−0.041.89+0.04−0.042.00+0.04−0.021.80+0.01−0.012.19+0.05−0.05
Norma1.04+0.05−0.060.68+0.04−0.050.13+0.04−0.040.56+0.16−0.020.57+0.06−0.071.60+0.07−0.06
Plasma temperature (keV)<8.5<18.6<8.3<8.4
Plasma optical depth τ<2.4<2.1<2.20.8+0.2−0.5
Fluxb1.05+0.05−0.050.59+0.03−0.030.98+0.02−0.020.22+0.03−0.03
▵χ2250240955194
Neutral Fe Kα line (keV)6.40+0.02−0.026.39+0.01−0.016.43+0.06−0.086.39+0.03−0.026.40+0.02−0.026.37+0.05−0.05
σNarrow (eV)<84<45<110<72<43<190
EWNarrow (eV)40+11−1290+9−1055+28−5160+38−2371+38−25<35
Flux (10−5 ph cm−2 s−1)1.49+0.41−0.462.50+0.26−0.270.29+0.15−0.271.00+0.63−0.381.92+1.02−0.69<1.49
Fe xxv line (keV)6.74+0.04−0.046.69+0.06−0.05
EW (eV)17+8−933+13−7
Flux (10−5 ph cm−2 s−1)0.46+0.22−0.240.47+0.18−0.10
▵χ2714
Fe xxvi line (keV)6.96+0.04−0.046.98+0.02−0.026.94+0.06−0.096.98+0.06−0.14
EW (eV)24+11−1130+8−10<5121+13−13
Flux (10−5 ph cm−2 s−1)0.70+0.32−0.320.70+0.19−0.23<0.220.24+0.07−0.06
▵χ2122438
Broad line (keV)6.36+0.08−0.096.16+0.20−0.206.33+0.21−0.246.27+0.13−0.176.32+0.06−0.116.32+0.24−0.25
σBroad (keV)0.32+0.11−0.090.37+0.21−0.130.28+0.28−0.130.50+0.13−0.110.15+0.07−0.030.39+0.45−0.26
EWBroad (eV)105+26−2453+27−3892+59−56134+42−3862+28−4471+53−48
Flux (10−5 ph cm−2 s−1)3.74+0.93−0.861.64+0.83−1.170.49+0.31−0.302.28+0.72−0.651.55+0.70−1.103.43+2.56−2.32
▵χ257283532412
Rfrac0.88+0.23−0.210.52+0.20−0.181.30+0.63−0.53<0.361.46+0.37−0.322.5+1.1−0.8
BAT constant1.04+0.12−0.110.83+0.11−0.101.25+0.28−0.281.24+0.16−0.140.78+0.13−0.12
NProb0.040.150.030.000.840.74
χ2ν713.4/648873.4/831606.1/540830.4/720765.7/809839.9/867
Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ1.99+0.04−0.041.87+0.04−0.041.89+0.04−0.042.00+0.04−0.021.80+0.01−0.012.19+0.05−0.05
Norma1.04+0.05−0.060.68+0.04−0.050.13+0.04−0.040.56+0.16−0.020.57+0.06−0.071.60+0.07−0.06
Plasma temperature (keV)<8.5<18.6<8.3<8.4
Plasma optical depth τ<2.4<2.1<2.20.8+0.2−0.5
Fluxb1.05+0.05−0.050.59+0.03−0.030.98+0.02−0.020.22+0.03−0.03
▵χ2250240955194
Neutral Fe Kα line (keV)6.40+0.02−0.026.39+0.01−0.016.43+0.06−0.086.39+0.03−0.026.40+0.02−0.026.37+0.05−0.05
σNarrow (eV)<84<45<110<72<43<190
EWNarrow (eV)40+11−1290+9−1055+28−5160+38−2371+38−25<35
Flux (10−5 ph cm−2 s−1)1.49+0.41−0.462.50+0.26−0.270.29+0.15−0.271.00+0.63−0.381.92+1.02−0.69<1.49
Fe xxv line (keV)6.74+0.04−0.046.69+0.06−0.05
EW (eV)17+8−933+13−7
Flux (10−5 ph cm−2 s−1)0.46+0.22−0.240.47+0.18−0.10
▵χ2714
Fe xxvi line (keV)6.96+0.04−0.046.98+0.02−0.026.94+0.06−0.096.98+0.06−0.14
EW (eV)24+11−1130+8−10<5121+13−13
Flux (10−5 ph cm−2 s−1)0.70+0.32−0.320.70+0.19−0.23<0.220.24+0.07−0.06
▵χ2122438
Broad line (keV)6.36+0.08−0.096.16+0.20−0.206.33+0.21−0.246.27+0.13−0.176.32+0.06−0.116.32+0.24−0.25
σBroad (keV)0.32+0.11−0.090.37+0.21−0.130.28+0.28−0.130.50+0.13−0.110.15+0.07−0.030.39+0.45−0.26
EWBroad (eV)105+26−2453+27−3892+59−56134+42−3862+28−4471+53−48
Flux (10−5 ph cm−2 s−1)3.74+0.93−0.861.64+0.83−1.170.49+0.31−0.302.28+0.72−0.651.55+0.70−1.103.43+2.56−2.32
▵χ257283532412
Rfrac0.88+0.23−0.210.52+0.20−0.181.30+0.63−0.53<0.361.46+0.37−0.322.5+1.1−0.8
BAT constant1.04+0.12−0.110.83+0.11−0.101.25+0.28−0.281.24+0.16−0.140.78+0.13−0.12
NProb0.040.150.030.000.840.74
χ2ν713.4/648873.4/831606.1/540830.4/720765.7/809839.9/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

b Flux for comptt quoted over the 0.5–10.0 keV range in units of 10−11 erg cm−2 s−1.

Table 4

Model A components for Suzaku XIS, HXD and BAT data from Swift.

Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ1.99+0.04−0.041.87+0.04−0.041.89+0.04−0.042.00+0.04−0.021.80+0.01−0.012.19+0.05−0.05
Norma1.04+0.05−0.060.68+0.04−0.050.13+0.04−0.040.56+0.16−0.020.57+0.06−0.071.60+0.07−0.06
Plasma temperature (keV)<8.5<18.6<8.3<8.4
Plasma optical depth τ<2.4<2.1<2.20.8+0.2−0.5
Fluxb1.05+0.05−0.050.59+0.03−0.030.98+0.02−0.020.22+0.03−0.03
▵χ2250240955194
Neutral Fe Kα line (keV)6.40+0.02−0.026.39+0.01−0.016.43+0.06−0.086.39+0.03−0.026.40+0.02−0.026.37+0.05−0.05
σNarrow (eV)<84<45<110<72<43<190
EWNarrow (eV)40+11−1290+9−1055+28−5160+38−2371+38−25<35
Flux (10−5 ph cm−2 s−1)1.49+0.41−0.462.50+0.26−0.270.29+0.15−0.271.00+0.63−0.381.92+1.02−0.69<1.49
Fe xxv line (keV)6.74+0.04−0.046.69+0.06−0.05
EW (eV)17+8−933+13−7
Flux (10−5 ph cm−2 s−1)0.46+0.22−0.240.47+0.18−0.10
▵χ2714
Fe xxvi line (keV)6.96+0.04−0.046.98+0.02−0.026.94+0.06−0.096.98+0.06−0.14
EW (eV)24+11−1130+8−10<5121+13−13
Flux (10−5 ph cm−2 s−1)0.70+0.32−0.320.70+0.19−0.23<0.220.24+0.07−0.06
▵χ2122438
Broad line (keV)6.36+0.08−0.096.16+0.20−0.206.33+0.21−0.246.27+0.13−0.176.32+0.06−0.116.32+0.24−0.25
σBroad (keV)0.32+0.11−0.090.37+0.21−0.130.28+0.28−0.130.50+0.13−0.110.15+0.07−0.030.39+0.45−0.26
EWBroad (eV)105+26−2453+27−3892+59−56134+42−3862+28−4471+53−48
Flux (10−5 ph cm−2 s−1)3.74+0.93−0.861.64+0.83−1.170.49+0.31−0.302.28+0.72−0.651.55+0.70−1.103.43+2.56−2.32
▵χ257283532412
Rfrac0.88+0.23−0.210.52+0.20−0.181.30+0.63−0.53<0.361.46+0.37−0.322.5+1.1−0.8
BAT constant1.04+0.12−0.110.83+0.11−0.101.25+0.28−0.281.24+0.16−0.140.78+0.13−0.12
NProb0.040.150.030.000.840.74
χ2ν713.4/648873.4/831606.1/540830.4/720765.7/809839.9/867
Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ1.99+0.04−0.041.87+0.04−0.041.89+0.04−0.042.00+0.04−0.021.80+0.01−0.012.19+0.05−0.05
Norma1.04+0.05−0.060.68+0.04−0.050.13+0.04−0.040.56+0.16−0.020.57+0.06−0.071.60+0.07−0.06
Plasma temperature (keV)<8.5<18.6<8.3<8.4
Plasma optical depth τ<2.4<2.1<2.20.8+0.2−0.5
Fluxb1.05+0.05−0.050.59+0.03−0.030.98+0.02−0.020.22+0.03−0.03
▵χ2250240955194
Neutral Fe Kα line (keV)6.40+0.02−0.026.39+0.01−0.016.43+0.06−0.086.39+0.03−0.026.40+0.02−0.026.37+0.05−0.05
σNarrow (eV)<84<45<110<72<43<190
EWNarrow (eV)40+11−1290+9−1055+28−5160+38−2371+38−25<35
Flux (10−5 ph cm−2 s−1)1.49+0.41−0.462.50+0.26−0.270.29+0.15−0.271.00+0.63−0.381.92+1.02−0.69<1.49
Fe xxv line (keV)6.74+0.04−0.046.69+0.06−0.05
EW (eV)17+8−933+13−7
Flux (10−5 ph cm−2 s−1)0.46+0.22−0.240.47+0.18−0.10
▵χ2714
Fe xxvi line (keV)6.96+0.04−0.046.98+0.02−0.026.94+0.06−0.096.98+0.06−0.14
EW (eV)24+11−1130+8−10<5121+13−13
Flux (10−5 ph cm−2 s−1)0.70+0.32−0.320.70+0.19−0.23<0.220.24+0.07−0.06
▵χ2122438
Broad line (keV)6.36+0.08−0.096.16+0.20−0.206.33+0.21−0.246.27+0.13−0.176.32+0.06−0.116.32+0.24−0.25
σBroad (keV)0.32+0.11−0.090.37+0.21−0.130.28+0.28−0.130.50+0.13−0.110.15+0.07−0.030.39+0.45−0.26
EWBroad (eV)105+26−2453+27−3892+59−56134+42−3862+28−4471+53−48
Flux (10−5 ph cm−2 s−1)3.74+0.93−0.861.64+0.83−1.170.49+0.31−0.302.28+0.72−0.651.55+0.70−1.103.43+2.56−2.32
▵χ257283532412
Rfrac0.88+0.23−0.210.52+0.20−0.181.30+0.63−0.53<0.361.46+0.37−0.322.5+1.1−0.8
BAT constant1.04+0.12−0.110.83+0.11−0.101.25+0.28−0.281.24+0.16−0.140.78+0.13−0.12
NProb0.040.150.030.000.840.74
χ2ν713.4/648873.4/831606.1/540830.4/720765.7/809839.9/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

b Flux for comptt quoted over the 0.5–10.0 keV range in units of 10−11 erg cm−2 s−1.

3.2 Model B – Self-consistent distant reflection

In an attempt to build a more self-consistent model, the pexrav and narrow 6.4 keV core are replaced with the reflionx ionized reflection model from Ross & Fabian (2005) which includes emission lines in addition to the reflection continuum. Model B provides fits of a reasonable quality to all six spectra; however, it is not expected to improve upon model A since it does not include a broad Gaussian to parametrize the red wing in an attempt to retain self-consistency and begin the construction of a more physically motivated model. Model B therefore acts as a ‘null hypothesis’ model in that all emission and reflection originate from distant matter. In most cases the fit statistic is somewhat worse than in model A, except for Fairall 9, in which the quality of the fit is actually improved to χ2ν = 864.7/834. The comptt model is still used here to model the soft excess in the four objects requiring it; however, the comptt parameters are not quoted in Table 5 (similarly in future tables) since the parameters obtained are consistent with those in model A.

Table 5

Model B components for Suzaku XIS, HXD and BAT data from Swift. The ionization parameter ξ is given in units of erg cm s−1.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.01−0.041.83+0.01−0.041.86+0.01−0.042.00+0.02−0.021.78+0.07−0.102.11+0.03−0.02
Norma0.95+0.06−0.060.64+0.02−0.040.12+0.04−0.030.51+0.01−0.020.54+0.04−0.061.52+0.04−0.05
ξ<22<20<2127+9−4<11<13
Fe/solar1.4+0.5−0.31.9+0.9−0.30.9+0.7−0.32.2+1.1−0.41.6+0.3−0.30.5+0.1−0.1
Normb0.89+0.42−0.511.11+0.12−0.620.27+0.06−0.110.16+0.05−0.051.18+0.10−0.202.05+0.33−1.00
Fe xxv line (keV)6.66+0.04−0.056.73+0.03−0.036.68+0.03−0.056.66+0.07−0.05
EW (eV)20+7−719+6−640+8−825+11−11
Flux (10−5 ph cm−2 s−1)0.72+0.25−0.250.52+0.16−0.160.64+0.13−0.130.94+0.41−0.41
▵χ220276716
Fe xxvi line (keV)6.95+0.03−0.036.98+0.02−0.026.94+0.05−0.116.96+0.06−0.166.98+0.09−0.26
EW (eV)31+9−931+7−739+24−2321+9−916+13−13
Flux (10−5 ph cm−2 s−1)0.92+0.27−0.270.71+0.16−−0.160.16+0.10−0.100.28+0.12−0.120.49+0.40−0.40
▵χ231387144
BAT constant1.07+0.12−0.120.83+0.10−0.101.03+0.23−0.231.03+0.10−0.090.70+0.10−0.10
NProb0.010.220.010.000.220.45
χ2ν741.8/649864.7/834623.6/543842.9/723840.2/812871.6/867
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.01−0.041.83+0.01−0.041.86+0.01−0.042.00+0.02−0.021.78+0.07−0.102.11+0.03−0.02
Norma0.95+0.06−0.060.64+0.02−0.040.12+0.04−0.030.51+0.01−0.020.54+0.04−0.061.52+0.04−0.05
ξ<22<20<2127+9−4<11<13
Fe/solar1.4+0.5−0.31.9+0.9−0.30.9+0.7−0.32.2+1.1−0.41.6+0.3−0.30.5+0.1−0.1
Normb0.89+0.42−0.511.11+0.12−0.620.27+0.06−0.110.16+0.05−0.051.18+0.10−0.202.05+0.33−1.00
Fe xxv line (keV)6.66+0.04−0.056.73+0.03−0.036.68+0.03−0.056.66+0.07−0.05
EW (eV)20+7−719+6−640+8−825+11−11
Flux (10−5 ph cm−2 s−1)0.72+0.25−0.250.52+0.16−0.160.64+0.13−0.130.94+0.41−0.41
▵χ220276716
Fe xxvi line (keV)6.95+0.03−0.036.98+0.02−0.026.94+0.05−0.116.96+0.06−0.166.98+0.09−0.26
EW (eV)31+9−931+7−739+24−2321+9−916+13−13
Flux (10−5 ph cm−2 s−1)0.92+0.27−0.270.71+0.16−−0.160.16+0.10−0.100.28+0.12−0.120.49+0.40−0.40
▵χ231387144
BAT constant1.07+0.12−0.120.83+0.10−0.101.03+0.23−0.231.03+0.10−0.090.70+0.10−0.10
NProb0.010.220.010.000.220.45
χ2ν741.8/649864.7/834623.6/543842.9/723840.2/812871.6/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization given in units of 10−5.

Table 5

Model B components for Suzaku XIS, HXD and BAT data from Swift. The ionization parameter ξ is given in units of erg cm s−1.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.01−0.041.83+0.01−0.041.86+0.01−0.042.00+0.02−0.021.78+0.07−0.102.11+0.03−0.02
Norma0.95+0.06−0.060.64+0.02−0.040.12+0.04−0.030.51+0.01−0.020.54+0.04−0.061.52+0.04−0.05
ξ<22<20<2127+9−4<11<13
Fe/solar1.4+0.5−0.31.9+0.9−0.30.9+0.7−0.32.2+1.1−0.41.6+0.3−0.30.5+0.1−0.1
Normb0.89+0.42−0.511.11+0.12−0.620.27+0.06−0.110.16+0.05−0.051.18+0.10−0.202.05+0.33−1.00
Fe xxv line (keV)6.66+0.04−0.056.73+0.03−0.036.68+0.03−0.056.66+0.07−0.05
EW (eV)20+7−719+6−640+8−825+11−11
Flux (10−5 ph cm−2 s−1)0.72+0.25−0.250.52+0.16−0.160.64+0.13−0.130.94+0.41−0.41
▵χ220276716
Fe xxvi line (keV)6.95+0.03−0.036.98+0.02−0.026.94+0.05−0.116.96+0.06−0.166.98+0.09−0.26
EW (eV)31+9−931+7−739+24−2321+9−916+13−13
Flux (10−5 ph cm−2 s−1)0.92+0.27−0.270.71+0.16−−0.160.16+0.10−0.100.28+0.12−0.120.49+0.40−0.40
▵χ231387144
BAT constant1.07+0.12−0.120.83+0.10−0.101.03+0.23−0.231.03+0.10−0.090.70+0.10−0.10
NProb0.010.220.010.000.220.45
χ2ν741.8/649864.7/834623.6/543842.9/723840.2/812871.6/867
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.01−0.041.83+0.01−0.041.86+0.01−0.042.00+0.02−0.021.78+0.07−0.102.11+0.03−0.02
Norma0.95+0.06−0.060.64+0.02−0.040.12+0.04−0.030.51+0.01−0.020.54+0.04−0.061.52+0.04−0.05
ξ<22<20<2127+9−4<11<13
Fe/solar1.4+0.5−0.31.9+0.9−0.30.9+0.7−0.32.2+1.1−0.41.6+0.3−0.30.5+0.1−0.1
Normb0.89+0.42−0.511.11+0.12−0.620.27+0.06−0.110.16+0.05−0.051.18+0.10−0.202.05+0.33−1.00
Fe xxv line (keV)6.66+0.04−0.056.73+0.03−0.036.68+0.03−0.056.66+0.07−0.05
EW (eV)20+7−719+6−640+8−825+11−11
Flux (10−5 ph cm−2 s−1)0.72+0.25−0.250.52+0.16−0.160.64+0.13−0.130.94+0.41−0.41
▵χ220276716
Fe xxvi line (keV)6.95+0.03−0.036.98+0.02−0.026.94+0.05−0.116.96+0.06−0.166.98+0.09−0.26
EW (eV)31+9−931+7−739+24−2321+9−916+13−13
Flux (10−5 ph cm−2 s−1)0.92+0.27−0.270.71+0.16−−0.160.16+0.10−0.100.28+0.12−0.120.49+0.40−0.40
▵χ231387144
BAT constant1.07+0.12−0.120.83+0.10−0.101.03+0.23−0.231.03+0.10−0.090.70+0.10−0.10
NProb0.010.220.010.000.220.45
χ2ν741.8/649864.7/834623.6/543842.9/723840.2/812871.6/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization given in units of 10−5.

The narrow ionized emission lines at 6.7 and 6.97 keV found in model A remain in model B and following the removal of the broad component, the significance of these lines is increased. In the case of SWIFT J2127.4+5654, for which there are no statistically significant ionized emission lines in model A, there are small excesses at these energies, indeed suggesting the possible presence of these lines. Emission lines of high significance in model A are also similarly significant in model B; for example, the Fe xxvi in Fairall 9 is particularly prominent in both A and B. This is contrary to the analysis of Fairall 9 by Schmoll et al. (2009) in which they did not detect the Fe xxvi line in the same Suzaku data, but they did however detect the He-like Fe xxv emission line which is also seen here in models A and B.

All objects show some indication of at least a small red wing (see Fig. 4), noting that in MCG −02-14-009 this feature is not particularly strong and in all objects the red wing does not extend below ∼5 keV. A smaller than expected excess at energies redwards of 6.4 keV is seen in Fairall 9, given the statistical significance of the broad Gaussian in model A. This could be due in part to a relatively high Fe abundance as measured by the Fe/solar parameter within reflionx at approximately twice the solar value. Some objects in the sample indicate a super-solar Fe abundance according to the reflionx component (see Table 5); however, the true abundance may in fact be lower than this measured value if part of the Fe Kα emission additionally arises from Compton-thin matter such as the broad-line region (BLR) or narrow-line region (NLR). A summary of the results obtained with model B is given in Table 5.

Ratio plots of model B (i.e. without including a broad Fe Kα) revealing excesses at energies redwards of 6.4 keV in some objects. Model B consisting of powerlaw +comptt (where required) + unblurred reflionx+ narrow ionized lines as required.
Figure 4

Ratio plots of model B (i.e. without including a broad Fe Kα) revealing excesses at energies redwards of 6.4 keV in some objects. Model B consisting of powerlaw +comptt (where required) + unblurred reflionx+ narrow ionized lines as required.

3.3 Model C –laor profile

The presence of emission from close-in to the BH has been indicated in the previous two models. In accordance with this, emission from the Fe K line region can be modelled with a relativistic component representing line emission from the accretion disc in addition to any appropriate emission from distant material (see Fig. 4). Model C advances upon model B by remodelling the Fe K region with the addition of a laor line profile operating under the assumption of a maximally spinning BH (Laor 1991). The line energy was restricted to 6.4–6.97 keV in the rest frame, with emissivity, inclination and the inner radius of emission allowed to vary. The outer radius of emission was fixed at 400 rg throughout. In the cases where the line emission reaches its lower limit, it is fixed at 6.4 keV in the rest frame. Whilst giving a more feasible interpretation of the broad emission in the Fe K region, it does not provide us with the most physically accurate representation given the assumption of a maximally rotating central BH with a = 0.998 within the laor model. Not all of the objects in the sample will feature a maximally rotating BH; indeed, Fairall 9 and SWIFT J2127.4+5654 have previously been found to have intermediate spin values (see Miniutti et al. 2009; Schmoll et al. 2009). The presence of narrow ionized emission lines was reassessed after accounting for the blue wing of the laor profile. Employment of this model therefore seeks to further parametrize the broad emission in the Fe K region, providing suitable and plausible parameters for use in the later kerrdisk models, whilst giving us an indication of the extent to which the spin of the central BH has an effect upon the observed spectrum.

The fit to all objects is improved over the purely distant emission in model B. The 6.97-keV line is found to be present in all but NGC 7469 and SWIFT J2127.4+5654 whilst the 6.7-keV line is only found in Mrk 335 (although relatively weak). Model C improves the fit to MCG −02-14-009 the least with only ▵χ2∼ 14 for three additional free parameters whereas the introduction of a laor profile offers a significant improvement for most other objects, particularly Fairall 9, Mrk 335 and SWIFT J2127.4+5654 (▵χ2∼ 33, 37 and 34, respectively) (see Table 6).

Table 6

Fit parameters from model C to Suzaku XIS, HXD and BAT data from Swift. Line energies are quoted in the rest frame. * denotes a frozen parameter for cases where the emissivity index is fixed at q = 3, it is unconstrained. The improvement ▵χ2 in the fit with the introduction of a laor profile is noted in comparison with the purely distant reflection as present in this model. Note that the Fe xxv emission line is no longer required in Ark 120, Fairall 9 and SWIFT J2127.4+5654.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.91+0.03−0.031.84+0.02−0.031.86+0.03−0.032.01+0.02−0.021.78+0.01−0.012.13+0.03−0.03
Norma0.97+0.03−0.050.64+0.03−0.020.12+0.04−0.030.56+0.02−0.020.55+0.04−0.061.54+0.04−0.05
laor profile
Line E (keV)6.4*<6.566.4*6.4*<6.446.4*
EW (eV)85+21−1779+21−22119+58−52126+33−3560+26−22140+41−41
q3.0+1.6−1.03.5+2.9−0.93.0*3.9+1.9−1.03.0*2.5+0.7−0.6
Rin (GM/c2)25+19−716+17−6>1332+12−1681+82−37<33
i°34+7−534+5−330+17−938+10−824+12−846+9−9
Flux (10−5 ph cm−2 s−1)3.85+0.95−0.771.79+0.50−0.500.59+0.30−0.262.28+0.62−0.642.56+1.02−0.845.39+1.58−1.58
▵χ2243314372034
ξ<19<21<1533+18−9<13<15
Fe/solar1.0*1.7+0.3−0.30.6+0.6−0.42.0+0.9−0.8<1.40.4+0.1−0.1
Normb1.06+0.15−0.570.97+0.23−0.320.23+0.07−0.110.11+0.06−0.051.07+0.08−0.481.96+0.24−0.92
Fe xxv line (keV)6.65+0.09−0.07
EW (eV)10+8−8
Flux (10−5 ph cm−2 s−1)0.17+0.13−0.13
▵χ22
Fe xxvi line (keV)6.96+0.04−0.046.99+0.03−0.036.94+0.04−0.067.02+0.05−0.06
EW (eV)33+12−1927−9−1542+23−2419+9−9
Flux (10−5 ph cm−2 s−1)0.98+0.34−0.570.63+0.21−0.350.18+0.10−0.100.25+0.12−0.12
▵χ21013511
BAT constant1.10+0.12−0.120.86+0.10−0.101.11+0.25−0.251.08+0.11−0.090.73+0.10−0.10
χ2ν721.6/648857.7/831609.6/540806.3/719820.3/808855.8/867
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.91+0.03−0.031.84+0.02−0.031.86+0.03−0.032.01+0.02−0.021.78+0.01−0.012.13+0.03−0.03
Norma0.97+0.03−0.050.64+0.03−0.020.12+0.04−0.030.56+0.02−0.020.55+0.04−0.061.54+0.04−0.05
laor profile
Line E (keV)6.4*<6.566.4*6.4*<6.446.4*
EW (eV)85+21−1779+21−22119+58−52126+33−3560+26−22140+41−41
q3.0+1.6−1.03.5+2.9−0.93.0*3.9+1.9−1.03.0*2.5+0.7−0.6
Rin (GM/c2)25+19−716+17−6>1332+12−1681+82−37<33
i°34+7−534+5−330+17−938+10−824+12−846+9−9
Flux (10−5 ph cm−2 s−1)3.85+0.95−0.771.79+0.50−0.500.59+0.30−0.262.28+0.62−0.642.56+1.02−0.845.39+1.58−1.58
▵χ2243314372034
ξ<19<21<1533+18−9<13<15
Fe/solar1.0*1.7+0.3−0.30.6+0.6−0.42.0+0.9−0.8<1.40.4+0.1−0.1
Normb1.06+0.15−0.570.97+0.23−0.320.23+0.07−0.110.11+0.06−0.051.07+0.08−0.481.96+0.24−0.92
Fe xxv line (keV)6.65+0.09−0.07
EW (eV)10+8−8
Flux (10−5 ph cm−2 s−1)0.17+0.13−0.13
▵χ22
Fe xxvi line (keV)6.96+0.04−0.046.99+0.03−0.036.94+0.04−0.067.02+0.05−0.06
EW (eV)33+12−1927−9−1542+23−2419+9−9
Flux (10−5 ph cm−2 s−1)0.98+0.34−0.570.63+0.21−0.350.18+0.10−0.100.25+0.12−0.12
▵χ21013511
BAT constant1.10+0.12−0.120.86+0.10−0.101.11+0.25−0.251.08+0.11−0.090.73+0.10−0.10
χ2ν721.6/648857.7/831609.6/540806.3/719820.3/808855.8/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization given in units of 10−5.

Table 6

Fit parameters from model C to Suzaku XIS, HXD and BAT data from Swift. Line energies are quoted in the rest frame. * denotes a frozen parameter for cases where the emissivity index is fixed at q = 3, it is unconstrained. The improvement ▵χ2 in the fit with the introduction of a laor profile is noted in comparison with the purely distant reflection as present in this model. Note that the Fe xxv emission line is no longer required in Ark 120, Fairall 9 and SWIFT J2127.4+5654.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.91+0.03−0.031.84+0.02−0.031.86+0.03−0.032.01+0.02−0.021.78+0.01−0.012.13+0.03−0.03
Norma0.97+0.03−0.050.64+0.03−0.020.12+0.04−0.030.56+0.02−0.020.55+0.04−0.061.54+0.04−0.05
laor profile
Line E (keV)6.4*<6.566.4*6.4*<6.446.4*
EW (eV)85+21−1779+21−22119+58−52126+33−3560+26−22140+41−41
q3.0+1.6−1.03.5+2.9−0.93.0*3.9+1.9−1.03.0*2.5+0.7−0.6
Rin (GM/c2)25+19−716+17−6>1332+12−1681+82−37<33
i°34+7−534+5−330+17−938+10−824+12−846+9−9
Flux (10−5 ph cm−2 s−1)3.85+0.95−0.771.79+0.50−0.500.59+0.30−0.262.28+0.62−0.642.56+1.02−0.845.39+1.58−1.58
▵χ2243314372034
ξ<19<21<1533+18−9<13<15
Fe/solar1.0*1.7+0.3−0.30.6+0.6−0.42.0+0.9−0.8<1.40.4+0.1−0.1
Normb1.06+0.15−0.570.97+0.23−0.320.23+0.07−0.110.11+0.06−0.051.07+0.08−0.481.96+0.24−0.92
Fe xxv line (keV)6.65+0.09−0.07
EW (eV)10+8−8
Flux (10−5 ph cm−2 s−1)0.17+0.13−0.13
▵χ22
Fe xxvi line (keV)6.96+0.04−0.046.99+0.03−0.036.94+0.04−0.067.02+0.05−0.06
EW (eV)33+12−1927−9−1542+23−2419+9−9
Flux (10−5 ph cm−2 s−1)0.98+0.34−0.570.63+0.21−0.350.18+0.10−0.100.25+0.12−0.12
▵χ21013511
BAT constant1.10+0.12−0.120.86+0.10−0.101.11+0.25−0.251.08+0.11−0.090.73+0.10−0.10
χ2ν721.6/648857.7/831609.6/540806.3/719820.3/808855.8/867
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.91+0.03−0.031.84+0.02−0.031.86+0.03−0.032.01+0.02−0.021.78+0.01−0.012.13+0.03−0.03
Norma0.97+0.03−0.050.64+0.03−0.020.12+0.04−0.030.56+0.02−0.020.55+0.04−0.061.54+0.04−0.05
laor profile
Line E (keV)6.4*<6.566.4*6.4*<6.446.4*
EW (eV)85+21−1779+21−22119+58−52126+33−3560+26−22140+41−41
q3.0+1.6−1.03.5+2.9−0.93.0*3.9+1.9−1.03.0*2.5+0.7−0.6
Rin (GM/c2)25+19−716+17−6>1332+12−1681+82−37<33
i°34+7−534+5−330+17−938+10−824+12−846+9−9
Flux (10−5 ph cm−2 s−1)3.85+0.95−0.771.79+0.50−0.500.59+0.30−0.262.28+0.62−0.642.56+1.02−0.845.39+1.58−1.58
▵χ2243314372034
ξ<19<21<1533+18−9<13<15
Fe/solar1.0*1.7+0.3−0.30.6+0.6−0.42.0+0.9−0.8<1.40.4+0.1−0.1
Normb1.06+0.15−0.570.97+0.23−0.320.23+0.07−0.110.11+0.06−0.051.07+0.08−0.481.96+0.24−0.92
Fe xxv line (keV)6.65+0.09−0.07
EW (eV)10+8−8
Flux (10−5 ph cm−2 s−1)0.17+0.13−0.13
▵χ22
Fe xxvi line (keV)6.96+0.04−0.046.99+0.03−0.036.94+0.04−0.067.02+0.05−0.06
EW (eV)33+12−1927−9−1542+23−2419+9−9
Flux (10−5 ph cm−2 s−1)0.98+0.34−0.570.63+0.21−0.350.18+0.10−0.100.25+0.12−0.12
▵χ21013511
BAT constant1.10+0.12−0.120.86+0.10−0.101.11+0.25−0.251.08+0.11−0.090.73+0.10−0.10
χ2ν721.6/648857.7/831609.6/540806.3/719820.3/808855.8/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization given in units of 10−5.

None of the objects in the sample require emission from within 6 GM/c2 indicating that a rotating BH, whilst possible, is not required to model the spectra. According to the fit parameters obtained with laor, the inner radius of emission lies at tens of rg for all six AGN. In most cases, the accretion disc is unlikely to be truncated at these distances from the BH and it is likely that these values arise from the assumption of a maximally rotating BH within the laor model. The emissivity indices also indicate that a high concentration of emission from very close to the BH is absent in the spectra from these AGN. An emissivity index of q > 5 would suggest that the emission from the accretion disc is very centrally concentrated (e.g. Miniutti et al. 2003), while for objects such as Mrk 335 and Fairall 9 this is within the error bars; one would also expect the inner radius of emission to be much closer to the BH if this interpretation were suitable for the spectra. For the cases where the emissivity index cannot be constrained, it is fixed to q = 3.0 which is consistent with the other objects in the sample.

Mrk 335 appears to have one of the strongest relativistic lines in the sample with a laor profile equivalent width of EW = 126+33−35 eV and an improvement of ▵χ2∼ 37. However, this is not particularly strong in comparison with previous studies of this AGN; for example, Longinotti et al. (2007) found an EW = 320+170−100 eV at a line energy of E = 6.93+0.77−0.27 keV in a 40-ks XMM–Newton observation of Mrk 335 in 2000, also using a laor profile. The width of the iron line emission from the disc and its suggested high ionization state could be due to the lack of an Fe xxvi line in the Longinotti spectra (possibly due to a low S/N), which has been accounted for in this analysis and is noted by O’Neill et al. (2007) in an analysis of the same XMM–Newton observation used here. Also in agreement with our results, O’Neill et al. found the equivalent width of the broad line in Mrk 335 to be EW = 115+14−14 eV in comparison with EW = 126+33−35 eV for the laor profile and EW = 113+46−52 eV for the broad Gaussian employed in model A.

The results obtained here are consistent with a recent analysis of Ark 120 by Nardini et al. (2010) who initially fitted the features in the Fe K region with a laor profile. In the Nardini et al. (2010) analysis, both the inclination of the accretion disc and emissivity are frozen at typical values of i = 40° and q = 3.0, respectively, finding an inner radius of emission of rin = 13+19−7rg. These values are consistent with those obtained here with model C, finding i° = 34+7−5, q = 3.0+1.6−1.0 and rin = 25+19−7rg. Also detected is the Fe xxvi narrow ionized emission line, again consistent with Nardini et al. (2010).

As in model A, there are no narrow ionized emission lines found in SWIFT J2127.4+5654 corresponding to Fe xxv and Fe xxvi. This implies that the excesses modelled as narrow components in model B may instead be due to relativistically broadened Fe Kα emission. When modelled as a single laor profile, all significant excesses in the Fe K band are removed with the profile centroid rest-frame energy at 6.4 keV and the blue wing peaking at ∼6.7 keV in the rest frame. Consequently, the observed excess at 6.7 keV observed in model A and more significantly in model B is more convincingly described as a blue wing. Additionally, the observed intermediate feature in MCG −02-14-009 is modelled well by the blue wing of the laor profile which occurs at ∼6.58 keV, prompting the removal of the previously employed narrow Gaussian.

3.4 Model D –kerrdisk profile

Replacing the laor model with a kerrdisk (Brenneman & Reynolds 2006) line profile further allows a more physically motivated fit to the data. The kerrdisk model allows the spin parameter to be varied between 0 < a < 0.998 with the aim of determining the extent to which we can rule in or out a non-rotating or maximally spinning BH. Similarly to model C, the outer radius of the disc is fixed at 400 rms, while assuming a disc of uniform emissivity. Throughout the fits using model D, the inner radius of emission is assumed to extend down to the ISCO (rISCO). The line energy is confined to 6.4–6.97 keV in the rest frame, being frozen at 6.4 keV if it reaches its lower limit.

Model D provides good fits to all six objects, producing the best-fitting physically motivated model for Mrk 335 and NGC 7469 prior to considering a blurred reflection component from the inner regions of the accretion disc (i.e. models E and F) (see Table 7 and Fig. 5). The presence of ionized emission lines is entirely consistent with model C, adding weight to the probability that these lines are present in the spectra. Further to this, the inclinations of the accretion discs to the observer are all comparable to those obtained in model C, as are the equivalent widths of the relativistic line profiles. Given that the spin of the central BH is a free parameter and the assumptions within this model, the emissivity indices are expected to vary from those obtained with the laor profile since the measured emissivity is degenerate to some extent with the inner radius of emission. Therefore as the spin parameter varies so does rISCO, consequently affecting the measured q.

Table 7

Fit parameters from model D to Suzaku XIS, HXD and BAT data from Swift. Line energies are quoted in the rest frame. * denotes a frozen parameter.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.03−0.031.84+0.02−0.031.86+0.01−0.012.02+0.03−0.021.78+0.01−0.012.13+0.04−0.03
Norma0.97+0.04−0.050.64+0.01−0.040.13+0.04−0.040.56+0.01−0.010.55+0.04−0.061.55+0.05−0.05
kerrdisk profile
Line E (keV)6.4*6.4*6.47+0.03−0.036.4*6.4*6.4*
EW (eV)95+32−2663+36−19142+47−46146+39−3991+9−8178+82−69
q2.3+0.4−0.32.7+0.7−0.42.0+0.4−0.42.6+0.5−0.31.7+0.4−0.62.6+1.0−0.4
a<0.940.44+0.04−0.11<0.880.70+0.12−0.010.69+0.09−0.090.70+0.10−0.14
i°33+2−1738+8−524+10−938+2−223+15−743+5−10
Flux (10−5 ph cm−2 s−1)3.62+1.22−0.991.78+1.02−0.546.64+2.20−2.152.52+0.67−0.672.23+0.22−0.206.35+2.91−2.46
▵χ2 (kerrdisk)233012402537
▵χ2 (zero spin)1404107
ξ<16<17<1641+7−16<11<13
Fe/solar1.0*1.6+0.3−0.40.8+0.6−0.41.9+0.9−0.60.8+0.4−0.20.3+0.2−0.2
Normb0.97+0.50−0.641.09+0.09−0.290.25+0.07−0.160.07+0.07−0.010.69+0.22−0.152.12+0.63−0.83
Fe xxv line (keV)6.67+0.06−0.05
EW (eV)16+8−8
Flux (10−5 ph cm−2 s−1)0.26+0.13−0.13
▵χ24
Fe xxvi line (keV)6.99+0.04−0.066.98+0.04−0.036.94+0.08−0.09formula
EW (eV)22+23−1530+8−1044+26−3822+10−10
Flux (10−5 ph cm−2 s−1)0.65+0.68−0.440.70+0.19−0.240.18+0.11−0.160.28+0.120.12
▵χ2712914
BAT constant1.05+0.12−0.120.85+0.10−0.101.05+0.23−0.231.06+0.11−0.090.67+0.09−0.09
χ2ν724.4/648861.2/832611.6/538803.0/719815.5/808852.3/867
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.03−0.031.84+0.02−0.031.86+0.01−0.012.02+0.03−0.021.78+0.01−0.012.13+0.04−0.03
Norma0.97+0.04−0.050.64+0.01−0.040.13+0.04−0.040.56+0.01−0.010.55+0.04−0.061.55+0.05−0.05
kerrdisk profile
Line E (keV)6.4*6.4*6.47+0.03−0.036.4*6.4*6.4*
EW (eV)95+32−2663+36−19142+47−46146+39−3991+9−8178+82−69
q2.3+0.4−0.32.7+0.7−0.42.0+0.4−0.42.6+0.5−0.31.7+0.4−0.62.6+1.0−0.4
a<0.940.44+0.04−0.11<0.880.70+0.12−0.010.69+0.09−0.090.70+0.10−0.14
i°33+2−1738+8−524+10−938+2−223+15−743+5−10
Flux (10−5 ph cm−2 s−1)3.62+1.22−0.991.78+1.02−0.546.64+2.20−2.152.52+0.67−0.672.23+0.22−0.206.35+2.91−2.46
▵χ2 (kerrdisk)233012402537
▵χ2 (zero spin)1404107
ξ<16<17<1641+7−16<11<13
Fe/solar1.0*1.6+0.3−0.40.8+0.6−0.41.9+0.9−0.60.8+0.4−0.20.3+0.2−0.2
Normb0.97+0.50−0.641.09+0.09−0.290.25+0.07−0.160.07+0.07−0.010.69+0.22−0.152.12+0.63−0.83
Fe xxv line (keV)6.67+0.06−0.05
EW (eV)16+8−8
Flux (10−5 ph cm−2 s−1)0.26+0.13−0.13
▵χ24
Fe xxvi line (keV)6.99+0.04−0.066.98+0.04−0.036.94+0.08−0.09formula
EW (eV)22+23−1530+8−1044+26−3822+10−10
Flux (10−5 ph cm−2 s−1)0.65+0.68−0.440.70+0.19−0.240.18+0.11−0.160.28+0.120.12
▵χ2712914
BAT constant1.05+0.12−0.120.85+0.10−0.101.05+0.23−0.231.06+0.11−0.090.67+0.09−0.09
χ2ν724.4/648861.2/832611.6/538803.0/719815.5/808852.3/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization given in units of 10−5.

Table 7

Fit parameters from model D to Suzaku XIS, HXD and BAT data from Swift. Line energies are quoted in the rest frame. * denotes a frozen parameter.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.03−0.031.84+0.02−0.031.86+0.01−0.012.02+0.03−0.021.78+0.01−0.012.13+0.04−0.03
Norma0.97+0.04−0.050.64+0.01−0.040.13+0.04−0.040.56+0.01−0.010.55+0.04−0.061.55+0.05−0.05
kerrdisk profile
Line E (keV)6.4*6.4*6.47+0.03−0.036.4*6.4*6.4*
EW (eV)95+32−2663+36−19142+47−46146+39−3991+9−8178+82−69
q2.3+0.4−0.32.7+0.7−0.42.0+0.4−0.42.6+0.5−0.31.7+0.4−0.62.6+1.0−0.4
a<0.940.44+0.04−0.11<0.880.70+0.12−0.010.69+0.09−0.090.70+0.10−0.14
i°33+2−1738+8−524+10−938+2−223+15−743+5−10
Flux (10−5 ph cm−2 s−1)3.62+1.22−0.991.78+1.02−0.546.64+2.20−2.152.52+0.67−0.672.23+0.22−0.206.35+2.91−2.46
▵χ2 (kerrdisk)233012402537
▵χ2 (zero spin)1404107
ξ<16<17<1641+7−16<11<13
Fe/solar1.0*1.6+0.3−0.40.8+0.6−0.41.9+0.9−0.60.8+0.4−0.20.3+0.2−0.2
Normb0.97+0.50−0.641.09+0.09−0.290.25+0.07−0.160.07+0.07−0.010.69+0.22−0.152.12+0.63−0.83
Fe xxv line (keV)6.67+0.06−0.05
EW (eV)16+8−8
Flux (10−5 ph cm−2 s−1)0.26+0.13−0.13
▵χ24
Fe xxvi line (keV)6.99+0.04−0.066.98+0.04−0.036.94+0.08−0.09formula
EW (eV)22+23−1530+8−1044+26−3822+10−10
Flux (10−5 ph cm−2 s−1)0.65+0.68−0.440.70+0.19−0.240.18+0.11−0.160.28+0.120.12
▵χ2712914
BAT constant1.05+0.12−0.120.85+0.10−0.101.05+0.23−0.231.06+0.11−0.090.67+0.09−0.09
χ2ν724.4/648861.2/832611.6/538803.0/719815.5/808852.3/867
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
Γ1.90+0.03−0.031.84+0.02−0.031.86+0.01−0.012.02+0.03−0.021.78+0.01−0.012.13+0.04−0.03
Norma0.97+0.04−0.050.64+0.01−0.040.13+0.04−0.040.56+0.01−0.010.55+0.04−0.061.55+0.05−0.05
kerrdisk profile
Line E (keV)6.4*6.4*6.47+0.03−0.036.4*6.4*6.4*
EW (eV)95+32−2663+36−19142+47−46146+39−3991+9−8178+82−69
q2.3+0.4−0.32.7+0.7−0.42.0+0.4−0.42.6+0.5−0.31.7+0.4−0.62.6+1.0−0.4
a<0.940.44+0.04−0.11<0.880.70+0.12−0.010.69+0.09−0.090.70+0.10−0.14
i°33+2−1738+8−524+10−938+2−223+15−743+5−10
Flux (10−5 ph cm−2 s−1)3.62+1.22−0.991.78+1.02−0.546.64+2.20−2.152.52+0.67−0.672.23+0.22−0.206.35+2.91−2.46
▵χ2 (kerrdisk)233012402537
▵χ2 (zero spin)1404107
ξ<16<17<1641+7−16<11<13
Fe/solar1.0*1.6+0.3−0.40.8+0.6−0.41.9+0.9−0.60.8+0.4−0.20.3+0.2−0.2
Normb0.97+0.50−0.641.09+0.09−0.290.25+0.07−0.160.07+0.07−0.010.69+0.22−0.152.12+0.63−0.83
Fe xxv line (keV)6.67+0.06−0.05
EW (eV)16+8−8
Flux (10−5 ph cm−2 s−1)0.26+0.13−0.13
▵χ24
Fe xxvi line (keV)6.99+0.04−0.066.98+0.04−0.036.94+0.08−0.09formula
EW (eV)22+23−1530+8−1044+26−3822+10−10
Flux (10−5 ph cm−2 s−1)0.65+0.68−0.440.70+0.19−0.240.18+0.11−0.160.28+0.120.12
▵χ2712914
BAT constant1.05+0.12−0.120.85+0.10−0.101.05+0.23−0.231.06+0.11−0.090.67+0.09−0.09
χ2ν724.4/648861.2/832611.6/538803.0/719815.5/808852.3/867

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization given in units of 10−5.

νFν plots of model D indicating the strength of any soft excess and relativistic line emission from the kerrdisk models. Consisting of powerlaw +comptt (where required) + unblurred reflionx+kerrdisk+ narrow ionized lines as required. XIS data are in black, HXD in red and BAT data are represented by blue circles.
Figure 5

νFν plots of model D indicating the strength of any soft excess and relativistic line emission from the kerrdisk models. Consisting of powerlaw +comptt (where required) + unblurred reflionx+kerrdisk+ narrow ionized lines as required. XIS data are in black, HXD in red and BAT data are represented by blue circles.

The results obtained from the spin parameter a for these six objects suggest that a maximally spinning central BH can be ruled out at the 90 per cent confidence level; however for some objects such as Ark 120 and MCG −02-14-009, the spin is unconstrained and only an upper limit can be measured (a < 0.94 and a < 0.88, suggesting that emission does not occur within 2.02 rg and 2.45 rg, respectively, consistent with the Nardini et al. (2010) analysis of Ark 120).

Here for Fairall 9 the best-fitting model gives a measured spin value of a = 0.44+0.04−0.11. A previous spin constraint for Fairall 9 by Schmoll et al. (2009) found a = 0.65+0.05−0.05 using a blurred reflector model (kerrconv (Brenneman & Reynolds 2006), convolved with reflionx). The measurement here also gives an intermediate value; however, it is only consistent with their findings when they ignore the spectra below 2 keV to ensure that the soft excess is not the component driving the main part of the fit. Schmoll et al. (2009) quote a worse spin constraint when these conditions are upheld a = 0.5+0.1−0.3. Within error bars at the 90 per cent level, these results are consistent with the value of the spin parameter found here. However, here we find that the emissivity index is constrained to q = 2.7+0.7−0.4 whereas Schmoll et al. only constrain this to q > 4.9.

Mrk 335 has been noted previously as an object with a broad relativistic Fe K line (Gondoin et al. 2002). Here we measure an intermediate spin a = 0.70+0.12−0.01, ruling out maximally spinning and non-rotating BHs only at the 95 per cent confidence level. The emissivity index is measured at a moderate value of q = 2.6+0.5−0.3 in agreement with previous studies of this object, for example by Longinotti et al. (2007). In agreement with model C, the kerrdisk component also features a relatively broad equivalent width EW = 146+39−39 with an improvement of ▵χ2∼ 40 (for four additional free parameters) over purely distant emission.

In NGC 7469 we measure a broad Fe line featuring a low emissivity index of q = 1.7+0.4−0.6, an inclination of i° = 23+15−7 and EW = 91+9−8 eV, providing good agreement with the previous measurements made within model C. We obtain a BH spin of a = 0.69+0.09−0.09, ruling out a Schwarzschild BH with 95 per cent confidence and a maximally rotating BH at the 99 per cent confidence level.

SWIFT J2127.4+5654 also shows a broad and statistically significant component, EW = 178+82−69 eV, improving the fit by ▵χ2∼ 37 (similarly to model C) for four additional free parameters. The emissivity and the inclination of the accretion disc are typical of the objects in this sample. The spin parameter measured here is consistent with that measured by Miniutti et al. (2009), 0.70+0.10−0.14 compared to 0.6+0.2−0.2 found using a blurred reflection model. Here we reject a maximally rotating central BH at greater than 99 per cent confidence and a non-rotating BH with 95 per cent confidence. The spin constraints obtained with model D for Mrk 335, NGC 7469 and SWIFT J2127.4+5654 can be seen in Fig. 6.

Confidence plots from model D for the spin parameter a for those objects in which it could be constrained. Dashed lines represent 90 per cent (black), 95 per cent (red), 99 per cent (green) and 99.5 per cent (blue) confidence levels. Note that the emissivity index q is a free parameter.
Figure 6

Confidence plots from model D for the spin parameter a for those objects in which it could be constrained. Dashed lines represent 90 per cent (black), 95 per cent (red), 99 per cent (green) and 99.5 per cent (blue) confidence levels. Note that the emissivity index q is a free parameter.

3.5 Model E – Blurred reflector (1)

In contrast to the previous models, model E does not include a contribution from the Comptonization of soft photons (i.e. the comptt model). Alternatively, the excess seen at soft energies is modelled using a blurred reflector in addition to reflection consistent with distant emission. The ionized reflection model reflionx is convolved with a kerrdisk kernel, modelling a relativistically blurred reflection spectrum from the accretion disc originating from very near to the BH instead of the kerrdisk model, i.e. the kerrconv model (Brenneman & Reynolds 2006). The reflection components consist of (reflionx+kerrconv*reflionx). The assumptions in using this convolution model are consistent with those used in model D; the inner radius of the accretion is assumed to extend down to rISCO and outer radius fixed at 400 rms for a disc of uniform emissivity. Narrow emission lines from ionized matter are also included where required, while the narrow 6.4-keV core is accounted for by the second unblurred reflionx component.

Model E is a significantly worse fit in all four AGN with a soft excess (by an average of ▵χ2∼ 23 compared to model D), requiring a significant amount of blurring to account for the excess at low energies (see Table 8). Given the amount of blurring required, parameters such as the emissivity index and the spin parameter approach more extreme values for all objects with a soft excess (e.g. q > 4 and a > 0.9). Whilst model E is a good fit to the data for Mrk 335 and NGC 7469, model D provides a better fit and without the use of such extreme parameters.

Table 8

Fit parameters from model E to Suzaku XIS, HXD and BAT data from Swift. (A) represents the unblurred reflionx and (B) represents the blurred reflionx. The ionization parameter ξ is given in units of erg cm s−1. * denotes a parameter frozen at the best-fitting value from model D. In some cases, the spin parameter a could not be constrained, denoted by –.

Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ2.12+0.01−0.022.01+0.01−0.011.90+0.03−0.032.15+0.01−0.021.80+0.01−0.012.20+0.05−0.05
Norma1.13+0.02−0.020.78+0.02−0.020.12+0.05−0.030.57+0.02−0.020.55+0.05−0.071.65+0.05−0.08
kerrconv
q4.1+0.8−0.8>6.15<3.56.6+2.0−1.0>3.02.2+0.6−0.4
a>0.970.98+0.01−0.01<0.960.87+0.05−0.06<0.97
i°48+3−670+5−240+11−1453+6−670+4−343+18−7
Fe/solar (A)1.4*1.9*0.4*2.0+0.8−0.41.6*>0.2
ξ (A)<13<11<543<11<11<26
Norm (A)b1.31+0.14−0.151.13+0.08−0.190.10+0.17−0.100.47+0.05−0.101.13+0.11−0.210.67+0.60−0.60
Fe/solar (B)1.5+0.3−0.30.8+0.2−0.3>1.11.0+0.1−0.1<0.41.0+0.8−0.3
ξ (B)56+12−324+17−9<14207+5−5<24<14
Norm (B)b0.38+0.21−0.130.24+0.04−0.030.27+0.19−0.200.06+0.01−0.010.30+0.36−0.123.68+1.50−1.99
Fe xxv line (keV)6.65+0.05−0.076.74+0.04−0.046.67+0.03−0.03
EW (eV)17+7−716+6−634+8−8
Flux (10−5 ph cm−2 s−1)0.62+0.26−0.270.44+0.17−0.170.55−0.13−0.13
▵χ281010
Fe xxvi line (keV)6.95+0.03−0.036.88+0.02−0.036.95+0.06−0.076.97+0.07−0.12
EW (eV)29+9−927+7−734+23−2316+9−8
Flux (10−5 ph cm−2 s−1)0.85+0.27−0.270.63+0.17−0.170.14+0.10−0.100.23+0.13−0.12
▵χ21625410
BAT constant1.11+0.11−0.110.87+0.10−0.101.06+0.23−0.230.93+0.08−0.080.68+0.10−0.09
χ2ν770.2/647879.6/832612.0/538811.6/722835.4/808830.4/865
Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ2.12+0.01−0.022.01+0.01−0.011.90+0.03−0.032.15+0.01−0.021.80+0.01−0.012.20+0.05−0.05
Norma1.13+0.02−0.020.78+0.02−0.020.12+0.05−0.030.57+0.02−0.020.55+0.05−0.071.65+0.05−0.08
kerrconv
q4.1+0.8−0.8>6.15<3.56.6+2.0−1.0>3.02.2+0.6−0.4
a>0.970.98+0.01−0.01<0.960.87+0.05−0.06<0.97
i°48+3−670+5−240+11−1453+6−670+4−343+18−7
Fe/solar (A)1.4*1.9*0.4*2.0+0.8−0.41.6*>0.2
ξ (A)<13<11<543<11<11<26
Norm (A)b1.31+0.14−0.151.13+0.08−0.190.10+0.17−0.100.47+0.05−0.101.13+0.11−0.210.67+0.60−0.60
Fe/solar (B)1.5+0.3−0.30.8+0.2−0.3>1.11.0+0.1−0.1<0.41.0+0.8−0.3
ξ (B)56+12−324+17−9<14207+5−5<24<14
Norm (B)b0.38+0.21−0.130.24+0.04−0.030.27+0.19−0.200.06+0.01−0.010.30+0.36−0.123.68+1.50−1.99
Fe xxv line (keV)6.65+0.05−0.076.74+0.04−0.046.67+0.03−0.03
EW (eV)17+7−716+6−634+8−8
Flux (10−5 ph cm−2 s−1)0.62+0.26−0.270.44+0.17−0.170.55−0.13−0.13
▵χ281010
Fe xxvi line (keV)6.95+0.03−0.036.88+0.02−0.036.95+0.06−0.076.97+0.07−0.12
EW (eV)29+9−927+7−734+23−2316+9−8
Flux (10−5 ph cm−2 s−1)0.85+0.27−0.270.63+0.17−0.170.14+0.10−0.100.23+0.13−0.12
▵χ21625410
BAT constant1.11+0.11−0.110.87+0.10−0.101.06+0.23−0.230.93+0.08−0.080.68+0.10−0.09
χ2ν770.2/647879.6/832612.0/538811.6/722835.4/808830.4/865

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization in units of 10−5.

Table 8

Fit parameters from model E to Suzaku XIS, HXD and BAT data from Swift. (A) represents the unblurred reflionx and (B) represents the blurred reflionx. The ionization parameter ξ is given in units of erg cm s−1. * denotes a parameter frozen at the best-fitting value from model D. In some cases, the spin parameter a could not be constrained, denoted by –.

Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ2.12+0.01−0.022.01+0.01−0.011.90+0.03−0.032.15+0.01−0.021.80+0.01−0.012.20+0.05−0.05
Norma1.13+0.02−0.020.78+0.02−0.020.12+0.05−0.030.57+0.02−0.020.55+0.05−0.071.65+0.05−0.08
kerrconv
q4.1+0.8−0.8>6.15<3.56.6+2.0−1.0>3.02.2+0.6−0.4
a>0.970.98+0.01−0.01<0.960.87+0.05−0.06<0.97
i°48+3−670+5−240+11−1453+6−670+4−343+18−7
Fe/solar (A)1.4*1.9*0.4*2.0+0.8−0.41.6*>0.2
ξ (A)<13<11<543<11<11<26
Norm (A)b1.31+0.14−0.151.13+0.08−0.190.10+0.17−0.100.47+0.05−0.101.13+0.11−0.210.67+0.60−0.60
Fe/solar (B)1.5+0.3−0.30.8+0.2−0.3>1.11.0+0.1−0.1<0.41.0+0.8−0.3
ξ (B)56+12−324+17−9<14207+5−5<24<14
Norm (B)b0.38+0.21−0.130.24+0.04−0.030.27+0.19−0.200.06+0.01−0.010.30+0.36−0.123.68+1.50−1.99
Fe xxv line (keV)6.65+0.05−0.076.74+0.04−0.046.67+0.03−0.03
EW (eV)17+7−716+6−634+8−8
Flux (10−5 ph cm−2 s−1)0.62+0.26−0.270.44+0.17−0.170.55−0.13−0.13
▵χ281010
Fe xxvi line (keV)6.95+0.03−0.036.88+0.02−0.036.95+0.06−0.076.97+0.07−0.12
EW (eV)29+9−927+7−734+23−2316+9−8
Flux (10−5 ph cm−2 s−1)0.85+0.27−0.270.63+0.17−0.170.14+0.10−0.100.23+0.13−0.12
▵χ21625410
BAT constant1.11+0.11−0.110.87+0.10−0.101.06+0.23−0.230.93+0.08−0.080.68+0.10−0.09
χ2ν770.2/647879.6/832612.0/538811.6/722835.4/808830.4/865
Soft excessArk 120Fairall 9MCG −02-14-009Mrk 335NGC 7469SWIFT J2127.4+5654
××
Γ2.12+0.01−0.022.01+0.01−0.011.90+0.03−0.032.15+0.01−0.021.80+0.01−0.012.20+0.05−0.05
Norma1.13+0.02−0.020.78+0.02−0.020.12+0.05−0.030.57+0.02−0.020.55+0.05−0.071.65+0.05−0.08
kerrconv
q4.1+0.8−0.8>6.15<3.56.6+2.0−1.0>3.02.2+0.6−0.4
a>0.970.98+0.01−0.01<0.960.87+0.05−0.06<0.97
i°48+3−670+5−240+11−1453+6−670+4−343+18−7
Fe/solar (A)1.4*1.9*0.4*2.0+0.8−0.41.6*>0.2
ξ (A)<13<11<543<11<11<26
Norm (A)b1.31+0.14−0.151.13+0.08−0.190.10+0.17−0.100.47+0.05−0.101.13+0.11−0.210.67+0.60−0.60
Fe/solar (B)1.5+0.3−0.30.8+0.2−0.3>1.11.0+0.1−0.1<0.41.0+0.8−0.3
ξ (B)56+12−324+17−9<14207+5−5<24<14
Norm (B)b0.38+0.21−0.130.24+0.04−0.030.27+0.19−0.200.06+0.01−0.010.30+0.36−0.123.68+1.50−1.99
Fe xxv line (keV)6.65+0.05−0.076.74+0.04−0.046.67+0.03−0.03
EW (eV)17+7−716+6−634+8−8
Flux (10−5 ph cm−2 s−1)0.62+0.26−0.270.44+0.17−0.170.55−0.13−0.13
▵χ281010
Fe xxvi line (keV)6.95+0.03−0.036.88+0.02−0.036.95+0.06−0.076.97+0.07−0.12
EW (eV)29+9−927+7−734+23−2316+9−8
Flux (10−5 ph cm−2 s−1)0.85+0.27−0.270.63+0.17−0.170.14+0.10−0.100.23+0.13−0.12
▵χ21625410
BAT constant1.11+0.11−0.110.87+0.10−0.101.06+0.23−0.230.93+0.08−0.080.68+0.10−0.09
χ2ν770.2/647879.6/832612.0/538811.6/722835.4/808830.4/865

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

breflionx normalization in units of 10−5.

In addition, the best fit obtained with model E for Ark 120 requires an inclination of i = 48°, whilst not unreasonable, this differs from the inclinations measured using models C and D. However, freezing the inclination at i = 35° worsens the fit further (χ2ν = 824.6/648 compared to χ2ν = 770.2/647). A similar dual-reflector fit as used here is considered by Nardini et al. (2010) in which the emissivity index is frozen at q = 5.0 compared to q = 4.1+0.9−0.9 measured here. Nardini et al. (2010) also noted that a relatively high inclination of i° = 57+5−12 obtained in their analysis is likely too large for an object such as Ark 120 although it is possible for the accretion disc and obscuring material to be misaligned. One area of disagreement between the results obtained with model E and the Nardini et al. (2010) fit is the determination of the spin parameter. Nardini et al. found 0.24 < a < 0.93 whereas here only a lower limit of a > 0.97 is found; this may be due to the degeneracies between the spin and emissivity index which are both allowed to vary in model E.

Similarly, NGC 7469 is best fitted with a kerrconv inclination parameter of i° = 70+4−3 whereas models C and D suggest that the disc is inclined at i∼ 23° to the observer, worsening the fit by ▵χ2∼ 80. These high inclinations are likely to be driven by the need to model the relatively smooth soft excess.

For the case of MCG −02-14-009, model E provides an approximately equal quality of fit compared to models C and D. No significant amount of blurring is required to model the spectrum, obtaining parameters similar to those found previously, i.e. low emissivity and an unconstrained spin parameter. The narrow Fe xxvi ionized emission line is also present, in line with findings from models C and D. The SWIFT J2127.4+5654 spectrum is fitted very well with model E (χ2ν = 830.4/865). Extreme parameters are not required to fit the data, which is in agreement with results from the previous disc line models. This may be due to the lack of any notable soft excess in these two objects. Since a significantly blurred spectrum simulates an excess at lower energies, a high emissivity and near maximally rotating BH would provide an over-excess at soft energies, inappropriate for objects such as MCG −02-14-009 and SWIFT J2127.4+5654. None the less, Miniutti et al. (2009) found q = 5.3+1.7−1.4 in SWIFT J2127.4+5654 compared to q = 2.2+0.6−0.4 here; this may be due to the high level of Galactic absorption and the additional small amount of intrinsic absorption at the redshift of the source. An increase in the emissivity index could be compensated for by increased absorption which is found in the Miniutti et al. (2009) analysis.

3.6 Model F – Blurred reflector (2)

Model F starts with model B as the base model but instead blurring the single reflionx reflection spectrum with the kerrconv convolution model to model any broad residuals present in the spectra. The Fe K line complex can still be modelled in this way since the reflionx model includes Fe Kα emission and blurring the spectrum emulates the resulting profile from a kerrdisk model. Since the narrow Fe Kα core included within reflionx is now relativistically blurred, a narrow Gaussian of a fixed width of 10 eV was added to ensure that the narrow 6.4 keV is still modelled. This may represent the case where the 6.4 keV line is observed from Compton-thin matter, such as the BLR or NLR. Contrary to model E, the soft excess (where present) is modelled using the comptt soft photon Comptonization model.

Any ionized emission due to Fe xxv and Fe xxvi required in model D is also found to be required in this model. In general a good fit to all objects is obtained, with model F clearly providing a better fit to the data compared to model E for Fairall 9, NGC 7469 and SWIFT J2127.4+5654. The kerrconv parameters are also consistent within error bars with the emissivity index, inclination and spin parameter found using the kerrdisk line profile previously, yielding typically slightly lower emissivity indices.

Note that the best-fitting spin parameter value with model F for Fairall 9 is a = 0.40+0.33−0.40 (quoted at the 75 per cent confidence level and is unconstrained at the 90 per cent confidence level) in agreement with a = 0.44+0.04−0.11 found in model D (at the 90 per cent confidence level) and with a = 0.5+0.1−0.3 found by Schmoll et al. (2009).

Given the independent modelling of the soft excess within this model, it is interesting to note that a = 0.72+0.18−0.17 is obtained for NGC 7469 at the 90 per cent confidence level. This value of the spin parameter is in agreement with that found in model D, although providing a slightly worse constraint.

3.7 Comparison with XMM–Newton results

Analysis of the data for these objects obtained with XMM–Newton (all objects other than SWIFT J2127.4+5654; see Table 2) yields no additional constraints upon the properties of the Fe K region; the reflection component is also difficult to constrain due to the lack of hard X-ray data. As with the Suzaku data, the Fe xxvi narrow ionized emission line is found throughout most objects and the Fe xxv line is relatively uncommon, observed only in Ark 120 and Fairall 9. However, the results obtained with a kerrdisk relativistic line profile are consistent with those obtained with model D in the Suzaku data. Similarly, a good fit is obtained with this model for all of the objects in the sample, with maximal spin ruled out from the fits. The 0.5–10.0 keV spectrum is generally well described by a combination of narrow distant emission and a small broad relativistic component, typically EW < 150 eV.

In particular, the XMM data for MCG −02-14-009 yields very similar results to the Suzaku data, the main difference being that a statistically significant Fe xxvi line is not present in the XMM–Newton data. All measurements obtained with a kerrdisk profile are consistent at the 90 per cent confidence level with the Suzaku data. Both data sets therefore reveal a relatively simple 0.5–10.0 keV spectrum; there is no observed soft excess and it is well described by distant reflection in addition to a moderately broad relativistic line profile from neutral Fe occurring at rin > 2.5 rg (obtained from the application of the laor model) with emissivity q = 2.8+2.0−0.8, as seen in Table 9.

Table 9

Fit parameters from model D to XMM–Newton EPIC-pn data. * denotes a frozen parameter. In some cases, the spin parameter a could not be constrained, denoted by –.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469
Γ1.58+0.01−0.011.89+0.01−0.011.80+0.02−0.021.64+0.06−0.051.62+0.02−0.04
Norma0.56+0.01−0.020.350.01−0.010.16+0.01−0.010.29+0.01−0.010.07+0.01−0.01
Plasma temperature (keV)<7.29.0+55.3−4.35.2+10.3−0.1
Plasma optical depth τ0.6+0.7−0.1<0.7<38.0
Fluxb2.90+0.21−0.212.15+0.18−0.180.014+0.003−0.002
▵χ2250595661
kerrdisk profile
Line E (keV)6.4*6.4*6.4*6.4*6.4*
Eqw (eV)114+10−17148+24−29149+76−73116+22−2257+16−20
q1.8+0.3−0.4<1.82.8+2.0−0.81.9+0.2−0.22.3+0.3−0.4
a<0.87<0.96<0.71<0.99
i°24+7−2<2637+12−443+6−1524+3−4
▵χ29717106212
ξ284+394−53<11<11233+24−131526+302−246
Fe/solar1.6+0.5−0.3<0.20.4+0.1−0.11.8+0.3−0.20.5+0.1−0.1
Normc<0.010.56+0.07−0.080.44+0.01−0.010.01+0.01−0.01<0.01
Fe xxv line (keV)6.69+0.03−0.036.78+0.06−0.05
EW (eV)21+6−640+18−18
Flux (10−5 ph cm−2 s−1)0.91+0.26−0.270.55+0.25−0.25
▵χ21112
Fe xxvi line (keV)7.01+0.03−0.037.09+0.06−0.067.00+0.05−0.067.03+0.03−0.03
EW31+7−749+22−2130+20−2034+10−10
Flux (10−5 ph cm−2 s−1)1.11+0.25−0.250.58+0.26−0.250.12+0.08−0.080.54+0.16−0.16
▵χ25414921
χ2ν1757.3/1663753.1/7611147.8/11511464.4/14111049.4/885
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469
Γ1.58+0.01−0.011.89+0.01−0.011.80+0.02−0.021.64+0.06−0.051.62+0.02−0.04
Norma0.56+0.01−0.020.350.01−0.010.16+0.01−0.010.29+0.01−0.010.07+0.01−0.01
Plasma temperature (keV)<7.29.0+55.3−4.35.2+10.3−0.1
Plasma optical depth τ0.6+0.7−0.1<0.7<38.0
Fluxb2.90+0.21−0.212.15+0.18−0.180.014+0.003−0.002
▵χ2250595661
kerrdisk profile
Line E (keV)6.4*6.4*6.4*6.4*6.4*
Eqw (eV)114+10−17148+24−29149+76−73116+22−2257+16−20
q1.8+0.3−0.4<1.82.8+2.0−0.81.9+0.2−0.22.3+0.3−0.4
a<0.87<0.96<0.71<0.99
i°24+7−2<2637+12−443+6−1524+3−4
▵χ29717106212
ξ284+394−53<11<11233+24−131526+302−246
Fe/solar1.6+0.5−0.3<0.20.4+0.1−0.11.8+0.3−0.20.5+0.1−0.1
Normc<0.010.56+0.07−0.080.44+0.01−0.010.01+0.01−0.01<0.01
Fe xxv line (keV)6.69+0.03−0.036.78+0.06−0.05
EW (eV)21+6−640+18−18
Flux (10−5 ph cm−2 s−1)0.91+0.26−0.270.55+0.25−0.25
▵χ21112
Fe xxvi line (keV)7.01+0.03−0.037.09+0.06−0.067.00+0.05−0.067.03+0.03−0.03
EW31+7−749+22−2130+20−2034+10−10
Flux (10−5 ph cm−2 s−1)1.11+0.25−0.250.58+0.26−0.250.12+0.08−0.080.54+0.16−0.16
▵χ25414921
χ2ν1757.3/1663753.1/7611147.8/11511464.4/14111049.4/885

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

b Flux for comptt quoted over the 0.5–10.0 keV range in units of 10−11 erg cm−2 s−1.

creflionx normalization given in units of 10−5.

Table 9

Fit parameters from model D to XMM–Newton EPIC-pn data. * denotes a frozen parameter. In some cases, the spin parameter a could not be constrained, denoted by –.

Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469
Γ1.58+0.01−0.011.89+0.01−0.011.80+0.02−0.021.64+0.06−0.051.62+0.02−0.04
Norma0.56+0.01−0.020.350.01−0.010.16+0.01−0.010.29+0.01−0.010.07+0.01−0.01
Plasma temperature (keV)<7.29.0+55.3−4.35.2+10.3−0.1
Plasma optical depth τ0.6+0.7−0.1<0.7<38.0
Fluxb2.90+0.21−0.212.15+0.18−0.180.014+0.003−0.002
▵χ2250595661
kerrdisk profile
Line E (keV)6.4*6.4*6.4*6.4*6.4*
Eqw (eV)114+10−17148+24−29149+76−73116+22−2257+16−20
q1.8+0.3−0.4<1.82.8+2.0−0.81.9+0.2−0.22.3+0.3−0.4
a<0.87<0.96<0.71<0.99
i°24+7−2<2637+12−443+6−1524+3−4
▵χ29717106212
ξ284+394−53<11<11233+24−131526+302−246
Fe/solar1.6+0.5−0.3<0.20.4+0.1−0.11.8+0.3−0.20.5+0.1−0.1
Normc<0.010.56+0.07−0.080.44+0.01−0.010.01+0.01−0.01<0.01
Fe xxv line (keV)6.69+0.03−0.036.78+0.06−0.05
EW (eV)21+6−640+18−18
Flux (10−5 ph cm−2 s−1)0.91+0.26−0.270.55+0.25−0.25
▵χ21112
Fe xxvi line (keV)7.01+0.03−0.037.09+0.06−0.067.00+0.05−0.067.03+0.03−0.03
EW31+7−749+22−2130+20−2034+10−10
Flux (10−5 ph cm−2 s−1)1.11+0.25−0.250.58+0.26−0.250.12+0.08−0.080.54+0.16−0.16
▵χ25414921
χ2ν1757.3/1663753.1/7611147.8/11511464.4/14111049.4/885
Ark 120Fairall 9MCG −02-14-009Mrk 335NGC 7469
Γ1.58+0.01−0.011.89+0.01−0.011.80+0.02−0.021.64+0.06−0.051.62+0.02−0.04
Norma0.56+0.01−0.020.350.01−0.010.16+0.01−0.010.29+0.01−0.010.07+0.01−0.01
Plasma temperature (keV)<7.29.0+55.3−4.35.2+10.3−0.1
Plasma optical depth τ0.6+0.7−0.1<0.7<38.0
Fluxb2.90+0.21−0.212.15+0.18−0.180.014+0.003−0.002
▵χ2250595661
kerrdisk profile
Line E (keV)6.4*6.4*6.4*6.4*6.4*
Eqw (eV)114+10−17148+24−29149+76−73116+22−2257+16−20
q1.8+0.3−0.4<1.82.8+2.0−0.81.9+0.2−0.22.3+0.3−0.4
a<0.87<0.96<0.71<0.99
i°24+7−2<2637+12−443+6−1524+3−4
▵χ29717106212
ξ284+394−53<11<11233+24−131526+302−246
Fe/solar1.6+0.5−0.3<0.20.4+0.1−0.11.8+0.3−0.20.5+0.1−0.1
Normc<0.010.56+0.07−0.080.44+0.01−0.010.01+0.01−0.01<0.01
Fe xxv line (keV)6.69+0.03−0.036.78+0.06−0.05
EW (eV)21+6−640+18−18
Flux (10−5 ph cm−2 s−1)0.91+0.26−0.270.55+0.25−0.25
▵χ21112
Fe xxvi line (keV)7.01+0.03−0.037.09+0.06−0.067.00+0.05−0.067.03+0.03−0.03
EW31+7−749+22−2130+20−2034+10−10
Flux (10−5 ph cm−2 s−1)1.11+0.25−0.250.58+0.26−0.250.12+0.08−0.080.54+0.16−0.16
▵χ25414921
χ2ν1757.3/1663753.1/7611147.8/11511464.4/14111049.4/885

a Powerlaw normalization given in units of 10−2 ph keV−1 cm−2 s−1.

b Flux for comptt quoted over the 0.5–10.0 keV range in units of 10−11 erg cm−2 s−1.

creflionx normalization given in units of 10−5.

These same observations of Ark 120 and Fairall 9 have previously been modelled with the blurred reflector approach, consisting of a combination of narrow distant emission and a smeared ionized reflector under the assumption of a Schwarzschild geometry (see Brenneman & Reynolds 2009). Consistent results are found for Fairall 9 for which narrow ionized emission lines are also found; however, the authors were unable to constrain the emissivity index, measured at q < 1.8 here and the spin value requires that emission occurs at rin > 1.8 rg. The results found here for Ark 120 suggest consistent measurements for the emissivity index with Brenneman & Reynolds (2009). In comparison with the Suzaku data for this object, both approaches are consistent with model D. Ionized emission lines are also found in Ark 120 throughout both analyses. The inner radius of emission is generally consistent, occurring at rin < 80 rg from Brenneman & Reynolds (2009) and rin > 2.5 rg as measured here from the spin parameter within model D (Table 9).

4 DISCUSSION

This small sample of six AGN includes typically bare Seyfert 1 galaxies featuring little or no intrinsic absorption. This property of these AGN is important since the spectra of these objects are simpler to model without complicating factors (such as absorption), allowing the observer to draw conclusions about the fundamentals of accretion disc properties and basic features of the Fe K region. These conclusions will therefore be less model dependent upon how the warm absorber is modelled (Turner & Miller 2009). It is important to effectively model the spectra of these AGN and assess the likely origin of various components of the spectrum before we can proceed to draw conclusions about more complicated AGN. Modelling of the broad-band continuum is also essential prior to analysis of the Fe K region; this has been achieved using data from both HXD and BAT hard X-ray detectors onboard Suzaku and Swift, respectively. This allows the spectrum spanning 0.5–100.0 keV to be modelled with better constraints upon distant reflection components modelled by pexrav and reflionx, whilst maintaining effective modelling of the soft excess. This approach furthers that taken by Nandra et al. (2007) in which EPIC-pn data from XMM–Newton are analysed in the 2.5–10.0 keV range.

The presence of narrow ionized emission lines due to Fe xxv and Fe xxvi assessed prior to modelling any broad residuals in the Fe K region has an important effect upon the parameters obtained with broad disc line profiles. In some spectra neglecting to model these lines (if present) can accentuate any apparent broad Fe K line, particularly those resulting from ionized species of Fe. Most commonly occurring in the objects analysed here is the 6.97 keV line from H-like Fe which is observed in all objects, except NGC 7469 and SWIFT J2127.4+5654. This is a surprising result given the rarity of such lines found by Nandra et al. (2007) in which only 2/26 sources showed evidence for these lines (although more common in work by Bianchi et al. 2009). This may be due to higher quality (longer exposure) data obtained with Suzaku allowing these lines to be more easily distinguished from broad residuals at these energies. Similarly, neglecting to include a narrow 6.97 keV emission line where present forces the relativistic line profile to have the emissivity index increased to particularly large values (q > 6) and the inclination is slightly increased to ∼45°–55° whereby the blue wing of the profile is forced to model the narrow excess at ∼6.97 keV.

In agreement with Nandra et al. (2007), however, is the rarity of the 6.7 keV line from He-like Fe which only features significantly here in Mrk 335. Models A and B earlier suggested the possible presence of Fe xxv in Fairall 9 and SWIFT J2127.4+5654; however, it was found that excesses at these energies coincided with the peak in the blue wing of the relativistic line profile. This suggests that high (i.e. calorimeter) resolution data are necessary to determine the presence of this particular emission line, particularly when employing a relativistic line profile, such as that which would be obtained with Astro-H (see Kelley et al. 2010).

The introduction of the kerrdisk model over a model considering emission purely from distant material (i.e. model B) improves the fit to the six objects by an average of ▵χ2∼ 20. This implies that broad residuals are a statistically significant feature in all of these objects (although less so in MCG −02-14-009).

4.1 Typical parameters of the accretion disc

The average inclination of the accretion disc to the observer as inferred by the kerrdisk line profiles in model D is i = 33°± 4°; this is consistent with Nandra et al. (2007) who found i = 38°± 6°. Also from model D we find an average emissivity index of q = 2.3 ± 0.2, much lower than those used in previous work for some of these objects, e.g. Fairall 9 (Schmoll et al. 2009) and SWIFT J2127.4+5654 (Miniutti et al. 2009), particularly when using blurred reflection models to model the whole continuum. The line profile produced within model D also suggests a broad-line profile with an average EW = 119 ± 19 eV consistent with an average EW = 91.3 ± 12.8 found by Nandra et al. (2007) using a Gaussian to model the broad residuals. Whilst here for models C and D the centroid rest energy of the line profile is not allowed to extend below 6.4 keV in the rest frame, the basic parametrization within model A suggests an average broad component line energy of 6.29 ± 0.03 keV compared to E = 6.27 ± 0.07 keV according to Nandra et al. (2007).

The low average emissivity index of the objects in this small sample may be due to a number of factors: independent modelling of the soft excess through a Comptonization of soft photons, modelling of ionized emission lines where present in the data and the assumption that emission extends down to rISCO (the kerrdisk model). In accordance with this, within models E and F facilitating the use of a blurred reflector, the average emissivity index is also dependent upon the way in which the soft excess is modelled. Model F (in which the soft excess is modelled with a comptt component, as in models A–D) suggests an average q = 2.0 ± 0.1 which in general agreement with model D. However, model E (in which purely the blurred reflection component is tasked with modelling the soft excess) suggests an average q = 5.0 ± 0.7 for those objects featuring a soft excess and an average q = 2.0 ± 0.4 for those without, namely MCG −02-14-009 and SWIFT J2127.4+5654. Thus, the high emissivity in model E for some objects is driven by the need to fit a featureless soft excess (e.g. see Fig. 7). A blurred reflection model with high emissivity therefore appears to be ruled out given model E produces a significantly worse fit in all four AGN with a soft excess.

νFν plots of models E and F for Mrk 335 showing the degree to which relativistic blurring of the reflection component is required according to differing interpretations of the soft excess. Note that the significant blurring required in model E to account for the soft excess reduces the accuracy of the fit to the Fe K region features. XIS data are in black, HXD in red and BAT is represented by blue circles.
Figure 7

νFν plots of models E and F for Mrk 335 showing the degree to which relativistic blurring of the reflection component is required according to differing interpretations of the soft excess. Note that the significant blurring required in model E to account for the soft excess reduces the accuracy of the fit to the Fe K region features. XIS data are in black, HXD in red and BAT is represented by blue circles.

Model C alternatively models reflection off a relativistic accretion disc under the assumption of a maximally rotating central BH, without the inner radius of emission fixed at the ISCO. This approach also suggests a similarly relatively small equivalent width of EW = 102 ± 14 eV and an average emissivity q = 3.2 ± 0.4 (for those objects for which the emissivity could be constrained), higher than that derived from model D and likely due to the fact that rin is not fixed at rISCO. The average inclination of the accretion disc is also very similar to model D above with i = 34°± 3°. Given the assumptions made within this model, the inner radius of emission is found to originate at tens of rg rather than <6 rg at an average rin = 39 ± 8 rg. Indeed, in Ark 120, Fairall 9 and MCG −02-14-009, model C and this interpretation of the accretion disc and SMBH provide a marginally better fit. However given the degeneracies between q, a and rin, it would be impossible to reasonably constrain any of these values without assumptions similar to those mentioned above.

4.2 Model-dependent black hole spin

The results for determining the spin of the central BHs of the AGN in this sample suggest that the spin derived is very much dependent upon which interpretation of the Fe K line region is followed. Modelling the soft excess through an independent model such as comptt tends to yield low-to-intermediate spin constraints for the objects in this sample, the exceptions being Ark 120 and MCG −02-14-009 in which only an upper limit could be placed. The employment of this interpretation also yields low-to-intermediate emissivity indices for the accretion discs (i.e. q∼ 2). Modelling the 0.5–100.0 keV spectrum with a blurred reflection component and no other modelling of the soft excess (i.e. model E) gives particularly high values of both emissivity index and spin parameter, typically q≳ 4 and a≳ 0.9, but these fits are statistically ruled out here.

This is, however, only in objects featuring a soft excess. In MCG −02-14-009, which has no obvious excess over a power law at low energies, the derived parameters are very similar throughout models D, E and F. Similarly for SWIFT J2127.4+5654, the parameters obtained with these three models are all consistent and indeed the resulting spin parameter is also consistent with previous findings by Miniutti et al. (2009) although they find that a much higher emissivity index is required (q = 5.3+1.7−1.4 compared to q = 2.2+0.3−0.9 from model E).

As discussed previously, the spin constraint obtained here for Fairall 9 agrees with that found by Schmoll et al. (2009), but only in the case where they ignore the spectrum below 2 keV (a = 0.5+0.1−0.3 compared to a = 0.44+0.04−0.11 found here in model D). According to model D, an intermediate spin of the central BH within NGC 7469 is found, a = 0.69+0.09−0.09. This is also consistent with that found within model F, a = 0.72+0.18−0.17, suggesting that the spin of this object is indeed a∼ 0.7.

These findings therefore suggest that the two differing interpretations give concurring spin constraints only when the blurred reflector is not required to model the soft excess, leading to the conclusion that the component responsible for the soft excess must be independently fitted before conclusive spin constraints can be made. Even so, the origin of the soft excess is not currently known (e.g. see Gierlinski & Done 2004). The relatively constant temperature of the soft excess versus BH mass suggests that it may not arise from direct thermal emission from the accretion disc, since the accretion temperature properties should scale with M−1/4BH in a standard accretion disc. An atomic origin of the soft excess has been suggested; however, no obvious spectral features are seen in high-resolution data meaning that if atomic emission is responsible, it must be significantly relativistically blurred such as here in model E (Fiore, Matt & Nicastro 1997; Ross, Fabian & Ballantyne 2002). Alternatively, it is suggested that the soft excess could arise from relativistically blurred absorption (Gierlinski & Done 2004) from a differentially rotating and outflowing disc wind (Murray & Chiang 1997). The smeared absorption creates a smooth hole in the spectrum resulting in the apparent soft excess and a hardening at high energies. Whilst in this paper we do not aim to determine the origin of the soft excess, we find that a simple parametrization of the soft excess continuum through a model such as comptt provides a better fit to the spectra than an atomic origin from a highly blurred reflection component.

4.3 Reflection off a disc wind?

One possible alternate origin of the observed broad component in AGN spectra is through reflection of primary continuum photons off a Compton thick disc wind, producing emission, absorption and broad features in the Fe K band (Sim et al. 2010). Fast outflows can add significant complexity to the Fe K region and can reproduce a wide range of spectral signatures owing to the many differing physical conditions. Sim et al. (2010) noted that such outflows can, however, have a significant effect upon the soft X-ray spectrum resulting in features such as highly blueshifted absorption lines, although these absorption lines are weaker than the features seen in the Fe K band. The lines from a disc wind may not be observed in objects such as those in this sample if they are observed relatively face-on; however even if we are not looking directly through the line of sight to the wind, it may be possible to observe broadened emission features in reflection. For instance, Sim et al. (2010) predict broadened Fe K emission down to ∼5 keV in their disc wind spectra, with equivalent widths typically several tens of eV.

5 CONCLUSIONS

Resulting from the above work on this sample of six ‘bare’ Seyfert AGN, we conclude the following.

  • Narrow emission from distant material is very important when modelling the Fe K region, and neglecting to include these components where present can have a significant effect upon the accretion disc parameters obtained with subsequent models. The narrow neutral 6.4 keV core is ubiquitous amongst these six AGN.

  • Ionized emission at ∼6.97 keV due to Fe xxvi is more common in these spectra than previously thought (e.g. Bianchi et al. 2004; Nandra et al. 2007) although Bianchi et al. (2009) found them to be more common. It is present here in 4/6 objects in the Suzaku data and it is also present in the same objects during the five XMM–Newton observations included here.

  • Emission at ∼6.70 keV due to Fe xxv is still relatively uncommon amongst these AGN, present in only 1/6 Suzaku and 2/5 XMM–Newton observations. However, it is likely that this emission line may be more common amongst AGN since this is also the energy at which the blue wing from a relativistically blurred neutral 6.4 keV line occurs with typical model parameters (e.g. i = 30°–40°, q = 2–3), making it difficult in most cases to distinguish between these two components in CCD resolution spectra.

  • The fit to all objects in the sample is significantly improved with the introduction of a relativistic line profile, such as that from a kerrdisk model. However, this broad emission is weaker after effective modelling of narrow components (mean EW = 119 ± 19 eV).

  • The average emissivity index of the accretion disc is a low to moderate q = 2.3 ± 0.2 when using comptt to model the soft excess. The emissivity index typically scales as R−2.3, consistent with an R−2 law for the illuminating continuum and therefore implying that strong GR effects (such as light bending) may not be required in these objects. The assumption that emission extends down to rISCO is still valid, however, since the flat emissivity of the disc indicates that emission is not centrally concentrated. This is consistent with emission occurring at typically tens of rg rather than within <6 g.

  • It is essential to effectively model the continuum and the soft excess prior to attempting to constrain parameters of the accretion disc and inner regions of AGN. Independent modelling of the soft excess with models such as comptt or a second soft powerlaw rather than blurring reflection from inner parts of the accretion disc as a way to model the soft excess yields consistent results with respect to AGN with and without an excess at lower energies.

  • Purely using a blurred reflection component to model the soft excess results in typically higher values of spin and emissivity index. This would lead to a skewed distribution amongst AGN towards the higher values of these parameters for those objects with a soft excess, i.e. a > 0.9 and q > 4. This would indicate that the regions of the accretion disc responsible for the soft excess (in an atomic origin of the soft excess) are not the same as those responsible for features in the Fe K region. In this small sample of six AGN, the fit is significantly worse using this interpretation for those objects with a soft excess, namely Ark 120, Fairall 9, Mrk 335 and NGC 7469.

  • New constraints suggest further intermediate spin values for Mrk 335 and NGC 7469 of a = 0.70+0.12−0.01 and 0.69+0.09−0.09, respectively. The spin value for SWIFT J2127.4+5654 found here is consistent with that found by Miniutti et al. (2009), whereas the spin found here for Fairall 9 is consistent with the Schmoll et al. (2009) analysis above 2 keV. Only upper bounds for the spin parameter can be placed for Ark 120 and MCG −02-14-009, requiring that emission does not occur within 2.02 rg and 2.45 rg, respectively.

  • Overall zero spin cannot be ruled out at a particularly high confidence level in all objects, as can be seen in Table 7. None of the objects in the sample show overwhelming evidence for a Kerr geometry and significantly centrally concentrated emission; maximal spin is ruled out in most cases (models D and F).

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). Data obtained with XMM–Newton have also been used within this paper, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA.

REFERENCES

Anders
E.
Grevesse
N.
,
1989
,
Geochimica Cosmochimica Acta
,
53
,
197

Arnaud
K. A.
,
1996
, in
Jacoby
G.
Barnes
J.
, eds, ASP Conf. Ser. Vol. 101,
Astronomical Data Analysis Software and Systems V
.
Astron. Soc. Pac.
, San Francisco, p.
17

Bianchi
S.
Matt
G.
Balestra
I.
Guainazzi
M.
Perola
G. C.
,
2004
,
A&A
,
422
,
65

Bianchi
S.
Guainazzi
M.
Matt
G.
Fonseca
Bonilla N.
Ponit
G.
,
2009
,
A&A
,
495
,
421

Blandford
R. D.
Znajek
R. L.
,
1977
,
MNRAS
,
179
,
433

Blustin
A. J.
et al.,
2007
,
A&A
,
466
,
107

Brenneman
L. W.
Reynolds
C. S.
,
2006
,
ApJ
,
652
,
1028

Brenneman
L. W.
Reynolds
C. S.
,
2009
,
ApJ
,
702
,
1367

Cusumano
G.
et al.,
2010
,
A&A
,
510
,
48

Fabian
A. C.
Rees
M. J.
Stella
L.
White
N. E.
,
1989
,
MNRAS
,
238
,
729

Fiore
F.
Matt
G.
Nicastro
F.
,
1997
,
MNRAS
,
284
,
731

Fukazawa
Y.
et al.,
2009
,
PASJ
,
61
,
17

Gallo
L. C.
Lehmann
I.
Pietsch
W.
Boller
Th.
Brinkmann
W.
Friedrich
P.
Grupe
D.
,
2006
,
MNRAS
,
365
,
688

George
I. M.
Fabian
A. C.
,
1991
,
MNRAS
,
249
,
352

Gierlinski
M.
Done
C.
,
2004
,
MNRAS
,
349
,
L7

Gondoin
P.
Orr
A.
Lumb
D.
Santos-Lleo
M.
,
2002
,
A&A
,
388
,
74

Gruber
D. E.
Matteson
J. L.
Peterson
L. E.
Jung
G. V.
,
1999
,
ApJ
,
520
,
124

Grupe
D.
Komossa
S.
Gallo
L. C.
Fabian
A. C.
Larsson
J.
Pradhan
A. K.
Xu
D.
Miniutti
G.
,
2008
,
ApJ
,
681
,
982

Kalberla
P. M. W.
Burton
W. B.
Hartmann
D.
Arnal
E. M.
Bajaja
E.
Morras
R.
Poppel
W. G. L.
,
2005
,
A&A
,
440
,
775

Kelley
R. L.
et al.,
2010
,
BAAS
,
41
,
738

King
A. R.
Pringle
J. E.
Hofmann
J. A.
,
2008
,
MNRAS
,
385
,
1621

Koyama
K.
et al.,
2007
,
PASJ
,
59
,
23

Laor
A.
,
1991
,
ApJ
,
376
,
90

Larsson
J.
Miniutti
G.
Fabian
A. C.
Miller
J. M.
Reynolds
C. S.
Ponti
G.
,
2008
,
MNRAS
,
384
,
1316

Longinotti
A. L.
Sim
S. A.
Nandra
K.
Cappi
M.
,
2007
,
MNRAS
,
374
,
237

Magdziarz
P.
Zdziarski
A. A.
,
1995
,
MNRAS
,
273
,
837

Malizia
A.
et al.,
2008
,
MNRAS
,
389
,
1360

Matt
G.
,
2002
,
MNRAS
,
337
,
147

Miniutti
G.
Fabian
A. C.
Goyder
R.
Lasenby
A. N.
,
2003
,
MNRAS
,
344
,
L22

Miniutti
G.
et al.,
2007
,
PASJ
,
59
,
315

Miniutti
G.
Panessa
F.
De Rosa
A.
Fabian
A. C.
Malizia
A.
Molina
M.
Miller
J. M.
Vaughan
S.
,
2009
,
MNRAS
,
398
,
255

Mitsuda
K.
et al.,
2007
,
PASJ
,
59
,
1

Morrison
R.
McCammon
D.
,
1983
,
ApJ
,
270
,
119

Murray
N.
Chiang
J.
,
1997
,
ApJ
,
474
,
91

Mushotzky
R. F.
Fabian
A. C.
Iwasawa
K.
Kunieda
H.
Matsuoka
M.
Nandra
K.
Tanaka
Y.
,
1995
,
MNRAS
,
272
,
L9

Nandra
K.
George
I. M.
Mushotzky
R. F.
Turner
T. J.
Yaqoob
T.
,
1997
,
ApJ
,
477
,
602

Nandra
K.
O’Neill
P. M.
George
I. M.
Reeves
J. N.
,
2007
,
MNRAS
,
382
,
194

Nardini
E.
Fabian
A. C.
Reis
R. C.
Walton
D. J.
,
2010
,
MNRAS
, in press (doi:10.1111/j.1365-2966.2010.17518.x)

O’Neill
P. M.
Nandra
K.
Cappi
M.
Longinotti
A. L.
Sim
S. A.
,
2007
,
MNRAS
,
381
,
94

Porquet
D.
,
2006
,
A&A
,
445
,
L5

Porquet
D.
Reeves
J. N.
O’Brien
P.
Brinkmann
W.
,
2004
,
A&A
,
422
,
85

Reeves
J. N.
et al.,
2007
,
PASJ
,
59
,
301

Ross
R. R.
Fabian
A. C.
,
2005
,
MNRAS
,
358
,
211

Ross
R. R.
Fabian
A. C.
Ballantyne
D. R.
,
2002
,
MNRAS
,
336
,
315

Schmoll
S.
et al.,
2009
,
ApJ
,
703
,
2171

Scott
J. E.
et al.,
2005
,
ApJ
,
634
,
193

Sim
S. A.
Miller
L.
Long
K. S.
Turner
T. J.
Reeves
J. N.
,
2010
,
MNRAS
,
404
,
1369

Takahashi
T.
et al.,
2007
,
PASJ
,
59
,
35

Tanaka
Y.
Inoue
H.
Holt
S. S.
,
1994
,
PASJ
,
46
,
L37

Tanaka
Y.
et al.,
1995
,
Nat
,
375
,
659

Thorne
K. S.
,
1974
,
ApJ
,
191
,
507

Titarchuk
L.
,
1994
,
ApJ
,
434
,
313

Tueller
J.
et al.,
2010
,
ApJS
,
186
,
378

Turner
T. J.
Miller
L.
,
2009
,
A&AR
,
17
,
47

Volonteri
M.
Sikora
M.
Lasota
J.-P.
,
2007
,
ApJ
,
667
,
704