The shape of dark matter haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in dark matter haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the peak-centroid offset. In particular, the correlation between halo ellipticity and halo concentration is nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multi-dimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-size halos. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields.

The clustering of active galactic nuclei (AGN) sheds light on their typical large (Mpc-scale) environments, which can constrain the growth and evolution of supermassive black holes. Here we measure the clustering of luminous X-ray-selected AGN in the Stripe 82X and XMM-XXL-North surveys around the peak epoch of black hole growth, in order to investigate the dependence of luminosity on large-scale AGN environment. We compute the auto-correlation function of AGN in two luminosity bins, $10^{43}\leq L_X<10^{44.5}$ erg s$^{-1}$ at $z\sim 0.8$ and $L_X\geq 10^{44.5}$ erg s$^{-1}$ at $z\sim 1.8$, and calculate the AGN bias taking into account the redshift distribution of the sources using three different methods. Our results show that while the less luminous sample has an inferred typical halo mass that is smaller than for the more luminous AGN, the host halo mass may be less dependent on luminosity than suggested in previous work. Focusing on the luminous sample, we calculate a typical host halo mass of $\sim 10^{13}$ M$_{\odot}~h^{-1}$, which is similar to previous measurements of moderate-luminosity X-ray AGN and significantly larger than the values found for optical quasars of similar luminosities and redshifts. We suggest that the clustering differences between different AGN selection techniques are dominated by selection biases, and not due to a dependence on AGN luminosity. We discuss the limitations of inferring AGN triggering mechanisms from halo masses derived by large-scale bias.

State-of-the-art models of massive black hole formation postulate that quasars at $z>6$ reside in extreme peaks of the cosmic density structure in the early universe. Even so, direct observational evidence of these overdensities is elusive, especially on large scales ($\gg$1 Mpc) as the spectroscopic follow-up of $z>6$ galaxies is observationally expensive. Here we present Keck / DEIMOS optical and IRAM / NOEMA millimeter spectroscopy of a $z\sim6$ Lyman-break galaxy candidate originally discovered via broadband selection, at a projected separation of 4.65 physical Mpc (13.94 arcmin) from the luminous $z$=6.308 quasar J1030+0524. This well-studied field presents the strongest indication to date of a large-scale overdensity around a $z>6$ quasar. The Keck observations suggest a $z\sim6.3$ dropout identification of the galaxy. The NOEMA 1.2mm spectrum shows a 3.5$\sigma$ line that, if interpreted as [CII], would place the galaxy at $z$=6.318 (i.e., at a line-of-sight separation of 3.9 comoving Mpc assuming that relative proper motion is negligible). The measured [CII] luminosity is $3\times10^8$ L$_\odot$, in line with expectations for a galaxy with a star formation rate $\sim15$ M$_\odot$ yr$^{-1}$, as inferred from the rest-frame UV photometry. Our combined observations place the galaxy at the same redshift as the quasar, thus strengthening the overdensity scenario for this $z>6$ quasar. This pilot experiment demonstrates the power of millimeter-wavelength observations in the characterization of the environment of early quasars.

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.

We report the discovery of a galaxy overdensity around a Compton-thick Fanaroff-Riley type II (FRII) radio galaxy at z=1.7 in the deep multiband survey around the z=6.3 QSO SDSS J1030+0524. Based on a 6hr VLT/MUSE and on a 4hr LBT/LUCI observation, we identify at least eight galaxy members in this structure with spectroscopic redshift z=1.687-1.699, including the FRII galaxy at z=1.699. Most of the identified overdensity members are blue, compact galaxies that are actively forming stars at rates of $\sim$8-60 $M_{\odot}$ yr$^{-1}$. Based on a 500ks Chandra ACIS-I observation we found that the FRII nucleus hosts a luminous QSO ($L_{2-10keV}=1.3\times10^{44}$ erg s$^{-1}$, intrinsic and rest-frame) that is obscured by Compton-thick absorption ($N_H=1.5\pm0.6\times 10^{24}$ cm$^{-2}$). Our Chandra observation, the deepest so far for a distant FRII within a galaxy overdensity, revealed significant diffuse X-ray emission within the region covered by the overdensity. In particular, X-ray emission extending for $\sim$240 kpc is found around the Eastern lobe of the FRII. Four out of the six MUSE star forming galaxies in the overdensity are distributed in an arc-like shape at the edge of this diffuse X-ray emission. The probability of observing by chance four out of the six $z=1.7$ sources at the edge of the diffuse emission is negligible. In addition, these four galaxies have the highest specific star formation rates of the MUSE galaxies in the overdensity and lie above the main sequence of field galaxies of equal stellar mass at z=1.7. We propose that the diffuse X-rays originate from an expanding bubble of gas that is shock-heated by the FRII jet, and that star formation is promoted by the compression of the cold interstellar medium of the galaxies around the bubble, which may be remarkable evidence of positive AGN feedback on cosmological scales. [shortened version]

We study the coherence of the near-infrared and X-ray background fluctuations and the X-ray spectral properties of the sources producing it. We use data from multiple Spitzer and Chandra surveys, including the UDS/SXDF surveys, the Hubble Deep Field North, the EGS/AEGIS field, the Chandra Deep Field South and the COSMOS surveys, comprising $\sim$2275 Spitzer/IRAC hours and $\sim$~16 Ms of Chandra data collected over a total area of $\sim$~1~deg$^2$. We report an overall $\sim$5$\sigma$ detection of a cross-power signal on large angular scales $>$ 20$''$ between the 3.6 and 4.5\mum\ and the X-ray bands, with the IR vs [1-2] keV signal detected at 5.2$\sigma$. The [0.5-1] and [2-4] keV bands are correlated with the infrared wavelengths at a $\sim$1$-$3$\sigma$ significance level. The hardest X-ray band ([4-7] keV) alone is not significantly correlated with any infrared wavelengths due to poor photon and sampling statistics. We study the X-ray SED of the cross-power signal. We find that its shape is consistent with a variety of source populations of accreting compact objects, such as local unabsorbed AGNs or high-z absorbed sources. We cannot exclude that the excess fluctuations are produced by more than one population. Because of poor statistics, the current relatively broad photometric bands employed here do not allow distinguishing the exact nature of these compact objects or if a fraction of the fluctuations have instead a local origin.

There exist hitherto unexplained fluctuations in the Cosmic Infrared Background (CIB) on arcminute scales and larger. These have been shown to cross-correlate with the Cosmic X-ray Background (CXB), leading several authors to attribute the excess to a high-redshift growing black hole population. In order to investigate potential sources that could explain this excess, in this paper, we develop a new framework to compute the power spectrum of undetected sources that do not have constant flux as a function of halo mass. In this formulation, we combine a semi-analytic model for black hole growth and their simulated spectra from hydrodynamical simulations. Revisiting the possible contribution of a high-redshift black hole population, we find that too much black hole growth is required at early epochs for z>6 accretion to explain these fluctuations. Examining a population of accreting black holes at more moderate redshifts, z\sim 2-3, we find that such models produce a poor fit to the observed fluctuations while simultaneously overproducing the local black hole mass density. Additionally, we rule out the hypothesis of a missing Galactic foreground of warm dust that produces coherent fluctuations in the X-ray via reflection of Galactic X-ray binary emission. Although we firmly rule out accreting massive black holes as the source of these missing fluctuations, additional studies will be required to determine their origin.

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.

While theoretical arguments predict that most of the early growth of supermassive black holes (SMBHs) happened during heavily obscured phases of accretion, current methods used for selecting $z>6$ quasars (QSOs) are strongly biased against obscured QSOs, thus considerably limiting our understanding of accreting SMBHs during the first Gyr of the Universe from an observational point of view. We report the $Chandra$ discovery of the first heavily obscured QSO candidate in the early universe, hosted by a close ($\approx5$ kpc) galaxy pair at $z=6.515$. One of the members is an optically classified type 1 QSO, PSO167-13. The companion galaxy was first detected as a [C II] emitter by ALMA. An X-ray source is significantly ($P=0.9996$) detected by $Chandra$ in the 2-5 keV band, with $<1.14$ net counts in the 0.5-2 keV band, although the current positional uncertainty does not allow a conclusive association with either PSO167-13 or its companion galaxy. From X-ray photometry and hardness-ratio arguments, we estimated an obscuring column density of $N_H>2\times10^{24}\,\mathrm{cm^{-2}}$ and $N_H>6\times10^{23}\,\mathrm{cm^{-2}}$ at $68\%$ and $90\%$ confidence levels, respectively. Thus, regardless of which of the two galaxies is associated with the X-ray emission, this source is the first heavily obscured QSO candidate at $z>6$.

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

In this white paper we explore the capabilities required to identify and study supermassive black holes formed from heavy seeds ($\mathrm{M_{\bullet}} \sim 10^4 - 10^6 \, \mathrm{M_{\odot}}$) in the early Universe. To obtain an unequivocal detection of heavy seeds we need to probe mass scales of $\sim 10^{5-6} \, \mathrm{M_{\odot}}$ at redshift $z \gtrsim 10$. From this theoretical perspective, we review the observational requirements and how they compare with planned/proposed instruments, in the infrared, X-ray and gravitational waves realms. In conclusion, detecting heavy black hole seeds at $z \gtrsim 10$ in the next decade will be challenging but, according to current theoretical models, feasible with upcoming/proposed facilities. Their detection will be fundamental to understand the early history of the Universe, as well as its evolution until now. Shedding light on the dawn of black holes will certainly be one of the key tasks that the astronomical community will focus on in the next decade.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The source-subtracted cosmic infrared background (CIB) fluctuations uncovered in deep Spitzer data cannot be explained by known galaxy populations and appear strongly coherent with unresolved cosmic X-ray background (CXB). This suggests that the source-subtracted CIB contains emissions from significantly abundant accreting black holes (BHs). We show that theoretically such populations would have the angular power spectrum which is largely independent of the epochs occupied by these sources, provided they are at z>~ 4, offering an important test of the origin of the new populations. Using the current measurements we reconstruct the underlying soft X-ray CXB from the new sources and show that its fluctuations, while consistent with a high-z origin, have an amplitude that cannot be reached in direct measurements with the foreseeable X-ray space missions. This necessitates application of the methods developed by the authors to future IR and X-ray datasets, which must cover large areas of the sky in order to measure the signal with high precision. The LIBRAE project within ESA's Euclid mission will probe source-subtracted CIB over ~1/2 the sky at three near-IR bands, and its cross-power with unresolved CXB can be measured then from the concurrent eROSITA mission covering the same areas of the sky. We discuss the required methodology for this measurement and evaluate its projected S/N to show the unique potential of this experimental configuration to accurately probe the CXB from the new BH sources and help identify their epochs.

We present the results of a search for galaxy clusters and groups in the $\sim2$ square degree of the COSMOS field using all available X-ray observations from the XMM-Newton and Chandra observatories. We reach an X-ray flux limit of $3\times10^{-16}\;ergs\;cm^{-2}\;s^{-1}$ in 0.5--2 keV range, and identify 247 X-ray groups with $M_{200c}=8\times10^{12}-3\times10^{14}\;M_{\odot}$ at a redshift range of $0.08\leq z<1.53$, using the multiband photometric redshift and the master spectroscopic redshift catalogues of the COSMOS. The X-ray centres of groups are determined using high-resolution Chandra imaging. We investigate the relations between the offset of the brightest group galaxies (BGGs) from halo X-ray centre and group properties and compare with predictions from semi-analytic models and hydrodynamical simulations. We find that BGG offset decreases with both increasing halo mass and decreasing redshift with no strong dependence on the X-ray flux and SNR. We show that the BGG offset decreases as a function of increasing magnitude gap with no considerable redshift dependent trend. The stellar mass of BGGs in observations extends over a wider dynamic range compared to model predictions. At $z<0.5$, the central dominant BGGs become more massive than those with large offsets by up to 0.3dex, in agreement with model prediction. The observed and predicted lognormal scatter in the stellar mass of both low- and large-offset BGGs at fixed halo mass is $\sim0.3$dex.

We present the X-UDS survey, a set of wide and deep Chandra observations of the Subaru-XMM Deep/UKIDSS Ultra Deep Survey (SXDS/UDS) field. The survey consists of 25 observations that cover a total area of 0.33 deg$^{2}$. The observations are combined to provide a nominal depth of ~600 ksec in the central 100 arcmin$^{2}$ region of the field that has been imaged with Hubble/WFC3 by the CANDELS survey and $\sim$200 ksec in the remainder of the field. In this paper, we outline the survey's scientific goals, describe our observing strategy, and detail our data reduction and point source detection algorithms. Our analysis has resulted in a total of 868 band-merged point sources detected with a false-positive Poisson probability of $<1\times10^{-4}$. In addition, we present the results of an X-ray spectral analysis and provide best-fitting neutral hydrogen column densities, $N_{\rm H}$, as well as a sample of 51 Compton-thick active galactic nucleus candidates. Using this sample, we find the intrinsic Compton-thick fraction to be 30-35% over a wide range in redshift ($z=0.1-3$), suggesting the obscured fraction does not evolve very strongly with epoch. However, if we assume that the Compton-thick fraction is dependent on luminosity, as is seen for Compton-thin sources, then our results are consistent with a rise in the obscured fraction out to $z\sim3$. Finally, an examination of the host morphologies of our Compton-thick candidates shows a high fraction of morphological disturbances, in agreement with our previous results. All data products described in this paper are made available via a public website.

We study the source-subtracted near-infrared and X-ray background fluctuations of the COSMOS field using data from the Spitzer SPLASH program ($\sim$1272 hours) and Chandra COSMOS Legacy Survey (4.6 Ms). The new auto power spectra of the cosmic infrared and X-ray background fluctuations reach maximum angular scales of $\sim$ 3000$''$ and $\sim$ 5000$''$, respectively. We measure the cross power spectra between each infrared and X-ray band and calculate the mean power above 20$''$. We find that the soft X-ray band is correlated with 3.6 and 4.5$\mu$m at $\sim$ 4 $\sigma$ significance level. The significance between hard X-ray and the 3.6$\mu$m (4.5$\mu$m) band is $\sim$2.2 $\sigma$ ($\sim$3.8$\sigma$). The combined infrared (3.6 + 4.5$\mu$m) data are correlated with the X-ray data in soft ([0.5-2] keV), hard ([2-7] keV) and broad ([0.5-7] keV) bands at $\sim$5.6$\sigma$, $\sim$4.4$\sigma$ and $\sim$6.6$\sigma$ level, respectively. We compare the new measurements with existing models for the contributions from known populations at $z$$<$7, which are not subtracted. The model predictions are consistent with the measurements but we cannot rule out contributions from other components, such as Direct Collapse Black Holes (DCBH). However, the stacked cross-power spectra, combining other available data, show excess fluctuations about an order of magnitude on average at $\sim$4$\sigma$ confidence at scales within $\sim$300$''$. By studying the X-ray SED of the cross-power signal, assuming no significant variation from the infrared, we find that its shape is consistent with DCBHs.

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.

The cosmic infrared background (CIB) contains emissions accumulated over the entire history of the Universe, including from objects inaccessible to individual telescopic studies. The near-IR (~1-10 mic) part of the CIB, and its fluctuations, reflects emissions from nucleosynthetic sources and gravitationally accreting black holes (BHs). If known galaxies are removed to sufficient depths the source-subtracted CIB fluctuations at near-IR can reveal sources present in the first-stars-era and possibly new stellar populations at more recent times. This review discusses the recent progress in this newly emerging field which identified, with new data and methodology, significant source-subtracted CIB fluctuations substantially in excess of what can be produced by remaining known galaxies. The CIB fluctuations further appear coherent with unresolved cosmic X-ray background (CXB) indicating a very high fraction of BHs among the new sources producing the CIB fluctuations. These observations have led to intensive theoretical efforts to explain the measurements and their properties. While current experimental configurations have limitations in decisively probing these theories, their potentially remarkable implications will be tested in the upcoming CIB measurements with the ESA's Euclid dark energy mission. We describe the goals and methodologies of LIBRAE (Looking at Infrared Background Radiation with Euclid), a NASA-selected project for CIB science with Euclid, which has the potential for transforming the field into a new area of precision cosmology.

We present the results from a $\sim500$ ks Chandra observation of the $z=6.31$ QSO SDSS J1030+0524. This is the deepest X-ray observation to date of a $z\sim6$ QSO. The QSO is detected with a total of 125 net counts in the full ($0.5-7$ keV) band and its spectrum can be modeled by a single power-law model with photon index of $\Gamma = 1.81 \pm 0.18$ and full band flux of $f=3.95\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$. When compared with the data obtained by XMM-Newton in 2003, our Chandra observation in 2017 shows a harder ($\Delta \Gamma \approx -0.6$) spectrum and a 2.5 times fainter flux. Such a variation, in a timespan of $\sim2$ yrs rest-frame, is unexpected for such a luminous QSO powered by a $> 10^9 \: M_{\odot}$ black hole. The observed source hardening and weakening could be related to an intrinsic variation in the accretion rate. However, the limited photon statistics does not allow us to discriminate between an intrinsic luminosity and spectral change, and an absorption event produced by an intervening gas cloud along the line of sight. We also report the discovery of diffuse X-ray emission that extends for 30"x20" southward the QSO with a signal-to-noise ratio of $\sim$6, hardness ratio of $HR=0.03_{-0.25}^{+0.20}$, and soft band flux of $f_{0.5-2 \: keV}= 1.1_{-0.3}^{+0.3} \times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$, that is not associated to a group or cluster of galaxies. We discuss two possible explanations for the extended emission, which may be either associated with the radio lobe of a nearby, foreground radio galaxy (at $z \approx 1-2$), or ascribed to the feedback from the QSO itself acting on its surrounding environment, as proposed by simulations of early black hole formation.

In this paper we report a systematic search for an emission line around 3.5 keV in the spectrum of the Cosmic X-ray Background using a total of $\sim$10 Ms Chandra observations towards the COSMOS Legacy and CDFS survey fields. We find a marginal evidence of a feature at an energy of $\sim$3.51 keV with a significance of 2.5-3 $\sigma$, depending on the choice of the statistical treatment. The line intensity is best fit at $8.8\ \pm\ {2.9}\times10^{-7}$ ph cm$^{-2}$s$^{-1}$ when using a simple $\Delta\chi^2$ or $10.2\ ^{+0.2}_{-0.4} \times10^{-7}$ ph cm$^{-2}$s$^{-1}$ when MCMC is used. Based on our knowledge of $Chandra$, and the reported detection of the line by other instruments, an instrumental origin for the line remains unlikely. We cannot though rule out a statistical fluctuation and in that case our results provide a 3$\sigma$ upper limit at 1.85$\times$10$^{-6}$ ph cm$^{-2}$s$^{-1}$. We discuss the interpretation of this observed line in terms of the iron line background; S {\sc XVI} charge exchange as well as potentially from sterile neutrino decay. We note that our detection is consistent with previous measurements of this line toward the Galactic center, and can be modeled as the result of sterile neutrino decay from the Milky Way for the dark matter distribution modeled as an NFW profile. For this case, we estimate a mass m$_{\nu}\sim$7.01 keV and a mixing angle sin$^2$(2$\theta$)= 0.83--2.75 $\times 10^{-10}$. These derived values are in agreement with independent estimates from galaxy clusters; the Galactic center and M31.

Super-massive black holes weighing up to $\sim 10^9 \, \mathrm{M_{\odot}}$ are in place by $z \sim 7$, when the age of the Universe is $\lesssim 1 \, \mathrm{Gyr}$. This implies a time crunch for their growth, since such high masses cannot be easily reached in standard accretion scenarios. Here, we explore the physical conditions that would lead to optimal growth wherein stable super-Eddington accretion would be permitted. Our analysis suggests that the preponderance of optimal conditions depends on two key parameters: the black hole mass and the host galaxy central gas density. In the high-efficiency region of this parameter space, a continuous stream of gas can accrete onto the black hole from large to small spatial scales, assuming a global isothermal profile for the host galaxy. Using analytical initial mass functions for black hole seeds, we find an enhanced probability of high-efficiency growth for seeds with initial masses $\gtrsim 10^4 \, \mathrm{M_{\odot}}$. Our picture suggests that a large population of high-$z$ lower-mass black holes that formed in the low-efficiency region, with low duty cycles and accretion rates, might remain undetectable as quasars, since we predict their bolometric luminosities to be $\lesssim 10^{41} \, \mathrm{erg \, s^{-1}}$. The presence of these sources might be revealed only via gravitational wave detections of their mergers.

Multi-wavelength surveys covering large sky volumes are necessary to obtain an accurate census of rare objects such as high luminosity and/or high redshift active galactic nuclei (AGN). Stripe 82X is a 31.3 deg$^2$ X-ray survey with $Chandra$ and $XMM$-Newton observations overlapping the legacy Sloan Digital Sky Survey (SDSS) Stripe 82 field, which has a rich investment of multi-wavelength coverage from the ultraviolet to the radio. The wide-area nature of this survey presents new challenges for photometric redshifts for AGN compared to previous work on narrow-deep fields because it probes different populations of objects that need to be identified and represented in the library of templates. Here we present an updated X-ray plus multi-wavelength matched catalog, including $Spitzer$ counterparts, and estimated photometric redshifts for 5961 (96% of a total of 6181) X-ray sources, which have a normalized median absolute deviation, $\sigma_{\rm nmad}$ = 0.06 and an outlier fraction, $\eta$ = 13.7%. The populations found in this survey, and the template libraries used for photometric redshifts, provide important guiding principles for upcoming large-area surveys such as $eROSITA$ and 3$XMM$ (in X-ray) and the Large Synoptic Survey Telescope (LSST; optical).

We present new measurements of the large scale clustering component of the cross-power spectra of the source-subtracted Spitzer-IRAC Cosmic Infrared Background (CIB) and Chandra-ACIS Cosmic X-ray Background (CXB) surface brightness fluctuations. Our investigation uses data from the Chandra Deep Field South (CDFS), Hubble Deep Field North (HDFN), EGS/AEGIS field and UDS/SXDF surveys, comprising 1160 Spitzer hours and $\sim$ 12 Ms of Chandra data collected over a total area of 0.3 deg$^2$. We report the first ($>$5$\sigma$) detection of a cross-power signal on large angular scales >20 arcsec between [0.5-2] keV and the 3.6 and 4.5$\mu$m bands, at $\sim$5$\sigma$ and 6.3$\sigma$ significance, respectively. The correlation with harder X-ray bands is marginally significant. Comparing the new observations with existing models for the contribution of the known unmasked source population at $z<$7, we find an excess of about an order of magnitude at 5$\sigma$ confidence. We discuss possible interpretations for the origin of this excess in terms of the contribution from accreting early black holes, including both direct collapse black holes and primordial black holes, as well as from scattering in the interstellar medium and intra-halo light.

Many cosmological studies predict that early supermassive black holes (SMBHs) can only form in the most massive dark matter halos embedded within large scale structures marked by galaxy over-densities that may extend up to 10 physical Mpc. This scenario, however, has not been confirmed observationally, as the search for galaxy over-densities around high-z quasars has returned conflicting results. The field around the z=6.28 quasar SDSSJ1030+0524 (J1030) is unique for multi-band coverage and represents an excellent data legacy for studying the environment around a primordial SMBH. In this paper we present wide-area (25x25 arcmin), Y- and J-band imaging of the J1030 field obtained with the near infrared camera WIRCam at the Canada-France-Hawaii Telescope (CFHT). We built source catalogues in the Y- and J-band, and matched those with our photometric catalogue in the r, z, i bands presented in Morselli et al. (2014). We used these new infrared data together with H and K and Spitzer/IRAC data to refine our selection of Lyman Break Galaxies (LBGs), extending our selection criteria to galaxies in the range 25.2<zAB<25.7. We selected 21 robust high-z candidates in the J1030 field with photometric redshift around 6 and colors i-z>=1.3. We found a significant asymmetry in the distribution of the high-z galaxies in J1030, supporting the existence of a coherent large-scale structure around the quasar. We compared our results with those of Bowler et al. (2015), who adopted similar LBGs selection criteria, and estimated an over-density of galaxies in the field of delta = 2.4, which is significant at >4 sigma. The over-density value and its significance are higher than those found in Morselli et al. (2014), and we interpret this as evidence of an improved LBG selection.

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.

We report the detection of morphology dependent stellar age in massive quenched galaxies (QGs) at z~1.2. The sense of the dependence is that compact QGs are 0.5-2 Gyr older than normal-sized ones. The evidence comes from three different age indicators, Dn4000, Hdelta and fits to spectral synthesis models, applied to their stacked optical spectra. All age indicators consistently show that the stellar populations of compact QGs are older than their normally-sized counterparts. We detect weak [OII] emission in a fraction of QGs, and the strength of the line, when present, is similar between the two samples; however, compact galaxies exhibit significantly lower frequency of [OII] emission than normal ones. A fraction of both samples are individually detected in 7 Ms Chandra X-ray images (luminosities$\sim10^{40}$-$10^{41}$ erg/sec). 7 Ms stacks of non-detected galaxies show similarly low luminosities in the soft band only, consistent with a hot gas origin for the X-ray emission. While both [OII] emitters and non-emitters are also X-ray sources among normal galaxies, no compact galaxy with [OII] emission is an X-ray source, arguing against an AGN powering the line in compact galaxies. We interpret the [OII] properties as further evidence that compact galaxies are older and further along into the process of quenching star-formation and suppressing gas accretion. Finally, we argue that the older age of compact QGs is evidence of progenitor bias: compact QGs simply reflect the smaller sizes of galaxies at their earlier quenching epoch, with stellar density most likely having nothing directly to do with cessation of star-formation.

Growing the supermassive black holes (~10^9 Msun) that power the detected luminous, highest redshift quasars (z > 6) from light seeds - the remnants of the first stars - within ~ 1 Gyr of the Big Bang poses a timing challenge for growth models. The formation of massive black hole seeds via direct collapse with initial masses ~ 10^4 - 10^5 Msun alleviates this problem. Physical conditions required to form these massive direct collapse black hole (DCBH) seeds are available in the early universe. These viable DCBH formation sites, satellite halos of star-forming galaxies, merge and acquire a stellar component. These produce a new, transient class of objects at high redshift, Obese Black hole Galaxies (OBGs), where the luminosity produced by accretion onto the black hole outshines the stellar component. Therefore, the OBG stage offers a unique way to discriminate between light and massive initial seeds. We predict the multi-wavelength energy output of OBGs and growing Pop III remnants at a fiducial redshift (z = 9), exploring both standard and slim disk accretion onto the growing central black hole for high and low metallicities of the associated stellar population. With our computed templates, we derive the selection criteria for OBGs, that comprise a pre-selection that eliminates blue sources; followed by color-color cuts ([F_{070W} - F_{220W}] > 0; -0.3 < [F_{200W} - F_{444W}] < 0.3) and when available, a high ratio of X-ray flux to rest-frame optical flux (F_X/F_{444W} >> 1) (Abridged).

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 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 report the discovery of a remarkable concentration of massive galaxies with extended X-ray emission at $z_{spec} = 2.506$, which contains 11 massive ($M_{*} \gtrsim 10^{11} M_{\odot}$) galaxies in the central 80kpc region (11.6$\sigma$ overdensity). We have spectroscopically confirmed 17 member galaxies with 11 from CO and the remaining ones from $H\alpha$. The X-ray luminosity, stellar mass content and velocity dispersion all point to a collapsed, cluster-sized dark matter halo with mass $M_{200c} = 10^{13.9\pm0.2} M_{\odot}$, making it the most distant X-ray-detected cluster known to date. Unlike other clusters discovered so far, this structure is dominated by star-forming galaxies (SFGs) in the core with only 2 out of the 11 massive galaxies classified as quiescent. The star formation rate (SFR) in the 80kpc core reaches $\sim$3400 $M_{\odot}$ yr$^{-1}$ with a gas depletion time of $\sim 200$ Myr, suggesting that we caught this cluster in rapid build-up of a dense core. The high SFR is driven by both a high abundance of SFGs and a higher starburst fraction ($\sim25\%$, compared to 3\%-5\% in the field). The presence of both a collapsed, cluster-sized halo and a predominant population of massive SFGs suggests that this structure could represent an important transition phase between protoclusters and mature clusters. It provides evidence that the main phase of massive galaxy passivization will take place after galaxies accrete onto the cluster, providing new insights into massive cluster formation at early epochs. The large integrated stellar mass at such high redshift challenges our understanding of massive cluster formation.

We exploit the 7 Ms \textit{Chandra} observations in the \chandra\,Deep Field-South (\mbox{CDF-S}), the deepest X-ray survey to date, coupled with CANDELS/GOODS-S data, to measure the total X-ray emission arising from 2076 galaxies at $3.5\leq z < 6.5$. This aim is achieved by stacking the \textit{Chandra} data at the positions of optically selected galaxies, reaching effective exposure times of $\geq10^9\mathrm{s}$. We detect significant ($>3.7\sigma$) X-ray emission from massive galaxies at $z\approx4$. We also report the detection of massive galaxies at $z\approx5$ at a $99.7\%$ confidence level ($2.7\sigma$), the highest significance ever obtained for X-ray emission from galaxies at such high redshifts. No significant signal is detected from galaxies at even higher redshifts. The stacking results place constraints on the BHAD associated with the known high-redshift galaxy samples, as well as on the SFRD at high redshift, assuming a range of prescriptions for X-ray emission due to X- ray binaries. We find that the X-ray emission from our sample is likely dominated by processes related to star formation. Our results show that low-rate mass accretion onto SMBHs in individually X-ray-undetected galaxies is negligible, compared with the BHAD measured for samples of X-ray detected AGN, for cosmic SMBH mass assembly at high redshift. We also place, for the first time, constraints on the faint-end of the AGN X-ray luminosity function ($\mathrm{logL_X\sim42}$) at $z>4$, with evidence for fairly flat slopes. The implications of all of these findings are discussed in the context of the evolution of the AGN population at high redshift.

Improving the capabilities of detecting faint X-ray sources is fundamental to increase the statistics on faint high-z AGN and star-forming galaxies. We performed a simultaneous Maximum Likelihood PSF fit in the [0.5-2] keV and [2-7] keV energy bands of the 4 Ms{\em Chandra} Deep Field South (CDFS) data at the position of the 34930 CANDELS H-band selected galaxies. For each detected source we provide X-ray photometry and optical counterpart validation. We validated this technique by means of a raytracing simulation. We detected a total of 698 X-ray point-sources with a likelihood $\mathcal{L}$$>$4.98 (i.e. $>$2.7$\sigma$). We show that the prior knowledge of a deep sample of Optical-NIR galaxies leads to a significant increase of the detection of faint (i.e. $\sim$10$^{-17}$ cgs in the [0.5-2] keV band) sources with respect to "blind" X-ray detections. By including previous X-ray catalogs, this work increases the total number of X-ray sources detected in the 4 Ms CDFS, CANDELS area to 793, which represents the largest sample of extremely faint X-ray sources assembled to date. Our results suggest that a large fraction of the optical counterparts of our X-ray sources determined by likelihood ratio actually coincides with the priors used for the source detection. Most of the new detected sources are likely star-forming galaxies or faint absorbed AGN. We identified a few sources sources with putative photometric redshift z$>$4. Despite the low number statistics and the uncertainties on the photo-z, this sample significantly increases the number of X--ray selected candidate high-z AGN.

The XMM-LSS, XMM-COSMOS, and XMM-CDFS surveys are complementary in terms of sky coverage and depth. Together, they form a clean sample with the least possible variance in instrument effective areas and PSF. Therefore this is one of the best samples available to determine the 2-10 keV luminosity function of AGN and its evolution. The samples and the relevant corrections for incompleteness are described. A total of 2887 AGN is used to build the LF in the luminosity interval 10^42-10^46 erg/s, and in the redshift interval 0.001-4. A new method to correct for absorption by considering the probability distribution for the column density conditioned on the hardness ratio is presented. The binned luminosity function and its evolution is determined with a variant of the Page-Carrera method, improved to include corrections for absorption and to account for the full probability distribution of photometric redshifts. Parametric models, namely a double power-law with LADE or LDDE evolution, are explored using Bayesian inference. We introduce the Watanabe-Akaike information criterion (WAIC) to compare the models and estimate their predictive power. Our data are best described by the LADE model, as hinted by the WAIC indicator. We also explore the 15-parameter extended LDDE model recently proposed by Ueda et al., and find that this extension is not supported by our data. The strength of our method is that it provides: un-absorbed non-parametric estimates; credible intervals for luminosity function parameters; model choice according to which one has more predictive power for future data.

We present multiwavelength photometric catalogues (HST, Spitzer and Hawk-I K band) for the first two of the Frontier Fields, Abell2744 and MACSJ0416 (plus their parallel fields). To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members as well as the ICL. We perform the detection on the HST H160 processed image to obtain a pure H-selected sample. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry in the other HST bands is obtained using SExtractor, performed again on residual images after the Galfit procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands has been obtained using our PSF-matching deconfusion code T-PHOT. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields. The procedure allows for the detection and the photometric measurements of ~2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of a complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K band from Hawk-I, and 3.6 - 4.5 {\mu}m from Spitzer. Full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. The multiwavelength catalogues are publicly available and are ready to be used for scientific purposes. Our procedures allows for the detection of outshined objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in a companion Paper II. [abridged]

We present the first public release of photometric redshifts, galaxy rest-frame properties and associated magnification values in the cluster and parallel pointings of the first two Frontier Fields, Abell-2744 and MACS-J0416. We exploit a multi-wavelength catalogue ranging from HST to ground-based K and Spitzer IRAC which is specifically designed to enable detection and measurement of accurate fluxes in crowded cluster regions. The multi-band information is used to derive photometric redshifts and physical properties of sources detected either in the H-band image alone or from a stack of four WFC3 bands. To minimize systematics median photometric redshifts are assembled from six different approaches to photo-z estimates. Their reliability is assessed through a comparison with available spectroscopic samples. State of the art lensing models are used to derive magnification values on an object-by-object basis by taking into account sources positions and redshifts. We show that photometric redshifts reach a remarkable ~3-5% accuracy. After accounting for magnification the H band number counts are found in agreement at bright magnitudes with number counts from the CANDELS fields, while extending the presently available samples to galaxies intrinsically as faint as H160~32-33 thanks to strong gravitational lensing. The Frontier Fields allow to probe the galaxy stellar mass distribution at 0.5-1.5 dex lower masses, depending on magnification, with respect to extragalactic wide fields, including sources at Mstar~ 10^7-10^8 Msun at z>5. Similarly, they allow the detection of objects with intrinsic SFRs>1dex lower than in the CANDELS fields reaching 0.1-1 Msun/yr at z~6-10. [abridged]

Star-forming galaxies are considered to be the leading candidate sources that dominate the cosmic reionization at z>7, and the search for analogs at moderate redshift showing Lyman continuum (LyC) leakage is currently a active line of research. We have observed a star-forming galaxy at z=3.2 with Hubble/WFC3 in the F336W filter, corresponding to the 730-890A rest-frame, and detect LyC emission. This galaxy is very compact and also has large Oxygen ratio [OIII]5007/[OII]3727 (>=10). No nuclear activity is revealed from optical/near-infrared spectroscopy and deep multi-band photometry (including the 6Ms X-ray, Chandra). The measured escape fraction of ionizing radiation spans the range 50-100\%, depending on the IGM attenuation. The LyC emission is detected at S/N=10 with m(F336W)=27.57+/-0.11 and it is spatially unresolved, with effective radius R_e<200pc. Predictions from photoionization and radiative transfer models are in line with the properties reported here, indicating that stellar winds and supernova explosions in a nucleated star-forming region can blow cavities generating density-bounded conditions compatible with optically thin media. Irrespective to the nature of the ionizing radiation, spectral signatures of these sources over the entire electromagnetic spectrum are of central importance for their identification during the epoch of reionization, when the LyC is unobservable. Intriguingly, the Spitzer/IRAC photometric signature of intense rest-frame optical emissions ([OIII]+Hbeta) observed recently at z~7.5-8.5 is similar to what is observed in this galaxy. Only the James Webb Space Telescope will measure optical line ratios at z>7 allowing a direct comparison with lower redshift LyC emitters, as reported here.

The active galactic nuclei X-ray luminosity function traces actively accreting supermassive black holes and is essential for the study of the properties of the active galactic nuclei (AGN) population, black hole evolution, and galaxy-black hole coevolution. Up to now, the AGN luminosity function has been estimated several times in soft (0.5-2 keV) and hard X-rays (2-10 keV). AGN selection in these energy ranges often suffers from identification and redshift incompleteness and, at the same time, photoelectric absorption can obscure a significant amount of the X-ray radiation. We estimate the evolution of the luminosity function in the 5-10 keV band, where we effectively avoid the absorbed part of the spectrum, rendering absorption corrections unnecessary up to NH=10^23 cm^-2. Our dataset is a compilation of six wide, and deep fields: MAXI, HBSS, XMM-COSMOS, Lockman Hole, XMM-CDFS, AEGIS-XD, Chandra-COSMOS, and Chandra-CDFS. This extensive sample of ~1110 AGN (0.01<z<4.0, 41<log L_x<46) is 98% redshift complete with 68% spectroscopic redshifts. We use Bayesian analysis to select the best parametric model from simple pure luminosity and pure density evolution to more complicated luminosity and density evolution and luminosity-dependent density evolution. We estimate the model parameters that describe best our dataset separately for each survey and for the combined sample. We show that, according to Bayesian model selection, the preferred model for our dataset is the luminosity-dependent density evolution (LDDE). Our estimation of the AGN luminosity function does not require any assumption on the AGN absorption and is in good agreement with previous works in the 2-10 keV energy band based on X-ray hardness ratios to model the absorption in AGN up to redshift three. Our sample does not show evidence of a rapid decline of the AGN luminosity function up to redshift four. [abridged]

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.

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 \mu m}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

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.

We present a new method to search for candidate z~>2 Herschel 500{\mu}m sources in the GOODS-North field, using a S500{\mu}m/S24{\mu}m "color deconfusion" technique. Potential high-z sources are selected against low-redshift ones from their large 500{\mu}m to 24{\mu}m flux density ratios. By effectively reducing the contribution from low-redshift populations to the observed 500{\mu}m emission, we are able to identify counterparts to high-z 500{\mu}m sources whose 24{\mu}m fluxes are relatively faint. The recovery of known z~4 starbursts confirms the efficiency of this approach in selecting high-z Herschel sources. The resulting sample consists of 34 dusty star-forming galaxies at z~>2. The inferred infrared luminosities are in the range 1.5x10^12-1.8x10^13 Lsun, corresponding to dust-obscured star formation rates (SFRs) of ~260-3100 Msun/yr for a Salpeter IMF. Comparison with previous SCUBA 850{\mu}m-selected galaxy samples shows that our method is more efficient at selecting high-z dusty galaxies with a median redshift of z=3.07+/-0.83 and 10 of the sources at z~>4. We find that at a fixed luminosity, the dust temperature is ~5K cooler than that expected from the Td-LIR relation at z<1, though different temperature selection effects should be taken into account. The radio-detected subsample (excluding three strong AGN) follows the far-infrared/radio correlation at lower redshifts, and no evolution with redshift is observed out to z~5, suggesting that the far-infrared emission is star formation dominated. The contribution of the high-z Herschel 500{\mu}m sources to the cosmic SFR density is comparable to that of SMG populations at z~2.5 and at least 40% of the extinction-corrected UV samples at z~4 (abridged).

We release the next installment of the Stripe 82 X-ray survey point-source catalog, which currently covers 31.3 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 Legacy field. In total, 6181 unique X-ray sources are significantly detected with {\it XMM-Newton} ($>5\sigma$) and {\it Chandra} ($>4.5\sigma$). This catalog release includes data from {\it XMM-Newton} cycle AO 13, which approximately doubled the Stripe 82X survey area. The flux limits of the Stripe 82X survey are $8.7\times10^{-16}$ erg s$^{-1}$ cm$^{-2}$, $4.7\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, and $2.1\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-10 keV), and full bands (0.5-10 keV), respectively, with approximate half-area survey flux limits of $5.4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, $2.9\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and $1.7\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. We matched the X-ray source lists to available multi-wavelength catalogs, including updated matches to the previous release of the Stripe 82X survey; 88\% of the sample is matched to a multi-wavelength counterpart. Due to the wide area of Stripe 82X and rich ancillary multi-wavelength data, including coadded SDSS photometry, mid-infrared {\it WISE} coverage, near-infrared coverage from UKIDSS and VHS, ultraviolet coverage from {\it GALEX}, radio coverage from FIRST, and far-infrared coverage from {\it Herschel}, as well as existing $\sim$30\% optical spectroscopic completeness, we are beginning to uncover rare objects, such as obscured high-luminosity AGN at high-redshift. The Stripe 82X point source catalog is a valuable dataset for constraining how this population grows and evolves, as well as for studying how they interact with the galaxies in which they live.

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

We present T-PHOT, a publicly available software aimed at extracting accurate photometry from low-resolution images of deep extragalactic fields, where the blending of sources can be a serious problem for the accurate and unbiased measurement of fluxes and colours. T-PHOT has been developed within the ASTRODEEP project and it can be considered as the next generation to TFIT, providing significant improvements above it and other similar codes. T-PHOT gathers data from a high-resolution image of a region of the sky, and uses it to obtain priors for the photometric analysis of a lower resolution image of the same field. It can handle different types of datasets as input priors: i) a list of objects that will be used to obtain cutouts from the real high-resolution image; ii) a set of analytical models; iii) a list of unresolved, point-like sources, useful e.g. for far-infrared wavelength domains. We show that T-PHOT yields accurate estimations of fluxes within the intrinsic uncertainties of the method, when systematic errors are taken into account (which can be done thanks to a flagging code given in the output). T-PHOT is many times faster than similar codes like TFIT and CONVPHOT (up to hundreds, depending on the problem and the method adopted), whilst at the same time being more robust and more versatile. This makes it an optimal choice for the analysis of large datasets. In addition we show how the use of different settings and methods significantly enhances the performance. Given its versatility and robustness, T-PHOT can be considered the preferred choice for combined photometric analysis of current and forthcoming extragalactic optical to far-infrared imaging surveys. [abridged]

In active galactic nuclei (AGN)-galaxy co-evolution models, AGN winds and outflows are often invoked to explain why super-massive black holes and galaxies stop growing efficiently at a certain phase of their lives. They are commonly referred to as the leading actors of feedback processes. Evidence of ultra-fast (v>0.05c) outflows in the innermost regions of AGN has been collected in the past decade by sensitive X-ray observations for sizable samples of AGN, mostly at low redshift. Here we present ultra-deep XMM-Newton and Chandra spectral data of an obscured (Nh~2x10^{23} cm^-2), intrinsically luminous (L2-10keV~4x10^{44} erg/s) quasar (named PID352) at z~1.6 (derived from the X-ray spectral analysis) in the Chandra Deep Field-South. The source is characterized by an iron emission and absorption line complex at observed energies of E~2-3 keV. While the emission line is interpreted as being due to neutral iron (consistent with the presence of cold absorption), the absorption feature is due to highly ionized iron transitions (FeXXV, FeXXVI) with an outflowing velocity of 0.14^{+0.02}_{-0.06}c, as derived from photoionization models. The mass outflow rate - ~2 Msun/yr - is similar to the source accretion rate, and the derived mechanical energy rate is ~9.5x10^{44} erg/s, corresponding to 9% of the source bolometric luminosity. PID352 represents one of the few cases where indications of X-ray outflowing gas have been observed at high redshift thus far. This wind is powerful enough to provide feedback on the host galaxy.

We study the properties of satellites in the environment of massive star-forming galaxies at z~1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of 10^11 Msun. At z>1.5, these galaxies typically trace halos of mass >10^13 Msun. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~6*10^9 Msun. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple beta-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to 68% of the central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr.

Highly obscured active galactic nuclei (AGN) are common in nearby galaxies, but are difficult to observe beyond the local Universe, where they are expected to significantly contribute to the black hole accretion rate density. Furthermore, Compton-thick (CT) absorbers (NH>10^24 cm^-2) suppress even the hard X-ray (2-10 keV) AGN nuclear emission, and therefore the column density distribution above 10^24 cm^-2 is largely unknown. We present the identification and multi-wavelength properties of a heavily obscured (NH>~10^25 cm^-2), intrinsically luminous (L(2-10keV)>10^44 erg s^-1) AGN at z=0.353 in the COSMOS field. Several independent indicators, such as the shape of the X-ray spectrum, the decomposition of the spectral energy distribution and X-ray/[NeV] and X-ray/6{\mu}m luminosity ratios, agree on the fact that the nuclear emission must be suppressed by a 10^25 cm^-2 column density. The host galaxy properties show that this highly obscured AGN is hosted in a massive star-forming galaxy, showing a barred morphology, which is known to correlate with the presence of CT absorbers. Finally, asymmetric and blueshifted components in several optical high-ionization emission lines indicate the presence of a galactic outflow, possibly driven by the intense AGN activity (L(Bol)/L(Edd) = 0.3-0.5). Such highly obscured, highly accreting AGN are intrinsically very rare at low redshift, whereas they are expected to be much more common at the peak of the star formation and BH accretion history, at z~2-3. We demonstrate that a fully multi-wavelength approach can recover a sizable sample of such peculiar sources in large and deep surveys such as COSMOS.

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 construct the rest-frame 2--10 keV intrinsic X-ray luminosity function of Active Galactic Nuclei (AGNs) from a combination of X-ray surveys from the all-sky Swift BAT survey to the Chandra Deep Field-South. We use ~3200 AGNs in our analysis, which covers six orders of magnitude in flux. The inclusion of the XMM and Chandra COSMOS data has allowed us to investigate the detailed behavior of the XLF and evolution. In deriving our XLF, we take into account realistic AGN spectrum templates, absorption corrections, and probability density distributions in photometric redshift. We present an analytical expression for the overall behavior of the XLF in terms of the luminosity-dependent density evolution, smoothed two power-law expressions in 11 redshift shells, three-segment power-law expression of the number density evolution in four luminosity classes, and binned XLF. We observe a sudden flattening of the low luminosity end slope of the XLF slope at z>~0.6. Detailed structures of the AGN downsizing have been also revealed, where the number density curves have two clear breaks at all luminosity classes above log LX>43. The two break structure is suggestive of two-phase AGN evolution, consisting of major merger triggering and secular processes.

Heavily obscured, Compton Thick (CT, NH>10^24 cm^-2) AGN may represent an important phase in AGN/galaxy co-evolution and are expected to provide a significant contribution to the cosmic X-ray background (CXB). Through direct X-ray spectra analysis, we selected 39 heavily obscured AGN (NH>3x10^23 cm^-2) in the 2 deg^2 XMM-COSMOS survey. After selecting CT AGN based on the fit of a simple absorbed two power law model to the XMM data, the presence of CT AGN was confirmed in 80% of the sources using deeper Chandra data and more complex models. The final sample of CT AGN comprises 10 sources spanning a large range of redshift and luminosity. We collected the multi-wavelength information available for all these sources, in order to study the distribution of SMBH and host properties, such as BH mass (M_BH), Eddington ratio (\lambda_Edd), stellar mass (M*), specific star formation rate (sSFR) in comparison with a sample of unobscured AGN. We find that highly obscured sources tend to have significantly smaller M_BH and higher \lambda_edd with respect to unobscured ones, while a weaker evolution in M* is observed. The sSFR of highly obscured sources is consistent with the one observed in the main sequence of star forming galaxies, at all redshift. We also present optical spectra, spectral energy distribution (SED) and morphology for the sample of 10 CT AGN: all the available optical spectra are dominated by the stellar component of the host galaxy, and a highly obscured torus component is needed in the SED of the CT sources. Exploiting the high resolution Hubble-ACS images available, we conclude that these highly obscured sources have a significantly larger merger fraction with respect to other X-ray selected samples of AGN. Finally we discuss implications in the context of AGN/galaxy co-evolutionary models, and compare our results with the predictions of CXB synthesis models.

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 present the hard-band ($2-10\,\mathrm{keV}$) X-ray luminosity function (HXLF) of $0.5-2\,\mathrm{keV}$ band selected AGN at high redshift. We have assembled a sample of 141 AGN at $3<z\lesssim5$ from X-ray surveys of different size and depth, in order to sample different regions in the $ L_X - z$ plane. The HXLF is fitted in the range $\mathrm{logL_X\sim43-45}$ with standard analytical evolutionary models through a maximum likelihood procedure. The evolution of the HXLF is well described by a pure density evolution, with the AGN space density declining by a factor of $\sim10$ from $z=3$ to 5. A luminosity-dependent density evolution model which, normally, best represents the HXLF evolution at lower redshift, is also consistent with the data, but a larger sample of low-luminosity ($\mathrm{logL_X}<44$), high-redshift AGN is necessary to constrain this model. We also estimated the intrinsic fraction of AGN obscured by a column density $\mathrm{logN_H}\geq23$ to be $0.54\pm0.05$, with no strong dependence on luminosity. This fraction is higher than the value in the Local Universe, suggesting an evolution of the luminous ($\mathrm{L_X>10^{44}\mathrm{erg\,s^{-1}}}$) obscured AGN fraction from $z=0$ to $z>3$.

We study the evolution of the total star formation (SF) activity, total stellar mass and halo occupation distribution in massive halos by using one of the largest X-ray selected sample of galaxy groups with secure spectroscopic identification in the major blank field surveys (ECDFS, CDFN, COSMOS, AEGIS). We provide an accurate measurement of SFR for the bulk of the star-forming galaxies using very deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations. For undetected IR sources, we provide a well-calibrated SFR from SED fitting. We observe a clear evolution in the level of SF activity in galaxy groups. The total SF activity in the high redshift groups (0.5<z<1.1) is higher with respect to the low redshift (0.15<z<0.5) sample at any mass by 0.8+/-0.12 dex. A milder difference (0.35+/-0.1 dex) is observed between the low redshift bin and the groups at z~0. We show that the level of SF activity is declining more rapidly in the more massive halos than in the more common lower mass halos. We do not observe any evolution in the halo occupation distribution and total stellar mass- halo mass relations in groups. The picture emerging from our findings suggests that the galaxy population in the most massive systems is evolving faster than galaxies in lower mass halos, consistently with a "halo downsizing" scenario.

We investigated the presence of galaxy overdensities around four $z\sim6$ QSOs, namely SDSS J1030+0524 (z = 6.28), SDSS J1148+5251 (z = 6.41), SDSS J1048+4637 (z = 6.20) and SDSS J1411+1217 (z = 5.95), through deep $r$-, $i$- and $z$- band imaging obtained with the wide-field ($\sim23'\times25'$) Large Binocular Camera (LBC) at the Large Binocular Telescope (LBT). We adopted color-color selections within the $i-z$ vs $r-z$ plane to identify samples of $i$-band dropouts at the QSO redshift and measure their relative abundance and spatial distribution in the four LBC fields, each covering $\sim8\times8$ physical Mpc at $z\sim6$. The same selection criteria were then applied to $z$-band selected sources in the $\sim$1 deg$^2$ Subaru-XMM Newton Deep Survey to derive the expected number of dropouts over a blank LBC-sized field ($\sim$0.14 deg$^2$). The four observed QSO fields host a number of candidates larger than what is expected in a blank field. By defining as $i$-band dropouts objects with $z_{AB}<25$, $i-z>1.4$ and undetected in the $r$-band, we found 16, 10, 9, 12 dropouts in SDSS J1030+0524, SDSS J1148+5251, SDSS J1048+4637, and SDSS J1411+1217, respectively, whereas only 4.3 such objects are expected over a 0.14 deg$^2$ blank field. This corresponds to overdensity significances of 3.3, 1.9, 1.7, 2.5$\sigma$, respectively. By considering the total number of dropouts in the four LBC fields and comparing it with what is expected in four blank fields of 0.14 deg$^2$ each, we find that high-z QSOs reside in overdense environments at the $3.7\sigma$ level. This is the first direct and unambiguous measurement of the large scale structures around $z\sim6$ QSOs. [shortened]

A spatial clustering signal has been established in Spitzer/IRAC measurements of the unresolved Cosmic near-Infrared Background (CIB) out to large angular scales, ~1 deg. This CIB signal, while significantly exceeding the contribution from the remaining known galaxies, was further found to be coherent at a highly statistically significant level with the unresolved soft Cosmic X-ray Background (CXB). This measurement probes the unresolved CXB to very faint source levels using deep near-IR source subtraction. We study contributions from extragalactic populations at low to intermediate redshifts to the measured positive cross-power signal of the CIB fluctuations with the CXB. We model the X-ray emission from AGN, normal galaxies and hot gas residing in virialized structures, calculating their CXB contribution including their spatial coherence with all infrared emitting counterparts. We use a halo model framework to calculate the auto and cross-power spectra of the unresolved fluctuations based on the latest constraints of the halo occupation distribution and the biasing of AGN, galaxies and diffuse emission. At small angular scales (<1'), the 4.5mic vs 0.5-2 keV coherence can be explained by shot noise from galaxies and AGN. However, at large angular scales (~10 arcmin) we find that the net contribution from the modeled populations is only able to account for ~3% of the measured CIBxCXB cross-power. The discrepancy suggests that the CIBxCXB signal originates from the same unknown source population producing the CIB clustering signal out to ~1 deg.

We took advantage of the observations carried out by XMM in the COSMOS field during 3.5 years, to study the long term variability of a large sample of AGN (638 sources), in a wide range of redshift (0.1<z<3.5) and X-ray luminosity ($10^{41}<$L(2-10)$<10^{45.5}$). Both a simple statistical method to asses the significance of variability, and the Normalized Excess Variance ($\sigma^{2}_{rms}$) parameter, where used to obtain a quantitative measurement of the variability. Variability is found to be prevalent in most AGN, whenever we have good statistic to measure it, and no significant differences between type-1 and type-2 AGN were found. A flat (slope -0.23+/-0.03) anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity is found, when significantly variable sources are considered all together. When divided in three redshift bins, the anti-correlation becomes stronger and evolving with z, with higher redshift AGN being more variable. We prove however that this effect is due to the pre-selection of variable sources: considering all the sources with available $\sigma^{2}_{rms}$ measurement, the evolution in redshift disappears. For the first time we were also able to study the long term X-ray variability as a function of $M_{\rm BH}$ and Eddington ratio, for a large sample of AGN spanning a wide range of redshift. An anti-correlation between $\sigma^{2}_{rms}$ and $M_{\rm BH}$ is found, with the same slope of the anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity, suggesting that the latter can be a byproduct of the former one. No clear correlation is found between $\sigma^{2}_{rms}$ and the Eddington ratio in our sample. Finally, no correlation is found between the X-ray $\sigma^{2}_{rms}$ and the optical variability.

We study the incidence of nuclear obscuration on a complete sample of 1310 AGN selected on the basis of their rest-frame 2-10 keV X-ray flux from the XMM-COSMOS survey, in the redshift range 0.3<z<3.5. We classify the AGN as obscured or un-obscured on the basis of either the optical spectral properties and the overall SED or the shape of the X-ray spectrum. The two classifications agree in about 70% of the objects, and the remaining 30% can be further subdivided into two distinct classes: at low luminosities X-ray un-obscured AGN do not always show signs of broad lines or blue/UV continuum emission in their optical spectra, most likely due to galaxy dilution effects; at high luminosities broad line AGN may have absorbed X-ray spectra, which hints at an increased incidence of small-scale (sub-parsec) dust-free obscuration. We confirm that the fraction of obscured AGN is a decreasing function of the intrinsic X-ray luminosity, while the incidence of absorption shows significant evolution only for the most luminous AGN, which appear to be more commonly obscured at higher redshift. We find no significant difference between the mean stellar masses and star formation rates of obscured and un-obscured AGN hosts. We conclude that the physical state of the medium responsible for obscuration in AGN is complex, and mainly determined by the radiation environment (nuclear luminosity) in a small region enclosed within the gravitational sphere of influence of the central black hole, but is largely insensitive to the wider scale galactic conditions.

We investigate the evolution of the star formation rate (SFR)-density relation in the Extended Chandra Deep Field South (ECDFS) and the Great Observatories Origin Deep Survey (GOODS) fields up to z~1.6. In addition to the "traditional method", in which the environment is defined according to a statistical measurement of the local galaxy density, we use a "dynamical" approach, where galaxies are classified according to three different environment regimes: group, "filament-like", and field. Both methods show no evidence of a SFR-density reversal. Moreover, group galaxies show a mean SFR lower than other environments up to z~1, while at earlier epochs group and field galaxies exhibit consistent levels of star formation (SF) activity. We find that processes related to a massive dark matter halo must be dominant in the suppression of the SF below z~1, with respect to purely density-related processes. We confirm this finding by studying the distribution of galaxies in different environments with respect to the so-called Main Sequence (MS) of star-forming galaxies. Galaxies in both group and "filament-like" environments preferentially lie below the MS up to z~1, with group galaxies exhibiting lower levels of star-forming activity at a given mass. At z>1, the star-forming galaxies in groups reside on the MS. Groups exhibit the highest fraction of quiescent galaxies up to z~1, after which group, "filament-like", and field environments have a similar mix of galaxy types. We conclude that groups are the most efficient locus for star-formation quenching. Thus, a fundamental difference exists between bound and unbound objects, or between dark matter haloes of different masses.

We present a detailed spectral analysis of new XMM-Newton data of the source CXOC J100043.1+020637, also known as CID-42, detected in the COSMOS survey at z = 0.359. Previous works suggested that CID-42 is a candidate recoiling supermassive black holes showing also an inverted P-Cygni profile in the X- ray spectra at ~6 keV (rest) with an iron emission line plus a redshifted absorption line (detected at 3sigma in previous XMM-Newton and Chandra observations). Detailed analysis of the absorption line suggested the presence of ionized material inflowing into the black hole at high velocity. In the new long XMM-Newton observation, while the overall spectral shape remains constant, the continuum 2-10 keV flux decreased of ~20% with respect to previous observation and the absorption line is undetected. The upper limit on the intensity of the absorption line is EW<162 keV. Extensive Monte Carlo simulations show that the non detection of the line is solely due to variation in the properties of the inflowing material, in agreement with the transient nature of these features, and that the intensity of the line is lower than the previously measured with a probability of 98.8%. In the scenario of CID-42 as recoiling SMBH, the absorption line can be interpreted as due to inflow of gas with variable density and located in the proximity of the SMBH and recoiling with it. New monitoring observations will be requested to further characterize this line.

We have analyzed the {\it XMM-Newton} and {\it Chandra} data overlapping $\sim$16.5 deg$^2$ of Sloan Digital Sky Survey Stripe 82, including $\sim$4.6 deg$^2$ of proprietary {\it XMM-Newton} data that we present here. In total, 3362 unique X-ray sources are detected at high significance. We derive the {\it XMM-Newton} number counts and compare them with our previously reported {\it Chandra} Log$N$-Log$S$ relations and other X-ray surveys. The Stripe 82 X-ray source lists have been matched to multi-wavelength catalogs using a maximum likelihood estimator algorithm. We discovered the highest redshift ($z=5.86$) quasar yet identified in an X-ray survey. We find 2.5 times more high luminosity (L$_x \geq 10^{45}$ erg s$^{-1}$) AGN than the smaller area {\it Chandra} and {\it XMM-Newton} survey of COSMOS and 1.3 times as many identified by XBo\"otes. Comparing the high luminosity AGN we have identified with those predicted by population synthesis models, our results suggest that this AGN population is a more important component of cosmic black hole growth than previously appreciated. Approximately a third of the X-ray sources not detected in the optical are identified in the infrared, making them candidates for the elusive population of obscured high luminosity AGN in the early universe.

In the local Universe, galaxy properties show a strong dependence on environment. In cluster cores, early type galaxies dominate, whereas star-forming galaxies are more and more common in the outskirts. At higher redshifts and in somewhat less dense environments (e.g. galaxy groups), the situation is less clear. One open issue is that of whether and how the star formation rate (SFR) of galaxies in groups depends on the distance from the centre of mass. To shed light on this topic, we have built a sample of X-ray selected galaxy groups at 0<z<1.6 in various blank fields (ECDFS, COSMOS, GOODS). We use a sample of spectroscopically confirmed group members with stellar mass M >10^10.3 M_sun in order to have a high spectroscopic completeness. As we use only spectroscopic redshifts, our results are not affected by uncertainties due to projection effects. We use several SFR indicators to link the star formation (SF) activity to the galaxy environment. Taking advantage of the extremely deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations, we have an accurate, broad-band measure of the SFR for the bulk of the star-forming galaxies. We use multi-wavelength SED fitting techniques to estimate the stellar masses of all objects and the SFR of the MIPS and PACS undetected galaxies. We analyse the dependence of the SF activity, stellar mass and specific SFR on the group-centric distance, up to z~1.6, for the first time. We do not find any correlation between the mean SFR and group-centric distance at any redshift. We do not observe any strong mass segregation either, in agreement with predictions from simulations. Our results suggest that either groups have a much smaller spread in accretion times with respect to the clusters and that the relaxation time is longer than the group crossing time.

The backbone of the large-scale structure of the Universe is determined by processes on a cosmological scale and by the gravitational interaction of the dominant dark matter. However, the mobile baryon population shapes the appearance of these structures. Theory predicts that most of the baryons reside in vast unvirialized filamentary structures that connect galaxy groups and clusters, but the observational evidence is currently lacking. Because the majority of the baryons are supposed to exist in a large-scale, hot and dilute gaseous phase, X-rays provide the ideal tool to progress our understanding. Observations with the Athena+ X-ray Integral Field Unit will reveal the location, chemical composition, physical state and dynamics of the active population of baryons.

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.

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.

In order to understand the nature of the sources producing the recently uncovered cosmic infrared background (CIB) fluctuations, we study cross-correlations between the fluctuations in the source-subtracted CIB from Spitzer/IRAC data and the unresolved cosmic X-ray background from deep Chandra observations. Our study uses data from the EGS/AEGIS field, where both data sets cover an 8'x45' region of the sky. Our measurement is the cross-power spectrum between the IR and X-ray data. The cross-power signal between the IRAC maps at 3.6 um and 4.5 um and the Chandra [0.5-2] keV data has been detected, at angular scales 20'', with an overall significance of 3.8 sigma and 5.6 sigma, respectively. At the same time we find no evidence of significant cross-correlations at the harder Chandra bands. The cross-correlation signal is produced by individual IR sources with 3.6 um and 4.5 um magnitudes m_AB 25-26 and [0.5-2] keV X-ray fluxes 7e-17 cgs.We determine that at least 15%-25% of the large scale power of the CIB fluctuations is correlated with the spatial power spectrum of the X-ray fluctuations. If this correlation is attributed to emission from accretion processes at both IR and X-ray wavelengths, this implies a much higher fraction of accreting black holes than among the known populations. We discuss the various possible origins for the cross-power signal and show that neither local foregrounds nor the known remaining normal galaxies and active galactic nuclei can reproduce the measurements. These observational results are an important new constraint on theoretical modeling of the near-IR CIB fluctuations.

X-ray spectroscopy of active galactic nuclei (AGN) offers the opportunity to directly probe the inner regions of the accretion disk. We present the results of our analysis of average AGN XMM-Newton X-ray spectra in the Chandra Deep Field South observation (XMM CDFS). We computed the average spectrum of a sample of 54 AGN with spectroscopic redshifts and signal-to-noise ratio S/N > 15 in the 2-12 keV rest-frame band in at least one EPIC camera. We have taken the effects of combining spectra from sources at different redshifts and both EPIC cameras into account, as well as their spectral resolution; we checked our results using thorough simulations. We explored the Fe line components of distant AGN focusing on the narrow core which arises from regions far from the central engine and on the putative relativistic component (from the accretion disk). The average spectrum shows a highly significant Fe feature. Our model-independent estimates of the equivalent width (EW) suggest a higher EW in a broader range. The line, modelled as an unresolved Gaussian, is significant at 6.8 sigma and has an EW=95 eV (full sample). The current data can be fitted equally well adding a relativistic profile to the narrow component (in the full sample, EW=140 eV and 67 eV respectively for the relativistic and narrow lines). Thanks to the high quality of the XMM CDFS spectra and to the detailed modelling of the continuum and instrumental effects, we have shown that the most distant AGN exhibit a highly significant Fe emission feature. It can be modelled both with narrow and broad lines. We found tantalising evidence for reflection by material both very close and far away from the central engine. The EW of both features are similar to those observed in individual nearby AGN, hence they must be a widespread characteristic of AGN, since otherwise the average values would be smaller than observed.

(abridged) The XMM-Newton survey in the Chandra Deep Field South (XMM-CDFS) aims at detecting and studying the spectral properties of a significant number of obscured and Compton-thick AGN. The large effective area of XMMin the 2--10 and 5--10 keV bands, coupled with a 3.45 Ms nominal exposure time, allows us to build clean samples in both bands, and makes the XMM-CDFS the deepest XMM survey currently published in the 5--10 keV band. The large multi-wavelength and spectroscopic coverage of the CDFS area allows for an immediate and abundant scientific return. We present the data reduction of the XMM-CDFS observations, the method for source detection in the 2--10 and 5--10keV bands, and the resulting catalogues. A number of 339 and 137 sources are listed in the above bands with flux limits of 6.6e-16 and 9.5e-16 erg/s/cm^2, respectively. The flux limits at 50% of the maximum sky coverage are 1.8e-15 and 4.0e-15 erg/s/cm^2, respectively. The catalogues have been cross-correlated with the Chandra ones: 315 and 130 identifications have been found with a likelihood-ratio method, respectively. A number of 15 new sources, previously undetected by Chandra, is found; 5 of them lie in the 4 Ms area. Redshifts, either spectroscopic or photometric, are available for ~92% of the sources. The number counts in both bands are presented and compared to other works. The survey coverage has been calculated with the help of two extensive sets of simulations, one set per band. The simulations have been produced with a newly-developed simulator, written with the aim of the most careful reproduction of the background spatial properties. We present a detailed decomposition of the XMM background into its components: cosmic, particle, and residual soft protons.

We present the first results of a wide area X-ray survey within the Sloan Digital Sky Survey (SDSS) Stripe 82, a 300 deg$^2$ region of the sky with a substantial investment in multi-wavelength coverage. We analyzed archival {\it Chandra} observations that cover 7.5 deg$^2$ within Stripe 82 ("Stripe 82 ACX"), reaching 4.5$\sigma$ flux limits of 7.9$\times10^{-16}$, 3.4$\times10^{-15}$ and 1.8$\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-7 keV) and full (0.5-7 keV) bands, to find 774, 239 and 1118 X-ray sources, respectively. Three hundred twenty-one sources are detected only in the full band and 9 sources are detected solely in the soft band. Utilizing data products from the {\it Chandra} Source Catalog, we construct independent Log$N$-Log$S$ relationships, detailing the number density of X-ray sources as a function of flux, which show general agreement with previous {\it Chandra} surveys. We compare the luminosity distribution of Stripe 82 ACX with the smaller, deeper CDF-S + E-CDFS surveys and with {\it Chandra}-COSMOS, illustrating the benefit of wide-area surveys in locating high luminosity AGN. We also investigate the differences and similarities of X-ray and optical selection to uncover obscured AGN in the local Universe. Finally, we estimate the population of AGN we expect to find with increased coverage of 100 deg$^2$ or 300 deg$^2$, which will provide unprecedented insight into the high redshift, high luminosity regime of black hole growth currently under-represented in X-ray surveys.

The Chandra Deep Field is the region of the sky with the highest concentration of X-ray data available: 4Ms of Chandra and 3Ms of XMM data, allowing excellent quality spectra to be extracted even for faint sources. We take advantage of this in order to compile a sample of heavily obscured Active Galactic Nuclei (AGN) using X-ray spectroscopy. We select our sample among the 176 brightest XMM sources, searching for either a) flat X-ray spectra (Photon index<1.4 at the 90% confidence level) suggestive of a reflection dominated continuum or b) an absorption turn-over suggestive of a column density higher than ~10^{24} cm-2. We find a sample of nine candidate heavily obscured sources satisfying the above criteria. Four of these show statistically significant FeKalpha lines with large equivalent widths (three out of four have EW consistent with 1 keV) suggesting that these are the most secure Compton-thick AGN candidates. Two of these sources are transmission dominated while the other two are most probably reflection dominated Compton-thick AGN. Although this sample of four sources is by no means statistically complete, it represents the best example of Compton-thick sources found at moderate-to-high redshift with three sources at z=1.2-1.5 and one source at z=3.7. Using Spitzer and Herschel observations, we estimate with good accuracy the X-ray to mid-IR (12 micron) luminosity ratio of our sources. These are well below the average AGN relation, independently suggesting that these four sources are heavily obscured.

We present the results of a search for extended X-ray sources and their corresponding galaxy groups from 800-ks Chandra coverage of the All-wavelength Extended Groth Strip International Survey (AEGIS). This yields one of the largest X-ray selected galaxy group catalogs from a blind survey to date. The red-sequence technique and spectroscopic redshifts allow us to identify 100$%$ of reliable sources, leading to a catalog of 52 galaxy groups. The groups span the redshift range $z\sim0.066-1.544$ and virial mass range $M_{200}\sim1.34\times 10^{13}-1.33\times 10^{14}M_\odot$. For the 49 extended sources which lie within DEEP2 and DEEP3 Galaxy Redshift Survey coverage, we identify spectroscopic counterparts and determine velocity dispersions. We select member galaxies by applying different cuts along the line of sight or in projected spatial coordinates. A constant cut along the line of sight can cause a large scatter in scaling relations in low-mass or high-mass systems depending on the size of cut. A velocity dispersion based virial radius can more overestimate velocity dispersion in comparison to X-ray based virial radius for low mass systems. There is no significant difference between these two radial cuts for more massive systems. Independent of radial cut, overestimation of velocity dispersion can be created in case of existence of significant substructure and also compactness in X-ray emission which mostly occur in low mass systems. We also present a comparison between X-ray galaxy groups and optical galaxy groups detected using the Voronoi-Delaunay method (VDM) for DEEP2 data in this field.

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.

(Abridged) We report the discovery of an X-ray group of galaxies located at a high redshift of z=1.61 in the Chandra Deep Field South. The group is first identified as an extended X-ray source. We use a wealth of deep multi-wavelength data to identify the optical counterpart -- our red sequence finder detects a significant over-density of galaxies at z~1.6 and the brightest group galaxy is spectroscopically confirmed at z=1.61. We measure an X-ray luminosity of L_{0.1-2.4 keV}= 1.8\pm0.6 \times 10^{43} erg/s, which then translates into a group mass of 3.2\pm0.8 \times 10^{13} M_sun. This is the lowest mass group ever confirmed at z>1.5. The deep optical-nearIR images from CANDELS reveal that the group exhibits a surprisingly prominent red sequence. A detailed analysis of the spectral energy distributions of the group member candidates confirms that most of them are indeed passive galaxies. Furthermore, their structural parameters measured from the near-IR CANDELS images show that they are morphologically early-type. The newly identified group at z=1.61 is dominated by quiescent early-type galaxies and the group appears similar to those in the local Universe. One possible difference is the high fraction of AGN (38^{+23}_{-20}%), which might indicate a role for AGN in quenching. But, a statistical sample of high-z groups is needed to draw a general picture of groups at this redshift. Such a sample will hopefully be available in near future surveys.

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 "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

We present results on a search of heavily obscured active galaxies z>1.7 using the rest-frame 9-20 keV excess for X-ray sources detected in the deep XMM-CDFS survey. Out of 176 sources selected with the conservative detection criteria (>8 sigma) in the first source catalogue of Ranalli et al., 46 objects lie in the redshift range of interest with the median redshift z~2.5. Their typical rest-frame 10-20 keV luminosity is 1e+44 erg/s, as observed. Among optically faint objects that lack spectroscopic redshift, four were found to be strongly absorbed X-ray sources, and the enhanced Fe K emission or absorption features in their X-ray spectra were used to obtain X-ray spectroscopic redshifts. Using the X-ray colour-colour diagram based on the rest-frame 3-5 keV, 5-9 keV, and 9-20 keV bands, seven objects were selected for their 9-20 keV excess and were found to be strongly absorbed X-ray sources with column density of nH > 0.6e+24 cm-2, including two possible Compton thick sources. While they are emitting at quasar luminosity, ~3/4 of the sample objects are found to be absorbed by nH > 1e+22 cm-2. A comparison with local AGN at the matched luminosity suggests an increasing trend of the absorbed source fraction for high-luminosity AGN towards high redshifts.

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

We present new IRAM PdBI 1.3mm continuum observations at ~1.5" resolution of 28 SMGs previously discovered with the 870um bolometer LABOCA at APEX within the central 0.7deg2 of the COSMOS field. 19 out of the 28 LABOCA sources were detected with the PdBI at a >~3sigma level of ~1.4mJy/b. A combined analysis of this new sample with existing interferometrically identified SMGs in the COSMOS field yields the following results: 1) >~15%, and possibly up to ~40% of single-dish detected SMGs consist of multiple sources, 2) statistical identifications of multi-wavelength counterparts to the single-dish SMGs yield that only ~50% of these single-dish SMGs have real radio or IR counterparts, 3) ~18% of interferometric SMGs have only radio or even no multi-wavelength counterpart at all, and 4) ~50-70% of z>~3 SMGs have no radio counterparts down to an rms of 7-12uJy at 1.4GHz. Using the exact interferometric positions to identify proper multi-wavelength counterparts allows us to determine accurate photometric redshifts for these sources. The redshift distributions of the combined and the individual 1.1mm and 870um selected samples have a higher mean and broader width than the redshift distributions derived in previous studies. Our sample supports the previous tentative trend that on average brighter and/or mm-selected SMGs are located at higher redshifts. There is a tentative offset between the mean redshift for the 1.1mm (<z>=3.1+/-0.4) and 870um (<z>=2.6+/-0.4) selected samples, with the 1.1mm sources lying on average at higher redshifts. Based on our nearly complete sample of AzTEC 1.1mm SMGs within a uniform 0.15deg2 area we infer a higher surface density of z>~4 SMGs than predicted by current cosmological models. In summary, our findings imply that (sub-)millimeter interferometric identifications are crucial to build statistically complete and unbiased samples of SMGs.

At present, 30-40 per cent of the baryons in the local Universe is still undetected. According to theoretical predictions, this gas should reside in filaments filling the large-scale structure (LSS) in the form of a Warm-Hot Intergalactic Medium (WHIM), at temperatures of 10^5 - 10^7 K, thus emitting in the soft X-ray energies via free-free interaction and line emission from heavy elements. In this work we characterize the properties of the X-ray emission of the WHIM, and the LSS in general, focusing on the influence of different physical mechanisms, namely galactic winds (GWs), black-hole feedback and star-formation, and providing estimates of possible observational constraints. To this purpose we use a set of cosmological hydrodynamical simulations that include a self-consistent treatment of star-formation and chemical enrichment of the intergalactic medium, that allows us to follow the evolution of different metal species. We construct a set of simulated light-cones to make predictions of the emission in the 0.3-10 keV energy range. We obtain that GWs increase by a factor of 2 the emission of both galaxy clusters and WHIM. The amount of oxygen at average temperature and, consequently, the amount of expected bright Ovii and Oviii lines is increased by a factor of 3 due to GWs and by 20 per cent when assuming a top-heavy IMF. We compare our results with current observational constraints and find that the emission from faint groups and WHIM should account from half to all of the unresolved X-ray background in the 1-2 keV band.

X-ray surveys contain sizable numbers of star forming galaxies, beyond the AGN which usually make the majority of detections. Many methods to separate the two populations are used in the literature, based on X-ray and multiwavelength properties. We aim at a detailed test of the classification schemes and to study the X-ray properties of the resulting samples. We build on a sample of galaxies selected at 1.4 GHz in the VLA-COSMOS survey, classified by Smolcic et al. (2008) according to their optical colours and observed with Chandra. A similarly selected control sample of AGN is also used for comparison. We review some X-ray based classification criteria and check how they affect the sample composition. The efficiency of the classification scheme devised by Smolcic et al. (2008) is such that ~30% of composite/misclassified objects are expected because of the higher X-ray brightness of AGN with respect to galaxies. The latter fraction is actually 50% in the X-ray detected sources, while it is expected to be much lower among X-ray undetected sources. Indeed, the analysis of the stacked spectrum of undetected sources shows, consistently, strongly different properties between the AGN and galaxy samples. X-ray based selection criteria are then used to refine both samples. The radio/X-ray luminosity correlation for star forming galaxies is found to hold with the same X-ray/radio ratio valid for nearby galaxies. Some evolution of the ratio may be possible for sources at high redshift or high luminosity, tough it is likely explained by a bias arising from the radio selection. Finally, we discuss the X-ray number counts of star forming galaxies from the VLA- and C-COSMOS surveys according to different selection criteria, and compare them to the similar determination from the Chandra Deep Fields. The classification scheme proposed here may find application in future works and surveys.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra program that has imaged the central 0.9 deg^2 of the COSMOS field down to limiting depths of 1.9 10^-16 erg cm^-2 s-1 in the 0.5-2 keV band, 7.3 10^-16 erg cm^-2 s^-1 in the 2-10 keV band, and 5.7 10^-16 erg cm^-2 s-1 in the 0.5-10 keV band. In this paper we report the i, K and 3.6micron identifications of the 1761 X-ray point sources. We use the likelihood ratio technique to derive the association of optical/infrared counterparts for 97% of the X-ray sources. For most of the remaining 3%, the presence of multiple counterparts or the faintness of the possible counterpart prevented a unique association. For only 10 X-ray sources we were not able to associate a counterpart, mostly due to the presence of a very bright field source close by. Only 2 sources are truly empty fields. Making use of the large number of X-ray sources, we update the "classic locus" of AGN and define a new locus containing 90% of the AGN in the survey with full band luminosity >10^42 erg/s. We present the linear fit between the total i band magnitude and the X-ray flux in the soft and hard band, drawn over 2 orders of magnitude in X-ray flux, obtained using the combined C-COSMOS and XMM-COSMOS samples. We focus on the X-ray to optical flux ratio (X/O) and we test its known correlation with redshift and luminosity, and a recently introduced anti-correlation with the concentration index (C). We find a strong anti-correlation (though the dispersion is of the order of 0.5 dex) between C and X/O, computed in the hard band, and that 90% of the obscured AGN in the sample with morphological information live in galaxies with regular morphology (bulgy and disky/spiral), suggesting that secular processes govern a significant fraction of the BH growth at X-ray luminosities of 10^43- 10^44.5 erg/s.

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.

X-ray selected galaxy group samples are usually generated by searching for extended X- ray sources that reflect the thermal radiation of the intragroup medium. On the other hand, large radio galaxies that regularly occupy galaxy groups also emit in the X-ray window, and their contribution to X-ray selected group samples is still not well understood. In order to investigate their relative importance, we have carried out a systematic search for non-thermal extended X-ray sources in the COSMOS field. Based on the morphological coincidence of X-ray and radio extensions, out of 60 radio galaxies, and \sim 300 extended X-ray sources, we find only one candidate where the observed extended X-ray emission arises from non- thermal processes related to radio galaxies. We present a detailed analysis of this source, and its environment. Our results yield that external Inverse Compton emission of the lobes is the dominant process that generates the observed X-ray emission of our extended X-ray candidate, with a minor contribution from the gas of the galaxy group hosting the radio galaxy. Finally, we show that finding only one potential candidate in the COSMOS field (in a redshift range 0 < z < 6 and with radio luminosity between 1025 and 1030 W/Hz) is consistent with expected X-ray-counts arising from synchrotron lobes. This implies that these sources are not a prominent source of contamination in samples of X-ray selected clusters/groups, but they could potentially dominate the z > 1 cluster counts at the bright end (S_X > 7 \cdot 10^-15 erg s^-1 cm^2).

We present a rest-frame spectral stacking analysis of ~1000 X-ray sources detected in the XMM-COSMOS field in order to investigate the iron K line properties of active galaxies beyond redshift z~1. In Type I AGN that have a typical X-ray luminosity of Lx~1.5e44 erg/s and z~1.6, the cold Fe K at 6.4 keV is weak (EW~0.05keV), in agreement with the known trend. In contrast, high-ionization lines of Fe XXV and Fe XXVI are pronounced. These high-ionization Fe K lines appear to have a connection with high accretion rates. While no broad Fe emission is detected in the total spectrum, it might be present, albeit at low significance, when the X-ray luminosity is restricted to the range below 3e44 erg/s, or when an intermediate range of Eddington ratio around 0.1 is selected. In Type II AGN, both cold and high-ionzation lines become weak with increasing X-ray luminosity. However, strong high-ionization Fe K (EW~0.3 keV) is detected in the spectrum of objects at z>2, while no 6.4 keV line is found. It is then found that the primary source of the high-ionization Fe K emission is those objects detected with Spitzer-MIPS at 24 micron. Given their median redshift of z=2.5, their bolometric luminosity is likely to reach 10^13 Lsun and the MIPS-detected emission most likely originates from hot dust heated by embedded AGN, probably accreting at high Eddington ratio. These properties match those of rapidly growing black holes in ultra-luminous infrared galaxies at the interesting epoch (z=2-3) of galaxy formation.

We explore the connection between black hole growth at the center of obscured quasars selected from the XMM-COSMOS survey and the physical properties of their host galaxies. We study a bolometric regime (<Lbol > 8 x 10^45 erg/s) where several theoretical models invoke major galaxy mergers as the main fueling channel for black hole accretion. We confirm that obscured quasars mainly reside in massive galaxies (Mstar>10^10 Msun) and that the fraction of galaxies hosting such powerful quasars monotonically increases with the stellar mass. We stress the limitation of the use of rest-frame color-magnitude diagrams as a diagnostic tool for studying galaxy evolution and inferring the influence that AGN activity can have on such a process. We instead use the correlation between star-formation rate and stellar mass found for star-forming galaxies to discuss the physical properties of the hosts. We find that at z ~1, ~62% of Type-2 QSOs hosts are actively forming stars and that their rates are comparable to those measured for normal star-forming galaxies. The fraction of star-forming hosts increases with redshift: ~71% at z ~2, and 100% at z ~3. We also find that the the evolution from z ~1 to z ~3 of the specific star-formation rate of the Type-2 QSO hosts is in excellent agreement with that measured for star-forming galaxies. From the morphological analysis, we conclude that most of the objects are bulge-dominated galaxies, and that only a few of them exhibit signs of recent mergers or disks. Finally, bulge-dominated galaxies tend to host Type-2 QSOs with low Eddington ratios (lambda<0.1), while disk-dominated or merging galaxies have at their centers BHs accreting at high Eddington ratios (lambda > 0.1).

Close encounters between galaxies are expected to be a viable mechanism, as predicted by numerical simulations, by which accretion onto supermassive black holes can be initiated. To test this scenario, we construct a sample of 562 galaxies (M*>2.5x10^10 M_sun) in kinematic pairs over the redshift range 0.25 < z < 1.05 that are more likely to be interacting than a well-matched control sample of 2726 galaxies not identified as being in a pair, both from the zCOSMOS 20k spectroscopic catalog. Galaxies that harbor an active galactic nucleus (AGN) are identified on the basis of their X-ray emission (L_x>2x10^42 erg s^-1) detected by Chandra. We find a higher fraction of AGN in galaxies in pairs relative to isolated galaxies of similar stellar mass. Our result is primarily due to an enhancement of AGN activity, by a factor of 1.9 (observed) and 2.6 (intrinsic), for galaxies in pairs of projected separation less than 75 kpc and line-of-sight velocity offset less than 500 km s^-1. This study demonstrates that close kinematic pairs are conducive environments for black hole growth either indicating a causal physical connection or an inherent relation, such as, to enhanced star formation. In the Appendix, we describe a method to estimate the intrinsic fractions of galaxies (either in pairs or the field) hosting an AGN with confidence intervals, and an excess fraction in pairs. We estimate that 17.8_{-7.4}^{+8.4}% of all moderate-luminosity AGN activity takes place within galaxies undergoing early stages of interaction that leaves open the question as to what physical processes are responsible for fueling the remaining ~80% that may include late-stage mergers.

With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.

As the mid-IR luminosity represents a good isotropic proxy of the AGN power, a low X-ray to mid-IR luminosity ratio is often claimed to be a reliable indicator for selecting Compton-thick (CT) AGN. We assess the efficiency of this diagnostic by examining the 12mu IRAS AGN sample for which high signal-to-noise XMM observations have been recently become available. We find that the vast majority (10/11) of the AGN that have been classified as CT on the basis the X-ray spectroscopy by Brightman & Nandra present a low Lx/L6 luminosity ratio, i.e. lower than a few percent of the average AGN ratio which is typical of reflection-dominated CT sources. At low Lx/L6 ratios we also find a comparable number of AGN, most of which are heavily absorbed but not CT. This implies that although most Compton-thick AGN present low Lx/L6 ratios, at least in the local, Universe, the opposite is not necessarily true. Next, we extend our analysis to higher redshifts. We perform the same analysis in the CDFS where excellent quality chandra (4 Ms) and xmm (3 Ms) X-ray spectra are available. We derive accurate X-ray luminosities for chandra sources using X-ray spectral fits, as well as 6mu luminosities from SED fits. We find 8 AGN with low Lx/L6 ratios in total, after excluding one source where the 6mu emission primarily comes from star-formation. One of these sources has been already demonstrated to host a CT nucleus, while for another one at a redshift of z=1.22 we argue it is most likely CT on the basis of its combined chandra and xmm spectrum. We find a large number of non CT contaminant with low Lx/L6 ratios. The above suggest that a low Lx/L6 ratio alone cannot ascertain the presence of a CT AGN, albeit the majority of low Lx/L6 AGN are heavily obscured. The two most reliable CT AGN in the high redshift Universe have high Lx/L6 ratios, showing that this method cannot provide complete CT AGN samples.

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.

We present the first results of the spectroscopy of distant, obscured AGN as obtained with the ultra-deep (~3.3 Ms) XMM-Newton survey in the Chandra Deep Field South (CDFS). One of the primary goals of the project is to characterize the X-ray spectral properties of obscured and heavily obscured Compton-thick AGN over the range of redhifts and luminosities that are relevant in terms of their contribution to the X-ray background. The ultra-deep exposure, coupled with the XMM detector's spectral throughput, allowed us to accumulate good quality X-ray spectra for a large number of X-ray sources and, in particular, for heavily obscured AGN at cosmological redshifts. Specifically we present the X-ray spectral properties of two high-redshift - z= 1.53 and z=3.70 - sources. The XMM spectra of both are very hard, with a strong iron Kalpha line at a rest-frame energy of 6.4 keV. A reflection-dominated continuum provides the best description of the X-ray spectrum of the z=1.53 source, while the intrinsic continuum of the z=3.70 AGN is obscured by a large column N_H ~ 10^24 cm-2 of cold gas. Compton-thick absorption, or close to it, is unambiguously detected in both sources. Interestingly, these sources would not be selected as candidate Compton thick AGN by some multiwavelength selection criteria based on the mid-infrared to optical and X-ray to optical flux ratios.

We present a study of the Arches cluster based on XMM-Newton observations performed over the past 8 years. Unexpectedly, we find that the X-ray emission associated with the cluster experienced a marked brightening in March/April 2007. We investigate the origin of both the X-ray continuum emission emanating from the star cluster and the flare. To study the time variability of the total X-ray flux, we stacked the PN and MOS data of observations performed within a time interval of a few days leading to the detection of the flaring episode. We then constructed two spectral datasets, one corresponding to the flare interval (March/April 2007) and another to the normal quiescent state of the source. The X-ray light curve of the Arches cluster shows, with high significance (8.6 sigma), a 70% increase in the X-ray emission in the March/April 2007 timeframe followed by a decline over the following year to the pre-flare level; the short-term duration of the flare is constrained to be longer than four days. The temperature and the line-of-sight column density inferred from the flare spectrum do not differ from those measured in the normal activity state of the cluster, suggesting that the flux enhancement is thermal in origin. We attribute the X-ray variability to in situ stellar activity: early-type stars may be responsible for the flare via wind collisions, whereas late-type stars may contribute by means of magnetic reconnection. These two possibilities are discussed.

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)<11.7 with high-resolution HST/ACS imaging from the COSMOS field. The visual analysis of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted morphologies in the AGN and control samples, i.e. quantifying the signature of recent mergers which might potentially be responsible for fueling/triggering the AGN. We find that (1) the vast majority (>85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

We report on the discovery of 10 additional galaxy clusters detected in the ongoing Swift/BAT all-sky survey. Among the newly BAT-discovered clusters there are: Bullet, Abell 85, Norma, and PKS 0745-19. Norma is the only cluster, among those presented here, which is resolved by BAT. For all the clusters we perform a detailed spectral analysis using XMM-Newton and Swift/BAT data to investigate the presence of a hard (non-thermal) X-ray excess. We find that in most cases the clusters' emission in the 0.3-200keV band can be explained by a multi-temperature thermal model confirming our previous results. For two clusters (Bullet and Abell 3667) we find evidence for the presence of a hard X-ray excess. In the case of the Bullet cluster, our analysis confirms the presence of a non-thermal, power-law like, component with a 20-100 keV flux of 3.4 \times 10-12 erg cm-2 s-1 as detected in previous studies. For Abell 3667 the excess emission can be successfully modeled as a hot component (kT=~13keV). We thus conclude that the hard X-ray emission from galaxy clusters (except the Bullet) has most likely thermal origin.

We present the clustering measurement of hard X-ray selected AGN in the local Universe. We used a sample of 199 sources spectroscopically confirmed detected by Swift-BAT in its 15-55 keV all-sky survey. We measured the real space projected auto-correlation function and detected a signal significant on projected scales lower than 200 Mpc/h. We measured a correlation length of r0=5.56+0.49-0.43 Mpc/h and a slope {\gamma}=1.64-0.08 -0.07. We also measured the auto-correlation function of Type I and Type II AGN and found higher correlation length for Type I AGN. We have a marginal evidence of luminosity dependent clustering of AGN, as we detected a larger correlation length of luminous AGN than that of low luminosity sources. The corresponding typical host DM halo masses of Swift-BAT are log(MDMH) 12-14 h^-1 M/M_sun, depending on the subsample. For the whole sample we measured log(MDMH)\sim 13.15 h-1 M/M_sun which is the typical mass of a galaxy group. We estimated that the local AGN population has a typical lifetime tau_AGN \sim 0.7 Gyr, it is powered by SMBH with mass MBH \sim 1-10x10^8 M_\odot and accreting with very low efficiency, log(epsilon)-2.0. We also conclude that local AGN host galaxies are typically red-massive galaxies with stellar mass of the order 2-80x10^10 h^-1 M_sun. We compared our results with clustering predictions of merger-driven AGN triggering models and found a good agreement.

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the core instrument on the Russian Spektrum-Roentgen-Gamma (SRG) mission which is scheduled for launch in late 2012. eROSITA is fully approved and funded by the German Space Agency DLR and the Max-Planck-Society. The design driving science is the detection of 50 - 100 thousands Clusters of Galaxies up to redshift z ~ 1.3 in order to study the large scale structure in the Universe and test cosmological models, especially Dark Energy. This will be accomplished by an all-sky survey lasting for four years plus a phase of pointed observations. eROSITA consists of seven Wolter-I telescope modules, each equipped with 54 Wolter-I shells having an outer diameter of 360 mm. This would provide and effective area at 1.5 keV of ~ 1500 cm2 and an on axis PSF HEW of 15" which would provide an effective angular resolution of 25"-30". In the focus of each mirror module, a fast frame-store pn-CCD will provide a field of view of 1 deg in diameter for an active FOV of ~ 0.83 deg^2. At the time of writing the instrument development is currently in phase C/D.

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the core instrument on the Russian Spektrum-Roentgen-Gamma (SRG) mission which is scheduled for launch in late 2012. eROSITA is fully approved and funded by the German Space Agency DLR and the Max-Planck-Society.The design driving science is the detection of 50-100 thousands Clusters of Galaxies up to redshift z~1.3 in order to study the large scale structure in the Universe and test cosmological models, especially Dark Energy. This will be accomplished by an all-sky survey lasting for four years plus a phase of pointed observations. eROSITA consists of seven Wolter-I telescope modules, each equipped with 54 Wolter-I shells having an outer diameter of 360 mm. This would provide and effective area at 1.5 keV of ~1500 cm2 and an on axis PSF HEW of 15" which would provide an effective angular resolution of 25"-30". In the focus of each mirror module, a fast frame-store pn-CCD will provide a field of view of 1deg in diameter for an active FOV of ~0.83 deg2. At the time of writing the instrument development is currently in phase C/D.

Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal SZ signal, and compare them to X-ray cluster detections. Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the simulation volume to identify true and spurious detections. Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations with the characteristics of ACT detect clusters with masses M>6-9e13 M_o/h, quite independent of redshift, reach 50% completeness at ~1e14 M_o/h and 100% completeness at ~2e14 M_o/h. Samples are contaminated by only a few per cent spurious detections. This is broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and X-ray maps confirms that the filters introduce at most a negligible bias.

We report the final optical identifications of the medium-depth (~60 ksec), contiguous (2 deg^2) XMM-Newton survey of the COSMOS field. XMM-Newton has detected ~800 X-ray sources down to limiting fluxes of ~5x10^{-16}, ~3x10^{-15}, and ~7x10^{-15} erg/cm2/s in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. The work is complemented by an extensive collection of multi-wavelength data from 24 micron to UV, available from the COSMOS survey, for each of the X-ray sources, including spectroscopic redshifts for ~50% of the sample, and high-quality photometric redshifts for the rest. The XMM and multiwavelength flux limits are well matched: 1760 (98%) of the X-ray sources have optical counterparts, 1711 (~95%) have IRAC counterparts, and 1394 (~78%) have MIPS 24micron detections. Thanks to the redshift completeness (almost 100%) we were able to constrain the high-luminosity tail of the X-ray luminosity function confirming that the peak of the number density of logL_X>44.5 AGN is at z~2. Spectroscopically-identified obscured and unobscured AGN, as well as normal and starforming galaxies, present well-defined optical and infrared properties. We devised a robust method to identify a sample of ~150 high redshift (z>1), obscured AGN candidates for which optical spectroscopy is not available. We were able to determine that the fraction of the obscured AGN population at the highest (L_X>10^{44} erg s^{-1}) X-ray luminosity is ~15-30% when selection effects are taken into account, providing an important observational constraint for X-ray background synthesis. We studied in detail the optical spectrum and the overall spectral energy distribution of a prototypical Type 2 QSO, caught in a stage transitioning from being starburst dominated to AGN dominated, which was possible to isolate only thanks to the combination of X-ray and infrared observations.

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.

Using Chandra observations we have identified a sample of seven off-nuclear X-ray sources, in the redshift range z=0.072-0.283, located within optically bright galaxies in the COSMOS Survey. Using the multi-wavelength coverage available in the COSMOS field, we study the properties of the host galaxies of these ULXs. In detail, we derived their star formation rate from H_alpha measurements and their stellar masses using SED fitting techniques with the aim to compute the probability to have an off-nuclear source based on the host galaxy properties. We divide the host galaxies in different morphological classes using the available ACS/HST imaging. We find that our ULXs candidates are located in regions of the SFR versus M$_star$ plane where one or more off-nuclear detectable sources are expected. From a morphological analysis of the ACS imaging and the use of rest-frame colours, we find that our ULXs are hosted both in late and early type galaxies. Finally, we find that the fraction of galaxies hosting a ULX ranges from ~0.5% to ~0.2% going from L[0.5-2 keV]=3 x 10^39 erg s^-1 to L[0.5-2 keV]= 2 x 10^40 erg s^-1.

VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5<z<5. At redshift z<1.4 we observe that radio-quiet AGN populate the locus of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- selected populations and show that it can reconcile apparently discrepant measurements presented in the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity.

We present a study of the X-ray to optical properties of a sample of 545 X-ray selected Type 1 AGN, from the XMM-COSMOS survey, over a wide range of redshifts ($0.04<\z<4.25$) and X-ray luminosities ($40.6 \leq \Log \Lhard \leq 45.3$). About 60% of them are spectroscopically identified Type 1 AGN, while the others have a reliable photometric redshift and are classified as Type 1 AGN on the basis of their multi-band Spectral Energy Distributions. We discuss the relationship between UV and X-ray luminosity, as parameterized by the $\alphaox$ spectral slope, and its dependence on redshift and luminosity. We compare our findings with previous investigations of optically selected broad-line AGN (mostly from SDSS). A highly significant correlation between $\alphaox$ and $\lo$ is found, in agreement with previous investigations of optically selected samples. We calculate bolometric corrections, $\kbol$, for the whole sample using hard X-ray luminosities ($\Lhard$), and the Eddington ratios for a subsample of 150 objects for which black hole mass estimates are available. We confirm the trend of increasing bolometric correction with increasing Eddington ratio as proposed in previous works. A tight correlation is found between $\alphaox$ and $\kbol$, which can be used to estimate accurate bolometric corrections using only optical and X-ray data. We find a significant correlation between $\alphaox$ and Eddington ratio, in which $\alphaox$ increases for increasing Eddington ratios.

We constrain the ratio of black hole (BH) mass to total stellar mass of type-1 AGN in the COSMOS survey at 1<z<2. For 10 AGN at mean redshift z~1.4 with both HST/ACS and HST/NICMOS imaging data we are able to compute total stellar mass M_(*,total), based on restframe UV-to-optical host galaxy colors which constrain mass-to-light ratios. All objects have virial BH mass-estimates available from the COSMOS Magellan/IMACS and zCOSMOS surveys. We find zero difference between the M_BH--M_(*,total)-relation at z~1.4 and the M_BH--M_(*,bulge)-relation in the local Universe. Our interpretation is: (a) If our objects were purely bulge-dominated, the M_BH--M_(*,bulge)-relation has not evolved since z~1.4. However, (b) since we have evidence for substantial disk components, the bulges of massive galaxies (logM_(*,total)=11.1+-0.25 or logM_BH~8.3+-0.2) must have grown over the last 9 Gyrs predominantly by redistribution of disk- into bulge-mass. Since all necessary stellar mass exists in the galaxy at z=1.4, no star-formation or addition of external stellar material is required, only a redistribution e.g. induced by minor and major merging or through disk instabilities. Merging, in addition to redistributing mass in the galaxy, will add both BH and stellar/bulge mass, but does not change the overall final M_BH/M_(*,bulge) ratio. Since the overall cosmic stellar and BH mass buildup trace each other tightly over time, our scenario of bulge-formation in massive galaxies is independent of any strong BH-feedback and means that the mechanism coupling BH and bulge mass until the present is very indirect.

Measurements of X-ray scaling laws are critical for improving cosmological constraints derived with the halo mass function and for understanding the physical processes that govern the heating and cooling of the intracluster medium. In this paper, we use a sample of 206 X-ray selected galaxy groups to investigate the scaling relation between X-ray luminosity (Lx) and halo mass (M00) where M200 is derived via stacked weak gravitational lensing. This work draws upon a broad array of multi-wavelength COSMOS observations including 1.64 square degrees of contiguous imaging with the Advanced Camera for Surveys (ACS) and deep XMM-Newton/Chandra imaging. The combined depth of these two data-sets allows us to probe the lensing signals of X-ray detected structures at both higher redshifts and lower masses than previously explored. Weak lensing profiles and halo masses are derived for nine sub-samples, narrowly binned in luminosity and redshift. The COSMOS data alone are well fit by a power law, M200 ~ Lx^a, with a slope of a=0.66+-0.14. These results significantly extend the dynamic range for which the halo masses of X-ray selected structures have been measured with weak gravitational lensing. As a result, tight constraints are obtained for the slope of the M-Lx relation. The combination of our group data with previously published cluster data demonstrates that the M-Lx relation is well described by a single power law, a=0.64+-0.03, over two decades in mass, 10^13.5-10^15.5 h72^-1 Msun. These results are inconsistent at the 3.7 level with the self-similar prediction of a=0.75. We examine the redshift dependence of the M-Lx relation and find little evidence for evolution beyond the rate predicted by self-similarity from z ~ 0.25 to z ~ 0.8.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra program, that covers the central contiguous ~0.92 deg^2 of the COSMOS field. C-COSMOS is the result of a complex tiling, with every position being observed in up to six overlapping pointings (four overlapping pointings in most of the central ~0.45 deg^2 area with the best exposure, and two overlapping pointings in most of the surrounding area, covering an additional ~0.47 deg^2). Therefore, the full exploitation of the C-COSMOS data requires a dedicated and accurate analysis focused on three main issues: 1) maximizing the sensitivity when the PSF changes strongly among different observations of the same source (from ~1 arcsec up to ~10 arcsec half power radius); 2) resolving close pairs; and 3) obtaining the best source localization and count rate. We present here our treatment of four key analysis items: source detection, localization, photometry, and survey sensitivity. Our final procedure consists of a two step procedure: (1) a wavelet detection algorithm, to find source candidates, (2) a maximum likelihood Point Spread Function fitting algorithm to evaluate the source count rates and the probability that each source candidate is a fluctuation of the background. We discuss the main characteristics of this procedure, that was the result of detailed comparisons between different detection algorithms and photometry tools, calibrated with extensive and dedicated simulations.

We present the results of a multifrequency campaign on the high-redshift (z = 3.1) blazar PKS 0537-286. The source was observed at different epochs from 2006 to 2008 with INTEGRAL and Swift, and nearly simultaneously with ground-based near-IR/optical telescopes. The SEDs are compatible with a model based on synchrotron radiation and external inverse Compton scattering. The campaign gives an insight into the physical environment of the blazar.

We present black hole masses and accretion rates for 182 Type 1 AGN in COSMOS. We estimate masses using the scaling relations for the broad Hb, MgII, and CIV emission lines in the redshift ranges 0.16<z<0.88, 1<z<2.4, and 2.7<z<4.9. We estimate the accretion rate using an Eddington ratio L_I/L_Edd estimated from optical and X-ray data. We find that very few Type 1 AGN accrete below L_I/L_Edd ~ 0.01, despite simulations of synthetic spectra which show that the survey is sensitive to such Type 1 AGN. At lower accretion rates the BLR may become obscured, diluted or nonexistent. We find evidence that Type 1 AGN at higher accretion rates have higher optical luminosities, as more of their emission comes from the cool (optical) accretion disk with respect to shorter wavelengths. We measure a larger range in accretion rate than previous works, suggesting that COSMOS is more efficient at finding low accretion rate Type 1 AGN. However the measured range in accretion rate is still comparable to the intrinsic scatter from the scaling relations, suggesting that Type 1 AGN accrete at a narrow range of Eddington ratio, with L_I/L_Edd ~ 0.1.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra} program that has imaged the central 0.5 sq.deg of the COSMOS field (centered at 10h, +02deg) with an effective exposure of ~160ksec, and an outer 0.4sq.deg. area with an effective exposure of ~80ksec. The limiting source detection depths are 1.9e-16 erg cm(-2) s(-1) in the Soft (0.5-2 keV) band, 7.3e(-16) erg cm^-2 s^-1 in the Hard (2-10 keV) band, and 5.7e(-16) erg cm(-2) s(-1) in the Full (0.5-10 keV) band. Here we describe the strategy, design and execution of the C-COSMOS survey, and present the catalog of 1761 point sources detected at a probability of being spurious of <2e(-5) (1655 in the Full, 1340 in the Soft, and 1017 in the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing positions with the ACIS-I imager, a remarkably uniform (to 12%) exposure across the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower flux limit. The widely different PSFs obtained in each exposure at each point in the field required a novel source detection method, because of the overlapping tiling strategy, which is described in a companion paper. (Puccetti et al. Paper II). This method produced reliable sources down to a 7-12 counts, as verified by the resulting logN-logS curve, with sub-arcsecond positions, enabling optical and infrared identifications of virtually all sources, as reported in a second companion paper (Civano et al. Paper III). The full catalog is described here in detail, and is available on-line.

Here we review some of the main issues related to multi-wavelength source identification and characterization, with particular emphasis on the field of X-ray surveys carried out over the last years. This complex and time-consuming process is going to represent one of the main difficulties over the coming years, when significantly larger surveys, both in area and depth, will be carried out with the new generations of space- and ground-based facilities like e.g. eROSITA, WISE, VISTA, Pan-STARRS, and LSST. The Virtual Observatory can offer a reliable way to approach to a new concept of data handling and multi-wavelength source characterization, provided that uniform and rigorous data analyses and extensive quality checks are performed.

The COSMOS survey is a multiwavelength survey aimed to study the evolution of galaxies, AGN and the large scale structure. The XMM-COSMOS is a deep X-ray survey over the full 2 deg2 of the COSMOS area. It consists of 55 XMM-Newton pointings for a total exposure of ~1.5 Ms with an average vignetting corrected depth of 40 ks across the field of view and a sky coverage of 2.13 deg2. We present the catalogue of point-like X-ray sources detected with the EPIC CCD cameras, the logN-logS relations and the X-ray colour-colour diagrams. The analysis was performed in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands. The completeness of the catalogue as well as logN-logS have been calibrated using Monte Carlo simulations. The catalogs contains a total of 1887 unique sources detected in at least one band. The survey, that shows unprecedented homogeneity, has a flux limit of ~1.7x10-15 erg cm-2 s-1, ~9.3x10-15 erg cm-2 s-1 and ~1.3x10-14 erg cm-2 s-1 over 90% of the area (1.92 deg2) in the 0.5-2 keV, 2-10 keV and 5-10 keV energy band, respectively. Thanks to the rather homogeneous exposure over a large area, the derived logN-logS relations are very well determined over the flux range sampled by XMM-COSMOS. These relations have been compared with XRB synthesis models, which reproduce the observations with an agreement of ~10% in the 5-10 keV and 2-10 keV band, while in the 0.5-2 keV band the agreement is of the order of ~20%. The hard X-ray colors confirmed that the majority of the extragalactic sources, in a bright subsample, are actually Type I or Type II AGN. About 20% of the sources have X-ray luminosity typical of AGN (L_X >1042 erg/s) although they do not show any clear signature of nuclear activity in the optical spectrum.

Tidal disruption events are possible sources of temporary nuclear activity in galactic nuclei and can be considered as good indicators of the existence of super massive black holes in the centers of galaxies. A new X-ray source has been serendipitously detected with ROSAT in a PSPC pointed observation of the galaxy cluster A3571. Given the strong flux decay of the object in subsequent detections, the tidal disruption scenario is investigated as possible explanationof the event. We followed up the evolution of the X-ray transient with ROSAT, XMM-Newton and Chandra for a total period of ~13 years. We also obtained 7-band optical/NIR photometry with GROND at the ESO/MPI 2.2m telescope. We report a very large decay of the X-ray flux of ROSAT source identified with the galaxy LEDA 095953 a member of the cluster Abell 3571. We measured a maximum 0.3-2.4 keV luminosity Log(L_X)=42.8 erg s^-1. The high state of the source lasted at least 150 ks, afterwards L_X declined as ~t^-2. The spectrum of the brightest epoch is consistent with a black body with temperature kT ~0.12 keV. The total energy released by this event in 10 yr was estimated to be Delta_E>2x10^50 erg. We interpret this event as a tidal disruption of a solar type star by the central supermassive black hole (i.e. ~10^7 M_sun) of the galaxy.

The impact of environment on AGN activity up to z~1 is assessed by utilizing a mass-selected sample of galaxies from the 10k catalog of the zCOSMOS spectroscopic redshift survey. We identify 147 AGN by their X-ray emission as detected by XMM-Newton from a parent sample of 7234 galaxies. We measure the fraction of galaxies with stellar mass M_*>2.5x10^10 Msun that host an AGN as a function of local overdensity using the 5th, 10th and 20th nearest neighbors that cover a range of physical scales (~1-4 Mpc). Overall, we find that AGNs prefer to reside in environments equivalent to massive galaxies with substantial levels of star formation. Specifically, AGNs with host masses between 0.25-1x10^11 Msun span the full range of environments (i.e., field-to-group) exhibited by galaxies of the same mass and rest-frame color or specific star formation rate. Host galaxies having M_*>10^11 Msun clearly illustrate the association with star formation since they are predominantly bluer than the underlying galaxy population and exhibit a preference for lower density regions analogous to SDSS studies of narrow-line AGN. To probe the environment on smaller physical scales, we determine the fraction of galaxies (M_*>2.5x10^10 Msun) hosting AGNs inside optically-selected groups, and find no significant difference with field galaxies. We interpret our results as evidence that AGN activity requires a sufficient fuel supply; the probability of a massive galaxy to have retained some sufficient amount of gas, as evidence by its ongoing star formation, is higher in underdense regions where disruptive processes (i.e., galaxy harrassment, tidal stripping) are lessened.

We study the spatial clustering of 538 X-ray selected AGN in the 2 deg^2 XMM-COSMOS field that are spectroscopically identified to I_{AB}<23 and span the redshift range z=0.2-3.0. The median redshift and luminosity of the sample are z = 0.98 and L_{0.5-10}=6.3 x 10^{43} erg/s, respectively. A strong clustering signal is detected at ~18sigma level, which is the most significant measurement obtained to date for clustering of X-ray selected AGN. By fitting the projected correlation function w(r_p) with a power law on scales of r_p=0.3-40 Mpc/h, we derive a best fit comoving correlation length of r_0 = 8.6 +- 0.5 Mpc/h and slope of gamma=1.88 +- 0.07 (Poissonian errors; bootstrap errors are about a factor of 2 larger). An excess signal is observed in the range r_p~5-15 Mpc/h, which is due to a large scale structure at z ~ 0.36 containing about 40 AGN. When removing the z ~ 0.36 structure, or computing w(r_p) in a narrower range around the peak of the redshift distribution (e.g. z=0.4-1.6), the correlation length decreases to r_0 ~ 5-6 Mpc/h, which is consistent with that observed for bright optical QSOs at the same redshift. We investigate the clustering properties of obscured and unobscured AGN separately. Within the statistical uncertainties, we do not find evidence that AGN with broad optical lines (BLAGN) cluster differently from AGN without broad optical lines (non-BLAGN). The correlation length measured for XMM-COSMOS AGN at z~1 is similar to that of massive galaxies (stellar mass M_*> 3 x 10^{10} M_sun) at the same redshift. This suggests that AGN at z~1 are preferentially hosted by massive galaxies, as observed both in the local and in the distant (z~2) Universe. (shortened)

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

We present a study of the host galaxies of AGN selected from the zCOSMOS survey to establish if accretion onto supermassive black holes and star formation are explicitly linked up to z~1. We identify 152 galaxies that harbor AGN, based on XMM observations of 7543 galaxies (i<22.5). Star formation rates (SFRs), including those weighted by stellar mass, are determined using the [OII]3727 emission-line, corrected for an AGN contribution. We find that the majority of AGN hosts have significant levels of star formation with a distribution spanning ~1-100 Msun yr^-1. The close association between AGN activity and star formation is further substantiated by an increase in the AGN fraction with the youthfulness of their stars as indicated by the rest-frame color (U-V) and spectral index Dn(4000); we demonstrate that mass-selection alleviates an artifical peak falling in the transition region when using luminosity-limited samples. We also find that the SFRs of AGN hosts evolve with cosmic time in a manner that closely mirrors the overall galaxy population and explains the low SFRs in AGNs (z<0.3) from the SDSS. We conclude that the conditions most conducive for AGN activity are a massive host galaxy and a large reservoir of gas. Furthermore, a direct correlation between mass accretion rate onto SMBHs and SFR is shown to be weak although the average ratio is constant with redshift, effectively shifting the evidence for a co-evolution scenario in a statistical manner to smaller physical scales. Our findings illustrate an intermittent scenario with an AGN lifetime substantially shorter than that of star formation and underlying complexities regarding fueling over vastly different physical scales yet to be determined [Abridged].

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.

We present a new measurement of the space density of high redshift (3.0<z<4.5), X-ray selected QSOs obtained by exploiting the deep and uniform multiwavelength coverage of the COSMOS survey. We have assembled a statistically large (40 objects), X-ray selected (F_{0.5-2 keV} >10^{-15} cgs), homogeneous sample of z>3 QSOs for which spectroscopic (22) or photometric (18) redshifts are available. We present the optical (color-color diagrams) and X-ray properties, the number counts and space densities of the z>3 X-ray selected quasars population and compare our findings with previous works and model predictions. We find that the optical properties of X-ray selected quasars are not significantly different from those of optically selected samples. There is evidence for substantial X-ray absorption (logN_H>23 cm^{-2}) in about 20% of the sources in the sample. The comoving space density of luminous (L_X >10^{44} erg s^-1) QSOs declines exponentially (by an e--folding per unit redshift) in the z=3.0-4.5 range, with a behavior similar to that observed for optically bright unobscured QSOs selected in large area optical surveys. Prospects for future, large and deep X-ray surveys are also discussed.

We present photometric redshifts and spectral energy distribution (SED) classifications for a sample of 1542 optically identified sources detected with XMM in the COSMOS field. Our template fitting classifies 46 sources as stars and 464 as non-active galaxies, while the remaining 1032 require templates with an AGN contribution. High accuracy in the derived photometric redshifts was accomplished as the result of 1) photometry in up to 30 bands with high significance detections, 2) a new set of SED templates including 18 hybrids covering the far-UV to mid-infrared, which have been constructed by the combination of AGN and non-active galaxies templates, and 3) multi-epoch observations that have been used to correct for variability (most important for type 1 AGN). The reliability of the photometric redshifts is evaluated using the sub-sample of 442 sources with measured spectroscopic redshifts. We achieved an accuracy of $\sigma_{\Delta z/(1+z_{spec})} = 0.014$ for i$_{AB}^*<$22.5 ($\sigma_{\Delta z/(1+z_{spec})} \sim0.015$ for i$_{AB}^*<$24.5). The high accuracies were accomplished for both type 2 (where the SED is often dominated by the host galaxy) and type 1 AGN and QSOs out to $z=4.5$. The number of outliers is a large improvement over previous photometric redshift estimates for X-ray selected sources (4.0% and 4.8% outliers for i$_{AB}^*<$22.5 and i$_{AB}^*<$24.5, respectively). We show that the intermediate band photometry is vital to achieving accurate photometric redshifts for AGN, whereas the broad SED coverage provided by mid infrared (Spitzer/IRAC) bands is important to reduce the number of outliers for normal galaxies.

We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies of the host galaxies of approximately 400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts with a sample of non-active galaxies drawn from the COSMOS field to match the magnitude and redshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies span a substantial range that peaks between those of early-type, bulge-dominated and late-type, disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. Unaccounted for, the nuclear point source can significantly bias results of these measured structural parameters, so we subtract the best-fit point source component to obtain images of the underlying host galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions are consistent with those of control galaxies. Combined with a lack of excess companion galaxies around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, but that this transition is not typically triggered by major mergers.

We report about the detection of 10 clusters of galaxies in the ongoing Swift/BAT all-sky survey. This sample, which comprises mostly merging clusters, was serendipitously detected in the 15--55 keV band. We use the BAT sample to investigate the presence of excess hard X-rays above the thermal emission. The BAT clusters do not show significant (e.g. >2 $\sigma$) non-thermal hard X-ray emission. The only exception is represented by Perseus whose high-energy emission is likely due to NGC 1275. Using XMM-Newton, Swift/XRT, Chandra and BAT data, we are able to produce upper limits on the Inverse Compton (IC) emission mechanism which are in disagreement with most of the previously claimed hard X-ray excesses. The coupling of the X-ray upper limits on the IC mechanism to radio data shows that in some clusters the magnetic field might be larger than 0.5 $\mu$G. We also derive the first logN-logS and luminosity function distribution of galaxy clusters above 15 keV.

The Lockman Hole field represents the sky area of lowest Galactic line-of-sight column density N_H=5.7X10^19 cm^-2. It was observed by the XMM-Newton X-ray observatory in 18 pointings for a total of 1.16 Msec (raw EPIC pn observing time) constituting the deepest XMM-Newton exposure so far. We present a catalogue of the X-ray sources detected in the central 0.196 deg^2 of the field and discuss the derived number counts and X-ray colours. In the 0.5--2.0 keV band, a sensitivity limit (defined as the faintest detectable source)of 1.9X10^-16 erg cm^-2 s^-1 was reached. The 2.0--10.0 keV band and 5.0--10.0 keV band sensitivity limits were 9X10^-16 erg cm^-2 s^-1 and 1.8X10^-15 erg cm^-2 s^-1, respectively.A total of 409 sources above a detection likelihood of 10 (3.9 sigma) were found within a radius of 15' off the field centre, of which 340, 266, and 98 sources were detected in the soft, hard, and very hard bands, respectively. The number counts in each energy band are in close agreement with results from previous surveys and with the synthesis models of the X-ray background. A 6% of Compton-thick source candidates have been selected from the X-ray colour-colour diagram. This fraction is consistent with the most recent predictions of X-ray background population synthesis models at our flux limits. We also estimated, for the first time, the logN-logS relation for Compton-thick AGN.

We present the analysis of the source counts in the XMM-COSMOS survey using data of the first year of XMM-Newton observations. The survey covers ~2 deg^2 within the region of sky bounded by 9^h57.5^m<R.A.<10^h03.5^m; 1^d27.5^m<DEC<2^d57.5^m with a total net integration time of 504 ks. Using a maximum likelihood algorithm we detected a total of 1390 sources at least in one band. Using Monte Carlo simulations to estimate the sky coverage we produced the logN-logS relations. These relations have been then derived in the 0.5--2 keV, 2--10 keV and 5--10 keV energy bands, down to flux limits of 7.2x10^-16 erg cm^-2 s^-1, 4.0x10^-15 erg cm^-2 s^-1 and 9.7x10^-15 erg cm^-2 s^-1, respectively. These relations have been compared to previous X-ray survey and to the most recent X-ray background model finding an excellent agreement. The slightly different normalizations observed in the source counts of COSMOS and previous surveys can be largely explained as a combination of low counting statistics and cosmic variance introduced by the large scale structure.

X-ray surveys facilitate investigations of the environment of AGNs. Deep Chandra observations revealed that the AGNs source surface density rises near clusters of galaxies. The natural extension of these works is the measurement of spatial clustering of AGNs around clusters and the investigation of relative biasing between active galactic nuclei and galaxies near clusters.The major aims of this work are to obtain a measurement of the correlation length of AGNs around clusters and a measure of the averaged clustering properties of a complete sample of AGNs in dense environments. We present the first measurement of the soft X-ray cluster-AGN cross-correlation function in redshift space using the data of the ROSAT-NEP survey. The survey covers 9x9 deg^2 around the North Ecliptic Pole where 442 X-ray sources were detected and almost completely spectroscopically identified. We detected a >3sigma significant clustering signal on scales s<50 h70^-1 Mpc. We performed a classical maximum-likelihood power-law fit to the data and obtained a correlation length s_0=8.7+1.2-0.3 h_70-1 Mpc and a slope gamma=1.7$^+0.2_-0.7 (1sigma errors). This is a strong evidence that AGNs are good tracers of the large scale structure of the Universe. Our data were compared to the results obtained by cross-correlating X-ray clusters and galaxies. We observe, with a large uncertainty, that the bias factor of AGN is similar to that of galaxies.

We present the data analysis and the X-ray source counts for the first season of XMM-Newton observations in the COSMOS field. The survey covers ~2 deg^2 within the region of sky bounded by 9^h57.5^m<R.A.<10^h03.5^m; 1^d27.5^m<DEC<2^d57.5^m with a total net integration time of 504 ks. A maximum likelihood source detection was performed in the 0.5-2 keV, 2-4.5 keV and 4.5-10 keV energy bands and 1390 point-like sources were detected in at least one band. Detailed Monte Carlo simulations were performed to fully test the source detection method and to derive the sky coverage to be used in the computation of the logN-logS relations. The 0.5--2 keV and 2--10 keV differential logN-logS were fitted with a broken power-law model which revealed a Euclidean slope (alpha~2.5) at the bright end and a flatter slope (alpha~1.5) at faint fluxes. In the 5--10 keV energy band a single power-law provides an acceptable fit to the observed source counts with a slope alpha~2.4. A comparison with the results of previous surveys shows good agreement in all the energy bands under investigation in the overlapping flux range. We also notice a remarkable agreement between our logN-logS relations and the most recent model of the XRB. The slightly different normalizations observed in the source counts of COSMOS and previous surveys can be largely explained as a combination of low counting statistics and cosmic variance introduced by the large scale structure.

We present the results of a search for galaxy clusters in the first 36 XMM-Newton pointings on the COSMOS field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 3x10-15 ergs cm-2 s-1, having one of the widest XMM-Newton contiguous raster surveys, covering an area of 2.1 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. Verification of the cluster candidates is done based on a galaxy concentration analysis in redshift slices of thickness of 0.1-0.2 in redshift, using the multi-band photometric catalog of the COSMOS field and restricting the search to z<1.3 and i_AB < 25. We identify 72 clusters and derive their properties based on the X-ray cluster scaling relations. A statistical description of the survey in terms of the cumulative log(N>S)-lg(S) distribution compares well with previous results, although yielding a somewhat higher number of clusters at similar fluxes. The X-ray luminosity function of COSMOS clusters matches well the results of nearby surveys, providing a comparably tight constraint on the faint end slope of alpha=1.93+/-0.04. For the probed luminosity range of 8x10+42 - 2x10+44 ergs s-1, our survey is in agreement with and adds significantly to the existing data on the cluster luminosity function at high redshifts and implies no substantial evolution at these luminosities to z=1.3.

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift <1.5, while the redshift distribution of the full BL AGN population peaks at z~1.5. Our analysis also suggests that the Type 2/Type 1 ratio decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral analysis and the most recent modeling of the X-ray luminosity function evolution.

We present a detailed spectral analysis of point-like X-ray sources in the XMM-COSMOS field. Our sample of 135 sources only includes those that have more than 100 net counts in the 0.3-10 keV energy band and have been identified through optical spectroscopy. The majority of the sources are well described by a simple power-law model with either no absorption (76%) or a significant intrinsic, absorbing column (20%).As expected, the distribution of intrinsic absorbing column densities is markedly different between AGN with or without broad optical emission lines. We find within our sample four Type-2 QSOs candidates (L_X > 10^44 erg/s, N_H > 10^22 cm^-2), with a spectral energy distribution well reproduced by a composite Seyfert-2 spectrum, that demonstrates the strength of the wide field XMM/COSMOS survey to detect these rare and underrepresented sources.

We present the first results of the measurements of angular auto-correlation functions (ACFs) of X-ray point sources detected in the XMM-Newton observations of the 2 deg2 COSMOS field (XMM-COSMOS). A significant positive signals have been detected in the 0.5-2 (SFT) band, in the angle range of 0.5-24 arcminutes, while the positive signals were at the 2 and 3 sigma levels in the 2-4.5 (MED) and 4.5-10 (UHD) keV bands respectively. Correctly taking integral constraints into account is a major limitation in interpreting our results. With power-law fits to the ACFs without the integral constraint term, we find correlation lengths of theta_c=1.9+-0.3, 0.8+0.5-0.4 and 6+-2 arcseconds for the SFT, MED, and UHD bands respectively for gamma=1.8. The inferred comoving correlation lengths, also taking into account the bias by the source merging due to XMM-Newton PSF, are 9.8+-0.7, 5.8+1.4-1.7 and 12+-2 h-1 Mpc at the effective redshifts of 1.1, 0.9, and 0.6 for the SFT, MED, and UHD bands respectively. If we include the integral constraint term in the fitting process,assuming that the power-law extends to the scale length of the entire XMM-COSMOSfield, the correlation lengths become larger by 20%-90%. Comparing the inferred rms fluctuations of the spatial distribution of AGNs with those of the underlying mass, the bias parameters of the X-ray source clustering at these effective redshifts are in the range b(AGN)=1.5-4. (Edited to simple ascii.)

We present the first set of XMM-Newton EPIC observations in the 2 square degree COSMOS field. The strength of the COSMOS project is the unprecedented combination of a large solid angle and sensitivity over the whole multiwavelength spectrum. The XMM-Newton observations are very efficient in localizing and identifying active galactic nuclei (AGN) and clusters as well as groups of galaxies. One of the primary goals of the XMM-Newton Cosmos survey is to study the co-evolution of active galactic nuclei as a function of their environment in the Cosmic web. Here we present the log of observations, images and a summary of first research highlights for the first pass of 25 XMM-Newton pointings across the field. In the existing dataset we have detected 1416 new X-ray sources in the 0.5-2, 2-4.5 and 4.5-10 keV bands to an equivalent 0.5-2 keV flux limit of 7x10-16 erg cm-2 s-1. The number of sources is expected to grow to almost 2000 in the final coverage of the survey. From an X-ray color color analysis we identify a population of heavily obscured, partially leaky or reflecting absorbers, most of which are likely to be nearby, Compton-thick AGN.

X-ray surveys facilitate investigations of the environment of AGNs. Deep Chandra observations revealed that the AGNs source surface density rises near clusters of galaxies. The natural extension of these works is the measurement of spatial clustering of AGNs around clusters and the investigation of relative biasing between active galactic nuclei and galaxies near clusters. The major aims of this work are to obtain a measurement of the correlation length of AGNs around clusters and a measure of the averaged clustering properties of a complete sample of AGNs in dense environments. We present the first measurement of the soft X-ray cluster-AGN cross-correlation function in redshift space using the data of the ROSAT-NEP survey. The survey covers 9x9 deg^2 around the North Ecliptic Pole where 442 X-ray sources were detected and almost completely spectroscopically identified. We detected a > 3 sigma significant clustering signal on scales s<50 h_70^-1 Mpc. We performed a classical maximum-likelihood power-law fit to the data and obtained a correlation length s_0=8.7^+1.2_-0.3 h70^-1 Mpc and a slope gamma=1.7^+0.2_-0.7 (1 sigma errors). This is a strong evidence that AGNs are good tracers of the large scale structure of the Universe. Our data were compared to the results obtained by cross-correlating X-ray clusters and galaxies. We observe, with a large uncertainty, a similar behaviour of the AGNs clustering around clusters similar to the clustering of galaxies around clusters.

With the superb angular resolution of the Chandra Observatory, it is now possible to detect X-ray point sources, either embedded in galaxy clusters or along the cluster line of sight, which could not be resolved by previous instruments. This now allows studies of source counts in distant cluster fields. A sample of 18 distant (0.25 < z < 1.01) galaxy clusters from the Chandra archive were used to analyze the inner region of clusters of galaxies to check for the presence of any overdensity of X-ray point sources embedded in the gas diffuse emission. We constructed the logN-logS, in both the soft and hard energy bands, for the X-ray point sources detected in the central cluster region (within 1 Mpc)to be compared with the counts of point sources detected in similarly deep fields without clusters. We find a ~ 2 sigma excess of cluster region sources at the bright end of the logN-logS. The radial distribution of the brightest X-ray point sources confirms this excess and indicates that it is confined to the inner 0.5 Mpc of the cluster region. The results suggest the possible existence of X-ray sources belonging to the cluster (most probably AGN, given their 0.5-10 keV luminosity ranging from 10^43 to 10^44 erg/s): on average one every three clusters. The X-ray source excess found here is much smaller than the excess of radio galaxies found recently in high-z X-ray selected clusters, possibly due to the higher sensitivity of the radio observations.

A comparative XMM-Newton archival data spectral study of the ultraluminous X-ray source (ULX) M33 X-8 has been performed by using both the standard disk model and the newly developed slim disk models. The results of this analysis tend to confirm the hypothesis that M33 X-8 is an X-ray binary with a stellar-mass black hole accreting at super-Eddington rate.

The first systematic study of the serendipitous X-ray source density around 10 high z (0.24<z<1.2) clusters performed with Chandra is presented. A factor \~2 overdensity has been found in 4 cluster fields. The result is statistically highly significant (at >2sigma per field, and >5sigma overall) only on scale of 8'X8' and peculiar to cluster fields. The blank field fluctuations, i.e. the most sensitive measurements to date of the X-ray cosmic variance on this scale, indicate (1 sigma) variations within ~15-25% around the average value. The first marginal evidence that these overdensities increase with cluster redshift is also presented. We speculate that the most likely explanation is that the sources in excess are AGNs tracing filaments connected to the clusters. If the association of the overdensities with large-scale structures and their positive relation with cluster redshift is confirmed, these studies could represent a new, direct, way of mapping the cosmic web and its cosmological parameters at high z.

The first systematic study of the serendipitous X-ray source density around 10 high z (0.24<z<1.2) clusters has been performed with Chandra. A factor ~2 overdensity has been found in 4 cluster fields, increasing to 11 the number of known source overdensities around high z clusters. The result is statistically highly significant (at >2 sigma per field, and >5 sigma overall) only on scales 8'X8' and peculiar to cluster fields. The blank field fluctuations, i.e. the most sensitive measurements to date of the X-ray cosmic variance on this scale, indicate variations within ~15-25% around the average value. The first marginal evidence that these overdensities increase with cluster redshift is also presented. We speculate that the most likely explanation for what has been observed is that the x-ray sources are AGNs which trace filaments connected to the clusters. Nevertheless other possible explanations such as an ``X-ray Butcher-Oemler'' effect are not ruled out by present results. If the association of the overdensities with large-scale structures and their positive relation with cluster redshift is confirmed, these studies could represent a new, direct, way of mapping the cosmic web and its parameters at high z.