Telescopes to be launched over the next decade-and-a-half, such as JWST, EUCLID, ATHENA and Lynx, promise to revolutionise the study of the high redshift Universe and greatly advance our understanding of the early stages of galaxy formation. We use a model that follows the evolution of the masses and spins of supermassive black holes (SMBHs) within a semi-analytic model of galaxy formation to make predictions for the Active Galactic Nucleus (AGN) luminosity function at $z\geq7$ in the broadband filters of JWST and EUCLID at near-infrared wavelengths, and ATHENA and Lynx at X-ray energies. The predictions of our model are relatively insensitive to the choice of seed black hole mass, except at the lowest luminosities ($L_{\mathrm{bol}}<10^{43}\mathrm{ergs^{-1}}$) and the highest redshifts ($z>10$). We predict that surveys with these different telescopes will select somewhat different samples of SMBHs, with EUCLID unveiling the most massive, highest accretion rate SMBHs, Lynx the least massive, lowest accretion rate SMBHs, and JWST and ATHENA covering objects inbetween. At $z=7$, we predict that typical detectable SMBHs will have masses, $M_{\mathrm{BH}}\sim10^{5-8}M_{\odot}$, and Eddington normalised mass accretion rates, $\dot{M}/\dot{M}_{\mathrm{Edd}}\sim0.6-2$. The SMBHs will be hosted by galaxies of stellar mass $M_{\star}\sim10^{8-10}M_{\odot}$, and dark matter haloes of mass $M_{\mathrm{halo}}\sim10^{11-12}M_{\odot}$. We predict that the detectable SMBHs at $z=10$ will have slightly smaller black holes, accreting at slightly higher Eddington normalised mass accretion rates, in slightly lower mass host galaxies compared to those at $z=7$, and reside in haloes of mass $M_{\mathrm{halo}}\sim10^{10-11}M_{\odot}$.

Understanding how Active Galactic Nuclei (AGN) evolve through cosmic time allows us to probe the physical processes that control their evolution. We use an updated model for the evolution of masses and spins of supermassive black holes (SMBHs), coupled to the latest version of the semi-analytical model of galaxy formation GALFORM using the Planck cosmology and a high resolution Millennium style dark matter simulation to make predictions for AGN and SMBH properties for $0 < z < 6$. We compare the model to the observed black hole mass function and the SMBH versus galaxy bulge mass relation at $z=0$, and compare the predicted bolometric, hard X-ray, soft X-ray and optical AGN luminosity functions to observations at $z < 6$, and find that the model is in good agreement with the observations. The model predicts that at $z<2$ and $L_{\mathrm{bol}} < 10^{43} \mathrm{ergs^{-1}}$, the AGN luminosity function is dominated by objects accreting in an Advection Dominated Accretion Flow (ADAF) disc state, while at higher redshifts and higher luminosities the dominant contribution is from objects accreting via a thin disc or at super-Eddington rates. The model also predicts that the AGN luminosity function at $z<3$ and $L_{\mathrm{bol}} < 10^{44} \mathrm{ergs^{-1}}$ is dominated by the contribution from AGN fuelled by quiescent hot halo accretion, while at higher luminosities and higher redshifts, the AGN luminosity function is dominated by the contribution from AGN fuelled by starbursts triggered by disc instabilities. We employ this model to predict the evolution of SMBH masses, Eddington ratios, and spins, finding that the median SMBH spin evolves very little for $0<z<6$.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

The northern tile of the wide-area and shallow XMM-XXL X-ray survey field is used to estimate the average dark matter halo mass of relatively luminous X-ray selected AGN [$\rm log\, L_X (\rm 2-10\,keV)= 43.6^{+0.4}_{-0.4}\,erg/s$] in the redshift interval $z=0.5-1.2$. Spectroscopic follow-up observations of X-ray sources in the XMM-XXL field by the Sloan telescope are combined with the VIPERS spectroscopic galaxy survey to determine the cross-correlation signal between X-ray selected AGN (total of 318) and galaxies (about 20,\,000). We model the large scales (2-25\,Mpc) of the correlation function to infer a mean dark matter halo mass of $\log M / (M_{\odot} \, h^{-1}) = 12.50 ^{+0.22} _{-0.30}$ for the X-ray selected AGN sample. This measurement is about 0.5\,dex lower compared to estimates in the literature of the mean dark matter halo masses of moderate luminosity X-ray AGN [$L_X (\rm 2-10\,keV)\approx 10^{42} - 10^{43}\,erg/s$] at similar redshifts. Our analysis also links the mean clustering properties of moderate luminosity AGN with those of powerful UV/optically selected QSOs, which are typically found in halos with masses few times $10^{12}\,M_{\odot}$. There is therefore evidence for a negative luminosity dependence of the AGN clustering. This is consistent with suggestions that AGN have a broad dark matter halo mass distribution with a high mass tail that becomes sub-dominant at high accretion luminosities. We further show that our results are in qualitative agreement with semi-analytic models of galaxy and AGN evolution, which attribute the wide range of dark matter halo masses among the AGN population to different triggering mechanisms and/or black hole fueling modes.

This is the fourth paper in a series that reports on our investigation of the clustering properties of active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the clustering of broad-line, luminous AGN at 0.16<z<0.36. We fit the H-alpha line profile in the SDSS spectra for all X-ray and optically-selected broad-line AGN, determine the mass of the super-massive black hole (SMBH), M_BH, and infer the accretion rate relative to Eddington (L/L_EDD). Since M_BH and L/L_EDD are correlated, we create AGN subsamples in one parameter while maintaining the same distribution in the other parameter. In both the X-ray and optically-selected AGN samples we detect a weak clustering dependence with M_BH and no statistically significant dependence on L/L_EDD. We find a difference of up to 2.7sigma when comparing the objects that belong to the 30% least and 30% most massive M_BH subsamples, in that luminous broad-line AGN with more massive black holes reside in more massive parent dark matter halos at these redshifts. These results provide evidence that higher accretion rates in AGN do not necessarily require dense galaxy environments in which more galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We also present semi-analytic models which predict a positive M_DMH dependence on M_BH, which is most prominent at M_BH ~ 10^{8-9} M_SUN.

We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.

PKS 2155-304 is one of the brightest extragalactic source in the X-ray and EUV bands, and is a prototype for the BL Lac class of objects. In this paper we investigate the large-scale environment of this source using new multi-object as well as long-slit spectroscopy, together with archival spectra and optical images. We find clear evidence of a modest overdensity of galaxies at z=0.11610, consistent with previous determinations of the BL Lac redshift. The galaxy group has a radial velocity dispersion of 250km/s and a virial radius of 0.22Mpc, yielding a role-of-thumb estimate of the virial mass of M(vir)~1.5x10$^{13}$Msun, i.e., one order of magnitude less than what observed in other similar objects. This result hints toward a relatively wide diversity in the environmental properties of BL Lac objects.

In a Universe where AGN feedback regulates star formation in massive galaxies, a strong correlation between these two quantities is expected. If the gas causing star formation is also responsible for feeding the central black hole, then a positive correlation is expected. If powerful AGNs are responsible for the star formation quenching, then a negative correlation is expected. Observations so far have mainly found a mild correlation or no correlation at all (i.e. a flat relation between star formation rate (SFR) and AGN luminosity), raising questions about the whole paradigm of "AGN feedback". In this paper, we report the predictions of the GALFORM semi-analytical model, which has a very strong coupling between AGN activity and quenching of star formation. The predicted SFR-AGN luminosity correlation appears negative in the low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive in the regime of the brightest AGN. Our predictions reproduce reasonably well recent observations by Rosario et al., yet there is some discrepancy in the normalisation of the correlation at low luminosities and high redshifts. Though this regime could be strongly influenced by observational biases, we argue that the disagreement could be ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the galaxies that dominate the regime where the observations imply a weak correlation are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these galaxies retain high enough molecular hydrogen contents to maintain relatively high star formation rates and strong infrared emission.

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.

We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of the coupled galaxy and black hole formation model GALFORM. We find good agreement when our calculation includes the hot-halo mode of accretion onto the central black hole. This type of accretion, which is additional to the common cold accretion during disk instabilities and galaxy mergers, is tightly coupled to the AGN feedback in the model. The hot-halo mode becomes prominent in DM haloes with masses greater than ~10^12.5 Msun, where AGN feedback typically operates, giving rise to a distinct class of moderate luminosity AGN that inhabit rich clusters and superclusters. Cold gas fuelling of the black hole cannot produce the observationally inferred DM halo masses of X-ray AGN. Switching off AGN feedback in the model results in a large population of luminous quasars (Lx(2-10 keV) > 10^44 erg/sec) in DM haloes with masses up to ~10^14 Msun, which is inconsistent with the observed clustering of quasars. The abundance of hot-halo AGN decreases significantly in the z~3-4 universe. At such high redshifts, the cold accretion mode is solely responsible for shaping the environment of moderate luminosity AGN. Our analysis supports two accretion modes (cold and hot) for the fuelling of supermassive black holes and strongly underlines the importance of AGN feedback in cosmological models both of galaxy formation and black hole growth.

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe.

We present the spatial clustering properties of 1466 X-ray selected AGN compiled from the Chandra CDF-N, CDF-S, eCDF-S, COSMOS and AEGIS fields in the 0.5-8 keV band. The X-ray sources span the redshift interval 0<z<3 and have a median value of Med{z}=0.976.We employ the projected two-point correlation function to infer the spatial clustering and find a clustering length of r0= 7.2+-0.6 h^{-1} Mpc and a slope of \gamma=1.48+-0.12, which corresponds to a bias of b=2.26+-0.16. Using two different halo bias models, we consistently estimate an average dark-matter host halo mass of Mh\sim 1.3 (+-0.3) x 10^{13} h^{-1} M_sun. The X-ray AGN bias and the corresponding dark-matter host halo mass, are significantly higher than the corresponding values of optically selected AGN (at the same redshifts). %indicating different populations of AGN. The redshift evolution of the X-ray selected AGN bias indicates, in agreement with other recent studies, that a unique dark-matter halo mass does not fit well the bias at all the different redshifts probed. Furthermore, we investigate if there is a dependence of the clustering strength on X-ray luminosity. To this end we consider only 650 sources around z~1 and we apply a procedure to disentangle the dependence of clustering on redshift. We find indications for a positive dependence of the clustering length on X-ray luminosity, in the sense that the more luminous sources have a larger clustering length and hence a higher dark-matter halo mass. In detail we find for an average luminosity difference of \delta\log_{10} L_x ~ 1 a halo mass difference of a factor of ~3. These findings appear to be consistent with a galaxy-formation model where the gas accreted onto the supermassive black hole in intermediate luminosity AGN comes mostly from the hot-halo atmosphere around the host galaxy.

We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.

We present the results of a narrow-band (H2S1, 2.121/0.021um) imaging search with WFCAM/UKIRT for H-alpha emitters around several potential signposts of rare (10^-7-10^-8 Mpc^-3) over-dense regions at z=2.23: an over-density of QSOs (2QZ cluster), a powerful, high-redshift radio galaxy (HzRG), and a concentration of submillimetre galaxies (SMGs) and optically faint radio galaxies (OFRGs). In total, we detect 137 narrow-band emitter candidates down to emission-line fluxes of 0.5-1 x 10^-16 erg s^-1 cm^-2, across a total area of 0.56 sq. degrees (2.1 x 10^5 comoving Mpc at z=2.23) in these fields. The BzK colours of the emitters suggest that at least 80% of our sample are likely to be H-alpha emitters (HAEs) at z=2.23. This is one of the largest HAE samples known at z>2. We find modest (~3-sigma) local over-densities of emitters associated with all the three targets. In the 2QZ cluster field, the emitters show a striking filamentary structure connecting four of the z=2.23 QSOs extending over 30Mpc (comoving). In the HzRG and SMG/OFRG fields, the structures appear to be smaller and seen only in the vicinities of the targets. The K-band magnitudes and the H-alpha equivalent widths of the emitters are weakly correlated with the over-density of the emitters: emitters in over-dense region are more evolved systems compared to those in under-dense regions at z=2.23. We find several examples of extended HAEs in our target fields, including a striking example with a spatial extent of 7.5 arcsec (60 kpc at z=2.23) in the 2QZ field, suggesting that these are relatively common in high-density regions. We conclude that narrow-band H-alpha surveys are efficient routes to map over-dense regions at high-z and thus to understand the relation between the growth of galaxies and their surrounding large-scale structures.

We use a coupled model of the formation and evolution of galaxies and black holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold dark matter universe. The model predicts the BH mass, spin and mass accretion history. BH mass grows via accretion triggered by discs becoming dynamically unstable or galaxy mergers (called the starburst mode) and accretion from quasi-hydrostatic hot gas haloes (called the hot-halo mode). By taking into account AGN obscuration, we obtain a very good fit to the observed luminosity functions (LF) of AGN (optical, soft and hard X-ray, and bolometric) for a wide range of redshifts (0<z<6). The model predicts a hierarchical build up of BH mass, with the typical mass of actively growing BHs increasing with decreasing redshift. Remarkably, despite this, we find downsizing in the AGN population, in terms of the differential growth with redshift of the space density of faint and bright AGN. This arises naturally from the interplay between the starburst and hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo mode. The bright end of the LF, on the other hand, is dominated by AGN which host BHs accreting close to or in excess of the Eddington limit during the starburst mode. The model predicts that the comoving space density of AGN peaks at z~3, similar to the star formation history. However, when taking into account obscuration, the space density of faint AGN peaks at lower redshift (z<2) than that of bright AGN (z~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.

We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical model which simulates the formation and evolution of galaxies in a cold dark matter (CDM) universe. During the evolution of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, and the resulting accretion power regulates the gas cooling and subsequent star formation. The accretion flow is assumed to form a geometrically thin cool disc when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick, radiatively inefficient hot flow when the accretion rate falls below this value. The resulting quasar optical luminosity function matches observations very well, and the mass of the SMBH correlates with the mass of the galaxy bulge as observed. The BH spin distribution depends strongly on whether the gas in any given accretion episode remains in the same plane (prolonged accretion) or whether, due to self-gravity, it fragments into multiple, randomly aligned accretion episodes (chaotic accretion). In the chaotic accretion model there is a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun) are hosted by giant elliptical galaxies and are rapidly spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Using the Blandford-Znajek mechanism for jet production to calculate the jet power, our model is able to reproduce the radio loudness of radio galaxies, LINERS and Seyferts. This is the first confirmation that a CDM galaxy formation model can reproduce the observed phenomenology of AGN.