A&A 472, 797-803 (2007)
DOI: 10.1051/0004-6361:20077394
F. La Franca1 -
S. Puccetti2, -
N. Sacchi1 -
C. Feruglio3 -
F. Fiore3 -
C. Gruppioni4 -
A. Lamastra1 -
I. Matute5 -
G. Melini1 -
F. Pozzi4,6
1 - Dipartimento di Fisica, Università Roma Tre,
via della Vasca Navale 84, 00146 Roma, Italy
2 -
ASI Science Data Center, via Galileo Galilei, 00044, Frascati, Italy
3 -
INAF - Osservatorio Astronomico di Roma, via Frascati 33, 00040 Monteporzio
Catone, Italy
4 -
INAF, Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
5 -
INAF, Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5,
50125 Firenze, Italy
6 -
Dipartimento di Astronomia, Università di Bologna,
via Ranzani 1, 40127 Bologna, Italy
Received 2 March 2007 / Accepted 2 June 2007
Abstract
Context. The counts of galaxies and AGN in the mid infra-red (MIR) bands are important instruments for studying their cosmological evolution. However, the classic spectral line ratios techniques can become misleading when trying to properly separate AGN from starbursts or even from apparently normal galaxies.
Aims. We use X-ray band observations to discriminate AGN activity in previously classified MIR-selected starburst galaxies and to derive updated AGN1 and (Compton thin) AGN2 counts at 15 m.
Methods. XMM observations of the ELAIS-S1 15 m sample down to flux limits
erg cm-2 s-1 (2-10 keV band) were used. We classified as AGN all those MIR sources with a unabsorbed 2-10 keV X-ray luminosity higher that
1042 erg s-1.
Results. We find that at least about % of the previously classified starburst galaxies harbor an AGN. According to these figures, we provide an updated estimate of the counts of AGN1 and (Compton thin) AGN2 at 15
m. It turns out that at least 24% of the extragalactic sources brighter than 0.6 my at 15
m are AGN (
13% contribution to the extragalactic background produced at fluxes brighter than 0.6 mJy).
Key words: cosmology: observations - infrared: galaxies - galaxies: active - surveys - galaxies: evolution - X-rays: galaxies
To measure the evolution of the density of active galactic nuclei (AGN) and galaxies at various wavelengths is one of the major goals in astrophysics, and extragalactic counts in the mid-infrared (MIR) are important instruments for this. MIR surveys, carried out with the ISOCAM instrument on board the Infrared Space Observatory ( ISO), indicates that the MIR galaxies have evolved significantly faster than deduced from optical surveys (e.g. Elbaz et al. 1999; Gruppioni et al. 2002; Metcalfe et al. 2003; Rodighiero et al. 2004; Pozzi et al. 2004). These results are now confirmed and extended by the data coming from the Spitzer Space Telescope (e.g. Gruppioni et al. 2005; Le Floc'h et al. 2005; Pérez-González et al. 2005; Babbedge et al. 2006; Caputi et al. 2007). On the contrary, the evolution of AGN at MIR wavelengths are similar to what has been measured in the optical and X-ray bands for type 1 AGN (AGN1) and even slower for type 2 AGN (Matute et al. 2002, 2006).
Proper analysis of the evolution of the sources emitting in the MIR is
based on a reliable classification that allows a separation of AGN
(more precisely type 2 AGN; AGN2) from the starburst
galaxies. Unfortunately there is growing evidence that the
classification of AGN based on their optical spectra alone provides an
incomplete and sometimes deceiving description of their true nature.
Indeed, some starbursts, or even normal, galaxies show signatures of a
hidden AGN at non-optical wavelengths (see e.g. Moorwood et al. 1996;
Marconi et al. 2000; Severgnini et al. 2003), while a fraction of
X-ray selected AGN show no signatures of AGN activity in their optical
spectra: the X-ray bright optically normal galaxies (XBONG, see
e.g. Fiore et al. 2000; Comastri et al. 2002; La Franca et al.
2005
and references therein). Combining Chandra and archival X-ray
data, Maiolino et al. (2003) have found that the local total density
of AGN may be a factor 2 higher than estimated from optical
spectroscopic surveys. Up to now the exact proportion of
AGN-dominated MIR galaxies is actually not known. This separation is
also complicated by the existence of mixed systems, where both star
formation and AGN activity significantly contribute to the mid-IR
emission (e.g. Fritz et al. 2006; Pozzi et al. 2007).
Mid-IR studies based on IRAS and ISO observations have
indirectly estimated a contribution of AGN to the cosmic infrared
background (CIRB) that is not larger than 5-10
(e.g. Franceschini
et al. 2001; Xu et al. 2001, 2003). This is the maximum room left in
their models by the strongly evolving starburst population. The total
contribution was, in any case, uncertain since mid-IR selected AGN2
and starburst galaxies were treated as a single population. A better
estimate of the AGN contribution comes from the X-ray band (0.5-10 keV), which offers a better wavelength regime for selecting and
identifying obscured sources (unless the sources are Compton-thick,
cm-2). A cross correlation of
X-ray and IR ISO sources detected in deep Chandra and
XMM-Newton observations allowed Fadda et al. (2002) to estimate that
the maximum fraction of the CIRB produced by AGN is
%. This
result is confirmed by Silva et al. (2004), who used the X-ray AGN
luminosity function (LF) and semi-empirical SEDs (linking the X-ray
to the Infrared) to derive that the contribution of AGN and their host galaxy to the CIRB is 10-20% in the 1-20
m range. These
results recently have been confirmed and extended by the analysis of
joint Spitzer and X-ray data, where a contribution of 3-11% is
derived (see e.g. Treister et al. 2006; Barmby et al. 2006; Brand et al. 2006).
On the other hand, it has been proposed that heavily X-ray obscured sources could be selected in the mid-infrared band where the dust heated by the AGN produces a spectral energy distribution (SED) that peaks at shorter wavelengths than the SED produced by dust heated by starbursts. Recently Spitzer data have been used to find highly obscured AGNs, both trough mid-infrared color selection (e.g. Stern et al. 2005; Polletta et al. 2006; Alonso-Herrero et al. 2006; Lacy et al. 2007) and joint radio-infrared selection (e.g. Donley et al. 2005; Weedman et al. 2006; Martinez-Sansigre et al. 2005, 2006). However, recent studies are showing that the situation is not simple. The selection of AGN with mid-IR colors is showing up as not simultaneously complete and reliable; and, moreover, when X-ray data are available, the AGN mid-IR SED do not directly correspond to the X-ray spectral type (absorbed/un-absorbed; e.g. Barmby et al. 2006; Donley et al. 2007).
In this framework it is interesting to work out more reliable AGN1 and
AGN2 counts in the MIR using the 2-10 keV X-ray luminosity as a
solid indication of nuclear activity. In this work we use XMM-Newton
observations of the S1 area of the European Large Area ISO
Survey (ELAIS hereafter; Oliver et al. 2000; Rowan-Robinson et al.
2004) to revise the optical spectroscopic classification of
the 15 m sample presented by La Franca et al. (2004) and then
to derive new, more reliable extragalactic counts of AGN1 and not heavily
absorbed (
cm-2, not Compton thick) AGN2 in the
0.5-10 mJy regime, where no other "class-separated'' data exist at
the moment.
The MIR and X-ray data are presented in Sect. 2. In Sect. 3, we describe the
X-ray properties of the MIR sources, while the revised AGN 15 m
counts are derived in Sect. 4. Section 5 contains our comments and
conclusions. Unless otherwise stated, all quoted errors are at the 68% confidence level. We assume H0= 70 km s-1 Mpc-1,
,
and
.
The mid-IR selected sources used in this work have been extracted from
the 15m sample selected in the central region 5 of the field
S1 of the ELAIS. The 15
m catalogue in the ELAIS S1 field has
been released by Lari et al. (2001). It covers an area of
4 deg2 centered at
,
and includes 462 mid-IR sources down to
a flux limit of 0.5 mJy. Mid-IR source counts based on this catalogue
have been presented and discussed by Gruppioni et al. (2002).
La Franca et al. (2004) presents the optical and spectroscopic
identifications and classification of a more reliable subsample of 406 sources from the Lari et al. (2001) catalogue. About 80
of these
sources have been optically identified on CCD exposures down to R
23. The spectral classification has been obtained for 90
of the optically identified sample. As discussed by La Franca et al. (2004), due to a different mid-IR flux limit coverage of the
ELAIS-S1 field, the total area has been divided into two regions: the
central and deepest region of S1 (S1-5) reaching mid-IR fluxes
(S15 hereafter) of 0.5 mJy, and the outer region (S1-rest) with a
0.9 mJy flux limit. The S1-5 area is spectroscopically complete at the
97
level down to R=21.6, while S1-rest completeness reaches the
98
level down to R=20.5. In total, 116 sources (29
of the
total mid-IR sample) do not have a spectroscopic identification due to
incompleteness of the follow-up or to the lack of optical counterpart
brighter than R=23. A detailed description of the optical
identification, spectroscopic classification, size, and completeness
function of the different areas used in the ELAIS-S1 sample are
presented and discussed by La Franca et al. (2004).
Classification for AGN-dominated sources in the ELAIS fields was based
on their optical spectral signatures. Sources showing broad emission
line profiles (rest frame
km s-1) were classified as
type-1. Type-2 sources were selected following classic diagnostic
diagrams (e.g. Tresse et al. 1996; Osterbrock 1989; Veilleux &
Osterbrock 1987) that included one or more of the following line
ratios: [NII]/H
,
SII/(H
+[NII]), OI/H
,
[OIII]/H
,
and [OII]/H
when available, depending on the
redshift of the source (e.g. log ([OIII]/H
) > 0.5 and
log ([NII]/H
) > -0.4).
Few objects have been re-observed in the framework of the follow-up
campaign of the Spitzer Wide-Area Infrared Extragalactic Survey
(SWIRE; see e.g. Lonsdale et al. 2003), which includes the ELAIS-S1
field. The ELAIS-S1 field has been observed in the B and V bands down
to 25 Vega mag and in the R band down to 24.5 mag (Berta et al.
2006). A spectroscopic follow-up campaign was carried out in the
period 2004-2006 at the 3.6 m and VLT ESO telescopes (La Franca et al.
in prep.). In comparison with La Franca et al. (2004), two more
15 m sources were spectroscopically identified in the ELAIS-S1-5
field: C15_J003317-431706 is a R=24.3 galaxy showing [OII] emission
at redshift 0.689, while C15_J003447-432447 is an AGN2 at redshift 1.076. Moreover, three sources with previously poor quality spectra
changed their spectroscopic classification: C15_J003330-431553, which
was wrongly classified as a starburst galaxy at z=0.473, showed broad
CIII and MgII emission at z=2.170 typical of AGN1 activity;
C15_J003603-433155, which was classified as AGN2, showed a broad MgII
emission typical of AGN1 activity; and C15_J003622-432826, which was
classified as a starburst galaxy, showed a clear [OIII]/H
ratio
typical of AGN2 activity. In summary, the MIR sample used in this
work contains 72 bona fide extragalactic sources, of which 67 sources
have been optically identified down to
.
Fifty-seven of them
have been spectroscopically identified: six are AGN1, three AGN2,
while almost all of the remaining 48 sources show the emission line
spectra typical of starburst activity. Most of the remaining 15 sources without spectroscopic identification (and for five sources
even optical identification) are so faint that it is very likely that
they are starburst galaxies at redshift higher than 0.5 (see La Franca
et al. 2004).
The region S1-5 of the ELAIS was observed with the XMM-Newton
telescope. The observations and data reduction are described by
Puccetti et al. (2006). The region was observed by four partially
overlapping pointings, covering a total area of 0.6 deg2 (see
Fig. 1). The observations were carried out from May to
July 2003, and each pointing lasted 84-100 ks. Source detection was
performed on Co-added PN+MOS1+MOS2 images accumulated in four energy
bands: 0.5-10 keV, 0.5-2 keV, 2-10 keV, and 5-10 keV. In total 478
sources were detected with a significance level
of
-5. The flux limits are
-16 erg cm-2 s-1 in the 0.5-2 keV band,
-15 erg cm-2 s-1 in the 2-10 keV band, and
-15 erg cm-2 s-1 in the 5-10 keV band.
![]() |
Figure 1:
Mosaic of the four XMM-Newton pointings in the ELAIS-S1 field
obtained by adding up the images of the PN and MOS cameras in
the energy range 0.5-10 keV (Puccetti et al. 2006). The circles are the
position of the 72 15 ![]() |
Open with DEXTER |
The source counts in each camera were obtained using the events files
in the energy range 0.5-10 keV for the PN and 0.3-10 keV for the
MOS. The counts of the two MOS cameras were eventually combined. The
counts of each source were extracted in a circular region with a
radius in the range
-
.
In general, the radius
value was chosen so that the S/N ratio was roughly optimized, but this
choice was limited in a few cases by the presence of nearby sources
or by a peculiar position of the source on the detector, for example
close to a gap in the CCD array. In some cases the source was
detected, and the corresponding counts extracted, only in either the
PN or in one or both of the MOS cameras, because PN and MOS do not
cover exactly the same sky regions and the position of the gaps
differs in the PN, MOS1, and MOS2 CCD arrays.
The background counts for each source were extracted from the nearest source-free region. In doing so, areas near gaps in the CCD array and near the edge of the telescope field of view have been excluded, as well as regions containing hot pixels and other CCD cosmetic defects.
The ancillary response files were generated for each source by means
of the tool ARFGEN, in order to properly correct for
energy-dependent vignetting and point spread function. The response
matrix files were generated with the tool RMFGEN (SAS 6.1.0).
The spectral counts, when higher than about 120, were first
accumulated in energy bins with 20 counts each, from 0.3 keV to 10 keV
in the MOS and from 0.5 to 10 keV in the PN. They were then fitted,
using XSPEC (version 11.3.1) and the
statistic with a model
comprising, in addition to the known galactic absorption
(
20 cm-2): (1) a power law, with two
parameters, normalization and photon spectral index
;
(2)
the absorption
at the redshift of the source. Sometimes a thin
plasma model with the temperature as a parameter was added. In the
case of degrees of freedom lower than 12, a fixed photon index
was assumed. When both PN and MOS data were available,
their relative normalization MOS/PN was left free to vary between 0.8 and 1.2. This interval was conservatively chosen wider than applicable
on-axis, because a fully reliable inter-calibration is still lacking
for sources off-axis. When the spectral counts were lower than about
120, the C statistic (Cash 1979) was used instead, as implemented in
XSPEC (Arnaud 2003
) after background subtraction (see
Alexander et al. 2003a, for a similar procedure) and with 5 counts in
each energy bin (the latter choice was made only for convenience; it
does not impair the correct use of the embedded statistics when using
the above-mentioned XSPEC implementation).
We wish to use the X-ray observations to identify possible
AGN activity among those 15 m galaxies whose optical spectroscopic
identification did no show any AGN signatures. We have thus
chosen a lower limit intrinsic un-absorbed X-ray luminosity of
1042 erg s-1 in order to classify a source as an AGN.
Indeed higher luminosities are typical of AGN activity (see e.g.
Ranalli et al. 2003; Maiolino et al. 2003), while the
most X-ray-luminous star-forming galaxy known, NGC 3256, has a
total X-ray luminosity of
erg s-1 (Lira et al. 2002).
In total, 13 out of the 72 XMM observed ELAIS 15
sources were
detected in the X-ray band. Table 1 shows the X-ray and optical
spectral characteristics for each X-ray detected source. For the
remaining 59 X-ray undetected ISO sources, three sigma-confidence
level upper limits have been measured in the 0.5-10 keV and 2-10 keV
energy bands. Table 2 shows a summary for all 72 ELAIS sources of the
MIR, optical and X-ray fluxes and luminosities. The redshift and
optical spectroscopic classification is given when available. The
X-ray luminosities of the detected sources have been corrected for
absorption.
![]() |
Figure 2:
Left. Redshift versus R-band magnitude distribution of the
15 ![]() ![]() ![]() |
Open with DEXTER |
Twelve of the 13 X-ray detected sources have an optical and spectroscopic identification, while the X-ray detected source C15_J003315-432829 has an optical counterpart with R=21.61 but no optical spectroscopic identification. Figure 2 shows the redshift of the spectroscopically identified ELAIS sources as a function of the optical or MIR fluxes. As it is possible to see, all sources with R-band magnitude fainter than R=20.5 have redshift higher than 0.45. Then, as discussed by La Franca et al. (2004), this (not classified) source (C15_J003315-432829) is likely to be a starburst galaxy or an AGN2 at redshift higher than 0.5. According to this guess, a lower limit of 1043.4 erg s-1 on the X-ray luminosity is derived. This result allowed us to classify this source as an AGN. Anyhow, for the sake of our analyses, even an assumed redshift as low as 0.2 would imply an X-ray luminosity higher than 1042 erg s-1, typical of AGN activity, and does not change our conclusions.
Four out of the six previously optically-classified AGN1 were detected
in the X-ray band and have 2-10 keV luminosities higher than
1043.2 erg s-1. Their X-ray spectral shape confirms this
classification, although two AGN1, C15_J003330-431553 and
C15_J003640-433925, show absorbed spectra with column densities of
cm-2. This result is not so unusual for AGN1 (see
La Franca et al. 2005, and references therein). The remaining two
undetected AGN1 and the three (previously optically classified) AGN2
have 2-10 keV X-ray-luminosity upper limits of
1043.2-1044.0 erg s-1, which are still consistent
with the typical X-ray luminosity of un-absorbed AGNs
(1042-1046 erg s-1).
All the 9 remaining X-ray detected sources have the optical spectra
typical of starburst activity. However, eight of them have X-ray
luminosities typical of AGN activity: five showed absorbed
(
cm-2) X-ray spectra, thus reinforcing the
classification as AGN2, while the other three sources have X-ray
spectra consistent with an unabsorbed power law. Finally, the last
source (C15_J003432-433922) has the X-ray luminosity typical of
starburst activity (
1040 erg s-1). This is also
confirmed by its X-ray spectrum, where a thin plasma model (MEKAL
with
)
has necessarily been added to the
power-law model.
There are 14 sources without X-ray detections and without measured redshift. According to either their R-band optical magnitudes or their R-band magnitude limits (R>24.5), most of these sources should lie at redshifts higher than 0.5 (see discussion in La Franca et al. 2004). In order to derive an upper limit on their X-ray luminosities, we made the conservative assumption that these sources lie at redshift 1.2. According to this assumption, the corresponding 2-10 keV unabsorbed X-ray luminosity upper limits are always greater than 1043 erg s-1 (see Fig. 3), and then no conclusion on their nature can be derived from the X-ray observations.
Figure 3 shows the 2-10 keV unabsorbed luminosity as a
function of the redshift and MIR luminosities. For 59 sources, we
have only upper limits on the X-ray fluxes and, then, on the un-absorbed X-ray luminosities (using the measured redshifts or, when
not available, a conservative estimate of z=1.2). However these
estimates do not take the effects of possible X-ray absorption into
account . This effect is more relevant at low redshifts (z<0.5-1)
because the absorption mainly affects the rest-frame, lower-energy
photons that are instead no longer observable in the 2-10 keV band at
higher redshift. In Fig. 3 (left) the expected observed X-ray
luminosity is shown as a function of absorption and the redshift of an
absorbed source having an intrinsic X-ray luminosity of 1042 erg s-1. We can expect 50% less emission at redshift
for
a Compton-thin source with a column density of
cm-2. In this case a
1041.7 erg s-1 absorbed
luminosity would be observed. The only solution for detecting such low
redshift-absorbed AGN would be to obtain deeper X-ray observations in
order to analyze their X-ray spectra and separate X-ray luminous
star-forming galaxies from absorbed AGN. As expected, the situation is
even worst for Compton-thick (
cm-2), absorbed AGN:
in this case the fluxes would be reduced by more than a factor of
ten. We thus conclude that our X-ray observations are only able to
detect a relevant, but nevertheless incomplete, fraction of the
Compton-thin AGN. Then our attempt to identify the population of AGN
missed by the optical spectroscopic identifications is partly
inefficient, limited to un-absorbed, or Compton-thin, absorbed AGN.
Therefore, the densities of AGN derived in the next sections are only
to be considered lower limits.
![]() |
Figure 3:
Left. 2-10 keV X-ray
luminosity versus redshift. Squares are AGN1, triangles are
optically classified AGN2, and circles are galaxies. Filled
symbols are X-ray detected sources. The dashed lines show our AGN
luminosity threshold, and the expected observed X-ray luminosity as a
function of absorption and redshift of an absorbed source having
an intrinsic X-ray luminosity of 1042 erg s-1.
Right. 2-10 keV X-ray luminosity versus 15 ![]() |
Open with DEXTER |
![]() |
Figure 4: Optical spectra of the seven objects showing AGN X-ray luminosity and previously classified as no-AGN galaxies. The position of the typical emission lines of AGN (even if not observed) are shown. |
Open with DEXTER |
We observed 72 ELAIS sources with XMM-Newton. Nine sources, previously
classified in the optical as AGN (6 AGN1 and 3 AGN2) showed X-ray
properties (either detections or upper limits) in agreement with their
optical classification. Eight out of the 63 previously classified (or
bona fide) starburst galaxies instead showed X-ray luminosities
typical of AGN activity. As our optical spectroscopic studies were
sensitive enough to detect broad emission lines (see La Franca et al. 2004), we classified them as new AGN2. As discussed in the
previous section, the object C15_J003315-432829 has no spectroscopic
identification, but its X-ray flux is so bright that will have an X-ray
luminosity higher than 1042 erg s-1 at redshifts higher than
z=0.2. We inspected the optical spectra of the remaining seven new
AGN2 (see Fig. 4). C15_J003317-431706 was observed
after the publication of the spectroscopic catalogue by La Franca et al. (2004) and is an R=24.3 galaxy at z=0.689 observed by VIMOS at VLT
(see Sect. 2.1). It has a noisy spectrum, just showing a clear
[OII]3727 emission line, which hampers a more detailed
classification. According to the classification scheme used by La
Franca et al. (2004; see also Dressler et al. 1999; Poggianti et al. 1999), the remaining six galaxies are: two e(a) dust-enshrouded
starburst galaxies, three e(c) star-forming spirals, and one k(e)
post-starburst galaxy with emission lines. C15_J003346-431942, which
was classified as an e(a) galaxy, has the line ratios typical of a
LINER, but could be compatible with an AGN2 classification if 10%
errors are taken into account. C15_J003318-431659, which was
classified as an e(c) galaxy (showing no detectable H
emission
or absorption) has a limit on the [OIII]5007/H
ratio which is
compatible with an AGN2. On the contrary, even taking the line ratio
uncertainties into account, the remaining five objects do not show any
signature in the optical spectrum compatible with an AGN activity.
As shown in Fig. 5, according to the large statistical
uncertainties, there is no significant dependence on the 15 flux by the fraction of the newly X-ray classified AGN2. We can then
conclude that at least a constant fraction of 12.7
+6.3-4.4%
of the previously optically-classified starburst galaxies are instead
AGN2.
According to our X-ray luminosities' upper limits or measurements, all
galaxies with
L15 < 109.5
are unlikely to harbor
an AGN (see Fig. 3). In contrast, all the newly classified
AGN2 have IR luminosities in the range 109.5
,
exactly in the region of the AGN2 LF
where a flattening has been observed (Matute et al. 2006). This result
confirms the hypothesis of Matute et al. (2006) that the observed
flattening of the AGN2 LF at redshift
0.35 was due to an
incompleteness in the AGN identification. As previously pointed out by
La Franca et al. (2004), the origin of this incompleteness has to be
ascribed to the fact that a relevant fraction of AGN2 do not show
clear AGN signatures in their optical spectra.
![]() |
Figure 5:
Up. Flux distribution of the 63 ELAIS-S1-5 sources that were
not optically classified as AGN. The shaded histogram shows the flux
distribution of the eight sources that showed AGN activity in the
X-ray band. Down. Fraction of the X-ray active AGN among the
galaxies as a function of the 15 ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
Revised AGN1 and (Compton thin) AGN2 integral counts at
15 ![]() |
Open with DEXTER |
The more IR luminous objects in our sample are AGN1, which show
15 m luminosities higher than 1011.5
.
As discussed
by Matute et al. (2006), the absence of AGN2 in this luminosity range
is mainly due to selection effects. The AGN1 SEDs have a fairly
constant infrared to optical ratio, while for AGN2 the
infrared-to-optical ratio increases with the increase in the
luminosity (Log (L15/LR)
0.5 Log L15). As a
consequence, AGN2 have on average fainter optical counterparts than
AGN1 at high luminosities and then miss the spectroscopic
identification (see Fig. 1 in Matute et al. 2006).
According to this above new classification, we have then corrected the
integral 15
counts, previously derived from the optical
spectroscopic classification of the ELAIS-S1 sources by La Franca et al. (2004). According to the estimates described in the previous
section, we have assumed that 12.7% of all ELAIS-S1 optically
classified or bona fide starburst galaxies harbor an AGN2. The results
are presented in Table 3 and shown in Fig. 6. The newly
derived counts are compared in Fig. 6 (left) with the
predictions from the AGN1 and AGN2 luminosity functions computed by
Matute et al. (2006), and in Fig. 6 (right) with the AGN
counts predicted by Silva et al. (2004) on the basis of the X-ray AGN
luminosity function (LF) and semi-empirical SEDs (linking the X-ray
to the infrared). It turns out that the AGN1 counts agree, within the
uncertainties, with the counts predicted by both the AGN1 LF of Matute
et al. (2006) and the un-absorbed (
cm-2)
AGN of Silva et al. (2004). Conversely, the new AGN2 counts are, as
expected, about a factor of 2-3 higher than predicted by the AGN2 LF
of Matute et al. (2006). At fluxes larger than 0.6 mJy, the density of
AGN2 is about 60 deg-2, while Matute et al. (2006) predict a
density of about 20 deg-2. On the other hand, the predicted
count of absorbed, Compton-thin AGN (10
cm-2) by Silva et al. (2004) are consistent with our estimates of
the AGN2 counts.
All the extragalactic sources, with 15
fluxes brighter than
0.6 mJy, produce a 15
background of
(
) = 0.343nW m-2 sm-1, while the AGNs produce
(
) = 0.046nW m-2 sm-1 (
13% contribution to the extragalactic
background produced at fluxes brighter than 0.6 mJy).
Table 3:
AGN1 and AGN2 integral counts at 15 m.
We have used XMM-Newton observations of the 15 m sample of the
ELAIS-S1-5 area to disentangle AGN activity among previously
spectroscopically classified starburst galaxies. We find that at least
a fraction of
% (not dependent on the flux) of the previously
spectroscopically classified starburst galaxies are instead bona fide
(Compton thin) AGN2. The fraction of newly identified AGN2 is not
large enough to significantly change the results previously found by
Pozzi et al. (2004) for the evolution of the LF of the starburst
galaxies obtained by using the ELAIS sample. On the other hand, as
the density of galaxies is about an order of magnitude higher than the
density of AGN at 15
m fluxes of
1 mJy, the new density of
AGN2 is significantly increased by a factor of
2-3. This result
agrees with the analysis of Matute et al. (2006), who ascribe the
observed flattening of the AGN2 LF at redshifts
0.5 to an
incompleteness of the optical spectroscopic identification techniques.
Here we wish to focus on the two results: a) that more than half of
the AGN2 could miss their classification if only the optical spectra
were used and b) that, according to our analysis and predictions
derived from the hard X-ray luminosity function, the surface density
of (Compton thin) AGN2 at fluxes larger than 0.6 mJy is in the range
50-70 deg-2. These are the first estimates of separated AGN1
and AGN2 counts in the MIR. This result corresponds to a fraction of
about 24% (17/72) of AGN (AGN1 plus AGN2) among all the extragalactic
sources at 15 m fluxes larger than 0.6 mJy (
13%
contribution to the extragalactic background produced at fluxes
brighter than 0.6 mJy).
According to our analysis, the ongoing studies of the density and evolution of absorbed/reddened AGN (AGN2) among the Spitzer sources, if only based on optical spectroscopic identifications, would miss about half of the true AGN2. A more efficient approach would require MIR color-selection criteria, or better if complemented with deep hard X-ray observations (see e.g. Barmby et al. 2006; Donley et al. 2007; Fiore et al. 2007; Daddi et al. 2007).
Acknowledgements
We thank Laura Silva for discussions and for having provided the predictions of the 15m AGN counts in machine-readable format. We thank the anonymous referee for the useful comments. This research has been partially supported by ASI, INAF, and MIUR grants.
Table 1: X-ray detected sources.
Table 2:
X-ray properties of the 15
sources.