Galactic outflows are known to consist of several gas phases, however, so far the connection between these multiple phases has been investigated little and only in a few objects. In this paper, we analyse MUSE/VLT data of 26 local (U)LIRGs and study their ionised and neutral atomic phases. We also include objects from the literature to obtain a total sample of 31 galaxies with spatially resolved multi-phase outflow information. We find that the ionized phase of the outflows has on average an electron density three times higher than the disc ($n_{\rm e, disc}$ $\sim$ 150 cm$^{-3}$ vs $n_{\rm e, outflow}$ $\sim$ 500 cm$^{-3}$), suggesting that cloud compression in the outflow is more important that cloud dissipation. We find that the difference in extinction between outflow and disc correlates with the outflow gas mass. Together with the analysis of the outflow velocities, this suggests that at least some of the outflows are associated with the ejection of dusty clouds from the disc. The presence of dust in outflows is relevant for potential formation of molecules inside them. We combine our data with millimetre data to investigate the molecular phase. We find that the molecular phase accounts for more than 60 $\%$ of the total mass outflow rate in most objects and this fraction is higher in AGN-dominated systems. The neutral atomic phase contributes of the order of 10 $\%$, while the ionized phase is negligible. The ionized-to-molecular mass outflow rate declines slightly with AGN luminosity, although with a large scatter.

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

We present ALMA band-7 data of the [CII] $\lambda157.74\,\mu{\rm m}$ emission line and underlying far-infrared (FIR) continuum for twelve luminous quasars at $z \simeq 4.8$, powered by fast-growing supermassive black holes (SMBHs). Our total sample consists of eighteen quasars, twelve of which are presented here for the first time. The new sources consists of six Herschel/SPIRE detected systems, which we define as "FIR-bright" sources, and six Herschel/SPIRE undetected systems, which we define as "FIR-faint" sources. We determine dust masses for the quasars hosts of $M_{dust} \le 0.2-25.0\times 10^8 M_{\odot}$, implying ISM gas masses comparable to the dynamical masses derived from the [CII] kinematics. It is found that on average the MgII line is blueshifted by $\sim 500\,{\rm km\,s}^{-1}$ with respect to the [CII] emission line, which is also observed when complementing our observations with data from the literature. We find that all of our "FIR-bright" subsample and most of the "FIR-faint" objects lie above the main sequence of star forming galaxies at $z \sim 5$. We detect companion sub-millimeter galaxies (SMGs) for two sources, both FIR-faint, with a range of projected distances of $\sim20-60$ kpc and with typical velocity shifts of $\left|\Delta v\right| \lesssim200\,{\rm km\,s}^{-1}$ from the quasar hosts. Of our total sample of eighteen quasars, 5/18 are found to have dust obscured starforming companions.

We present the discovery of PSO J030947.49+271757.31, the radio brightest (23.7 mJy at 1.4 GHz) active galactic nucleus (AGN) at z>6.0. It was selected by cross-matching the NRAO VLA Sky Survey and the Panoramic Survey Telescope and Rapid Response System PS1 databases and its high-z nature was confirmed by a dedicated spectroscopic observation at the Large Binocular Telescope. A pointed Neil Gehrels $Swift$ Observatory XRT observation allowed us to measure a flux of $\sim$3.4$\times$10$^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the [0.5-10] keV energy band, which also makes this object the X-ray brightest AGN ever observed at z>6.0. Its flat radio spectrum ($\alpha_{\nu r}$<0.5), very high radio loudness (R>10$^3$), and strong X-ray emission, compared to the optical, support the hypothesis of the blazar nature of this source. %i.e. a radio-loud (RL) AGN with the relativistic jet pointed toward us. Assuming that this is the only blazar at this redshift in the surveyed area of sky, we derive a space density of blazars at z$\sim$6 and with M$_{1450 \mbox{\AA}}$ < -25.1 of 5.5$^{+11.2}_{-4.6}$$\times$10$^{-3}$ Gpc$^{-3}$. From this number, and assuming a reasonable value of the bulk velocity of the jet ($\Gamma$=10), we can also infer a space density of the entire radio-loud AGN population at z$\sim$6 with the same optical/UV absolute magnitude of 1.10$^{+2.53}_{-0.91}$ Gpc$^{-3}$. Larger samples of blazars will be necessary to better constrain these estimates.

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

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.

To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity.

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)

Dynamical black hole mass measurements in some gas-rich galaxy mergers indicate that they are overmassive relative to their host galaxy properties. Overmassive black holes in these systems present a conflict with the standard progression of galaxy merger - quasar evolution; an alternative explanation is that a nuclear concentration of molecular gas driven inward by the merger is affecting these dynamical black hole mass estimates. We test for the presence of such gas near the two black holes in NGC 6240 using long-baseline ALMA Band 6 observations (beam size 0"06 $\times$ 0"03 or 30 pc$\times$15 pc). We find (4.2-9.8) $\times10^{7}$ M$_{\odot}$ and (1.2-7.7) $\times10^{8}$ M$_{\odot}$ of molecular gas within the resolution limit of the original black hole mass measurements for the north and south black holes, respectively. In the south nucleus, this measurement implies that 6-89% of the original black hole mass measurement actually comes from molecular gas, resolving the tension in the original black hole scaling relations. For the north, only 5% to 11% is coming from molecular gas, suggesting the north black hole is actually overmassive. Our analysis provides the first measurement of significant molecular gas masses contaminating dynamical black hole mass measurements. These high central molecular gas densities further present a challenge to theoretical black hole accretion prescriptions, which often assume accretion proceeds rapidly through the central 10 pc.

We present the discovery and properties of DESJ014132.4-542749.9 (DES0141-54), a new powerful radio-loud active galactic nucleus (AGN) in the early Universe (z=5.0). It was discovered by cross-matching the first data release of the Dark Energy Survey (DES DR1) with the Sidney University Molonglo Survey (SUMSS) radio catalog at 0.843 GHz. This object is the first radio-loud AGN at high redshift discovered in the DES. The radio properties of DES0141-54, namely its very large radio-loudness (R>10$^{4}$), the high radio luminosity (L$_{0.8 GHz}$=1.73$\times$10$^{28}$ W Hz$^{-1}$), and the flatness of the radio spectrum ($\alpha$=0.35) up to very high frequencies (120 GHz in the source's rest frame), classify this object as a blazar, meaning, a radio-loud AGN observed along the relativistic jet axis. However, the X--ray luminosity of DESJ0141-54 is much lower compared to those of the high redshift (z$\geq$4.5) blazars discovered so far. Moreover its X-ray-to-radio luminosity ratio (log($\frac{L_{[0.5-10]keV}}{L_{1.4GHz}}$)=9.96$\pm$0.30 Hz) is small also when compared to lower redshift blazars: only 2\% of the low-z population has a similar ratio. By modeling the spectral energy distribution we found that this peculiar X--ray weakness and the powerful radio emission could be related to a particularly high value of the magnetic field. Finally, the mass of the central black hole is relatively small (M$_{BH}$ = 3-8 $\times$10$^8$ M$_{\odot}$) compared to other confirmed blazars at similar redshift, making DES0141-54 the radio-loud AGN that host the smallest supermassive black hole ever discovered at z$\geq$5.

We report the first Atacama large millimeter/submillimeter array observations of MCG-03-58-007, a local ($z=0.03236\pm0.00002$, this work) AGN ($L_{AGN}\sim10^{45}~\rm erg~s^{-1}$), hosting a powerful X-ray ultra-fast ($v=0.1c$) outflow (UFO). The CO(1-0) line emission is observed across $\sim18\,$kpc scales with a resolution of $\sim 1\,\rm kpc$. About 78\% of the CO(1-0) luminosity traces a galaxy-size rotating disk. However, after subtracting the emission due to such rotating disk, we detect with a S/N=20 a residual emission in the central $\sim 4\,$kpc. Such residuals may trace a low velocity ($v_{LOS}=170\,\rm km\,s^{-1}$) outflow. We compare the momentum rate and kinetic power of such putative molecular outflow with that of the X-ray UFO and find $\dot{P}_{out}/\dot{P}_{UFO}=0.3\pm0.2$ and $\dot{E}_{mol}/\dot{E}_{UFO}\sim4\cdot10^{-3}$. This result is at odds with the energy-conserving scenario suggested by the large momentum boosts measured in some other molecular outflows. An alternative interpretation of the residual CO emission would be a compact rotating structure, distinct from the main disk, which would be a factor of $\sim10-100$ more extended and massive than typical circumnuclear disks revealed in Seyferts. However, in both scenarios, our results rule out the hypothesis of a momentum-boosted molecular outflow in this AGN, despite the presence of a powerful X-ray UFO. [Abridged]

The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.

State-of-the-art rest-frame UV and FIR photometric and spectroscopic observations are now pushing the redshift frontiers of galaxy formation studies up to $z\sim9-11$ and beyond. Recent HST observations unveiled the presence of a star-forming galaxy exhibiting the Lyman break at $\lambda_{\rm obs}=1.47\pm0.01$ $\mu$m, i.e., a $z=11.09^{+0.08}_{-0.12}$ galaxy with a stellar mass of $\sim10^9 M_\odot$, demonstrating that galaxy build-up was well underway early in the epoch of reionization (EoR) at $z>10$. Targeted spectroscopy of a lensed Lyman break galaxy uncovers the earliest metals known to date up to $z=9.1096\pm0.0006$ by detecting the bright [OIII] 88~$\mu$m nebular line, indicating the onset of star formation 250 million years after the Big Bang, i.e., corresponding to a redshift of $z\sim15$. These latest findings lead us to a number of key questions: How and when metal enrichment happened in the EoR? What was the nature of the earliest-epoch star-forming galaxies at $z=10-15$? What was the spatial distribution of such galaxies, and what was the relation to the putative large-scale ionization bubbles during the EoR? What were the dark-halo masses of such earliest-epoch star-forming galaxies? To address all these questions, we need to uncover a statistically large number of $z=10-15$ galaxies in the pre-reionization era. Here we argue two possible pathways: (1) a wide-area, sensitive blind spectroscopic survey of [OIII] 88 $\mu$m line-emitting galaxies at submillimeter wavelengths, and (2) an ultra-wide-area, high-cadence photometric survey of transient sources at radio-to-(sub)millimeter wavelengths, together with the immediate follow-up spectroscopy with an ultra-wide-band spectrograph, to catch the pop-III $\gamma$-ray bursts.

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.

Based on new near-infrared spectroscopic data from the instrument EMIR on the 10.4 m Gran Telescopio Canarias (GTC) we report the presence of an ionized and warm molecular outflow in the luminous type-2 quasar J150904.22+043441.8 (z = 0.1118). The ionized outflow is faster than its molecular counterpart, although the outflow sizes that we derive for them are consistent within the errors (1.34$\pm$0.18 kpc and 1.46$\pm$0.20 kpc respectively). We use these radii, the broad emission-line luminosities and in the case of the ionized outflow, the density calculated from the trans-auroral [OII] and [SII] lines, to derive mass outflow rates and kinetic coupling efficiencies. Whilst the ionized and warm molecular outflows represent a small fraction of the AGN power ($\leq$0.033% and 0.0001% of L$_{bol}$ respectively), the total molecular outflow, whose mass is estimated from an assumed warm-to-cold gas mass ratio of 6$\times10^{-5}$, has a kinetic coupling efficiency of $\sim$1.7%L$_{bol}$. Despite the large uncertainty, this molecular outflow represents a significant fraction of L$_{bol}$ and it could potentially have a significant impact on the host galaxy. In addition, the quasar spectrum reveals bright and patchy narrow Pa$\alpha$ emission extending out to 4 arcsec (8 kpc) South-East and North-West from the active nucleus.

Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black hole growth lies in a thoroughly-studied census of dual AGN across cosmic history, which will be enabled by next-generation observational capabilities, theoretical advances, and simulations. This white paper outlines the key questions in galaxy mergers, dual and offset AGN, and proposes multiwavelength solutions using future high-resolution observatories in the X-rays (AXIS, Lynx), near and mid-infrared (30 meter class telescopes, JWST), and submillimeter (ALMA).

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

We study molecular outflows in a sample of 45 local galaxies, both star forming and AGN, primarily by using CO data from the ALMA archive and from the literature. For a subsample we also compare the molecular outflow with the ionized and neutral atomic phases. We infer an empirical analytical function relating the outflow rate simultaneously to the SFR, $L_{\rm AGN}$, and galaxy stellar mass; this relation is much tighter than the relations with the individual quantities. The outflow kinetic power shows a larger scatter than in previous, more biased studies, spanning from 0.1 to 5 per cent of $L_{\rm AGN}$, while the momentum rate ranges from 1 to 30 times $L_{\rm AGN}/c$, indicating that these outflows can be both energy-driven, but with a broad range of coupling efficiencies with the ISM, and radiation pressure-driven. For about 10 per cent of the objects the outflow energetics significantly exceed the maximum theoretical values; we interpret these as 'fossil outflows' resulting from activity of a past strong AGN, which has now faded. We estimate that, in the stellar mass range probed here ($>$ 10$^{10}~\rm M_{\odot}$), less than 5 per cent of the outflowing gas escapes the galaxy. The molecular gas depletion time associated with the outflow can be as short as a few million years in powerful AGN, however, the total gas (H$_2$+HI) depletion times are much longer. Altogether, our findings suggest that even AGN-driven outflows might be relatively ineffective in clearing galaxies of their entire gas content, although they are likely capable of clearing and quenching the central region.

High redshift blazars are an important class of Active Galactic Nuclei (AGN) that can provide an independent estimate of the supermassive black-hole mass function in high redshift radio-loud AGN without the bias due to absorption along the line-of-sight. Using the Cosmic Lens All Sky Survey (CLASS) we built a complete radio flux-limited sample of high redshift (z>4) blazars suitable for statistical studies. By combining dedicated optical observations and the SDSS spectroscopic database, we obtained a sample of 26 blazar candidates with a spectroscopic redshift above 4. On the basis of their radio spectrum we distinguish between blazars and QSO with a Gigahertz Peaked Spectrum (GPS) like spectrum. Out of the 18 confirmed blazars 14 constitute a completely identified, flux-limited sample down to a magnitude of 21 (AB). Using this complete sample we derive a space density of blazars with 4<z<5.5 of rho=0.13 (+0.05,-0.03) Gpc^-3. This is the first actual estimate of the blazar space density in this range of redshift. This value is in good agreement with the extrapolation of the luminosity function and cosmological evolution based on a sample of flat-spectrum radio quasars selected at lower redshifts and it is consistent with a cosmological evolution peaking at z$\sim$2 similar to radio-quiet QSO. We do not confirm, instead, the presence of a peak at z~4 in the space density evolution, recently suggested using an X-ray selected sample of blazars. It is possible that this extreme peak of the evolution is present only among the most luminous blazars.

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 results of a detailed analysis of a deep simultaneous $130\,\rm ks$ \textit{XMM-Newton & NuSTAR} observation of the nearby ($z=0.0315$) and bright ($L_{\rm bol}\sim3\times10^{45}\,\rm erg\,s^{-1}$) starburst-AGN Seyfert\,2 system: MCG--03--58--007. From the broadband fitting we show that most of the obscuration needs to be modeled with a toroidal type reprocessor such as \texttt{MYTorus} \citep{MurphyYaqoob09}. Nonetheless the signature of a powerful disc-wind is still apparent at higher energies and the observed rapid short-term X-ray spectral variability is more likely caused by a variable zone of highly ionized fast wind rather than by a neutral clumpy medium. We also detect X-ray emission from larger scale gas as seen from the presence of several soft narrow emission lines in the RGS, originating from a contribution of a weak star forming activity together with a dominant photoionized component from the AGN.

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 study the molecular gas properties of two star-forming galaxies separated by 6 kpc in the projected space and belonging to a galaxy cluster selected from the Irac Shallow Cluster Survey, at a redshift $z=1.2$, i.e., $\sim2$ Gyr after the cosmic star formation density peak. This work describes the first CO detection from $1<z<1.4$ star forming cluster galaxies with no reported clear evidence of AGN. We exploit observations taken with the NOEMA interferometer at $\sim3$ mm to detect CO(2-1) line emission from the two selected galaxies, unresolved by our observations. Based on the CO(2-1) spectrum we estimate a total molecular gas mass $M({\rm H_2})=(2.2^{+0.5}_{-0.4})\times10^{10}$ $M_\odot$ and dust mass $M_{\rm dust}<4.2\times10^8~M_\odot$ for the two blended sources. The two galaxies have similar stellar masses and a large relative velocity of $\sim$400 km/s estimated from the CO(2-1) line width. These findings tend to privilege a scenario where both sources contribute to the observed CO(2-1). By using the archival Spitzer MIPS flux at 24$\mu$m we estimate an ${\rm SFR(24\mu m)}=(28^{+12}_{-8})~M_\odot$/yr for each of the two galaxies. Assuming that the two sources equally contribute to the observe CO(2-1) our analysis yields a depletion time scale $\tau_{\rm dep}=(3.9^{+1.4}_{-1.8})\times10^8$ yr, and a molecular gas to stellar mass ratio $0.17\pm0.13$ for each of two sources, separately. Our results are in overall agreement with those of other distant cluster galaxies. The two target galaxies have molecular gas mass and depletion time that are marginally compatible with, but smaller than those of main sequence field galaxies, suggesting that the molecular gas has not been refueled enough. Higher resolution and higher frequency observations will enable us to spatially resolve the two sources and possibly distinguish between different gas processing mechanisms.

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 discovery of a new candidate for a fast disk wind, in the nearby Seyfert 2 galaxy MCG-03-58-007. This wind is discovered in a deep Suzaku observation that was performed in 2010. Overall the X-ray spectrum of MCG-03-58-007 is highly absorbed by a neutral column density of NH~10^23 cm^-2, in agreement with the optical classification as a type 2 AGN. In addition, this observation unveiled the presence of two deep absorption troughs at E = 7.4 +- 0.1 keV and E = 8.5 +- 0.2 keV. If associated with blue-shifted FeXXVI, these features can be explained with the presence of two highly ionised (log \xi/(erg cm/s)~ 5.5) and high column density (NH~5-8 x 10^23cm^-2) outflowing absorbers with v_out1~ -0.1c and v_out2~ -0.2c. The disk wind detected during this observation is most likely launched from within a few hundreds gravitational radii from the central black and has a kinetic output that matches the prescription for significant feedback. The presence of the lower velocity component of the disk wind is independently confirmed by the analysis of a follow-up XMM-Newton & NuSTAR observation. A faster (v_out~ -0.35 c) component of the wind is also seen in this second observation. During this observation we also witnessed an occultation event lasting \Delta t ~ 120 ksec, which we ascribe to an increase of the opacity of the disk wind (\Delta NH~1.4x10^24 cm^-2). Our interpretation is that the slow zone (v_out~ -0.1c) of the wind is the most stable but inhomogeneous component, while the faster zones could be associated with two different inner streamlines of the wind.

Dual/binary Supermassive Black Hole (SMBH) systems are the inevitable consequence of the current Lambda Cold Dark Matter cosmological paradigm. In this context, we discuss here the properties of MCG+11-11-032, a local (z=0.0362) Seyfert 2 galaxy. This source was proposed as a dual AGN candidate on the basis of the presence of double-peaked [OIII] emission lines in its optical spectrum. MCG+11-11-032 is also an X-ray variable source and was observed several times by the Swift X-ray Telescope (XRT) on time scales from days to years. In this work, we analyze the SDSS-DR13 spectrum and find evidence for double-peaked profiles in all the strongest narrow emission lines. We also study the XRT light curve and unveil the presence of an alternating behavior of the intrinsic 0.3-10 keV flux, while the 123-month Swift BAT light curve supports the presence of almost regular peaks and dips almost every 25 months. In addition, the XRT spectrum suggests for the presence of two narrow emission lines with rest-frame energies of E~6.16 keV and E~6.56 keV. Although by considering only the optical emission lines, different physical mechanisms may be invoked to explain the kinematical properties, the X-ray results are most naturally explained by the presence of a binary SMBH in the center of this source. In particular, we evidence a remarkable agreement between the putative SMBH pair orbital velocity derived from the BAT light curve and the velocity offset derived by the rest-frame Delta_E between the two X-ray line peaks in the XRT spectrum (i.e. Delta_v~0.06c).

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.

Observations and simulations show that outflows in active galactic nuclei (AGN) contain gas in different phases. To understand their true impact on galaxy evolution, we advocate consistent and unbiased investigation of these multiphase winds in large AGN samples.

We present a long-term, multi-wavelength project to understand the epoch of fastest growth of the most massive black holes by using a sample of 40 luminous quasars at z~4.8. These quasars have rather uniform properties, with typical accretion rates and black hole masses of L/L_Edd~0.7 and M_BH~10^9 M_sol. The sample consists of "FIR-bright" sources with a previous Herschel/SPIRE detection, suggesting SFR > 1000 M_sol/yr, as well as of "FIR-faint" sources for which Herschel stacking analysis implies a typical SFR of ~400 M_sol/yr. Six of the quasars have been observed by ALMA in [CII]{\lambda}157.74 micron line emission and adjacent rest-frame 150 micron continuum, to study the dusty cold ISM. ALMA detected companion, spectroscopically confirmed sub-mm galaxies (SMGs) for three sources - one FIR-bright and two FIR-faint. The companions are separated by ~14-45 kpc from the quasar hosts, and we interpret them as major galaxy interactions. Our ALMA data therefore clearly support the idea that major mergers may be important drivers for rapid, early SMBH growth. However, the fact that not all high-SFR quasar hosts are accompanied by interacting SMGs, and their ordered gas kinematics observed by ALMA, suggest that other processes may be fueling these systems. Our analysis thus demonstrates the diversity of host galaxy properties and gas accretion mechanisms associated with early and rapid SMBH growth.

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 investigate quasar outflows at $z \geq 6$ by performing zoom-in cosmological hydrodynamical simulations. By employing the SPH code GADGET-3, we zoom in the $2 R_{200}$ region around a $2 \times 10^{12} M_{\odot}$ halo at $z = 6$, inside a $(500 ~ {\rm Mpc})^3$ comoving volume. We compare the results of our AGN runs with a control simulation in which only stellar/SN feedback is considered. Seeding $10^5 M_{\odot}$ BHs at the centers of $10^{9} M_{\odot}$ halos, we find the following results. BHs accrete gas at the Eddington rate over $z = 9 - 6$. At $z = 6$, our most-massive BH has grown to $M_{\rm BH} = 4 \times 10^9 M_{\odot}$. Fast ($v_{r} > 1000$ km/s), powerful ($\dot{M}_{\rm out} \sim 2000 M_{\odot}$/yr) outflows of shock-heated low-density gas form at $z \sim 7$, and propagate up to hundreds kpc. Star-formation is quenched over $z = 8 - 6$, and the total SFR (SFR surface density near the galaxy center) is reduced by a factor of $5$ ($1000$). We analyse the relative contribution of multiple physical process: (i) disrupting cosmic filamentary cold gas inflows, (ii) reducing central gas density, (iii) ejecting gas outside the galaxy; and find that AGN feedback has the following effects at $z = 6$. The inflowing gas mass fraction is reduced by $\sim 12 \%$, the high-density gas fraction is lowered by $\sim 13 \%$, and $\sim 20 \%$ of the gas outflows at a speed larger than the escape velocity ($500$ km/s). We conclude that quasar-host galaxies at $z \geq 6$ are accreting non-negligible amount of cosmic gas, nevertheless AGN feedback quenches their star formation dominantly by powerful outflows ejecting gas out of the host galaxy halo.

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.

We present the final data release of the APEX low-redshift legacy survey for molecular gas (ALLSMOG), comprising CO(2-1) emission line observations of 88 nearby, low-mass (10^8.5<M* [M_Sun]<10^10) star-forming galaxies carried out with the 230 GHz APEX-1 receiver on the APEX telescope. The main goal of ALLSMOG is to probe the molecular gas content of more typical and lower stellar mass galaxies than have been studied by previous CO surveys. We also present IRAM 30m observations of the CO(1-0) and CO(2-1) emission lines in nine galaxies aimed at increasing the M*<10^9 M_Sun sample size. In this paper we describe the observations, data reduction and analysis methods and we present the final CO spectra together with archival HI 21cm line observations for the entire sample of 97 galaxies. At the sensitivity limit of ALLSMOG, we register a total CO detection rate of 47%. Galaxies with higher M*, SFR, nebular extinction (A_V), gas-phase metallicity (O/H), and HI gas mass have systematically higher CO detection rates. In particular, the parameter according to which CO detections and non-detections show the strongest statistical differences is the gas-phase metallicity, for any of the five metallicity calibrations examined in this work. We investigate scaling relations between the CO(1-0) line luminosity and galaxy-averaged properties using ALLSMOG and a sub-sample of COLD GASS for a total of 185 sources that probe the local main sequence (MS) of star-forming galaxies and its +-0.3 dex intrinsic scatter from M* = 10^8.5 M_Sun to M* = 10^11 M_Sun. L'_{CO(1-0)} is most strongly correlated with the SFR, but the correlation with M* is closer to linear and almost comparably tight. The relation between L'_{CO(1-0)} and metallicity is the steepest one, although deeper CO observations of galaxies with A_V<0.5 mag may reveal an as much steep correlation with A_V. [abridged]

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 new ALMA band-7 data for a sample of six luminous quasars at z~4.8, powered by fast-growing supermassive black holes (SMBHs) with rather uniform properties: the typical accretion rates and black hole masses are L/L_Edd~0.7 and M_BH~10^9 M_sol. Our sample consists of three "FIR-bright" sources, which were individually detected in previous Herschel/SPIRE observations, with star formation rates of SFR>1000 M_sol/yr, and three "FIR-faint" sources for which Herschel stacking analysis implies a typical SFR of ~400 M_sol/yr. The dusty interstellar medium in the hosts of all six quasars is clearly detected in the ALMA data and resolved on scales of ~2 kpc, in both continuum (\lambda_rest ~150 um) and [C II]157.74 um line emission. The continuum emission is in good agreement with the expectations from the Herschel data, confirming the intense SF activity in the quasar hosts. Importantly, we detect companion sub-millimeter galaxies (SMGs) for three sources -- one FIR-bright and two FIR-faint, separated by ~14-45 kpc and <450 km/s from the quasar hosts. The [C II]-based dynamical mass estimates for the interacting SMGs are within a factor of ~3 of the quasar hosts' masses, while the continuum emission implies SFR(quasar)~(2-11)xSFR(SMG). Our ALMA data therefore clearly support the idea that major mergers are important drivers for rapid early SMBH growth. However, the fact that not all high-SFR quasar hosts are accompanied by interacting SMGs and the gas kinematics as observed by ALMA suggest that other processes may be fueling these systems. Our analysis thus demonstrates the diversity of host galaxy properties and gas accretion mechanisms associated with early and rapid SMBH growth.

In this work we present an analysis of the behaviour of galaxies in a four-dimensional parameter space defined by stellar mass, metallicity, star formation rate, and molecular gas mass. We analyse a combined sample of 227 galaxies, which draws from a number of surveys across the redshift range 0 < z < 2 (> 90% of the sample at z~0), and covers > 3 decades in stellar mass.Using Principle Component Analysis, we demonstrate that galaxies in our sample lie on a 2-dimensional plane within this 4D parameter space, indicative of galaxies that exist in an equilibrium between gas inflow and outflow. Furthermore, we find that the metallicity of galaxies depends only on stellar mass and molecular gas mass. In other words, gas-phase metallicity has a negligible dependence on star formation rate, once the correlated effect of molecular gas content is accounted for. The well-known `fundamental metallicity relation', which describes a close and tight relationship between metallicity and SFR (at fixed stellar mass) is therefore entirely a by-product of the underlying physical relationship with molecular gas mass (via the Schmidt-Kennicutt relation).

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 investigate the properties of star formation-driven outflows by using a large spectroscopic sample of ~160,000 local "normal" star forming galaxies, drawn from the SDSS, spanning a wide range of star formation rates and stellar masses. The galaxy sample is divided into a fine grid of bins in the M_*-SFR parameter space, for each of which we produce a composite spectrum by stacking together the SDSS spectra of the galaxies contained in that bin. We exploit the high signal-to-noise of the stacked spectra to study the emergence of faint features of optical emission lines that may trace galactic outflows and would otherwise be too faint to detect in individual galaxy spectra. We adopt a novel approach that relies on the comparison between the line-of-sight velocity distribution (LoSVD) of the ionised gas (as traced by the [OIII]5007 and Halpha+[NII]6548,6583 emission lines) and the LoSVD of the stars, which are used as a reference tracing virial motions. Significant deviations of the gas kinematics from the stellar kinematics in the high velocity tail of the LoSVDs are interpreted as a signature of outflows. Our results suggest that the incidence of ionised outflows increases with SFR and sSFR. The outflow velocity (v_out) correlates tightly with the SFR for SFR>1 M_Sun/yr, whereas at lower SFRs the dependence of v_out on SFR is nearly flat. The outflow velocity, although with a much larger scatter, increases also with the stellar velocity dispersion, and we infer velocities as high as v_out~(6-8)*sigma_stars. Strikingly, we detect the signature of ionised outflows only in galaxies located above the main sequence (MS) of star forming galaxies in the M_*-SFR diagram, and the incidence of such outflows increases sharply with the offset from the MS. Our complementary analysis of the stellar kinematics reveals the presence of blue asymmetries of a few 10 km/s in the stellar LoSVDs. [abridged]

There has been much recent work dedicated to exploring secondary correlations in the mass-metallicity relation, with significant dependence on both the star formation rate and HI content being demonstrated. Previously, a paucity of molecular gas data (combined with sample selection bias) hampered the investigation of any such relation with molecular gas content. In this work, we assemble a sample of 221 galaxies from a variety of surveys in the redshift range 0 < z < 2, to explore the connection between molecular gas content and metallicity. We explore the effect of gas mass on the mass-metallicity relation, finding that the offset from the relation is negatively correlated against both molecular and total gas mass. We then employ a principle component analysis technique to explore secondary dependences in the mass-metallicity relation, finding that the secondary dependence with gas mass is significantly stronger than with star formation rate, and as such the underlying `Fundamental Metallicity Relation' is between stellar mass, metallicity, and gas mass. In particular, the metallicity dependence on SFR is simply a byproduct of the dependence on the molecular gas content, via the Schmidt-Kennicutt relation. Finally, we note that our principle component analysis finds essentially no connection between gas-phase metallicity and the efficiency of star formation.

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.

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.

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.

We present ALLSMOG, the APEX Low-redshift Legacy Survey for MOlecular Gas. ALLSMOG is a survey designed to observe the CO(2-1) emission line with the APEX telescope, in a sample of local galaxies (0.01 < z < 0.03), with stellar masses in the range 8.5 < log(M*/Msun) < 10. This paper is a data release and initial analysis of the first two semesters of observations, consisting of 42 galaxies observed in CO(2-1). By combining these new CO(2-1) emission line data with archival HI data and SDSS optical spectroscopy, we compile a sample of low-mass galaxies with well defined molecular gas masses, atomic gas masses, and gas-phase metallicities. We explore scaling relations of gas fraction and gas consumption timescale, and test the extent to which our findings are dependent on a varying CO/H2 conversion factor. We find an increase in the H2/HI mass ratio with stellar mass which closely matches semi-analytic predictions. We find a mean molecular gas fraction for ALLSMOG galaxies of MH2/M* = (0.09 - 0.13), which decreases with stellar mass. We measure a mean molecular gas consumption timescale for ALLSMOG galaxies of 0.4 - 0.7 Gyr. We also confirm the non-universality of the molecular gas consumption timescale, which varies (with stellar mass) from ~100 Myr to ~2 Gyr. Importantly, we find that the trends in the H2/HI mass ratio, gas fraction, and the non-universal molecular gas consumption timescale are all robust to a range of recent metallicity-dependent CO/H2 conversion factors.

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.

Recent results have suggested that the well known mass-metallicity relation has a strong dependence on the star formation rate, to the extent that a three dimensional `fundamental metallicity relation' exists which links the three parameters with minimal scatter. In this work, we use a sample of 4253 local galaxies observed in atomic hydrogen from the ALFALFA survey to demonstrate, for the first time, that a similar fundamental relation (the HI-FMR) also exists between stellar mass, gas-phase metallicity, and HI mass. This latter relation is likely more fundamental, driving the relation between metallicity, SFR and mass. At intermediate masses, the behaviour of the gas fundamental metallicity relation is very similar to that expressed via the star formation rate. However, we find that the dependence of metallicity on HI content persists to the highest stellar masses, in contrast to the `saturation' of metallicity with SFR. It is interesting to note that the dispersion of the relation is very low at intermediate stellar masses (9< log(M*/Msun) <11), suggesting that in this range galaxies evolve smoothy, in an equilibrium between gas inflow, outflow and star formation. At high and low stellar masses, the scatter of the relation is significantly higher, suggesting that merging events and/or stochastic accretion and star formation may drive galaxies outside the relation. We also assemble a sample of galaxies observed in CO. However, due to a small sample size, strong selection bias, and the influence of a metallicity-dependent CO/H2 conversion factor, the data are insufficient to test any influence of molecular gas on metallicity.

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.

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.

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.