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.

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.

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.

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.

The two-point correlation function of the galaxy distribution is a key cosmological observable that allows us to constrain the dynamical and geometrical state of our Universe. To measure the correlation function we need to know both the galaxy positions and the expected galaxy density field. The expected field is commonly specified using a Monte-Carlo sampling of the volume covered by the survey and, to minimize additional sampling errors, this random catalog has to be much larger than the data catalog. Correlation function estimators compare data-data pair counts to data-random and random-random pair counts, where random-random pairs usually dominate the computational cost. Future redshift surveys will deliver spectroscopic catalogs of tens of millions of galaxies. Given the large number of random objects required to guarantee sub-percent accuracy, it is of paramount importance to improve the efficiency of the algorithm without degrading its precision. We show both analytically and numerically that splitting the random catalog into a number of subcatalogs of the same size as the data catalog when calculating random-random pairs, and excluding pairs across different subcatalogs provides the optimal error at fixed computational cost. For a random catalog fifty times larger than the data catalog, this reduces the computation time by a factor of more than ten without affecting estimator variance or bias.

We study the spatial clustering through the projected two-point correlation function of $632$ $(1130)$ XMM-COSMOS Active Galactic Nuclei (AGNs) with known spectroscopic (spectroscopic or photometric) redshifts in the range $z = [0.1 - 2.5]$ in order to measure the AGN bias and estimate the typical mass of the hosting dark matter (DM) halo as a function of AGN host galaxy properties. We create AGN subsamples in terms of stellar mass $M_*$ and specific black hole accretion rate $L_X/M_*$, to probe how AGN environment depends on these quantities. For the full spectroscopic AGN sample, we measure a typical DM halo mass of $\log (M_\mathrm{halo} / h^{-1}\mathrm{M}_\odot)= 12.79_{-0.43}^{+0.26}$, similar to galaxy group environments and in line with previous studies for moderate-luminosity X-ray selected AGN. We find no significant dependence on $L_X/M_*$, with $\log (M_\mathrm{halo} / h^{-1}\mathrm{M}_\odot) = 13.06_{-0.38}^{+0.23}$ ($12.97_{-1.26}^{+0.39}$) for the low (high) $L_X/M_*$ subsample. We also find no difference in the hosting halos in terms of $M_*$ with $\log (M_\mathrm{halo} / h^{-1}\mathrm{M}_\odot) = 12.93_{-0.62}^{+0.31}$ ($12.90_{-0.62}^{+0.30}$) for the low (high) $M_*$ subsample. By comparing the $M_*-M_\mathrm{halo}$ relation derived for XMM-COSMOS AGN subsamples with what is expected for normal non-active galaxies by abundance matching and clustering results, we find that the typical DM halo mass of our high $M_*$ AGN subsample is similar to that of non-active galaxies. However, AGNs in our low $M_*$ subsample are found in more massive halos than non-active galaxies. By excluding AGNs in galaxy groups from the clustering analysis, we find evidence that the result for low $M_*$ may be due a larger fraction of AGNs as satellites in massive halos.

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 characterize the environments of local accreting supermassive black holes by measuring the clustering of AGN in the Swift/BAT Spectroscopic Survey (BASS). With 548 AGN in the redshift range 0.01<z<0.1 over the full sky from the DR1 catalog, BASS provides the largest, least biased sample of local AGN to date due to its hard X-ray selection (14-195 keV) and rich multiwavelength/ancillary data. By measuring the projected cross-correlation function between the AGN and 2MASS galaxies, and interpreting it via halo occupation distribution (HOD) and subhalo-based models, we constrain the occupation statistics of the full sample, as well as in bins of absorbing column density and black hole mass. We find that AGN tend to reside in galaxy group environments, in agreement with previous studies of AGN throughout a large range of luminosity and redshift, and that on average they occupy their dark matter halos similar to inactive galaxies of comparable stellar mass. We also find evidence that obscured AGN tend to reside in denser environments than unobscured AGN, even when samples were matched in luminosity, redshift, stellar mass, and Eddington ratio. We show that this can be explained either by significantly different halo occupation distributions or statistically different host halo assembly histories. Lastly, we see that massive black holes are slightly more likely to reside in central galaxies than black holes of smaller mass.

We present results of a ground-based near-infrared campaign with Palomar TripleSpec, Keck NIRSPEC, and Gemini GNIRS to target two samples of reddened active galactic nucleus (AGN) candidates from the 31 deg$^2$ Stripe 82 X-ray survey. One sample, which is $\sim$89\% complete to $K<16$ (Vega), consists of eight confirmed AGNs, four of which were identified with our follow-up program, and is selected to have red $R-K$ colors ($>4$, Vega). The fainter sample ($K>17$, Vega) represents a pilot program to follow-up four sources from a parent sample of 34 that are not detected in the single-epoch SDSS catalog and have {\it WISE} quasar colors. All twelve sources are broad-line AGNs (at least one permitted emission line has a FWHM exceeding 1300 km s$^{-1}$) and span a redshift range $0.59 < z < 2.5$. Half the ($R-K$)-selected AGNs have features in their spectra suggestive of outflows. When comparing these sources to a matched sample of blue Type 1 AGNs, we find the reddened AGNs are more distant ($z > 0.5$) and a greater percentage have high X-ray luminosities ($L_{\rm X,full} > 10^{44}$ erg s$^{-1}$). Such outflows and high luminosities may be consistent with the paradigm that reddened broad-line AGNs represent a transitory phase in AGN evolution as described by the major merger model for black hole growth. Results from our pilot program demonstrate proof-of-concept that our selection technique is successful in discovering reddened quasars at $z > 1$ missed by optical surveys.

We study the X-ray variability properties of distant AGNs in the Chandra Deep Field-South region over 17 years, up to $z\sim 4$, and compare them with those predicted by models based on local samples. We use the results of Monte Carlo simulations to account for the biases introduced by the discontinuous sampling and the low-count regime. We confirm that variability is an ubiquitous property of AGNs, with no clear dependence on the density of the environment. The variability properties of high-z AGNs, over different temporal timescales, are most consistent with a Power Spectral Density (PSD) described by a broken (or bending) power-law, similar to nearby AGNs. We confirm the presence of an anti-correlation between luminosity and variability, resulting from the dependence of variability on BH mass and accretion rate. We explore different models, finding that our acceptable solutions predict that BH mass influences the value of the PSD break frequency, while the Eddington ratio $\lambda_{Edd}$ affects the PSD break frequency and, possibly, the PSD amplitude as well. We derive the evolution of the average $\lambda_{Edd}$ as a function of redshift, finding results in agreement with measurements based on different estimators. The large statistical uncertainties make our results consistent with a constant Eddington ratio, although one of our models suggest a possible increase of $\lambda_{Edd}$ with lookback time up to $z\sim 2-3$. We conclude that variability is a viable mean to trace the accretion history of supermassive BHs, whose usefulness will increase with future, wide-field/large effective area X-ray missions.

We investigate the star formation properties of a large sample of ~2300 X-ray-selected Type 2 Active Galactic Nuclei (AGNs) host galaxies out to z~3 in the Chandra COSMOS Legacy Survey in order to understand the connection between the star formation and nuclear activity. Making use of the existing multi-wavelength photometric data available in the COSMOS field, we perform a multi-component modeling from far-infrared to near-ultraviolet using a nuclear dust torus model, a stellar population model and a starburst model of the spectral energy distributions (SEDs). Through detailed analysis of SEDs, we derive the stellar masses and the star formation rates (SFRs) of Type 2 AGN host galaxies. The stellar mass of our sample is in the range 9 < log M_{stellar}/M_{\odot} < 12 with uncertainties of ~0.19 dex. We find that Type 2 AGN host galaxies have, on average, similar SFRs compared to the normal star-forming galaxies with similar M_{stellar} and redshift ranges, suggesting no significant evidence for enhancement or quenching of star formation. This could be interpreted in a scenario, where the relative massive galaxies have already experienced substantial growth at higher redshift (z>3), and grow slowly through secular fueling processes hosting moderate-luminosity AGNs.

Using {\em Chandra} observations in the 2.15 deg$^{2}$ COSMOS legacy field, we present one of the most accurate measurements of the Cosmic X-ray Background (CXB) spectrum to date in the [0.3-7] keV energy band. The CXB has three distinct components: contributions from two Galactic collisional thermal plasmas at kT$\sim$0.27 and 0.07 keV and an extragalactic power-law with photon spectral index $\Gamma$=1.45$\pm{0.02}$. The 1 keV normalization of the extragalactic component is 10.91$\pm{0.16}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. Removing all X-ray detected sources, the remaining unresolved CXB is best-fit by a power-law with normalization 4.18$\pm{0.26}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$ and photon spectral index $\Gamma$=1.57$\pm{0.10}$. Removing faint galaxies down to i$_{AB}\sim$27-28 leaves a hard spectrum with $\Gamma\sim$1.25 and a 1 keV normalization of $\sim$1.37 keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. This means that $\sim$91\% of the observed CXB is resolved into detected X-ray sources and undetected galaxies. Unresolved sources that contribute $\sim 8-9\%$ of the total CXB show a marginal evidence of being harder and possibly more obscured than resolved sources. Another $\sim$1\% of the CXB can be attributed to still undetected star forming galaxies and absorbed AGN. According to these limits, we investigate a scenario where early black holes totally account for non source CXB fraction and constrain some of their properties. In order to not exceed the remaining CXB and the $z\sim$6 accreted mass density, such a population of black holes must grow in Compton-thick envelopes with N$_{H}>$1.6$\times$10$^{25}$ cm$^{-2}$ and form in extremely low metallicity environments $(Z_\odot)\sim10^{-3}$.

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.

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$.

We present the X-ray spectral analysis of the 1855 extragalactic sources in the Chandra COSMOS-Legacy survey catalog having more than 30 net counts in the 0.5-7 keV band. 38% of the sources are optically classified Type 1 active galactic nuclei (AGN), 60% are Type 2 AGN and 2% are passive, low-redshift galaxies. We study the distribution of AGN photon index and of the intrinsic absorption N(H,z) based on the sources optical classification: Type 1 have a slightly steeper mean photon index than Type 2 AGN, which on the other hand have average intrinsic absorption ~3 times higher than Type 1 AGN. We find that ~15% of Type 1 AGN have N(H,z)>1E22 cm^(-2), i.e., are obscured according to the X-ray spectral fitting; the vast majority of these sources have L(2-10keV)>$1E44 erg/s. The existence of these objects suggests that optical and X-ray obscuration can be caused by different phenomena, the X-ray obscuration being for example caused by dust-free material surrounding the inner part of the nuclei. ~18% of Type 2 AGN have N(H,z)<1E22 cm^(-2), and most of these sources have low X-ray luminosities (L(2-10keV)<$1E43 erg/s). We expect a part of these sources to be low-accretion, unobscured AGN lacking of broad emission lines. Finally, we also find a direct proportional trend between N(H,z) and host galaxy mass and star formation rate, although part of this trend is due to a redshift selection effect.

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.

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 square degrees of the COSMOS field, combination of 56 new overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 micron identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 micron information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift towards faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGN and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.

Ultra-deep observations of ECDF-S with Chandra and XMM-Newton enable a search for extended X-ray emission down to an unprecedented flux of $2\times10^{-16}$ ergs s$^{-1}$ cm$^{-2}$. We present the search for the extended emission on spatial scales of 32$^{\prime\prime}$ in both Chandra and XMM data, covering 0.3 square degrees and model the extended emission on scales of arcminutes. We present a catalog of 46 spectroscopically identified groups, reaching a redshift of 1.6. We show that the statistical properties of ECDF-S, such as logN-logS and X-ray luminosity function are broadly consistent with LCDM, with the exception that dn/dz/d$\Omega$ test reveals that a redshift range of $0.2<z<0.5$ in ECDF-S is sparsely populated. The lack of nearby structure, however, makes studies of high-redshift groups particularly easier both in X-rays and lensing, due to a lower level of clustered foreground. We present one and two point statistics of the galaxy groups as well as weak-lensing analysis to show that the detected low-luminosity systems are indeed low-mass systems. We verify the applicability of the scaling relations between the X-ray luminosity and the total mass of the group, derived for the COSMOS survey to lower masses and higher redshifts probed by ECDF-S by means of stacked weak lensing and clustering analysis, constraining any possible departures to be within 30% in mass. Abridged.

We present the first measurement of the projected correlation function of 485 gamma-ray selected Blazars, divided in 175 BLLacertae (BL Lacs) and 310 Flat Spectrum Radio Quasars (FSRQs) detected in the 2-year all-sky survey by Fermi-Large Area Telescope. We find that Fermi BL Lacs and FSRQs reside in massive dark matter halos (DMHs) with logMh=13.35+0.20/-0.14 and logMh = 13.40+0.15/-0.19 Msun/h, respectively, at low (z=0.4) and high (z =1.2) redshift. In terms of clustering properties, these results suggest that BL Lacs and FSRQs are similar objects residing in the same dense environment typical of galaxy groups, despite their different spectral energy distribution, power and accretion rate. We find no difference in the typical bias and hosting halo mass between Fermi Blazars and radio-loud AGN, supporting the unifcation scheme simply equating radio-loud objects with misaligned Blazar counterparts. This similarity in terms of typical environment they preferentially live in, suggests that Blazars preferentially occupy the centre of DMHs, as already pointed out for radio-loud AGN. This implies, in light of several projects looking for the gamma-ray emission from DM annihilation in galaxy clusters, a strong contamination from Blazars to the expected signal from DM annihilation.

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 the statistical properties of the Normalized Excess Variance of variability process characterized by a red-noise power spectral density (PSD), as the case of Active Galactic Nuclei (AGN). We perform Monte Carlo simulations of lightcurves, assuming both a continuous and a sparse sampling pattern and various signal-to-noise (S/N) ratios. We show that the normalized excess variance is a biased estimate of the variance even in the case of continuously sampled lightcurves. The bias depends on the PSD slope and on the sampling pattern, but not on the S/N ratio. We provide a simple formula to account for the bias, which yields unbiased estimates with an accuracy better than 15%. We show that the normalized excess variance estimates based on single lightcurves (especially for sparse sampling and S/N less than 3) are highly uncertain (even if corrected for bias) and we propose instead the use of an ensemble estimate, based on multiple lightcurves of the same object, or on the use of lightcurves of many objects. These estimates have symmetric distributions, known errors, and can also be corrected for biases.We use our results to estimate the ability to measure the intrinsic source variability in current data, and show that they could also be useful in the planning of the observing strategy of future surveys such as those provided by X-ray missions studying distant and/or faint AGN populations and, more in general, in the estimation of the variability amplitude of sources that will result from future surveys such as Pan-STARRS, and LSST.

In this paper we investigate the power spectrum of the unresolved 0.5-2 keV CXB with deep Chandra 4 Ms observations in the CDFS. We measured a signal which, on scales >30", is significantly higher than the Shot-Noise and is increasing with the angular scale. We interpreted this signal as the joint contribution of clustered undetected sources like AGN, Galaxies and Inter-Galactic-Medium (IGM). The power of unresolved cosmic sources fluctuations accounts for \sim 12% of the 0.5-2 keV extragalactic CXB. Overall, our modeling predicts that \sim 20% of the unresolved CXB flux is made by low luminosity AGN, \sim 25% by galaxies and \sim 55% by the IGM (Inter Galactic Medium). We do not find any direct evidence of the so called Warm Hot Intergalactic Medium (i.e. matter with 10^5K<T<10^7K and density contrast {\delta} <1000), but we estimated that it could produce about 1/7 of the unresolved CXB. We placed an upper limit to the space density of postulated X-ray-emitting early black hole at z>7.5 and compared it with SMBH evolution models.

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.

The study of the angular and spatial structure of the X-ray sky has been under investigation since the times of the Einstein X-ray Observatory. This topic has fascinated more than two generations of scientists and slowly unveiled an unexpected scenario regarding the consequences of the angular and spatial distribution of X-ray sources. It was first established from the clustering of sources making the CXB that the source spatial distribution resembles that of optical QSO. It then it became evident that the distribution of X-ray AGN in the Universe was strongly reflecting that of Dark Matter. In particular one of the key result is that X-ray AGN are hosted by Dark Matter Halos of mass similar to that of galaxy groups. This result, together with model predictions, has lead to the hypothesis that galaxy mergers may constitute the main AGN triggering mechanism. However detailed analysis of observational data, acquired with modern telescopes, and the use of the new Halo Occupation formalism has revealed that the triggering of an AGN could also be attributed to phenomena like tidal disruption or disk instability, and to galaxy evolution. This paper reviews results from 1988 to 2011 in the field of X-ray selected AGN clustering.

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.

The Wide Field X-ray Telescope (WFXT) is a proposed mission with a high survey speed, due to the combination of large field of view (FOV) and effective area, i.e. grasp, and sharp PSF across the whole FOV. These characteristics make it suitable to detect a large number of variable and transient X-ray sources during its operating lifetime. Here we present estimates of the WFXT capabilities in the time domain, allowing to study the variability of thousand of AGNs with significant detail, as well as to constrain the rates and properties of hundreds of distant, faint and/or rare objects such as X-ray Flashes/faint GRBs, Tidal Disruption Events, ULXs, Type-I bursts etc. The planned WFXT extragalactic surveys will thus allow to trace variable and transient X-ray populations over large cosmological volumes.

Recent work on nearby AGNs has shown that X-ray variability is correlated with the mass and accretion rate onto the central SMBH. Here we present the application of the variability-luminosity relation to high redshift AGNs in the CDFS, making use of XMM-Newton observations. We use Monte Carlo simulations in order to properly account for bias and uncertainties introduced by the sparse sampling and the very low statistics. Our preliminary results indicate that BH masses span over the range from 10^5 to 10^9 solar mass while accretion rates range from 10^-3 up to values greater than 1, in unit of Eddington accretion rate.