Issue |
A&A
Volume 507, Number 3, December I 2009
|
|
---|---|---|
Page(s) | 1277 - 1289 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/200912261 | |
Published online | 01 October 2009 |
A&A 507, 1277-1289 (2009)
Black hole growth and starburst activity at z = 0.6-4 in the Chandra Deep Field South
Host galaxies properties of obscured AGN![[*]](/icons/foot_motif.png)
M. Brusa1 - F. Fiore2 - P. Santini2 - A. Grazian2 - A. Comastri3 - G. Zamorani3 - G. Hasinger1,4 - A. Merloni1,5 - F. Civano6 - A. Fontana2 - V. Mainieri7
1 - Max-Planck-Institute for extraterrestrial Physik, Giessenbachstrasse, 1, 85748 Garching bei München, Germany
2 -
INAF - Osservatorio Astronomico di Roma, via Frascati 33, Monteporzio-Catone (Roma) 00040, Italy
3 -
INAF - Osservatorio Astronomico di Bologna, via Ranzani 1 40127 Bologna, Italy
4 -
Max-Planck-Institute for Plasma Physics, Boltzmannstrasse 2, 85748 Garching bei München, Germany
5 -
Excellence Cluster Universe, Boltzmannstrasse 2, 85748 Garching bei München, Germany
6 -
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
7 -
European Southern Observatory, Karl-Schwarzschild-str. 2, 85748 Garching bei München, Germany
Received 2 April 2009 / Accepted 22 September 2009
Abstract
Aims. The co-evolution of host galaxies and the active black
holes which reside in their centre is one of the most important topics
in modern observational cosmology. Here we present a study of the
properties of obscured active galactic nuclei (AGN) detected in the
CDFS 1 Ms observation and their host galaxies.
Methods. We limited the analysis to the MUSIC area, for which deep K-band
observations obtained with ISAAC@VLT are available, ensuring accurate
identifications of the counterparts of the X-ray sources as well as
reliable determination of photometric redshifts and galaxy parameters,
such as stellar masses and star formation rates. In particular, we:
1) refined the X-ray/infrared/optical association of
179 sources in the MUSIC area detected in the Chandra observation;
2) studied the host galaxies observed and rest frame colors and
properties.
Results. We found that X-ray selected (
erg s-1)
AGN show Spitzer colors consistent with both AGN and starburst
dominated infrared continuum; the latter would not have been selected
as AGN from infrared diagnostics. The host galaxies of X-ray selected
obscured AGN are all massive (
)
and, in 50% of the cases, are also actively forming stars (1/SSFR <
)
in dusty environments. The median
value of the active nucleus is between 2% and 10% depending on the assumed
/M* ratio. Finally, we found that the X-ray selected AGN fraction increases with the stellar mass up to a value of
% at z>1 and M*>3
,
a fraction significantly higher than in the local Universe for AGN of similar luminosities.
Key words: galaxies: active - galaxies: starburst - X-rays: galaxies
1 Introduction
Galaxy interactions, and more in general the large scale structure (LSS) galaxy environment, are thought to play a major role in regulating both star-formation and accretion onto nuclear super-massive black holes (SMBHs). This implies that the full understanding of galaxy evolution requires a good knowledge of the SMBH census through cosmic time. In particular, feedback between the central SMBHs in their active phases and the interstellar medium is likely to affect strongly the evolution of their host galaxies.
A short, powerful but highly obscured growth phase of both SMBHs and their host galaxies is predicted by many models for the co-evolution of galaxies and AGN (Silk & Rees 1998; Fabian 1999; Granato et al. 2004; Di Matteo et al. 2005; Menci et al. 2008). This phase ends when strong AGN winds and shocks heat the interstellar medium, blowing away the dust and gas and inhibiting further star-formation in the AGN host galaxies. According to this ``evolutionary sequence'' (e.g. Hopkins et al. 2008), highly obscured AGN should be associated to young galaxies in the process of assembling most of their stellar mass through significant episodes of star formation. On the contrary, unobscured AGN should be associated to galaxies with low or absent episodes of star-formation, given that most of the gas and dust responsible for the star formation has been blown away by the effect of AGN feedback. Many observational evidences (see Alexander et al. 2005; Page et al. 2004; Stevens et al. 2006) and theoretical arguments (Menci et al. 2008, and references therein) in favor of the evolutionary sequence do exist. These results challenge our 20-years old AGN view, in which the differences we see in different classes of sources - especially between ``obscured'' and ``unobscured'' ones - are simply due to orientations effects (Antonucci et al. 1985; Antonucci 2003; Urry & Padovani 1995).
A correct and complete identification of obscured and highly obscured AGN at high redshift is therefore crucial because they may represent the first, still little explored, phase of the common growth of both SMBHs and their host galaxies.
Moderately and highly obscured AGN at ,
where most of the accretion and star-formation processes are on-going,
offer an ideal tool for a direct test of feedback models. In these
objects the nuclear light is completely blocked or strongly reduced and
does not overshine the galaxy
optical and near-infrared light. This gives us the possibility to study
host galaxy morphologies,
colors and spectral energy distributions without the difficulty of
disentangling star-light from nuclear light. We can therefore study
galaxy properties during active phases, e.g. when nuclear feedback should be in action.
Obscured AGN at cosmological distances are however usually faint in the
observed optical bands, because the UV rest frame is strongly reduced
by dust extinction. On the other hand, moderately obscured AGN
(or Compton thin, 1022 < < 1024 cm-2) are common in
X-ray images (Bauer et al. 2004; Comastri & Brusa 2008, and reference therein), making up
50% of the full X-ray population at fluxes
<10-14 erg cm-2 s-1 (Gilli et al. 2007) once the contribution from the normal starforming galaxy population is removed
(Ranalli et al. 2005). Compton thick AGN (CT,
cm-2)
are faint also in the X-ray band, because photoelectric absorption and
Compton scattering strongly reduce the X-ray flux up to
>10 keV, and over the entire X-ray range if
cm-2.
Indeed, only a dozen CT AGN is present in the deepest images of
the X-ray sky, the Chandra Deep Fields (Norman et al. 2002; Tozzi et al. 2006; Georgantopoulos et al. 2008).
Identification of the correct counterparts of obscured AGN at
cosmological redshift is not trivial. These objects are faint in the
optical band, because 1) the intrinsic AGN emission is
absorbed by the surrounding material; and 2) the host galaxy light
is strongly reduced by cosmological dimming (L* galaxies at z=1-2 have
).
The probability to find by chance a galaxy of these magnitudes in the
Chandra error boxes is not negligible. The identification process is
made easier by using deep near infrared images, because the surface
density of these sources is smaller than that of optical ones, bringing
the probability of finding a galaxy by chance in the Chandra error
boxes to comfortably small values. Moreover, the K-band
flux is more tightly correlated with the X-ray flux than the optical
(obscured) one. For this reason, we have reanalyzed the identifications
of the Chandra Deep Field South (CDFS) 1 Ms sources (Alexander
et al. 2003) in the area covered by sensitive K band and IRAC 3.6
m and 4.5
m
observations (the GOODS MUSIC area, Grazian et al. 2006).
The excellent multiband photometry available in this area allows the determination of a reliable photometric redshift for the faint sources not reachable by optical spectroscopy. Once identified the counterparts of the X-ray sources and their redshift, it is possible to proceed to a detailed study of the properties of their host galaxies. We present here a study of the properties of the host galaxies of X-ray obscured AGN, such as their morphology and close environment, optical and infrared rest frame colors, mid infrared to optical spectral energy distributions. Galaxy masses and star-formation rates are derived from the fit of galaxy templates to the optical to near infrared spectral energy distributions. We finally compare the masses and star-formation rates of obscured AGN host galaxies to those of inactive galaxies in the field selected in both optical and infrared bands.
The paper is organized as follows: Sect. 2 presents the
CDFS datasets, the X-ray to optical/infrared association and the
obscured AGN sample. Section 3 presents the observed frame
colors of the obscured AGN sample. Section 4 presents the host
galaxies properties (masses and
star-formation rates, SFRs) of the sample and our estimates of the
fraction of AGN in mass selected samples. Section 5 presents a
discussion of the results and Sect. 6 outlines a summary of the
most important points. A cosmology with H0=70 km s-1 Mpc-1,
= 0.3,
is adopted throughout (Spergel et al. 2003).
2 Obscured AGN selection in the CDFS
2.1 X-ray data: catalogs and multiband identifications
X-ray data for a total exposure time of 1 Ms
in the GOODS field have been obtained by the CDFS consortium and
are publicly available. The source catalog has been published by
Giacconi et al. (2002),
including the basic X-ray properties in three different X-ray bands
(source counts, fluxes, exposure times). Subsequently, Alexander
et al. (2003, hereinafter
A03) published a re-analysis of the CDFS data that led mainly to
an improved astrometry of the X-ray
positions with some differences in the number of detected sources at
the very faint counts level (see A03 for more details).
Optical identifications for the sources detected by Giacconi et al. (2002), including the position of the candidate counterparts, optical magnitudes and the spectroscopic redshifts from a dedicated VLT campaign, are reported in Szokoly et al. (2004). Photometric redshifts based on good quality optical photometry (VLT and ACS/HST) and deep infrared (ISAAC/VLT) observations are available for all but four CDFS sources without spectroscopic redshift (Zheng et al. 2004; Mainieri et al. 2005), leading to a virtually complete catalog of identified X-ray sources. Using this catalog and redshifts information, Tozzi et al. (2006) derived the rest-frame X-ray properties for the 1 Ms CDFS sources. IRAC photometry, photometric redshifts and radio properties of different subsamples of the CDFS X-ray sources have been published in others works appeared in the literature in the past few years, e.g. Rigby et al. (2005), Alonso-Herrero et al. (2006), Georgantopoulos et al. (2007), Rovilos & Georgantopoulos (2007) and Tozzi et al. (2009).
These published works rely on the association of the X-ray sources to a candidate optical counterpart. However, given that the optical flux of the X-ray counterparts can be strongly suppressed by dust absorption, an identification process that uses as a prior a catalog generated at a redder wavelength (K-band or IRAC bands) can be more efficient in isolating the correct association to the X-ray source (see, e.g. Brusa et al. 2007). The lower surface density of IRAC sources with respect to faint optical objects is widely used as an argument to indicate that, in case of ambiguous association between two sources, the most likely candidate is the one detected at the brightest IRAC fluxes (whatever the optical flux).
In order to quantify how (and if) the use of catalogs selected
in different bands leads to significantly different X-ray to optical
associations, we decided to perform again the identification of the
Chandra sources, using a new version of the multi-color GOODS-MUSIC
sample (GOODS MUlticolor Southern Infrared Catalog; Grazian et al.
2006; De Santis et al. 2007), extracted from the public data of the GOODS-South survey (Giavalisco
et al. 2004), and updated as described in Santini et al. (2009). The 15-bands multi-wavelength coverage ranges from 0.35 to 24 m,
as a result of the combination of images coming from different
instruments (2.2ESO, VIMOS/VLT, ACS/HST, ISAAC/VLT, IRAC/Spitzer and
MIPS/Spitzer). In the following we will refer to the catalog with full
coverage in all bands. Such a catalog covers an area of
143.2 arcmin2
and contains 15 208 sources. After culling Galactic stars,
it contains 14 999, extragalactic objects selected in the z band, in the Ks band and
at 4.5
m (see Santini et al. 2009, for further details).
A total of 183 X-ray sources in the A03 catalog are found within
the GOODS-MUSIC area. Of these, 3 objects are associated with
stars (ID #120,279,289 in A03) and an additional source is most
likely a Ultra Lumimous X-ray source (ULX) in a nearby
galaxy (ID #213). Therefore, the number of extragalactic
sources (AGN or galaxies) in the MUSIC area is 179, 154 in
common with Giacconi et al. (2002) and an additional 25 detected only by A03.
167 of the 179 AO3 sources in the MUSIC area are in common with the catalog obtained from the full 2 Ms exposure (Luo et al. 2008). We have visually inspected and carefully checked all the 12 sources missing in Luo et al. (2008). One of these sources, #225 is most likely a spurious detection on the wings of the Chandra PSF of a brighter source (#224). Other three sources, #183, #215 and #239, lie in problematic areas, close to a CCD gap and are probably spurious detection in AO3. These four sources are flagged as likely spurious in Table 1. The other 8 sources present in the AO3 catalog but not in the Luo et al. (2008) catalog are probably real sources, which have not been detected in the full 2 Ms exposure because of variability or a more conservative detection threshold in Luo et al. (2008), or both.
We proceeded for the association of the counterparts to the
X-ray sources as follows. First, we used a statistical method
(the ``likelihood ratio technique'', following the approach
described by
Brusa et al. (2007) to isolate the most likely z-band, K-band and IRAC counterparts from the z-, K-
and 4.5 micron selected catalogs in the MUSIC area. The
method calculates the probability that a source is the correct
association by weighing the information on the X-ray to optical
distance, the surface density of (possible) false coincidence
background objects and the brightness of the chosen counterpart. With
this method we were able to isolate the most obvious associations, and,
at the same time to pick up the objects with 2 or more different
counterparts with comparable likelihood in the same or different bands ( of
the entire sample). Then, we visually inspected all the proposed
179 associations in order to resolve possible ambiguous cases.
The result of our work is summarized in Table 1.
For each X-ray source we report the X-ray identificator from A03 (AID,
Col. 1), the X-ray identificator from G02/S04 (XID, Col. 2),
the MUSIC identificator (OID, Col. 3), the positions of the
proposed counterpart (Cols. 4, 5), the distance of the chosen
counterpart from the X-ray centroid (Col. 6), a flag for the
identification (Col. 7), the observed X-ray fluxes in the hard
(2-10 keV) and soft (0.5-2 keV) bands (Cols. 8
and 9), the redshift of the source (Col. 10, see next
section) and the rest frame 2-10 keV luminosity
(Col. 11). The X-ray fluxes (in erg cm-2 s-1) have been derived from the counts reported in A03, using the following conversion factors:


for a power law spectrum with energy index equal to 0.4, i.e. the average value for deep fields sources (see, e.g., Tozzi et al. 2001)
![[*]](/icons/foot_motif.png)
Of the 154 published counterparts (first and second part of the table, mostly from Szokoly et al. 2004, but see also notes in Table 1), we were able to confirm 144 X-ray to optical associations (93.5%, sources with flag = 0), we changed 9 X-ray to optical associations (5.8%, flag = 1)), while for the remaining case (AID 186) we were not able to assign a reliable IRAC or optical counterpart (0.6%, i.e. neither to confirm the published counterpart nor to propose a new counterpart). The fraction of previously misidentified sources rises up to 15% when the optically faint (z>24) sample is considered (9/61).
Figure 1 shows the Chandra, ACS (z-band), K and IRAC-4.5 micron cutouts for the 9 sources we identify with a counterpart different from that proposed in the literature (see again notes in Table 1). The cyan circle marks the new counterpart, while the yellow circle marks the old counterpart. In almost all the cases the new proposed counterpart is absent in the optical images.
For the remaining 25 sources detected only in A03, we were able to unambiguously identify the counterpart for 23 of them. The optical counterparts of these 23 sources are reported in the third part of Table 1.
![]() |
Figure 1:
Chandra, ACS (z-band), K-band and IRAC-4.5 micron cutouts (5'' |
Open with DEXTER |
![]() |
Figure 2: Redshift distribution of the X-ray sources. Open histogram: all sources; right diagonal (blue) shaded histogram: spectroscopic sample; left diagonal (red) shaded histogram: photometric sample. |
Open with DEXTER |
We computed the rest frame, absorption corrected 2-10 keV luminosities from the 2-10 keV
fluxes, when available, and from the 0.5-2 keV fluxes in the remaining cases. We assumed a power law spectrum with
reduced at low energy by rest frame photoelectric
absorption. We used the
values reported by Tozzi et al. (2006) for the sources with redshift consistent between this and the previous work (
% of the sources have
,
see also next section). For the remaining 32% of the sources we evaluated
from the hardness ratio between the hard and soft bands.
2.2 Optical+Infrared photometry and photometric redshifts
Table 2 shows the compilation of optical b,v,i,z bands (Cols. 2-5) from GOODS data (Santini et al. 2009), infrared J,H,K
(Cols. 6-8) from ESO/VLT, IRAC and MIPS (Cols. 9-13)
photometry. Of the 179 X-ray sources, 174 are detected in the
z-band, 172 in the 4.5 m IRAC channel and 133 in the 24
m MIPS survey.
The spectroscopic or photometric redshifts available for the proposed optical counterparts
of the 179 CDFS-GOODS sources are listed in Col. 10 of Table 1.
When a spectroscopic redshift is present, the relative source catalog
is reported, as follows: (1) CXO-CDFS (S04); (2) K20 (Mignoli
et al. 2005); (3) GOODS (Vanzella et al. 2005, 2006, 2008; Popesso et al. 2008); (4) COMBO-17 (Wolf et al. 2003); (5) VVDS (Le Fevre et al. 2004); (6) other (Daddi et al. 2005; Cristiani et al. 2000). Otherwise, the redshift is from SED fitting (Santini et al. 2009) and is reported
with its 1 confidence
interval. We have spectroscopic or photometric redshifts for all but
one (AID 282) of the identified X-ray sources. However, in eleven
cases, only a lower limit on the
photometric redshift solution is available from the SED fitting,
i.e. at higher redshifts the
curve is relatively flat and an upper limit on the basis of a
could not be
determined.
Figure 2 shows the redshift distribution of the entire sample (open black histogram) and the distributions derived for the spectroscopic and photometric sub-samples (blue and red histograms, respectively). About half of the redshifts at z>1 come from the photometric sample. For about 25% of the sources the redshifts listed in Table 2 are different (i.e. not consistent within the errors) from the ones published in previous works (as compiled by Zheng et al. 2004; Tozzi et al. 2006), for the following reasons: 1) for 9 objects the counterparts have changed from the ones reported in the literature (see previous section); 2) for about 30 sources the availability of longer wavelength (IRAC) data increases the reliability of the photometric redshifts. The photometric redshifts for the sources in the MUSIC sample are instead consistent with a new analysis of the sources in the 2 Ms CDFS (Luo et al. 2008) and the Extended Chandra Deep Field South (E-CDFS, Lehmer et al. 2005) surveys obtained by Rafferty et al. (2009, in preparation).
2.3 Obscured AGN sample
We extracted from our sample of 179 sources, the sample of ``obscured''
AGN defined as outlined below. First, we considered all the sources
with rest frame (2-10 keV) >
41.8 in order to minimize the contribution of non-AGN objects among the
X-ray selected sample. Indeed, following Ranalli et al. (2003), the observed X-ray emission can be due to stellar processes only in objects with SFR larger than
100
yr-1
(see also Sects. 4.2 and 4.4 for further discussion on the
luminosity threshold). We also removed one source that shows extended
X-ray emission, and the three sources without a unique/secure
identification (see Sect. 2.1). We were left with
142 sources. Secondly, we isolated all the 19 sources
spectroscopically identified as Broad Line (BL) AGN in Szokoly
et al. (2004), and additional 7 objects still spectroscopically unidentified but with a sg parameter
(stellarity, from SExtractor, Bertin & Arnouts 1996) larger
than 0.9. We named these 26 sources as BL/pointlike sample.
Following Silverman et al. (2008), we computed the ratio of the flux in the z850-band
in two fixed apertures (a small aperture with a radius of
3 pixels (0.09'') and a larger aperture with a radius of
25 pixels (0.75'')). This quantity is similar to the inverse
of the concentration index (Abraham et al. 1994). Fluxes were
determined by extracting counts from background-subtracted cutouts of
ACS images centered aroud the position of the X-ray sources.
Figure 3 shows the distribution of the z-band
counts ratio for the BL/pointlike sample (blue histogram), and the
remaining sources (116 objects, red histogram). All the objects
with ratio values >0.3 are from the BL/pointlike sample, as expected
in sources for which the
pointlike and/or nuclear emission dominates overall the host-galaxy
light. The contribution of pointlike/unobscured AGN at ratio values
<0.3 is
% (being zero in the two lowest ratio bins, 0-0.1), therefore we expect a negligible contamination from this class of
sources in the sample with ratio <0.3. We call ``obscured'' AGN the 116 objects with the ratio of z-band counts <0.3 and not classified as BL AGN or pointlike in the ACS images. Ninety of these ``obscured'' AGN (
%) are also detected at MIPS 24
m.
![]() |
Figure 3: Distribution of the z-band count ratio for the ``obscured'' AGN sample (116 objects, red long dashed histogram) and the BL/pointlike sample (26 objects, blue dotted histogram). |
Open with DEXTER |
X-ray column densities from the 1 Ms sources in the CDFS derived
from an accurate spectral fitting analysis have been published by Tozzi
et al. (2006). As discussed in the previous section, a fraction of the X-ray to optical associations changed from the Tozzi et al. (2006)
analysis to the present analysis (9 sources), and in addition some
more robust photometric redshifts have been obtained using the Spitzer
photometry. Among the sixty-nine sources in the ``obscured''
AGN sample in common with the Tozzi et al. sample and with a
from the redshift
values reported in Tozzi et al. (2006), fifty-four (
%) show X-ray absorption in excess of 1022 cm-2, further suggesting that our sample is representative of the X-ray obscured AGN population.
We visually inspected all the sources in the z-band ACS images and divided them in three main classes: isolated sources (35), interacting/disturbed/clumpy galaxies (53) and objects too faint to tell something on the morphology (26). All of these 26 sources populate the first bins in the histogram in Fig. 3.
![]() |
Figure 4: ACS (z-band) images of obscured AGN counterparts with disturbed morphology or interacting systems (AID 88, 107, 226 and 84). |
Open with DEXTER |
Figure 4 shows four examples of obscured AGN counterparts with disturbed morphology or interacting systems.
![]() |
Figure 5:
Left panel: 24 |
Open with DEXTER |
2.4 Host galaxy properties
For the host galaxies of the obscured AGN we retrieved from the GOODS-MUSIC catalog the rest frame U-V colors, the absolute V and K-band magnitudes, the SFRs and the stellar masses M*
derived from the SED fitting. More specifically, stellar masses
are derived from the comparison of SED and libraries of synthetic
spectra, assuming a Salpeter (1955) Initial Mass Function (IMF) and simple stellar population (SPP) models (Bruzual & Charlot 2003), as described at length in Fontana et al. (2006). SFRs have been estimated from SED fitting, using Bruzual & Charlot (2003)
synthetic models and fitting the UV rest frame photometric bands.
We parametrize the star formation histories with a variety of
exponentially declining
laws (with timescales ranging from 0.1 to 15 Gyr),
metallicities (from Z = 0.02
to Z = 2.5
)
and dust extinctions (0 < E(B-V) <
1.1, with a Calzetti or Small Magellanic extinction curve). These SFR
estimates have been compared with the ones derived from the
observed IR (MIPS) emission, following the approach presented by
Papovich et al. (2007):
despite the majority of galaxies are assigned a consistent SFR (with
overall median ratio around unity, see details in Santini et al. 2009),
in the following, we will use the SFR from the SED fitting for all the
obscured AGN in our sample, in order to minimize the possible
AGN contribution in the computation of SFRs from MIPS data in
AGN samples (see Daddi et al. 2007a).
As a word of caution, we note here that the derivation of host galaxy properties is affected by both statistical and systematics undertainties.
Statistical errors for both masses and SFR have been calculated from the 1 convidence level in each parameter estimate by scanning the
levels, allowing the redshift to change within the errors in case of objects with photometric redshifts
(see Fontana et al. 2006; Santini et al. 2009).
The masses and SFR estimates are, however, limited by the systematics
related to model degeneracy and, in particular, the ones associated to
the adopted IMF, the SSP, the metallicity and the star formation
history (see Santini et al. 2009; and the discussion in Ilbert et al. 2009).
The largest, systematic, uncertainty on the estimate of the stellar
mass is the assumption on the IMF. The masses obtained using a
Chabrier IMF are on average 0.24 dex lower than those
obtained adopting a Salpeter IMF, and this effect is totally
independent from the range of redshift and mass probed (Salimbeni
et al. 2009). The uncertainties related to the adopted
SSP models (Bruzual & Charlot 2003; Charlot & Bruzual 2007; Maraston 2005) to construct
libraries of synthetic spectra, are also of the order of
= 0.2-0.3, but the effect is less evident in the high-mass part and at high-z (see Salimbeni et al. 2009, for a proper discussion on these effects in the GOODS-MUSIC sample).
Similar systematics effects are present also in the SFR estimates, and are expected to be of the order of
)
0.2-0.4 (see, e.g., discussion in Ferreras et al. 2005, for the variation due to different IMF).
Since most of these effects are systematics and since we are interested at general trends over several decades of mass and not in the detailed description of the mass of the host galaxies of each X-ray source, we are confident that the overall results discussed in the rest of the paper remain unchanged.
3 Observed frame colors
A proper classification of source nature (AGN vs. starburst vs. normal evolved galaxy) should ideally be obtained via a complete analysis of its emission over the electromagnetic spectrum, using both spectroscopic (e.g. emission line widths and ratios) and photometric (all-wavelengths sampled SED) observables. However, in reality, a complete source characterization is not obtainable, even in fields where the best and deepest imaging and spectroscopic campaigns have been performed, such as the GOODS/MUSIC field. Partial, but still reliable, information (at least in a statistical sense) can be obtained through the analysis of the emission in bands where differences between nuclear and star formation emission are emphasized, comparing the observed SEDs with galaxy templates over a range as broad as possible. In particular, the combination of observed-frame mid-infrared, near infrared and optical flux ratios has been exploited recently in the literature to isolate obscured AGN (see, e.g., Martinez-Sansigre et al. 2005, 2006; Fiore et al. 2008; Donley et al. 2008; Dey et al. 2008).
Figure 5, left panel shows the 24 m to optical (I-band) flux ratio vs. the I-[3.6] color
for the X-ray sources in our sample. The green contours show the
isodensities of the field objects in the MIPS selected sample. The
sizes and shape (colors) of the symbols reflect the
AGN 2-10 keV
luminosity: (blue) small triangles are objects with
,
(red) large squares
objects with
.
Open circles mark all the other X-ray sources (both BL/pointlike AGN and sources with
<
41.8). The well-established correlation between the
MIR/O color and the total source luminosity traced in this case by
the X-ray emission (see Fiore et al. 2008)
is also present in this plot, where the most luminous sources (red
squares)
cluster at average values of MIR/O ratio higher than those of less
luminous sources (blue triangles), and on average towards redder
optical to near infrared colors.
The right panel of Fig. 5
shows an IRAC color-color diagram where the [3.6]-[4.5] color is
plotted against the [5.8]-[8.0] color. The dotted lines isolate a
wedge for the selection of bright, optically unobscured quasars,
originally proposed in Stern et al. (2005). The
sizes and shapes (colors) of the symbols reflect the AGN X-ray luminosity, as in the left panel of Fig. 5.
Open small symbols are unobscured AGN and all the other X-ray sources.
While the most luminous ``obscured'' AGN occupy the region in the
color-color diagram expected from a power-law SED, i.e. within the
wedge marked by the dashed lines, the moderate luminous ``obscured''
AGN population mainly clusters in two regions with roughly the
same [5.8]-[8.0] color but with a bimodal distribution in the
[3.6]-[4.5] color (see also results presented by Cardamone
et al. 2008). About half of the moderate luminous obscured AGN in
the sample are clustered around [3.6]-[4.5]
-0.25. Most of these objects are in
the redshift range z = 0.6-0.8 (i.e. in the structures present in the CDFS field, see Gilli et al. 2005; Silverman et al. 2008).
In this region also fall the 60% (4/7) of the radio detected
sources in the obscured AGN sample, extracted from the work by
Tozzi et al. (2009, see also the results presented by Hickox et al. 2009). The other half of the obscured AGN sample clusters at redder [3.6]-[4.5] colors (
), where both
star-forming and normal galaxies at z
1-3 are expected. The isodensity curves from the GOODS-MUSIC sample (Santini et al. 2009), overlaid as green contours in the plot, show the same bimodal distribution as observed for the X-ray sources.
In Fig. 5, right panel, we also plot as magenta asterisks the population of the most obscured, possibly Compton Thick AGN at
discussed by Fiore et al. (2008) and selected on the basis of their high 24
m to optical flux ratio and extremely red colors (R-K > 4.5; see also Daddi et al. 2007a; Dey et al. 2008). These X-ray undetected, Compton Thick AGN candidates occupy the same region in which X-ray detected, z>1 obscured AGN lie (see also Georgantopoulos et al. 2008).
![]() |
Figure 6:
Spitzer color color diagram used to isolate star-formation and AGN
dominated high-z galaxies (Pope et al. 2008). The small green points represent the distribution of field galaxies in the 24 |
Open with DEXTER |
Figure 6 shows the combined MIPS-IRAC color-color diagram introduced by Pope et al. (2008; see also Ivison et al. 2004) to distinguish between AGN and star-formation dominated sources in high-z (z > 1.5) infrared selected sources, where the 24 m to 8.0
m flux ratio is plotted against the 8.0
m to 4.5
m flux ratio. Indeed, at z>1.5, the observed 24
m emission samples the rest frame
m
emission, where strong starburst PAH lines are commonly present
and can boost the measured IR flux, giving a higher 24
m to 8.0
m flux ratio. Similarly, the 8.0
m to 4.5
m flux ratio is a measure of the rest-frame near infrared continuum slope. Pope et al. (2008) proposed that dust obscured galaxies with large values of the 8.0
m to 4.5
m flux ratio (>2; dotted line in Fig. 6)
are AGN dominated sources, according to the presence of
AGN continuum in Spitzer IRS spectra, while objects with
lower values of the 8.0
m to 4.5
m flux ratio (<2) are ``star-formation dominated'', in agreement with the fact that they span a wide range of
24
m to 8.0
m
flux ratio. In both panels of this figure, we plot the field population
extracted from the MUSIC sample as small crosses. In the left
panel, we report as squares the CDFS obscured AGN with a MIPS
24
m detection, and in the redshift range z=1.5-4
(38 objects in total). In the right panel, we mark as triangles
the sources satisfying, among the underlying population, the dust
obscured galaxies selection used by Pope et al. (MIPS/O ratio
larger than 1000 and R-K > 4.5, see also Fiore et al. 2008; Dey et al. 2008)
and as circles the X-ray selected obscured QSO from the
Hellas2XMM survey for which deep Spitzer imaging has been obtained
(Pozzi et al. 2007).
According to the classification introduced by Pope et al. (2008),
the subsample of X-ray selected, ``obscured'' AGN with
MIPS detection show Spitzer colors consistent with both an AGN
dominated continuum and a starburst dominated continuum in the
mid-infrared. Actually, a large fraction of the most X-ray luminous AGN
lie in the Pope et al. (2008)
``star-formation dominated'' part of the diagram, similarly to the
X-ray selected obscured quasars from the Hellas2XMM survey. In
conclusion, for at least half of the obscured AGN with
MIPS detection in our sample, the accretion activity is revealed
thanks to the presence of the X-ray emission and not by a large
8.0 m to 4.5
m
flux ratio. Moreover, about half (13/25) of the sources classified
as AGN in the ``starburst'' part of the diagram have 24
m to 8.0
m flux ratio lower than 5, while in the original Pope et al. diagram this region is almost empty.
4 Rest-frame colors: host galaxy properties
4.1 Reddening
![]() |
Figure 7: Left panel: U-V (rest frame) vs. the absolute V-band magnitude for X-ray obscured AGN (squares) and the underlying optically selected galaxy population (contours). The sizes and colors of the symbols reflect the AGN X-ray luminosity (as in Figs. 5 and 6). Right panel: [3.6]-[4.5] color vs. the U-V (rest frame) for obscured X-ray sources (squares) and the underlying optically selected galaxy population (contours). The upper panels show the comparison of the X-ray selected sources with the optically-selected (z-band) galaxy population in three different redshift bins: z = 0.6-1.0, z = 1-2, and z = 2.4 from left to right, respectively. The lower panels show the comparison for the subsample of objects detected also at 24 micron, in the same redshift bins. The sizes and colors of the symbols reflect the AGN X-ray luminosity as in Figs. 5 and 6. Unobscured AGN are not shown here. |
Open with DEXTER |
The left panel of Fig. 7 shows the U-V (rest) vs. the absolute V-band magnitude for X-ray selected ``obscured'' AGN (squares) and the underlying optically selected galaxy population (contours) in the redshift interval z = 0.6-4.0. The sizes and colors of the symbols reflect the AGN X-ray luminosity (as in Figs. 5 and 6). X-ray selected, obscured AGN show a (U-V) rest color redder than the average color of the optically selected galaxy population and are hosted preferentially by luminous (MV < -20), red (U-V > 0.5) galaxies, populating mostly the red sequence and the green valley of the bimodal color-magnitude distribution observed for field objects. We do not find any trend with the X-ray luminosity: the U-V color does not depend on the strength of the X-ray emission. This result is in agreement with previous findings reported in the literature (Nandra et al. 2007; Rovilos & Georgantoupolos 2007; Silverman et al. 2008; Hasinger 2008); in the next section we will further discuss it and its consequences for our fully understanding of AGN hosts properties.
The right panel of Fig. 7 shows the [3.6]-[4.5] color vs. the U-V
(rest frame) color for the sample of obscured AGN (squares) and the
underlying optically selected galaxy population (contours), in three
different redshift bins:
z=0.6-1.0, z=1-2, and z=2-4 from left to right,
respectively. The upper panels show the comparison of X-ray selected
obscured sources with the optically-selected (z-band) galaxy population; the lower panels show the comparison for the subsample of objects detected also at 24 m.
We do not find a significant evolution in the (U-V) rest frame color with redshift for both the overall galaxy population and the obscured AGN sample. Instead, the average value of the [3.6]-[4.5] color increases with redshift and the average [3.6]-[4.5] color of the obscured AGN is consistent with that of the underlying galaxy population (see also Fig. 5). The only difference among the field and obscured AGN population is, therefore, the fact that obscured AGN are redder, in U-V color, than optically selected sources (as outlined in the left panel), while sharing the same [3.6]-[4.5] color. This is probably due to the fact that the observed [3.6]-[4.5] color strongly correlates with redshift and it is not related to any intrinsic property of the sources (e.g. star formation or AGN activity). If we consider the MIPS selected sample (lower panels) the X-ray selected, obscured AGN closely follow the distributions of the underlying star-forming population with redder U-V rest frame colors than in the z-selected sample.
![]() |
Figure 8: Left panel: U-V (rest frame) vs. the absolute K-band magnitude for the obscured AGN (squares) and the underlying galaxy population (contours). The upper panels show the comparison of the X-ray selected sources with the optically-selected (z-band) galaxy population in three different redshift bins: z=0.6-1.0, z=1-2, and z=2.4 from left to right, respectively. The lower panels show the comparison for the subsample of objects detected also at 24 micron, in the same redshift bins. The shapes and colors of the symbols reflect the AGN X-ray luminosity (see Fig. 5). Right panel: stellar mass M* vs. redshift for field galaxies (crosses) and our ``obscured'' AGN (squares and circles, same notation as in Fig. 5). The vertical lines mark the redshift range we chose for our analysis. |
Open with DEXTER |
![]() |
Figure 9:
(Left panel): the distribution of the SFR for both the optically ( upper) and MIPS ( lower) selected sample of field galaxies (black open histogram) and obscured AGN (red filled histogram) in the redshift interval z=1-4 and in a specific stellar mass range (
|
Open with DEXTER |
Figure 8 (left panel), shows the U-V (rest frame) vs. the absolute K-band magnitude for the obscured AGN in our sample (squares) and the isodensity contours of the galaxy population (green curves), in the same redshift bins as in Fig. 7, and for both the optically (z-band) and MIPS selected samples.
At all redshifts, most of the obscured AGN host galaxies show absolute magnitude brighter than MK<-23, much brighter than that of the bulk of the optically selected population (upper panel): the large majority of optically selected galaxies cluster around lower values of MK and bluer values of the (U-V) rest frame colors.
The distribution of sources in the U-V vs. MK plane is different when the MIPS selected sample is considered: obscured AGN hosts share the same U-V colors
of the mid infrared selected galaxy population. In this case, a color
bimodality is not observed at any redshift, and the U-V color
correlates smoothly with the galaxy mass (also see Aussel et al.
2009, in prep.). AGN hosts are found preferably at the peak
of the galaxy isodensity contours in Fig. 8 lower panels. In particular, the comparison of the upper and lower panels of Fig. 8
suggests a scenario in which star forming galaxies hosted in low mass
systems are not particularly dust obscured and therefore the
reprocessed emission is not revealed at IR bands (at least at
the limiting fluxes of the 24 m GOODS survey), consistent with the results by Daddi et al. (2007b).
Given that MK is a good tracer of the
stellar mass of the galaxies, this suggests that obscured AGN are
preferentially hosted in high-mass systems. This is shown in the right
panel of Fig. 8, where the best fit
is plotted against the redshift (see also Alonso-Herrero et al. 2008)
for both the field population from the MUSIC catalog (crosses) and the
obscured AGN (red circles). As discussed in Sect. 2.4, the
determination of the stellar mass suffers from severe systematic
uncertainties. However, these uncertainties are present for both the
AGN and the underlying galaxy population, with the net results of a
likely shift in the Y-axis (
= -0.3).
4.2 Star formation activity
The black open histogram in Fig. 9
left panel, shows a normalized distribution of the SFR (see
Sect. 2.4) for both the entire MUSIC catalog (upper panel) and the
MIPS selected sample (lower panel). The red filled histogram shows the
distribution of SFR for the ``obscured'' AGN in our sample. In all
cases the SFR have been restricted to objects in the redshift interval z=0.6-4, and in the stellar mass range
.
We chose this stellar mass range because most of the obscured AGN in
our sample have masses larger than this value (see Fig. 8
right panel) and the completeness of both the entire MUSIC catalog and
the MIPS selected sub-sample have been studied in details (see
discussion in Fontana et al. 2009). The total sample considered here includes 103 obscured AGN (82 in the MIPS detected subsample).
In % of the cases the host galaxies of these AGN show substantial (
yr-1)
amount of SFR. The distribution of SFR is significantly different from
the one of optically selected galaxies (the probability that the two
distributions shown in the upper-left panel of Fig. 9 are
drawn from the same parent population is only 0.1%, according to a
Kolmogorov-Smirnov statistics) suggesting an enhancement of star
formation activity in objects hosting obscured AGN. This result is
solid despite the large systematics uncertainties (dominating the SFR
measurements, see discussion in Sect. 2.4), given that they are
expected to be present in both the optically and X-ray selected population.
An enhancement of the SFR activity in obscured AGN is qualitatively in agreement with the results at
obtained from an X-ray selected sample of AGN from the XMM-COSMOS survey (Silverman et al. 2009). As already pointed out in Fig. 8, the distribution of the SFR in the AGN sample mirrors the one of the MIPS selected sample (Fig. 9, lower-left panel) suggesting that AGN activity and star formation are indeed tightly related phenomena.
For each AGN host galaxy we could then calculate the Specific
Star-Formation Rate (SSFR), defined as the instantaneous ratio of SFR
to the total stellar mass (
), and its inverse,
,
the so called ``growth time'', i.e. the time it would take a
galaxy to double
its stellar mass if forming stars at the observed instantaneous rate.
Objects with growth time lower than the age of the Universe at their
redshift (
)
are actively forming stars (see, e.g. Fontana et al. 2009), while objects with growth time larger
than the age of the Universe (
)
can be considered dominated by a quiescent population. We will name
hereafter ``active host galaxies'' the former, and ``inactive host
galaxies'' the latter, when referring to the AGN hosts.
The right panel of Fig. 9
shows the growth time as a function of the stellar mass of the galaxies
in three different redshift bins for the overall field population
(green crosses) and for the obscured AGN sample (filled symbols,
color-coded on the basis of the X-ray luminosity, see Fig. 5).
The horizontal lines mark the age of the Universe at the two redshift
boundaries of the chosen intervals. Taking the estimates of the SSFR at
their face values, the majority (%) of obscured AGN in each redshift interval is hosted in ``active host galaxies'' confirming the trend outlined in Fig. 9
(left panel). However, we should note that in the highest redshift bin,
the SSFR estimates suffer also from large statistical uncertainties
(see error bars in the upper-left part of each panel in Fig. 9,
right), and, as a consequence, the distinction between ``active'' and
``inactive'' host galaxies is rather fuzzy. Therefore, the results for
this bin should be taken only as indicative.
4.3 AGN fraction
In the following we derive the AGN fraction as a function of the
stellar mass of the underlying galaxy population. The AGN fraction
has been computed as the ratio of X-ray selected AGN over the overall
field population in each mass bin. Given that at z>2 we can sample only objects with
erg s-1,
the fractions have been computed further imposing this threshold to the
obscured AGN sample (i.e. considering only objects marked as
squares in the figures
of the paper) in all redshifts bins and, therefore, they should be
considered representative of this AGN population only
.
Figure 10 shows the AGN fraction as a function of M*
in three different redshift bins. We find that the fraction of objects
hosting AGN increases with the stellar mass of the galaxy, being up to
30% for
,
while it is <1% for
M*=1010-3
.
We should note that the MUSIC sample cannot be considered complete for galaxies with masses
already at
.
Indeed, different levels of completeness affect different kind of
galaxies. The total sample is expected to be complete against passive
galaxies with
up to z=4. In particular, the stellar mass limit of star-forming galaxies, which have M/L lower than passive galaxies, is lower than that of passive galaxies at all redshifts (Fontana et al. 2006). The incompleteness of the MUSIC sample for
galaxies implies that the plotted AGN fractions are somewhat
higher than the real fractions at those masses, and that the actual
trends of increasing AGN fraction as a function of the stellar
mass is probably steeper than those in Fig. 10.
![]() |
Figure 10:
Fraction of obscured AGN with
|
Open with DEXTER |
The observed trend of the AGN fraction increasing with the stellar mass
means that, for a given mass-selected sample, AGN will populate the
high-mass tail of the parent sample mass distribution.
Qualitatively similar results are obtained by Bundy et al. (2008) and Yamada et al. (2009). This therefore also explains the position of the AGN in the left panel of Fig. 8, i.e. among the most luminous, massive and reddest sources. Best et al. (2005)
found a similar behaviour for the AGN fraction as a function of
the stellar mass in SDSS AGN samples, selected on the basis of
line ratios diagnostic diagrams (Baldwin et al. 1981; Kauffmann et al. 2003). In Fig. 10 we also plot the results from Best et al. (2005)
corresponding to AGN with L([OIII]) luminosity larger than 108.0 and
,
as labeled. Using the relation between the observed L(OIII) and the X-ray luminosity given by Netzer et al. (2006, Eq. (1)), these L([OIII]) luminosities
represent a fiducial range for local objects with X-ray luminosities
similar to those in our sample (
erg s-1) when an average extinction E(B-V) = 0.4
0.2 is taken into account (see discussion in Netzer et al. 2006).
It is interesting to note, however, that the increase with mass
of the fraction of radio selected AGN (i.e. objects in which the
radio emission is not related to star formation processes) is much
steeper than what we find for our sample of AGN. This is true both
locally (Best et al. 2005) and up to
(Smolcic et al. 2009). This difference in the mass-dependence of the fraction of radio versus X-ray or optically selected AGN is likely
to be due to different accretion modes in these different classes of AGN (Smolcic et al. 2009).
4.4 Eddington ratio distribution function
From our multiwavelength database we were able to estimate another
important parameter for the characterization of AGN emission,
namely its Eddington ratio (
), which gives the ratio at which the AGN is emitting compared to its expected Eddington luminosity (
).
Indeed, the ratio between the X-ray luminosity and the host galaxy
stellar mass is proportional to the AGN Eddington ratio, the
proportionality factor depending on the bolometric correction between
the X-ray and the bolometric luminosity, and on the
/M* ratio. We derived the Eddington ratio distribution function for our sources as explained below.
First, we calculated the ratio between the X-ray luminosity and the stellar mass for each obscured AGN with
in our sample (
). Figure 11 shows the distribution of this quantity for the sources in two different redshifts bins at z>1 (solid histograms in
the upper left and bottom panels of Fig. 11). The median (seminterquartiles) values of the distributions of
are similar,
= 32.21 (0.61) and
= 32.73 (0.42) for the two redshift bins, respectively. When only the sources with
erg s-1 are considered, the median values of the
distribution are almost identical,
= 32.87 and
= 32.83, respectively (dashed histograms).
In the z=1-2 redshift bin, where the uncertainties on
the SSFR estimates are reasonably small and therefore we are able to
divide the sample in ``active'' and ``inactive'' galaxies (see
Sect. 4.2), we also plotted the distribution of
for the most luminous sources (
erg s-1) separately for these two populations (upper right panel of Fig. 11). We found that AGN in ``inactive'' host galaxies tend to have lower
values (dashed histogram in the upper right panel of Fig. 11)
than those in ``active'' host galaxies although the statistical
significance of this result is only marginal, the probability that the
two distributions are drawn from the same parent population
being 5%, according to a Kolmogorov-Smirnov statistics.
Then, we assumed the local relation between
and M* (e.g. Marconi & Hunt 2003,
M*) to convert the stellar mass in a BH mass
and calculate
,
and we applied a bolometric correction of 20 (e.g. Marconi et al. 2004) to derive the bolometric luminosity from the X-ray luminosity. With these assumptions, the median values of
correspond to
0.02-0.08 (
0.1 for the
subsample).
We note that, given the systematics uncertainties affecting the
derivation of the stellar masses (see discussion in 2.4), and the
fact that total stellar masses are derived and not bulge
masses, the derived BH massses are most likely upper limits, this
translating in conservative (lower-limits) estimate of the
.
On the other hand, recently, it has been suggested that the average
-host galaxy mass ratio at z=1-2 is higher by a factor
than the local value (McLure et al. 2006; Peng et al. 2006; Merloni et al. 2009). Should this be the case, our estimates of
should be reduced by a factor of
at z=1-2 and possibly more in the higher redshift bin if the average
ratio further increases to higher redshift.
![]() |
Figure 11: The distribution of the ratio between the X-ray luminosity and the stellar mass for the obscured AGN in our sample (empty histograms), in two different redshifts bins as labeled. |
Open with DEXTER |
5 Discussions and conclusions
We presented a new catalog of the counterparts of the 179 extragalactic X-ray sources detected in the 1 Ms Chandra observation of the MUSIC/CDFS/GOODS field and an extensive analysis of the host galaxies properties of obscured AGN.
We quantified the bias in the determination of the counterparts of
X-ray selected sources when the match is limited to optical catalogs
only, with respect to the combined use of optical and near infrared
(deep K-band
and IRAC) data. We estimate that the fraction of misidentified X-ray
sources previously reported in the literature is of the order of ,
and rises up to
% when optically faint (z>24) sources are considered (see Fig. 1); the use of an optically-based catalog biases the identification against the most extreme, obscured sources
therefore preventing the exact knowledge of the multiwavelength properties of the X-ray counterparts.
In order to study the host-galaxies of obscured AGN, we defined a sample of 116 ``bona fide'' obscured AGN, by selecting sources without broad lines in the optical spectra and with small optical nuclear emission with respect to the host galaxy optical emission (see Fig. 3). From eyeball inspection of the host galaxy morphology, we found a variety of cases, and a disturbed morphology (due to activity/merging/star formation) in more than half of the sub-sample for which a morphological classification could be made (see Fig. 4).
We investigated the optical to infrared colors of these obscured AGN.
The most striking result is that half of the X-ray selected obscured
AGN in the redshift range
z=1.5-4.0 have a 8.0 m to 4.5
m flux ratio <2 and according to Pope et al. (2008) would have been classified as ``star-formation'' dominated objects (see Fig. 6); 50% of them in addition have 24.0
m to 8.0
m flux ratio <5, where the original Pope et al. (2008) diagram is almost empty. Moreover, previous analysis based on Chandra stacking analysis of sources
with high MIR/O (Daddi et al. 2007a; Fiore et al. 2008, 2009) claimed a large contribution (
%) of heavily obscured (Compton Thick) sources among the stacked population (but see also discussion in Donley et al. 2008).
This suggests that 1) the accretion activity in high-redshift
sources is unambiguously revealed thanks to the presence of a strong
X-ray emission; and 2) the star-formation region as defined by
Pope et al. (2008)
contains not only objects in which the bolometric luminosity is
dominated by the star-formation processes, but also a not-negligible
number of objects hosting candidate obscured/Compton Thick AGN. In
conclusion, our obscured AGN show Spitzer colors consistent with both
an AGN dominated continuum and a starburst dominated continuum in
the MIR.
We found that the hosts of obscured AGN are redder in (U-V) rest frame than the overall galaxy population at the same redshift: in particular, obscured AGN mainly populate the red sequence and the green valley in the color-magnitude plots (Fig. 7, left panel), in agreement with the results of Nandra et al. (2007), and Silverman et al. (2007). For the MUSIC sample the U-V galaxy colors are strongly correlated with the K band absolute magnitude (Fig. 8, left panel), and therefore with the galaxy stellar mass, with the most massive systems having a redder color. The hosts of the obscured AGN are therefore found in the red-massive tail of the distribution of optically selected galaxies in all three redshift bins considered. AGN feedback is often invoked as one of the main responsible for the observed galaxy colors (Nandra et al. 2007; Hasinger 2008). However, it is well known that the main ingredient for nuclear activity is the presence of a SMBH in galaxy nuclei, and that SMBHs are found nearly exclusively in massive galaxies (e.g. Magorrian et al. 1998). Therefore, it is not truly surprising to find AGN hosted in massive galaxies and the simple presence of AGN in massive red galaxies is not enough to argue for a significant feedback effect on the observed colors, because of the strong color-mass correlation. Were AGN feedback responsible for the observed red colors, since galaxy colors are strongly correlated with the galaxy mass and AGN are found preferably in massive galaxies, then AGN feedback should be also considered as one of the main players in the building of the galaxy mass-color correlation.
We found that about 2/3 of the obscured AGN hosts at all redshifts show substantial (
yr-1) star formation activity (Fig. 9 left panel) and about half live in galaxies which are still actively forming stars with respect to their mass (Fig. 9,
right panel). For these sources, the observed red colors are likely due
to dust extinction rather than evolved stellar population. We then
conclude that a significant fraction of obscured AGN live in massive,
dusty star-forming galaxies with red optical colors. This result is in
qualitative agreement with the morphological analysis. Higher
luminosity X-ray selected AGN are not systematically found in objects
with the highest SSFR (see Fig. 9 right panel), in agreement with Alonso-Herrero et al. (2008).
We compared the number of obscured AGN and of all X-ray selected AGN to
the number of field galaxies in broad bins of galaxy stellar mass (
)
and redshifts (z=0.6-1, z=1-2, z=2-4). We find that the AGN fraction increases with the host galaxy stellar mass, from
1% at
to
30% at
(see also Yamada et al. 2009),
and the actual trend of increasing AGN fraction as a function of
the stellar mass is probably steeper given the uncompleteness of the
MUSIC sample at
(see Sect. 4.4). The uncertainties on the stellar mass
estimate from the SED fitting have the effect of shifting at lower
masses (of
dex) the datapoints, leaving the total fraction and the slope unchanged.
We compared this trend with that observed in the local Universe (Best et al. 2005) for AGN with luminosity above similar thresholds. While the observed trend is the same, in all the investigated redshift bins the AGN fraction is higher than that observed in the local Universe, and it could likely be even higher. In fact, we are comparing here AGN selected with two different methods: forbidden emission line luminosity (SDSS) and X-ray emission (GOODS). The latter sample does not contain most Compton thick AGN. On the other hand, Compton thick AGN may well be present in [OIII] selected AGN samples. The fraction of Compton thick AGN not directly detected in deep Chandra surveys is estimated between 40% and 100% of the X-ray selected AGN, using infrared selection or other techniques (see Donley et al. 2008; Fiore et al. 2009, and references therein). Therefore, under the simplest assumption that this Compton Thick AGN fraction is constant with the galaxy mass, the discrepancy observed in Fig. 10 can increase by up of a factor of 2.
The fraction of active galaxies to the total galaxy population is proportional to the AGN duty cycle. Our results would thus suggest higher AGN duty cycles at z=2-4 than at z=0, in agreement with expectations from most recent semi-analytic models (e.g. Menci et al. 2008), in which at higher redshift the AGN activity is present in a larger number of galaxies than locally.
The fact that the most luminous obscured AGN are found in the most
massive galaxies at all investigated redshifts may suggest that the
of the obscured AGN is similar, particularly in the case of the most luminous sources (
> 43 erg s-1),
for which the threshold in luminosity introduces a bias against the
sources accreting at lower rates in the lowest redshift bin. Assuming
the local Magorrian relation between
and M* (e.g. Marconi & Hunt 2003) and a bolometric correction of 20 (e.g. Marconi et al. 2004) the median values of
= 32.83 and
= 32.87 correspond to
.
Although suffering from large uncertainties associated with the stellar
mass and BH mass estimates (see Sect. 4.4), this value
can be considered as representative of the accretion state of the most
luminous, obscured AGN in the present sample.
Similar results are obtained from a comprehensive analysis of host
galaxies properties of Chandra Deep Field North X-ray sources sources
at z=2-4 (Yamada et al. 2009)
and are also typical of unobscured type 1 AGN at z>1 (Merloni et al. 2009; Trump et al. 2009).
It is instructive to compare these findings with similar estimates in the local Universe. Kauffmann & Heckman (2009) using SDSS data found an average
value of
and a log normal distribution for this parameter, for AGN hosted in
galaxies with significant on-going star formation, while AGN in
inactive host galaxies follow a power-law distribution.
As discussed in Sect. 4.4, and shown in the upper right panel
of Fig. 11,
we also find marginal evidence of different accretion rates
distributions for the populations of AGN in ``active'' and ``inactive''
host galaxies.
6 Summary
Taking advantage of the 1Ms CDFS Chandra observations and the associated deep multiwavelength follow-up, we have studied the host galaxies properties of obscured AGN at z>1.
The most important results of our analysis can be summarized as follows:
- X-ray selected (
erg s-1) AGN show Spitzer colors consistent with both AGN and starburst dominated infrared continuum and only the combination of the two selections (X-ray + Mid-infrared) can fully characterize the growth of SMBH at z > 1;
- the host galaxies of X-ray selected obscured AGN are all massive and in 50% of the cases are also actively forming stars. Their red colors are due to dust extinction;
- the X-ray selected AGN fraction increases with the stellar mass up to a value of
% at z>1 and M*>3
, a fraction significantly higher than what is observed in the local Universe (
%) for AGN of similar luminosities;
- the median
value is between 2% and 10% depending on the assumed
ratio, the X-ray luminosity threshold and the host galaxies properties of the obscured AGNs.
The authors are grateful to P. N. Best for deriving the SDSS OIII source fractions as a function of stellar mass with an optical luminosity cut that matches the present X-ray data. M.B. and F.F. acknowledge useful discussions with Roberto Gilli, Paolo Tozzi, John Silverman, Cristian Vignali, G. Cesare Perola, Silvia Bonoli, Niel Brandt, Bin Luo, Hagai Netzer and Jochen Greiner. We thank the anonymous referee for insightful and detailed comments which helped in the presentation of the results. M.B. acknowledges support from the XMM-Newton DLP grant 50-)G-0502; F.F. and A.C. acknowledge support from ASI contracts I/088/06/0 and I/016/07/0 and PRIN/MIUR 2006-02-5203.
References
- Alexander, D. M., Bauer, F. E., Brandt, W. N., et al. 2003, AJ, 126, 539 [NASA ADS] [CrossRef]
- Alexander, D. M., Bauer, F. E., Chapman, S. C., et al. 2005, ApJ, 632, 736 [NASA ADS] [CrossRef]
- Alonso-Herrero, A., Perez-Gonzalez, P. G., Rieke, G. H., et al. 2008, ApJ, 677, 127 [NASA ADS] [CrossRef]
- Antonucci, R. R. 1993, ARA&A, 31, 473 [NASA ADS]
- Antonucci, R. R., & Miller, J. S. 1985, ApJ, 297, 621 [NASA ADS] [CrossRef]
- Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [NASA ADS] [CrossRef]
- Bauer, F. E., Alexander, D. M., Brandt, W. N., et al. 2004, AJ, 128, 2048 [NASA ADS] [CrossRef]
- Best, P. N., Kauffmann, G., Heckman, T. M., et al. 2005, MNRAS, 362, 35 [NASA ADS]
- Brusa, M., Comastri, A., Daddi, E., et al. 2005, A&A, 432, 69 [NASA ADS] [CrossRef] [EDP Sciences]
- Brusa, M., Zamorani, G., Comastri, A., et al. 2007, ApJS, 172, 353 [NASA ADS] [CrossRef]
- Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef]
- Bundy, K., Georgakakis, A., Nandra, K., et al. 2008, ApJ, 681, 931 [NASA ADS] [CrossRef]
- Cardamone, C. N., Urry, C. M., Damen, M., et al. 2008, ApJ, 680, 130 [NASA ADS] [CrossRef]
- Cristiani, S., Appenzeller, I., Arnouts, S., et al. 2000, A&A, 359, 489 [NASA ADS]
- Comastri, A., & Brusa, M. 2008, AN, 329, 122 [NASA ADS]
- Daddi, E., Renzini, A., Pirzkal, N., et al. 2005, ApJ, 626, 680 [NASA ADS] [CrossRef]
- Daddi, E., Alexander, D. M., Dickinson, M., et al. 2007a, ApJ, 670, 173 [NASA ADS] [CrossRef]
- Daddi, E., Dickinson, M., Morrison, G., et al. 2007b, ApJ, 670, 156 [NASA ADS] [CrossRef]
- De Santis, C., Grazian, A., Fontana, A., & Santini, P. 2007, New Astron., 12, 271 [NASA ADS] [CrossRef]
- Dey, A., Soifer, B. T., Desai, V., et al. 2008, ApJ, 677, 943 [NASA ADS] [CrossRef]
- Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 [NASA ADS] [CrossRef]
- Donley, J. L., Rieke, G. H., Perez-Gonzalez, P. G., & Barro, G. 2008, ApJ, 687, 111 [NASA ADS] [CrossRef]
- Fabian, A. C. 1999, MNRAS, 308, L39 [NASA ADS] [CrossRef]
- Ferreras, I., Lisker, T., Carollo, C. M., Lilly, S. J., & Mobasher, B., 2005, ApJ, 635, 243 [NASA ADS] [CrossRef]
- Fiore, F., Grazian, A., Santini, P., et al. 2008, ApJ, 692, 74
- Fiore, F., Puccetti, S., Brusa, M., et al. 2009, ApJ, 693, 447 [NASA ADS] [CrossRef]
- Fontana, A., Salimbeni, S., Grazian, A., et al. 2006, 459, 745
- Fontana, A., Santini, P., Grazian, A., et al. 2009, A&A, 501, 15 [NASA ADS] [CrossRef] [EDP Sciences]
- 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]
- Giacconi, R., Zirm, A., Wang, J., et al. 2002, ApJS, 139, 369 [NASA ADS] [CrossRef]
- Giavalisco, M., Ferguson, H. C., Koekemoer, A. M., et al. 2004, ApJ, 600, L93 [NASA ADS] [CrossRef]
- Gilli, R., Daddi, E., Zamorani, G., et al. 2005, A&A, 430, 811 [NASA ADS] [CrossRef] [EDP Sciences]
- Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 [NASA ADS] [CrossRef] [EDP Sciences]
- Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 2004, ApJ, 600, 580 [NASA ADS] [CrossRef]
- Grazian, A., Fontana, A., de Santis, C., et al. 2006, A&A, 449, 951 [NASA ADS] [CrossRef] [EDP Sciences]
- Hasinger, G. 2008, A&A, 490, 905 [NASA ADS] [CrossRef] [EDP Sciences]
- Hickox, R., Jones, C., Forman, W. R., et al. 2009, ApJ, 696, 891 [NASA ADS] [CrossRef]
- Hopkins, P. F., Hernquist, L., Cox, T. J., & Keres, D. 2008, ApJS, 175, 365 [NASA ADS]
- Kauffmann, G., & Heckman, T. M. 2009, MNRAS, 397, 153 [CrossRef]
- Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055 [NASA ADS] [CrossRef]
- Koekemoer, A. M., Alexander, D. M., Bauer, F. B., et al. 2004, ApJ, 600, L123 [NASA ADS] [CrossRef]
- Ilbert, O., capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236 [NASA ADS] [CrossRef]
- Ivison, R. J., Greve, T. R., Serjeant, S., et al. 2004, ApJS, 154, 124 [NASA ADS] [CrossRef]
- Le Fevre, O., Vettolani, G., Paltani, S., et al. 2004, A&A, 428, 1043 [NASA ADS] [CrossRef] [EDP Sciences]
- Lehmer, B., Brandt, W. N., Alexander, D. M., et al. 2005, ApJS, 161, 21 [NASA ADS] [CrossRef]
- Magorrian, J., Tremaine, S., Richstone, D., et al. 1998, AJ, 115, 2285 [NASA ADS] [CrossRef]
- Mainieri, V., Rosati, P., Tozzi, P., et al. 2005, A&A, 437, 805 [NASA ADS] [CrossRef] [EDP Sciences]
- Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 [NASA ADS] [CrossRef]
- Marconi, A., Risaliti, G., Gilli, R., et al. 2004, MNRAS, 351, 169 [NASA ADS] [CrossRef]
- Martinez-Sansigre, A., Rawlings, S., Lacy, M., et al. 2005, Nature, 436, 666 [NASA ADS] [CrossRef]
- Martinez-Sansigre, A., Rawlings, S., Lacy, M., et al. 2006, MNRAS, 370, 1479 [NASA ADS] [CrossRef]
- McLure, R. J., Jarvis, M. J., Targett, T. A., Dunlop, J. S., & Best, P. N. 2006, MNRAS, 368, 1395 [NASA ADS] [CrossRef]
- Menci, N., Fontana, A., Giallongo, E., Grazian, A., & Salimbeni, S. 2006, ApJ, 647, 753 [NASA ADS] [CrossRef]
- Menci, N., Fiore, F., Puccetti, S., & Cavaliere, A. 2008, ApJ, 686, 219 [NASA ADS] [CrossRef]
- Merloni, A., & Heinz, S. 2008, MNRAS, 338, 1011 [NASA ADS]
- Merloni, A., Bongiorno, A., Bolzonella, M., et al. 2009, ApJ, in press,[arXiv:0910.4970]
- Mignoli, M., Cimatti, A., Zamorani, G., et al. 2005, A&A, 437, 883 [NASA ADS] [CrossRef] [EDP Sciences]
- Netzer, H., Mainieri, V., Rosati, P., & Trakhtenbrot, B. 2006, A&A, 453, 525 [NASA ADS] [CrossRef] [EDP Sciences]
- Nandra, K., Georgakakis, A., Willmer, C. N. A., et al. 2007, ApJ, 660, L11 [NASA ADS] [CrossRef]
- Norman, C., Hasinger, G., Giacconi, R., et al. 2002, ApJ, 571, 218 [NASA ADS] [CrossRef]
- Page, M. J., Stevens, J. A., Ivison, R. J., & Carrera, F. J. 2004, ApJ, 611, L85 [NASA ADS] [CrossRef]
- Papovich, C., Rudnick, G., Le Floch, E., et al. 2007, ApJ, 668, 45 [NASA ADS] [CrossRef]
- Peng, C. Y., Impey, C. D., Rix, H.-W., et al. 2006, ApJ, 649, 616 [NASA ADS] [CrossRef]
- Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81 [NASA ADS] [CrossRef]
- Pope, A., Bussman, R. S., Dey, A., et al. 2008, ApJ, 689, 127 [NASA ADS] [CrossRef]
- Popesso, P., Dickinson, M., Nonino, M., et al. 2009, A&A, 494, 443 [NASA ADS] [CrossRef] [EDP Sciences]
- Pozzi, F., Vignali, C., Comastri, A., et al. 2007, A&A, 468, 603 [NASA ADS] [CrossRef] [EDP Sciences]
- Ranalli, P., Comastri, A., & Setti, G. 2003, A&A, 399, 39 [NASA ADS] [CrossRef] [EDP Sciences]
- Ranalli, P., Comastri, A., & Setti, G. 2005, A&A, 440, 23 [NASA ADS] [CrossRef] [EDP Sciences]
- Rigby, J. R., Rieke, G. H., Perez-Gonzalez, P. G., et al. 2005, ApJ, 627, 134 [NASA ADS] [CrossRef]
- Rodighiero, G., Cimatti, A., Franceschini, A., et al. 2007, A&A, 470, 21 [NASA ADS] [CrossRef] [EDP Sciences]
- Rovilos, E., & Georgantopoulos, I. 2007, A&A, 475, 115 [NASA ADS] [CrossRef] [EDP Sciences]
- Salpeter, E. E. 1955, ApJ, 121, 161 [NASA ADS] [CrossRef]
- Santini, P., Fontana, A., Grazian, A., et al. 2009, A&A, 504, 751 [CrossRef] [EDP Sciences]
- Silk, J., & Rees, M. J. 1998, A&A, 311, L1 [NASA ADS]
- Silverman, J. D., Green, P. J., Barkhouse, W. A., et al. 2008, ApJ, 679, 118 [NASA ADS] [CrossRef]
- Silverman, J. D., et al. 2009, ApJ, in press
- Smolcic, V., Zamorani, G., Schinnerer, E., et al. 2009, ApJ, 696, 24 [NASA ADS] [CrossRef]
- Stern, D., Eisenhardt, P., Gorjian, V., et al. 2005, ApJ, 631, 163 [NASA ADS] [CrossRef]
- Stevens, J. A., et al. 2006, MNRAS, 360, 610 [NASA ADS] [CrossRef]
- Szokoly, G. P., Bergeron, J., Hasinger, G., et al. 2004, ApJS, 155, 271 [NASA ADS] [CrossRef]
- Tozzi, P., Rosati, P., Nonino, M., et al. 2001, ApJ, 562, 42 [NASA ADS] [CrossRef]
- Tozzi, P., Gilli, R., Mainieri, V., et al. 2006, A&A, 451, 457 [NASA ADS] [CrossRef] [EDP Sciences]
- Tozzi, P., Mainieri, V., Rosati, P., et al. 2009, ApJ, 698, 740 [NASA ADS] [CrossRef]
- Trump, J. R., Impey, C. D., Kelly, B. C., et al. 2009, ApJ, 700, 49 [NASA ADS] [CrossRef]
- Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 [NASA ADS] [CrossRef]
- Vanzella, E., Cristiani, S., Dickinson, M., et al. 2005, A&A, 434, 53 [NASA ADS] [CrossRef] [EDP Sciences]
- Vanzella, E., Cristiani, S., Dickinson, M., et al. 2006, A&A, 454, 23 [NASA ADS] [CrossRef]
- Vanzella, E., Cristiani, S., Dickinson, M., et al. 2008, A&A, 478, 83 [NASA ADS] [CrossRef] [EDP Sciences]
- Wolf, C., Meisenheimer, K., Kleinheinrich, M., et al. 2004, A&A, 421, 913 [NASA ADS] [CrossRef] [EDP Sciences]
- Yamada, T., Kajisawa, M., Akyiama, M., et al. 2009, ApJ, 699, 1354 [NASA ADS] [CrossRef]
- Yan, L., Sajina, A., Fadda, D., et al. 2007, ApJ, 658, 778 [NASA ADS] [CrossRef]
- Zheng, W., Mikles, V. J., Mainieri, V., et al. 2004, ApJS, 155, 73 [NASA ADS] [CrossRef]
Online Material
Table 1: Counterparts to CDFS X-ray sources: positions, fluxes, redshifts and luminosities.
Table 2: Multiwavelength photometry for counterparts to CDFS X-ray sources.
Footnotes
- ... AGN
- Tables 1 and 2 are only available in electronic form at http://www.aanda.org
- ... by A03
- We do not discuss here the 10 sources in the MUSIC area present only in Giacconi et al. (2002).
- ...2001)
- We haven't used the fluxes tabulated by A03 because they have been calculated with different (source-by-source) energy indeces and make it difficult to correct for the absorption and derive rest-frame, unobscured luminosity (see 2.2).
- ... sample
- We recall that our sample contains 16 sources detected only in the A03 sample and not present at all in the Tozzi et al. (2006) spectral analysis.
- ... redshifts
- We note here that photo-z accuracy, in general, is not the leading source of error in deriving the galaxies and AGN parameters (see also discussion in Ilbert et al. 2009).
- ... only
- This threshold further assures that the observed X-ray
emission is mostly due to AGN activity, with a negligible
contribution from star-formation, even for the objects with the highest
observed SFR (
yr-1, Ranalli et al. 2003).
All Tables
Table 1: Counterparts to CDFS X-ray sources: positions, fluxes, redshifts and luminosities.
Table 2: Multiwavelength photometry for counterparts to CDFS X-ray sources.
All Figures
![]() |
Figure 1:
Chandra, ACS (z-band), K-band and IRAC-4.5 micron cutouts (5'' |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Redshift distribution of the X-ray sources. Open histogram: all sources; right diagonal (blue) shaded histogram: spectroscopic sample; left diagonal (red) shaded histogram: photometric sample. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Distribution of the z-band count ratio for the ``obscured'' AGN sample (116 objects, red long dashed histogram) and the BL/pointlike sample (26 objects, blue dotted histogram). |
Open with DEXTER | |
In the text |
![]() |
Figure 4: ACS (z-band) images of obscured AGN counterparts with disturbed morphology or interacting systems (AID 88, 107, 226 and 84). |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Left panel: 24 |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Spitzer color color diagram used to isolate star-formation and AGN
dominated high-z galaxies (Pope et al. 2008). The small green points represent the distribution of field galaxies in the 24 |
Open with DEXTER | |
In the text |
![]() |
Figure 7: Left panel: U-V (rest frame) vs. the absolute V-band magnitude for X-ray obscured AGN (squares) and the underlying optically selected galaxy population (contours). The sizes and colors of the symbols reflect the AGN X-ray luminosity (as in Figs. 5 and 6). Right panel: [3.6]-[4.5] color vs. the U-V (rest frame) for obscured X-ray sources (squares) and the underlying optically selected galaxy population (contours). The upper panels show the comparison of the X-ray selected sources with the optically-selected (z-band) galaxy population in three different redshift bins: z = 0.6-1.0, z = 1-2, and z = 2.4 from left to right, respectively. The lower panels show the comparison for the subsample of objects detected also at 24 micron, in the same redshift bins. The sizes and colors of the symbols reflect the AGN X-ray luminosity as in Figs. 5 and 6. Unobscured AGN are not shown here. |
Open with DEXTER | |
In the text |
![]() |
Figure 8: Left panel: U-V (rest frame) vs. the absolute K-band magnitude for the obscured AGN (squares) and the underlying galaxy population (contours). The upper panels show the comparison of the X-ray selected sources with the optically-selected (z-band) galaxy population in three different redshift bins: z=0.6-1.0, z=1-2, and z=2.4 from left to right, respectively. The lower panels show the comparison for the subsample of objects detected also at 24 micron, in the same redshift bins. The shapes and colors of the symbols reflect the AGN X-ray luminosity (see Fig. 5). Right panel: stellar mass M* vs. redshift for field galaxies (crosses) and our ``obscured'' AGN (squares and circles, same notation as in Fig. 5). The vertical lines mark the redshift range we chose for our analysis. |
Open with DEXTER | |
In the text |
![]() |
Figure 9:
(Left panel): the distribution of the SFR for both the optically ( upper) and MIPS ( lower) selected sample of field galaxies (black open histogram) and obscured AGN (red filled histogram) in the redshift interval z=1-4 and in a specific stellar mass range (
|
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Fraction of obscured AGN with
|
Open with DEXTER | |
In the text |
![]() |
Figure 11: The distribution of the ratio between the X-ray luminosity and the stellar mass for the obscured AGN in our sample (empty histograms), in two different redshifts bins as labeled. |
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.