Issue |
A&A
Volume 498, Number 1, April IV 2009
|
|
---|---|---|
Page(s) | 67 - 81 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/200811282 | |
Published online | 24 February 2009 |
Revealing X-ray obscured quasars in SWIRE sources with extreme mid-IR/optical flux ratios
G. Lanzuisi1,2 - E. Piconcelli2 - F. Fiore2 - C. Feruglio3 - C. Vignali4,5 - M. Salvato6 - C. Gruppioni5
1 - Dipartimento di Fisica, Università di Roma La Sapienza, P.le A. Moro 2, 00185 Roma, Italy
2 -
Osservatorio Astronomico di Roma (INAF), via Frascati 33, 00040
Monteporzio Catone, Italy
3 -
CEA, Irfu, Service d'Astrophysique, Centre de Saclay, 91191 Gif-sur-Yvette, France
4 -
Dipartimento di Astronomia, Università degli Studi di Bologna, via Ranzani 1, 40127 Bologna, Italy
5 -
INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
6 -
California Institute of Technology, MC 105-24, 1200 East California Boulevard, Pasadena, CA 91125, USA
Received 3 November 2008 / Accepted 6 February 2009
Abstract
Recent works have suggested that selection criteria based on mid-IR properties, i.e.
extreme colors and bright flux levels, can
be used to reveal a population of dust-enshrouded, extremely-luminous quasars at
-2.
However, the X-ray spectral properties of these intriguing objects still remain
largely unexplored.
We have performed an X-ray study of a large sample of bright mid-IR
(
mJy) galaxies showing an extreme MIR/Optical flux
ratio (
/
)
in order to confirm the presence of a luminous active
nucleus in these very red objects.
Sampling of a large area is required to pick up objects at the
highest luminosities given their low surface density.
Accordingly, we have applied our selection criteria to an area of
6 deg2 covered by XMM-Newton/Chandra observations within the
50 deg2 SWIRE
survey, resulting in a final sample of 44 objects.
The vast majority of the source redshifts, both spectroscopic and photometric, are in the range
.
The X-ray coverage of the sample is highly inhomogeneous (from
snap-shot 5 ks Chandra observations to medium-deep XMM-Newton exposures of 70 ks)
and, consequently, a sizable fraction of them (
43%)
remains undetected in the 0.5-10 keV band.
Using spectral or hardness information we were able to estimate the value of
the absorbing column density in 23 sources. 95% of
them are consistent with being obscured by neutral gas with an intrinsic column
density of
22 2. Remarkably, we also find that
55% of these
sources can be classified as type 2 quasars on the basis of their absorption
properties and X-ray luminosity.
Moreover, most of the X-ray undetected sources show extreme mid-IR colors,
consistent with being luminous AGN-powered objects, suggesting they might host
heavily obscured (possibly Compton-thick) quasars in X-rays.
This demonstrates that our selection criteria applied to a wide area survey is very efficient
in finding a large number of type 2 quasars at
.
The existence of this class
of very powerful, obscured quasars at high z could have important
implications in the context of the formation and cosmological evolution of
accreting supermassive black holes and their host galaxies.
Key words: galaxies: active - galaxies: high-redshift - galaxies: nuclei - infrared: galaxies - X-rays: galaxies
1 Introduction
The similarity in the anti-hierarchical
evolution of both star formation rate and AGN activity at z<2 (La
Franca et al. 2005; Arnouts et al. 2007) reveals that the assembly
and evolution of galaxies and those of the super massive black holes (SMBHs) at their centers are
intimately connected (e.g. Page et al. 2004; Alexander et al. 2005; Hopkins et al. 2006; Best et al. 2005; Croton et al. 2006).
However, the census of accreting SMBHs through different cosmic epochs is still
highly incomplete and represents one of the major challenges in modern
extragalactic astronomy.
Clear-cut evidence of the existence of a large number of missing AGNs arises by the fact
that the sources
detected in the recent Chandra/XMM-Newton surveys (e.g. Brandt & Hasinger 2005)
are able to resolve only the 50-60% of the cosmic X-ray background (CXB) in the 5-10 keV band (Worsley et al. 2005),
and their integrated radiation significantly fails to account for the CXB emission peak at 30 keV (Gilli et al. 2007; Comastri 2004). A population of heavily obscured AGNs showing a spectral shape consistent with the unresolved CXB component is usually invoked to fill this gap.
Intriguingly, as stressed by Fiore et al. (2003), about 20-30% of
the sources from Chandra/XMM-Newton surveys are optically-faint
(i.e. ), and with a very high X-ray-to-optical flux ratio
(F2-10/
). Optical spectroscopy of these extreme F2-10/
objects
is therefore extremely challenging even with the largest telescopes.
Results from deep exposures and accurate multiwavelength studies
indicate that the vast majority of these sources could be obscured quasars at
high-z (i.e. z>1, e.g. Mainieri et al. 2005; Maiolino et al. 2006). We have therefore started to uncover the long-sought
missing AGN population of the X-ray absorbed, i.e. type 2, quasars (QSO2s hereafter), that
is a key ingredient for both the CXB synthesis models as well as the
prediction of the SMBH/galaxy formation and evolution models
(e.g. Silk & Rees 1998; Fabian 1999; Granato et al. 2004; Hopkins et al. 2006). In particular, it has been hypothesized that before shining
in the optical as bright point-like blue objects, quasars undergo an
early dust enshrouded phase being associated with the rapid
SMBH growth. This obscured phase is likely
triggered by multiple galaxy encounters and mergers that favor the
infall and the concentration in the nuclear region of large reservoirs
of gas feeding the SMBH as well as the assembly of the massive
spheroids at high z. The link between the growth of the SMBHs and
their host galaxies is believed to be strong as witnessed by the tight
correlations between SMBH masses and bulge luminosity, mass and
velocity dispersion (Merritt & Ferrarese 2001). A full understanding
of AGN feedback processes is therefore a key question to be addressed
in the context of galaxy evolution (e.g. Di Matteo et al. 2005; Alexander et al. 2005; Menci et al. 2008; Somerville et al. 2008).
In the latest years there have been many efforts aimed at unveiling
the missing population of heavily obscured quasars at high zusing alternative approaches with respect to the pure X-ray
selection.
Multiwavelength selection criteria based on the IR band seem to be very promising,
since absorbed AGN optical/UV/X-ray light should be isotropically
re-emitted by the obscuring material at these frequencies.
In particular, multiple independent lines of evidence point towards selection
criteria using mid-IR (3-30 m; MIR hereafter) photometry or combination of MIR and multiwavelenght data, to efficiently collect large samples
of optically-faint quasars candidates at
(e.g. Houck et al. 2005; Yan et al. 2005; Weedman et al. 2006; Polletta et al. 2006; Alonso-Herrero et al. 2006; Daddi et al. 2007; Hickox et al. 2007; Donley et al. 2008, and references therein).
Combining MIR and radio data, Martinez-Sansigre et al. (2005) claimed the
discovery of a distant, optically-faint QSO2 population whose number
density is comparable with that observed for classical optically
bright quasars. This might explain the paucity of obscured objects in
the high luminosity range of the X-ray selected AGN population that
has been reported in many works (Ueda et al. 2003; Steffen et al. 2003; La Franca et al. 2005; but see also Dwelly & Page 2006).
Furthermore, recent studies have provided strong evidence that a large fraction of the dust-obscured galaxies (DOGs hereafter) with extreme
/
ratios (
1000) and bright at 24
m (i.e.
-1 mJy), harbor
a high-z obscured active nucleus (e.g. Dey et al. 2008; Polletta et al. 2008a; Georgantopoluos et al. 2008; Fiore et al. 2008a,b).
However, an accurate determination of the X-ray column density of the absorber in these sources is hampered by the faintness of these sources in the X-ray band (i.e. Polletta et al. 2006; Martinez-Sansigre et al. 2007). Therefore, most of the X-ray information has been derived by the stacking of Chandra/XMM-Newton data or, for small number of sources detected with only a handful of counts, by the hardness-ratio analysis.
![]() |
Figure 1:
a) Left panel: 24 |
Open with DEXTER |
Evidence of X-ray spectroscopically-confirmed heavily-absorbed QSO2s at cosmological distance still remains ambiguous (i.e. Norman et al. 2002; Alexander et al. 2008; Erlund et al. 2008). Indeed, the Compton-thick (i.e. with
2) nature of the best Compton-thick quasar candidate,
IRAS 09104+4109, has recently been questioned (and not confirmed) by Piconcelli et al. (2007a) using XMM-Newton data.
One of the reasons of this difficulty is that building up large and complete samples of bright (i.e. with a 2-10 keV flux
erg cm-2 s-1) QSO2s by homogeneous selection criteria is complex and time-consuming,
because of the large area that must be covered with deep multi-band observations. This is necessary due to the low space density of the AGNs belonging to the bright end of the X-ray luminosity function.
The Spitzer Wide-area InfraRed Extragalactic (SWIRE) survey (Lonsdale et al. 2003) has imaged nearly 50 square degrees in 6 fields: ELAIS-S1, ELAIS-N1, ELAIS-N2, Lockman Hole, XMM-LSS and CDFS with medium-deep MIPS and IRAC photometry, providing a unique opportunity to build up a large sample of luminous QSO2s.
In this paper we report results on an X-ray spectroscopic study of a sample of SWIRE
sources selected on the basis of their extreme
mid-IR/optical flux ratios (
/
)
and 24
m flux density (i.e.
mJy).
These selection cuts have been applied in order to preferentially select sources harboring a luminous, dust-enshrouded AGN (e.g. Polletta et al. 2008a; Dey et al. 2008). Our aim is to confirm the presence of the active nucleus in these extreme DOGs and to directly measure the column density of the absorber via the characterization of their X-ray spectral properties.
The optical magnitudes used throughout the paper are Vega magnitudes.
The cosmological parameters used are H0=70 km s-1 Mpc-1,
and q0 = 0.
Table 1: Optical-IR parameters and redshift of the sample.
2 The sample: selection criteria and datasets
2.1 Optical and IR observations
Infrared and optical photometric data are available from the public archive of the SWIRE survey (Lonsdale et al. 2003, and references therein).
The catalog contains medium-deep photometric data in MIPS 24 m band down to 400
Jy, 5
,
and in the four IRAC channels from 3.6 to 8.0
m (4.1
Jy at 3.6
m, 5
,
Surace et al. 2004).
Optical photometry is also available in different bands for selected areas in each of the 5 SWIRE fields used in this work (ELAIS-S1, ELAIS-N1, ELAIS-N2, Lockman Hole, XMM-Newton SXDS). We did not consider the CDFS here since these data were already studied by Fiore et al. (2008a) using a
/
selection technique similar to our approach.
ELAIS-N1 and N2 have optical photometric information in 5 bands (
U'g'r'i'Z') from the Wide Field survey (WFS, McMahon et al. 2001). The
Lockman Hole sky region is covered with 4-band Ug'r'i' photometry from the SWIRE photometry program.
XMM-LSS has very deep (
RAB<27.7) 5 band BVRi'z' photometry in 1.12 deg2 from the Subaru XMM Deep Survey (SXDS, e.g. Sekiguchi et al. 2005; Furusawa et al. 2008).
Full details on the optical and IR photometry for the ELAIS-S1 field can be found in Feruglio et al. (2008).
We first selected sources showing a 24 m flux density F(24
m
Jy and
/
in order to have an initial sample of DOGs (560 sources e.g. Fig. 1a).
Several sources in ELAIS-S1, ELAIS-N1, ELAIS-N2 and Lockman Hole do not show an R magnitude value in the catalogs mostly since they are faint objects, i.e. with a magnitude above the detection limit. For these sources we performed aperture photometry at the position of the MIPS source.
The photometry was performed using the IRAF apphot tool, with an 1.2 arcsec aperture radius (consistent with the released optical aperture photometry), and a 10 arcsec wide annulus for background estimate. For undetected sources we computed a 3
limiting magnitude using the median value of
at the positions of the source and the relation
![]() |
(1) |
where A is the area corresponding to the aperture.
2.2 Selection criteria
We selected all the sources with available X-ray data from the initial sample of DOGs (i.e. F(24
Jy and
/
sources). The portion of sky covered by XMM-Newton and/or Chandra observations (for a total of
6 deg2) significantly changes across the five SWIRE fields discussed here, varying from
100% for the XMM-SXDS field to
5% for the ELAIS-N2 field. Nonetheless, since such a selection is completely random and not based on specific target observations, it should not introduce any selection effect in the final sample.
We started selecting sources with
/
at the highest
flux density levels since bright 24
m objects with large
/
are predominantly MIR power-law (i.e. AGN-dominated) sources (e.g. Dey et al. 2008; Donley et al. 2008).
We then gradually decreased the
flux density threshold until a statistically significant sample of objects was collected.
In fact, fainter sources (
mJy) selected with the same extreme
/
colors, turn out to be dominated by star-formation activity both in the X-ray and MIR bands (Donley et al. 2008; Pope et al. 2008a).
Accordingly we selected all the sources showing a
/
and a 24
m flux density F(24
m)
1.3 mJy (see Fig. 1a) obtaining in this way a sample of 50 objects, referred as Extreme DOGs (i.e. EDOGs) hereafter.
Since in four cases the same optical object has a couple of IRAC-MIPS counterparts, we retained as the most likely counterpart the MIR source located at the smallest distance from the optical position.
In this way we got a sample of 46 EDOGs within the 6 deg2 of SWIRE survey covered by X-ray observations. This implies a number density of 7.7 deg-2 that is in very good agreement with that obtained applying the same selection criteria to the sample of DOGs in the 8.6 deg2 of the NOAO Deep Wide-Field Survey
field analyzed in Dey et al. (2008) (i.e. 7.6 sources deg-2). The fraction of sources with
/
among all the sources with
m)
1.3 mJy is
10%.
Two sources (namely SWIRE2 J160734.22+544217.3 and SWIRE2 J160945.43+534944.6) fall in a CCD gap in the X-ray image and therefore were excluded from the sample since there is no useful X-ray information for them.
The final sample therefore consists of 44 sources and is listed in Table 1.
Figure 1b shows the MIR color-color diagram (e.g. Lacy et al. 2004) for the sources with
mJy in the 5 SWIRE fields considered in this work (SWIRE
sample hereafter) obtained using the flux densities measured in the four Spitzer-IRAC channels
.
Two distinct populations can be readily separated in this plot. One has blue colors in F5.8/F3.6 and red colors in F8.0/F4.5, most likely consisting of a mixture of starlight-dominated galaxies and
star-forming galaxies with MIR 5-12
m spectra dominated by PAH emission features at low-redshift (
-0.4).
The other population has redder colors and can be mainly identified with galaxies with a MIR power-law spectrum that is a typical signature of an AGN-powered emission.
The thick lines mark the region of the diagram typically occupied by AGNs (Lacy et al. 2004; Alonso-Herrero et al. 2006; Polletta et al. 2006; Lacy et al. 2007; Barmby et al. 2008; see also Donley et al. 2008, for a detailed discussion about the level of contamination by star-forming galaxies in this region).
It is worth noting that all but one sources in our sample lie within the AGN region, and most of them show extreme red colors i.e.
and
similarly to the luminous obscured quasars at z=1.3-3 in the Polletta et al. (2008a) sample. The only source falling outside the box has the lowest spectroscopic redshift in the sample, i.e. source No. 14 at z=0.22, and therefore it is likely that the contribution from the host galaxy stellar light can explain the low value of the F5.8/F3.6 ratio. The range of MIR colors observed for our EDOGs can be likely ascribed to the relative strength of the contribution of star formation to the total MIR luminosity (e.g. Dey et al. 2008; Sajina et al. 2005).
Interestingly, Donley et al. (2007) and Cardamone et al. (2008) found that the most luminous X-ray sources with
erg s-1 in their large samples of MIR sources fall in this box of the IRAC color-color diagram.
2.3 Redshifts
![]() |
Figure 2: a) Left panel: near IR spectrum of source No. 1 at z=1.007. b) Right panel: near IR spectrum of source No. 6 at z=0.914. Both spectra were obtained by MOIRCS/SUBARU observations. |
Open with DEXTER |
Spectroscopic redshifts are available for 7 objects in the sample (0.22
2.54; e.g. Table 1).
In particular, two of them are obtained from new observations performed with the Subaru telescope (see Sect. 2.3.1 for details of data reduction and analysis).
For the remaining 37 sources we were able to derive photometric redshifts thanks to the available optical and IR photometric data. We used the code described in Fontana et al. (2000) and Fontana (2001).
Photometric redshifts were calculated using a
minimization technique on spectral libraries which include starburst, passive galaxy and AGN semi-empirical template spectra from Polletta et al. (2007), Fiore et al. (2008a) and Pozzi at al. (2007).
All the available photometric data points, from the optical to the MIPS 24
m band, were used in the computation (e.g. Feruglio et al. 2008).
The photo-z are reported in Table 1 with errors corresponding to a confidence level of 1.6
and their distribution spans in the range 0.54
3.77.
In four cases the
minimization technique was not able to properly constrain the z value. For these sources (namely Nos. 16, 18, 20 and 44) we decided to assume a fixed value of z=1 which is close to the mean value of the z distribution in our sample (
1.2). It is worth noting that this is a conservative assumption since the typical value of z found for large samples of DOGs is
2 (e.g. Polletta et al. 2008a; Dey et al. 2008, and references therein).
For 19 of 37 sources photometric redshifts from the SWIRE Photometric Redshift Catalogue published in Rowan-Robinson et al. (2008) are available.
Their approach is based on galaxy and quasar templates applied to data at 0.36-4.5
m, and on a set of 4 IR emission templates fitted to IR excess data at 3.6-170
m.
The photo-z distribution from the SWIRE Photometric Redshift Catalogue spans in the range 0.754
2.236.
![]() |
Figure 3:
a) Left panel: |
Open with DEXTER |
2.3.1 Spectroscopic redshifts
We obtained the spectroscopic redshifts for two sources (namely source Nos. 1 and 6) via near-IR spectroscopy.
These observations were performed with MOIRCS (Multi-Object
InfraRed Camera and Spectrograph) mounted on the Subaru telescope
(Hawaii) on 2007 Dec. 22, as backup program (program No. S07B0087N, P.I.
M. Salvato).
We used the zJ500 grism (optimized for the observation of 0.9-1.78
m region) with a R500 filter blocking the contamination of the
second-order spectra at
1.5
m.
Each observation consists of 2 frames dithered along the slit by 5 arcsec.
After standard reduction the spectra have been subtracted by each other in order to improve sky subtraction and flat-fielding and then combined. SWIRE2 J021559.46-045517.5 has been observed for a total of 1200 s and SWIRE2 J021749.00-052306.9 for a total of 1800 s. Sky lines have been used for the wavelength calibration. Results can be summarized as follows:
- SWIRE2 J021559.46-045517.5 (No. 1):
we obtained a redshift of z=1.007 on the basis of the identification of H
( FWHM < 1100 km s-1) [NII]
and [NII]
narrow emission lines (see Fig. 2a).
- SWIRE2 J021749.00-052306.9 (No. 6):
we obtained a redshift of z=0.914 on the basis of the identification of strong [OIII]
, H
and [NII]
narrow emission lines (see Fig. 2b). This value does not match with the tentative redshift of z=0.987 proposed by Lacy et al. (2007), using very low S/N optical spectral data. The FWHM of the H
line is FWHM<1000 km s-1. This support a narrow line quasar classification for this EDOG.
Table 2: X-ray coverage and 0.5-10 keV band flux of the sample.
2.4 MIR properties of EDOGs
We computed 5.8 m monochromatic luminosities of our sources using a linear interpolation between the observed 8
m and 24
m flux densities at the wavelength corresponding to 5.8
m at the source rest-frame (e.g. Fiore et al. 2008a, for further details).
The
/
values as a function of the
luminosity for the SWIRE
mJy sample and for our EDOG sample are compared in Fig. 3a.
It can be seen that our selection criteria (
/
and
mJy) are well suited to find highly obscured, high-luminosity
(
erg s-1) sources. Furthermore, the combination of large
and extreme MIR colors guarantees that the vast majority of our EDOGs are AGN-dominated in the MIR band (e.g. Polletta et al. 2008a; Dey et al. 2008; Donley et al. 2008).
Our finding is therefore in agreement with the correlation between
/
and MIR luminosity reported by Fiore et al. (2008a) for obscured AGNs (see their Figs. 2a and b).
On the other hand, highly luminous (i.e.
> 1044 erg s-1) MIR sources
with a
/
value in the range 10-100 predominantly show a continuum-dominated MIR spectrum with PAH
emission features, suggesting they are composite ULIRG systems containing both a
dominant AGN and a starburst component (Sajina et al. 2007; Yan et al. 2007). MIR-selected,
luminous type 1 quasars typically fall in this region of the
/
-
diagram (e.g. Brand et al. 2007; Fiore et al. 2008b; Polletta et al. 2008a).
Finally, sources at lower luminosities (i.e.
< 1044 erg s-1) can be likely associated with Seyfert-like AGNs (obscured and unobscured) and starburst galaxies.
Table 3: X-ray spectral properties of the 23 sources showing a detection in at least one of the two energy bands, S and H.
In Fig. 3b is shown the distribution of the
luminosities as a function of the redshift for our EDOGs sample.
For comparison we also plot the sources from the CDFS and COSMOS fields (Fiore et al. 2008a,b) with
/
.
The combination of our selection criteria with the large area of the SWIRE survey allows us to collect the most IR luminous sources at all redshifts: objects in our EDOG sample are, on average,
0.5 (
1.5-2) order of magnitude brighter than those selected in the 2 (<0.1) deg2 COSMOS (CDFS) survey.
Assuming that EDOGs host an AGN at their center, it is possible to give a rough estimate of the bolometric luminosity (
)
using the relationship log
= log
+ 1.16 from Elvis et al. (1994). In particular, Polletta et al. (2008a) point out that this relationship likely provides only a lower limit on real value of
.
By proceeding in this way we obtained values of bolometric luminosities in the range 10
.
EDOGs are therefore placed at the brightest end of the bolometric luminosity function of galaxies in the Universe.
3 X-ray observations
3.1 Data reduction
We searched for all Chandra and XMM-Newton archival observations available for the SWIRE fields listed in Sect. 2.1 as of February 2008.
The X-ray coverage of the sources in our sample is highly inhomogeneous, ranging from shallow 5 ks Chandra observations of the
field, up to medium-deep observations in the XMM-Newton SXDS field (
50 ks) and in the Chandra Lockman Hole region (
70 ks).
Twenty-one sources fall within the FOV of an XMM-Newton observation, while the sources with Chandra data are twenty-three.
Event files of the Chandra observations were retrieved from the Chandra X-ray Center (CXC) via the Web ChaSeR (Chandra Archive Search and Retrieval) interface.
The data reduction was performed following the standard procedures outlined in the Science Analysis Threads for ACIS data at the CIAO Web site.
Standard XMM-Newton SAS tasks epproc and emproc (SAS v. 7.1.0) were used to linearize the PN and MOS event files.
The event files were processed using the latest calibration files
and cleaned up from effects of hot pixels and cosmic rays contamination.
X-ray events corresponding to patterns 0-12 (0-4) for the MOS(PN) cameras were selected.
XMM-Newton observations sometimes experience high background flaring periods due to the crossing of clouds of soft (i.e. E<1 MeV) protons, produced by the sun and projected towards the Earth in the solar wind, that are accelerated by magnetic reconnections in the magnetosphere.
In order to produce a cleaned event file we removed all the high background periods identified in
light curves at energies >10 keV, where the contribution from the emission from real X-ray sources is minimized.
3.2 X-ray fluxes
The count-rate in the 0.5-10 keV band (indicated as T band hereafter) was obtained by the command sosta in the Ximage package.
Using this command, we derived the number of source counts within a box region with a 5 (15) arcsec side
for Chandra (XMM-Newton) data centered on the position of the optical counterpart of the X-ray source and
also corrected them for the effects due to the mirror vignetting and the off-axis point spread function degradation.
We repeated the same procedure to estimate the background by the extraction of the background counts from a source-free region around the target. We were thus able to infer the background-subtracted source count rate and its statistical significance. A 3
upper limit is automatically calculated by sosta if the detection significance is less than 10-3. According to this procedure, we obtained 25 positive detections and 19 upper limits in the T band for the 44 sources in our sample.
The count rate was subsequently converted into flux using the PIMMS tool. As underlying spectral model we assumed a power law (with a photon index fixed to
= 1.9) modified by absorption from neutral gas with a rest-frame column density of
= 1023 cm-2. This is an intermediate
value for
/
selected DOGs (see below).
Table 2 lists the resulting T band count rate and flux along with the effective exposure time for each source in the sample. The latter has been obtained from the nominal exposure time of the observation after correction for the vignetting factor at the source position on the detector.
We found that 45% of the sources have an effective exposure time that is
15 ks.
Fluxes in the T band span from
to
erg s-1 cm-2with a distribution mean value of
erg cm-2 s-1.
4 X-ray spectral properties
In order to characterize their spectral properties (with particular emphasis on the absorbing column density ), we selected the sources showing a detection in at least one of the following energy bands: soft (i.e. 0.5-2 keV, S band hereafter) or hard (i.e. 2-10 keV, H band).
This criterion indeed allows to derive the hardness ratio (HR) and, therefore, provide an estimate
of the absorption column density even if just a handful of counts are collected during an observation.
This condition was satisfied by 23 out of 44 sources, i.e.
50% of the sample and,
for the 12 brightest sources, i.e. with
40 net counts in the T band, we were also able to perform a basic spectral analysis (we refer hereafter to this sample of 23 EDOGs as the ``X-ray sample'').
4.1 Spectral analysis
All sources detected with a number of T band counts 40 are observed with XMM-Newton. Spectra were acquired from the cleaned event files of each EPIC camera, i.e. PN, MOS1 and MOS2.
Source counts were extracted from circular regions with a radius in the range of 10-20 arcsec. The radius was chosen in order to optimize the signal-to-noise ratio (S/N) for each source, but in same cases it was limited by the presence of nearby sources or
CCD gaps.
Appropriate response and ancillary files for all the EPIC cameras were created using RMFGEN and ARFGEN tasks in the XMM-Newton SAS, respectively.
In case of multiple observations of the same source, we combined the source and the background spectra from different exposures using the tool mathpha in order to improve the photon statistics.
In doing so we also computed the averages of Response Matrix Files and Ancillary Region Files, weighted for exposure times, using the tools addrmf and addarf.
We performed the spectral analysis for 12 sources having net T band counts in the range 40-600 (with a median value of 80 counts). All the fits discussed in this paper include the absorption due to the line-of-sight Galactic column density (Dickey & Lockman 1990).
For the nine sources with a number of counts <150 we used the C-statistic minimization method (Cash 1979),
being particularly well-suited for low-count spectra.
This is a maximum likelihood method that allows the use of unbinned data. Since it assumes that the counts collected in a given channel follow a Poisson distribution, it cannot be applied to background-subtracted data. Therefore, before performing the spectral modeling we accurately parametrize the background spectrum, extracted from a large source-free region for each source, with a model consisting of a broken power-law plus Gaussian lines, as suggested by Lumb et al. (2002). Then we included this background model in the source spectral fit (fixing the values of the parameters and rescaling the normalizations to the area of the source extraction region).
We assumed an absorbed power-law as spectral fitting model.
Due to the low number of counts we decided to fix the power-law photon index
to 1.9, i.e. the average value found for large samples of quasars in the 2-10 keV band (e.g. Reeves & Turner 2000; Piconcelli et al. 2005, and references therein)
For the remaining three sources (namely source Nos. 6, 39, 44) with >150 net spectral counts we were able to perform the spectral analysis using the
minimization method.
In particular, the large number of collected counts for the source No. 6 (SWIRE2 J021749.00-052306.9) allowed to significantly detect an Fe K
emission line at
6.4 keV in this high-luminosity (L2-10
erg s-1) obscured quasar at
1 (see Sect. 6.2 for further details).
Table 3 reports X-ray counts and spectral parameters for these 12 sources (labeled with SP). We also report the observed L2-10, the absorption corrected L2-10 and the ratio
/L2-10 value for each of them.
![]() |
Figure 4:
a) Left panel: comparison of the best fit values of |
Open with DEXTER |
4.2 Hardness ratio analysis
For the 11 sources detected with <40 T band counts we calculated the background-subtracted count-rate in the S and H band.
We computed the HR value, defined by the formula -
/
+
,
where
(
)
is the net count-rate in S (H) band.
The HR provides an indication about the flatness of an X-ray spectrum and, therefore, it can be used for an estimate of the strength of X-ray obscuration, since the photoelectric absorption mainly affects the soft X-ray portion of the spectrum (i.e. photoelectric absorption cross-section
).
To evaluate the absorbing column density ,
we created a simulated spectrum adopting as model a power law (with photon index
)
modified by a rest-frame absorption. We varied the column density
in the range
(
) < 24 2 by steps of
log(
) = 0.2, and the redshift in the range
by steps of
.
By computing the HR of the simulated spectrum for each value of
and z, we created a set of curves corresponding to different values of
in the plane z versus HR that allowed us to evaluate the absorption
for each source.
Thus, we were able to compute both the observed and the absorption corrected L2-10.
For 2 source (Nos. 8 and 31) with a positive detection only in the H band, we obtained a lower limit on HR and, hence,
.
The latter was used in the computation of the absorption corrected luminosity
Table 3 reports X-ray counts and spectral parameters for these 11 source (labeled HR).
Previous works (e.g. Perola et al. 2004) suggested that the values obtained through the HR technique tend to be overestimated.
In order to check the reliability of the HR technique in case of a small number of counts, we compared the
inferred by spectral analysis with the
value obtained using the method outlined above for the subsample of 12 sources (see Sect. 4.1). The values of
determined for each sources are in good agreement within errors as shown in Fig. 4a. Notably, the only source whose
estimate does not match the one-to-one relationship is the Compton-thick QSO2 candidate SWIRE2 J022003.95-045220.4 (i.e. source No. 12, see Appendix for further details).
This suggests that, due to the intrinsically different spectral shape of Compton-thick sources, the HR analysis can lead to underestimate the
value in these sources.
![]() |
Figure 5:
a) Left panel: F2-10/ |
Open with DEXTER |
4.3 Absorption properties
All the sources in the X-ray sample show a hard X-ray luminosity value typical of AGNs, i.e. L2-10
erg s-1.
We found that the vast majority (i.e.
95%) of the 23 EDOGs with X-ray spectral data are consistent with being obscured by neutral gas with an intrinsic column density of
22 2 (Fig. 4b). Unfortunately, the combination of large z values and, mostly, the limited photon statistics do not allow us to test the presence of warmly ionized absorbing gas in our sources.
In particular the fraction of X-ray obscured EDOGs showing a cold absorption column density strictly larger than 22 2 is
55%. The uncertainties on this fraction is due to the large errors on the
value for the sources detected with a small number of counts.
Table 4 lists the source breakdown according to the level of absorption. This table, as well as Fig. 4b, also reveals the common presence (
55%) of QSO2s among the EDOGs in the X-ray sample. This finding highlights the efficiency of our MIR selection criteria in collecting a large number of these elusive AGNs, which are easily missed in optical and X-ray surveys. This makes promising their application to large samples of X-ray well-exposed EDOGs with the goal of studying the properties of the absorption at the brightest end of the AGN X-ray luminosity distribution.
It is also worth noting that the sources with the highest absorption column densities fall within the QSO2 locus in the
versus 2-10 keV luminosity plane.
To further test the real capability of our MIR selection criteria in gathering QSO2s, we compared these results with those obtained for
a control sample consisting of 20 randomly-selected sources from the same SWIRE fields taken into account in our work.
We picked out these sources for having a hard X-ray flux range (
10-15-10-13 erg cm-2 s-1) and a redshift distribution similar to those of our X-ray sample, regardless of their MIR and optical properties.
Photometric redshifts for these sources are taken from Rowan-Robinson et al. (2008) (13 sources) or have been derived from optical and IR photometric data (7 sources).
Spectral parameters were derived by the individual X-ray spectral analysis of each source.
The fraction of sources in the control sample with
cm-2 is about 15%. Finally, the fraction of QSO2s is
5%.
Table 4: The source breakdown according to the absorption properties. See Sect. 4.3 for details.
4.4 X-ray to optical flux ratio
Optically-selected AGNs typically show an F2-10/
flux ratio, in the range 0.1-10, while
20% of AGNs selected in hard X-ray deep surveys has an F2-10/
(see Alexander et al. 2001; Fiore et al. 2003). Follow-up studies of
these EXOs found that most of them are dust-enshrouded AGN at
(e.g. Maiolino et al. 2006, and references therein).
Fiore et al. (2003) discovered the existence of a striking correlation between F2-10/
and 2-10 keV luminosity for non type 1 AGNs whereby
the most luminous X-ray sources tend to show the highest F2-10/
values.
This is due to the fact that Compton-thin AGNs have their nuclear optical/UV light strongly attenuated differently from the X-ray light and, therefore, their F2-10/
represents the ratio between the X-ray flux from the AGN and the host galaxy star light, which spans a small luminosity range.
Furthermore, the F2-10/
value is a proxy for the absorption in AGN (i.e. the higher is F2-10/
,
the higher is
)
and, therefore, can be used to select obscured AGNs.
![]() |
Figure 6:
a) Left panel: |
Open with DEXTER |
In order to investigate the effects of our MIR/Optical selection on the distribution of F2-10/,
we computed this ratio for the sources in our sample
(see Table 2).
In Fig. 5a we plot the distribution of F2-10/
as a function of the 2-10 keV luminosity for the X-ray sample.
A clear correlation between F2-10/
and hard X-ray luminosity is found, even if characterized by a quite large scatter, in particular for the most X-ray absorbed sources. The solid line in the plot represents the best log(F2-10/
)-log(
)
linear regression obtained by Fiore et al. (2003) from a large sample of X-ray selected, optically-obscured AGNs.
The F2-10/
value as a function of the column density
for our X-ray sample is displayed in Fig. 5b. Notably, all but two sources with
show X-ray absorption column densities consistent with being
2. We also report in this diagram the values obtained for a sample of the brightest 99 hard X-ray sources detected in the XMM-LSS survey (e.g. Garcet et al. 2007). This survey covers
6 deg2 and, therefore, is one of the surveys with the widest coverage area in X-rays. This sample of 99 optically identified sources (half of sources show a z>0.5) covers roughly the same range of 2-10 keV flux (
erg cm-2 s-1) of our X-ray sample, and the
values were obtained from individual spectral analysis. From this figure it can be seen that EDOGs have higher values of
and that also the
values are significantly higher, as expected on the basis of Fiore et al. (2003) and Cocchia et al. (2007) results.
5 Discussion
5.1 Efficiency of our MIR selection criteria to collect QSO2s
Our study was aimed at testing the capability of selection criterion (
/
and F(24
)
mJy)
to isolate luminous, X-ray absorbed quasars at
.
The discovery of a QSO2 in
55% of EDOGs in the X-ray sample
is the most interesting result of our work (e.g. Table 4).
It is worth stressing that all the
EDOGs with column densities consistent with being Compton-thick
fall in the QSO regime, i.e. with a L2-10 >1044 erg s-1.
We have thus tested the efficiency of our MIR technique
to recover a large number of these elusive objects missed in traditional
surveys. They remained poorly-sampled even by
blind-searches performed in the deepest X-ray surveys
(e.g. Fabian 2002; Tozzi et al. 2006). Our choice of sampling at
m)
mJy appears crucial to identify very luminous
AGN-powered objects (as confirmed by the L2-10 > 1043 erg s-1 measured for the sources in the X-ray sample) and minimize any
strong contamination by non-AGN emission (e.g. Weedman et al. 2006;
Dey et al. 2008; Sacchi et al., in prep.).
Our findings also highlights the excellent opportunity
represented by this simple approach in enabling many important studies
concerning the nature of this dust-enshrouded, bright AGN population at high
redshifts. In particular, the systematic application of this method to wide-field MIR surveys with appropriately deep and
homogeneous X-ray coverage will provide
useful insights on largely unexplored issues such as: (i) the demography
and surface density of highly obscured quasars; (ii) the multiwavelength
SED of distant QSO2s; (iii) the distribution at the largest X-ray luminosities; (iv) the luminosity
function of Compton-thin and Compton-thick quasars at z> 1-2,
i.e. near the peak of the AGN activity. Our work may be, therefore, a pathfinder study
for the next-generation observatories, when large areas of the sky will be surveyed at faint hard X-ray
fluxes.
Our study represents the first attempt to characterize in details the X-ray spectral properties of EDOGs. As expected on the basis of previous works based on analogous (although less severe) MIR selection criteria (Fiore et al. 2008a,b; Georgantopoulos et al. 2008; Polletta et al. 2008a) these extremely red sources are found to harbor X-ray type 2 AGNs. However, this MIR bright galaxy population deserves further and deep investigations in X-rays since it constitutes by far one of the most intriguing recently-discovered class of extragalactic objects. Extraordinary bolometric luminosities, high dust content and redshift distribution of these sources lead to believe that the EDOG phenomenon can be potentially associated with the early stages of quasars evolution as envisaged by Silk & Rees (1998).
Interestingly, as stressed by Dey et al. (2008), the space density of the F(24 m)
1 mJy and
/
DOGs is comparable to that of unobscured quasars derived from the MIR luminosity function at z=2 (Brown et al. 2006).
The brightest DOGs likely are AGN-powered sources deeply obscured both in the optical and soft X-rays and, therefore, cannot be revealed by traditional optical/soft-X surveys, which miss a sizable fraction of the accretion history of the Universe.
In this paper we focus on the most extreme DOGs (representing 10% of the
mJy and
/
DOGs)
as a first step in the exploration of the X-ray spectral properties of this class of MIR bright, dust enshrouded galaxies.
We indeed found that the vast majority of the EDOGs harbor an obscured, intrinsically luminous X-ray AGN. It would be very interesting to extend this study to the whole DOG population, in order to provide an accurate estimate of the AGN space density at high z.
5.2 MIR luminosity versus hard X-ray luminosity
In Fig. 6a we compare the 2-10 keV luminosity with the 5.8 m luminosity for sources at the high-luminosity end of the AGN population.
Objects from our X-ray sample are indicated with filled circles, while the shaded area represents the region of L2-10-
plain derived for AGNs in the local Universe by Lutz et al. (2004). In the figure are also plotted data taken from other samples of QSO2s available in the literature. In particular, we have included: (i) three QSO2s from the Polletta et al. (2008a) sample of z>1 objects with F24>1 mJy,
/
,
> 1012
and 10
1023 2 (namely LH_02, LH_A5 and LH_A6; filled squares); (ii) six Compton-thick
2 QSO2 candidates from Alexander et al. (2008) which are extremely weak (or even undetected) in deep X-ray exposures but show evidence for strong AGN activity in the optical and/or MIR band (filled pentagons); and (iii) eight X-ray selected QSO2s (
L0.5-10>1044 erg s-1 and
> 1022 2) from the Sturm et al. (2006) sample (filled triangles).
We have also included in the plot two well-known, very powerful (
>
1012
)
dust-enshrouded QSO2s, i.e. IRAS 09104+4109 (empty star, z=0.442; e.g.
Hines & Wills 1993; Piconcelli et al. 2007a) and 4C 23.56
(empty triangle, z=2.48; e.g. Johnson et al. 2006; Seymour et al. 2008).
Sources absorbed by a column density of Log(
)
24 2 are marked with empty squares.
The solid line in this plot represents the linear regression derived for our sample:
once 3 sources with L2-10 <1043 erg s-1 were excluded, since their weakness in the hard X-ray band (see Table No. 3) suggests that their












As expected, we found a correlation between the two luminosities. In particular this relationship is not linear, and the
/L2-10 ratio increases with increasing
.
The shift at the highest MIR luminosities from the local Universe relationship suggests a possible evolution of the
/L2-10 with redshift.
Nonetheless it is evident from this plot that there is a significant scatter in this relationship.
This probably means that we have put together a mix of different classes of sources, characterized by different properties and dominant energy mechanisms.
The line resulting from Eq. (2) divides the
-L2-10 plane in two regions.
The region above this line includes all the X-ray selected QSO2s from Sturm et al. (2006) and the locus typical of AGNs in the local Universe,
which was derived by Lutz et al. (2004) from a sample of optically/X-ray bright Seyfert 2 galaxies.
Selection effects therefore favor this region to be inhabited by the X-ray/optically ``loud'' objects. All the sources from the Alexander et al. (2008) fall in this part of the diagram reflecting the fact that they are optically bright enough to be robustly identified via optical spectroscopy even if they likely are Compton-thick QSO2s. The intrinsic X-ray luminosities of these sources was indeed estimated on the basis of the luminosity of optical emission lines. Two out of three objects taken from the Polletta et al. (2008a) sample showing moderate absorption in X-rays (
1023 cm-1) are located in this area.
Interestingly, Sturm et al. (2006) emphasize the absence of dust emission and absorption features in the MIR spectra of their QSO2s unlike ULIRGs, which are very dusty systems with powerful circumnuclear starbursts. This points out that these highly luminous X-ray quasars were able to sweep away dust envelopes from their nuclear environment.
In the context of the evolutionary sequence described in the previous section, it seems possible that most of the objects falling in this region of the
L2-10-
plane have already passed through the initial dust-enshrouded phase of quasar activity/galaxy assembly (e.g. Silk & Rees 1998) and are going to evolve into classical ``blue'' QSOs traditionally collected in optical/UV surveys. IRAS 09104+4109 and 4C +24.36, which are respectively IR- and radio-selected, lies outside the locus of local Universe AGNs.
On the contrary, the majority of the sources located below the solid line in Fig. 6a might be objects with dust-absorbed powerful (ULIRG-like) nuclei in which the star-formation activity could significantly contribute to their MIR emission or characterized by intrinsic extreme MIR/X luminosity ratios (e.g. Polletta et al. 2008b; Piconcelli et al. 2009, in preparation).
Accordingly, these sources may be the tip of a larger population representing the link between X-ray-selected QSO2s with featureless MIR spectra (e.g. Sturm et al. 2006) and the submillimeter galaxies (SMGs) (e.g. Alexander et al. 2005; Pope et al. 2008b) which are well-known sites of strong star-formation activity.
In particular, these EDOGs could belong to the minority fraction of SMGs with optical-IR emission dominated by a dust-obscured quasar (e.g. Vignali et al. 2009) with MIR colors overlapping those of MIR power-law continuum AGNs (Yun et al. 2008; Pope et al. 2008a). Thus, even if the position of a source in Fig. 6a strongly depends on the criterion of selection adopted, EDOGs may therefore represent a mixed-bag population that approximately come in two flavors, i.e. ``X-ray loud'' QSO2-like and SMGs-like, which essentially indicate different stages of quasar evolution.
However, caution must be used before drawing any firm conclusion from these plots because of the uncertainties on the intrinsic and the AGN luminosity due to the low quality of X-ray data and the use of photometric redshifts for most of the sources considered here.
Figure 6b shows the observed 2-10 keV luminosity (
)
as a function of
for the same sources plotted in Fig. 6a.
The solid line represents the relationship obtained from the X-ray sample (i.e. Eq. (2)) which can be used to estimate the intrinsic hard X-ray luminosity. The dotted and dashed line represent the observed
-
ratio derived from Eq. (2) in case of a column density of
= 1023 and
= 1024 cm-2, respectively.
Downward-pointing arrows in the plot indicate the 3
upper limit on
for the X-ray undetected sources in our sample lying in the AGN region of the IRAC color-color diagram (Fig. 1b), i.e. sources showing a power law spectrum and very red MIR colors (log(F5.8/
F3.6) > 0.2 and log(F8/
F4.5) > 0.2). Hereafter we refer to this group of 15 sources as the Red Sample.
For these 15 X-ray undetected AGNs, we converted the upper limit on the 2-10 keV flux (assuming a photon index of
)
into the 2-10 keV observed luminosity
.
The intrinsic L2-10 for these sources was determined from their
,
using the relationship in Eq. (2).
By the comparison between the observed and the intrinsic values of hard X-ray luminosity we were able to provide a basic estimate of the column density distribution in the Red Sample. We found that: (i) 93% of these sources show a
cm-2; (ii) 27% of them show a column density
>1024 cm-2 and are Compton-Thick AGN candidates; and (iii) QSO2s-like objects account for 27% of the total (see Table 5).
Since 75
of the Red Sample sources have X-ray effective exposures of
5 ks, the
upper limits are weakly constrained. Accordingly, the real values of the column density may likely be significantly higher than the
values derived here, which should be considered as lower limits.
A handful of objects in the Red Sample, namely Nos. 2, 3 and 38, were not
detected in X-ray medium/deep (i.e. >15 ks) observations (big downward-pointing arrows in Fig. 6b).
This fact makes them very interesting objects since they may likely be
more obscured than those found in the X-ray sample.
In the MIR color-color diagram presented in Fig. 1b,
these EDOGs indeed occupy the area with the extreme red colors, i.e.
and
,
that Polletta et al. (2008a) consider as the typical one of highly luminous and obscured quasars at
.
In particular, we derived an L2-10
erg s-1 for sources Nos. 2 and 3
and L2-10
erg s-1 for the No. 38, and the estimated
value is
1024 cm-2 for all these three sources, making them very likely Compton-Thick QSO2 candidates.
Table 5: The source breakdown for the X-ray sample, the Red sample, and the total sample. See Sect. 5.2 for details.
Merging the results for the Red Sample with those previously found for the X-ray sample, we obtained a total sample of 38 sources for which has been possible to infer information on and hard X-ray luminosity (see Table 5).
This allows us to improve the statistics and get sounder insights on the X-ray properties
of EDOGs, which basically turn out to host heavily obscured AGNs (
20% are consistent to be Compton-thick) with high X-ray luminosity (
40% are QSO2s).
6 Summary and conclusions
We exploited the multiwavelength coverage of the large area SWIRE survey to
collect low-surface density/high luminosity (i.e. 10
)
objects
on the basis of their extreme MIR/Optical flux ratio (
/
)
and bright MIR flux density, i.e.
mJy (referred to here as EDOGs).
Cross-correlating with publicly available XMM-Newton/Chandra observations of 5 SWIRE fields we defined a final sample of 44 EDOGs with X-ray coverage (see Table 1 and Sect. 2).
Our study has been focused on the understanding of largely unexplored X-ray spectral properties of this recently-discovered class of sources.
The main results of this paper can be summarized as follows:
- All but one sources fall within the AGN region of the MIR color-color diagram proposed by Lacy et al. (2004). In particular, the vast majority of EDOGs lies in the obscured luminous QSO locus (Polletta et al. 2008a) characterized by very red colors (Fig. 1b).
- The very bulk of the redshifts, regardless of they are spectroscopic or photometric, are in the range 0.7
z
2.5 (e.g. Sect. 2.3).
- Thanks to the large area sampled (
6 deg2), our EDOGs are on average more luminous in the MIR band by a factor of
0.5-2 orders of magnitude than their equivalents in the CDFS and COSMOS surveys (see Fig. 3b).
- Twenty-three sources (the so-called X-ray sample, e.g. Sect. 4) have a positive detection either in the soft or in the hard band. The non-detection of the remaining sources in the sample can be mainly attributed to the inadequate X-ray coverage (most of the undetected sources indeed have an exposure
5 ks). All the X-ray detected sources have a 2-10 keV luminosity L2-10 > 1043 erg s-1, whereby they fall well within the AGN X-ray luminosity range (Fig. 4b).
- Our analysis reveals that 50% of the sources show a column density strictly larger than 1022 2. Furthermore
50% of them can be classified as QSO2s on the basis of their absorption properties and X-ray luminosity (i.e. L2-10 > 1044 erg s-1).
- At least three EDOGs at z >1.4 were found to be very likely Compton-thick QSO2 candidates (e.g. Sect. 5.2). Moreover such deeply buried AGNs may be present in a sizable fraction (
30%, e.g. Table 5) of the X-ray undetected sources with extreme red colors (log(F5.8/ F3.6) > 0.2 and log(F8/ F4.5) > 0.2), as suggested by their large MIR luminosity (e.g. Sect. 5.2). They therefore deserve further investigation by future properly-planned, deeper X-ray observations.

It seems probable that EDOGs might be an early dust-embedded phase of the AGN activity, possibly associated with merger events
(e.g. Georgantopoulos et al. 2008; Brodwin et al. 2008). Accordingly, an interesting aspect that should be examined in the next future is the FIR/sub-mm emission of these intriguing sources in order to evaluate the contribution of the star-formation activity to their large bolometric luminosity, i.e. using an approach similar to that by Pope et al. (2008a) for a sample of DOGs at low
levels (
Jy).
Acknowledgements
The authors are very grateful to the referee for the exhaustive list of helpful comments and her/his prompt response. We thank for support the Italian Space Agency (contracts ASI-INAF I/023/05/0, ASI I/088/06/0 and PRIN MIUR 2006-02-5203). The authors kindly thank Simonetta Puccetti, Mari Polletta and Roberto Maiolino for many useful comments. Based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. This research has made use of software provided by the Chandra X-ray Center (CXC) in the application package CIAO. This research has made use of the NASA/IPA Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.
Appendix A: Notes on individual sources
- Source No. 6. This is a narrow-line quasar at z=0.914, as shown in on Sect. 2.3.1 on the basis of a Subaru/MOIRCS spectrum.
This source is the brightest in our sample, with a 2-10 keV flux of
10-13 erg cm-2 s-1. We fitted the EPIC data with a continuum model consisting of a power law modified by an intrinsic neutral absorption component, plus a power law with the same slope and Galactic
at soft energies (
2 keV). The latter accounts for the ``soft-excess'' component, i.e. a mixture of photoionized/scattered emission superimposed upon the primary obscured continuum typically observed in most X-ray obscured AGNs (Turner et al. 1997; Piconcelli et al. 2007b; Cappi et al. 2006). Given the redshift of this EDOG, we also included in the best-fit model a Compton reflection component with R = 1 (where R is the solid angle in units of 2
subtended by the reflecting material) as commonly detected in the rest-frame 10-50 keV spectra of AGN (e.g. Risaliti 2002; Reeves et al. 2006). We obtained
and NHz
10.8
22 2, implying a 2-10 keV luminosity of L2-10 = 5.1
1044 erg s-1. The most remarkable result of our analysis is, however, the discovery of an Fe K
emission line in the X-ray spectrum of this QSO2 (Fig. A.1a). This line (
keV) is required by the data at >99.9% confidence level and shows an equivalent width of
eV. It represents one of the best examples of Fe K
line in a high-z QSO2s detected so far.
- Sources Nos. 12, 32 and 33. These sources (with z=1.443, 2.54 and 2.31 respectively) show most of their X-ray counts at E> 2 keV (e.g. Figs. A.1b, A.2a and A.2b) and they are therefore promising deeply buried (possibly Compton-thick) QSO2 candidates.
In particular, source No. 12 was detected with
40(30) net counts in the 0.5-10 (2-10 keV) keV band, i.e. at the limit of Cash statistic reliability for the spectral fitting procedure (Sect. 4.1). Fitting these data with a simple spectral model consisting of a power-law with fixed photon index (
) modified by an intrinsic absorption component, we obtained a lower limit on the column density of
cm-2. The corresponding X-ray absorption-corrected luminosity is L2-10
erg s-1, making this EDOG one of the best Compton-thick QSO2 candidates in our sample.
Sources Nos. 32 and 33 have only
20 net counts, even if they were detected in deep (
70 ks) Chandra exposures, and almost all of them fall in the hard X-ray band. The HR of these sources are therefore very hard (
0.5, e.g. Table 3) and suggest that they are also highly obscured quasars. We derived column density values of
2 for both sources (but keep in mind the large uncertainties on the HR value), confirming thus the results reported in Weedman et al. (2006, sources A4 and A8 in their paper). Beyond providing the spectroscopic redshifts, these authors carried out a detailed multiwavelength study for both sources that allowed to classify them as AGN-powered ULIRGs.
- Source No. 27. This source was undetected in a
70 ks Chandra exposure, being thus the only EDOG without useful X-ray information in such a deep pointing and, possibly, the weakest X-ray sources in our sample. Its IRAC colors suggest that a sizable part of the IR emission in this source can be powered by star formation (e.g. Sect. 2.2 and Yun et al. 2008). According to Rowan-Robinson et al. (2008), the MIR Spectral Energy Distribution (SED) is indeed dominated by the starburst component. This hypothesis is also supported by the large F24/F8 ratio (
20) measured for this EDOG. Such a value is similar to those reported for the infrared luminous galaxies at
in Yan et al. (2007), for which both deeply embedded starburst and AGN are believed to significantly contribute to the dust heating (e.g. Polletta et al. 2008a).
![]() |
Figure A.1: a) Left panel: XMM-Newton spectrum of source No. 6 (i.e. SWIRE2 J021749.00-052306.9). The background subtracted spectra of EPIC PN (black), MOS1 (red) and MOS2 (green) are shown. b) Right panel: XMM-Newton spectrum of source No. 12 (i.e. SWIRE2 J022003.95-045220.4). The source spectrum (black stars) is superimposed to the background spectrum (red empty squares). The best fit models for the source (black solid line) and the background (red dotted line) are also shown. |
Open with DEXTER |
![]() |
Figure A.2: a) Left panel: Chandra spectrum of source No. 32 (i.e. SWIRE2 J104409.95+585224.7). b) Right panel: Chandra spectrum of source No. 33 (i.e. SWIRE2 J104528.29+591326.6). Black stars (red empty squares) represent the source (background) binned counts, with 2 (15) counts per bin. |
Open with DEXTER |
References
- Alexander, D. M., Brandt, W. N., Hornschemeier, A. E., et al. 2001, AJ, 122, 2156 [NASA ADS] [CrossRef] (In the text)
- Alexander, D. M., Smail, I., Bauer, F. E., et al. 2005, Nature, 434, 738 [NASA ADS] [CrossRef] (In the text)
- Alexander, D. M., Chary, R. R., Pope, A., et al. 2008, ApJ, in press,[arXiv:0803.0636] (In the text)
- Alonso-Herrero, A., Perez-González, P. G., Alexander, D. M., et al., 2006, ApJ, 640, 167 [NASA ADS] [CrossRef] (In the text)
- Antonucci, R. 1993, ARA&A, 31, 473 [NASA ADS]
- Arnouts, S., Walcher, C. J., Le Fevre, O., et al. 2007, A&A, 476, 137 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Barmby, P., Alonso-Herrero, A., Donley, J. L., et al., 2006, ApJ, 642, 126 [NASA ADS] [CrossRef]
- Best, P. N., Kauffmann, G., Heckman, T. M, et al., 2005, MNRAS, 362, 25 [NASA ADS] [CrossRef] (In the text)
- Brand, K., Dey, A., Weedman, D., et al. 2006, ApJ, 644, 143 [NASA ADS] [CrossRef]
- Brand, K., Dey, A., Desai, V., et al. 2007, ApJ, 663, 204 [NASA ADS] [CrossRef] (In the text)
- Brandt W. N., & Hasinger, G. 2005, ARA&A, 43, 827 [NASA ADS] [CrossRef] (In the text)
- Brodwin, M., Dey, A., Brown, M. J., et al. 2008, ApJ, 687, L65 [NASA ADS] [CrossRef] (In the text)
- Brown, M. J. I., Brand, K., Dey, A., et al. 2006, ApJ, 638, 88 [NASA ADS] [CrossRef] (In the text)
- Cappi, M., Panessa, F., Bassani, L., et al. 2006, A&A, 446, 459 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Cardamone, C. N., Urry, C. M., Damen, M., et al. 2008, ApJ, 680, 130 [NASA ADS] [CrossRef] (In the text)
- Cash, W. 1979, ApJ, 228, 939 [NASA ADS] [CrossRef] (In the text)
- Cocchia, F., Fiore, F., Vignali, C., et al. 2007, A&A, 466, 31 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Comastri, A., 2004, Multiwavelength AGN Surveys, Proceedings of the Guillermo Haro Conference, ed. R. Mujica, & R. Maiolino (Singapore: World Scientific Publishing Company), 323 (In the text)
- Cowie, L. L., Barger, A. J., Bautz, M. W., Brandt, W. N., & Garmire, G. P. 2003, ApJ, 584, 57 [NASA ADS] [CrossRef]
- Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 11 [NASA ADS] [CrossRef] (In the text)
- Daddi, E., Alexander, D. M., Dickinson, M., et al., 2007, ApJ, 670, 173 [NASA ADS] [CrossRef] (In the text)
- Dey, A., Soifer, B. T., Desai, V., et al. 2008, ApJ, 677, 943 [NASA ADS] [CrossRef] (In the text)
- Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 [NASA ADS] [CrossRef] (In the text)
- Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 [NASA ADS] [CrossRef] (In the text)
- Donley, J. L., Rieke, G. H., Pérez-González, P. G., Rigby, J. R., & Alonso-Herrero, A. 2007, ApJ, 660, 167 [NASA ADS] [CrossRef] (In the text)
- Donley, J. L., Rieke, G. H., Pérez-González, P. G., & Barro, G. 2008, ApJ, 687, 111 [NASA ADS] [CrossRef] (In the text)
- Dwelly, T., & Page, M. J. 2006, MNRAS, 372, 1755 [NASA ADS] [CrossRef] (In the text)
- Elvis, M., Wilkes, B. J., McDowell, J. C., et al. 1994, ApJS, 95, 1 [NASA ADS] [CrossRef] (In the text)
- Erlund, M. C., Fabian, A. C., Blundell, K. M., & Crawford, C. S. 2008, MNRAS, 385, L125 [NASA ADS] (In the text)
- Fabian, A. C. 1999, MNRAS, 308, L39 [NASA ADS] [CrossRef] (In the text)
- Fabian, A. C. 2002, ASPC, 258, 185 [NASA ADS] (In the text)
- Feruglio, C., Fiore, F., La Franca, F., et al. 2008, A&A, 488, 417 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Fiore, F., Brusa, M., Cocchia, F., et al. 2003, A&A, 409, 79 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Fiore, F., Grazian, A., Santini, P., et al. 2008a, ApJ, 672, 94 [NASA ADS] [CrossRef] (In the text)
- Fiore, F., Puccetti, S., Brusa, M., et al. 2008b, ApJ, in press,[arXiv:0810.0720] (In the text)
- Fontana, A., D'Odorico, S., Poli, S., et al. 2000, AJ, 120, 2206 [NASA ADS] [CrossRef] (In the text)
- Fontana, A. 2001, Astrophys. Space Sci. Suppl., 277, 535 [NASA ADS] [CrossRef] (In the text)
- Furusawa. H., Kosugi, G., Akiyama, M., et al., 2008, ApJS, 176, 1 [NASA ADS] [CrossRef] (In the text)
- Garcet, O., Gandhi, P., Gosset, E., et al. 2007, A&A, 474, 473 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Georgantopoulos, I., Georgakakis, A., & Akylas, A. 2007, A&A, 466, 823 [NASA ADS] [CrossRef] [EDP Sciences]
- Georgantopoulos, I., Georgakakis, A., Rowan-Robinson, M., & Rovilos, E. 2008, A&A, 484, 671 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Gilli, R., Comastri, A., & Hasinger, G., 2007, A&A, 463, 79 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 2004, ApJ, 600, 580 [NASA ADS] [CrossRef] (In the text)
- Hasinger, G. 2008, A&A, 490, 905 [NASA ADS] [CrossRef] [EDP Sciences]
- Hickox, R., Jones, C., Forman, W. R., et al. 2007, ApJ, 671, 1365 [NASA ADS] [CrossRef] (In the text)
- Hines, D.C., & Wills, B. J. 1993, ApJ, 415, 82 [NASA ADS] [CrossRef] (In the text)
- Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., & Springel, V. 2006, ApJS, 163, 1 [NASA ADS] [CrossRef] (In the text)
- Houck, J. R., Soifer, B. T., Weedman, D., et al. 2005, ApJ, 622, L105 [NASA ADS] [CrossRef] (In the text)
- Johnson, O., Almaini, O., Best, P. N., & Dunlop, J. 2007, MNRAS, 376, 151 [NASA ADS] [CrossRef] (In the text)
- Lacy, M., Petric, A. O., Sajina, A., et al. 2007, AJ, 133, 186 [NASA ADS] [CrossRef] (In the text)
- Lacy, M., Storrie-Lombardi, L. J., Sajina, A., et al., 2004, ApJS, 154, 166 [NASA ADS] [CrossRef] (In the text)
- La Franca, F., Fiore, F., Comastri, A., et al. 2005, ApJ, 635, 864 [NASA ADS] [CrossRef] (In the text)
- Lonsdale, C. J., Smith, H. E., Rowan-Robinson, M., et al. 2003, PASP, 115, 897 [NASA ADS] [CrossRef] (In the text)
- Lumb, D. H., Warwick, R. S., Page, M., & De Luca, A. 2002, A&A, 389, 93 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Lutz, D., Maiolino, R., Spoon, H. W. W., & Moorwood, A. F. M. 2004, A&A, 418, 465 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Mainieri, V., Rosati, P., Tozzi, P., et al. 2005, A&A, 437, 805 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Maiolino, R., Comastri, A., Gilli, R., et al., 2003, MNRAS, 344, 59 [NASA ADS] [CrossRef]
- Maiolino, R., Mignoli, M., Pozzetti, L., et al. 2006, A&A, 445, 457 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Maiolino, R., Shemmer, O., Imanishi, M., et al. 2007, A&A, 468, 979 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Martínez-Sansigre, A., Rawlings, S., Lacy, M., et al., 2005, Nature, 436, 666 [NASA ADS] [CrossRef] (In the text)
- Martínez-Sansigre, A., Rawlings, S., Bonfield, D. G., et al. 2007, MNRAS, 379, L6 [NASA ADS]
- McMahon, R. G., Walton, N. A., Irwin, M. J., et al. 2001, New Astron. Rev., 45, 97 [NASA ADS] [CrossRef] (In the text)
- Menci, N., Fiore, F., Puccetti, S., & Cavaliere, A. 2008, ApJ, 686, 219 [NASA ADS] [CrossRef] (In the text)
- Merritt, D., & Ferrarese, L., 2001, The Central Kiloparsec of Starbursts and AGN: The La Palma Connection, ASP Conf. Proc., 249, 335 (In the text)
- Norman, C., Hasinger, G., Giacconi, R., et al. 2002, ApJ, 571, 218 [NASA ADS] [CrossRef] (In the text)
- Page, M. J., Stevens, J. A., Ivison, R. J., & Carrera, F. J. 2004, ApJ, 611, L85 [NASA ADS] [CrossRef] (In the text)
- Perola, G. C., Puccetti, S., Fiore, F., et al. 2004, A&A, 421, 491 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Piconcelli, E., Jimenez-Bailón, E., Guainazzi, et al. 2005, A&A, 432, 15 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Piconcelli, E., Fiore, F., Nicastro, et al. 2007a, A&A, 473, 85 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Piconcelli, E., Bianchi, S., Guainazzi, M., Fiore, F., & Chiaberge, M. 2007b, A&A, 466, 855 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Polletta, M., Wilkes, B. J., Siana, B., et al., 2006, ApJ, 642, 673 [NASA ADS] [CrossRef] (In the text)
- Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81 [NASA ADS] [CrossRef] (In the text)
- Polletta, M., Weedman, D., Hönig, S., et al. 2008a, ApJ, 675, 960 [NASA ADS] [CrossRef] (In the text)
- Polletta, M., Omont, A., Berta, S., et al. 2008b, A&A, 492, 81 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Pope, A., Bussmann, R. S., Dey, A., et al. 2008a, ApJ, 689, 127 [NASA ADS] [CrossRef] (In the text)
- Pope, A., Chary, R. R., Alexander, D. M., et al. 2008b, ApJ, 675, 1171 [NASA ADS] [CrossRef] (In the text)
- Pozzi, F., Vignali, C., Comastri, A., et al. 2007, A&A, 468, 603 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Reeves, J. N., & Turner, M. J. L. 2000, MNRAS, 316, 234 [NASA ADS] [CrossRef] (In the text)
- Reeves, J. N., Fabian, A. C., Kataoka, J., et al., 2006, AN, 327, 1079 [NASA ADS] (In the text)
- Risaliti, G. 2002, A&A, 386, 379 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Rosati, P., Tozzi, P., Giacconi, R., et al. 2002, ApJ, 566, 667 [NASA ADS] [CrossRef]
- Rowan-Robinson, M., Babbedge, T., Oliver, S., et al. 2008, MNRAS, 386, 697 [NASA ADS] [CrossRef] (In the text)
- Sajina, A., Lacy, M., & Scott, D. 2005, ApJ, 621, 256 [NASA ADS] [CrossRef] (In the text)
- Sajina, A., Yan, L., & Armus, L. 2007, ApJ, 664, 713 [NASA ADS] [CrossRef] (In the text)
- Seymour, N., Ogle, P., De Breuck, C., et al. 2008, ApJ, 681, L1 [NASA ADS] [CrossRef] (In the text)
- Sekiguchi, K., Akiyama, M., Furusawa, H., et al. 2005, Proceedings of the ESO Workshop, Multiwavelength mapping of galaxy formation and evolution, ed. A. Renzini, & R. Bender, 82 (In the text)
- Silk, J., & Rees, M. J. 1998, A&A, 331, L1 [NASA ADS] (In the text)
- Simpson, C. 2005, MNRAS, 360, 565 [NASA ADS] [CrossRef]
- Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hernquist, L. 2008, MNRAS, 391, 481 [NASA ADS] [CrossRef] (In the text)
- Steffen, A. T., Barger, A. J., Cowie, L. L., Mushotzky, R. F., & Yang, Y. 2003, ApJ, 596, L23 [NASA ADS] [CrossRef] (In the text)
- Sturm, E., Hasinger, G., Lehmann, I., et al., 2006, ApJ, 642, 81 [NASA ADS] [CrossRef] (In the text)
- Surace, J. A., et al. 2004, BAAS, 36, 1450 [NASA ADS] (In the text)
- Tozzi, P., Gilli, R., Mainieri, V., et al., 2006, A&A, 451, 457 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Turner, T. J., George, I. M., Nandra, K., & Mushotzky, R. F. 1997, ApJS, 113, 23 [NASA ADS] [CrossRef] (In the text)
- Ueda, Y., Akiyama, M., Ohta, K., & Miyaji, T. 2003, ApJ, 598, 886 [NASA ADS] [CrossRef] (In the text)
- Vignali et al. 2008, submitted
- Weedman, D. W., Le Floc'h, E., Higdon, S. J. U., et al. 2006, ApJ, 638, 613 [NASA ADS] [CrossRef] (In the text)
- Worsley, M. A., Fabian, A. C., Bauer, F. E., et al. 2005, MNRAS, 357, 1281 [NASA ADS] [CrossRef] (In the text)
- Yan, L., Chary, R., Armus, L., et al., 2005, ApJ, 628, 604 [NASA ADS] [CrossRef] (In the text)
- Yan, L., Sajina, A., Fadda, D., et al., 2007, ApJ, 658, 778 [NASA ADS] [CrossRef] (In the text)
- Yun, M. S., Aretxaga, I., Ashby, M. L. N., et al. 2008, MNRAS, 389, 333 [NASA ADS] [CrossRef] (In the text)
Footnotes
- ... channels
- Source Nos. 8, 9, 25 and 43 are undetected in one (or two) IRAC bands and, therefore, they are not plotted in Fig. 1b.
- ... available
- For 15 sources we were able to compare the photo-z from the SWIRE Photometric Redshift Catalogue with those computed with our method. We found that for
50% (
85%) of the sources, the photo-z values are consistent within
(
).
- ... site
- http://cxc.harvard.edu/ciao
- ... 7.1.0)
- http://xmm.esac.esa.int/sas/
- ... Ximage
- http://heasarc.gsfc.nasa.gov/xanadu/ximage/ximage.html
- ... PIMMS
- http://heasarc.nasa.gov/Tools/w3pimms.html
- ... sample
- For the sources which are undetected in the 2-10 keV band we derived an upper limit in the F2-10/
value extrapolating the H band flux using as underlying spectral model an absorbed (
= 1023 2) power law with
, e.g. Sect. 3.2.
- ... 23.56
- We re-analyzed the XMM-Newton data of 4C 23.56. Our best-fit (an absorbed power law with
cm-2 and
fixed to 1.9 plus a soft (
) unabsorbed power law) is consistent with those reported by Johnson et al. (2006), confirming the QSO2 nature of this high-z radiogalaxy. We derived a 2-10 keV intrinsic(observed) luminosity of L2-10 =
erg s-1.
- ... properties
- For the X-ray sample we considered the percentages derived using the lowest value in 90% confidence level for the
parameter obtained for each source.
All Tables
Table 1: Optical-IR parameters and redshift of the sample.
Table 2: X-ray coverage and 0.5-10 keV band flux of the sample.
Table 3: X-ray spectral properties of the 23 sources showing a detection in at least one of the two energy bands, S and H.
Table 4: The source breakdown according to the absorption properties. See Sect. 4.3 for details.
Table 5: The source breakdown for the X-ray sample, the Red sample, and the total sample. See Sect. 5.2 for details.
All Figures
![]() |
Figure 1:
a) Left panel: 24 |
Open with DEXTER | |
In the text |
![]() |
Figure 2: a) Left panel: near IR spectrum of source No. 1 at z=1.007. b) Right panel: near IR spectrum of source No. 6 at z=0.914. Both spectra were obtained by MOIRCS/SUBARU observations. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
a) Left panel: |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
a) Left panel: comparison of the best fit values of |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
a) Left panel: F2-10/ |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
a) Left panel: |
Open with DEXTER | |
In the text |
![]() |
Figure A.1: a) Left panel: XMM-Newton spectrum of source No. 6 (i.e. SWIRE2 J021749.00-052306.9). The background subtracted spectra of EPIC PN (black), MOS1 (red) and MOS2 (green) are shown. b) Right panel: XMM-Newton spectrum of source No. 12 (i.e. SWIRE2 J022003.95-045220.4). The source spectrum (black stars) is superimposed to the background spectrum (red empty squares). The best fit models for the source (black solid line) and the background (red dotted line) are also shown. |
Open with DEXTER | |
In the text |
![]() |
Figure A.2: a) Left panel: Chandra spectrum of source No. 32 (i.e. SWIRE2 J104409.95+585224.7). b) Right panel: Chandra spectrum of source No. 33 (i.e. SWIRE2 J104528.29+591326.6). Black stars (red empty squares) represent the source (background) binned counts, with 2 (15) counts per bin. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.