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

Hyperluminous quasars ($L_{\rm bol}\gtrsim 10^{47}$ erg s$^{-1}$) are ideal laboratories to study the interaction and impact of extreme radiative field and the most powerful winds in the AGN nuclear regions. They typically exhibit low coronal X-ray luminosity ($L_{\rm X}$) compared to the UV and MIR radiative outputs ($L_{\rm UV}$ and $L_{\rm MIR}$) with a non-negligible fraction of them reporting even $\sim$1 dex weaker $L_{\rm X}$ compared to the prediction of the well established $L_{\rm X}$-$L_{\rm UV}$ and $L_{\rm X}$-$L_{\rm MIR}$ relations followed by the bulk of the AGN population. We report in our WISE/SDSS-selected Hyperluminous (WISSH) $z=2-4$ broad-line quasar sample, the discovery of a dependence between the intrinsic 2-10 keV luminosity ($L_{\rm 2-10}$) and the blueshifted velocity of the CIV emission line ($v_{\rm CIV}$) indicative of accretion disc winds. In particular, sources with fastest winds ($v_{\rm CIV}\gtrsim 3000~\rm km s^{-1}$) possess $\sim$0.5-1 dex lower $L_{\rm 2-10}$ than sources with negligible $v_{\rm CIV}$. No similar dependence is found on $L_{\rm UV}$, $L_{\rm MIR}$, $L_{\rm bol}$, photon index and absorption column density. We interpret these findings in the context of accretion disc wind models. Both magnetohydrodynamic and line-driven models can qualitatively explain the reported relations as a consequence of X-ray shielding from the inner wind regions. In case of line-driven winds, the launch of fast winds is favoured by a reduced X-ray emission, and we speculate that these winds may play a role in directly limiting the coronal hard X-ray production.

During the last years, Ly$\alpha$ nebulae have been routinely detected around high-z, radio-quiet quasars (RQQs) thanks to the advent of sensitive integral field spectrographs. Constraining the physical properties of the Ly$\alpha$ nebulae is crucial for a full understanding of the circum-galactic medium (CGM), which is a venue of feeding and feedback processes. The most luminous quasars are privileged test-beds to study these processes, given their large ionizing fluxes and dense CGM environments in which they are expected to be embedded. We aim at characterizing the rest-frame UV emission lines in the CGM around a hyper-luminous, broad emission line, RQQ at z~3.6, that exhibits powerful outflows at both nuclear and host galaxy scales. We analyze VLT/MUSE observations of the quasar J1538+08 and perform a search for extended UV emission lines to characterize its morphology, emissivity, kinematics and metal content. We report the discovery of a very luminous ($\sim2 \times10^{44}~erg~s^{-1}$), giant Ly$\alpha$ nebula and a likely associated extended CIV nebula. The Ly$\alpha$ nebula emission exhibits moderate blueshift compared with the quasar systemic redshift and large average velocity dispersion ($\sigma_{v}$ ~700 $km~s^{-1}$) across the nebula, while the CIV nebula shows $\sigma_{v}$~$350~km~s^{-1}$. The Ly$\alpha$ line profile exhibits a significant asymmetry towards negative velocity values at 20-30 kpc south of the quasar and is well parameterized by two Gaussian components: a narrow ($\sigma$~$470~km~s^{-1}$) systemic one plus a broad ($\sigma$~1200 $km~s^{-1}$), blueshifted (~1500 $km~s^{-1}$) one. Our analysis of the MUSE observation of J1538+08 reveals metal-enriched CGM around this hyper-luminous quasar and our detection of blueshifted emission in the line profile of the Ly$\alpha$ nebula suggests that powerful nuclear outflows can propagate through the CGM over tens of kpc.

We present Atacama Large Millimiter/submillimiter Array (ALMA) observations of eight highly excited CO (J$_{\rm up}>8$) lines and continuum emission in two $z\sim6$ quasars: SDSS J231038.88+185519.7 (hereafter J2310), for which CO(8-7), CO(9-8), and CO(17-16) lines have been observed, and ULAS J131911.29+095951.4 (J1319), observed in the CO(14-13), CO(17-16) and CO(19-18) lines. The continuum emission of both quasars arises from a compact region ($< 0.9$ kpc). By assuming a modified black-body law, we estimate dust masses of Log$(M_{\rm dust}/M_{\odot})=8.75\pm0.07$ and Log$(M_{\rm dust}/M_{\odot})=8.8\pm0.2$ and dust temperatures of $T_{\rm dust}=76\pm3~{\rm K}$ and $T_{\rm dust}=66^{+15}_{-10}~{\rm K}$, respectively for J2310 and J1319. Only CO(8-7) and CO(9-8) in J2310 are detected, while $3\sigma$ upper limits on luminosities are reported for the other lines of both quasars. The CO line luminosities and upper limits measured in J2310 and J1319 are consistent with those observed in local AGN and starburst galaxies, and other $z\sim 6$ quasars, except for SDSS J1148+5251 (J1148), the only quasar at $z=6.4$ with a previous CO(17-16) line detection. By computing the CO SLEDs normalised to the CO(6-5) line and FIR luminosities for J2310, J1319, and J1149, we conclude that different gas heating mechanisms (X-ray radiation and/or shocks) may explain the different CO luminosities observed in these $z\sim6$ quasar. Future J$_{\rm up}>8$ CO observations will be crucial to understand the processes responsible for molecular gas excitation in luminous high-$z$ quasars.

We present interferometric observations of the CN(1-0) line emission in Mrk231 and combine them with previous observations of CO and other H$_2$ gas tracers to study the physical properties of the massive molecular outflow. We find a strong boost of the CN/CO(1-0) line luminosity ratio in the outflow, which is unprecedented compared to any other known Galactic or extragalactic source. For the dense gas phase in the outflow traced by the HCN and CN emissions, we infer $\rm X_{\rm CN}\equiv [CN]/[H_2] > X_{\rm HCN}$ by at least a factor of three, with H$_2$ gas densities of $n_{\rm H_2}\sim10^{5-6}$ cm$^{-3}$. In addition, for the first time, we resolve narrow spectral features in the HCN(1-0) and HCO$^+$(1-0) high-velocity line wings tracing the dense phase of the outflow. The velocity dispersions of these spectral features, $\sigma_v\sim7-20$ km s$^{-1}$, are consistent with those of massive extragalactic giant molecular clouds detected in nearby starburst nuclei. The H$_2$ gas masses inferred from the HCN data are quite high, $M_{mol}\sim0.3-5\times10^8$ $M_{\odot}$. Our results suggest that massive, denser molecular gas complexes survive embedded into the more diffuse H$_2$ phase of the outflow, and that the chemistry of such outflowing dense clouds is affected by enhanced UV radiation.

We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged)

The WISSH quasars project aims at studying the nuclear and host galaxy properties of the most luminous quasars ($L_{bol}>10^{47}$ erg/s, $1.8<z<4.6$). Nuclear winds are manifested as UV broad ($\geq$ 2,000 km/s) absorption lines (BAL) in $\sim$ 15\% of quasars. We aim at studying the incidence and properties of such winds in the WISSH sample, to investigate possible differences with respect to lower luminosity AGN regimes. We collected optical spectra from the SDSS data release 12, and identified those showing absorption troughs in the region between \siiv and \civ emission lines. We find a higher observed fraction of \civ BAL quasars in the WISSH sample (24\%) with respect to previous catalogues (10-15\%). These WISSH BAL quasars are also characterized by a larger average BI ($\sim$4,000 km/s) and maximum velocity ($\sim$17,000 km/s). Moreover, for two objects we discovered BAL features bluewards of the \siiv peak, which can be associated to \civ absorption with velocity of 0.15c. Finally, we updated previous studies on the dependence of maximum outflow velocity upon bolometric luminosity, showing that BAL winds have intermediate properties compared to molecular/ionized winds and ultra fast outflows (UFOs). Finally, the radio properties of the WISSH BAL quasars as a whole are in line with those of samples at lower luminosities from previous studies. Our results suggest that the higher $L_{\rm bol}$ of the WISSH quasars likely favours the acceleration of BAL outflows and that their most likely driving mechanism is radiation pressure. Furthermore we estimate that the kinetic power associated to these winds in hyperluminous quasars is sufficient, for the highest column density and fastest winds, to provide efficient feedback onto the host galaxy.

Starburst galaxies and star-forming active galactic nuclei (AGN) are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hours of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very high energy band. We did not detect significant gamma-ray emission, and set upper limits at 95\% confidence level to the gamma-ray flux above 200 GeV f<5.1x10^{-13} cm^{-2} s ^{-1} . This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr^{-1}.

We report on the first ALMA observation of the CO(3$-$2) and rest-frame ~340 GHz continuum emission in PDS 456, which is the most luminous, radio-quiet QSO in the local Universe ($z$~0.18), with a bolometric luminosity $L_{\rm Bol}\sim10^{47}$ erg s$^{-1}$. ALMA angular resolution allowed us to map scales as small as ~700 pc. The molecular gas reservoir, traced by the core of the very bright CO(3$-$2) emission line, is distributed in a compact rotating disk, with size of ~1.3 kpc, seen close to face-on ($i$~25 deg). Fast CO(3$-$2) emission in the velocity range $v\in[-1000,+500]$ km s$^{-1}$ is also present. Specifically, we detect several blue-shifted clumps out to ~5 kpc from the nucleus, in addition to a compact ($R\lesssim1.2$ kpc), broad emission component. These components reveal a galaxy-wide molecular outflow, with a total mass $M_{\rm mol}^{\rm out}\sim2.5\times10^8$ $M_{\odot}$ and a mass outflow rate $\dot{M}_{\rm mol}\sim290$ $M_{\odot}$ yr$^{-1}$. The corresponding depletion time is ~8 Myr, shorter than the rate at which the molecular gas is converted into stars, indicating that the detected outflow is potentially able to quench star-formation in the host. The momentum flux of the molecular outflow normalised to the radiative momentum output (i.e. $L_{\rm Bol}/c$) is $\lesssim1$, comparable to that of the X-ray ultra-fast outflow (UFO) detected in PDS 456. This is at odds with the expectations for an energy-conserving expansion suggested for most of the large-scale outflows detected in low-luminosity AGN so far. We suggest three possible scenarios that may explain this observation: (i) in very luminous AGN such as our target the molecular gas phase is tracing only a fraction of the total outflowing mass; (ii) a small coupling between the shocked gas by the UFO and the host-galaxy ISM (iii) AGN radiation pressure may play an important role in driving the outflow.

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

We present ALMA rest-frame 230 GHz continuum and CO(2-1) line observations of the nearby Compton-thick Seyfert galaxy ESO428-G14, with angular resolution 0.7 arcsec (78 pc). We detect CO(2-1) emission from spiral arms and a circum-nuclear ring with 200 pc radius, and from a transverse gas lane with size of $\sim100$ pc, which crosses the nucleus and connects the two portions the circumnuclear ring. The molecular gas in the host galaxy is distributed in a rotating disk with intrinsic circular velocity $v_{rot}=135$ km/s, inclination $i=57$ deg, and dynamical mass $M_{dyn }=5\times 10^9~\rm M_{\odot}$ within a radius of $\sim 1$ kpc. In the inner 100 pc region CO is distributed in a equatorial bar, whose kinematics is highly perturbed and consistent with an inflow of gas towards the AGN. This inner CO bar overlaps with the most obscured, Compton-thick region seen in X-rays. We derive a column density of $\rm N(H_2) \approx 2\times10^{23}~ cm^{-2}$ in this region, suggesting that molecular gas may contribute significantly to the AGN obscuration. We detect a molecular outflow with a total outflow rate $\rm \dot M_{of}\approx 0.8~M_{\odot}/yr$, distributed along a bi-conical structure with size of $700$ pc on both sides of the AGN. The bi-conical outflow is also detected in the $\rm H_2$ emission line at 2.12 $\mu$m, which traces a warmer nuclear outflow located within 170 pc from the AGN. This suggests that the outflow cools with increasing distance from the AGN. We find that the hard X-ray emitting nuclear region mapped with Chandra is CO-deprived, but filled with warm molecular gas traced by $\rm H_2$ - thus confirming that the hard (3-6 keV) continuum and Fe K$\alpha$ emission are due to scattering from dense neutral clouds in the ISM.

Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.

We investigated the interstellar medium properties of the disc and outflowing gas in the central regions of nine nearby Seyfert galaxies, all characterised by prominent outflows. These objects are part of the Measuring Active Galactic Nuclei Under MUSE Microscope survey, which aims to probe their physical conditions and ionisation mechanism by exploiting MUSE unprecedented sensitivity. We studied the different properties of the gas in the disc and outflow with spatially and kinematically resolved maps by dividing the strongest emission lines in velocity bins. We associated the core of the lines with the disc, consistent with the stellar velocity, and the redshifted and the blueshifted wings with the outflow. We find that the outflowing gas is characterised by higher values of density and ionisation parameter than the disc, which presents a higher dust extinction. Moreover, we distinguish high- and low-ionisation regions across the portion of spatially resolved narrow-line region traced by the outflowing gas. The high-ionisation regions characterised by the lowest [NII]/H{\alpha} and [SII]/H{\alpha} line ratios generally trace the innermost parts along the axis of the emitting cones where the [SIII]/[SII] line ratio is enhanced, while the low-ionisation regions follow the cone edges and/or the regions perpendicular to the axis of the outflows, also characterised by a higher [OIII] velocity dispersion. A possible scenario to explain these features relies on the presence of two distinct populations of line emitting clouds: one is optically thin to the radiation and is characterised by the highest excitation, while the other is optically thick and is impinged by a filtered, and thus harder, radiation field which generates strong low-excitation lines. The highest values of [NII]/H{\alpha} and [SII]/H{\alpha} line ratios may be due to shocks and/or a hard filtered radiation field from the AGN.

We report the discovery of Serenity-18, a galaxy at z=5.939 for which we could measure the content of molecular gas, M(H_2)~ 5 x10^9 M_sun, traced by the CO(6-5) emission, together with the metal-poor ([Fe/H]=-3.08 +- 0.12, [Si/H]=-2.86 +- 0.14) gas clump/filament which is possibly feeding its growth. The galaxy has an estimated star formation rate of ~100 M_sun yr^{-1}, implying that it is a typical main sequence galaxy at these redshifts. The metal-poor gas is detected through a damped Lyman-alpha absorber (DLA) observed at a spatial separation of 40 kpc and at the same redshift of Serenity-18, along the line of sight to the quasar SDSS J2310+1855 (z_em = 6.0025). The chemical abundances measured for the damped Lyman-alpha system are in very good agreement with those measured for other DLAs discovered at similar redshifts, indicating an enrichment due to massive PopII stars. The galaxy/Damped system we discovered is a direct observational evidence of the assembly of a galaxy at the edge of the reionization epoch.

We report on the coexistence of powerful gas outflows observed in millimeter and X-ray data of the Radio-Loud Narrow Line Seyfert 1 Galaxy IRAS 17020+4544. Thanks to the large collecting power of the Large Millimeter Telescope, a prominent line arising from the 12CO(1-0) transition was revealed in recent observations of this source. The complex profile is composed by a narrow double-peak line and a broad wing. While the double-peak structure may be arising in a disk of molecular material, the broad wing is interpreted as the signature of a massive outflow of molecular gas with an approximate bulk velocity of -660 km/s. This molecular wind is likely associated to a multi-component X-ray Ultra-Fast Outflow with velocities reaching up to ~0.1c and column densities in the range 10^{21-23.9} cm^-2 that was reported in the source prior to the LMT observations. The momentum load estimated in the two gas phases indicates that within the observational uncertainties the outflow is consistent with being propagating through the galaxy and sweeping up the gas while conserving its energy. This scenario, which has been often postulated as a viable mechanism of how AGN feedback takes place, has so far been observed only in ULIRGs sources. IRAS 17020+4544 with bolometric and infrared luminosity respectively of 5X10^{44} erg/s and 1.05X10^{11} L_sun appears to be an example of AGN feedback in a NLSy1 Galaxy (a low power AGN). New proprietary multi-wavelength data recently obtained on this source will allow us to corroborate the proposed hypothesis.

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1<z<3) has been poorly explored so far, and current observations in this redshift range are mostly limited to targets with high chances to be in an outflowing phase. This paper introduces SUPER (a SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback), an ongoing ESO's VLT/SINFONI Large Programme. SUPER will perform the first systematic investigation of ionized outflows in a sizeable and blindly-selected sample of 39 X-ray AGN at z~2, which reaches high spatial resolutions (~2 kpc) thanks to the adaptive optics-assisted IFU observations. The outflow morphology and star formation in the host galaxy will be mapped through the broad component of [OIII] and the narrow component of Ha emission lines. The main aim of our survey is to infer the impact of outflows on the on-going star formation and to link the outflow properties to a number of AGN and host galaxy properties. We describe here the survey characteristics and goals, as well as the selection of the target sample. Moreover, we present a full characterization of its multi-wavelength properties: we measure, via spectral energy distribution fitting of UV-to-FIR photometry, stellar masses (4x10^9-2x10^11 Msun), star formation rates (25-680 Msun yr^-1) and AGN bolometric luminosities (2x10^44-8x10^47 erg s^-1), along with obscuring column densities (<2x10^24 cm^-2) and 2-10 keV luminosities (2x10^43-6x10^45 erg s^-1) derived through X-ray spectral analysis. Finally, we classify our AGN as jetted or non-jetted according to their radio and FIR emission.

Ionized outflows, revealed by broad asymmetric wings of the [OIII] line, are commonly observed in AGN but the low intrinsic spatial resolution of observations has generally prevented a detailed characterization of their properties. The MAGNUM survey aims at overcoming these limitations by focusing on the nearest AGN, including NGC 1365, a nearby Seyfert galaxy (D~17 Mpc), hosting a low-luminosity AGN (Lbol ~ 2x10^43 erg/s). We want to obtain a detailed picture of the ionized gas in the central ~5 kpc of NGC 1365 in terms of physical properties, kinematics, and ionization mechanisms. We also aim to characterize the warm ionized outflow as a function of distance from the nucleus and its relation with the nuclear X-ray wind. We employed VLT/MUSE optical integral field spectroscopic observations to investigate the warm ionized gas and Chandra ACIS-S X-ray data for the hot highly-ionized phase. We obtained flux, kinematic, and diagnostic maps of the optical emission lines, which we used to disentangle outflows from disk motions and measure the gas properties down to a spatial resolution of ~70 pc. [OIII] emission mostly traces an AGN-ionized kpc-scale biconical outflow with velocities up to ~200 km/s. H{\alpha} emission traces instead star formation in a circumnuclear ring and along the bar, where we detect non-circular motions. Soft X-rays are mostly due to thermal emission from the star-forming regions, but we could isolate the AGN photoionized component which matches the [OIII] emission. The mass outflow rate of the extended ionized outflow matches that of the nuclear X-ray wind and then decreases with radius. However, the hard X-ray emission from the circumnuclear ring suggests that star formation might contribute to the outflow. The integrated mass outflow rate, kinetic energy rate, and outflow velocity are broadly consistent with the typical relations observed in more luminous AGN.

We present ALMA observations of the CO(6-5) and [CII] emission lines and the sub-millimeter continuum of the $z\sim6$ quasi-stellar object (QSO) SDSS J231038.88+185519.7. Compared to previous studies, we have analyzed a synthetic beam that is ten times smaller in angular size, we have achieved ten times better sensitivity in the CO(6-5) line, and two and half times better sensitivity in the [CII] line, enabling us to resolve the molecular gas emission. We obtain a size of the dense molecular gas of $2.9\pm0.5$ kpc, and of $1.4\pm0.2$ kpc for the 91.5 GHz dust continuum. By assuming that CO(6-5) is thermalized, and by adopting a CO--to--$H_2$ conversion factor $\rm \alpha_{CO} = 0.8~ M_{\odot}~K^{-1}~ (km/s)^{-1} ~pc^{2}$, we infer a molecular gas mass of $\rm M(H_2)=(3.2 \pm0.2) \times 10^{10}\rm M_{\odot}$. Assuming that the observed CO velocity gradient is due to an inclined rotating disk, we derive a dynamical mass of $\rm M_{dyn}~sin^2(i) = (2.4\pm0.5) \times 10^{10}~ M_{\odot}$, which is a factor of approximately two smaller than the previously reported estimate based on [CII]. Regarding the central black hole, we provide a new estimate of the black hole mass based on the C~IV emission line detected in the X-SHOOTER/VLT spectrum: $\rm M_{BH}=(1.8\pm 0.5) \times 10^{9}~ M_{\odot}$. We find a molecular gas fraction of $\rm \mu=M(H_2)/M^*\sim4.4$, where $\rm M^*\approx M_{dyn} - M(H_2)-M(BH)$. We derive a ratio $v_{rot}/\sigma \approx 1-2$ suggesting high gas turbulence, outflows/inflows and/or complex kinematics due to a merger event. We estimate a global Toomre parameter $Q\sim 0.2-0.5$, indicating likely cloud fragmentation. We compare, at the same angular resolution, the CO(6-5) and [CII] distributions, finding that dense molecular gas is more centrally concentrated with respect to [CII]. We find that the current BH growth rate is similar to that of its host galaxy.

Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming, disk galaxies and quiescent spheroidal systems. Yet, our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work we use new deep, high-spatial resolution ALMA CO(3-2) and archival VLT/SINFONI H$\alpha$ observations to study the molecular and ionized components of the AGN-driven outflow in zC400528 ---a massive, main sequence galaxy at z=2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient and extends for at least ~10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale at least twice shorter than that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star-formation activity in zC400528.

We present a 155ks NuSTAR observation of the $z\sim2$ hot dust-obscured galaxy (hot DOG) W1835+4355. We extracted spectra from the two NuSTAR detectors and analyzed them jointly with the archival XMM PN and MOS spectra. We performed a spectroscopic analysis based on both phenomenological and physically motivated models employing toroidal and spherical geometry for the obscurer. In all the modelings, the source exhibits a Compton-thick column density $N_{\rm H} \gtrsim 10^{24}$ cm$^{-2}$, a 2-10 keV luminosity $L_{2-10}\approx2\times10^{45}$ erg s$^{-1}$ , and a prominent soft excess ($\sim5-10$ % of the primary radiative output), which translates into a luminosity $\sim10^{44}$ erg s$^{-1}$. We modeled the spectral energy distribution from 1.6 to 850 $\mu m$ using a clumpy two-phase dusty torus model plus a modified blackbody to account for emission powered by star formation in the far-infrared. We employed several geometrical configurations consistent with those applied in the X-ray analysis. In all cases we obtained a bolometric luminosity $L_{\rm bol}\approx3-5\times10^{47}$ erg s$^{-1}$, which confirms the hyperluminous nature of this active galactic nucleus. Finally, we estimate a prodigious star formation rate of $\sim$3000 $M_{\odot}\,yr^{-1}$, which is consistent with the rates inferred for $z\approx2-4$ hyperluminous type I quasars. The heavily obscured nature, together with $L_{\rm bol}$, the ratio of X-ray to mid-infrared luminosity, the rest-frame optical morphology, and the host star formation rate are indicative of its evolutionary stage. We can interpret this as a late-stage merger event in the transitional, dust-enshrouded, evolutionary phase eventually leading to an optically bright AGN.

We present ALMA and ACA [CI]$^{3}P_{1}-^{3}P_{0}$ ([CI](1-0)) observations of NGC6240, which we combine with ALMA CO(2-1) and IRAM PdBI CO(1-0) data to study the physical properties of the massive molecular (H$_2$) outflow. We discover that the receding and approaching sides of the H$_2$ outflow, aligned east-west, exceed 10 kpc in their total extent. High resolution ($0.24"$) [CI](1-0) line images surprisingly reveal that the outflow emission peaks between the two AGNs, rather than on either of the two, and that it dominates the velocity field in this nuclear region. We combine the [CI](1-0) and CO(1-0) data to constrain the CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$) in the outflow, which is on average $2.1\pm1.2~\rm M_{\odot} (K~km~s^{-1}~pc^2)^{-1}$. We estimate that $60\pm20$ % of the total H$_2$ gas reservoir of NGC6240 is entrained in the outflow, for a resulting mass-loss rate of $\dot{M}_{\rm out}=2500\pm1200~M_{\odot}~yr^{-1}$ $\equiv 50\pm30$ SFR. This energetics rules out a solely star formation-driven wind, but the puzzling morphology challenges a classic radiative-mode AGN feedback scenario. For the quiescent gas we compute $\langle\alpha_{\rm CO}\rangle = 3.2\pm1.8~\rm M_{\odot} (K~km~s^{-1}~pc^2)^{-1}$, which is at least twice the value commonly employed for (U)LIRGs. We observe a tentative trend of increasing $r_{21}\equiv L^{\prime}_{\rm CO(2-1)}/L^{\prime}_{\rm CO(1-0)}$ ratios with velocity dispersion and measure $r_{21}>1$ in the outflow, whereas $r_{21}\simeq1$ in the quiescent gas. We propose that molecular outflows are the location of the warmer, strongly unbound phase that partially reduces the opacity of the CO lines in (U)LIRGs, hence driving down their global $\alpha_{\rm CO}$ and increasing their $r_{21}$ values.

We present the broad-band soft and hard X-ray spectral analysis of 8 Narrow Line Seyfert 1 galaxies extracted from the IBISCO Sample. The study also focuses on the properties of the NLS1 in our sample in relation to those of the IBISCO parent Seyfert population. The IBISCO sample comprises 57 AGN selected from the INTEGRAL IBIS AGN catalogue (in the 20-100 keV band), with z$<$0.05 and covering a wide range of luminosities, BH masses and absorption. All AGN have also measurements of the molecular gas (H$_2$) content of their host galaxies, through the detection of CO emission lines. The main goals of this analysis are to accurately determine the X-ray continuum emission, investigate the presence of absorption features around 7 keV (indicative of the presence of outflows) and measure the bolometric luminosity in order to study the accretion parameters of the eight IBISCO NLS1, and study the accretion mechanisms and investigate the feeding and feedback cycle in these peculiar AGN. Preliminary results show that NLS1 tend to have higher Eddington ratios and larger molecular gas fractions than their parent Seyfert population in the IBISCO sample. Nuclear (AGN) vs. host galaxy properties scaling relations of NLS1 in relation to the parent Seyfert population are also explored.

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 an ALMA high-resolution observation of the 840 um continuum and [CII] line emission in the WISE-SDSS selected hyper-luminous (WISSH) QSO J1015+0020 at z~4.4. Our analysis reveals an exceptional overdensity of [CII]-emitting companions with a very small (<150 km/s) velocity shift with respect to the QSO redshift. We report the discovery of the closest companion observed so far in submillimetre observations of high-z QSOs. It is only 2.2 kpc distant and merging with J1015+0020, while two other [CII] emitters are found at 8 and 17 kpc. Two strong continuum emitters are also detected at <3.5 arcsec. They are likely associated to the same overdense structure of J1015+0020, as they exceed by a factor of 100 the number of expected sources, considering the Log(N)-Log(S) at 850 um. The host galaxy of J1015+0020 shows a SFR of about 100 Msun/yr while the total SFR of the QSO and its companion galaxies is a factor of 10 higher, indicating that substantial stellar mass assembly at early epochs may have taken place in the QSO satellites. For J1015+0020 we compute a SMBH mass MBH~6E9 Msun and a dynamical mass Mdyn~4E10 Msun . This translates into an extreme ratio Mdyn/MBH~7. The total stellar mass of the QSO host galaxy plus the [CII] emitters already exceeds 1E11 Msun at z~4.4. These sources will likely merge and develop into a giant galaxy of 1.3E12 Msun. Under the assumption of constant mass or Eddington accretion rate equal to the observed values, we find that the growth timescale of the host galaxy is comparable or even shorter than that inferred for the SMBH.

We present the analysis of the rest-frame optical-to-UV spectrum of APM 08279+5255, a well-known lensed broad absorption line (BAL) quasar at $z = 3.911$. The spectroscopic data are taken with the optical DOLoRes and near-IR NICS instruments at TNG, and include the previously unexplored range between C III] $\lambda$1910 and [O III] $\lambda\lambda$4959,5007. We investigate the possible presence of multiple BALs by computing "balnicity" and absorption indexes (i.e. BI, BI$_0$ and AI) for the transitions Si IV $\lambda$1400, C IV $\lambda$1549, Al III $\lambda$1860 and Mg II $\lambda$2800. No clear evidence for the presence of absorption features is found in addition to the already known, prominent BAL associated to C IV, which supports a high-ionization BAL classification for APM 08279+5255. We also study the properties of the [O III], H$\beta$ and Mg II emission lines. We find that [O III] is intrinsically weak ($F_{\rm [OIII]}/F_{\rm H\beta} \lesssim 0.04$), as it is typically found in luminous quasars with a strongly blueshifted C IV emission line ($\sim$2500 km s$^{-1}$ for APM 08279+5255). We compute the single-epoch black hole mass based on Mg II and H$\beta$ broad emission lines, finding $M_{\rm BH} = (2 \div 3) \times 10^{10}\mu^{-1}$ M$_\odot$, with the magnification factor $\mu$ that can vary between 4 and 100 according to CO and rest-frame UV-to-mid-IR imaging respectively. Using a Mg II equivalent width (EW)-to-Eddington ratio relation, the EW$_{\rm MgII} \sim 27$ \AA\ measured for APM 08279+5255 translates into an Eddington ratio of $\sim$0.4, which is more consistent with $\mu=4$. This magnification factor also provides a value of $M_{\rm BH}$ that is consistent with recent reverberation-mapping measurements derived from C IV and Si IV.

We have undertaken a multi-band observing program aimed at obtaining a complete census of winds in a sample of WISE/SDSS selected hyper-luminous (WISSH) QSOs at z~2-4. We have analyzed the rest-frame optical (LBT/LUCI and VLT/SINFONI) and UV (SDSS) spectra of 18 randomly selected WISSH QSOs to measure the SMBH mass and study the properties of winds both in the NLR and BLR traced by blueshifted/skewed [OIII] and CIV emission lines, respectively. These WISSH QSOs are powered by SMBH with masses $\ge$10$^9$ Msun accreting at 0.4<$\lambda_{Edd}$<3.1. We have found the existence of two sub-populations characterized by the presence of outflows at different distances from the SMBH. One population ([OIII] sources) exhibits powerful [OIII] outflows, rest-frame EW (REW) of the CIV emission REW$_{CIV}\approx$20-40 A and modest CIV velocity shift (v$_{CIV}^{peak}$) with respect to the systemic redshift (<=2000 km/s). The second population (Weak [OIII] sources), representing ~70% of the analyzed WISSH QSOs, shows weak/absent [OIII] emission and an extremely large v$_{CIV}^{peak}$ (up to ~8000 km/s and REW$_{CIV}$<=20 A). We propose two explanations for the observed behavior of the strength of the [OIII] emission in terms of orientation effects of the line of sight and ionization cone. The dichotomy in the presence of BLR and NLR winds could be likely due to inclination effects considering a polar geometry scenario for the BLR winds. We find a strong correlation with L$_{Bol}$ and an anti-correlation with $\alpha_{ox}$, whereby the higher L$_{Bol}$, the steeper $\alpha_{ox}$ and the larger is the v$_{CIV}^{peak}$. Finally, the observed dependence v$_{CIV}^{peak}\propto L_{Bol}^{0.28\pm0.04}$ is consistent with radiatively driven winds scenario, where strong UV continuum is necessary to launch the wind and a weakness of the X-ray emission is fundamental to prevent overionization of the wind itself.

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.

We report the detection of the CO(3-2) emission line from the host galaxy of Gamma Ray Burst (GRB) 080207 at $z$ = 2.086. This is the first detection of molecular gas in emission from a GRB host galaxy beyond redshift 1. We find this galaxy to be rich in molecular gas with a mass of $1.1 \times 10^{11}\,\rm M_{\odot}$ assuming $\alpha_{\rm CO}=$ 4.36 $\rm M_{\odot}(\rm K\,km\,s^{-1}\,pc^2)^{-1}$. The molecular gas mass fraction of the galaxy is $\sim$ 0.5, typical of star forming galaxies (SFGs) with similar stellar masses and redshifts. With a $\rm SFR_{FIR}$ of 260 $\rm M_{\odot}\,yr^{-1}$, we measure a molecular-gas-depletion timescale of 0.43 Gyr, near the peak of the depletion timescale distribution of SFGs at similar redshifts. Our findings are therefore in contradiction with the proposed molecular gas deficiency in GRB host galaxies. We argue that the reported molecular gas deficiency for GRB hosts could be the artifact of improper comparisons or neglecting the effect of the typical low metallicities of GRB hosts on the CO-to-molecular-gas conversion factor. We also compare the kinematics of the CO(3-2) emission line to that of the H$\alpha$ emission line from the host galaxy. We find the H$\alpha$ emission to have contributions from two separate components, a narrow and a broad one. The narrow component matches the CO emission well in velocity space. The broad component, with a FWHM of $\sim$ 1100 $\rm km\,s^{-1}$, is separated by $+390$ $\rm km\,s^{-1}$ in velocity space from the narrow component. We speculate this broad component to be associated with a powerful outflow in the host galaxy or in an interacting system.

This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

We present the first attempt to detect outflows from galaxies approaching the Epoch of Reionization (EoR) using a sample of 9 star-forming ($\rm SFR=31\pm 20~M_{\odot}~yr^{-1}$) $z\sim 5.5$ galaxies for which the [CII]158$\mu$m line has been previously obtained with ALMA. We first fit each line with a Gaussian function and compute the residuals by subtracting the best fitting model from the data. We combine the residuals of all sample galaxies and find that the total signal is characterised by a flux excess of $\sim 0.5$ mJy extended over $\sim 1000$ km~s$^{-1}$. Although we cannot exclude that part of this signal is due to emission from faint satellite galaxies, we show that the most probable explanation for the detected flux excess is the presence of broad wings in the [CII] lines, signatures of starburst-driven outflows. We infer an average outflow rate of $\rm \dot{M}=54\pm23~ M_{\odot}~yr^{-1}$, providing a loading factor $\eta=\rm \dot{M}/SFR=1.7\pm1.3$ in agreement with observed local starbursts. Our interpretation is consistent with outcomes from zoomed hydro-simulations of {\it Dahlia}, a $z\sim 6$ galaxy ($\rm SFR\sim 100~\rm M_{\odot}~yr^{-1}$) whose feedback-regulated star formation results into an outflow rate $\rm \dot{M}\sim 30~ M_{\odot}~yr^{-1}$. The quality of the ALMA data is not sufficient for a detailed analysis of the [CII] line profile in individual galaxies. Nevertheless, our results suggest that starburst-driven outflows are in place in the EoR and provide useful indications for future ALMA campaigns. Deeper observations of the [CII] line in this sample are required to better characterise feedback at high-$z$ and to understand the role of outflows in shaping early galaxy formation.

Galactic winds driven by active galactic nuclei (AGN) have been invoked to play a fundamental role in the co-evolution between supermassive black holes and their host galaxies. Finding observational evidence of such feedback mechanisms is of crucial importance and it requires a multi-wavelength approach in order to compare winds at different scales and phases. In Tombesi et al. (2015) we reported the detection of a powerful ultra-fast outflow (UFO) in the Suzaku X-ray spectrum of the ultra-luminous infrared galaxy IRAS F11119$+$3257. The comparison with a galaxy-scale OH molecular outflow observed with Herschel in the same source supported the energy-conserving scenario for AGN feedback. The main objective of this work is to perform an independent check of the Suzaku results using the higher sensitivity and wider X-ray continuum coverage of NuSTAR. We clearly detect a highly ionized Fe K UFO in the 100ks NuSTAR spectrum with parameters $N_H = (3.2\pm1.5)\times 10^{24}$ cm$^{-2}$, log$\xi$$=$$4.0^{+1.2}_{-0.3}$ erg s$^{-1}$ cm, and $v_{\text{out}} = 0.253^{+0.061}_{-0.118}$ c. The launching radius is likely at a distance of $r \ge 16$ $r_s$ from the black hole. The mass outflow rate is in the range $\dot{M}_{out}$$\simeq$0.5-2 $M_{\odot}$ yr$^{-1}$. The UFO momentum rate and power are $\dot{P}_{out} \simeq$0.5-2 $L_{AGN}/c$ and $\dot{E}_{out} \simeq$7-27\% $L_{AGN}$, respectively. The UFO parameters are consistent between the 2013 Suzaku and the 2015 NuSTAR observations. Only the column density is found to be variable, possibly suggesting a clumpy wind. The comparison with the energetics of molecular outflows estimated in infrared and millimeter wavelengths support a connection between the nuclear and galaxy-scale winds in luminous AGN.

We have performed a high sensitivity observation of the UFO/BAL quasar APM 08279+5255 at z=3.912 with NOEMA at 3.2 mm, aimed at detecting fast moving molecular gas. We report the detection of blueshifted CO(4-3) with maximum velocity (v95\%) of $-1340$ km s$^{-1}$, with respect to the systemic peak emission, and a luminosity of $L' = 9.9\times 10^9 ~\mu^{-1}$ K km s$^{-1}$ pc$^{-2}$ (where $\mu$ is the lensing magnification factor). We discuss various scenarios for the nature of this emission, and conclude that this is the first detection of fast molecular gas at redshift $>3$. We derive a mass flow rate of molecular gas in the range $\rm \dot M=3-7.4\times 10^3$ M$_\odot$/yr, and momentum boost $\dot P_{OF} / \dot P_{AGN} \sim 2-6$, therefore consistent with a momentum conserving flow. For the largest $\dot P_{OF}$ the scaling is also consistent with a energy conserving flow with an efficiency of $\sim$10-20\%. The present data can hardly discriminate between the two expansion modes. The mass loading factor of the molecular outflow $\eta=\dot M_{OF}/SFR$ is $>>1$. We also detect a molecular emission line at a frequency of 94.83 GHz, corresponding to a rest frame frequency of 465.8 GHz, which we tentatively identified with the cation molecule $\rm N_2H^+$(5-4), which would be the first detection of this species at high redshift. We discuss the alternative possibility that this emission is due to a CO emission line from the, so far undetected, lens galaxy. Further observations of additional transitions of the same species with NOEMA can discriminate between the two scenarios.

We present new ALMA observations aimed at mapping molecular gas reservoirs through the CO(3-2) transition in three quasars at $z\simeq2.4$, LBQS 0109+0213, 2QZ J002830.4-281706, and [HB89] 0329-385. Previous [OIII]5007 observations of these quasars showed evidence for ionised outflows quenching star formation in their host galaxies. Systemic CO(3-2) emission has been detected only in one quasar, LBQS 0109+0213, where the CO(3-2) emission is spatially anti-correlated with the ionised outflow, suggesting that most of the molecular gas may have been dispersed or heated in the region swept by the outflow. In all three sources, including the one detected in CO, our constraints on the molecular gas mass indicate a significantly reduced reservoir compared to main-sequence galaxies at the same redshift, supporting a negative feedback scenario. In the quasar 2QZ J002830.4-281706, we tentatively detect an emission line blob blue-shifted by $v\sim-2000$ km/s with respect to the galaxy systemic velocity and spatially offset by 0.2 arcsec (1.7 kpc) with respect to the ALMA continuum peak. Interestingly, such emission feature is coincident in both velocity and space with the ionised outflow as seen in [OIII]5007. This tentative detection must be confirmed with deeper observations but, if real, it could represent the molecular counterpart of the ionised gas outflow driven by the AGN. Finally, in all ALMA maps we detect the presence of serendipitous line emitters within a projected distance $\sim 160$ kpc from the quasars. By identifying these features with the CO(3-2) transition, the serendipitous line emitters would be located within |$\Delta v$|$<$500 km/s from the quasars, hence suggesting an overdensity of galaxies in two out of three quasars.

Studying the coupling between the energy output produced by the central quasar and the host galaxy is fundamental to fully understand galaxy evolution. Quasar feedback is indeed supposed to dramatically affect the galaxy properties by depositing large amounts of energy and momentum into the ISM. In order to gain further insights on this process, we study the SEDs of sources at the brightest end of the quasar luminosity function, for which the feedback mechanism is supposed to be at its maximum. We model the rest-frame UV-to-FIR SEDs of 16 WISE-SDSS Selected Hyper-luminous (WISSH) quasars at 1.8 < z < 4.6 disentangling the different emission components and deriving physical parameters of both the nuclear component and the host galaxy. We also use a radiative transfer code to account for the contribution of the quasar-related emission to the FIR fluxes. Most SEDs are well described by a standard combination of accretion disk+torus and cold dust emission. However, about 30% of them require an additional emission component in the NIR, with temperatures peaking at 750K, which indicates the presence of a hotter dust component in these powerful quasars. We measure extreme values of both AGN bolometric luminosity (LBOL > 10^47 erg/s) and SFR (up to 2000 Msun/yr). A new relation between quasar and star-formation luminosity is derived (LSF propto LQSO^(0.73)) by combining several Herschel-detected quasar samples from z=0 to 4. Future observations will be crucial to measure the molecular gas content in these systems, probe the impact between quasar-driven outflows and on-going star-formation, and reveal the presence of merger signatures in their host galaxies.

We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field <U> with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

Similarly to the cosmic star formation history, the black hole accretion rate density of the Universe peaked at 1<z<3. This cosmic epoch is hence best suited for investigating the effects of radiative feedback from AGN. Observational efforts are underway to quantify the impact of AGN feedback, if any, on their host galaxies. Here we present a study of the molecular gas content of AGN hosts at z~1.5 using CO[2-1] line emission observed with ALMA for a sample of 10 AGNs. We compare this with a sample of galaxies without an AGN matched in redshift, stellar mass, and star formation rate. We detect CO in 3 AGNs with $\mathrm{L_{CO} \sim 6.3-25.1\times 10^{9} L_{\odot}}$ which translates to a molecular hydrogen gas mass of $\mathrm{2.5-10\times 10^{10} M_{\odot}}$ assuming conventional conversion factor of $\mathrm{\alpha_{CO}}\sim3.6$. Our results indicate a >99% probability of lower depletion time scales and lower molecular gas fractions in AGN hosts with respect to the non-AGN comparison sample. We discuss the implications of these observations on the impact that AGN feedback may have on star formation efficiency of z>1 galaxies.

We present the 78-ks Chandra observations of the $z=6.4$ quasar SDSS J1148+5251. The source is clearly detected in the energy range 0.3-7 keV with 42 counts (with a significance $\gtrsim9\sigma$). The X-ray spectrum is best-fitted by a power-law with photon index $\Gamma=1.9$ absorbed by a gas column density of $\rm N_{\rm H}=2.0^{+2.0}_{-1.5}\times10^{23}\,\rm cm^{-2}$. We measure an intrinsic luminosity at 2-10 keV and 10-40 keV equal to $\sim 1.5\times 10^{45}~\rm erg~s^{-1}$, comparable with luminous local and intermediate-redshift quasar properties. Moreover, the X-ray to optical power-law slope value ($\alpha_{\rm OX}=-1.76\pm 0.14$) of J1148 is consistent with the one found in quasars with similar rest-frame 2500 \AA ~luminosity ($L_{\rm 2500}\sim 10^{32}~\rm erg~s^{-1}$\AA$^{-1}$). Then we use Chandra data to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assuming that X-ray emission is attenuated by intervening, metal-rich ($Z\geq \rm Z_{\odot}$) molecular clouds distributed on $\sim$kpc scales in the host galaxy. Our analysis favors a black hole mass $M_{\rm BH} \sim 3\times 10^9 \rm M_\odot$ and a molecular hydrogen mass $M_{\rm H_2}\sim 2\times 10^{10} \rm M_\odot$, in good agreement with estimates obtained from previous studies. We finally discuss strengths and limits of our analysis.

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

We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.

Models and observations suggest that both power and effects of AGN feedback should be maximised in hyper-luminous (L_Bol>10^47 erg/s) quasars, i.e. objects at the brightest end of the AGN luminosity function. We present the first results of a multi-wavelength observing program, focusing on a sample of WISE/SDSS selected hyper-luminous (WISSH) broad-line quasars at z~1.5-5. The WISSH quasars project has been designed to reveal the most energetic AGN-driven outflows, estimate their occurrence at the peak of quasar activity and extend the study of correlations between outflows and nuclear properties up to poorly-investigated extreme AGN luminosities (L_Bol~10^47 -10^48 erg/s). We present NIR, long-slit LBT/LUCI1 spectroscopy of five WISSH quasars at z~2.3-3.5 showing prominent [OIII] emission lines with broad (FWHM~1200-2200 km/s) and skewed profiles. The luminosities of the broad [OIII] wings are the highest measured so far (L_[OIII]^broad >~5x10^44 erg/s) and reveal the presence of powerful ionised outflows with mass outflow rates Mdot >~1700 M_Sun/yr and kinetic powers Edot >~10^45 erg/s. Although these estimates are affected by large uncertainties, due to the use of [OIII] as tracer of ionized outflows and the very basic outflow model we assume, these results suggest that the AGN is highly efficient in pushing outwards large amounts of ionised gas in our targets. The mechanical outflow luminosities for WISSH quasars correspond to higher fractions (~1-3%) of L_Bol than those derived for lower L_Bol AGN. Our targets host very massive (M_BH>~2x10^9 M_Sun) black holes which are still accreting at a high rate (i.e. a factor of ~0.4-3 of the Eddington limit). These findings demonstrate that WISSH quasars offer the opportunity of probing the extreme end of both luminosity and SMBH mass functions and revealing powerful ionised outflows able to affect the evolution of their host galaxies.

We study the redshift evolution of the quasar UV Luminosity Function (LF) for 0.5 < z < 6.5, by collecting the most up to date observational data and, in particular, the recently discovered population of faint AGNs. We fit the QSO LF using either a double power-law or a Schechter function, finding that both forms provide good fits to the data. We derive empirical relations for the LF parameters as a function of redshift and, based on these results, predict the quasar UV LF at z=8. From the inferred LF evolution, we compute the redshift evolution of the QSO/AGN comoving ionizing emissivity and hydrogen photoionization rate. If faint AGNs are included, the contribution of quasars to reionization increases substantially. However, their level of contribution critically depends on the detailed shape of the QSO LF, which can be constrained by efficient searches of high-z quasars. To this aim, we predict the expected (i) number of z>6 quasars detectable by ongoing and future NIR surveys (as EUCLID and WFIRST), and (ii) number counts for a single radio-recombination line observation with SKA-MID (FoV = 0.49 deg^2) as a function of the Hnalpha flux density, at 0<z<8. These surveys (even at z<6) will be fundamental to better constrain the role of quasars as reionization sources.

We present ALMA observations of cold dust and molecular gas in four high-luminosity, heavily reddened (A$_{\rm{V}} \sim 2.5-6$ mag) Type 1 quasars at $z\sim2.5$ with virial M$_{\rm{BH}} \sim 10^{10}$M$_\odot$, to test whether dusty, massive quasars represent the evolutionary link between submillimetre bright galaxies (SMGs) and unobscured quasars. All four quasars are detected in both the dust continuum and in the $^{12}$CO(3-2) line. The mean dust mass is 6$\times$10$^{8}$M$_\odot$ assuming a typical high redshift quasar spectral energy distribution (T=41K, $\beta$=1.95 or T=47K, $\beta$=1.6). The implied star formation rates are very high - $\gtrsim$1000 M$_\odot$ yr$^{-1}$ in all cases. Gas masses estimated from the CO line luminosities cover $\sim$1-5$\times10^{10}$($\alpha_{\rm{CO}} / 0.8$)M$_\odot$ and the gas depletion timescales are very short - $\sim5-20$Myr. A range of gas-to-dust ratios is observed in the sample. We resolve the molecular gas in one quasar - ULASJ2315$+$0143 ($z=2.561$) - which shows a strong velocity gradient over $\sim$20 kpc. The velocity field is consistent with a rotationally supported gas disk but other scenarios, e.g. mergers, cannot be ruled out at the current resolution of these data. In another quasar - ULASJ1234+0907 ($z=2.503$) - we detected molecular line emission from two millimetre bright galaxies within 200 kpc of the quasar, suggesting that this quasar resides in a significant over-density. The high detection rate of both cold dust and molecular gas in these sources, suggests that reddened quasars could correspond to an early phase in massive galaxy formation associated with large gas reservoirs and significant star formation.

Negative feedback from active galactic nuclei (AGN) is considered a key mechanism in shaping galaxy evolution. Fast, extended outflows are frequently detected in the AGN host galaxies at all redshifts and luminosities, both in ionised and molecular gas. However, these outflows are only "potentially" able to quench star formation and we are still missing a decisive evidence of negative feedback in action. Here we present Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) H- and K-band integral-field spectroscopic observations of two quasars at $z\sim$2.4 characterised by fast, extended outflows detected through the [OIII]$\lambda$5007 line (Carniani et al. 2015). The high signal-to-noise ratio of our observations allows us to identify faint narrow (FWHM $< 500$ km/s), and spatially extended components in [OIII]$\lambda$5007 and H$\alpha$ emission associated with star formation in the host galaxy. Such star-formation powered emission is spatially anti-correlated with the fast outflows. The ionised outflows therefore appear to be able to suppress star formation in the region where the outflow is expanding. However the detection of narrow, spatially-extended H$\alpha$ emission indicates star formation rates of at least $\sim 50-90$ M$_{\odot}$/yr, suggesting either that AGN feedback does not affect the whole galaxy or that many feedback episodes are required before star formation is completely quenched. On the other hand, the narrow H$\alpha$ emission extending along the edges of the outflow cone may lead also to a positive feedback interpretation. Our results highlight the possible twofold role of galaxy-wide outflows in host galaxy evolution.

We explore the possibility of detecting hydrogen radio recombination lines from 0 < z < 10 quasars. We compute the expected Hnalpha flux densities as a function of absolute magnitude and redshift by considering (i) the range of observed AGN spectral indices from UV to X-ray bands, (ii) secondary ionizations from X-ray photons, and (iii) stimulated emission due to nonthermal radiation. All these effects are important to determine the line fluxes. We find that the combination of slopes: alpha_X,hard = -1.11, alpha_X,soft = -0.7, alpha_EUV = -1.3, alpha_UV = -1.7, maximizes the expected flux, f_Hnalpha = 10 microJy for z = 7 quasars with M_AB = -27 in the n = 50 lines; allowed SED variations produce variations by a factor of 3 around this value. Secondaries boost the line intensity by a factor of 2 to 4, while stimulated emission in high-z quasars with M_AB = -26 provides an extra boost to RRL flux observed at nu = 1 GHz if recombinations arise in HII regions with T_e = 10^3-5 K, n_e = 10^3-5 cm^-3. We compute the sensitivity required for a 5sigma detection of Hnalpha lines using the SKA, finding that the SKA-MID could detect sources with M_AB < -27 (M_AB < -26) at z < 8 (z < 3) in less than 100 hrs of observing time. These observations could open new paths to searches for obscured SMBH progenitors, complementing X-ray, optical/IR and sub-mm surveys.

Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star-forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure the molecular gas content of seven high-redshift (z ~ 1.6) starburst galaxies with the Atacama Large (sub-)Millimeter Array and IRAM/Plateau de Bure Interferometer. Our targets are selected from the sample of Herschel far-infrared detected galaxies having star formation rates (~300-800 Msolar/yr) elevated (>4x) above the star-forming `main sequence', and included in the FMOS-COSMOS near-infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. We detect CO emission in all cases at high levels of significance, indicative of high gas fractions (~30-50%). Even more compelling, we firmly establish with a clean and systematic selection that starbursts, identified as main-sequence outliers, at high redshift generally have a lower ratio of CO to total infrared luminosity as compared to typical main-sequence star-forming galaxies, although with a smaller offset than expected based on past studies of local starbursts. We put forward a hypothesis that there exists a continuous increase in star formation efficiency with elevation from the main sequence with galaxy mergers as a possible physical driver. Along with a heightened star formation efficiency, our high-redshift sample is similar in other respects to local starbursts such as being metal rich and having a higher ionization state of the interstellar medium.

We present the best sensitivity and angular resolution maps of the molecular disk and outflow of Mrk 231, as traced by CO observations obtained with IRAM/PdBI, and we analyze archival Chandra and NuSTAR observations. We constrain the physical properties of both the molecular disk and outflow, the presence of a highly-ionized ultra-fast nuclear wind, and their connection. The molecular outflow has a size of ~1 kpc, and extends in all directions around the nucleus, being more prominent along the south-west to north-east direction, suggesting a wide-angle biconical geometry. The maximum projected velocity of the outflow is nearly constant out to ~1 kpc, thus implying that the density of the outflowing material decreases from the nucleus outwards as $r^{-2}$. This suggests that either a large part of the gas leaves the flow during its expansion or that the bulk of the outflow has not yet reached out to ~1 kpc, thus implying a limit on its age of ~1 Myr. We find $\dot M_{OF}=[ 500-1000]~ M_{\odot}~yr^{-1}$ and $\dot E_{kin,OF}=[7-10]\times 10^{43}$ erg s$^{-1}$. Remarkably, our analysis of the X-ray data reveals a nuclear ultra-fast outflow (UFO) with velocity -20000 km s$^{-1}$, $\dot M_{UFO}=[0.3- 2.1] ~M_\odot yr^{-1}$, and momentum load $\dot P_{UFO}/\dot P_{rad}=[0.2-1.6]$.We find $\dot E_{kin,UFO}\sim \dot E_{kin,OF}$ as predicted for outflows undergoing an energy conserving expansion. This suggests that most of the UFO kinetic energy is transferred to mechanical energy of the kpc-scale outflow, strongly supporting that the energy released during accretion of matter onto super-massive black holes is the ultimate driver of giant massive outflows. We estimate a momentum boost $\dot P_{OF}/\dot P_{UFO}\approx [30-60]$. The ratios $\dot E_{kin, UFO}/L_{bol,AGN} =[ 1-5]\%$ and $\dot E_{kin,OF}/L_{bol,AGN} = [1-3]\%$ agree with the requirements of the most popular models of AGN feedback.

We study the ionization and kinematics of the ionized gas in the nuclear region of the barred Seyfert 2 galaxy NGC~5643 using MUSE integral field observations in the framework of the MAGNUM (Measuring Active Galactic Nuclei Under MUSE Microscope) survey. The data were used to identify regions with different ionization conditions and to map the gas density and the dust extinction. We find evidence for a double sided ionization cone, possibly collimated by a dusty structure surrounding the nucleus. At the center of the ionization cone, outflowing ionized gas is revealed as a blueshifted, asymmetric wing of the [OIII] emission line, up to projected velocity v(10)~-450 km/s. The outflow is also seen as a diffuse, low luminosity radio and X-ray jet, with similar extension. The outflowing material points in the direction of two clumps characterized by prominent line emission with spectra typical of HII regions, located at the edge of the dust lane of the bar. We propose that the star formation in the clumps is due to `positive feedback' induced by gas compression by the nuclear outflow, providing the first candidate for outflow induced star formation in a Seyfert-like radio quiet AGN. This suggests that positive feedback may be a relevant mechanism in shaping the black hole-host galaxy coevolution.

Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $\mu$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($\Delta I/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($\nu_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/\mu{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4\sigma$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,\mu$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.

AGN-driven outflows are invoked by galaxy evolutionary models to quench star formation and to explain the origin of the relations observed locally between super massive black holes and their host galaxies. This work aims to detect the presence of extended ionised outflows in luminous quasars where we expect the maximum activity both in star formation and in black hole accretion. Currently, there are only a few studies based on spatially resolved observations of outflows at high redshift, $z>2$. We analyse a sample of six luminous (${\rm L>10^{47} \ erg/s}$) quasars at $z\sim2.4$, observed in H-band using the near-IR integral field spectrometer SINFONI at VLT. We perform a kinematic analysis of the [OIII] emission line at $\lambda = 5007\AA$. [OIII] has a complex gas kinematic, with blue-shifted velocities of a few hundreds of km/s and line widths up to 1500 km/s. Using the spectroastrometric method we infer size of the ionised outflows of up to $\sim$2 kpc. The properties of the ionised outflows, mass outflow rate, momentum rate and kinetic power, are correlated with the AGN luminosity. The increase in outflow rate with increasing AGN luminosity is consistent with the idea that a luminous AGN pushes away the surrounding gas through fast outflows driven by radiation pressure, which depends on the emitted luminosity. We derive mass outflow rates of about 6-700 M$_{\odot}$/yr for our sample, which are lower than those observed in molecular outflows. Indeed physical properties of ionised outflows show dependences on AGN luminosity which are similar to those of molecular outflows but indicating that the mass of ionised gas is smaller than that of the molecular one. Alternatively, this discrepancy between ionised and molecular outflows could be explained with different acceleration mechanisms.

Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching.

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

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4<S/N<=4.5 if a multiwavelength counterpart within 1.5" is also present, the total number of detections in our survey is 22. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or multiple components. Using optical to NIR photometric redshifts, available spectroscopic redshifts, and redshifts estimated from the radio-to-submm spectral index we infer a median redshift of $\tilde{z}=3.20\pm0.25$ for our sample. To study the overall multiplicity and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest AzTEC SMGs detected at 1.1 mm, AzTEC1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. We find that the median redshift for the 15 brightest AzTEC SMGs ($\tilde{z}=3.05\pm0.44$) is consistent with that for AzTEC16-30. This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux density. For the combined AzTEC1-30 sample we derive a median redshift of $\tilde{z}=3.17\pm0.27$, consistent with previous results based on mm-selected samples. We further infer that within the combined AzTEC1-30 sample $\sim25\pm9\%$ of sources separate into multiple components.

Even 10 billion years ago, the cores of the first galaxy clusters are often found to host a characteristic population of massive galaxies with already suppressed star formation. Here we search for distant cluster candidates at z~2 using massive passive galaxies as tracers. With a sample of ~40 spectroscopically confirmed passive galaxies at 1.3<z<2.1, we tune photometric redshifts of several thousands passive sources in the full 2 sq.deg. COSMOS field. This allows us to map their density in redshift slices, probing the large scale structure in the COSMOS field as traced by passive sources. We report here on the three strongest passive galaxy overdensities that we identify in the redshift range 1.5<z<2.5. While the actual nature of these concentrations is still to be confirmed, we discuss their identification procedure, and the arguments supporting them as candidate galaxy clusters (likely mid-10^13 M_sun range). Although this search approach is likely biased towards more evolved structures, it has the potential to select still rare, cluster-like environments close to their epoch of first appearance, enabling new investigations of the evolution of galaxies in the context of structure growth.

We present the first X-ray spectrum of a Hot dust-obscured galaxy (DOG), namely W1835+4355 at z ~ 2.3. Hot DOGs represent a very rare population of hyperluminous (>= 10^47 erg/s), dust-enshrouded objects at z > 2 recently discovered in the WISE All Sky Survey. The 40 ks XMM-Newton spectrum reveals a continuum as flat (Gamma ~ 0.8) as typically seen in heavily obscured AGN. This, along with the presence of strong Fe Kalpha emission, clearly suggests a reflection-dominated spectrum due to Compton-thick absorption. In this scenario, the observed luminosity of L(2-10 keV) ~ 2 x 10^44 erg/s is a fraction (<10%) of the intrinsic one, which is estimated to be >~ 5 x 10^45 erg/s by using several proxies. The Herschel data allow us to constrain the SED up to the sub-mm band, providing a reliable estimate of the quasar contribution (~ 75%) to the IR luminosity as well as the amount of star formation (~ 2100 Msun/yr). Our results thus provide additional pieces of evidence that associate Hot DOGs with an exceptionally dusty phase during which luminous quasars and massive galaxies co-evolve and a very efficient and powerful AGN-driven feedback mechanism is predicted by models.

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

Past observations of QSO host galaxies at z >6 have found cold gas and star formation on compact scales of a few kiloparsecs. We present new high sensitivity IRAM PdBI follow-up observations of the [CII] 158micron emission line and FIR continuum in the host galaxy of SDSS J1148+5152, a luminous QSO at redshift 6.4189. We find that a large fraction of the gas traced by [CII] is at high velocities, up to ~1400 km/s relative to the systemic velocity, confirming the presence of a major quasar-driven outflow indicated by previous observations. The outflow has a complex morphology and reaches a maximum projected radius of ~30 kpc. The extreme spatial extent of the outflow allows us, for the first time in an external galaxy, to estimate mass-loss rate, kinetic power and momentum rate of the outflow as a function of the projected distance from the nucleus and the dynamical time-scale. These trends reveal multiple outflow events during the past 100 Myr, although the bulk of the mass, energy and momentum appear to have been released more recently, within the past ~20 Myr. Surprisingly, we discover that also the quiescent gas at systemic velocity is extremely extended. More specifically, we find that, while 30% of the [CII] within v\in(-200, 200) km/s traces a compact component that is not resolved by our observations, 70% of the [CII] emission in this velocity range is extended, with a projected FWHM size of 17.4+-1.4 kpc. We detect FIR continuum emission associated with both the compact and the extended [CII] components, although the extended FIR emission has a FWHM of 11+-3 kpc, thus smaller than the extended [CII] source. Overall, our results indicate that the cold gas traced by [CII] is distributed up to r~30 kpc. A large fraction of extended [CII] is likely associated with star formation on large scales, but the [CII] source extends well beyond the FIR continuum.

Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, preventing massive galaxies to over-grow and producing the red colors of ellipticals. On the other hand, some models are also requiring `positive' AGN feedback, inducing star formation in the host galaxy through enhanced gas pressure in the interstellar medium. However, finding observational evidence of the effects of both types of feedback is still one of the main challenges of extragalactic astronomy, as few observations of energetic and extended radiatively-driven winds are available. Here we present SINFONI near infrared integral field spectroscopy of XID2028, an obscured, radio-quiet z=1.59 QSO detected in the XMM-COSMOS survey, in which we clearly resolve a fast (1500 km/s) and extended (up to 13 kpc from the black hole) outflow in the [OIII] lines emitting gas, whose large velocity and outflow rate are not sustainable by star formation only. The narrow component of Ha emission and the rest frame U band flux from HST-ACS imaging enable to map the current star formation in the host galaxy: both tracers independently show that the outflow position lies in the center of an empty cavity surrounded by star forming regions on its edge. The outflow is therefore removing the gas from the host galaxy (`negative feedback'), but also triggering star formation by outflow induced pressure at the edges (`positive feedback'). XID2028 represents the first example of a host galaxy showing both types of feedback simultaneously at work.

Understanding the relationship between the formation and evolution of galaxies and their central super massive black holes (SMBH) is one of the main topics in extragalactic astrophysics. Links and feedback may reciprocally affect both black hole and galaxy growth. Observations of the CO line at redshifts of 2-4 are crucial to investigate the gas mass, star formation activity and accretion onto SMBHs, as well as the effect of AGN feedback. Potential correlations between AGN and host galaxy properties can be highlighted by observing extreme objects. Despite their luminosity, hyper-luminous QSOs at z=2-4 are still little studied at mm wavelengths. We targeted CO(3-2) in ULAS J1539+0557, an hyper-luminos QSO (Lbol> 10^48 erg/s) at z=2.658, selected through its unusual red colors in the UKIDSS Large Area Survey (ULAS). We find a molecular gas mass of 4.1+-0.8 10^10 Msun, and a gas fraction of 0.4-0.1, depending mostly on the assumed source inclination. We also find a robust lower limit to the star-formation rate (SFR=250-1600 Msun/yr) and star-formation efficiency (SFE=25-350 Lsun/(K km s-1 pc2) by comparing the observed optical-near-infrared spectral energy distribution with AGN and galaxy templates. The black hole gas consumption timescale, M(H_2)/dM(accretion)/dt, is ~160 Myr, similar or higher than the gas consumption timescale. The gas content and the star formation efficiency are similar to those of other high-luminosity, highly obscured QSOs, and at the lower end of the star-formation efficiency of unobscured QSOs, in line with predictions from AGN-galaxy co-evolutionary scenarios. Further measurements of the (sub)-mm continuum in this and similar sources are mandatory to obtain a robust observational picture of the AGN evolutionary sequence.

We study the properties of massive, galactic-scale outflows of molecular gas and investigate their impact on galaxy evolution. We present new IRAM PdBI CO(1-0) observations of local ULIRGs and QSO hosts: clear signature of massive and energetic molecular outflows, extending on kpc scales, is found in the CO(1-0) kinematics of four out of seven sources, with measured outflow rates of several 100 Msun/yr. We combine these new observations with data from the literature, and explore the nature and origin of massive molecular outflows within an extended sample of 19 local galaxies. We find that starburst-dominated galaxies have an outflow rate comparable to their SFR, or even higher by a factor of ~ 2-4, implying that starbursts can indeed be effective in removing cold gas from galaxies. Nevertheless, our results suggest that the presence of an AGN can boost the outflow rate by a large factor, which is found to increase with the L_AGN/L_bol ratio. The gas depletion time-scales due to molecular outflows are anti-correlated with the presence and luminosity of an AGN in these galaxies, and range from a few hundred million years in starburst galaxies, down to just a few million years in galaxies hosting powerful AGNs. In quasar hosts the depletion time-scales due to the outflow are much shorter than the depletion time-scales due to star formation. We estimate the outflow kinetic power and find that, for galaxies hosting powerful AGNs, it corresponds to about 5% of the AGN luminosity, as expected by models of AGN feedback. Moreover, we find that momentum rates of about 20 L_AGN/c are common among the AGN-dominated sources in our sample. For "pure" starburst galaxies our data tentatively support models in which outflows are mostly momentum-driven by the radiation pressure from young stars onto dusty clouds.

We present high spatial resolution CO(1-0) mapping of the luminous infrared galaxy NGC6240 obtained with the IRAM - Plateau de Bure Interferometer. This source is a well-known early-stage merging system hosting two AGN. We find a broad CO(1-0) line profile, with maximum velocity 800 km/s and FWZI=1400 km/s, and displaying several kinematic components, revealing the complexity of the gas dynamics in this system. We detect a blueshifted CO emission with velocity in the range -200 and -500 km/s, peaked around the southern AGN, at the same position where the H2 outflow is located, and with a mass loss rate of ~500 Msun/yr. We interpret this blueshifted component as a outflow, originating from the southern nucleus. The spatial and spectral match strongly suggests that the CO outflow is connected to the H2 superwind located around the southern nucleus, and to the large scale CO outflow, with similar velocities, extended on scales of 10 kpc. The large mass loading factor (dM/dt /SFR ~ 10) of the molecular gas suggests that the outflow is likely driven by both SNa winds and the radiation of the southern AGN. We discovered a nuclear, redshifted CO emission peaking in the mid point of the two nuclei, as it is the case for the CO emission at the systemic velocity. The large velocity dispersion, which reaches its maximum (~500 Km/s) in the mid point between the two nuclei, suggests that the gas might be highly turbulent in this region, although the presence of an unresolved rotation component cannot be ruled out.

Understanding the astrophysics of feedback in active galactic nuclei (AGN) is key to understanding the growth and co-evolution of supermassive black holes and galaxies. AGN-driven winds/outflows are potentially the most effective way of transporting energy and momentum from the nuclear scales to the host galaxy, quenching star formation by sweeping away the gas reservoir. Key questions in this field are: 1) how do accretion disks around black holes launch winds/outflows, and how much energy do these carry? 2) How are the energy and metals accelerated in winds/outflows transferred and deposited into the circumgalactic medium? X-ray observations are a unique way to address these questions because they probe the phase of the outflows which carries most of the kinetic energy. We show how a high throughput, high spectral resolution instrument like the X-ray Integral Field Unit (X-IFU) on Athena+ will allow us to address these questions by determining the physical parameters (ionization state, density, temperature, abundances, velocities, geometry, etc.) of the outflows on a dynamical time-scale, in a broad sample of nearby bright AGN. The X-IFU will also allow direct spectral imaging of the impact of these winds on the host galaxy for local AGN, forming a template for understanding AGN at higher redshifts where wind shocks cannot be resolved.

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

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

We present deep CO observations of NGC6240 performed with the IRAM Plateau de Bure Interferometer (PdBI). NGC6240 is the prototypical example of a major galaxy merger in progress, caught at an early stage, with an extended, strongly-disturbed butterfly-like morphology and the presence of a heavily obscured active nucleus in the core of each progenitor galaxy. The CO line shows a skewed profile with very broad and asymmetric wings detected out to velocities of -600 km/s and +800 km/s with respect to the systemic velocity. The PdBI maps reveal the existence of two prominent structures of blueshifted CO emission. One extends eastward, i.e. approximately perpendicular to the line connecting the galactic nuclei, over scales of ~7 kpc and shows velocities up to -400 km/s. The other extends southwestward out to ~7 kpc from the nuclear region, and has a velocity of -100 km/s with respect to the systemic one. Interestingly, redshifted emission with velocities 400 to 800 km/s is detected around the two nuclei, extending in the east-west direction, and partly overlapping with the eastern blue-shifted structure, although tracing a more compact region of size ~1.7 kpc. The overlap between the southwestern CO blob and the dust lanes seen in HST images, which are interpreted as tidal tails, indicates that the molecular gas is deeply affected by galaxy interactions. The eastern blueshifted CO emission is co-spatial with an Halpha filament that is associated with strong H2 and soft X-ray emission. The analysis of Chandra X-ray data provides strong evidence for shocked gas at the position of the Halpha emission. Its association with outflowing molecular gas supports a scenario where the molecular gas is compressed into a shock wave that propagates eastward from the nuclei. If this is an outflow, the AGN are likely the driving force.

We report the results from a 2011 Suzaku observation of the nearby low-ionization BAL quasar/ULIRG Markarian 231. These data reveal that the X-ray spectrum has undergone a large variation from the 2001 XMM-Newton and BeppoSAX observations. We interpret this finding according to a scenario whereby the X-ray continuum source is obscured by a two-component partial-covering absorber with NH ~10^22 and ~10^24 cm^-2, respectively. The observed spectral change is mostly explained by a progressive appearance of the primary continuum at <10 keV due to the decrease of the covering fraction of the denser absorption component. The properties of the X-ray obscuration in Mrk 231 match well with those of the X-ray shielding gas predicted by the theoretical models for an efficient radiatively-driven acceleration of the BAL wind. In particular, the X-ray absorber might be located at the extreme base of the outflow. We measure a 2-10 keV luminosity of L(2-10) = 3.3 x 10^43 erg s^-1 for the 2011 data set, i.e. an increase of 30% with respect to the 2001 value.

We present observations of the CO[3-2] emission towards two massive and infrared luminous Lyman Break Galaxies at z = 3.21 and z = 2.92, using the IRAM Plateau de Bure Interferometer, placing first constraints on the molecular gas masses (Mgas) of non-lensed LBGs. Their overall properties are consistent with those of typical (Main-Sequence) galaxies at their redshifts, with specific star formation rates ~1.6 and ~2.2 Gyr^(-1), despite their large infrared luminosities L_IR ~2-3 x 10^12 Lsun derived from Herschel. With one plausible CO detection (spurious detection probability of 10^(-3)) and one upper limit, we investigate the evolution of the molecular gas-to-stellar mass ratio (Mgas/M*) with redshift. Our data suggest that the steep evolution of Mgas/M* of normal galaxies up to z~2 is followed by a flattening at higher redshifts, providing supporting evidence for the existence of a plateau in the evolution of the specific star formation rate at z > 2.5.

Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.

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

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

We present a long (~76 ks) Chandra observation of IRAS 09104+4109, a hyper-luminous galaxy, optically classified as a Type 2 AGN hosted in a cD galaxy in a cluster at z=0.442. We also report on the results obtained by fitting its broad-band spectral energy distribution. The Compton-thick nature of this source (which has been often referred to as an "archetype" of Compton-thick Type 2 quasars) was formerly claimed on the basis of its marginal detection in the PDS instrument onboard BeppoSAX, being then disputed using XMM-Newton data. Both Chandra analysis and optical/mid-IR spectral fitting are consistent with the presence of heavy (~1-5 10^{23} cm^{-2}) but not extreme (Compton-thick) obscuration. However, using the mid-IR and the [OIII] emission as proxies of the nuclear hard X-ray luminosity suggests the presence of heavier obscuration. The 54-month Swift BAT map shows excess hard X-ray emission likely related to a nearby (z=0.009) Type 2 AGN, close enough to IRAS 09104+4109 to significantly enhance and contaminate its emission in the early BeppoSAX PDS data.

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

Tracing the formation and evolution of all supermassive black holes, including the obscured ones, understanding how black holes influence their surroundings and how matter behaves under extreme conditions, are recognized as key science objectives to be addressed by the next generation of instruments. These are the main goals of the COSPIX proposal, made to ESA in December 2010 in the context of its call for selection of the M3 mission. In addition, COSPIX, will also provide key measurements on the non thermal Universe, particularly in relation to the question of the acceleration of particles, as well as on many other fundamental questions as for example the energetic particle content of clusters of galaxies. COSPIX is proposed as an observatory operating from 0.3 to more than 100 keV. The payload features a single long focal length focusing telescope offering an effective area close to ten times larger than any scheduled focusing mission at 30 keV, an angular resolution better than 20 arcseconds in hard X-rays, and polarimetric capabilities within the same focal plane instrumentation. In this paper, we describe the science objectives of the mission, its baseline design, and its performances, as proposed to ESA.

We report the detection of the 6.4 keV Iron K{\alpha} emission line in two infrared-luminous, massive, star-forming BzK galaxies at z = 2.578 and z = 2.90 in the CDF-S. The Chandra 4 Ms spectra of BzK4892 and BzK8608 show a reflection dominated continuum with strong Iron lines, with rest-frame equivalent widths EW-2.3 keV and 1.2 keV, respectively, demonstrating Compton thick obscuration of the central AGN. For BzK8608 the line identification closely matches the existing photometric redshift derived from the stellar emission. We use the observed luminosities of the Iron K{\alpha} line, of the rest-frame mid-IR continuum and of the UV rest-frame narrow emission lines to infer intrinsic L(2-10 keV)>1e44 erg/s, about 1.0-2.5 dex larger than the observed ones, hence confirming the presence of an absorber with N(H) > 1e24 cm-2. The two BzK galaxies have stellar masses of 5e10 M\odot and, based on VLA 1.4 GHz and submm 870{\mu}m observations, they appear to host vigorous starburst activity with SFR - 300-700 M\odot/yr that is also optically thick. We estimate that the AGN might also conceivably account for an important fraction of the bolometric far-IR emission of the galaxies. The implied volume density of Compton thick (CT) AGN with L(2-10 keV)>1e44 erg/s is in agreement with predictions from X-ray background synthesis models. These sources provide one of the first clearcut observations of the long-sought phase of simultaneous, heavily obscured quasar and star formation activity, predicted by models of massive galaxy evolution at high redshifts.

In the standard scenario for galaxy evolution young star-forming galaxies transform into red bulge-dominated spheroids, where star formation has been quenched. To explain such a transformation, a strong negative feedback generated by accretion onto a super-massive black hole is often invoked. The depletion of gas resulting from quasar-driven outflows should eventually stop star-formation in the host galaxy and lead the black hole to "suicide" for starvation. Direct observational evidence for a major quasar feedback onto the host galaxy is still missing, since outflows previously observed in quasars are generally associated with the ionized component of the gas, which only accounts for a minor fraction of the total gas content, and typically occurring in the central regions. We used the IRAM PdB Interferometer to observe CO(1-0) in Mrk 231, the closest QSO known. Thanks to the wide band we detect broad wings of the CO line, with velocities up to 750 km/s and spatially resolved on the kpc scale. Such broad CO wings trace a giant molecular outflow of about 700 MSun/yr, far larger than the ongoing star-formation rate (~200 MSun/yr) observed in the host galaxy. This wind will totally expel the cold gas reservoir in Mrk 231 in about 1e7 yrs, therefore halting the star-formation activity on the same timescale. The inferred kinetic energy in the molecular outflow is ~1.2e44 erg/s, corresponding to a few percent of the AGN bolometric luminosity, very close to the fraction expected by models ascribing quasar feedback to highly supersonic shocks generated by radiatively accelerated nuclear winds. The direct observational evidence for quasar feedback reported here provides solid support to the scenarios ascribing the observed properties of local massive galaxies to quasar-induced large scale winds.

We investigate the effects of the environment on star-formation in a sample of massive luminous and ultra-luminous infrared galaxies (LIRGs and ULIRGs) with S(24 micron)>80 uJy and i+<24 in the COSMOS field. We exploit the accurate photometric redshifts in COSMOS to characterize the galaxy environment and study the evolution of the fraction of LIRGs and ULIRGs in different environments in the redshift range z=0.3-1.2 and in bins of stellar mass. We find that the environment plays a role in the star formation processes and evolution of LIRGs and ULIRGs. We find an overall increase of the ULIRG+LIRG fraction in an optically-selected sample with increasing redshift, as expected from the evolution of the star formation rate density. We find that the ULIRG+LIRG fraction decreases with increasing density up to z~1, and that the dependence on density flattens with increasing redshift. We do not observe the reversal of the star-formation rate density relation up to z=1 in massive LIRGs and ULIRGs, suggesting that such reversal might occur at higher redshift in this infrared luminosity range.

We present the first results obtained from the identification of ~30,000 sources in the 24mic observations of the COSMOS field at S24>~80muJy. Using accurate photo-z's and extrapolations of the counts at faint fluxes we resolve with unprecedented detail the build-up of the mid-IR background across cosmic ages. We find that ~50% and ~80% of the 24mic background intensity originate from galaxies at z<~1 and z<~2 respectively, supporting the scenario where highly obscured sources at very high redshifts (z>~2) contribute only marginally to the CIB. Assuming flux-limited selections at optical wavelengths, we also find that the fraction of i-band sources with 24mic detection strongly increases up to z~2 as a consequence of the rapid evolution that star-forming galaxies have undergone with lookback time. Nonetheless this rising trend shows a clear break at z~1.3, probably due to k-correction effects implied by the complexity of SEDs in the mid-IR. Finally, we compare our results with the predictions from different models of galaxy formation. We note that semi-analytical formalisms fail to reproduce the redshift distributions observed at 24mic. Furthermore the simulated galaxies exhibit R-K colors much bluer than observed and the predicted K-band fluxes are systematically underestimated at z>0.5. Unless these discrepancies result from an incorrect treatment of extinction they may reflect an underestimate of the predicted density of high-z massive sources with strong star formation, which would point to more fundamental processes that are still not fully controlled in the simulations. The most recent backward evolution scenarios reproduce reasonably well the flux/redshift distribution of 24mic sources up to z~3, although none of them is able to exactly match our results at all redshifts. [Abridged]

We present a catalog of optical spectroscopic identifications of sources detected by Spitzer at 3.6 or 24 micron down to 10 and 280 microJy, respectively, in the SWIRE/XMM-Newton/ELAIS-S1 field and classified via line width analysis and diagnostic diagrams. A total of 1376 sources down to R~24.2 mag have been identified (1362 detected at 3.6 micron, 419 at 24 micron, and 405 at both) by low-resolution optical spectroscopy carried out with FORS2, VIMOS, and EFOSC2 at the Very Large Telescope and 3.6m ESO telescopes. The spectroscopic campaigns have been carried out over the central 0.6 square degrees area of ELAIS-S1 which, in particular, has also been observed by XMM-Newton and Chandra. We find the first direct optical spectroscopic evidence that the fraction of active galactic nuclei (AGN; mostly AGN2) increases with increasing F(24 micron)/F(R) ratio, reaching values of 70(+/-20)% in the range 316<F(24 micron)F(R)<1000. We present an IRAC-MIPS color-color diagram able to separate AGN1 from obscured AGN2 candidates. After having corrected for the spectroscopic incompleteness of our sample, it results that the AGN fraction at F(24 micron)=0.8 mJy is ~22(+/-7)% and decreases slowly to ~19(+/-5)% down to F(24 micron)=0.3 mJy.

Recent works have suggested that selection criteria based on MIR colors can be used to reveal a population of dust-enshrouded, extremely luminous quasars at z>1. However the X-ray spectral properties of these intriguing sources still remain largely unexplored. We report on an X-ray spectroscopic study of a sample of 44 very bright mid-IR galaxies with extreme mid-IR to optical flux ratios (MIR/O>2000). The X-ray coverage of the sample is highly inhomogeneous (from snap-shot 5 ks Chandra observations to medium-deep XMM exposures of 70 ks) and, consequently, a sizable fraction of them (~43%) remains undetected in the 0.5-10 keV band. The vast majority (95%) of the detected sources (23) show an absorption column density NH>10e22 cm-2 and, remarkably, we also find that 50% of them can be classified as Type 2 quasars on the basis of their absorption properties and X-ray luminosity. Moreover, most of the X-ray undetected sources show extreme mid-IR colors, consistent with being luminous AGN-powered objects, suggesting they might host heavily obscured (possibly Compton-thick) quasars in X-rays. This demonstrates that our selection criteria applied to a wide area survey is very efficient in finding a large number of Type 2 quasars at z > 1. The existence of this class of very powerful, obscured quasars at high z could have important implications in the context of the formation and cosmological evolution of accreting supermassive black holes and their host galaxies.

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

We present optical identifications and a multi-band catalogue of a sample of 478 X-ray sources in the XMM and Chandra surveys of the central 0.6 deg^2 of the ELAIS-S1 field. The optical/infrared counterpart of each X-ray source was identified using R and IRAC 3.6 um bands. This method was complemented by the precise positions obtained through Chandra observations. Approximately 94% of the counterparts are detected in the R band, while the remaining are blank fields in the optical down to R~24.5, but have a near-infrared counterpart detected by IRAC within 6 arcsec from the XMM centroid. The multi-band catalogue contains photometry in ten photometric bands (B to the MIPS 24 um). We determined redshift and classification for 237 sources (~50% of the sample) brighter than R=24. We classified 47% of the sources with spectroscopic redshift as broad-line active galactic nuclei (BL AGNs) with z=0.1-3.5, while sources without broad-lines are about 46% of the spectroscopic sample and are found up to z=2.6. We identified 11 type 2 QSOs among the sources with X/O>8, with z=0.9-2.6, high 2-10 keV luminosity (log(L2-10keV)>=43.8 erg/s) and hard X-ray colors suggesting large absorbing columns at the rest frame (logN_H up to 23.6 cm^-2). BL AGNs show on average blue optical-to-near-infrared colors, softer X-ray colors and X-ray-to-optical colors typical of optically selected AGNs. Conversely, narrow-line sources show redder optical colors, harder X-ray flux ratio and span a wider range of X-ray-to-optical colors. On average the SEDs of high-luminosity BL AGNs resemble the power-law typical of unobscured AGNs. The SEDs of NOT BL AGNs are dominated by the galaxy emission in the optical/near-infrared, and show a rise in the mid-infrared which suggests the presence of an obscured active nucleus.

(abridged) The ESIS survey is the optical follow up of the SWIRE/Spitzer in the ELAIS-S1 region of the sky. In the era of observational cosmology, the main efforts are focused on the study of galaxy evolution and its environmental dependence. Wide area, multiwavelength, extragalactic surveys are needed in order to probe sufficiently large volumes, minimize cosmic variance and find significant numbers of rare objects. We present VIMOS I and z band imaging belonging to the ESIS survey. A total of ~4 deg2 were targeted in I and ~1 deg2 in z. More than 300000 galaxies have been detected in the I band and ~50000 in the z band. Object coordinates are defined within an uncertainty of ~0.2 arcsec r.m.s., with respect to GSC 2.2. We reach a 90% average completeness at 23.1 and 22.5 mag (Vega) in the I and z bands, respectively. On the basis of IRAC colors, we identified galaxies having the 1.6 um stellar peak shifted to z=1-3. The new I, z band data provide reliable constraints to avoid low-redshift interlopers and reinforce this selection. Roughly 1000 galaxies between z=2-3 were identified over the ESIS ~4 deg2, at the SWIRE 5.8 um depth (25.8 uJy at 3sigma). These are the best galaxy candidates to dominate the massive tail (M>1e11 Msun) of the z>2 mass function.

Simbol-X is a French-Italian mission, with a participation of German laboratories, for X-ray astronomy in the wide 0.5-80 keV band. Taking advantage of emerging technology in mirror manufacturing and spacecraft formation flying, Simbol-X will push grazing incidence imaging up to ~80 keV, providing an improvement of roughly three orders of magnitude in sensitivity and angular resolution compared to all instruments that have operated so far above 10 keV. This will open a new window in X-ray astronomy, allowing breakthrough studies on black hole physics and census and particle acceleration mechanisms. We describe briefly the main scientific goals of the Simbol-X mission, giving a few examples aimed at highlighting key issues of the Simbol-X design.

(Abridged) We selected high-z massive galaxies at 5.8 microns, in the SWIRE ELAIS-S1 field (1 sq. deg.). Galaxies with the 1.6 microns stellar peak redshifted into the IRAC bands (z~1-3, called ``IR-peakers'') were identified. Stellar masses were derived by means of spectro-photometric fitting and used to compute the stellar mass function (MF) at z=1-2 and 2-3. A parametric fit to the MF was performed, based on a Bayesian formalism, and the stellar mass density of massive galaxies above z=2 determined. We present the first systematic study of the very-massive tail of the galaxy stellar mass function at high redshift. A total of 326 sources were selected. The majority of these galaxies have stellar masses in excess of 1e11 Msun and lie at z>1.5. The availability of mid-IR data turned out to be a valuable tool to constrain the contribution of young stars to galaxy SEDs, and thus their M(stars)/L ratio. The influence of near-IR data and of the chosen stellar library on the SED fitting are also discussed. A significant evolution is found not only for galaxies with M~1e11 Msun, but also in the highest mass bins considered. The comoving number density of these galaxies was lower by more than a factor of 10 at z=2-3, with respect to the local estimate. SWIRE 5.8 micron peakers more massive than 1.6x1e11 Msun provide 30-50% of the total stellar mass density in galaxies at z=2-3.

A large population of heavily obscured, Compton Thick AGNs is predicted by models of galaxy formation, models of Cosmic X-ray Background and by the ``relic'' super-massive black-hole mass function measured from local bulges. However, so far only a handful of Compton thick AGNs have been possibly detected using even the deepest Chandra and XMM surveys. Compton-thick AGNs can be recovered thanks to the reprocessing of the AGN UV emission in the infrared by selecting sources with AGN luminosity's in the mid-infrared and faint near-infrared and optical emission. To this purpose, we make use of deep HST, VLT, Spitzer and Chandra data on the Chandra Deep Field South to constrain the number of Compton thick AGN in this field. We show that sources with high 24$\mu$m to optical flux ratios and red colors form a distinct source population, and that their infrared luminosity is dominated by AGN emission. Analysis of the X-ray properties of these extreme sources shows that most of them (80$\pm15%$) are indeed likely to be highly obscured, Compton thick AGNs. The number of infrared selected, Compton thick AGNs with 5.8$\mu$m luminosity higher than $10^{44.2}$ erg s$^{-1}$ turns out to be similar to that of X-ray selected, unobscured and moderately obscured AGNs with 2-10 keV luminosity higher than $10^{43}$ erg s$^{-1}$ in the redshift bin 1.2-2.6. This ``factor of 2'' source population is exactly what it is needed to solve the discrepancies between model predictions and X-ray AGN selection.

Context: The counts of galaxies and AGN in the mid infra-red (MIR) bands are important instruments for studying their cosmological evolution. However, the classic spectral line ratios techniques can become misleading when trying to properly separate AGN from starbursts or even from apparently normal galaxies. Aims: We use X-ray band observations to discriminate AGN activity in previously classified MIR-selected starburst galaxies and to derive updated AGN1 and (Compton thin) AGN2 counts at 15 um. Methods: XMM observations of the ELAIS-S1 15um sample down to flux limits ~2x10^-15 erg cm^-2 s^-1 (2-10 keV band) were used. We classified as AGN all those MIR sources with a unabsorbed 2-10 keV X-ray luminosity higher that ~10^42 erg/s. Results: We find that at least about 13(+/-6) per cent of the previously classified starburst galaxies harbor an AGN. According to these figures, we provide an updated estimate of the counts of AGN1 and (Compton thin) AGN2 at 15 um. It turns out that at least 24% of the extragalactic sources brighter than 0.6 my at 15 um are AGN (~13% contribution to the extragalactic background produced at fluxes brighter than 0.6 mJy).

(Abridged) We present the results of the photometric and spectroscopic identification of 110 hard X-ray selected sources (6e-15<F(2-10 keV)<4e-14 cgs) from 5 additional XMM-Newton fields, nearly doubling the original HELLAS2XMM sample. We spectroscopically identified 59 new sources, bringing the completeness of the full HELLAS2XMM sample to ~70%. We found optical counterparts for 214 out of the 232 X-ray sources of the full sample down to R~25. We measure the flux and luminosity of the [OIII](5007) emission line for 59 such sources. Results. We use the full HELLAS2XMM and the CDF samples to estimate the LogN-LogS and the obscured QSO surface densities down to F(2-10 keV)~1e-14 and 1e-15 cgs. At these flux limits the fraction of X-ray selected obscured and unobscured QSO turns out to be similar. Since X-ray selection misses most Compton thick AGN, the number of obscured QSO may well outnumber that of unobscured QSOs. We find that hard X-ray selected AGNs with detected [OIII] emission span a wide range of L(2-10 keV)/L[OIII] with a median higher than that of optically selected AGNs, suggesting that optically selected samples are at least partly incomplete, and/or [OIII] emission is not a perfect isotropic indicator of the nuclear power. Differently from narrow-line AGNs in the same z interval, the 7 XBONGs candidates in the sample have L(2-10 keV)/L[OIII]> 1000, while their Lx, L(R) and Nh are similar. This suggests that while the central engine of narrow-line AGNs and XBONGs looks similar, the narrow-line region in XBONGs could be strongly inhibited or obscured.

We have surveyed with XMM-Newton the central ~0.6 deg2 region of the ELAIS-S1 field down to flux limits of ~5.5X10-16 cgs (0.5-2 keV, S band), ~2X10-15 cgs (2-10 keV, H band), and ~4X10-15 cgs (5-10 keV, HH band). We detect a total of 478 sources, 395 and 205 of which detected in the S and H bands respectively. We identified 7 clearly extended sources and estimated their redshift through X-ray spectral fits with thermal models. In four cases the redshift is consistent with z=0.4. We have computed the angular correlation function of the sources in the S and H bands, finding best fit correlation angles theta_0=5.2+/-3.8 arcsec and theta_0=12.8+/-7.8 arcsec respectively. A rough estimate of the present-day correlation length r_0 can be obtained inverting the Limber equation and assuming an appropriate redshift distribution dN/dz. The results range between 12.8 and 9.8 h-1 Mpc in the S band and between 17.9 and 13.4 h-1 Mpc in the H band, with 30-40% statistical errors, assuming either smooth redshift distributions or redshift distributions with spikes accounting for the presence of a structure at z=0.4. The relative density of the S band sources is higher near the clusters and groups at z~0.4 and extends toward East and toward South/West. This suggests that the structure is complex, with a size comparable to the full XMM-Newton field. Conversely, the highest relative source densities of the H band sources are located in the central-west region of the field.

We have determined the cosmological evolution of the density of AGN and of their Nh distribution as a function of the un-absorbed 2-10 keV luminosity up to redshift 4. We used the HELLAS2XMM sample combined with other published catalogs, yielding a total of 508 AGN. Our best fit is obtained with a luminosity-dependent density evolution (LDDE) model where low luminosity (Lx~10^43 erg s^-1) AGN peak at z~0.7, while high luminosity AGN (Lx>10^45 erg s^-1) peak at z~2.0. A pure luminosity evolution model (PLE) can instead be rejected. There is evidence that the fraction of absorbed (Nh>10^22 cm^-2) AGN decreases with the intrinsic X-ray luminosity, and increases with the redshift. Our best fit solution provides a good fit to the observed counts, the cosmic X-ray background, and to the observed fraction of absorbed AGN as a function of the flux in the 10^-15<S(2-10)<10^-10 erg s^-1 cm^-2 range. We find that the absorbed, high luminosity (Lx>10^44 erg s^-1) AGN have a density of 267 deg^-2 at fluxes S(2-10)>10^-15 erg s^-1 cm^-2. Using these results, we estimate a density of supermassive black holes in the local Universe of 3.2 (h^2_70) x 10^5 M_sol Mpc^-3, which is consistent with the recent measurements of the black hole mass function in the local galaxies.