A&A 470, 497-503 (2007)
DOI: 10.1051/0004-6361:20067036
M. Krumpe - G. Lamer - A. D. Schwope - B. Husemann
Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
Received 28 December 2006 / Accepted 2 April 2007
Abstract
Aims. We present the analysis and results of a 20 ks XMM-Newton observation of RBS1423.
Methods. X-ray spectral analysis is used to establish a significantly broadened relativistic iron K
line from a highly ionised disk.
Results. A QSO at z=2.262 was considered to be the optical counterpart of this ROSAT Bright Survey X-ray source. Based on the improved XMM-Newton source position we identified a z=0.208 QSO as optical counterpart to RBS1423. The 0.2-12 keV X-ray luminosity of this radio-quiet QSO is
.
The XMM-EPIC spectra are well described by a power law with a significantly broadened iron K
line. Disk line models for both Schwarzschild and Kerr black holes require hydrogen-like iron ions to fit the measured line profile. Significant ionisation of the reflection disk is confirmed by model fits with ionised disk models, resulting in an ionisation parameter
.
Key words: galaxies: active - quasar: general - X-rays: individuals: RBS1423
The X-ray spectra of Seyfert type I galaxies and radio quiet
quasi stellar objects (QSOs)
are generally dominated by a power law component which
is often accompanied by a softer component below 1 keV
("soft excess'', Arnaud et al. 1985).
An iron
fluorescence line at 6.4 keV was first observed
in OSO-8 spectra of Centaurus A (Mushotzky et al. 1978). It was found to be a common
feature in the X-ray spectra of AGN by the GINGA satellite
(Nandra & Pounds 1994).
With the higher spectral resolution of the ASCA satellite, very broad emission features
were discovered in the X-ray spectra of a few nearby Seyfert galaxies
(e.g. MCG-6-30-15, Tanaka et al. 1995; NGC 3783, George et al. 1998).
These features are usually interpreted as iron K
fluorescence
originating from reflection by the innermost parts of the AGN accretion
disk, where they are broadened by gravitational redshift and
relativistic motion.
With the launch of XMM-Newton and CHANDRA higher quality spectra
became available, which in many cases confirmed the relativistic
reflection scenario (e.g. for MCG-6-30-15, Wilms et al. 2001;
NGC 3516, Turner et al. 2002).
While the relativistically broadened iron lines are now
found in large fraction among the local, low luminosity Seyfert galaxies
(Nandra et al. 2006; Guainazzi et al. 2006),
detections in higher luminosity Seyfert galaxies and QSOs e.g., E1821+643
(Jimenéz-Bailón et al. 2007), 3C109 (Miniutti et al. 2006), or Q0056-363
(Porquet & Reeves 2003) are still rare.
The low number of detections in luminous QSOs is partly
due the small number of bright, low redshift, luminous QSOs.
On the other hand, it has been shown that the equivalent width of narrow
iron K
emission lines is anti-correlated with X-ray luminosity
(Iwasawa & Taniguchi 1993). This so-called X-ray Baldwin effect has also been
found for the equivalent widths of relativistically broadened
iron lines (Nandra et al. 1997).
In a sample of 38 PG QSOs and (radio-loud) quasars, only three XMM-Newton spectra
showed significant detections of relativistically broadened iron lines
(Jimenéz-Bailón et al. 2005).
Here we present a 20 ks XMM-Newton observation of the
ROSAT X-ray source RBS1423 (RA = 1444
14
,
Dec = 06
32
30
).
In the ROSAT bright survey (RBS, Schwope et al. 2000)
this source was identified with a QSO at z=2.262.
Based on the XMM-Newton position, we
assign a new optical counterpart to RBS1423,
a QSO with z=0.208 and an 0.2-12 keV X-ray luminosity
.
The paper is organised as follows. In Sects. 2-4 we describe the XMM-Newton data, the optical counterpart identification to RBS1423, and optical observations. In Sect. 5 we analyse the X-ray data by fitting different models. Finally, our conclusions are discussed in Sect. 6.
RBS1423 was observed by XMM-Newton on February 11, 2005
(ObsID 0207130401, orbit 948, exposure time 20.2 ks) with
the European Photon Imaging Cameras (EPIC) as primary instruments.
The EPIC-PN camera was operated in standard Full Window mode with
thin filter, while both EPIC-MOS cameras observed in Full Window mode
with medium filters. The data were processed with SAS version 7.0
(Science Analysis Sofware) package, including the corresponding calibration
files. The epchain and emchain tasks were used for generating
linearised event lists from the raw PN and MOS data.
After cleaning the data of observation time with high background,
we obtained net exposure times of 14.1 ks for the EPIC-PN camera and
16.7 ks for both EPIC-MOS cameras.
The count rates
were
and
per camera, respectively.
Therefore, the effects of photon pile-up are negligible for both cameras.
![]() |
Figure 1: Left: SDSS R-band image of RBS1423 with overplotted XMM-Newton EPIC-PN X-ray contours. The former optical counterpart of RBS1423 (QSO at z=2.262) is clearly not associated with the XMM-Newton source. Instead, the QSO at z=0.208 (marked as Sy1) is the related optical counterpart. The cross indicates the ROSAT X-ray position of RBS1423. The X-ray source in the upper right corner is unrelated to RBS1423. Right: optical spectrum of the QSO at z=0.208 (optical counterpart to RBS1423). The SDSS spectrum (green solid line) shows a significantly increased flux compared to the Keck spectrum (black solid line). XMM-Newton Optical Monitor magnitudes are plotted as triangles, while SDSS u-, g-, r-, i-, z-magitudes are plotted as crosses. The most prominent quasar emission lines are labelled. |
Open with DEXTER |
The RBS-catalogue originally identified a QSO at z=2.262(Schwope et al. 2000)
as X-ray source of RBS1423, but mentioned a possible contribution of
a Seyfert I galaxy at 54 arcsec distance (z=0.208).
Considering the measured 0.5-2.0 keV ROSAT flux
of
erg cm-2 s-1, RBS1423 would have been one of the most
luminous X-ray sources in the sky with
.
Based on the linearised event list, we produced 0.2-10 keV broad band X-ray
images for all EPIC cameras. A combined source detection was performed
(edetect_chain) with a minimum detection likelihood ml=5. The source
detection list was astrometrically corrected for a systematic XMM-Newton
position error (eposcorr). Furthermore, we generated an
astrometrically corrected EPIC-PN X-ray image, which was used to plot the
X-ray contours in Fig. 1. Our XMM-Newton X-ray image
clearly reveals the Seyfert I galaxy at z=0.208 as the only counterpart to
RBS1423. According to the conventional dividing line of
between Seyfert galaxies and QSOs, we classify the
(see Sect. 4.3) object as a low luminosity QSO.
The QSO at redshift z=2.262 was not
detected as an X-ray source above the detection likelihood limit.
The offset between the XMM-Newton and the ROSAT position
(the latter indicated by a cross in Fig. 1) is 41 arcsec.
The angular distance between the z=2.26 QSO and the
z=0.2 QSO is 54 arcsec. Hence, the off-axis
location of the X-ray source is negligible for the further data analysis.
In total, we collected 40 000 X-ray photons from RBS1423 allowing proper
spectral analysis (Sect. 5).
More than four years later, on June 15, 2001, the Sloan Digital Sky Survey
(SDSS, see e.g. Adelman-McCarthy et al. 2006) obtained photometric magnitudes in the u-, g-, r-,
i-,
z-SDSS-bands for the QSO at z=0.2. On May 3, 2005, an SDSS spectrum was
taken with an exposure time of 3600 s. The SDSS spectrum is shown in
Fig. 1, right panel as a green line. The line widths are
FWHM
km s-1, FWHM
km s-1, and
FWHM
km s-1.
Optical and near-ultraviolet images of RBS1423 were obtained with the optical
monitor (OM, Mason et al. 2001) onboard XMM-Newton through U and UVW1 filters with central wavelengths
3440
and 2910
,
respectively. The exposure time in the U-filter was 5000 s
and 3174 s in the UVW1.
The data were processed with the OM reduction pipeline omichain.
The resulting count rates were corrected for detector dead time and for the time-dependent degradation factor, as computed by the SAS task ommag. Zero point and flux conversion factors from the XMM-Newton calibration files allowed us to calculate optical fluxes. The U- and UVW1-fluxes of the OM are marked in Fig. 1 (right panel) as triangles.
We converted the SDSS g- and r-magnitudes to Johnson B-magnitudes
following Smith et al. (2002) also applying a k-correction.
Based on the derived absolute B-band magnitude of
,
we re-classify RBS1423 as a low luminosity QSO.
Comparing the different optical data sets of this QSO
(Fig. 1), it is obvious that
RBS1423 is highly variable in the optical.
The SDSS spectrum shows a fivefold increased continuum flux and
a tenfold increased H-
line flux compared to the Keck spectrum.
The OM photometric data show a pronounced rise of the spectral flux at short
wavelength.
We searched for X-ray variability in our XMM-Newton data, but found no
significant deviation from the mean count rate.
The measured XMM-Newton flux of
erg cm-2 s-1 is comparable to the ROSAT flux.
Multiple observations covering longer time scales are required to observe a possible
X-ray variability.
![]() |
Figure 2:
Line profile analysis of H-![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
PN (black) and MOS (green, red) cameras spectra of RBS1423 (in observer frame).
A single power law has been fitted to the 0.2-12.0 keV data. A broad
and significant positive residual is seen between ![]() ![]() |
Open with DEXTER |
X-ray spectra of RBS1423 were extracted from calibrated photon event lists
with SAS task especget which also
creates the appropriate redistribution matrices and exposure-corrected
effective area files.
Instead of using the published effective area file shipped with SAS 7.0
for EPIC-PN we used the more recent version
XRT3_XAREAEF_0011.CCF kindly made available by Haberl (MPE).
This version implements a slight reduction of the XRT3 effective areas
(with respect to XRT1 and XRT2) in the 6-8 keV range,
which had been measured during the pre-flight calibration of the telescopes.
The new version improves the concordance of the MOS and PN spectra
at the energy of the Fe K line.
X-ray events corresponding to patterns 0-12 and 0-4 events (single and double
pixels) were selected for MOS and PN. For the EPIC-PN camera
only X-ray events with flag 0 were included. The MOS X-ray events were
extracted by using the recommended flag selection XMMEA_SM. A low-energy
cutoff was set to 0.2 keV, while the high-energy cutoff was fixed to 12.0 keV.
Identical source regions were used for all EPIC cameras in order to extract
spectra. All spectra were binned to a minimum of 30 counts per bin to apply
the
minimisation technique in the spectral fitting, which was
perfomed by using the XSPEC package version 12.0.
All quoted errors are 68% limits unless otherwise mentioned.
The galactic hydrogen column density in the line of sight is
(Dickey & Lockman 1990).
A single power law model with galactic absorption was
fitted to the broad-band 0.2-12.0 keV PN and MOS spectra. The best fitting
power law (/d.o.f. = 1051/754 with
)
showed significant positive residuals around
5 keV, typical for iron
K
emission (Fig. 3).
Below 0.5 keV, a discrepancy between PN and both MOS
cameras is noticable (
15% higher count rates in MOS cameras relative to PN).
The PN data are broadly consistent with a power law emission spectrum and
foreground galactic absorption,
while the MOS data indicate a possible soft excess.
Due to this discrepancy we use only the spectra above 0.5 keV for the analysis
of the possible Fe-K
feature.
Moreover, we fit the PN and MOS data separately.
Table 1: Results of the spectral analysis for the PN detector (0.5-12.0 keV). Notation of the models is according to XSPEC. Frozen parameters are between exclamation marks. Energies and equivalent width are given in the quasar frame.
Table 2: Results of the spectral analysis for the MOS detectors (0.5-12.0 keV). Notation of the models is according to XSPEC. Frozen parameters are between exclamation marks. Energies and equivalent width are given in the quasar frame.
All fit results are summarised in Table 1 for the PN detector and in Table 2 for the MOS cameras. Line energies and equivalent widths presented in these tables are given in the source frame.
The line feature is also apparent, when we fit single power laws to the PN and MOS data separately in the energy range 0.5-12.0 keV. We then tried to fit the spectra by introducing a second power law (labelled as wabs(power1 + power2)). However, the remaining excess still indicates the presence of a broad emission line feature.
As the next simple approach, we tried to fit the data using
a single power law with a Gaussian line superposed.
The line was detected at
keV with a
) keV
in the PN detector and at
keV (
keV) in the
MOS detectors. This corresponds to a rest-frame energy of 5.6 keV and 6.3 keV,
respectively. Considering the statistical errors, the line energies
measured by the MOS and PN cameras differ only marginally. However, some
systematic influence of the effective area calibrations on the line energies
cannot be ruled out.
The significant width of the line suggests
relativistic broadening of the line.
We therefore fitted iron K disk line
models to the spectra: The XSPEC model diskline (Fabian et al. 1989)
calculates the line broadening in an accretion disk around a non-rotating
(Schwarzschild) black hole (BH)
and the model laor (Laor 1991) around a spinning (Kerr) BH.
The line emissivity is
parameterised by a power law as a function of radius, R-q with q
as a free parameter.
Since q is unconstrained we used an accretion disk
emissivity law
for the Schwarzschild BH model
and q=3 for the Kerr BH model (R-3).
The relevant model parameters and the equivalent widths of the emission line
fits are given in Tables 1 and 2.
For both cameras the laor fit results are statistically the better
fits, favouring the rotating black hole scenario.
Figure 4 shows that for the laor model
the emission line is required with high significance and that it
is significantly broadened in the relativistic regime.
The rest-frame line energy for a simultaneous fit of the PN and MOS data is
(
) keV.
Even if a broken power law plus a relativistic disk model laor
is used, the broad-line feature is dectected with a significance of more
than 3
in the PN and MOS data.
In none of the cases is the source frame energy of the line consistent with
the K
energy expected for neutral iron (6.4 keV).
Instead, the fits suggest that the iron feature is due to fluorescence
from hydrogen-like iron with a K
energy of 6.9 keV.
Fixing the source frame line energy to 6.4 keV and fitting the
accretion disk inclination angle in the model laor as a
free parameter formally leads to acceptable
values, but leaves
visible residuals in the range beyond 6.4 keV.
![]() |
Figure 4: Confidence contours (68%, 90%, 99%) of the parameters Ri and line normalisation for a single line with laor profile. Both PN and MOS spectra have been used to calculate the contours. |
Open with DEXTER |
Since broad iron K
features in AGN are connected with
reflection by the accretion disk, other effects modifying
the reflected X-ray spectrum have to be taken into account.
We therefore applied the self-consistent ionised disk reflection models
by Ross & Fabian (2005) to model the X-ray spectrum of RBS1423.
Their calculations include a range of relevant fluorescence lines, absorption
edges, Compton-reflection, and Compton-broadening of the atomic features.
The models are available in tabulated form for fitting with XSPEC (reflion)
and can be convolved with relativistic blurring using the code by
Laor 1991 (model kdblur).
The observed X-ray spectrum is then interpreted as superposition of the illuminated
power law and the reflected component, described by the XSPEC models
wabs(kdblur(power+atabreflion.mod
)).
Since the data quality is not sufficient to constrain all parameters, we set the
outer radius of the relativistically blurring model to 400 Schwarzschild
radii and the inclination to 30 degrees. This choice is not critical to our results.
Furthermore, the power law of the
reflected component is chosen to be equal to the power law of the illuminating
radiation. Again, we fitted the PN and MOS data separately. Within the
calculated errors, the fit parameters agreed very well for both data sets.
The values of the reflion fits (see Fig. 5)
are slightly better
than those of the laor and diskline fits.
The iron abundance was found to be nearly solar (Fe
,
Fe
).
This value seems to be contradictory to the large equivalent widths
of the iron line resulting from the fits including a power law and
a gaussian or relativistically blurred line.
Indeed, if we modify the best fit reflion model by setting the
iron abundance to the lowest possible value (A=0.11) and replace the
iron emission feature by a Gaussian line, we get an estimate of 170 eV for the
equivalent width of the iron line.
It is likely that the iron line fluxes from the power law
based fits are overestimated, since these models do not account for
the effect of the iron K absorption edge and as a result the slope of the power law
becomes too steep.
This explanation is confirmed by the fact that the iron
equivalent width is largest for the fits with the steepest
power law models (see Tables 1 and 2).
Therefore, we consider the results of the physically consistent reflion
model as more reliable.
The best fitting ionisation parameters are
and
.
This result confirms that the reflecting parts of the accretion disk in
RBS1423 are highly ionised. Figure 6 shows that the data
cannot be fitted with a lower ionization parameter value and a different
accretion disk inclination angle. Independent of the inclination angle, a
highly ionised disk is needed.
Figure 7 shows the E2 f(E)-model spectrum for the reflion fit to the PN data. The best fit model is reflection-dominated with a very weak direct power law component. However, the errors in the normalisations of both the direct and reflected components are rather large.
The 0.2-12.0 keV X-ray flux of RBS1423 according to the model is
erg cm-2 s-1. The
rest-frame X-ray luminosity (0.2-12.0 keV) is
.
This makes RBS1423 one of the few X-ray luminous AGN with a detected
broad iron fluorescence line.
Assuming
is representing the
bolometric luminosity, the accretion rate
can be calculated by using
![]() |
(1) |
![]() |
Figure 5:
PN spectrum of RBS1423 (observer frame).
A relativistically blurred self-consitent ionised disk
reflection model ( kdblur, reflion) has
been fitted to the 0.5-12.0 keV data. The comparison with
Fig. 3 shows that the broad positive deviation in
the ![]() |
Open with DEXTER |
![]() |
Figure 6: PN confidence level contours (68, 90, and 99 per cent) for the relativistically blurred disk reflection. The ionization parameter is plotted in dependence on the accretion disk inclination angle. |
Open with DEXTER |
![]() |
Figure 7: E2 f(E) - model spectrum for the RBS1423 PN data with a wabs(kdblur(power + reflion))-fit. The lower dotted line represents the illuminating power law radiation, while the upper dotted line shows the reflected component by the accretion disk. Both components were gravitationally blurred and absorbed by the galactic hydrogen column density. The solid line indicates the superposition of both components. Note that the normalisations of the direct and reflected components are highly degenerate. |
Open with DEXTER |
As an alternative to the disk refelection models we fitted a partial covering model
(XSPEC model zpcfabs) to the
EPIC spectra (see Tables 1 and 2) .
The best fit absorbing column density is
cm-2 with a covering fraction of
40% (PN).
With
values of
/d.o.f. = 420/376 (PN) and
/d.o.f. = 293/278 (MOS) the partial covering model fits the data equally as
well as the ionised reflection model.
Since the X-ray spectrum of RBS1423 does not show a strong
soft excess, models with partial coverage by an ionised absorber do not fit the data
well. For the PN spectrum the 90% upper limit for the ionisation parameter is
.
We analysed an XMM-Newton observation (PN net time exposure
14.1 ks) of RBS1423, a suspected high-luminosity QSO.
We showed that the corresponding optical
counterpart is a broad-line AGN at z=0.208, thus correcting the
identification in Schwope et al. (2000).
The absolute magnitude in the B-band of
and the X-ray
luminosity of
(0.5-12 keV) classify the objects as a QSO.
Several optical observations clearly indicate high optical
variability and a strong UV excess.
In the X-ray spectrum of RBS1423 we detected
a broad emission feature around 5 keV,
which is well described by iron K
fluorescence
in the innermost parts of the accretion disk.
This makes RBS1423 one of the most luminous radio-quiet AGN
with a detection of relativistic disk reflection.
Alternatively, the X-ray spectrum can be modelled equally well
with a partially covered power law model.
The required column density is
with a covering fraction of
.
Since the shape of a highly absorbed component is similar to the shape of
a relativistically broadened iron line, even much better signal-to-noise
EPIC spectra of Seyfert galaxies can be explained by both reflection
and partial covering models.
Fabian et al. (2004) successfully fitted the variable XMM-Newton
spectra of the Seyfert galaxy 1H 0707-495
with ionised reflection models, while Gallo et al. (2004) modelled the
same data with partial absorption.
The spectrum of RBS1423 does not show a strong soft excess. Therefore, the partial absorber must be neutral to avoid a strong contribution to the spectrum below 1 keV.
For a rather luminous AGN like RBS1423 a two-component model, where one component is not absorbed at all and the second is absorbed by a dense column of neutral gas, is hardly feasible (see also the discussion in Miniutti et al. (2006) on the possibility of a partial absorber in 3C 109). Therefore, we clearly favour disk reflection models to interpret the RBS1423 spectrum.
The energy of the line resulting from disk line fits is not
consistent with K
from neutral iron, but suggests
the presence of highly ionised iron.
This is confirmed by ionised disk model fits (reflion, Ross & Fabian 2005)
resulting in an ionisation parameter
.
We compared the self-consistent ionised disk model fits with models
including a power law and a broadened emission line.
We conclude that
the large equivalent line widths of the power law models are
likely to be overestimated, since the underlying power law tends to be
modelled too steeply when the iron edge is not taken into account.
Our detection adds a new piece of information to the ongoing debate about whether
the accretion disks in high luminosity AGN
are generally more highly ionised than in low
luminosity objects. Nandra et al. (1997) find evidence for a change in the
iron K line profile with increasing luminosity
and also find an anti-correlation of broad-line equivalent width
and luminosity (X-ray Baldwin effect).
They attribute this correlation to the fact that ionisation
increases with luminosity.
Interestingly, the line profile in RBS1423 is similar to the
average profile of the AGN in the
range in Nandra et al. (1997), nearly symmetric with a strong blue wing.
On the other hand, no evidence for strong ionisation has been found in the disk reflection spectra of the luminous quasars PG 1425+267 (Miniutti & Fabian 2006) and 3C 109 (Miniutti et al. 2006). However, these quasars are radio-loud and therefore not fully comparable with RBS1423.
Porquet & Reeves 2003 found significant broad K
reflection from neutral
iron in the radio-quiet QSO Q0056-363, which is about twice as X-ray luminous
as RBS1423. They conclude that the accretion rate of the object has to be less
than
of the Eddington accretion rate to avoid ionisation of the disk. This is at odds
with other estimates of the accretion rate for this object.
In another high luminosity radio-quiet QSO, E1821+643,
a double peaked Fe K
profile was detected in the EPIC spectrum (Jimenéz-Bailón et al. 2007).
The first peak is consistent with neutral iron, the second with FE XXV or
Fe XXVI.
However, the ionised iron line in the object could also be due to contamination
by emission from a surrounding cluster of galaxies.
In summary, we strongly favour the relativistic disk line reflection over the partial covering model. The physically self-consistent ionised disk reflection model (Ross & Fabian 2005) yields a robust detection of disk ionisation. Therefore, we conclude that RBS1423 is one of the first radio-quiet QSO for which ionised disk reflection is established.
Acknowledgements
Mirko Krumpe is supported by the Deutsches Zentrum für Luft- und Raumfahrt (DLR) GmbH under contract No. FKZ 50 OR 0404. Georg Lamer acknowledges support by the Deutsches Zentrum für Luft- und Raumfahrt (DLR) GmbH under contract no. FKZ 50 OX 0201.