The X-ray luminosity function of Active Galactic Nuclei in the redshift interval z=3-5

We combine deep X-ray survey data from the Chandra observatory and the wide-area/shallow XMM-XXL field to estimate the AGN X-ray luminosity function in the redshift range z=3-5. The sample consists of nearly 340 sources with either photometric (212) or spectroscopic (128) redshift in the above range. The combination of deep and shallow survey fields provides a luminosity baseline of three orders of magnitude, Lx(2-10keV)~1e43-1e46erg/s at z>3. We follow a Bayesian approach to determine the binned AGN space density and explore their evolution in a model-independent way. Our methodology accounts for Poisson errors in the determination of X-ray fluxes and uncertainties in photometric redshift estimates. We demonstrate that the latter is essential for unbiased measurement of space densities. We find that the AGN X-ray luminosity function evolves strongly between the redshift intervals z=3-4 and z=4-5. There is also suggestive evidence that the amplitude of this evolution is luminosity dependent. The space density of AGN with Lx<1e45erg/s drops by a factor of 5 between the redshift intervals above, while the evolution of brighter AGN appears to be milder. Comparison of our X-ray luminosity function with that of UV/optical selected QSOs at similar redshifts shows broad agreement at bright luminosities, Lx>1e45erg/s. The faint-end slope of UV/optical luminosity functions however, is steeper than for X-ray selected AGN. This implies that the type-I AGN fraction increases with decreasing luminosity at z>3, opposite to trends established at lower redshift. We also assess the significance of AGN in keeping the hydrogen ionised at high redshift. Our X-ray luminosity function yields ionising photon rate densities that are insufficient to keep the Universe ionised at redshift z>4. A source of uncertainty in this calculation is the escape fraction of UV photons for X-ray selected AGN.

Large scale clustering measurements with photometric redshifts: comparing the dark matter halos of X-ray AGN, star-forming and passive galaxies at z~1

We combine multiwavelength data in the AEGIS-XD and C-COSMOS surveys to measure the typical dark matter halo mass of X-ray selected AGN [Lx(2-10keV)>1e42 erg/s] in comparison with far-infrared selected star-forming galaxies detected in the Herschel/PEP survey (PACS Evolutionary Probe; Lir>1e11 solar) and quiescent systems at z~1. We develop a novel method to measure the clustering of extragalactic populations that uses photometric redshift Probability Distribution Functions in addition to any spectroscopy. This is advantageous in that all sources in the sample are used in the clustering analysis, not just the subset with secure spectroscopy. The method works best for large samples. The loss of accuracy because of the lack of spectroscopy is balanced by increasing the number of sources used to measure the clustering. We find that X-ray AGN, far-infrared selected star-forming galaxies and passive systems in the redshift interval 0.6<z<1.4 are found in halos of similar mass, $log M_{DMH}/(M_{odot},h^{-1})approx13.0$. We argue that this is because the galaxies in all three samples (AGN, star-forming, passive) have similar stellar mass distributions, approximated by the J-band luminosity. Therefore all galaxies that can potentially host X-ray AGN, because they have stellar masses in the appropriate range, live in dark matter haloes of $log M_{DMH}/(M_{odot},h^{-1})approx13.0$ independent of their star-formation rates. This suggests that the stellar mass of X-ray AGN hosts is driving the observed clustering properties of this population. We also speculate that trends between AGN properties (e.g. luminosity, level of obscuration) and large scale environment may be related to differences in the stellar mass of the host galaxies.

Investigating evidence for different black hole accretion modes since redshift z~1

Chandra data in the COSMOS, AEGIS-XD and 4Ms CDFS are combined with optical/near-IR photometry to determine the rest-frame U-V vs V-J colours of X-ray AGN hosts at mean redshifts 0.40 and 0.85. This combination of colours (UVJ) provides an efficient means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasises differences between AGN split by their UVJ colours. AGN in quiescent galaxies are dominated by spheroids, while star-forming hosts are split between bulges and disks. The UVJ diagram of AGN hosts is then used to set limits on the accretion density associated with evolved and star-forming systems. Most of the black hole growth since z~1 is associated with star-forming hosts. Nevertheless, ~15-20% of the X-ray luminosity density since z~1, is taking place in the quiescent region of the UVJ diagram. For the z~0.40 subsample, there is tentative evidence (2sigma significance), that AGN split by their UVJ colours differ in Eddington ratio. AGN in star-forming hosts dominate at high Eddington ratios, while AGN in quiescent hosts become increasingly important as a fraction of the total population toward low Eddington ratios. At higher redshift, z~0.8, such differences are significant at the 2sigma level only at Eddington ratios >1e-3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of SMBHs at the centres of galaxies. We compare our results with the GALFORM semi-analytic model, which postulates two black hole fuelling modes, the first linked to star-formation and the second occuring in passive galaxies. GALFORM predicts a larger fraction of black hole growth in quiescent galaxies at z<1, compared to the data. Relaxing the strong assumption of the model that passive AGN hosts have zero star-formation rate could reconcile this disagreement.

Constraining the fraction of Compton-thick AGN in the Universe by modelling the diffuse X-ray background spectrum

This paper investigates what constraints can be placed on the fraction of Compton-thick (CT) AGN in the Universe from the modeling of the spectrum of the diffuse X-ray background (XRB). We present a model for the synthesis of the XRB that uses as input a library of AGN X-ray spectra generated by the Monte Carlo simulations described by Brightman & Nandra. This is essential to account for the Compton scattering of X-ray photons in a dense medium and the impact of that process on the spectra of obscured AGN. We identify a small number of input parameters to the XRB synthesis code which encapsulate the minimum level of uncertainty in reconstructing the XRB spectrum. These are the power-law index and high energy cutoff of the intrinsic X-ray spectra of AGN, the level of the reflection component in AGN spectra and the fraction of CT AGN in the Universe. We then map the volume of the space allowed to these parameters by current observations of the XRB spectrum in the range 3-100 keV. One of the least constrained parameters is the fraction of CT AGN. Statistically acceptable fits to the XRB spectrum at the 68% confidence level can be obtained for CT fractions in the range 5-50%. This is because of degeneracies among input parameters to the XRB synthesis code and uncertainties in the modeling of AGN spectra (e.g. reflection). The most promising route for constraining the fraction of CT AGN in the Universe is via the direct detection of those sources in high energy (>10keV) surveys. It is shown that the observed fraction of CT sources identified in the SWIFT/BAT survey, limits the intrinsic fraction of CT AGN, at least at low redshift, to 10-20% (68% confidence level). We also make predictions on the number density of CT sources that current and future X-ray missions are expected to discover. Testing those predictions will constrain the intrinsic fraction of CT AGN as a function of redshift.

The radio spectra of reddened 2MASS QSOs: evidence for young radio jets

Multifrequency radio continuum observations (1.4-22 GHz) of a sample of reddened QSOs are presented. We find a high incidence (13/16) of radio spectral properties, such as low frequency turnovers, high frequency spectral breaks or steep power-law slopes, similar to those observed in powerful compact steep spectrum (CSS) and gigahertz-peaked spectrum (GPS) sources. The radio data are consistent with relatively young radio jets with synchotron ages <1e6-1e7yr. This calculation is limited by the lack of high resolution (milli-arcsec) radio observations. For the one source in the sample that such data are available a much younger radio age is determined, <2e3yr, similar to those of GPS/CSS sources. These findings are consistent with claims that reddened QSOs are young systems captured at the first stages of the growth of their supermassive black holes. It also suggests that expanding radio lobes may be an important feedback mode at the early stages of the evolution of AGN.

Observational constraints on the physics behind the evolution of AGN since z ~ 1

We explore the evolution with redshift of the rest-frame colours and space densities of AGN hosts (relative to normal galaxies) to shed light on the dominant mechanism that triggers accretion onto supermassive black holes as a function of cosmic time. Data from serendipitous wide-area XMM surveys of the SDSS footprint (XMM/SDSS, Needles in the Haystack survey) are combined with Chandra deep observations in the AEGIS, GOODS-North and GOODS-South to compile uniformly selected samples of moderate luminosity X-ray AGN [L_X(2-10keV) = 1e41-1e44erg/s] at redshifts 0.1, 0.3 and 0.8. It is found that the fraction of AGN hosted by red versus blue galaxies does not change with redshift. Also, the X-ray luminosity density associated with either red or blue AGN hosts remains nearly constant since z=0.8. X-ray AGN represent a roughly fixed fraction of the space density of galaxies of given optical luminosity at all redshifts probed by our samples. In contrast the fraction of X-ray AGN among galaxies of a given stellar mass decreases with decreasing redshift. These findings suggest that the same process or combination of processes for fueling supermassive black holes are in operation in the last 5 Gyrs of cosmic time. The data are consistent with a picture in which the drop of the accretion power during that period (1dex since z=0.8) is related to the decline of the space density of available AGN hosts, as a result of the evolution of the specific star-formation rate of the overall galaxy population. Scenarios which attribute the evolution of moderate luminosity AGN since z approx 1 to changes in the suppermassive black hole accretion mode are not favored by our results.

Searching for an Intermediate Mass Black Hole in the Blue Compact Dwarf galaxy MRK 996

The possibility is explored that accretion on an intermediate mass black hole contributes to the ionisation of the interstellar medium of the Compact Blue Dwarf galaxy MRK996. Chandra observations set tight upper limits (99.7 per cent confidence level) in both the X-ray luminosity of the posited AGN, Lx(2-10keV)<3e40erg/s, and the black hole mass, <1e4/lambda Msolar, where lambda, is the Eddington ratio. The X-ray luminosity upper limit is insufficient to explain the high ionisation line [OIV]25.89mu m, which is observed in the mid-infrared spectrum of the MRK996 and is proposed as evidence for AGN activity. This indicates that shocks associated with supernovae explosions and winds of young stars must be responsible for this line. It is also found that the properties of the diffuse X-ray emission of MRK996 are consistent with this scenario, thereby providing direct evidence for shocks that heat the galaxys interstellar medium and contribute to its ionisation.

Infrared Excess Sources: Compton Thick QSOs, low luminosity Seyferts or starbursts?

We explore the nature of Infrared Excess sources (IRX), which are proposed as candidates for luminous L_X(2-10keV)>1e43erg/s Compton Thick (N_H>2e24cm^{-2}$) QSOs at z~2. Lower redshift, z~1, analogues of the distant IRX population are identified by firstly redshifting to z=2 the SEDs of all sources with secure spectroscopic redshifts in the AEGIS (6488) and the GOODS-North (1784) surveys and then selecting those that qualify as IRX sources at that redshift. A total of 19 galaxies are selected. The mean redshift of the sample is $zapprox1$. We do not find strong evidence for Compton Thick QSOs in the sample. For 9 sources with X-ray counterparts, the X-ray spectra are consistent with Compton Thin AGN. Only 3 of them show tentative evidence for Compton Thick obscuration. The SEDs of the X-ray undetected population are consistent with starburst activity. There is no evidence for a hot dust component at the mid-infrared associated with AGN heated dust. If the X-ray undetected sources host AGN, an upper limit of L_X(2-10keV) =1e43erg/s is estimated for their intrinsic luminosity. We propose that a large fraction of the $zapprox2$ IRX population are not Compton Thick QSOs but low luminosity [L_X(2-10keV)<1e43erg/s], possibly Compton Thin, AGN or dusty starbursts. It is shown that the decomposition of the AGN and starburst contribution to the mid-IR is essential for interpreting the nature of this population, as star-formation may dominate this wavelength regime.

Host galaxy morphologies of X-ray selected AGN: assessing the significance of different black hole fueling mechanisms to the accretion density of the Universe at z~1

We use morphological information of X-ray selected AGN hosts to set limits on the fraction of the accretion density of the Universe at z~1 that is not likely to be associated with major mergers. Deep X-ray observations are combined with high resolution optical data from the Hubble Space Telescope in the AEGIS, GOODS North and GOODS South fields to explore the morphological breakdown of X-ray sources in the redshift interval 0.5<z<1.3. The sample is split into disks, early-type bulge dominated galaxies, peculiar systems and point-sources in which the nuclear source outshines the host galaxy. The X-ray luminosity function and luminosity density of AGN at z~1 are then calculated as a function of morphological type. We find that disk-dominated hosts contribute 30pm9 per cent to the total AGN space density and 23pm6 per cent to the luminosity density at z~1. We argue that AGN in disk galaxies are most likely fueled not by major merger events but by minor interactions or internal instabilities. We find evidence that these mechanisms may be more efficient in producing luminous AGN (L_X>1e44 erg/s) compared to predictions for the stochastic fueling of massive black holes in disk galaxies.

X-ray selected AGN in groups at redshifts z~1

We explore the role of the group environment in the evolution of AGN at the redshift interval 0.7<z<1.4, by combining deep Chandra observations with extensive optical spectroscopy from the All-wavelength Extended Groth strip International Survey (AEGIS). The sample consists of 3902 optical sources and 71 X-ray AGN. Compared to the overall optically selected galaxy population, X-ray AGN are more frequently found in groups at the 99% confidence level. This is partly because AGN are hosted by red luminous galaxies, which are known to reside, on average, in dense environments. Relative to these sources, the excess of X-ray AGN in groups is significant at the 91% level only. Restricting the sample to 0.7<z<0.9 and M_B<-20mag in order to control systematics we find that X-ray AGN represent (4.7pm1.6) and (4.5pm1.0)% of the optical galaxy population in groups and in the field respectively. These numbers are consistent with the AGN fraction in low redshift clusters, groups and the field. The results above, although affected by small number statistics, suggest that X-ray AGN are spread over a range of environments, from groups to the field, once the properties of their hosts (e.g. colour, luminosity) are accounted for. There is also tentative evidence, significant at the 98% level, that the field produces more X-ray luminous AGN compared to groups, extending similar results at low redshift to z~1. This trend may be because of either cold gas availability or the nature of the interactions occurring in the denser group environment (i.e. prolonged tidal encounters).