A statistical study of intermediate Palomar Transient Factory supernovae (SNe) in Type 1 AGN has shown a major deficit of supernovae around Type 1 AGN host galaxies, with respect to Type 2 AGN hosts. The aim of this work is to test whether there is any preference for Type 1 AGN to host SN of a specific kind. Through the analysis of SN occurrence and their type (thermonuclear vs core-collapse), we can directly link the type of stars producing the SN events, thus this is an indirect way to study host galaxies in Type 1 AGN. We examine the detection fractions of SNe, the host galaxies and compare the sample properties to typical host galaxies in the Open Supernova Catalog (OSC; Guillochon et al. 2017). The majority of the host galaxies in the AGN sample are late-type, similar to typical galaxies hosting SN within the OSC. The findings are supportive of a deficiency of SNe near Type 1 AGN, although we cannot with certainty assess the overall detection fractions of SNe in Type 1 AGN relative to other SN host galaxies. We can state that Type 1 AGN has equal detection fractions of thermonuclear vs core-collapse SNe. However, we note the possibility of a higher detection rate of core-collapse supernovae in Type-1 AGN with insecure AGN classifications.

We have recently used the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey to show that red quasars have fundamentally different radio properties to typical blue quasars: a significant (factor $\sim3$) enhancement in the radio-detection fraction, which arises from systems around the radio-quiet threshold with compact ($<5''$) radio morphologies. To gain greater insight into these physical differences, here we use the DR14 Sloan Digital Sky Survey (SDSS) and more sensitive, higher resolution radio data from the Very Large Array (VLA) Stripe 82 (S82) and VLA-COSMOS 3 GHz (C3GHz) surveys. With the S82 data, we perform morphological analyses at a resolution and depth three times that of the FIRST radio survey, and confirm an enhancement in radio-faint and compact red quasars over typical quasars; we now also find tentative evidence for an enhancement in red quasars with slightly extended radio structures ($16-43$ kpc at $z=1.5$). These analyses are complemented by C3GHz, which is deep enough to detect radio emission from star-formation processes. From our data we find that the radio enhancement from red quasars is due to AGN activity on compact scales ($< 43$ kpc) for radio-intermediate-radio-quiet sources ($-5<R<-3.4$, where $R=L_{1.4GHz}/L_{6 \mu m}$), which decreases at $R<-5$ as the radio emission from star-formation starts to dilute the AGN component. Overall our results argue against a simple orientation scenario and are consistent with red quasars representing a younger, earlier phase in the overall evolution of quasars.

Red quasi-stellar objects (QSOs) are a subset of the luminous end of the cosmic population of active galactic nuclei (AGN), most of which are reddened by intervening dust along the line-of-sight towards their central engines. In recent work from our team, we developed a systematic technique to select red QSOs from the Sloan Digital Sky Survey (SDSS), and demonstrated that they have distinctive radio properties using the Faint Images of the Radio Sky at Twenty centimeters (FIRST) radio survey. Here we expand our study using low-frequency radio data from the LOFAR Two-metre Sky Survey (LoTSS). With the improvement in depth that LoTSS offers, we confirm key results: compared to a control sample of normal "blue" QSOs matched in redshift and accretion power, red QSOs have a higher radio detection rate and a higher incidence of compact radio morphologies. For the first time, we also demonstrate that these differences arise primarily in sources of intermediate radio-loudness: radio-intermediate red QSOs are $\times 3$ more common than typical QSOs, but the excess diminishes among the most radio-loud and the most radio-quiet systems in our study. We develop Monte-Carlo simulations to explore whether differences in star formation could explain these results, and conclude that, while star formation is an important source of low-frequency emission among radio-quiet QSOs, a population of AGN-driven compact radio sources is the most likely cause for the distinct low-frequency radio properties of red QSOs. Our study substantiates the conclusion that fundamental differences must exist between the red and normal blue QSO populations.

Our group has developed a technique that makes use of the observed non-linear relation between the ultraviolet and the X-ray luminosity in quasars to provide an independent measurement of their distances, thus turning quasars into standardizable candles. This technique, at present, it is mostly based upon quasar samples with data from public catalogues both in the X-rays and in the optical/ultraviolet and extends the Hubble diagram of supernovae to a redshift range still poorly explored (z>2). From the X-ray perspective, we are now on the eve of a major change, as the upcoming mission eROSITA is going to provide us with up to ~3 millions of active galactic nuclei across the entire sky. Here we present predictions for constraining cosmological parameters such as the amount of dark matter ($\Omega_{\rm m}$), dark energy ($\Omega_\Lambda$) and the evolution of the equation of state of dark energy ($w$) through the Hubble diagram of quasars, based on the 4-year eROSITA all-sky survey. Our simulations show that the eROSITA quasars, complemented by redshift and broad-band photometric information, will supply the largest quasar sample at z<2, but with very few objects available for cosmology at higher redshifts that survives the cut for the Malmquist bias, as eROSITA will sample the brighter end of the X-ray luminosity function. The power of the quasar Hubble diagram for precision cosmology lies in the high-redshift regime, where quasars can be observed up to redshift ~7.5, essential to discriminate amongst different model extrapolations. Therefore, to be competitive for cosmology, the eROSITA quasar Hubble diagram must be complemented with the already available quasar samples and dedicated (deep) large programmes at redshift z>3.

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

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

(abridged) We present the X-ray analysis of a sample of 30 luminous quasars at $z\simeq3.0-3.3$ with deep XMM-Newton observations, selected from the SDSS-DR7 to be representative of the most luminous, intrinsically blue quasar population. By construction, the sample boasts a unique degree of homogeneity in terms of optical/UV properties. In the X-rays, only four sources are too faint for a detailed spectral analysis. Neglecting a radio-loud object, the other 25 quasars are, as a whole, the most X-ray luminous ever observed, with rest-frame 2-10 keV luminosities of $0.5-7\times10^{45}$ erg/s. The continuum photon index distribution, centred at $\Gamma\sim1.85$, is in excellent agreement with those in place at lower redshift, luminosity and black-hole mass, confirming the universal nature of the X-ray emission mechanism in quasars. Even so, when compared against the well-known $L_{\rm X}-L_{\rm UV}$ correlation, our quasars unexpectedly split into two distinct subsets. About 2/3 of the sources are clustered around the relation with a minimal scatter of 0.1 dex, while the remaining 1/3 appear to be X-ray underluminous by factors of $>3-10$. Such a large incidence ($\approx25\%$) of X-ray weakness has never been reported in radio-quiet, non-BAL quasar samples. Several factors could contribute to enhance the X-ray weakness fraction among our $z\simeq3$ blue quasars. However, the X-ray weak objects also have, on average, flatter spectra, with no clear evidence of absorption. Indeed, column densities in excess of a few $\times10^{22}$ cm$^{-2}$ can be ruled out for most of the sample. We suggest that, at least in some of our X-ray weak quasars, the corona might experience a radiatively inefficient phase due to the presence of a powerful accretion-disc wind, which substantially reduces the accretion rate through the inner disc and so the availability of seed photons for Compton up-scattering.

We present the design, methods, and first results of the MUSE Analysis of Gas around Galaxies (MAGG) survey, a large programme on the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) which targets 28 z > 3.2 quasars to investigate the connection between optically-thick gas and galaxies at z~3-4. MAGG maps the environment of 52 strong absorption line systems at z > 3, providing the first statistical sample of galaxies associated with gas-rich structures in the early Universe. In this paper, we study the galaxy population around a very metal poor gas cloud at z~3.5 towards the quasar J124957.23-015928.8. We detect three Lyman alpha emitters within <200km/s of the cloud redshift, at projected separations <185 kpc (physical). The presence of star-forming galaxies near a very metal-poor cloud indicates that metal enrichment is still spatially inhomogeneous at this redshift. Based on its very low metallicity and the presence of nearby galaxies, we propose that the most likely scenario for this LLS is that it lies within a filament which may be accreting onto a nearby galaxy. Taken together with the small number of other LLSs studied with MUSE, the observations to date show a range of different environments near strong absorption systems. The full MAGG survey will significantly expand this sample and enable a statistical analysis of the link between gas and galaxies to pin down the origin of these diverse environments at z~3-4.

A tight non-linear relation between the X-ray and the optical-ultraviolet (UV) emission has been observed in Active Galactic Nuclei (AGN) over a wide range of redshift and several orders of magnitude in luminosity, suggesting the existence of an ubiquitous physical mechanism regulating the energy transfer between the accretion disc and the X-ray emitting corona. Recently, our group developed a method to use this relation in the observational cosmology, turning quasars into standardizable candles. This work has the main aim to investigate the potential evolution of this correction at high redshifts. We thus studied the $L_{\rm X}-L_{\rm UV}$ relation for a sample of quasars in the redshift range 4<$z$<7, adopting the selection criteria proposed in our previous work regarding their spectral properties. The resulting sample consists of 53 Type 1 (unobscured) quasars, observed either with Chandra or XMM-Newton, for which we performed a full spectral analysis, determining the rest-frame 2 keV flux density, as well as more general X-ray properties such as the estimate of photon index, and the soft (0.5-2 keV) and hard (2-10 keV) unabsorbed luminosities. We find that the relation shows no evidence for evolution with redshift. The intrinsic dispersion of the L$_X$-L$_{UV}$ for a sample free of systematics/contaminants is of the order of 0.22 dex, which is consistent with previous estimates from our group on quasars at lower redshift.

We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect.

In the current framework, the standard parametrization of our Universe is the so-called Lambda Cold Dark Matter ({\Lambda}CDM) model. Recently, Risaliti & Lusso (2019) have shown a ~4{\sigma} tension with the {\Lambda}CDM model through a model-independent parametrization of a Hubble Diagram of supernovae Ia (SNe Ia) from the JLA survey and quasars. Model-independent approaches and independent samples over a wide redshift range are key to testing this tension and any possible systematics. Here we present an analysis of a combined Hubble Diagram of SNe Ia, quasars, and gamma-ray bursts (GRBs) to check the agreement of the quasar and GRB cosmological parameters at high redshifts (z>2) and to test the concordance flat {\Lambda}CDM model with improved statistical accuracy. We build a Hubble diagram with SNe Ia from the Pantheon sample (Scolnic et al. 2018), quasars from the Risaliti & Lusso (2019) sample, and GRBs from the Demianski et al. (2017a) sample, where quasars are standardised through the observed non-linear relation between their ultraviolet and X-ray emission and GRBs through the correlation between the spectral peak energy and the isotropic-equivalent radiated energy (the so-called "Amati relation"). We fit the data with cosmographic models consisting of a fourth-order logarithmic polynomial and a fifth-order linear polynomial, and compare the results with the expectations from a flat {\Lambda}CDM model. We confirm the tension between the best fit cosmographic parameters and the {\Lambda}CDM model at ~4{\sigma} with SNe Ia and quasars, at ~2{\sigma} with SNe Ia and GRBs, and at >4{\sigma} with the whole SNe Ia+quasars+GRB data set. The completely independent high-redshift Hubble diagrams of quasars and GRBs are fully consistent with each other, strongly suggesting that the deviation from the standard model is not due to unknown systematic effects but to new physics.

A minority of the optically selected quasar population are red at optical wavelengths due to the presence of dust along the line-of-sight. A key focus of many red quasar studies is to understand their relationship with the overall quasar population: are they blue quasars observed at a (slight) inclination angle or do they represent a transitional phase in the evolution of quasars? Identifying fundamental differences between red and blue quasars is key to discriminate between these two paradigms. To robustly explore this, we have uniformly selected quasars from the Sloan Digital Sky Survey with mid-infrared counterparts, carefully controlling for luminosity and redshift effects. We take a novel approach to distinguish between colour-selected quasars in the redshift range of $0.2 < z < 2.4$ by constructing redshift-sensitive $g^* - i^*$ colour cuts. From cross-matching this sample to the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey, we have found a factor $\approx$ 3 larger fraction of radio-detected red quasars with respect to that of blue quasars. Through a visual inspection of the FIRST images and an assessment of the radio luminosities (rest-frame ${L_{\rm 1.4~GHz}}$ and ${L_{\rm 1.4~GHz}}/{L_{\rm 6\mu m}}$) we find that the radio-detection excess for red quasars is primarily due to compact and radio-faint systems (around the radio quiet--radio loud threshold). We show that our results rule out orientation as the origin for the differences between red and blue quasars and argue that they provide broad agreement with an evolutionary model.

We present the X-ray properties of 108 Dust-Obscured Galaxies (DOGs; F$_{24 \mu m}$/F$_{R} >$ 1000) in the COSMOS field, all of which detected in at least three far-infrared bands with the Herschel Observatory. Out of the entire sample, 22 are individually detected in the hard 2-8 keV X-ray band by the Chandra COSMOS Legacy survey, allowing us to classify them as AGN. Of them, 6 (27%) are Compton Thick AGN candidates with column densities N$_{H}$$>$10$^{24}$ cm$^{-2}$ while 15 are moderately obscured AGNs with 10$^{22}$ $<$ N$_{H}$ $<$ 10$^{24}$ cm$^{-2}$. Additionally, we estimate AGN contributions to the IR luminosity (8-1000$\mu$m rest-frame) greater than 20% for 19 DOGs based on SED decomposition using Spitzer/MIPS 24$\mu$m and the five Herschel bands (100-500 $\mu$m). Only 7 of these are detected in X-rays individually. We performed a X-ray stacking analysis for the 86 undetected DOGs. We find that the AGN fraction in DOGs increases with 24$\mu$m flux and that it is higher than that of the general 24$\mu$m population. However, no significant difference is found when considering only X-ray detections. This strongly motivates the combined use of X-ray and far-IR surveys to successfully probe a wider population of AGNs, particularly for the most obscured ones.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations.

We recently proposed the equivalent width (EW) of the narrow [OIII]5007{\AA} emission line as an orientation indicator for active galactic nuclei. We tested this method on about $12,300$ optically selected broad line quasars from the SDSS 7th Data Release at redshift $z < 0.8$ and with full width at half-maximum values of broad emission lines (H$\alpha$, H$\beta$, and MgII) larger than $2000$ km/s. We now examine their spectral energy distributions (SEDs) using broad band photometry from the near-infrared to the ultraviolet to look for variations in the overall shape as a function of the EW[OIII]. We find that quasars with low EW[OIII] values (close to face-on position) have flatter near-infrared SEDs with respect to sources with high EW[OIII] values (almost edge-on). Moreover, quasars with high EW[OIII] values show a factor of $\sim 2$ lower emission in the UV to quasars with low EW[OIII] values. Our findings indicate that the torus is clumpy and, on average, co-axial with the accretion disc and broad line region, in agreement with the most recent theoretical models for the obscuring torus.

We investigate the multi-wavelength properties of host galaxies of 3701 X-ray-selected active galactic nuclei (AGNs) out to z~5 in the Chandra-COSMOS Legacy Survey. Thanks to the extensive multi-wavelength photometry available in the COSMOS field, we derive AGN luminosities, host stellar masses, and star formation rates (SFRs) via a multi-component SED fitting technique. Type 1 and Type 2 AGNs follow the same intrinsic Lx-L6um relation, suggesting that mid-infrared emission is a reasonably good measure of the AGN accretion power regardless of obscuration. We find that there is a strong increase in Type 1 AGN fraction toward higher AGN luminosity, possibly due to the fact that Type 1 AGNs tend to be hosted by more massive galaxies. The AGN luminosity and SFR are consistent with an increase toward high stellar mass, while both the Mstellar-dependence is weaker towards the high-mass end, which could be interpreted as a consequence of quenching both star formation and AGN activity in massive galaxies. AGN host galaxies tend to have SFRs that are consistent with normal star-forming galaxies, independent of AGN luminosities. We confirm that black hole accretion rate and SFR are correlated up to z~5, when forming stars. The majority (~73%) of our AGN sample are faint in the far-infrared, implying that the moderate-luminosity AGNs seem to be still active after the star formation is suppressed. It is not certain whether AGN activity plays a role in quenching the star formation. We conclude that both AGN activity and star formation might be more fundamentally related to host stellar mass.

We have recently demonstrated that the non-linear relation between ultraviolet and X-ray luminosity in quasars is very tight (with an intrinsic dispersion of ~0.2 dex), once contaminants (e.g. dust reddening, X-ray absorption), variability, and differences in the active galactic nuclei (AGN) physical properties are taken into account. This relation has thus the great potential to advance our understanding in both supermassive black hole accretion physics and observational cosmology, by targeting a single class of objects. Here we focus on the various contributions to the observed X-ray variability in a homogenous sample of 791 quasars selected from SDSS-DR7 with X-ray data from the 3XMM-DR7 source catalogue. The 250 quasars in this cleaned data set with at least two X-ray observations typically vary with a standard deviation of fractional variation of 15-30% on timescales of weeks/years. Yet, when the count rates are computed at progressively smaller off-axis values, the same quantity is reduced to roughly 10-25%. This suggests that, when estimating variability indicators, part of the quoted variability amplitude could be due to instrumental/calibration issues rather than true variations in the quasar emission.

The concordance (LambdaCDM) model reproduces the main current cosmological observations assuming the validity of general relativity at all scales and epochs, the presence of cold dark matter, and of a cosmological constant, equivalent to a dark energy with constant density in space and time. However, the LambdaCDM model is poorly tested in the redshift interval between the farthest observed Type Ia supernovae5 and that of the Cosmic Microwave background (CMB). We present new measurements of the expansion rate of the Universe in the range 0.5<z<5.5 based on a Hubble diagram of quasars. The quasar distances are estimated from their X-ray and ultraviolet emission, following a method developed by our group. The distance modulus-redshift relation of quasars at z<1.4 is in agreement with that of supernovae and with the concordance model. Yet, a deviation from the LambdaCDM model emerges at higher redshift, with a statistical significance of ~4 sigma. If an evolution of the dark energy equation of state is allowed, the data suggest a dark energy density increasing with time.

We report on the detection of a blueshifted Fe K absorption feature in two consecutive XMM-Newton observations of the luminous blue quasar Ton 28, at the 4$\sigma$ cumulative significance. The rest energy of 9.2 keV implies the presence of an accretion disc wind with bulk outflow velocity of $\sim$0.28$c$, while the kinetic power is most likely a few per cent of the quasar luminosity. Remarkably, Ton 28 had been specifically selected as an optimal target to reveal an ultra-fast X-ray wind based on its total luminosity ($L_{\rm bol} > 10^{46}$ erg s$^{-1}$) and [OIII] $\lambda$5007\AA\ equivalent width (EW < 6 \AA), suggestive of high accretion rate and low inclination, respectively. Other peculiar optical/UV emission-line properties include narrow H$\beta$, strong FeII and blueshifted CIV. These are key parameters in the Eigenvector 1 formalism, and are frequently found in active galaxies with ongoing accretion disc winds, hinting at a common physical explanation. Provided that the effectiveness of our selection method is confirmed with similar sources, this result could represent the first step towards the characterization of black-hole winds through multiwavelength indicators in the absence of high-quality X-ray spectra.

Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission). These quasars have a median redshfit of z=3.17, absolute $i$ magnitude in the range $-29.67\leq M_i(z=2)\leq-27.03$, and different levels of radio-loudness. This sample unveils diverse specimens of Ly$\alpha$ nebulosities extending for tens of kiloparsecs around these quasars (on average out to a maximum projected distance of 80 kpc) above a surface brightness SB$>8.8\times10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ ($2\sigma$). The bulk of the extended Ly$\alpha$ emission is within R<50 kpc, and is characterized by relatively quiescent kinematics, with average velocity dispersions of $\langle \sigma_{\rm Ly\alpha}\rangle < 400$~km~s$^{-1}$. Therefore, the motions within all these Ly$\alpha$ nebulosities have amplitudes consistent with gravitational motions expected in dark matter halos hosting quasars at these redshifts, possibly reflecting the complexity in propagating a fast wind on large scales. Our current data suggest a combination of photoionization and resonant scattering as powering mechanisms of the Ly$\alpha$ emission. We discover the first $z\sim3$ ELAN, which confirms a very low probability ($\sim1\%$) of occurrence of such extreme systems at these cosmic epochs. Finally, we discuss the redshift evolution currently seen in extended Ly$\alpha$ emission around radio-quiet quasars from $z\sim3$ to $z\sim2$, concluding that it is possibly linked to a decrease of cool gas mass within the quasars' CGM from $z\sim3$ to $z\sim2$, and thus to the balance of cool vs hot media.

The standard AGN-galaxy co-evolutionary scenario predicts a phase of deeply buried supermassive black hole growth coexisting with a starburst (SB) before feedback phenomena deplete the cold molecular gas reservoir of the galaxy and an optically luminous QSO is revealed ('SB-QSO evolutionary sequence'). The aim of this work is to measure the cold gas reservoir of three highly obscured QSOs to test if their gas fraction is similar to that of sub-millimeter galaxies (SMGs), as expected by some models, and place these measurements in the context of the SB-QSO framework. We target CO(1-0) transition in BzK4892, a Compton Thick (CT) QSO at z=2.6, CO(1-0) in BzK8608 and CO(2-1) in CDF153, two highly obscured QSOs at z=2.5 and z=1.5, respectively. For all these targets, we place 3$\sigma$ upper limits on the CO, with $L'_{CO} < (1.5\div 2.8)\times 10^{10}$ K km/s pc$^2$. We also compare the molecular gas conditions of our targets with those of other systems at z>1, considering normal star forming galaxies and SMGs, unobscured and obscured AGN from the literature. For the AGN samples, we provide an updated and (almost) complete collection of targets with CO follow-up. BzK4892 displays a high star formation efficiency (SFE$=L_{IR}/L'_{CO}>410$ L$_{\odot}$/(K km s$^{-1}$ pc$^2$)) and a gas fraction $f_{gas}<0.1$. Less stringent constraints are derived for the other two targets ($f_{gas}<0.5$ and SFE$>10$). From the comparison with literature data, we found that a) obscured AGN at z>1 are associated with higher SFE and lower $f_{gas}$ with respect to star forming galaxies; b) mildly and highly obscured active galaxies have comparable gas fractions; c) the SFE of CT and obscured AGN are similar to those of unobscured AGN. Within the SB-QSO framework, these findings could be consistent with a scenario where feedback can impact the host galaxy already from the early phases of the SB-QSO sequence.

We present a search for nuclear X-ray emission in the brightest cluster galaxies (BCGs) of a sample of groups and clusters of galaxies extracted from the Chandra archive. The exquisite angular resolution of Chandra allows us to obtain robust photometry at the position of the BCG, and to firmly identify unresolved X-ray emission when present, thanks to an accurate characterization of the extended emission at the BCG position. We consider two redshift bins (0.2<z<0.3 and 0.55<z<0.75) and analyze all the clusters observed by Chandra with exposure time larger than 20 ks. Our samples have 81 BCGs in 73 clusters and 51 BCGs in 49 clusters in the low- and high-redshift bin, respectively. X-ray emission in the soft (0.5-2 keV) or hard (2-7 keV) band is detected only in 14 and 9 BCGs ($\sim 18$% of the total samples), respectively. The X-ray photometry shows that at least half of the BCGs have a high hardness ratio, compatible with significant intrinsic absorption. This is confirmed by the spectral analysis with a power law model plus intrinsic absorption. We compute the fraction of X-ray bright BCGs above a given hard X-ray luminosity, considering only sources with positive photometry in the hard band (12/5 sources in the low/high-z sample). In the 0.2<z<0.3 interval the hard X-ray luminosity ranges from $10^{42}$ to $7 \times 10^{43}$ erg s$^{-1}$, with most sources found below $10^{43}$ erg s$^{-1}$. In the $0.55<z<0.75$ range, we find a similar distribution of luminosities below $\sim 10^{44}$ erg s$^{-1}$, plus two very bright sources of a few $10^{45}$ erg s$^{-1}$ associated with two radio galaxies. We also find that X-ray luminous BCGs tend to be hosted by cool-core clusters, despite the majority of cool cores do not host nuclear X-ray emission (Abridged).

We present the first results from our survey of intervening and proximate Lyman limit systems (LLSs) at $z$$\sim$2.0-2.5 using the Wide Field Camera 3 on-board the Hubble Space Telescope. The quasars in our sample are projected pairs with proper transverse separations $R_\perp$$\leq$150 kpc and line of sight velocity separations $\lesssim$11,000 km/s. We construct a stacked ultraviolet (rest-frame wavelengths 700-2000\AA) spectrum of pairs corrected for the intervening Lyman forest and Lyman continuum absorption. The observed spectral composite presents a moderate flux excess for the most prominent broad emission lines, a $\sim$30% decrease in flux at $\lambda$=800-900\AA\ compared to a stack of brighter quasars not in pairs at similar redshifts, and lower values of the mean free path of the HI ionizing radiation for pairs ($\lambda_{\rm mfp}^{912}=140.7\pm20.2~h_{70}^{-1}$Mpc) compared to single quasars ($\lambda_{\rm mfp}^{912}=213.8\pm28~h_{70}^{-1}$Mpc) at the average redshift $z\simeq2.44$. From the modelling of LLS absorption in these pairs, we find a higher ($\sim$20%) incidence of proximate LLSs with $\log N_{\rm HI}\geq17.2$ at $\delta v$$<$5,000 km/s compared to single quasars ($\sim$6%). These two rates are different at the 5$\sigma$ level. Moreover, we find that optically-thick absorbers are equally shared between foreground and background quasars. Based on these pieces of evidence, we conclude that there is a moderate excess of gas absorbing Lyman continuum photons in our closely-projected quasar pairs compared to single quasars. We argue that this gas arises mostly within large-scale structures or partially neutral regions inside the dark matter haloes where these close pairs reside.

We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.

The formal division of the distance indicators into primary and secondary leads to difficulties in description of methods which can actually be used in two ways: with, and without the support of the other methods for scaling. Thus instead of concentrating on the scaling requirement we concentrate on all methods of distance determination to extragalactic sources which are designated, at least formally, to use for individual sources. Among those, the Supernovae Ia is clearly the leader due to its enormous success in determination of the expansion rate of the Universe. However, new methods are rapidly developing, and there is also a progress in more traditional methods. We give a general overview of the methods but we mostly concentrate on the most recent developments in each field, and future expectations.

We propose a modified version of the observed non-linear relation between the X-ray (2 keV) and the ultraviolet (2500 \AA) emission in quasars (i.e. $L_{\rm X}\propto L_{\rm UV}^{\gamma}$) which involves the full width at half-maximum, FWHM, of the broad emission line, i.e. $L_{\rm X}\propto L_{\rm UV}^{\hat\gamma}~FWHM^{\hat\beta}$. By analysing a sample of 550 optically selected non-jetted quasars in the redshift range of 0.36-2.23 from the Sloan Digital Sky Survey cross matched with the XMM-Newton catalogue 3XMM-DR6, we found that the additional dependence of the observed $L_{\rm X}-L_{\rm UV}$ correlation on the FWHM of the MgII broad emission line is statistically significant. Our statistical analysis leads to a much tighter relation with respect to the one neglecting FWHM, and it does not evolve with redshift. We interpret this new relation within an accretion disc corona scenario where reconnection and magnetic loops above the accretion disc can account for the production of the primary X-ray radiation. For a broad line region size depending on the disc luminosity as $R_{\rm blr}\propto L_{\rm disc}^{0.5}$, we find that $L_{\rm X}\propto L_{\rm UV}^{4/7} ~FWHM^{4/7}$, which is in very good agreement with the observed correlation.

The cosmological model is at present not tested between the redshift of the farthest observed supernovae (z ~ 1.4) and that of the Cosmic Microwave Background (z ~ 1,100). Here we introduce a new method to measure the cosmological parameters: we show that quasars can be used as "standard candles" by employing the non-linear relation between their intrinsic UV and X-ray emission as an absolute distance indicator. We built a sample of ~ 1,900 quasars with available UV and X-ray observations, and produced a Hubble Diagram up to z ~ 5. The analysis of the quasar Hubble Diagram, when used in combination with supernovae, provides robust constraints on the matter and energy content in the cosmos. The application of this method to forthcoming, larger quasar samples, will also provide tight constraints on the dark energy equation of state and its possible evolution with time.

We present an analysis of the average spectral properties of 12,000 SDSS quasars as a function of accretion disc inclination, as measured from the equivalent width of the [O III] 5007{\AA} line. The use of this indicator on a large sample of quasars from the SDSS DR7 has proven the presence of orientation effects on the features of UV/optical spectra, confirming the presence of outflows in the NLR gas and that the geometry of the BLR is disc-like. Relying on the goodness of this indicator, we are now using it to investigate other bands/components of AGN. Specifically, the study of the UV/optical/IR SED of the same sample provides information on the obscuring "torus". The SED shows a decrease of the IR fraction moving from face-on to edge-on sources, in agreement with models where the torus is co-axial with the accretion disc. Moreover, the fact we are able to observe the broad emission lines also in sources in an edge-on position, suggests that the torus is rather clumpy than smooth as in the Unified Model. The behaviour of the SED as a function of EW[OIII] is in agreement with the predictions of the clumpy torus models as well.

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.

A tight non-linear relation exists between the X-ray and UV emission in quasars (i.e. $L_{\rm X}\propto L_{\rm UV}^{\gamma}$), with a dispersion of $\sim$0.2~dex over $\sim$3~orders of magnitude in luminosity. Here we propose a modified version of this relation which involves the broad emission lines' full-width at half maximum, $L_{\rm X}\propto L_{\rm UV}^{\hat\gamma}\upsilon_{\rm fwhm}^{\hat\beta}$. We interpret this new relation through a simple, {\it ad-hoc} model of accretion disc corona, derived from previous works where it is assumed that reconnection and magnetic loops above the accretion disc can account for the production of the primary X-ray radiation. We find that the monochromatic optical-UV (2500 \AA) and X-ray (2 keV) luminosities depend on the black hole mass and accretion rate as $L_{\rm UV}\propto M_{\rm BH}^{4/3} (\dot{M}/\dot{M}_{\rm Edd})^{2/3}$ and $L_{\rm X}\propto M_{\rm BH}^{19/21} (\dot{M}/\dot{M}_{\rm Edd})^{5/21}$, respectively. Assuming a broad line region size function of the disc luminosity $R_{\rm blr}\propto L_{\rm disc}^{0.5}$ we finally have that $L_{\rm X}\propto L_{\rm UV}^{4/7} \upsilon_{\rm fwhm}^{4/7}$. Such relation is remarkably consistent with the slopes and the normalization obtained from a fit of a sample of 545 optically selected quasars from SDSS DR7 cross matched with the latest XMM-Newton catalogue 3XMM-DR6. The homogeneous sample used here has a dispersion of 0.21 dex, which is much lower than previous works in the literature and suggests a tight physical relation between the accretion disc and the X-ray emitting corona. The proposed new relation does not evolve with time, and thus it can be employed as a cosmological indicator to robustly estimate cosmological parameters.

The non-linear relation between X-ray and UV luminosity in quasars can be used to estimate their distance. Recently, we have shown that despite the large dispersion of the relation, a Hubble Diagram made of large samples of quasars can provide unique constraints on cosmology at high redshift. Furthermore, the dispersion of the relation is heavily affected by measurement errors: until now we have used serendipitous X-ray observations, but dedicated observations would significantly increase the precision of the distance estimates. We discuss the future role of XMM in this new field, showing (1) the fundamental contribution of the Serendipitous Source Catalogue and of large surveys, and (2) the breakthrough advancements we may achieve with the observation of a large number of SDSS quasars at high redshift: every ~10 quasars observed at z$\sim$3 would be equivalent to discovering a supernova at that redshift.

We present AGNfitter, a publicly available open-source algorithm implementing a fully Bayesian Markov Chain Monte Carlo method to fit the spectral energy distributions (SEDs) of active galactic nuclei (AGN) from the sub-mm to the UV, allowing one to robustly disentangle the physical processes responsible for their emission. AGNfitter makes use of a large library of theoretical, empirical, and semi-empirical models to characterize both the nuclear and host galaxy emission simultaneously. The model consists of four physical emission components: an accretion disk, a torus of AGN heated dust, stellar populations, and cold dust in star forming regions. AGNfitter determines the posterior distributions of numerous parameters that govern the physics of AGN with a fully Bayesian treatment of errors and parameter degeneracies, allowing one to infer integrated luminosities, dust attenuation parameters, stellar masses, and star formation rates. We tested AGNfitter's performace on real data by fitting the SEDs of a sample of 714 X-ray selected AGN from the XMM-COSMOS survey, spectroscopically classified as Type1 (unobscured) and Type2 (obscured) AGN by their optical-UV emission lines. We find that two independent model parameters, namely the reddening of the accretion disk and the column density of the dusty torus, are good proxies for AGN obscuration, allowing us to develop a strategy for classifying AGN as Type1 or Type2, based solely on an SED-fitting analysis. Our classification scheme is in excellent agreement with the spectroscopic classification, giving a completeness fraction of $\sim 86\%$ and $\sim 70\%$, and an efficiency of $\sim 80\%$ and $\sim 77\%$, for Type1 and Type2 AGNs, respectively.

A Hubble diagram (HD) has recently been constructed in the redshift range 0<z<6.5 using a non-linear relation between the ultraviolet and X-ray luminosities of QSOs. The Type Ia SN HD has already provided a high-precision test of cosmological models, but the fact that the QSO distribution extends well beyond the supernova range (z<1.8), in principle provides us with an important complementary diagnostic whose significantly greater leverage in z can impose tighter constraints on the distance versus redshift relationship. In this paper, we therefore perform an independent test of nine different cosmological models, among which six are expanding, while three are static. Many of these are disfavoured by other kinds of observations (including the aforementioned Type Ia SNe). We wish to examine whether the QSO HD confirms or rejects these earlier conclusions. We find that four of these models (Einstein-de Sitter, the Milne universe, the Static Universe with simple tired light and the Static universe with plasma tired light) are excluded at the >99% C.L. The Quasi-Steady State Model is excluded at >95% C.L. The remaining four models (Lambda-CDM/wCDM, the R_h=ct Universe, the Friedmann open universe and a Static universe with a linear Hubble law) all pass the test. However, only Lambda-CDM/wCDM and $R_{\rm h}=ct$ also pass the Alcock-Paczynski (AP) test. The optimized parameters in Lambda-CDM/wCDM are Omega_m=0.20^{+0.24}_{-0.20} and w_{de}=-1.2^{+1.6}_{-infinity} (the dark-energy equation-of-state). Combined with the AP test, these values become Omega_m=0.38^{+0.20}_{-0.19} and w_{de}=-0.28^{+0.52}_{-0.40}. But whereas this optimization of parameters in Lambda-CDM/wCDM creates some tension with their concordance values, the $R_{\rm h}=ct$ Universe has the advantage of fitting the QSO and AP data without any free parameters.

The observed relation between the soft X-ray and the optical-ultraviolet emission in active galactic nuclei (AGN) is non-linear and it is usually parametrized as a dependence between the logarithm of the monochromatic luminosity at 2500 {\AA} and at 2 keV. Previous investigations have found that the dispersion of this relation is rather high (~0.35-0.4 in log units), which may be caused by measurement uncertainties, variability, and intrinsic dispersion due to differences in the AGN physical properties (e.g. different accretion modes). We show that, once optically-selected quasars with homogeneous SED and X-ray detection are selected, and dust reddened and/or gas obscured objects are not included, the measured dispersion drops to significantly lower values (i.e. ~0.21-0.24 dex). We show that the residual dispersion is due to some extent to variability, and to remaining measurement uncertainties. Therefore, the real physical intrinsic dispersion should be <0.21 dex. Such a tight relation, valid over 4 decades in luminosity, must be the manifestation of an intrinsic (and universal) physical relation between the disk, emitting the primary radiation, and the hot electron corona emitting X-rays.

We present a new method to test the cosmological model, based on the non-linear relation between ultraviolet and X-ray luminosity of quasars. We built a data set of ~1,138 quasars by merging several literature samples with X-ray measurements at 2 keV and SDSS photometry, which was used to estimate the extinction-corrected 2500 Angstrom flux. We obtained three main results: (1) we checked the relation between X-ray and UV luminosities in small redshift bins up to z~6, confirming that it holds at all redshifts with the same slope; (2) we built a Hubble diagram for quasars up to z~6, which is well matched to that of supernovae in the common z=0-1.4 interval, and extends the test of the cosmological model up to z~6; (3) we showed that this relation is a powerful tool to estimate cosmological parameters. Assuming a LambdaCDM model, we obtain Omega_M=0.22+0.08-0.10 and Omega_Lambda=0.92+0.18-0.30 (Omega_M=0.28+-0.04 and Omega_\Lambda=0.73+-0.08 from a joint quasar-SNe fit). However, much more precise measurements will be achieved in the future. A few thousands SDSS quasars already have serendipitous X-ray observations with Chandra or XMM-Newton, and at least 100,000 quasars with UV and X-ray data will be available from the eROSITA all-sky survey in a few years. Euclid, LSST, and Athena surveys will further increase the sample size to at least several hundred thousands. Our simulations show that these samples will provide tight constraints on the cosmological parameters, and will allow to test possible deviations from the standard model with higher precisions than available today.

The aim of this paper is to test the basic model of negative AGN feedback. According to this model, once the central black hole accretes at the Eddington limit and reaches a certain critical mass, AGN driven outflows blow out gas, suppressing star formation in the host galaxy and self-regulating black hole growth. We consider a sample of 224 quasars selected from the SDSS at z<1 observed in the infrared band by Herschel. We evaluate the star formation rate in relation to several outflow signatures traced by the [OIII]4959,5007 and [OII]3726,3729 emission lines in about half of the sample with high quality spectra. Most of the quasars show asymmetric and broad wings in [OIII], which we interpret as outflow signatures. We separate the quasars in two groups, ``weakly'' and ``strongly'' outflowing, using three different criteria. When we compare the mean star formation rate in five redshift bins in the two groups, we find that the SFRs are comparable or slightly larger in the strongly outflowing quasars. We estimate the stellar mass from SED fitting and the quasars are distributed along the star formation main sequence, although with a large scatter. The scatter from this relation is uncorrelated with respect to the kinematic properties of the outflow. Moreover, for quasars dominated in the infrared by starburst or by AGN emission, we do not find any correlation between the star formation rate and the velocity of the outflow, a trend previously reported in the literature for pure starburst galaxies. We conclude that the basic AGN negative feedback scenario seems not to agree with our results. Although we use a large sample of quasars, we did not find any evidence that the star formation rate is suppressed in the presence of AGN driven outflows on large scale. A possibility is that feedback is effective over much longer timescales than those of single episodes of quasar activity.

The ionising continuum from active galactic nuclei (AGN) is fundamental for interpreting their broad emission lines and understanding their impact on the surrounding gas. Furthermore, it provides hints on how matter accretes onto supermassive black holes. Using HST's Wide Field Camera 3 we have constructed the first stacked ultraviolet (rest-frame wavelengths 600-2500\AA) spectrum of 53 luminous quasars at z=2.4, with a state-of-the-art correction for the intervening Lyman forest and Lyman continuum absorption. The continuum slope ($f_\nu \propto \nu^{\alpha_\nu}$) of the full sample shows a break at ~912\AA\ with spectral index $\alpha_\nu=-0.61\pm0.01$ at $\lambda>912$\AA\ and a softening at shorter wavelengths ($\alpha_\nu=-1.70 \pm 0.61$ at $\lambda\leq 912$\AA). Our analysis proves that a proper intergalactic medium absorption correction is required to establish the intrinsic continuum emission of quasars. We interpret our average ultraviolet spectrum in the context of photoionisation, accretion disk models, and quasar contribution to the ultraviolet background. We find that observed broad line ratios are consistent with those predicted assuming an ionising slope of $\alpha_\mathrm{ion}=$-2.0, similar to the observed ionising spectrum in the same wavelength range. The continuum break and softening are consistent with accretion disk plus X-ray corona models when black hole spin is taken into account. Our spectral energy distribution yields a 30% increase to previous estimates of the specific quasar emissivity, such that quasars may contribute significantly to the total specific Lyman limit emissivity estimated from the Ly$\alpha$ forest at z<3.2.

Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host

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.

Aims. The co-evolution of galaxies and super massive black holes (SMBHs) requires that some sort of feedback mechanism is operating during the active galactic nuclei (AGN) phases. AGN driven winds are the most likely candidates for such feedback mechanism, but direct observational evidence of their existence and of their effects on the host galaxies are still scarce and their physical origin is still hotly debated. Methods. X-shooter observations of a sample of X-ray selected, obscured quasars at z$\sim$1.5, selected on the basis of their observed red colors and X-ray-to-optical flux ratio, have shown the presence of outflowing ionized gas identified by broad [OIII] emission lines in 6 out of 8 objects, confirming the efficiency of the selection criteria. Here we present slit-resolved spectroscopy for the two brightest sources, XID2028 and XID5321, to study the complex emission and absorption line kinematics. Results. We detect outflow extended out to $\sim$ 10 kpc from the central black hole (BH), both as blueshifted and redshifted emission. Interestingly, we also detect kpc scale outflows in the [OII] emission lines and in the neutral gas component, traced by the sodium D and magnesium absorption lines, confirming that a substantial amount of the outflowing mass is in the form of neutral gas. Conclusions. The measured gas velocities and the outflow kinetic powers, inferred under reasonable assumptions on the geometry and physical properties of these two systems, favor an AGN origin for the observed winds.

We present X-shooter at the Vewry Large Telescope observations of a sample of 10 luminous, X-ray obscured QSOs at z$\sim1.5$ from the XMM-COSMOS survey, expected to be caught in the transitioning phase from starburst to active galactic nucleus (AGN)-dominated systems. The main selection criterion is X-ray detection at bright fluxes (L$>=10^{44}$ erg s$^{-1}$) coupled to red optical-to-NIR-to-MIR colors. Thanks to its large wavelength coverage, X-shooter allowed us to determine accurate redshifts from the presence of multiple emission lines for five out of six targets for which we had only a photometric redshift estimate, with an 80 percent success rate, significantly larger than what is observed in similar programs of spectroscopic follow-up of red QSOs. We report the detection of broad and shifted components in the [OIII]$\lambda\lambda$5007,4959 complexes for six out of eight sources with these lines observable in regions free from strong atmospheric absorptions. The full width half-maximum (FWHM) associated with the broad components are in the range FWHM$\sim900-1600$ km s$^{-1}$, larger than the average value observed in SDSS Type 2 AGN samples at similar observed [OIII] luminosity, but comparable to those observed for QSO/ultraluminous infrared galaxies systems for which the presence of kpc scale outflows have been revealed through integral field unit spectroscopy. Although the total outflow energetics (inferred under reasonable assumptions) may be consistent with winds accelerated by stellar processes, we favour an AGN origin for the outflows given the high outflow velocities oberved (v$>1000$ km s$^{-1}$) and the presence of strong winds also in objects undetected in the far infrared.

Recoiling supermassive black holes (SMBHs) are considered one plausible physical mechanism to explain high velocity shifts between narrow and broad emission lines sometimes observed in quasar spectra. If the sphere of influence of the recoiling SMBH is such that only the accretion disc is bound, the dusty torus would be left behind, hence the SED should then present distinctive features (i.e. a mid-infrared deficit). Here we present results from fitting the Spectral Energy Distributions (SEDs) of 32 Type-1 AGN with high velocity shifts between broad and narrow lines. The aim is to find peculiar properties in the multi-wavelength SEDs of such objects by comparing their physical parameters (torus and disc luminosity, intrinsic reddening, and size of the 12$\mu$m emitter) with those estimated from a control sample of $\sim1000$ \emph{typical} quasars selected from the Sloan Digital Sky Survey in the same redshift range. We find that all sources, with the possible exception of J1154+0134, analysed here present a significant amount of 12~$\mu$m emission. This is in contrast with a scenario of a SMBH displaced from the center of the galaxy, as expected for an undergoing recoil event.

We present AGNfitter: a Markov Chain Monte Carlo algorithm developed to fit the spectral energy distributions (SEDs) of active galactic nuclei (AGN) with different physical models of AGN components. This code is well suited to determine in a robust way multiple parameters and their uncertainties, which quantify the physical processes responsible for the panchromatic nature of active galaxies and quasars. We describe the technicalities of the code and test its capabilities in the context of X-ray selected obscured AGN using multiwavelength data from the XMM-COSMOS survey.

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.

Quasi-stellar objects (QSOs) occur in galaxies in which supermassive black holes (SMBHs) are growing substantially through rapid accretion of gas. Many popular models of the co-evolutionary growth of galaxies and SMBHs predict that QSOs are also sites of substantial recent star formation, mediated by important processes, such as major mergers, which rapidly transform the nature of galaxies. A detailed study of the star-forming properties of QSOs is a critical test of such models. We present a far-infrared Herschel/PACS study of the mean star formation rate (SFR) of a sample of spectroscopically observed QSOs to z~2 from the COSMOS extragalactic survey. This is the largest sample to date of moderately luminous AGNs studied using uniform, deep far-infrared photometry. We study trends of the mean SFR with redshift, black hole mass, nuclear bolometric luminosity and specific accretion rate (Eddington ratio). To minimize systematics, we have undertaken a uniform determination of SMBH properties, as well as an analysis of important selection effects within spectroscopic QSO samples that influence the interpretation of SFR trends. We find that the mean SFRs of these QSOs are consistent with those of normal massive star-forming galaxies with a fixed scaling between SMBH and galaxy mass at all redshifts. No strong enhancement in SFR is found even among the most rapidly accreting systems, at odds with several co-evolutionary models. Finally, we consider the qualitative effects on mean SFR trends from different assumptions about the star-forming properties of QSO hosts and redshift evolution of the SMBH-galaxy relationship. While limited currently by uncertainties, valuable constraints on AGN-galaxy co-evolution can emerge from our approach.

The fraction of AGN luminosity obscured by dust and re-emitted in the mid-IR is critical for understanding AGN evolution, unification, and parsec-scale AGN physics. For unobscured (Type-1) AGN, where we have a direct view of the accretion disk, the dust covering factor can be measured by computing the ratio of re-processed mid-IR emission to intrinsic nuclear bolometric luminosity. We use this technique to estimate the obscured AGN fraction as a function of luminosity and redshift for 513 Type-1 AGN from the XMM-COSMOS survey. The re-processed and intrinsic luminosities are computed by fitting the 18-band COSMOS photometry with a custom SED-fitting code, which jointly models emission from: hot-dust in the AGN torus, the accretion disk, and the host-galaxy. We find a relatively shallow decrease of the luminosity ratio as a function of Lbol, which we interpret as a corresponding decrease in the obscured fraction. In the context of the receding torus model, where dust sublimation reduces the covering factor of more luminous AGN, our measurements require a torus height which increases with luminosity as h ~ Lbol^{0.3-0.4}. Our obscured fraction-luminosity relation agrees with determinations from SDSS censuses of Type-1 and Type-2 quasars, and favors a torus optically thin to mid-IR radiation. We find a much weaker dependence of obscured fraction on 2-10 keV luminosity than previous determinations from X-ray surveys, and argue that X-ray surveys miss a significant population of highly obscured Compton-thick AGN. Our analysis shows no clear evidence for evolution of obscured fraction with redshift.

Theoretically, bound binaries of massive black holes are expected as the natural outcome of mergers of massive galaxies. From the observational side, however, massive black hole binaries remain elusive. Velocity shifts between narrow and broad emission lines in quasar spectra are considered a promising observational tool to search for spatially unresolved, dynamically bound binaries. In this series of papers we investigate the nature of such candidates through analyses of their spectra, images and multi-wavelength spectral energy distributions. Here we investigate the properties of the optical spectra, including the evolution of the broad line profiles, of all the sources identified in our previous study. We find a diverse phenomenology of broad and narrow line luminosities, widths, shapes, ionization conditions and time variability, which we can broadly ascribe to 4 classes based on the shape of the broad line profiles: 1) Objects with bell-shaped broad lines with big velocity shifts (>1000 km/s) compared to their narrow lines show a variety of broad line widths and luminosities, modest flux variations over a few years, and no significant change in the broad line peak wavelength. 2) Objects with double-peaked broad emission lines tend to show very luminous and broadened lines, and little time variability. 3) Objects with asymmetric broad emission lines show a broad range of broad line luminosities and significant variability of the line profiles. 4) The remaining sources tend to show moderate to low broad line luminosities, and can be ascribed to diverse phenomena. We discuss the implications of our findings in the context of massive black hole binary searches.

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

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.

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

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)

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

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.

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

We present a study of the multi-wavelength properties, from the mid-infrared to the hard X-rays, of a sample of 255 spectroscopically identified X-ray selected Type-2 AGN from the XMM-COSMOS survey. Most of them are obscured the X-ray absorbing column density is determined by either X-ray spectral analyses (for the 45% of the sample), or from hardness ratios. Spectral Energy Distributions (SEDs) are computed for all sources in the sample. The average SEDs in the optical band is dominated by the host-galaxy light, especially at low X-ray luminosities and redshifts. There is also a trend between X-ray and mid-infrared luminosity: the AGN contribution in the infrared is higher at higher X-ray luminosities. We calculate bolometric luminosities, bolometric corrections, stellar masses and star formation rates (SFRs) for these sources using a multi-component modeling to properly disentangle the emission associated to stellar light from that due to black hole accretion. For 90% of the sample we also have the morphological classifications obtained with an upgraded version of the Zurich Estimator of Structural Types (ZEST+). We find that on average Type-2 AGN have lower bolometric corrections than Type-1 AGN. Moreover, we confirm that the morphologies of AGN host-galaxies indicate that there is a preference for these Type-2 AGN to be hosted in bulge-dominated galaxies with stellar masses greater than 10^10 solar masses.

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 a sizable class of type 1 active galactic nuclei (AGNs) with unusually weak near-infrared (1-3{\mu}m) emission in the XMM-COSMOS type 1 AGN sample. The fraction of these "hot-dust-poor" AGNs increases with redshift from 6% at lowredshift (z < 2) to 20% at moderate high redshift (2 < z < 3.5). There is no clear trend of the fraction with other parameters: bolometric luminosity, Eddington ratio, black hole mass, and X-ray luminosity. The 3{\mu}m emission relative to the 1{\mu}m emission is a factor of 2-4 smaller than the typical Elvis et al. AGN spectral energy distribution (SED), which indicates a "torus" covering factor of 2%-29%, a factor of 3-40 smaller than required by unified models. The weak hot dust emission seems to expose an extension of the accretion disk continuum in some of the source SEDs. We estimate the outer edge of their accretion disks to lie at $(0.3-2.0) /times 10^4$ Schwarzschild radii, ~10-23 times the gravitational stability radii. Formation scenarios for these sources are discussed.

We discuss some of the main open issues in the evolution of Active Galactic Nuclei which can be solved by the sensitive, wide area surveys to be performed by the proposed Wide Field X-ray Telescope mission.

By means of a one-zone evolutionary model, we study the co-evolution of supermassive black holes and their host galaxies, as a function of the accretion radiative efficiency, dark matter content, and cosmological infall of gas. In particular, the radiation feedback is computed by using the self-regulated Bondi accretion. The models are characterized by strong oscillations when the galaxy is in the AGN state with a high accretion luminosity. We found that these one-zone models are able to reproduce two important phases of galaxy evolution, namely an obscured-cold phase when the bulk of star formation and black hole accretion occur, and the following quiescent hot phase in which accretion remains highly sub-Eddington. A Compton-thick phase is also found in almost all models, associated with the cold phase. An exploration of the parameter space reveals that the closest agreement with the present-day Magorrian relation is obtained, independently of the dark matter halo mass, for galaxies with a low-mass seed black hole, and the accretion radiative efficiency ~0.1.

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.

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.

(Abriged) We report on the measurement of the rest frame K-band luminosity and total stellar mass of the hosts of 89 broad line Active Galactic Nuclei detected in the zCOSMOS survey in the redshift range 1<z<2.2. The unprecedented multiwavelength coverage of the survey field allows us to disentangle the emission of the host galaxy from that of the nuclear black hole in their Spectral Energy Distributions. We derive an estimate of black hole masses through the analysis of the broad Mg II emission lines observed in the medium-resolution spectra taken with VIMOS/VLT as part of the zCOSMOS project. We found that, as compared to the local value, the average black hole to host galaxy mass ratio appears to evolve positively with redshift, with a best fit evolution of the form (1+z)^{0.68 \pm0.12 +0.6 -0.3}, where the large asymmetric systematic errors stem from the uncertainties in the choice of IMF, in the calibration of the virial relation used to estimate BH masses and in the mean QSO SED adopted. A thorough analysis of observational biases induced by intrinsic scatter in the scaling relations reinforces the conclusion that an evolution of the MBH-M* relation must ensue for actively growing black holes at early times: either its overall normalization, or its intrinsic scatter (or both) appear to increase with redshift. This can be interpreted as signature of either a more rapid growth of supermassive black holes at high redshift, a change of structural properties of AGN hosts at earlier times, or a significant mismatch between the typical growth times of nuclear black holes and host galaxies.