Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations.

Strong scaling relations between host galaxy properties (such as stellar mass, bulge mass, luminosity, effective radius etc) and their nuclear supermassive black hole's mass point towards a close co-evolution. In this work, we first review previous efforts supporting the fundamental importance of the relation between supermassive black hole mass and stellar velocity dispersion ($M_{\rm BH}$-$\sigma_{\rm e}$). We then present further original work supporting this claim via analysis of residuals and principal component analysis applied to some among the latest compilations of local galaxy samples with dynamically measured supermassive black hole masses. We conclude with a review of the main physical scenarios in favour of the existence of a $M_{\rm BH}$-$\sigma_{\rm e}$ relation, with a focus on momentum-driven outflows.

In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey.

The AGN bolometric correction is a key element to understand BH demographics and compute accurate BH accretion histories from AGN luminosities. However, current estimates still differ from each other by up to a factor of two to three, and rely on extrapolations at the lowest and highest luminosities. Here we revisit this fundamental issue presenting general hard X-ray ($K_{X}$) and optical ($K_{O}$) bolometric corrections, computed combining several AGN samples spanning the widest (about 7 dex) luminosity range ever used for this kind of studies. We analysed a total of $\sim 1000$ type 1 and type 2 AGN for which a dedicated SED-fitting has been carried out. We provide a bolometric correction separately for type 1 and type 2 AGN; the two bolometric corrections results to be in agreement in the overlapping luminosity range and therefore, for the first time, a universal bolometric correction for the whole AGN sample (both type 1 and type 2) has been computed. We found that $K_{X}$ is fairly constant at $log(L_{BOL}/L_{\odot}) < 11$, while it increases up to about one order of magnitude at $log(L_{BOL}/L_{\odot}) \sim 14.5$. A similar increasing trend has been observed when its dependence on either the Eddington ratio or the BH mass is considered, while no dependence on redshift up to $z\sim3.5$ has been found. On the contrary, the optical bolometric correction appears to be fairly constant (i.e. $K_{O} \sim 5$) whatever is the independent variable. We also verified that our bolometric corrections correctly predict the AGN bolometric luminosity functions. According to this analysis, our bolometric corrections can be applied to the whole AGN population in a wide range of luminosity and redshift.

There exist discrepancies in measurements of the number and evolution of galaxy pairs. The pair fraction appears to be sensitive to both the criteria used to select pair fraction and the methods used to analyze survey data. This paper explores the connection between stellar mass estimation and the pair fraction of galaxies making use of STEEL, the Statistical sEmi-Emprical modeL. Previous results have found the pair fraction is sensitive to choices made when selecting what qualifies as a pair, for example luminosity or stellar mass selections. We find that different estimations of stellar mass such as photometric choice mass-to-light ratio or IMF that effect the stellar mass function also significantly affect the derived galaxy pair fraction. By making use of the galaxy halo connection we investigate these systematic affects on the pair fraction. We constrain the galaxy halo connection using the stellar-mass-halo-mass relationship for two observed stellar mass functions, and the Illustris TNG stellar mass function. Furthermore, we also create a suite of toy models where the stellar-mass-halo-mass relationship is manually changed. For each stellar-mass-halo-mass relation the pair fraction, and its evolution, are generated. We find that enhancements to the number density of high mass galaxies cause steepening of the stellar-mass-halo mass relation, resulting in a reduction of the pair fraction. We argue this is a considerable cause of bias that must be accounted for when comparing pair fractions.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings.

We investigate the cosmic evolution of the ratio between black hole mass (MBH) and host galaxy total stellar mass (Mstellar) out to z~2.5 for a sample of 100 X-ray-selected moderate-luminosity, broad-line active galactic nuclei (AGNs) in the Chandra-COSMOS Legacy Survey. By taking advantage of the deep multi-wavelength photometry and spectroscopy in the COSMOS field, we measure in a uniform way the galaxy total stellar mass using a SED decomposition technique and the black hole mass based on broad emission line measurements and single-epoch virial estimates. Our sample of AGN host galaxies has total stellar masses of 10^10-12Msun, and black hole masses of 10^7.0-9.5Msun. Combining our sample with the relatively bright AGN samples from the literature, we find no significant evolution of the MBH-Mstellar relation with black hole-to-host total stellar mass ratio of MBH/Mstellar~0.3% at all redshifts probed. We conclude that the average black hole-to-host stellar mass ratio appears to be consistent with the local value within the uncertainties, suggesting a lack of evolution of the MBH-Mstellar relation up to z~2.5.

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.

A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes.

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $\epsilon \sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).

Aims. We perform clustering measurements of 800 X-ray selected Chandra COSMOS Legacy (CCL) Type 2 AGN with known spectroscopic redshift to probe the halo mass dependence on AGN host galaxy properties, such as galaxy stellar mass Mstar, star formation rate (SFR) and specific black hole accretion rate (BHAR), in the redshift range z = [0 - 3]. Methods. We split the sample of AGN with known spectroscopic redshits according to Mstar, SFR and specific BHAR, while matching the distributions in terms of the other parameters, including redshift. We measure the projected two-point correlation function wp(rp) and model it with the 2-halo term to derive the large-scale bias b and the corresponding typical mass of the hosting halo, for the different subsamples. Results. We found no significant dependence of the large-scale bias and typical halo mass on galaxy stellar mass and specific BHAR for CCL Type 2 AGN at mean z~1, while a negative dependence on SFR is observed, with lower SFR AGN residing in richer environment. Mock catalogs of AGN matched to have the same X-ray luminosity, stellar mass, specific BHAR and SFR of CCL Type 2 AGN, almost reproduce the observed Mstar, specific BHAR and SFR-Mh relations, when assuming a fraction of satellite AGN of 15%, which corresponds to a ratio between the probabilities of satellite and central AGN of being active Q = 2. Mock matched normal galaxies follow a slightly steeper Mstar -Mh relation with low mass mock galaxies residing in less massive halos than mock AGN of similar mass, and are less biased than mock AGN with similar specific BHAR and SFR, at least for Q > 1.

Recent advances in observations have provided a wealth of measurements of the expansions of outflows in galactic discs out to large radii in a variety of galactic hosts. To provide an updated baseline for the interpretation of such data, and to assess to what extent the present status of the modeling is consistent with the existing observations, we provide a compact two-dimensional description for the expansion of AGN-driven shocks in realistic galactic discs with exponential gas density profiles in a disc geometry. We derive solutions for the outflow expansion and the mass outflow rates in different directions with respect to the plane of the disc. These are expressed in terms of the global properties of the host galaxy and of the central AGN to allow for an easy and direct comparison with existing observations in a variety of galactic hosts with measured properties, and out to distances $\sim 10$ kpc from the centre. The results are compared with a state-of-the-art compilation of observed outflows in 19 galaxies with different measured gas and dynamical mass, allowing for a detailed, one-by-one comparison with the model predictions. The agreement we obtain for a wide range of host galaxy gas mass and AGN bolometric luminosity provides a quantitative systematic test for the modeling of AGN-driven outflows in galactic discs. We also consider a larger sample of galaxies with no reliable measurements of the gas and dynamical mass. In this case we perform a comparison of the model predictions for different bins of AGN luminosities assuming values for the gas mass and dynamical mass derived from scaling relations. The encouraging, quantitative agreement of the model predictions with a wide set of existing observations constitutes a baseline for the interpretation of forthcoming data, and for a more detailed treatment of AGN feedback in galaxy formation models.

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models.

We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass functions the satellite mass function is progressively built in time via abundance matching techniques constrained by number densities of centrals in the local Universe. By enforcing dynamical merging timescales as predicted by high-resolution N-body simulations, we obtain satellite distributions as a function of stellar mass and halo mass consistent with current data. We show that stellar stripping, star formation, and quenching play all a secondary role in setting the number densities of massive satellites above $M_*\gtrsim 3\times 10^{10}\, M_{\odot}$. We further show that observed star formation rates used in our empirical model over predict low-mass satellites below $M_*\lesssim 3\times 10^{10}\, M_{\odot}$, whereas, star formation rates derived from a continuity equation approach yield the correct abundances similar to previous results for centrals.

Which mechanism(s) are mainly driving nuclear activity in the centres of galaxies is a major unsettled question. In this study, we investigate the statistical relevance of galaxy mergers for fuelling gas onto the central few kpc of a galaxy, potentially resulting in an active galactic nucleus (AGN). To robustly address that, we employ large-scale cosmological hydrodynamic simulations from the Magneticum Pathfinder set, including models for BH accretion and AGN feedback. Our simulations predict that for luminous AGN ($L_{\rm AGN} > 10^{45} {\rm erg/s}$) at $z = 2$, more than 50 per cent of their host galaxies have experienced a merger in the last 0.5~Gyr. These high merger fractions, however, merely reflect the intrinsically high merger fractions of massive galaxies at $z=2$, in which luminous AGN preferentially occur. Apart from that, our simulations suggest that merger events are not the statistically dominant fuelling mechanism for nuclear activity over a redshift range $z=0-2$: irrespective of AGN luminosity, less than 20 per cent of AGN hosts have on average undergone a recent merger, in agreement with a number of observational studies. The central ISM conditions required for inducing AGN activity can be, but are not necessarily caused by a merger. Despite the statistically minor relevance of mergers, at a given AGN luminosity and stellar mass, the merger fractions of AGN hosts can be by up to three times higher than that of inactive galaxies. Such elevated merger fractions still point towards an intrinsic connection between AGN and mergers, consistent with our traditional expectation.

We analyze 168 Swift monitoring observations of the nearest broad absorption line quasar Mrk 231 in the UV and X-ray bands, where we detect significant variability in the UV ($\sim$2246\AA) light curve with a null probability of $4.3\times10^{-10}$ for a constant model. Separately, from an archival sample of Swift observed active galactic nuclei (AGN), we measure the relation between UV excess variance and luminosity, finding that the normalized UV excess variance decreases with luminosity. Comparing to this mean relation, the normalized UV excess variance of Mrk 231 is smaller, however within the scatter characterising the full population. The upper limit of the X-ray excess variance is consistent with other AGN. The power spectrum density of the UV light curve can be well fit by a power law model with a slope of $1.82\pm0.14$ between $10^{-7.5}$ and $10^{-6}$ Hz, consistent with those for typical AGN, with no obvious quasi-periodical oscillation peaks. The UV variability and its power spectrum suggest that a significant amount of the UV emission of Mrk 231 is from the accretion disk. The consistencies in the normalized UV variability and the shape of the power spectrum density between Mrk 231 and other normal AGN suggest that the origin of UV variability of broad absorption line quasars is similar to other AGN, and dust scattering at large scales such as the torus is not a dominating process for the UV emission of Mrk 231. Significant scattering, if present, is constrained to smaller than $\sim$10 light days. We perform lagged correlation analysis between the UV and X-ray light curves and find the correlation insignificant within the present data.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

The active galaxy NGC 1275 lies at the center of the Perseus cluster of galaxies, being an archetypal BH-galaxy system that is supposed to fit well with the M_BH-host scaling relations obtained for quiescent galaxies. Since it harbours an obscured AGN, only recently our group has been able to estimate its black hole mass. Here our aim is to pinpoint NGC 1275 on the less dispersed scaling relations, namely the M_BH-sigma_star and M_BH-L_bul planes. Starting from our previous work Ricci et al. 2017b, we estimate that NGC 1275 falls well outside the intrinsic dispersion of the M_BH-sigma_star plane being ~1.2 dex (in black hole mass) displaced with respect to the scaling relations. We then perform a 2D morphological decomposition analysis on Spitzer/IRAC images at 3.6 mic and find that, beyond the bright compact nucleus that dominates the central emission, NGC 1275 follows a de Vaucouleurs profile with no sign of significant star formation nor clear merger remnants. Nonetheless, its displacement on the M_BH-L_(3.6,bul) plane with respect to the scaling relation is as high as observed in the M_BH-sigma_star.

Within a galaxy the stellar mass-to-light ratio $\Upsilon_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $\Upsilon_*$ gradients than if the IMF is held fixed. If $\Upsilon_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $\Upsilon_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $\Upsilon_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $\Upsilon_*$ gradients changes the high-mass slope of the stellar mass function $\phi(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $\Upsilon_*$ gradients in larger samples.

[ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early fast evolution comprises: biased collapse of the low angular momentum gaseous baryons located in the inner regions of the host dark matter halo; cooling, fragmentation, and infall of the gas down to the radius set by the centrifugal barrier; further rapid compaction via clump/gas migration toward the galaxy center, where strong heavily dust-enshrouded star-formation takes place and most of the stellar mass is accumulated; ejection of substantial gas amount from the inner regions by feedback processes, which causes a dramatic puffing up of the stellar component. In the late slow evolution, passive aging of stellar populations and mass additions by dry merger events occur. We describe these processes relying on prescriptions inspired by basic physical arguments and by numerical simulations, to derive new analytical estimates of the relevant sizes, timescales, and kinematic properties for individual galaxies along their evolution. Then we obtain quantitative results as a function of galaxy mass and redshift, and compare them to recent observational constraints on half-light size $R_e$, on the ratio $v/\sigma$ between rotation velocity and velocity dispersion (for gas and stars) and on the specific angular momentum $j_\star$ of the stellar component; we find good consistency with the available multi-band data in average values and dispersion, both for local ETGs and for their $z\sim 1-2$ star-forming and quiescent progenitors.

We present the first data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 103 nearby galaxies from the Palomar sample, observed with the eMERLIN array, as part of the LeMMINGs survey. This sample includes galaxies which are active (LINER and Seyfert) and quiescent (HII galaxies and Absorption line galaxies, ALG), which are reclassified based upon revised emission-line diagrams. We detect radio emission $\gtrsim$ 0.2 mJy for 47/103 galaxies (22/34 for LINERS, 4/4 for Seyferts, 16/51 for HII galaxies and 5/14 for ALGs) with radio sizes typically of $\lesssim$100 pc. We identify the radio core position within the radio structures for 41 sources. Half of the sample shows jetted morphologies. The remaining half shows single radio cores or complex morphologies. LINERs show radio structures more core-brightened than Seyferts. Radio luminosities of the sample range from 10$^{32}$ to 10$^{40}$ erg s$^{-1}$: LINERs and HII galaxies show the highest and the lowest radio powers respectively, while ALGs and Seyferts have intermediate luminosities. We find that radio core luminosities correlate with black hole (BH) mass down to $\sim$10$^{7}$ M$_{\odot}$, but a break emerges at lower masses. Using [O III] line luminosity as a proxy for the accretion luminosity, active nuclei and jetted HII galaxies follow an optical fundamental plane of BH activity, suggesting a common disc-jet relationship. In conclusion, LINER nuclei are the scaled-down version of FR I radio galaxies; Seyferts show less collimated jets; HII galaxies may host weak active BHs and/or nuclear star-forming cores; and recurrent BH activity may account for ALG properties.

Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $\sigma_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $\sigma$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm \alpha_{JAM}}$, from $M_{*}$ and $\sigma_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,\sigma_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(\sigma_R/\sigma_e) = (R/R_e)^{-\gamma(n)}$. The resulting stellar mass functions, based on $M_*^{\rm \alpha_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm \alpha_{JAM}}/M_* - \sigma_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $\phi(M_*^{\rm \alpha_{JAM}})$ is an estimate of the stellar mass function in which $\sigma_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $\sigma_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form.

We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and S\'{e}rsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be S\'{e}rsic and the input S\'{e}rsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing S\'{e}rsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution.

We present QSFit, a new software package to automatically perform the analysis of Active Galactic Nuclei (AGN) optical spectra. The software provides luminosity estimates for the AGN continuum, the Balmer continuum, both optical and UV iron blended complex, host galaxy and emission lines, as well as width, velocity offset and equivalent width of 20 emission lines. Improving on a number of previous studies on AGN spectral analysis, QSFit fits all the components simultaneously, using an AGN continuum model which extends over the entire available spectrum, and is thus a probe of the actual AGN continuum whose estimates are scarcely influenced by localized features (e.g. emission lines) in the spectrum. We used QSFit to analyze 71,251 optical spectra of Type 1 AGN at z < 2 (obtained by the Sloan Digital Sky Survey, SDSS) and to produce a publicly available catalog of AGN spectral properties. Such catalog allowed us (for the first time) to estimate the AGN continuum slope and the Balmer continuum luminosity on a very large sample, and to show that there is no evident correlation between these quantities the redshift. All data in the catalog, the plots with best fitting model and residuals, and the IDL code we used to perform the analysis, are available on a dedicated website. The whole fitting process is customizable for specific needs, and can be extended to analyze spectra from other data sources. The ultimate purpose of QSFit is to allow astronomers to run standardized recipes to analyze the AGN data, in a simple, replicable and shareable way.

Type 2 active galactic nuclei (AGN) represent the majority of the AGN population. However, due to the difficulties in measuring their black hole (BH) masses, it is still unknown whether they follow the same BH mass-host galaxy scaling relations valid for quiescent galaxies and type 1 AGN. Here we present the locus of type 2 AGN having virial BH mass estimates in the $M_{BH}-\sigma_\star$ plane. Our analysis shows that the BH masses of type 2 AGN are $\sim0.9$ dex smaller than type 1 AGN at $\sigma_\star\sim 185$ km s$^{-1}$, regardless of the (early/late) AGN host galaxy morphology. Equivalently, type 2 AGN host galaxies have stellar velocity dispersions $\sim 0.2$ dex higher than type 1 AGN hosts at $M_{BH}\sim10^7$ M$_\odot$.

Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by using a comprehensive galaxy evolution semi-analytic model, which self-consistently evolves supermassive black holes from high-redshift seeds via gas accretion and mergers, and also includes AGN feedback. Once tuned to reproduce the (mean) correlation of black-hole mass with velocity dispersion, the model is unable to also account for the correlation with stellar mass. This behaviour is independent of the model's parameters, thus suggesting an internal inconsistency in the data. The predicted distributions, especially at the low-mass end, are also much broader than observed. However, if selection effects are included, the model's predictions tend to align with the observations. We also demonstrate that the correlations between the residuals of the local scaling relations are more effective than the scaling relations themselves at constraining AGN feedback models. In fact, we find that our semi-analytic model, while in apparent broad agreement with the scaling relations when accounting for selection biases, yields very weak correlations between their residuals at fixed stellar mass, in stark contrast with observations. This problem persists when changing the AGN feedback strength, and is also present in the $z\sim 0$ outputs of the hydrodynamic cosmological simulation Horizon-AGN, which includes state-of-the-art treatments of AGN feedback. This suggests that current AGN feedback models may be too weak or are simply not capturing the effect of the black hole on the stellar velocity dispersion.

The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from different treatments of the stellar population remain of order < 0.5 dex, systematics in photometry are now about 0.1 dex, despite recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey (SDSS) pipeline photometry leads to underestimates of rho_*(> M_*) by factors of 3-10 in the mass range 10^11 - 10^11.6 M_Sun, but up to a factor of 100 at higher stellar masses. This impacts studies which match massive galaxies to dark matter halos. Although systematics which arise from different treatments of the stellar population remain of order < 0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low and high redshift samples are compared.

Feedback from accreting SMBHs is often identified as the main mechanism responsible for regulating star-formation in AGN host galaxies. However, the relationships between AGN activity, radiation, winds, and star-formation are complex and still far from being understood. We study scaling relations between AGN properties, host galaxy properties and AGN winds. We then evaluate the wind mean impact on the global star-formation history, taking into account the short AGN duty cycle with respect to that of star-formation. We first collect AGN wind observations for 94 AGN with detected massive winds at sub-pc to kpc spatial scales. We then fold AGN wind scaling relations with AGN luminosity functions, to evaluate the average AGN wind mass-loading factor as a function of cosmic time. We find strong correlations between the AGN molecular and ionised wind mass outflow rates and the AGN bolometric luminosity. The power law scaling is steeper for ionised winds (slope 1.29+/-0.38) than for molecular winds (0.76+/-0.06), meaning that the two rates converge at high bolometric luminosities. The molecular gas depletion timescale and the molecular gas fraction of galaxies hosting powerful AGN winds are 3-10 times shorter and smaller than those of main-sequence galaxies with similar SFR, stellar mass and redshift. These findings suggest that, at high AGN bolometric luminosity, the reduced molecular gas fraction may be due to the destruction of molecules by the wind, leading to a larger fraction of gas in the atomic ionised phase. The AGN wind mass-loading factor $\eta=\dot M_{OF}/SFR$ is systematically higher than that of starburst driven winds. Our analysis shows that AGN winds are, on average, powerful enough to clean galaxies from their molecular gas only in massive systems at z<=2, i.e. a strong form of co-evolution between SMBHs and galaxies appears to break down for the least massive galaxies.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and S\'ersic indices of the bulges of their host galaxies, have dynamical mass and S\'ersic index distributions consistent with those of SDSS early-type galaxies of similar bulge stellar mass. The host galaxies of local black hole samples thus do not appear structurally different from SDSS galaxies, sharing similar dynamical masses, light profiles and light distributions. Analysis of the residuals reveals that velocity dispersion is more fundamental than S\'ersic index n in the scaling relations between black holes and galaxies. Indeed, residuals with S\'ersic index could be ascribed to the (weak) correlation with bulge mass or even velocity dispersion. Finally, targetted Monte Carlo simulations that include the effects of the sphere of influence of the black hole, and tuned to reproduce the observed residuals and scaling relations in terms of velocity dispersion and stellar mass, show that, at least for galaxies with Mbulge > 1e10 Msun and n>5, the observed mean black hole mass at fixed S\'ersic index is biased significantly higher than the intrinsic value.

We compare the set of local galaxies having dynamically measured black holes with a large, unbiased sample of galaxies extracted from the Sloan Digital Sky Survey. We confirm earlier work showing that the majority of black hole hosts have significantly higher velocity dispersions sigma than local galaxies of similar stellar mass. We use Monte-Carlo simulations to illustrate the effect on black hole scaling relations if this bias arises from the requirement that the black hole sphere of influence must be resolved to measure black hole masses with spatially resolved kinematics. We find that this selection effect artificially increases the normalization of the Mbh-sigma relation by a factor of at least ~3; the bias for the Mbh-Mstar relation is even larger. Our Monte Carlo simulations and analysis of the residuals from scaling relations both indicate that sigma is more fundamental than Mstar or effective radius. In particular, the Mbh-Mstar relation is mostly a consequence of the Mbh-sigma and sigma-Mstar relations, and is heavily biased by up to a factor of 50 at small masses. This helps resolve the discrepancy between dynamically-based black hole-galaxy scaling relations versus those of active galaxies. Our simulations also disfavour broad distributions of black hole masses at fixed sigma. Correcting for this bias suggests that the calibration factor used to estimate black hole masses in active galaxies should be reduced to values of fvir~1. Black hole mass densities should also be proportionally smaller, perhaps implying significantly higher radiative efficiencies/black hole spins. Reducing black hole masses also reduces the gravitational wave signal expected from black hole mergers.

Galaxy interactions are thought to be one of the main triggers of Active Galactic Nuclei (AGN), especially at high luminosities, where the accreted gas mass during the AGN lifetime is substantial. Evidence for a connection between mergers and AGN, however, remains mixed. Possible triggering mechanisms remain particularly poorly understood for luminous AGN, which are thought to require triggering by major mergers, rather than secular processes. We analyse the host galaxies of a sample of 20 optically and X-ray selected luminous AGN (log($L_{bol}$ [erg/s]) $>$ 45) at z $\sim$ 0.6 using HST WFC3 data in the F160W/H band. 15/20 sources have resolved host galaxies. We create a control sample of mock AGN by matching the AGN host galaxies to a control sample of non-AGN galaxies. Visual signs of disturbances are found in about 25% of sources in both the AGN hosts and control galaxies. Using both visual classification and quantitative morphology measures, we show that the levels of disturbance are not enhanced when compared to a matched control sample. We find no signs that major mergers play a dominant role in triggering AGN at high luminosities, suggesting that minor mergers and secular processes dominate AGN triggering up to the highest AGN luminosities. The upper limit on the enhanced fraction of major mergers is $\leqslant$20%. While major mergers might increase the incidence of (luminous AGN), they are not the prevalent triggering mechanism in the population of unobscured AGN.

The cosmological process of hydrogen (HI) reionization in the intergalactic medium is thought to be driven by UV photons emitted by star-forming galaxies and ionizing active galactic nuclei (AGN). The contribution of QSOs to HI reionization at $z>4$ has been traditionally believed to be quite modest. However, this view has been recently challenged by new estimates of a higher faint-end UV luminosity function (LF). To set firmer constraints on the emissivity of AGN at $z<6$, we here make use of complete X-ray selected samples including deep Chandra and new COSMOS data, capable to efficiently measure the 1 ryd comoving AGN emissivity up to $z\sim5-6$ and down to five magnitudes fainter than probed by current optical surveys, without any luminosity extrapolation. We find good agreement between the logN$\rm{_H}\lesssim21-22$ cm$^{-2}$ X-ray LF and the optically-selected QSO LF at all redshifts for $M_{1450}\leq -23$. The full range of the logN$\rm{_H}\lesssim21-22$ cm$^{-2}$ LF ($M_{1450} \leq -17$) was then used to quantify the contribution of AGN to the photon budget critical value needed to keep the Universe ionized. We find that the contribution of ionizing AGN at $z = 6$ is as small as 1\% - 7\%, and very unlikely to be greater than 30\%, thus excluding an AGN-dominated reionization scenario.

The time lag between optical and near-infrared continuum emission in active galactic nuclei (AGN) shows a tight correlation with luminosity and has been proposed as a standardisable candle for cosmology. In this paper, we explore the use of these AGN hot-dust time lags for cosmological model fitting under the constraints of the new VISTA Extragalactic Infrared Legacy Survey VEILS. This new survey will target a 9 deg^2 field observed in J- and Ks-band with a 14-day cadence and will run for three years. The same area will be covered simultaneously in the optical griz bands by the Dark Energy Survey, providing complementary time-domain optical data. We perform realistic simulations of the survey setup, showing that we expect to recover dust time lags for about 450 objects out of a total of 1350 optical type 1 AGN, spanning a redshift range of 0.1 < z < 1.2. We use the lags recovered from our simulations to calculate precise distance moduli, establish a Hubble diagram, and fit cosmological models. Assuming realistic scatter in the distribution of the dust around the AGN as well as in the normalisation of the lag-luminosity relation, we are able to constrain {\Omega}_{\Lambda} in {\Lambda}CDM with similar accuracy as current supernova samples. We discuss the benefits of combining AGN and supernovae for cosmology and connect the present work to future attempts to reach out to redshifts of z > 4.

We present the measurement of the projected and redshift space 2-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy AGN at 2.9$\leq$z$\leq$5.5 ($\langle L_{bol} \rangle \sim$10$^{46}$ erg/s) using the generalized clustering estimator based on phot-z probability distribution functions (Pdfs) in addition to any available spec-z. We model the projected 2pcf estimated using $\pi_{max}$ = 200 h$^{-1}$ Mpc with the 2-halo term and we derive a bias at z$\sim$3.4 equal to b = 6.6$^{+0.60}_{-0.55}$, which corresponds to a typical mass of the hosting halos of log M$_h$ = 12.83$^{+0.12}_{-0.11}$ h$^{-1}$ M$_{\odot}$. A similar bias is derived using the redshift-space 2pcf, modelled including the typical phot-z error $\sigma_z$ = 0.052 of our sample at z$\geq$2.9. Once we integrate the projected 2pcf up to $\pi_{max}$ = 200 h$^{-1}$ Mpc, the bias of XMM and \textit{Chandra} COSMOS at z=2.8 used in Allevato et al. (2014) is consistent with our results at higher redshift. The results suggest only a slight increase of the bias factor of COSMOS AGN at z$\gtrsim$3 with the typical hosting halo mass of moderate luminosity AGN almost constant with redshift and equal to logM$_h$ = 12.92$^{+0.13}_{-0.18}$ at z=2.8 and log M$_h$ = 12.83$^{+0.12}_{-0.11}$ at z$\sim$3.4, respectively. The observed redshift evolution of the bias of COSMOS AGN implies that moderate luminosity AGN still inhabit group-sized halos at z$\gtrsim$3, but slightly less massive than observed in different independent studies using X-ray AGN at z$\leq2$.

There is still much debate surrounding how the most massive, central galaxies in the local universe have assembled their stellar mass, especially the relative roles of in-situ growth versus later accretion via mergers. In this paper, we set firmer constraints on the evolutionary pathways of the most massive central galaxies by making use of empirical estimates on their abundances and stellar ages. The most recent abundance matching and direct measurements strongly favour that a substantial fraction of massive galaxies with Mstar>3x10^11 Msun reside at the centre of clusters with mass Mhalo>3x10^13 Msun. Spectral analysis supports ages >10 Gyrs, corresponding to a formation redshift z_form >2. We combine these two pieces of observationally-based evidence with the mass accretion history of their host dark matter haloes. We find that in these massive haloes, the stellar mass locked up in the central galaxy is comparable to, if not greater than, the total baryonic mass at z_form. These findings indicate that either only a relatively minor fraction of their present-day stellar mass was formed in-situ at z_form, or that these massive, central galaxies form in the extreme scenario where almost all of the baryons in the progenitor halo are converted into stars. Interestingly, the latter scenario would not allow for any substantial size growth since the galaxy's formation epoch either via mergers or expansion. We show our results hold irrespective of systematic uncertainties in stellar mass, abundances, galaxy merger rates, stellar initial mass function, star formation rate and dark matter accretion histories.

We present the largest high-redshift (3<z<6.85) sample of X-ray-selected active galactic nuclei (AGN) on a contiguous field, using sources detected in the Chandra COSMOS Legacy survey. The sample contains 174 sources, 87 with spectroscopic redshift, the other 87 with photometric redshift (z_phot). In this work we treat z_phot as a probability weighted sum of contributions, adding to our sample the contribution of sources with z_phot<3 but z_phot probability distribution >0 at z>3. We compute the number counts in the observed 0.5-2 keV band, finding a decline in the number of sources at z>3 and constraining phenomenological models of X-ray background. We compute the AGN space density at z>3 in two different luminosity bins. At higher luminosities (logL(2-10 keV) > 44.1 erg/s) the space density declines exponentially, dropping by a factor ~20 from z~3 to z~6. The observed decline is ~80% steeper at lower luminosities (43.55 erg/s < logL(2-10 keV) < 44.1 erg/s), from z~3 to z~4.5. We study the space density evolution dividing our sample in optically classified Type 1 and Type 2 AGN. At logL(2-10 keV) > 44.1 erg/s, unobscured and obscured objects may have different evolution with redshift, the obscured component being three times higher at z~5. Finally, we compare our space density with predictions of quasar activation merger models, whose calibration is based on optically luminous AGN. These models significantly overpredict the number of expected AGN at logL(2-10 keV) > 44.1 erg/s with respect to our data.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

Recent studies have established that extreme dwarf galaxies --whether satellites or field objects-- suffer from the so called "too big to fail" (TBTF) problem. Put simply, the TBTF problem consists of the fact that it is difficult to explain both the measured kinematics of dwarfs and their observed number density within the LCDM framework. The most popular proposed solutions to the problem involve baryonic feedback processes. For example, reionization and baryon depletion can decrease the abundance of halos that are expected to host dwarf galaxies. Moreover, feedback related to star formation can alter the dark matter density profile in the central regions of low-mass halos. In this article we assess the TBTF problem for field dwarfs, taking explicitly into account the baryonic effects mentioned above. We find that 1) reionization feedback cannot resolve the TBTF problem on its own, because the halos in question are too massive to be affected by it, and that 2) the degree to which profile modification can be invoked as a solution to the TBTF problem depends on the radius at which galactic kinematics are measured. Based on a literature sample of about 90 dwarfs with interferometric observations in the 21cm line of atomic hydrogen (HI), we conclude that the TBTF problem persists despite baryonic effects. However, the preceding statement assumes that the sample under consideration is representative of the general population of field dwarfs. In addition, the unexplained excess of dwarf galaxies in LCDM could be as small as a factor of ~ 1.8, given the current uncertainties in the measurement of the galactic velocity function. Both of these caveats highlight the importance of upcoming uniform surveys with HI interferometers for advancing our understanding of the issue.

Supermassive black hole -- host galaxy relations are key to the computation of the expected gravitational wave background (GWB) in the pulsar timing array (PTA) frequency band. It has been recently pointed out that standard relations adopted in GWB computations are in fact biased-high. We show that when this selection bias is taken into account, the expected GWB in the PTA band is a factor of about three smaller than previously estimated. Compared to other scaling relations recently published in the literature, the median amplitude of the signal at $f=1$yr$^{-1}$ drops from $1.3\times10^{-15}$ to $4\times10^{-16}$. Although this solves any potential tension between theoretical predictions and recent PTA limits without invoking other dynamical effects (such as stalling, eccentricity or strong coupling with the galactic environment), it also makes the GWB detection more challenging.

We analyzed a large sample of radio-loud and radio-quiet quasar spectra at redshift 1.0 < z < 1.2 to compare the inferred underlying quasar continuum slopes (after removal of the host galaxy contribution) with accretion disk models. The latter predict redder (decreasing) alpha_3000 continuum slopes (L_\nu~\nu^alpha at 3000Ang) with increasing black hole mass, bluer alpha_3000 with increasing luminosity at 3000Ang, and bluer alpha_3000 with increasing spin of the black hole, when all other parameters are held fixed. We find no clear evidence for any of these predictions in the data. In particular we find that: (i) alpha_3000 shows no significant dependence on black hole mass or luminosity. Dedicated Monte Carlo tests suggest that the substantial observational uncertainties in the black hole virial masses can effectively erase any intrinsic dependence of alpha_3000 on black hole mass, in line with some previous studies. (ii) The mean slope alpha_3000 of radio-loud sources, thought to be produced by rapidly spinning black holes, is comparable to, or even redder than, that of radio-quiet quasars. Indeed, although quasars appear to become more radio loud with decreasing luminosity, we still do not detect any significant dependence of alpha_3000 on radio loudness. The predicted mean alpha_3000 slopes tend to be bluer than in the data. Disk models with high inclinations and dust extinction tend to produce redder slopes closer to empirical estimates. Our mean alpha_3000 values are close to the ones independently inferred at z<0.5 suggesting weak evolution with redshift, at least for moderately luminous quasars.

The triggering mechanisms for Active Galactic Nuclei (AGN) are still debated. Some of the most popular ones include galaxy interactions (IT) and disk instabilities (DI). Using an advanced semi analytic model (SAM) of galaxy formation, coupled to accurate halo occupation distribution modeling, we investigate the imprint left by each separate triggering process on the clustering strength of AGN at small and large scales. Our main results are as follows: i) DIs, irrespective of their exact implementation in the SAM, tend to fall short in triggering AGN activity in galaxies at the center of halos with $M_h>10^{13.5} h^{-1}M_{\odot}$. On the contrary, the IT scenario predicts abundance of active, central galaxies that generally agrees well with observations at every halo mass. ii) The relative number of satellite AGN in DIs at intermediate-to-low luminosities is always significantly higher than in IT models, especially in groups and clusters. The low AGN satellite fraction predicted for the IT scenario might suggest that different feeding modes could simultaneously contribute to the triggering of satellite AGN. iii) Both scenarios are quite degenerate in matching large-scale clustering measurements, suggesting that the sole average bias might not be an effective observational constraint. iv) Our analysis suggests the presence of both a mild luminosity and a more consistent redshift dependence in the AGN clustering, with AGN inhabiting progressively less massive dark matter halos as the redshift increases. We also discuss the impact of different observational selection cuts in measuring AGN clustering, including possible discrepancies between optical and X-ray surveys.

We describe the luminosity function, based on Sersic fits to the light profiles, of CMASS galaxies at z ~ 0.55. Compared to previous estimates, our Sersic-based reductions imply more luminous, massive galaxies, consistent with the effects of Sersic- rather than Petrosian or de Vaucouleur-based photometry on the Sloan Digital Sky Survey (SDSS) main galaxy sample at z ~ 0.1. This implies a significant revision of the high mass end of the correlation between stellar and halo mass. Inferences about the evolution of the luminosity and stellar mass functions depend strongly on the assumed, and uncertain, k+e corrections. In turn, these depend on the assumed age of the population. Applying k+e corrections taken from fitting the models of Maraston et al. (2009) to the colors of both SDSS and CMASS galaxies, the evolution of the luminosity and stellar mass functions appears impressively passive, provided that the fits are required to return old ages. However, when matched in comoving number- or luminosity-density, the SDSS galaxies are less strongly clustered compared to their counterparts in CMASS. This rules out the passive evolution scenario, and, indeed, any minor merger scenarios which preserve the rank ordering in stellar mass of the population. Potential incompletenesses in the CMASS sample would further enhance this mismatch. Our analysis highlights the virtue of combining clustering measurements with number counts.

We present a catalog of visual like H-band morphologies of $\sim50.000$ galaxies ($H_{f160w}<24.5$) in the 5 CANDELS fields (GOODS-N, GOODS-S, UDS, EGS and COSMOS). Morphologies are estimated with Convolutional Neural Networks (ConvNets). The median redshift of the sample is $<z>\sim1.25$. The algorithm is trained on GOODS-S for which visual classifications are publicly available and then applied to the other 4 fields. Following the CANDELS main morphology classification scheme, our model retrieves the probabilities for each galaxy of having a spheroid, a disk, presenting an irregularity, being compact or point source and being unclassifiable. ConvNets are able to predict the fractions of votes given a galaxy image with zero bias and $\sim10\%$ scatter. The fraction of miss-classifications is less than $1\%$. Our classification scheme represents a major improvement with respect to CAS (Concentration-Asymmetry-Smoothness)-based methods, which hit a $20-30\%$ contamination limit at high z. The catalog is released with the present paper via the $\href{http://rainbowx.fis.ucm.es/Rainbow_navigator_public}{Rainbow\,database}$

Recent inspections of local available data suggest that the almost linear relation between the stellar mass of spheroids ($M_{\rm sph}$) and the mass of the super massive Black Holes (BHs) residing at their centres, shows a break below $M_{\rm sph} \sim 10^{10}\ {\rm M}_\odot$, with a steeper, about quadratic relation at smaller masses. We investigate the physical mechanisms responsible for the change in slope of this relation, by comparing data with the results of the semi-analytic model of galaxy formation MORGANA, which already predicted such a break in its original formulation. We find that the change of slope is mostly induced by effective stellar feedback in star-forming bulges. The shape of the relation is instead quite insensitive to other physical mechanisms connected to BH accretion such as disc instabilities, galaxy mergers, Active Galactic Nucleus (AGN) feedback, or even the exact modelling of accretion onto the BH, direct or through a reservoir of low angular momentum gas. Our results support a scenario where most stars form in the disc component of galaxies and are carried to bulges through mergers and disc instabilities, while accretion onto BHs is connected to star formation in the spheroidal component. Therefore, a model of stellar feedback that produces stronger outflows in star-forming bulges than in discs will naturally produce a break in the scaling relation. Our results point to a form of co-evolution especially at lower masses, below the putative break, mainly driven by stellar feedback rather than AGN feedback.

[abridged] We investigate the coevolution of galaxies and hosted supermassive black holes throughout the history of the Universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source term to reconstruct the supermassive black hole (BH) mass function from the AGN luminosity functions. Such an approach includes physically-motivated AGN lightcurves tested on independent datasets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disc conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the UV and far-IR luminosity functions associated to star formation in galaxies. These results imply that the buildup of stars and BHs in galaxies occurs via in-situ processes, with dry mergers playing a marginal role at least for stellar masses < 3 10^11 M_sun and BH masses < 10^9 M_sun, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. Finally, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, so further validating our results from the continuity equation.

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.

The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth and assembly of spheroids. However, direct comparison between models and observations is usually prevented by the so-called "progenitor bias", i.e., new galaxies entering the observational selection at later epochs, thus eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z=0 along their main progenitors, and assign as their "descendants" local SDSS central galaxies matched in host halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in stellar mass of a factor of >2 since z~1, and of >2.5 in mean effective radius. We then compare the inferred stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass and size increase of a factor ~2-4 observed in cluster galaxies since z ~1. The predicted, average (major) merger rate since z~1 is in good agreement with the latest observational estimates.

We use the Arecibo legacy fast ALFA (ALFALFA) 21cm survey to measure the number density of galaxies as a function of their rotational velocity, $V_\mathrm{rot,HI}$ (as inferred from the width of their 21cm emission line). Based on the measured velocity function we statistically connect galaxies with their host halo, via abundance matching. In a lambda cold dark matter ($\Lambda$CDM) cosmology, dwarf galaxies are expected to be hosted by halos that are significantly more massive than indicated by the measured galactic velocity; if smaller halos were allowed to host galaxies, then ALFALFA would measure a much higher galactic number density. We then seek observational verification of this predicted trend by analyzing the kinematics of a literature sample of gas-rich dwarf galaxies. We find that galaxies with $V_\mathrm{rot,HI} \lesssim 25$ $\mathrm{km} \, \mathrm{s}^{-1}$ are kinematically incompatible with their predicted $\Lambda$CDM host halos, in the sense that hosts are too massive to be accommodated within the measured galactic rotation curves. This issue is analogous to the "too big to fail" problem faced by the bright satellites of the Milky Way, but here it concerns extreme dwarf galaxies in the field. Consequently, solutions based on satellite-specific processes are not applicable in this context. Our result confirms the findings of previous studies based on optical survey data and addresses a number of observational systematics present in these works. Furthermore, we point out the assumptions and uncertainties that could strongly affect our conclusions. We show that the two most important among them -namely baryonic effects on the abundances of halos and on the rotation curves of halos- do not seem capable of resolving the reported discrepancy.

We present a census of the active black hole population at 1<z<2, by constructing the bivariate distribution function of black hole mass and Eddington ratio, employing a maximum likelihood fitting technique. The study of the active black hole mass function (BHMF) and the Eddington ratio distribution function (ERDF) allows us to clearly disentangle the active galactic nuclei (AGN) downsizing phenomenon, present in the AGN luminosity function, into its physical processes of black hole mass downsizing and accretion rate evolution. We are utilizing type-1 AGN samples from three optical surveys (VVDS, zCOSMOS and SDSS), that cover a wide range of 3 dex in luminosity over our redshift interval of interest. We investigate the cosmic evolution of the AGN population as a function of AGN luminosity, black hole mass and accretion rate. Compared to z = 0, we find a distinct change in the shape of the BHMF and the ERDF, consistent with downsizing in black hole mass. The active fraction or duty cycle of type-1 AGN at z~1.5 is almost flat as a function of black hole mass, while it shows a strong decrease with increasing mass at z=0. We are witnessing a phase of intense black hole growth, which is largely driven by the onset of AGN activity in massive black holes towards z=2. We finally compare our results to numerical simulations and semi-empirical models and while we find reasonable agreement over certain parameter ranges, we highlight the need to refine these models in order to match our observations.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

We investigate, for the first time at z~3, the clustering properties of 189 Type 1 and 157 Type 2 X-ray active galactic nuclei (AGN) of moderate luminosity (log<Lbol> = 45.3 erg/s), with photometric or spectroscopic redshifts in the range 2.2<z<6.8. These samples are based on Chandra and XMM-Newton data in COSMOS. We find that Type 1 and Type 2 COSMOS AGN at z=3 inhabit DMHs with typical mass of logMh = 12.84+0.10/-0.11 and 11.73+0.39/-0.45 Msun/h, respectively. This result requires a drop in the halo masses of Type 1 and 2 COSMOS AGN at z~3 compared to z<2 XMM COSMOS AGN with similar luminosities. Additionally, we infer that unobscured COSMOS AGN at z~3 reside in 10 times more massive halos compared to obscured COSMOS AGN, at 2.6sigma level. This result extends to z~3 that found in COSMOS at z<2, and rules out the picture in which obscuration is purely an orientation effect. A model which assumes that the AGN activity is triggered by major mergers is quite successful in predicting both the low halo mass of COSMOS AGN and the typical mass of luminous SDSS quasars at z~3, with the latter inhabiting more massive halos respect to moderate luminosity AGN. Alternatively we can argue, at least for Type 1 COSMOS AGN, that they are possibly representative of an early phase of fast (i.e. Eddington limited) BH growth induced by cosmic cold flows or disk instabilities. Given the moderate luminosity, these new fast growing BHs have masses of e7-8 Msun at z~3 which might evolve into e8.5-9 Msun mass BHs at z=0. Following our clustering measurements, we argue that this fast BH growth at z~3 in AGN with moderate luminosity occurs in DMHs with typical mass of 6 times e12 Msun/h.

We study of the role of major mergers (mass ratios >1:4) in driving size growth in high-redshift (1<z<2) spheroidal galaxies (SGs) with stellar masses between 10^9.5 MSun and 10^10.7 MSun. This is a largely unexplored mass range at this epoch, containing the progenitors of more massive SGs on which the bulk of the size-evolution literature is based. We visually split our SGs into systems that are relaxed and those that exhibit tidal features indicative of a recent merger. Numerical simulations indicate that, given the depth of our images, only tidal features due to major mergers will be detectable at these epochs (features from minor mergers being too faint), making the disturbed SGs a useful route for estimating major-merger-driven size growth. The disturbed SGs are offset in size from their relaxed counterparts, lying close to the upper envelope of the local size -- mass relation. The mean size ratio of the disturbed SGs to their relaxed counterparts is ~2. Combining this observed size growth with empirical major-merger histories from the literature suggests that the size evolution of a significant fraction (around two-thirds) of SGs in this mass range could be driven by major mergers. If, as is likely, our galaxies are progenitors of more massive M > 10^10.7 MSun SGs at z<1, then major mergers are also likely to play an important role in the size growth of at least some massive SGs in this mass range.

We quantify the systematics in the size-luminosity relation of galaxies in the SDSS main sample which arise from fitting different 1- and 2-component model profiles to the images. In objects brighter than L*, fitting a single Sersic profile to what is really a two-component SerExp system leads to biases: the half-light radius is increasingly overestimated as n of the fitted single component increases; it is also overestimated at B/T ~ 0.6. However, the net effect on the R-L relation is small, except for the most luminous tail, where it curves upwards towards larger sizes. We also study how this relation depends on morphological type. Our analysis is one of the first to use Bayesian-classifier derived weights, rather than hard cuts, to define morphology. Crudely, there appear to be only two relations: one for early-types (Es, S0s and Sa's) and another for late-types (Sbs and Scds). However, closer inspection shows that within the early-type sample S0s tend to be 15% smaller than Es of the same luminosity, and, among faint late types, Sbs are more than 25% smaller than Scds. Neither the early- nor the late-type relations are pure power-laws: both show significant curvature, which we quantify. However, the R-L relations of the bulges of early-types are almost pure power laws. Our analysis confirms that two mass scales are special for early-type galaxies: M* = 3e10 and 2e11 Msun. These same mass scales are also special for late types: there is almost no correlation between R and M* below the former, and almost no late-types above the latter. We also show that the intrinsic scatter around the relation decreases at large luminosity and/or stellar mass; this should provide additional constraints on models of how the most massive galaxies formed.

We show that recently documented trends in galaxy sizes with mass and redshift can be understood in terms of the influence of underlying cosmic evolution; a holistic view which is complimentary to interpretations involving the accumulation of discreet evolutionary processes acting on individual objects. Using standard cosmology theory, supported with results from the Millennium simulations, we derive expected size trends for collapsed cosmic structures, emphasising the important distinction between these trends and the assembly paths of individual regions. We then argue that the observed variation in the stellar mass content of these structures can be understood to first order in terms of natural limitations of cooling and feedback. But whilst these relative masses vary by orders of magnitude, galaxy and host radii have been found to correlate linearly. We explain how these two aspects will lead to galaxy sizes that closely follow observed trends and their evolution, comparing directly with the COSMOS and SDSS surveys. Thus we conclude that the observed minimum radius for galaxies, the evolving trend in size as a function of mass for intermediate systems, and the observed increase in the sizes of massive galaxies, may all be considered an emergent consequence of the cosmic expansion.

There is growing evidence of a higher intrinsic fraction of broad absorption line quasars (BALQSOs) than that obtained in optical surveys, on which most previous X-ray studies of BALQSOs have focused. Here we present Chandra observations of 18 BALQSOs at $z\sim2$, selected from a near-infrared (2MASS) sample, where the BALQSO fraction is likely to be close to the intrinsic fraction. We measure photon indices using the stacked spectra of the optically-faint ($i-K_s\geq 2.3$ mag) and optically-bright ($i-K_s < 2.3$ mag) samples to be $\Gamma \simeq 1.5$--$2.1$. We constrain their intrinsic column density by modelling the X-ray fractional hardness ratio, finding a mean column density of $3.5\times10^{22}$ \cmsq\ assuming neutral absorption. We incorporate SDSS optical measurements (rest frame UV) to study the broadband spectral index between the X-ray and UV bands, and compare this to a large sample of normal quasars. We estimate that the optically-faint BALQSOs are X-ray weaker than the optically-bright ones, and the entire sample of BALQSOs are intrinsically X-ray weak when compared to normal AGN. Correcting for magnification of X-ray emission via gravitational lensing by the central black hole viewed at large inclination angles makes these BALQSOs even more intrinsically X-ray weak. Finally, we estimate AGN kinetic feedback efficiencies of a few percent for an X-ray wind of $0.3c$ in high-ionization BALQSOs. Combined with energy carried by low-ionization BALQSOs and UV winds, the total kinetic energy in BALQSOs can be sufficient to provide AGN kinetic feedback required to explain the co-evolution between black holes and host galaxies.

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

We present an investigation of the scaling relations between Supermassive Black Hole (SMBH) masses (Mbh), and their host galaxies' K-band bulge (Lbul) and total (Ltot) luminosities. The wide-field WIRCam imager at the Canada-France-Hawaii-Telescope (CFHT) was used to obtain the deepest and highest resolution near infrared images available for a sample of 35 galaxies with securely measured Mbh, selected irrespective of Hubble type. For each galaxy, we derive bulge and total magnitudes using a two-dimensional image decomposition code that allows us to account, if necessary, for large- and small-scale disks, cores, bars, nuclei, rings, envelopes and spiral arms. We find that the present-day Mbh-Lbul and Mbh-Ltot relations have consistent intrinsic scatter, suggesting that Mbh correlates equally well with bulge and total luminosity of the host. Our analysis provides only mild evidence of a decreased scatter if the fit is restricted to elliptical galaxies. The log-slopes of the Mbh-Lbul and Mbh-Ltot relations are 0.75+/-0.10 and 0.92+/-0.14, respectively. However, while the slope of the Mbh-Lbul relation depends on the detail of the image decomposition, the characterization of Mbh-Ltot does not. Given the difficulties and ambiguities of decomposing galaxy images into separate components, our results indicate that Ltot is more suitable as a tracer of SMBH mass than Lbul, and that the Mbh-Ltot relation should be used when studying the co-evolution of SMBHs and galaxies.

The early--type galaxy (ETG) mass--size relation has been largely studied to understand how these galaxies have assembled their mass. One key observational result of the last years is that massive galaxies increased their size by a factor of a few at fixed stellar mass from $z\sim2$. Hierarchical models favor minor mergers as a plausible driver of this size growth. Some of these models, predict a significant environmental dependence in the sense that galaxies residing in more massive halos tend to be larger than galaxies in lower mass halos, at fixed stellar mass and redshift. At present, observational results of this environmental dependence have been contradictory. In this paper we revisit this issue in the local Universe, by investigating how the sizes of massive ETGs depend on large-scale environment using an updated and accurate sample of ETGs in different environments - field, group, clusters - from the Sloan Digital Sky Survey DR7. Our analysis does not show any significant environmental dependence of the sizes of central and satellites ETGs at fixed stellar mass at $z\sim0$. The size-mass relation of early-type galaxies at $z\sim0$ seems to be universal, i.e., independent of the mass of the host halo and of the position of the galaxy in that halo (central or satellite). The result is robust to different galaxy selections based on star formation, morphology or central density. In fact, considering our observational errors and the size of the sample, any size ratio larger than $30-40%$ between massive galaxies ($log(M_*/M_\odot)>11$) living in clusters and in the field can be ruled out at $3\sigma$ level.

In addition to the large systematic differences arising from assumptions about the stellar mass-to-light ratio, the massive end of the stellar mass function is rather sensitive to how one fits the light profiles of the most luminous galaxies. We quantify this by comparing the luminosity and stellar mass functions based on SDSS cmodel magnitudes, and PyMorph single-Sersic and Sersic-Exponential fits to the surface brightness profiles of galaxies in the SDSS. The PyMorph fits return more light, so that the predicted masses are larger than when cmodel magnitudes are used. As a result, the total stellar mass density at z~0.1 is about 1.2x larger than in our previous analysis of the SDSS. The differences are most pronounced at the massive end, where the measured number density of objects having M* > 6 x 10^{11} Msun is ~5x larger. Alternatively, at number densities of 10^{-6} Mpc^{-3}, the limiting stellar mass is 2x larger. The differences with respect to fits by other authors, typically based on Petrosian-like magnitudes, are even more dramatic, although some of these differences are due to sky-subtraction problems, and are sometimes masked by large differences in the assumed $M_*/L$ (even after scaling to the same IMF). Our results impact studies of the growth and assembly of stellar mass in galaxies, and of the relation between stellar and halo mass, so we provide simple analytic fits to these new luminosity and stellar mass functions and quantify how they depend on morphology, as well as the binned counts in electronic format.

In this contributed paper I review our current knowledge of the local Black Hole (BH) scaling relations, and their impact on the determination of the local BH mass function. I particularly emphasize the remaining systematic uncertainties impinging upon a secure determination of the BH mass function and how progress can be made. I then review and discuss the evidence for a different time evolution for separate BH-galaxy scaling relations, and how these independent empirical evidences can be reconciled with the overall evolution of the structural properties of the host galaxies. I conclude discussing BH demography in the context of semi-empirical continuity accretion models, as well as more complex evolutionary models, emphasizing the general constraints we can set on them.

[abridged] The mass-size relation of early-type galaxies (ETGs) has been largely studied in the last years to probe the mass assembly of the most massive objects in the Universe. In this paper, we focus on the mass-size relation of quiescent massive ETGs (Mstar/Msol > 3*10^10) living in massive clusters (M200 ~ 10^14 Mstar) at 0.8< z <1.5, as compared to those living in the field at the same epoch. Our sample contains ~ 400 ETGs in clusters and the same number in the field. Therefore, our sample is approximately an order of magnitude larger than previous studies in the same redshift range for galaxy clusters. We find that ETGs living in clusters are between ~30-50% larger than galaxies with the same stellar mass residing in the field. We parametrize the size using the mass-normalized size, gamma=Re/Mstar^0.57. The gamma distributions in both environments peak at the same position but the distributions in clusters are more skewed towards larger sizes. Since this size difference is not observed in the local Universe, the size evolution at fixed stellar mass from z~1.5 of cluster galaxies is less steep ((1+z)-0.53pm0.04) than the evolution of field galaxies ((1+z)-0.92pm0.04). The size difference seems to be essentially driven by the galaxies residing in the clusters cores (R<0.5*R200). If part of the size evolution is due to mergers, the difference we see between cluster and field galaxies could be due to higher merger rates in clusters at higher redshift, probably during the formation phase of the clusters when velocity dispersions are lower. We cannot exclude however that the difference is driven by newly quenched galaxies which are quenched more efficiently in clusters. The implications of these results for the hierarchical growth of ETGs will be discussed in a companion paper.

A crucial challenge in astrophysics over the coming decades will be to understand the origins of supermassive black holes (SMBHs) that lie at the centres of most, if not all, galaxies. The processes responsible for the initial formation of these SMBHs and their early growth via accretion - when they are seen as Active Galactic Nuclei (AGN) - remain unknown. To address this challenge, we must identify low luminosity and obscured z>6 AGNs, which represent the bulk of early SMBH growth. Sensitive X-ray observations are a unique signpost of accretion activity, uncontaminated by star formation processes, which prevent reliable AGN identification at other wavelengths (e.g. optical, infrared). The Athena+ Wide Field Imager will enable X-ray surveys to be carried out two orders of magnitude faster than with Chandra or XMM-Newton, opening a new discovery space and identifying over 400 z>6 AGN, including obscured sources. Athena+ will also play a fundamental role to enhance the scientific return of future multiwavelength facilities that will probe the physical conditions within the host galaxies of early SMBHs, which is vital for understanding how SMBHs form, what fuels their subsequent growth, and to assess their impact on the early Universe. Follow-up of samples of z>6 galaxies with the Athena+ X-ray Integral Field Unit could also reveal the presence of highly obscured AGNs, thanks to the detection of strong iron lines. Thus, Athena+ will enable the first quantitative measurements of the extent and distribution of SMBH accretion in the early Universe.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We present a Special Issue on the interplay of galaxies and Supermassive Black Holes (SMBHs) recently published in Advances in Astronomy. This is the introductory paper containing the motivation for this Special Issue together with a brief description of the articles which are part of the manuscript and the link to the entire book (http://www.hindawi.com/journals/aa/si/610485/). We hope this Special Issue will be useful for many astronomers who want to get an update on the current status of the AGN-Galaxy coevolution topic.

[abridged] We study the dependence of the galaxy size evolution on morphology, stellar mass and large scale environment for a sample of 298 group and 384 field quiescent early-type galaxies from the COSMOS survey, selected from z~1 to the present, and with masses $log(M/M_\odot)>10.5$. The galaxy size growth depends on galaxy mass and early-type galaxy morphology, e.g., elliptical galaxies evolve differently than lenticular galaxies. At the low mass end -$10.5<Log(M/M_\odot)<11$, ellipticals do not show strong size growth from $z\sim1$ to the present (10% to 30% depending on the morphological classification). On the other end, massive ellipticals -log(M/M_\odot)>11.2$- approximately doubled their size. Interestingly, lenticular galaxies display different behavior: they appear more compact on average and they do show a size growth of \sim60% since z=1 independent of stellar mass. We compare our results with state-of-the art semi-analytic models. While major and minor mergers can account for most of the galaxy size growth, we find that with present data and the theoretical uncertainties in the modeling we cannot state clear evidence favoring either merger or mass loss via quasar and/or stellar winds as the primary mechanism driving the evolution. The galaxy mass--size relation and the size growth do not depend on environment in the halo mass range explored in this work (field to group mass $log(M_h/M_\odot)<14$), i.e., group and field galaxies follow the same trends, which is at variance with predictions from current hierarchical models that show a clear dependence of size growth on halo mass for massive ellipticals -$log(M_*/M_\odot)>11.2$.

We develop semi-empirical models of the supermassive black hole and active galactic nucleus (AGN) populations, which incorporate the black hole growth implied by the observed AGN luminosity function assuming a radiative efficiency \epsilon, and a distribution of Eddington ratios \lambda. By generalizing these continuity-equation models to allow a distribution P(\lambda|mbh,z) we are able to draw on constraints from observationally estimated P(\lambda) distributions and active galaxy fractions while accounting for the luminosity thresholds of observational samples. We consider models with a Gaussian distribution of log \lambda, and Gaussians augmented with a power-law tail to low \lambda. Within our framework, reproducing the high observed AGN fractions at low redshift requires a characteristic Eddington ratio \lambda_c that declines at late times, and matching observed Eddington ratio distributions requires P(\lambda) that broadens at low redshift. To reproduce the observed increase of AGN fraction with black hole or galaxy mass, we also require a \lambda_c that decreases with increasing black hole mass, reducing the AGN luminosity associated with the most massive black holes. Finally, achieving a good match to the high mass end of the local black hole mass function requires an increased radiative efficiency at high black hole mass. We discuss the potential impact of black hole mergers or a \lambda-dependent bolometric correction, and we compute evolutionary predictions for black hole and galaxy specific accretion rates. Despite the flexibility of our framework, no one model provides a good fit to all the data we consider.

We present basic predictions of an updated version of the Munich semi-analytic hierarchical galaxy formation model that grows bulges via mergers and disk instabilities. Overall, we find that while spheroids below Ms ~ 10^11 Msun grow their sizes via a mixture of disk instability and mergers, galaxies above it mainly evolve via mergers. Including gas dissipation in major mergers, efficiently shrinks galaxies, especially those with final mass Ms < 10^11 Msun that are the most gas-rich, improving the match with different observables. We find that the predicted scatter in sizes at fixed stellar mass is still larger than the observed one by up to <40%. Spheroids are, on average, more compact at higher redshifts at fixed stellar mass, and at fixed redshift and stellar mass larger galaxies tend to be more starforming. More specifically, while for bulge-dominated galaxies the model envisages a nearly mass-independent decrease in sizes, the predicted size evolution for intermediate-mass galaxies is more complex. The z=2 progenitors of massive galaxies with mass around Ms and B/T>0.7 at z=0, are found to be mostly disc-dominated galaxies with a median B/T ~ 0.3, with only ~20% remaining bulge-dominated. The model also predicts that central spheroids living in more massive haloes tend to have larger sizes at fixed stellar mass. Including host halo mass dependence in computing velocity dispersions, allows the model to properly reproduce the correlations with stellar mass. We also discuss the fundamental plane, the correlations with galaxy age, the structural properties of pseudobulges, and the correlations with central black holes.

We present the first direct measurement of the mean Halo Occupation Distribution (HOD) of X-ray selected AGN in the COSMOS field at z < 1, based on the association of 41 XMM and 17 C-COSMOS AGN with member galaxies of 189 X-ray detected galaxy groups from XMM and Chandra data. We model the mean AGN occupation in the halo mass range logM_200[Msun] = 13-14.5 with a rolling-off power-law with the best fit index alpha = 0.06(-0.22;0.36) and normalization parameter f_a = 0.05(0.04;0.06). We find the mean HOD of AGN among central galaxies to be modelled by a softened step function at logMh > logMmin = 12.75 (12.10,12.95) Msun while for the satellite AGN HOD we find a preference for an increasing AGN fraction with Mh suggesting that the average number of AGN in satellite galaxies grows slower (alpha_s < 0.6) than the linear proportion (alpha_s = 1) observed for the satellite HOD of samples of galaxies. We present an estimate of the projected auto correlation function (ACF) of galaxy groups over the range of r_p = 0.1-40 Mpc/h at <z> = 0.5. We use the large-scale clustering signal to verify the agreement between the group bias estimated by using the observed galaxy groups ACF and the value derived from the group mass estimates. We perform a measurement of the projected AGN-galaxy group cross-correlation function, excluding from the analysis AGN that are within galaxy groups and we model the 2-halo term of the clustering signal with the mean AGN HOD based on our results.

Using the wide multi-band photometry available in the COSMOS field we explore the host galaxy properties of a large sample of Active Galactic Nuclei (AGN) obtained by combining X-ray and optical spectroscopic selections. Based on a careful study of their Spectral Energy Distribution (SED), which has been parametrized using a 2-component (AGN+galaxy) model fit, we derived dust-corrected rest-frame magnitudes, colors, stellar masses and star formation rates (SFRs). We find that AGN hosts span a large range of stellar masses and SFRs. No color-bimodality is seen at any redshift in the AGN hosts, which are found to be mainly massive, red galaxies. Once accounting for the color-mass degeneracy in well defined mass-matched samples, we find a residual marginal enhancement of AGN incidence in redder galaxies with lower specific star formation rates, and we argue that this result might emerge because of our ability to properly account for AGN light contamination and dust extinction. Interestingly, we find that the probability for a galaxy to host a black hole growing at any given "specific accretion rate" (i.e. the ratio of X-ray luminosity to the host stellar mass) is almost independent of the host galaxy mass, while it decreases as a power-law with Lx/M. By analyzing the normalization of such probability distribution, we show how the incidence of AGN increases with redshift as rapidly as (1+z)^4, in close resemblance with the overall evolution of the specific star formation rate of the entire galaxy population. Although AGN activity and star formation in galaxies do appear to have a common triggering mechanism, at least in a statistical sense, within the COSMOS sample we do not find strong evidence of any 'smoking gun' signaling powerful AGN influence on the star-forming properties of their hosts galaxies.

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are presented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% and 52% respectively), while accurate photometric redshifts are available for the rest of the sample. The study of such a large sample of X-ray selected AGN with a high quality multi-wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 micron and 160 micron) to the optical-UV allows us to obtain accurate estimates of bolometric luminosities, bolometric corrections and Eddington ratios. The kbol-Lbol relations derived in the present work are calibrated for the first time against a sizable AGN sample, and rely on observed redshifts, X-ray luminosities and column density distributions. We find that kbol is significantly lower at high Lbol with respect to previous estimates by Marconi et al. (2004) and Hopkins et al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while black hole masses for Type-2 AGN are computed for 481 objects using the black hole mass-stellar mass relation and the morphological information. We confirm a trend between kbol and lambda_Edd, with lower hard X-ray bolometric corrections at lower Eddington ratios for both Type-1 and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given Lbol.

We confirm the detection of 3 groups in the Lynx supercluster, at z~1.3, and give their redshifts and masses. We study the properties of the group galaxies as compared to the central clusters, RXJ0849+4452 and RXJ0848+4453, selecting 89 galaxies in the clusters and 74 galaxies in the groups. We morphologically classify galaxies by visual inspection, noting that our early-type galaxy (ETG) sample would have been contaminated at the 30% -40% level by simple automated classification methods (e.g. based on Sersic index). In luminosity selected samples, both clusters and groups show high fractions of Sa galaxies. The ETG fractions never rise above ~50% in the clusters, which is low compared to the fractions observed in clusters at z~1. However, ETG plus Sa fractions are similar to those observed for ETGs in clusters at z~1. Bulge-dominated galaxies visually classified as Sas might also be ETGs with tidal features or merger remnants. They are mainly red and passive, and span a large range in luminosity. Their star formation seems to have been quenched before experiencing a morphological transformation. Because their fraction is smaller at lower redshifts, they might be the spiral population that evolves into ETGs. For mass-selected samples, the ETG fraction show no significant evolution with respect to local clusters, suggesting that morphological transformations occur at lower masses and densities. The ETG mass-size relation shows evolution towards smaller sizes at higher redshift in both clusters and groups, while the late-type mass-size relation matches that observed locally. The group ETG red sequence shows lower zero points and larger scatters than in clusters, both expected to be an indication of a younger galaxy population. The estimated age difference is small when compared to the difference in age at different galaxy masses.

There is mounting evidence that a significant fraction of Black Holes (BHs) today live in late-type galaxies, including bulge-less galaxies and those hosting pseudobulges, and are significantly undermassive with respect to the scaling relations followed by their counterpart BHs in classical bulges of similar stellar (or even bulge) mass. Here we discuss the predictions of two state-of-the-art hierarchical galaxy formation models in which BHs grow via mergers and, in one, also via disk instability. Our aim is to understand if the wealth of new data on local BH demography is consistent with standard models. We follow the merger trees of representative subsamples of BHs and compute the fractional contributions of different processes to the final BH mass. We show that the model in which BHs always closely follow the growth of their host bulges, also during late disk instabilities (i.e., bars), produces too narrow a distribution of BHs at fixed stellar mass to account for the numerous low-mass BHs now detected in later-type galaxies. Models with a looser connection between BH growth and bar instability instead predict the existence of a larger number of undermassive BHs, in better agreement with the observations. The scatter in the updated local BH-bulge mass relation (with no restriction on galaxy type) appears to be quite large when including later-type systems, but it can still be managed to be reproduced within current hierarchical models. However, the fuelling of BHs during the late bar-instability mode needs to be better quantified/improved to properly fit the data. We conclude discussing how the possibly large number of BHs in later type galaxies demands for an in-depth revision of the local BH mass function and its modelling.

We present the Kormendy and mass-size relations for early-type galaxies (ETGs) as a function of environment at z~1.3. Our sample includes 76 visually classified ETGs with masses 10^10 < M/Msun < 10^11.5, selected in the Lynx supercluster and in the GOODS/CDF-S field, 31 ETGs in clusters, 18 in groups and 27 in the field, all with multi-wavelength photometry and HST/ACS observations. The Kormendy relation, in place at z~1.3, does not depend on the environment. The mass-size relation reveals that ETGs overall appear to be more compact in denser environments: cluster ETGs have sizes on average around 30-50% smaller than those of the local universe, and a distribution with a smaller scatter, whereas field ETGs show a mass-size relation with a similar distribution than the local one. Our results imply that (1) the mass-size relation in the field did not evolve overall from z ~ 1.3 to present; this is interesting and in contrast to the trend found at higher masses from previous works; (2) in denser environments, either ETGs have increased their size by 30-50%, on average, and spread their distributions, or more ETGs have been formed within the dense environment from not ETG progenitors or larger galaxies have been accreted to a pristine compact population to reproduce the mass-size relation observed in the local Universe. Our results are driven by galaxies with masses M<2*10^11Msun and those with masses M~10^11Msun follow the same trends that the entire sample. Following Valentinuzzi et al. definition of superdense ETGs, around 35-45% of our cluster sample is made of superdense ETGs.

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.

We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 Msec Chandra observation and adopting a new detection algorithm, based on a targeted search at the position of known high-z galaxies. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 42.75<log L(2-10 keV)<44.5. We join this data with the luminous AGN luminosity functions from optical surveys and find that the evolution of the high-z, wide luminosity range luminosity function can be best modeled by pure luminosity evolution with L* decreasing from 6x10^44 ergs/s at z=3 to L*=2x10^44 ergs/s at z=6. We compare the high-z luminosity function with the prediction of theoretical models using galaxy interactions as AGN triggering mechanism. We find that these models are broadly able to reproduce the high-z AGN luminosity functions. A better agreement is found assuming a minimum Dark Matter halo mass for black hole formation and growth. We compare our AGN luminosity functions with galaxy mass functions to derive high-z AGN duty cycle using observed Eddington ratio distributions to derive black hole masses. We find that the duty cycle increases with galaxy stellar mass and redshift by a factor 10-30 from z=0.25 to z=4-5. We also report on the detection of a large fraction of highly obscured, Compton thick AGN at z>3 (18+17-10%). Their optical counterparts are not strongly reddened and we thus conclude that the size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly star-forming galaxy at z=3.47. If confirmed, this would be one of the farthest objects in which stellar sources are detected in X-rays.

We present the high redshift (3<z<5.3) 0.5-2 keV number counts and the 2-10 keV (rest frame) space density of X-ray selected AGNs detected in the Chandra COSMOS survey. The sample comprises 81 X-ray detected sources with available spectroscopic (31) and photometric (50) redshifts plus 20 sources with a formal z_phot<3 but with a broad photometric redshift probability distribution, such that z_phot+1sigma>3. 81 sources are selected in the 0.5-2 keV band, 14 are selected in the 2-10 keV and 6 in the 0.5-10 keV bands. We sample the high luminosity (log L_(2-10 keV)>44.15 erg/s) space density up to z~5 and a fainter luminosity range (43.5<log L(2-10 keV)<44.15 erg/s) than previous studies, up to z=3.5. We weighted the contribution to the number counts and the space density of the sources with photometric redshift by using their probability of being at z>3. We find that the space density of high-luminosity AGNs declines exponentially at all the redshifts, confirming the trend observed for optically selected quasars. At lower luminosity, the measured space density is not conclusive, and a larger sample of faint sources is needed. Comparisons with optical luminosity functions and BH formation models are presented together with prospects for future surveys.

We extend a chemical evolution model relating galaxy stellar mass and gas-phase oxygen abundance (the mass-metallicity relation) to explicitly consider the mass-dependence of galaxy gas fractions and outflows. Using empirically derived scalings of galaxy mass with halo virial velocity in conjunction with the most recent observations of z~0 total galaxy cold gas fractions and the mass-metallicity relation, we place stringent global constraints on the magnitude and scaling of the efficiency with which star forming galaxies expel metals. We demonstrate that under the assumptions that metal accretion is negligible and the stellar initial mass function does not vary, efficient outflows are required to reproduce the mass-metallicity relation; without winds, gas-to-stellar mass ratios >~ 0.3 dex higher than observed are needed. Moreover, z=0 gas fractions are low enough that while they have some effect on the magnitude of outflows required, the slope of the gas fraction--stellar mass relation does not strongly affect our conclusions on how the wind efficiencies must scale with galaxy mass. Despite systematic uncertainties in the normalization and slope of the mass-metallicity relation, we show that the metal expulsion efficiency zetaw=(Zw/Zg)etaw (where Zw is the wind metallicitiy and Zg is the interstellar medium metallicity) must be both high and scale steeply with mass. Specifically, we show that zetaw >> 1 and zetaw proportional to vvir^-3 or steeper. In contrast, momentum- or energy-driven outflow models suggest that etaw should scale as vvir^-1 or vvir^-2, respectively, implying that the Zw-Mstar relation should be shallower than the Zg-Mstar relation. [abridged]

We develop a new empirical methodology to study the relation between the stellar mass of galaxies and the mass of their host subhaloes. Our approach is similar to abundance matching, and is based on assigning a stellar mass to each subhalo within a large cosmological N-body simulation. However, we significantly extend the method to include a different treatment for central and satellite galaxies, allowing the stellar mass of satellite galaxies to depend on both the host subhalo mass, and on the halo mass. Unlike in previous studies, our models are constrained by two observations: the stellar mass function of galaxies, and their auto-correlation functions (CFs). We present results for ~10^6 different successful models, showing that the uncertainty in the host subhalo mass reaches a factor of ~10 for a given stellar mass. Our results thus indicate that the host subhalo mass of central and satellite galaxies is poorly constrained by using information coming solely from the abundance and CFs of galaxies. In addition, we demonstrate that the fraction of stellar mass relative to the universal baryon fraction can vary between different models, and can reach ~0.6 for a specific population of low mass galaxies. We conclude that additional observational constraints are needed, in order to better constrain the mass relation between haloes and galaxies. These might include weak lensing, satellite kinematics, or measures of clustering other than the CFs. Moreover, CFs based on larger surveys with a better estimate of the systematic uncertainties are needed.

We present Suzaku observations of three iron low-ionization broad absorption line quasars (FeLoBALs). We detect J1723+5553 (3\sigma) in the observed 2--10 keV band, and constrain its intrinsic nh column density to nh > 6\times10^{23} \cmsq by modeling its X-ray hardness ratio. We study the broadband spectral index, aox, between the X-ray and UV bands by combining the X-ray measurements and the UV flux extrapolated from 2MASS magnitudes, assuming a range of intrinsic column densities, and then comparing the aox values for the three FeLoBALs with those from a large sample of normal quasars. We find that the FeLoBALs are consistent with the spectral energy distribution (SED) of normal quasars if the intrinsic nh column densities are nh > 7\times10^{23} \cmsq for J0943+5417, nh > 2\times10^{24} \cmsq for J1352+4293, and 6\times10^{23} < \nh < 3\times10^{24} \cmsq for J1723+5553. At these large intrinsic column densities, the optical depth from Thomson scattering can reach \sim6, which will significantly modulate the UV flux. Our results suggest that the X-ray absorbing material is located at a different place from the UV absorbing wind, likely between the X-ray and UV emitting regions. We find a significant kinetic feedback efficiency for FeLoBALs, indicating that the outflows are an important feedback mechanism in quasars.

We present a study of the redshift evolution of the projected correlation function of 593 X-ray selected AGN with I_AB<23 and spectroscopic redshifts z<4, extracted from the 0.5-2 keV X-ray mosaic of the 2.13 deg^2 XMM-COSMOS survey. We introduce a method to estimate the average bias of the AGN sample and the mass of AGN hosting halos, solving the sample variance using the halo model and taking into account the growth of the structure over time. We find evidence of a redshift evolution of the bias factor for the total population of XMM-COSMOS AGN from b(z=0.92)=2.30 +/- 0.11 to b(z=1.94)=4.37 +/- 0.27 with an average mass of the hosting DM halos logM [h^-1 M_sun] ~ 13.12 +/- 0.12 that remains constant at all z < 2. Splitting our sample into broad optical lines AGN (BL), AGN without broad optical lines (NL) and X-ray unobscured and obscured AGN, we observe an increase of the bias with redshift in the range z=0.7-2.25 and z=0.6-1.5 which corresponds to a constant halo mass logM [h^-1 M_sun] ~ 13.28 +/- 0.07 and logM [h^-1 M_sun] ~ 13.00 +/- 0.06 for BL /X-ray unobscured AGN and NL/X-ray obscured AGN, respectively. The theoretical models which assume a quasar phase triggered by major mergers can not reproduce the high bias factors and DM halo masses found for X-ray selected BL AGN with L_BOL ~ 2e45 erg s^-1. Our work extends up to z ~ 2.2 the z <= 1 statement that, for moderate luminosity X-ray selected BL AGN, the contribution from major mergers is outnumbered by other processes, possibly secular such as tidal disruptions or disk instabilities.

Low-ionization (MgII, FeII, FeIII) broad absorption line quasars (LoBALs) probe a relatively obscured quasar population, and could be at an early evolutionary stage for quasars. We study the intrinsic fractions of LoBALs using the SDSS, 2MASS, and FIRST surveys. We find that the LoBAL fractions of the near infra-red (NIR) and radio samples are approximately 5--7 times higher than those measured in the optical sample. This suggests that the fractions measured in the NIR and radio bands are closer to the intrinsic fractions of the populations, and that the optical fractions are significantly biased due to obscuration effects, similar to high-ionization broad absorption line quasars (HiBALs). We also find that the LoBAL fractions decrease with increasing radio luminosities, again, similar to HiBALs. In addition, we find tentative evidence for high fractions of LoBALs at high NIR luminosities, especially for FeLoBALs with a fraction of ~18 per cent at M_K_s < -31 mag. This population of NIR luminous LoBALs may be at an early evolutionary stage of quasar evolution. We use a two-component model of LoBALs including a pure geometric component and a luminosity dependent component at high NIR luminosities, and obtain better fits than those from a pure geometric model. Therefore, the LoBAL population can be modelled as a hybrid of both the geometric and evolutionary models, where the geometric component constitutes 3.4\pm0.3, 5.8\pm0.4, and 1.5\pm0.3 per cent of the quasar population for BI-LoBALs, AI-LoBALs, and FeLoBALs, respectively. Considering a population of obscured quasars that do not enter the SDSS survey, which could have a much higher LoBAL fraction, we expect that intrinsic fraction of LoBALs could be even higher.

For early-type galaxies, the correlations between stellar mass and size, velocity dispersion, surface brightness, color, axis ratio and color-gradient all indicate that two mass scales, M* = 3 x 10^10 Msun and M* = 2 x 10^11 Msun, are special. The smaller scale could mark the transition between wet and dry mergers, or it could be related to the interplay between SN and AGN feedback, although quantitative measures of this transition may be affected by morphological contamination. At the more massive scale, mean axis ratios and color gradients are maximal, and above it, the colors are redder, the sizes larger and the velocity dispersions smaller than expected based on the scaling at lower M*. In contrast, the color-sigma relation, and indeed, most scaling relations with sigma, are not curved: they are well-described by a single power law, or in some cases, are almost completely flat. When major dry mergers change masses, sizes, axis ratios and color gradients, they are expected to change the colors or velocity dispersions much less. Therefore, the fact that scaling relations at sigma > 150 km/s show no features, whereas the size-M*, b/a-M*, color-M* and color gradient-M* relations do, suggests that M* = 2 x 10^11 Msun is the scale above which major dry mergers dominate the assembly histories of early-type galaxies.

The color-magnitude relation of early-type galaxies differs slightly but significantly from a pure power-law, curving downwards at low and upwards at large luminosities (Mr>-20.5 and Mr<-22.5). This remains true of the color-size relation, and is even more apparent with stellar mass (M* < 3x10^10 Msun and M* > 2x10^11 Msun). The upwards curvature at the massive end does not appear to be due to stellar population effects. In contrast, the color-sigma relation is well-described by a single power law. Since major dry mergers change neither the colors nor sigma, but they do change masses and sizes, the clear features observed in the scaling relations with M*, but not with sigma > 150 km/s, suggest that M* > 2x10^11 Msun is the scale above which major dry mergers dominate the assembly history. We discuss three models of the merger histories since z ~ 1 which are compatible with our measurements. In all three models, dry mergers are responsible for the flattening of the color-M* relation at M* > 3x10^10 Msun - wet mergers only matter at smaller masses. At M* > 2 x 10^11 Msun, the merger histories in one model are dominated by major rather than minor dry mergers, as suggested by the axis ratio and color gradient trends. In another, although both major and minor mergers occur at the high mass end, the minor mergers contribute primarily to the formation of the ICL, rather than to the mass growth of the central massive galaxy. A final model assumes that the reddest objects were assembled by a mix of major and minor dry mergers.

We investigate the characteristic radiative efficiency \epsilon, Eddington ratio \lambda, and duty cycle P_0 of high-redshift active galactic nuclei (AGN), drawing on measurements of the AGN luminosity function at z=3-6 and, especially, on recent measurements of quasar clustering at z=3-4.5 from the Sloan Digital Sky Survey. The free parameters of our models are \epsilon, \lambda, and the normalization, scatter, and redshift evolution of the relation between black hole mass \mbh and halo virial velocity V_vir. We compute the luminosity function from the implied growth of the black hole mass function and the quasar correlation length from the bias of the host halos. We test our adopted formulae for the halo mass function and halo bias against measurements from the large N-body simulation developed by the MICE collaboration. The strong clustering of AGNs observed at z=3 and, especially, at z=4 implies that massive black holes reside in rare, massive dark matter halos. Reproducing the observed luminosity function then requires high efficiency \epsilon and/or low Eddington ratio \lambda, with a lower limit (based on 2\sigma agreement with the measured z=4 correlation length) \epsilon> 0.7\lambda/(1+0.7\lambda), implying \epsilon > 0.17 for \lambda > 0.25. Successful models predict high duty cycles, P_0~0.2, 0.5, and 0.9 at z=3.1, 4.5 and 6, respectively, and they require that the fraction of halo baryons locked in the central black hole is much larger than the locally observed value. The rapid drop in the abundance of the massive and rare host halos at z>7 implies a proportionally rapid decline in the number density of luminous quasars, much stronger than simple extrapolations of the z=3-6 luminosity function would predict. (abridged)

We use Shen et al.'s (2009) measurements of luminosity-dependent clustering in the SDSS Data Release 5 Quasar Catalog, at redshifts 0.4 < z < 2.5, to constrain the relation between quasar luminosity and host halo mass and to infer the duty cycle f_opt, the fraction of black holes that shine as optically luminous quasars at a given time. We assume a monotonic mean relation between quasar luminosity and host halo mass, with log-normal scatter \Sigma. For specified f_opt and \Sigma, matching the observed quasar space density determines the normalization of the luminosity-halo mass relation, from which we predict the clustering bias. The data show no change of bias between the faint and bright halves of the quasar sample but a modest increase in bias for the brightest 10%. At the mean redshift z=1.45 of the sample, the data can be well described either by models with small intrinsic scatter (\Sigma=0.1 dex) and a duty cycle f_opt=6*10^(-4) or by models with much larger duty cycles and larger values of the scatter. "Continuity equation" models of the black hole mass population imply f_opt > 2*10^(-3) in this range of masses and redshifts, and the combination of this constraint with the clustering measurements implies scatter \Sigma > 0.4 dex. These findings contrast with those inferred from the much stronger clustering of high-luminosity quasars at z~4, which require minimal scatter between luminosity and halo mass and duty cycles close to one.

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

We provide fits to the distribution of galaxy luminosity, size, velocity dispersion and stellar mass as a function of concentration index C_r and morphological type in the SDSS. We also quantify how estimates of the fraction of `early' or `late' type galaxies depend on whether the samples were cut in color, concentration or light profile shape, and compare with similar estimates based on morphology. Our fits show that Es account for about 20% of the r-band luminosity density, rho_Lr, and 25% of the stellar mass density, rho_*; including S0s and Sas increases these numbers to 33% and 40%, and 50% and 60%, respectively. Summed over all galaxy types, we find rho_* ~ 3 * 10^8 M_Sun Mpc^{-3} at z ~ 0. This is in good agreement with expectations based on integrating the star formation history. However, compared to most previous work, we find an excess of objects at large masses, up to a factor of ~ 10 at M_* ~ 5*10^{11} M_Sun. The stellar mass density further increases at large masses if we assume different IMFs for Es and spiral galaxies, as suggested by some recent chemical evolution models, and results in a better agreement with the dynamical mass function. We also show that the trend for ellipticity to decrease with luminosity is primarily because the E/S0 ratio increases at large L. However, the most massive galaxies, M_* > 5 * 10^{11} M_Sun, are less concentrated and not as round as expected if one extrapolates from lower L, and they are not well-fit by pure deVaucouleur laws. This suggests formation histories with recent radial mergers. Finally, we show that the age-size relation is flat for Es of fixed dynamical mass, but, at fixed M_dyn, S0s and Sas with large sizes tend to be younger. Explaining this difference between E and S0 formation is a new challenge for models of early-type galaxy formation.

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

We cross-correlate the SDSS DR3 quasar sample with FIRST and the Vestergaard et al. black hole (BH) mass sample to compare the mean accretion histories of optical and radio quasars. We find significant statistical evidence that radio quasars have a higher mean Eddington ratio Lambda at z > 2 with respect to optical quasars, while the situation is clearly reverse at z < 1. At z > 2 radio quasars happen to be less massive than optical quasars; however, as redshift decreases radio quasars appear in increasingly more massive BHs with respect to optical quasars. These two trends imply that radio sources are not a mere random subsample of optical quasars. No clear correlation between radio activity and BH mass and/or accretion rate is evident from our data, pointing to other BH properties, possibly the spin, as the driver of radio activity. We have checked that our main results do not depend on any evident bias. We perform detailed modelling of reasonable accretion histories for optical and radio quasars, finding that radio quasars grow by a factor of a few, at the most, since z ~ 4. The comparison between the predicted mass function of active radio quasars and the observed optical luminosity function of radio quasars, implies a significantly lower probability for lower mass BHs to be radio loud at all epochs, in agreement with what is observed in the local universe.

We use the large catalogues of haloes available for the Millennium Simulation to test whether recently merged haloes exhibit stronger large-scale clustering than other haloes of the same mass. This effect could help to understand the very strong clustering of quasars at high redshift. However, we find no statistically significant excess bias for recently merged haloes over the redshift range 2 < z < 5, with the most massive haloes showing an excess of at most ~5%. We also consider galaxies extracted from a semianalytic model built on the Millennium Simulation. At fixed stellar mass, we find an excess bias of ~ 20-30% for recently merged objects, decreasing with increasing stellar mass. The fact that recently-merged galaxies are found in systematically more massive haloes than other galaxies of the same stellar mass accounts for about half of this signal, and perhaps more for high-mass galaxies. The weak merger bias of massive systems suggests that objects of merger-driven nature, such as quasars, do not cluster significantly differently than other objects of the same characteristic mass. We discuss the implications of these results for the interpretation of clustering data with respect to quasar duty cycles, visibility times, and evolution in the black hole-host mass relation.

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

Supermassive black holes (BHs) appear to be ubiquitous at the center of all galaxies which have been observed at high enough sensitivities and resolution with the Hubble Space Telescope. Their masses are found to be tightly linked with the masses and velocity dispersions of their host galaxies. On the other hand, BHs are widely held to constitute the central engines of quasars and active galactic nuclei (AGN) in general. It is however still unclear how BHs have grown, and whether they have co-evolved with their hosts. In this Review I discuss how, in ways independent of specific models, constraints on the growth history of BHs and their host galaxies have been set by matching the statistics of local BHs to the emissivity, number density, and clustering properties of AGNs at different cosmological epochs. I also present some new results obtained through a novel numerical code which evolves the BH mass function and clustering adopting broad distributions of Eddington ratios. I finally review BH evolution in a wider cosmological context, connecting BH growth to galaxy evolution.

We use a sample of about 48,000 SDSS early-type galaxies to show that older galaxies have smaller half-light radii re and larger velocity dispersions sigma than younger ones of the same stellar mass Mstar. We use the age-corrected luminosity Lcorr as a proxy for Mstar to minimize biases: below Lcorr~1e11 Lsun, galaxies with age ~11 Gyrs have re smaller by 40% and sigma larger by 25%, compared to galaxies that are ~4 Gyrs younger. The sizes and velocity dispersions of more luminous galaxies vary by less than 15%, whatever their age, a challenge for current galaxy formation models. A closer check reveals that the lowering in the dispersion is caused by older galaxies that show a significant departure from the re--Lcorr and sigma--Lcorr relations at high Lcorr. Such features might find an explanation in models where more massive galaxies undergo more minor mergers than less massive galaxies at late times, thus causing a break in the homology. In terms of the Fundamental Plane of early-type galaxies, the data indicate that all galaxies show a significant and similar increase in the dynamical-to-stellar mass ratio with increasing mass, independent of their age. However, older galaxies have smaller Mdyn/Mstar ratios than objects which formed more recently. These findings may suggest that lower mass galaxies and, at fixed stellar mass, higher redshift galaxies, formed from gas-richer progenitors, thus underwent more dissipation and contraction.

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2\sigma$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

We construct evolutionary models of the populations of AGN and supermassive black holes, in which the black hole mass function grows at the rate implied by the observed luminosity function, given assumptions about the radiative efficiency and the Eddington ratio. We draw on a variety of recent X-ray and optical measurements to estimate the bolometric AGN luminosity function and compare to X-ray background data and the independent estimate of Hopkins et al. (2007) to assess remaining systematic uncertainties. The integrated AGN emissivity closely tracks the cosmic star formation history, suggesting that star formation and black hole growth are closely linked at all redshifts. Observational uncertainties in the local black hole mass function remain substantial, with estimates of the integrated black hole mass density \rho_BH spanning the range 3-5.5x10^5 Msun/Mpc^3. We find good agreement with estimates of the local mass function for a reference model where all active black holes have efficiency \eps=0.065 and L_bol/L_Edd~0.4. In this model, the duty cycle of 10^9 Msun black holes declines from 0.07 at z=3 to 0.004 at z=1 and 0.0001 at z=0. The decline is shallower for less massive black holes, a signature of "downsizing" evolution in which more massive black holes build their mass earlier. The predicted duty cycles and AGN clustering bias in this model are in reasonable accord with observational estimates. If the typical Eddington ratio declines at z<2, then the "downsizing" of black hole growth is less pronounced. Matching the integrated AGN emissivity to the local black hole mass density implies \eps=0.075 (\rho_BH/4.5x10^5 Msun/Mpc^3)^{-1} for our standard luminosity function estimate (25% higher for Hopkins et al.'s), lower than the values \eps=0.16-0.20 predicted by MHD simulations of disk accretion.

We find that the fraction of classical Broad Absorption Line quasars (BALQSOs) among the FIRST radio sources in the Sloan Data Release 3, is 20.5^{+7.3}_{-5.9}% at the faintest radio powers detected (L_{\rm 1.4 GHz}~10^{32} erg/s), and rapidly drops to <8% at L_{\rm 1.4 GHz}~3*10^{33} erg/s. Similarly, adopting the broader Absorption Index (AI) definition of Trump et al. (2006) we find the fraction of radio BALQSOs to be 44^{+8.1}_{-7.8}% reducing to 23.1^{+7.3}_{-6.1}% at high luminosities. While the high fraction at low radio power is consistent with the recent near-IR estimates by Dai et al. (2008), the lower fraction at high radio powers is intriguing and confirms previous claims based on smaller samples. The trend is independent of the redshift range, the optical and radio flux selection limits, or the exact definition of a radio match. We also find that at fixed optical magnitude, the highest bins of radio luminosity are preferentially populated by non-BALQSOs, consistent with the overall trend. We do find, however, that those quasars identified as AI-BALQSOs but \emph{not} under the classical definition, do not show a significant drop in their fraction as a function of radio power, further supporting independent claims for which these sources, characterized by lower equivalent width, may represent an independent class with respect to the classical BALQSOs. We find the balnicity index, a measure of the absorption trough in BALQSOs, and the mean maximum wind velocity to be roughly constant at all radio powers. We discuss several plausible physical models which may explain the observed fast drop in the fraction of the classical BALQSOs with increasing radio power, \emph{although no one is entirely satisfactory}. (abridged).

We utilize the local velocity dispersion function (VDF) of spheroids, together with their inferred age--distributions, to predict the VDF at higher redshifts (0<z<6), under the assumption that (i) most of the stars in each nearby spheroid formed in a single episode, and (ii) the velocity dispersion sigma remained nearly constant afterward. We assume further that a supermassive black hole (BH) forms concurrently with the stars, and within ~1 Gyr of the formation of the potential well of the spheroid, and that the relation between the mass of the BH and host velocity dispersion maintains the form M_BH ~ sigma^{beta} with beta~4, but with the normalization allowed to evolve with redshift as ~(1+z)^{alpha}. We compute the BH mass function associated with the VDF at each redshift, and compare the accumulated total BH mass density with that inferred from the integrated quasar luminosity function (LF; the so--called Soltan argument). This comparison is insensitive to the assumed duty cycle or Eddington ratio of quasar activity, and we find that the match between the two BH mass densities favors a relatively mild redshift evolution, with alpha ~ 0.26, with a positive evolution as strong as alpha>1.3 excluded at the 99% confidence level. A direct match between the characteristic BH mass in the VDF--based and quasar LF--based BH mass functions also yields a mean Eddington ratio of lambda ~ 0.5-1 that is roughly constant within 0<z<3. A strong positive evolution in the M_BH-sigma relation is still allowed by the data if galaxies increase, on average, their velocity dispersions since the moment of formation, due to dissipative processes. If we assume that the mean velocity dispersion of the host galaxies evolves as sigma(z)=sigma(0)*(1+z)^{-gamma}, we find a lower limit of gamma>0.23 for alpha>1.5. abridged

We constrain the mean kinetic efficiency of radio-loud active galactic nuclei by using an optically selected sample for which both the optical and the radio luminosity functions (LFs) have been determined; the former traces the bolometric luminosity L, while the latter traces the kinetic power L_k, empirically correlated to the radio emission. Thus in terms of the ratio g_k=L_k/L, we can convert the optical LF of the sample into a radio one. This computed LF is shown to match the directly observed LF for the same sample if g_k=0.10^{+0.05}_{-0.01} holds, with a scatter \sigma=0.38^{+0.04}_{-0.09} dex; with these values we also match a number of independent correlations between L_k, L and radio emission, that we derive through Monte Carlo simulations. We proceed to translate the value of g_k into a constraint on the kinetic efficiency for the production of radio jets or winds, namely, \epsilon_k=L_k/(Mdot*c^2)~0.01 in terms of the rate Mdot of mass accretion onto the central black hole. Then, on assuming that on average the radio sources share the same kinetic efficiency, we compute a solid lower limit of about 25% on the contribution of radio sources to the local black hole mass density.

The intrinsic fraction of broad absorption line quasars (BALQSOs) is important in constraining geometric and evolutionary models of quasars. We present the fraction of BALQSOs in 2MASS detected quasars within the SDSS DR3 sample in the redshift range of 1.7 < z < 4.38. The fraction of BALQSOs is 40.4^{+3.4}_{-3.3}% in the 2MASS 99% database K_s band completeness sample, and 38.5^{+1.7}_{-1.7}% in the larger 2MASS sample extending below the completeness limit. These fractions are significantly higher than the 26% reported in the optical bands for the same parent sample. We also present the fraction of BALQSOs as functions of apparent magnitudes, absolute magnitudes, and redshift in the 2MASS and SDSS bands. The 2MASS fractions are consistently higher than the SDSS fractions in every comparison, and the BALQSO fractions steadily increase with wavelength from the SDSS u to the 2MASS K_s bands. Furthermore, the i - K_s color distributions of BALQSOs and non-BALQSOs indicate that BALQSOs are redder than non-BALQSOs, with a K-S test probability of 2e-12. These results are consistent with the spectral difference between BALQSOs and non-BALQSOs including both the absorption troughs and dust extinction in BALQSOs, which leads to significant selection biases against BALQSOs in the optical bands. Using a simple simulation incorporating the luminosity function of quasars and the amount of obscuration for BALQSOs, we simultaneously fit the BALQSO fractions in the SDSS and 2MASS bands. We obtain a true BALQSO fraction of 43\pm2% for luminous quasars (M_{K_s} \lesssim -30.1 mag).

Using the results of recent optical surveys we conclude that the {\it non}-detection of quasars down to faint magnitudes implies a significant flattening of the high redshift (z~6) optical active galactic nuclei (AGN) luminosity function for M_{1450}>-24.7. We find that all the data are consistent with a faint-end slope for the optical AGN luminosity function of \beta=-2.2 and \beta=-2.8, at the 90% and 99% confidence level respectively, flatter than the bright-end slope of \beta'~ -3.2. We also show that X-ray deep surveys have probed even fainter magnitudes than the optical ones yielding more significant constraints on the shallow faint-end slope of the optical luminosity function. The inclusion of Type II AGN candidates, detected in the Chandra deep fields, hints towards an higher normalization for the total AGN luminosity function, if these sources lie at 5<z <6.5. We then discuss simple theoretical models of AGN formation and evolution in the context of cold dark matter cosmology. The comparison with the total AGN luminosity function favors a redshift-dependent relation between black hole and dark matter halo masses of the type M_BH~M_ halo^\alpha, with 1.3< \alpha< 1.7, compatible with independent studies from statistical analysis and rotation curve measurements. Finally we compute the quasar contribution to reionization to be <9% at z~6, up to ~30% when integrated within 5.5< z<6.5, significantly smaller than that from galaxies.

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of black holes and host galaxies predicts quasar luminosity functions at redshifts 1.5<z<6 and local properties in nice agreement with observations. In our model the quasar activity marks and originates the transition between an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while the latter accounts for the high number density of massive galaxies at substantial redshifts z>1.5, the population of Extremely Red Objects, and the properties of local ellipticals.

We show that the ABC scenario we proposed for the co-evolution of spheroids and QSOs predicts accretion rates and masses of supermassive black holes in sub-mm galaxies in keeping with recent X-ray determinations. These masses are well below the local values, and those predicted by alternative models. The observed column densities may be mostly due to ISM in the galaxy. The contribution of the associated nuclear activity to the X-ray background is likely negligible, while they may contribute a sizeable fraction $\sim 10 %$ to hard-X cumulative counts at the faintest observed fluxes.

At variance with most semi-analytic models, in the Anti-hierarchical Baryon Collapse scenario (Granato et al. 2001, 2004) the main driver of the galaxy formation and evolution is not the merging sequence but are baryon processes. This approach emphasizes, still in the framework of the hierarchical clustering paradigm for dark matter halos, feedback processes from supernova explosions and from active nuclei, that tie together star formation in spheroidal galaxies and the growth of black holes at their centers. We review some recent results showing the remarkably successful predictive power of this scenario, which allows us to account for the evolution with cosmic time of a broad variety of properties of galaxies and active nuclei, which proved to be very challenging for competing models.

We present new relationships between halo masses M_h and several galaxy properties, including r*-band luminosities L_r, stellar (M_star) and baryonic masses, stellar velocity dispersions (sigma), and black hole masses (M_BH). Approximate analytic expressions are given. In the galaxy halo mass range 3*10^(10) M_sun <= M_h <= 3*10^(13) M_sun the M_h-L_r, M_star-M_h, and M_BH-M_h are well represented by a double power law, with a break at M_h_break~3*10^(11) M_sun, corresponding to a mass in stars M_star~ 1.2*10^(10) M_sun, to a r*-band luminosity L_r~ 5*10^9 L_sun, to a stellar velocity dispersion sigma ~ 88 km/s, and to a black hole mass M_BH~ 9*10^(6) M_sun. The sigma-M_h relation can be approximated by a single power law, though a double power law is a better representation. We interpret these relations in terms of the effect of feedback from supernovae and from the active nucleus on the interstellar medium. We argue that the break of the power laws occurs at a mass which marks the transition between the dominance of the stellar and the AGN feedback.

We compute the angular momentum, the spin parameter and the related distribution function for Dark Matter halos hosting a spiral galaxy. We base on scaling laws, inferred from observations, that link the properties of the galaxy to those of the host halo; we further assume that the Dark Matter has the same total specific angular momentum of the baryons. Our main results are: (i) we find that the gas component of the disk significantly contributes to the total angular momentum of the system; (ii) by adopting for the Dark Matter the observationally supported Burkert profile, we compute the total angular momentum of the disk and its correlation with the rotation velocity; (iii) we find that the distribution function of the spin parameter $\lambda$ peaks at a value of about 0.03, consistent with a no-major-merger scenario for the late evolution of spiral galaxies.

We investigate how accretion episodes onto massive black holes power quasars and AGN while they accumulate mass into the holes. We compute both the trend and the stochastic component to the trigger of the accretion events, as provided by structure buildup after the hierarchical paradigm. We base on host galaxy evolution proceeding from the protogalactic era at redshifts z > 2.5 dominated by major merging events in high density regions, to the subsequent era marked by galaxy-galaxy interactions in newly forming groups. These dynamical events perturb the gravitational equilibrium of the gas reservoir in the hosts, and trigger recurrent accretion episodes controlled by energy feedback from the very source emission onto the surrounding gas. Depletion of the latter by these events (adding to quiescent star formation) concurs with the slowing down of the clustering to cause a fast drop of the activity in dense regions. Meanwhile, in the "field" later and rarer events are triggered by interactions of still gas-rich galaxies, and eventually by captures of dwarf satellite galaxies; these are also included in our analytic model. Thus we compute the quasar and AGN luminosity functions; we find these to brighten and rise from z ~ 6 to z ~ 2.5, and then toward z ~ 0 to dim and fall somewhat, in detailed agreement with the observations. We predict that for z < 2.5 the mass distribution of the holes progressively rises and shifts rightwards; we compare our results with the local data. We also find that downward of z ~ 2.5 the Eddington ratios related to emitting, most massive holes drift below unity on average, with a widening scatter; meanwhile, some smaller holes flare up closer to the Eddington limit.

We show that the observed Velocity Dispersion Function of E/S0 galaxies matches strikingly well the distribution function of virial velocities of massive halos virializing at z > 1.5, as predicted by the standard hierarchical clustering scenario in a \LambdaCDM cosmology, for a constant ratio sigma/V_vir = 0.55 \pm 0.05, close to the value expected at virialization if it typically occurred at z > 3. This strongly suggests that dissipative processes and later merging events had little impact on the matter density profile. Adopting the above sigma/V_vir ratio, the observed relationships between photometric and dynamical properties which define the fundamental plane of elliptical galaxies, such as the luminosity-sigma (Faber-Jackson) and the luminosity-effective radius relations, as well as the M_BH-sigma relation, are nicely reproduced. Their shapes turn out to be determined by the mutual feedback of star-formation (and supernova explosions)and nuclear activity, along the lines discussed by Granato et al. (2004). To our knowledge, this is the first semi-analytic model for which simultaneous fits of the fundamental plane relations and of the epoch-dependent luminosity function of spheroidal galaxies have been presented.

We present a physical model for the coevolution of massive spheroidal galaxies and active nuclei at their centers. Supernova heating is increasingly effective in slowing down the star formation and in driving gas outflows in smaller and smaller dark matter halos. Thus the more massive protogalaxies virializing at early times are the sites of faster star formation. The correspondingly higher radiation drag causes a faster angular momentum loss by the gas and induces a larger accretion rate onto the central black hole. In turn, the kinetic energy of the outflows powered by the active nuclei can unbind the residual gas in a time shorter for larger halos. The model accounts for a broad variety of dynamical, photometric and metallicity properties of early-type galaxies, for the M_BH -- \sigma relation and for the local supermassive black-hole mass function.

We have performed a detailed analysis of the local super-massive black-hole (SMBH) mass function based on both kinematic and photometric data and derived an accurate analytical fit in the range 10^6 <= (M_BH/M_sun) <= 5*10^9. We find a total SMBH mass density of (4.2+/-1.1)*10^5 M_sun/Mpc^3, about 25% of which is contributed by SMBHs residing in bulges of late type galaxies. Exploiting up-to-date luminosity functions of hard X-ray and optically selected AGNs, we have studied the accretion history of the SMBH population. If most of the accretion happens at constant \dot{M_BH}/M_BH the local SMBH mass function is fully accounted for by mass accreted by X-ray selected AGNs, with bolometric corrections indicated by current observations and a standard mass-to-light conversion efficiency \epsilon ~10%. The analysis of the accretion history highlights that the most massive BHs accreted their mass faster and at higher redshifts (z>1.5), while the lower mass BHs responsible for most of the hard X-ray background have mostly grown at z<1.5. The accreted mass function matches the local SMBH mass function if \epsilon ~0.09(+0.04,-0.03) and the Eddington ratio \lambda=L/L_Edd \~0.3(+0.3,-0.1) (68% confidence errors). The visibility time, during which AGNs are luminous enough to be detected by the currently available X-ray surveys, ranges from ~0.1 Gyr for present day BH masses M_BH(z=0) ~10^6 M_sun to ~0.3 Gyr for M_BH(z=0) >= 10^9 M_sun. The mass accreted during luminous phases is >= 25-30% even if we assume extreme values of \epsilon (\epsilon \~0.3-0.4). An unlikely fine tuning of the parameters would be required to account for the local SMBH mass function accomodating a dominant contribution from 'dark' BH growth (due, e.g., to BH coalescence).