We study quenching in seven green valley galaxies on kpc scales by resolving their molecular gas content using \textsuperscript{12}CO(1-0) observations obtained with NOEMA and ALMA, and their star-formation rate using spatially resolved optical spectroscopy from the MaNGA survey. We perform radial stacking of both datasets to increase the sensitivity to molecular gas and star formation, thereby avoiding biases against strongly quenched regions. We find that both spatially resolved gas fraction ($\rm {f_{gas}}$) and star formation efficiency ($\rm {SFE}$) are responsible for quenching green valley galaxies at all radii: both quantities are suppressed with respect to typical star-forming regions. $\rm {f_{gas}}$ and $\rm {SFE}$ have roughly equal influence in quenching the outer disc. We are, however, unable to identify the dominant mechanism in the strongly quenched central regions. We find that $\rm {f_{gas}}$ is reduced by $\rm \sim 1~dex$ in the central regions, but the star formation rate is too low to be measured, leading to upper limits for the $\rm {SFE}$. Moving from the outer disc to central regions, the reduction in $\rm {f_{gas}}$ is driven by an increasing $\rm \Sigma_{\star}$ profile rather than a decreasing $\rm \Sigma_{H_{2}}$ profile. The reduced $\rm {f_{gas}}$ may therefore be caused by a decrease in the gas supply rather than molecular gas ejection mechanisms, such as winds driven by active galactic nuclei. We warn more generally that studies investigating $\rm {f_{gas}}$ may be deceiving in inferring the cause of quenching, particularly in the central (bulge-dominated) regions of galaxies.

ALMA observations have revealed the presence of dust in galaxies in the Epoch of Reionization (redshift $z>6$). However, the dust temperature, $T_d$, remains unconstrained, and this introduces large uncertainties, particularly in the dust mass determinations. Using an analytical and physically-motivated model, we show that dust in high-$z$, star-forming giant molecular clouds (GMC), largely dominating the observed far-infrared luminosity, is warmer ($T_d > 60\ \mathrm{K}$) than locally. This is due to the more compact GMC structure induced by the higher gas pressure and turbulence characterizing early galaxies. The compactness also delays GMC dispersal by stellar feedback, thus $\sim 40\%$ of the total UV radiation emitted by newly born stars remains obscured. A higher $T_d$ has additional implications: it (a) reduces the tension between local and high-$z$ IRX-$\beta$ relation, (b) alleviates the problem of the uncomfortably large dust masses deduced from observations of some EoR galaxies.

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.

ALMA observations have revealed that [CII] 158$\mu$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($\Delta^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations.

We present ALMA observations of a merging system at z ~ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII] 158 micron and far-infrared continuum data with multi-wavelength ancillary data we find that the system is composed of two massive (Mstar >~ 10^10 Msun) star-forming galaxies experiencing a major merger (stellar mass ratio r_mass ~ 0.9) at close spatial (~13 kpc; projected) and velocity (delta_v < 300 km/s) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger-scale protocluster environment and is coincident to one of its overdensity peaks. ALMA reveals also the presence of [CII] emission arising from a circumgalactic gas structure, extending out to ~30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterized by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-z galaxies, and thus play a key role in the galaxy baryon cycle at early epochs.

Far-infrared (FIR) emission lines are a powerful tool to investigate the properties of the interstellar medium, especially in high-redshift galaxies, where ALMA observations have provided unprecedented information. Interpreting such data with state-of-the-art cosmological simulations post-processed with CLOUDY, has provided insights on the internal structure and gas dynamics of these systems. However, no detailed investigation of the consistency and uncertainties of this kind of analysis has been performed to date. Here, we compare different approaches to estimate FIR line emission from state-of-the-art cosmological simulations, either with CLOUDY or with on-the-fly non-equilibrium chemistry. We find that [CII]$_{158\mu}$ predictions are robust. [OI] emission lines are instead model-dependent, as these lines are strongly affected by the thermodynamic state of the gas and non-equilibrium photoionisation effects. For the same reasons, [OI] lines represent an excellent ISM diagnostic. Future observations targeting these lines will be also crucial to constrain models.

We present a novel method to simultaneously characterise the star formation law and the interstellar medium properties of galaxies in the Epoch of Reionization (EoR) through the combination of [CII] 158$\mu$m (and its known relation with star formation rate) and CIII]$\lambda$1909{\AA} emission line data. The method, based on a Markov Chain Monte Carlo algorithm, allows to determine the target galaxy average density, $n$, gas metallicity, $Z$, and burstiness parameter, $\kappa_s$, quantifying deviations from the Kennicutt-Schmidt relation. As an application, we consider COS-3018 (z=6.854), the only EoR Lyman Break Galaxy so far detected in both [CII] and CIII]. We show that COS-3018 is a moderate starburst ($\kappa_s \approx 3$), with $Z\approx 0.4\, Z_{\odot}$, and $n \approx 500\, {\rm cm^{-3}}$. Our method will be optimally applied to joint ALMA and JWST targets.

This work aims at studying the $M_{BH}-M_{dyn}$ relation of a sample of $2<z<7$ quasars by constraining their host galaxy masses through a full kinematical modelling of the cold gas kinematics thus avoiding all possible biases and effects introduced by the rough virial estimates usually adopted so far. For this purpose we retrieved public observations of $72$ quasar host galaxies observed in ${\rm [CII]_{158\mu m}}$ or ${\rm CO}$ transitions with the Atacama Large Millimeter Array (ALMA). We then selected those quasars whose line emission is spatially resolved and performed a kinematic analysis on ALMA observations. We estimated the dynamical mass of the systems by modelling the gas kinematics with a rotating disc taking into account geometrical and instrumental effects. Our dynamical mass estimates, combined with $M_{BH}$ obtained from literature and our own new ${\rm CIV}\lambda1550$ observations have allowed us to investigate the $ M_{BH}/M_{dyn}$ in the early Universe. Overall we obtained a sample of $10$ quasars at $z\sim2-7$ in which line emission is detected with high SNR ($> 5-10$) and the gas kinematics is spatially resolved and dominated by ordered rotation. The estimated dynamical masses place $6$ out of $10$ quasars above the local relation yielding to a $M_{BH}/M_{dyn}$ ratios $\sim10\times$ higher than those estimated in low-$z$ galaxies. On the other hand we found that $4$ quasars at $z\sim 4-6$ have dynamical-to-BH mass ratios consistent with what is observed in early type galaxies in the local Universe.

While massive black holes (MBHs) are known to inhabit all massive galaxies, their ubiquitous presence in dwarf galaxies has not been confirmed yet, with only a limited number of sources detected so far. Recently, some studies proposed infrared emission as an alternative way to identify MBHs in dwarfs, based on a similar approach usually applied to quasars. In this study, by accurately combining optical and infrared data taking into account resolution effects and source overlapping, we investigate in detail the possible limitations of this approach with current ground-based facilities, finding a quite low ($\sim$0.4 per cent) fraction of active MBH in dwarfs that are luminous in mid-infrared, consistent with several previous results. Our results suggest that the infrared selection is strongly affected by several limitations that make the identification of MBHs in dwarf galaxies currently prohibitive, especially because of the very poor resolution compared to optical surveys, and the likely contamination by nearby sources, although we find a few good candidates worth further follow-ups. Optical, X-ray and radio observations, therefore, still represent the most secure way to search for MBH in dwarfs.

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 report the discovery of a high-redshift, massive molecular outflow in the starburst galaxy SPT0346-52 ($z=5.656$) via the detected absorption of high-excitation water transitions (H$_2$O $4_{2,3}-4_{1,4}$ and H$_2$O $3_{3,0}-3_{2,1}$) with the Atacama Large Millimeter/submillimeter Array (ALMA). The host galaxy is one of the most powerful starburst galaxies at high redshift (star formation rate; SFR $\sim3,600$M$_\odot$year$^{-1}$), with an extremely compact ($\sim320$pc) star formation region and a star formation rate surface density ($\Sigma_{SFR}\sim5,500$M$_{\odot}~$year$^{-1}~$kpc$^{-2}$) five times higher than `maximum' (i.e. Eddington-limited) starbursts, implying a highly transient phase. The estimated outflow rate is $\sim500$M$_{\odot}$year$^{-1}$, which is much lower than the SFR, implying that in this extreme starburst the outflow capabilities saturate and the outflow is no longer capable of regulating star formation, resulting in a runaway process in which star formation will use up all available gas in less than 30Myr. Finally, while previous kinematic investigations of this source revealed possible evidence for an ongoing major merger, the coincidence of the hyper-compact starburst and high-excitation water absorption indicates that this is a single starburst galaxy surrounded by a disc.

We study the impact of active galactic nuclei (AGN) on the CO Spectral Line Energy Distribution (SLED) of high-$z$ galaxies. In particular, we want to assess if the CO SLED can be used as a probe of AGN activity. To this purpose, we develop a semi-analytical model that takes into account the radiative transfer and the clumpy structure of giant molecular clouds where the CO lines are excited, their distribution in the galaxy disk, and the torus obscuration of the AGN radiation. We study the joint effect on the CO SLED excitation of (i) the X-ray luminosity from the AGN ($L_{X}$), (ii) the size of the molecular disk, (iii) the inclination angle ($\Omega$) of the torus with respect to the molecular disk, and (iv) the GMC mean density. We also discuss the possibility of an enhanced Cosmic Ray Ionization Rate (CRIR). We find that the X-ray Dominated Region (XDR) generated by the AGN in every case enhances the CO SLED for $J>5$, with increasing excitation of high-$J$ CO lines for increasing X-ray luminosity. Because high-$z$ galaxies are compact, the XDR region typically encloses the whole disk, thus its effect can be more important with respect to lower redshift objects. The impact of the XDR can be disentangled from an enhanced CRIR either if $L_X>10^{44} \rm \,erg\, s^{-1}$, or if $\Omega \geq 60^{\circ}$. We finally provide predictions on the CO(7-6)/[CII] and CO(17-16)/[CII] ratios as a function of $L_X$, which can be relevant for ALMA follow up of galaxies and quasars previously detected in [CII].

The Sunyaev-Zel'dovich (SZ) effect can potentially be used to investigate the heating of the circumgalactic medium and subsequent suppression of cold gas accretion onto the host galaxy caused by quasar feedback. We use a deep ALMA observation of HE0515-4414 in band 4, the most luminous quasar known at the peak of cosmic star formation (z=1.7), to search for the SZ signal tracing the heating of the galaxy's halo. ALMA's sensitivity to a broad range of spatial scales enables us to disentangle emitting compact sources from the negative, extended SZ signal. We obtain a marginal S-Z detection (~3.3$\sigma$) on scales of about 300 kpc (30-40 arcsec), at the 0.2 mJy level, 0.5 mJy after applying a correction factor for primary beam attenuation and flux that is resolved out by the array. We show that our result is consistent with a simulated ALMA observation of a similar quasar in the FABLE cosmological simulations. We emphasise that detecting an SZ signal is more easily achieved in the visibility plane than in the (inferred) images. We also confirm a marginal detection (3.2$\sigma$) of a potential SZ dip on smaller scales (<100 kpc) already claimed by other authors, possibly highlighting the complex structure of the halo heating. Finally, we use SZ maps from the FABLE cosmological simulations, convolved with ALMA simulations, to illustrate that band 3 observations are much more effective in detecting the SZ signal with higher significance, and discuss the optimal observing strategy.

We report the first detection obtained with ALMA of the [N II] 122$\mu$m line emission from a galaxy group BRI 1202-0725 at $z=4.69$ consisting of a QSO and a submilimeter-bright galaxy (SMG). Combining with a detection of [N II] 205$\mu$m line in both galaxies, we constrain the electron densities of the ionized gas based on the line ratio of [NII]122/205. The derived electron densities are $26^{+12}_{-11}$ and $134^{+50}_{-39}$ cm$^{-3}$ for the SMG and the QSO, respectively. The electron density of the SMG is similar to that of the Galactic Plane and to the average of the local spirals. Higher electron densities by up to a factor of three could, however, be possible for systematic uncertainties of the line flux estimates. The electron density of the QSO is comparable to high-$z$ star-forming galaxies at $z=1.5-2.3$, obtained using rest-frame optical lines and with the lower limits suggested from stacking analysis on lensed starbursts at $z=1-3.6$ using the same tracer of [NII]. Our results suggest a large scatter of electron densities in global scale at fixed star formation rates for extreme starbursts. The success of the [N II] 122$\mu$m and 205$\mu$m detections at $z=4.69$ demonstrates the power of future systematic surveys of extreme starbursts at $z>4$ for probing the ISM conditions and the effects on surrounding environments.

A tight relation between the [CII]158$\mu$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($\Sigma_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [CII]-deficiency we have developed an analytical model that fits local [CII] data, and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of $\Sigma_{\rm CII}$ with the surface star formation rate ($\Sigma_*$). However, for $\Sigma_* > 1 M_\odot~{\rm yr}^{-1}~{\rm kpc}^{-2}$, $\Sigma_{\rm CII}$ saturates. We conclude that underluminous [CII] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation ($\kappa_s \gg 1$), parameterized by the "burstiness" parameter $\kappa_s$; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [CII] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [OIII]88$\mu$m, CIII]1909A, CO lines). Simple formulae are given to interpret available data for low and high-$z$ galaxies.

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.

We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$\mu$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use "Althaea", a highly-resolved ($30\, \rm pc$) simulated prototypical Lyman Break Galaxy, in the redshift range $z = 6-7$, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: I) Merger , II) Spiral Disk, and III) Disturbed Disk. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in I), II), and III), respectively. We derive a generalised dynamical mass vs. [CII]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor $2$ ($4$). A Tully-Fisher relation is found for the observed high-$z$ galaxies, i.e. $L_{\rm[CII]}\propto (FWHM)^{1.80\pm 0.35}$ for which we provide a simple, physically-based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [CII]. When seen face-on, Althaea is always detected at $> 5\sigma$; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects.

We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([CII], [NII], and [OIII]). At $z=8$, the most massive galaxy, `Freesia', has an age $t_\star \simeq 409\,{\rm Myr}$, stellar mass $M_{\star} \simeq 4.2\times 10^9 {\rm M}_{\odot}$, and a star formation rate ${\rm SFR} \simeq 11.5\,{\rm M}_{\odot}{\rm yr}^{-1}$, due to a recent burst. Freesia has two stellar components (A and B) separated by $\simeq 2.5\, {\rm kpc}$; other 11 galaxies are found within $56.9 \pm 21.6 \, {\rm kpc}$. The mean ISRF in the Habing band is $G = 7.9\, G_0$ and is spatially uniform; in contrast, the ionisation parameter is $U = 2^{+20}_{-2} \times 10^{-3}$, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionising escape fraction from Freesia is $f_{\rm esc}\simeq 2\%$. While [CII] emission is extended (radius 1.54 kpc), [OIII] is concentrated in Freesia-A (0.85 kpc), where the ratio $\Sigma_{\rm [OIII]}/\Sigma_{\rm [CII]} \simeq 10$. As many high-$z$ galaxies, Freesia lies below the local [CII]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) which disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [CII]-SFR deviations.

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.

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.

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.

The James Webb Space Telescope (JWST) will enable the detection of optical emission lines in galaxies spanning a broad range of luminosities out to redshifts z>10. Measurements of key galaxy properties, such as star formation rate and metallicity, through these observations will provide unique insight into, e.g., the role of feedback from stars and active galactic nuclei (AGNs) in regulating galaxy evolution, the co-evolution of AGNs and host galaxies, the physical origin of the 'main sequence' of star-forming galaxies and the contribution by star-forming galaxies to cosmic reionization. We present an original framework to simulate and analyse observations performed with the Near Infrared Spectrograph (NIRSpec) on board JWST. We use the BEAGLE tool (BayEsian Analysis of GaLaxy sEds) to build a semi-empirical catalogue of galaxy spectra based on photometric spectral energy distributions (SEDs) of dropout galaxies in the Hubble Ultra Deep Field (HUDF). We demonstrate that the resulting catalogue of galaxy spectra satisfies different types of observational constraints on high redshift galaxies, and use it as input to simulate NIRSpec/prism (R~100) observations. We show that a single 'deep' (~100 ks) NIRSpec/prism pointing in the HUDF will enable S/N>3 detections of multiple optical emission lines in ~30 (~60) galaxies at z>6 (z~4-6) down to m_F160W<30 AB mag. Such observations will allow measurements of galaxy star formation rates, ionization parameters and gas-phase metallicities within factors of 1.5, mass-to-light ratios within a factor of 2, galaxy ages within a factor of 3 and V-band attenuation optical depths with a precision of 0.3.

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.

We present deep spectroscopic follow-up observations of the Bremer Deep Field (BDF) where the two $z\sim$7 bright Ly$\alpha$ emitters (LAE) BDF521 and BDF3299 were previously discovered by Vanzella et al. (2011) and where a factor of $\sim$3-4 overdensity of faint LBGs has been found by Castellano et al. (2016). We confirm a new bright Ly$\alpha$ emitter, BDF2195, at the same redshift of BDF521, $z=7.008$, and at only $\sim$90 kpc physical distance from it, confirming that the BDF area is likely an overdense, reionized region. A quantitative assessment of the Ly$\alpha$ fraction shows that the number of detected bright emitters is much higher than the average found at z$\sim$7, suggesting a high Ly$\alpha$ transmission through the inter-galactic medium (IGM). However, the line visibility from fainter galaxies is at odds with this finding, since no Ly$\alpha$ emission is found in any of the observed candidates with $M_{UV}>$-20.25. This discrepancy can be understood either if some mechanism prevents Ly$\alpha$ emission from fainter galaxies within the ionized bubbles from reaching the observer, or if faint galaxies are located outside the reionized area and bright LAEs are solely responsible for the creation of their own HII regions. A thorough assessment of the nature of the BDF region and of its sources of re-ionizing radiation will be made possible by JWST spectroscopic capabilities.

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 an original phenomenological model to describe the evolution of galaxy number counts, morphologies, and spectral energy distributions across a wide range of redshifts (0.2<z<15) and stellar masses [Log10 M/Msun >6]. Our model follows observed mass and luminosity functions of both star-forming and quiescent galaxies, and reproduces the redshift evolution of colors, sizes, star-formation and chemical properties of the observed galaxy population. Unlike other existing approaches, our model includes a self-consistent treatment of stellar and photoionized gas emission and dust attenuation based on the BEAGLE tool. The mock galaxy catalogs generated with our new model can be used to simulate and optimize extragalactic surveys with future facilities such as the James Webb Space Telescope (JWST), and to enable critical assessments of analysis procedures, interpretation tools, and measurement systematics for both photometric and spectroscopic data. As a first application of this work, we make predictions for the upcoming JWST Advanced Deep Extragalactic Survey (JADES), a joint program of the JWST/NIRCam and NIRSpec Guaranteed Time Observations teams. We show that JADES will detect, with NIRCam imaging, thousands of galaxies at z>6, and tens at z>10 at m_AB<30 (5-sigma) within the 236 arcmin^2 of the survey. The JADES data will enable accurate constraints on the evolution of the UV luminosity function at z>8, and resolve the current debate about the rate of evolution of galaxies at z>8. Ready to use mock catalogs and software to generate new realizations are publicly available as the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) package.

In this work we report deep MUSE observations of a Broad Absorption Line (BAL) quasar at z ~ 5, revealing a Lyman alpha nebula with a maximum projected linear size of ~ 60 kpc around the quasar (down to our 2-sigma SB limit per layer of ~ 9e-19 erg/s/cm^2/arcsec^2 for a 1 arcsec^2 aperture). After correcting for the cosmological surface brightness dimming, we find that our nebula, at z ~ 5, has an intrinsically less extended Lyman alpha emission than nebulae at lower redshift. However, such a discrepancy is greatly reduced when referring to comoving distances, which take into account the cosmological growth of dark matter (DM) haloes, suggesting a positive correlation between the size of Lyman alpha nebulae and the sizes of DM haloes/structures around quasars. Differently from the typical nebulae around radio-quiet non-BAL quasars, in the inner regions (~ 10 kpc) of the circumgalactic medium (CGM) of our source, the velocity dispersion of the Lyman alpha emission is very high (FWHM > 1000 km/s), suggesting that in our case we may be probing outflowing material associated with the quasar.

We present the first spatially-resolved rest-frame UV study of the gravitationally lensed galaxy, the 'Cosmic Horseshoe' (J1148+1930) at z=2.38. Our gravitational lens model shows that the system is made up of four star-forming regions, each ~4-8 kpc^2 in size, from which we extract four spatially exclusive regional spectra. We study the interstellar and wind absorption lines, along with CIII] doublet emission lines, in each region to investigate any variation in emission/absorption line properties. The mapped CIII] emission shows distinct kinematical structure, with velocity offsets of ~+/-50 km/s between regions suggestive of a merging system, and a variation in equivalent width that indicates a change in ionisation parameter and/or metallicity between the regions. Absorption line velocities reveal a range of outflow strengths, with gas outflowing between -200<v(km/s)<-50 relative to the systemic velocity of that region. Interestingly, the strongest gas outflow appears to emanate from the most diffuse star-forming region. The star-formation rates remain relatively constant (~8-16 M_sol/yr), mostly due to large uncertainties in reddening estimates. As such, the outflows appear to be 'global' rather than 'locally' sourced. We measure electron densities with a range of log(Ne)=3.92-4.36 cm^-3, and point out that such high densities may be common when measured using the CIII] doublet due to its large critical density. Overall, our observations demonstrate that while it is possible to trace variations in large scale gas kinematics, detecting inhomogeneities in physical gas properties and their effects on the outflowing gas may be more difficult. This study provides important lessons for the spatially-resolved rest-frame UV studies expected with future observatories, such as JWST.

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

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang.

AGN feedback, acting through strong outflows accelerated in the nuclear region of AGN hosts, is invoked as a key ingredient for galaxy evolution by many models to explain the observed BH-galaxy scaling relations. Recently, some direct observational evidence of radiative mode feedback in action has been finally found in quasars at $z$>1.5. However, it is not possible to study outflows in quasars at those redshifts on small scales ($\lesssim$100 pc), as spatial information is limited by angular resolution. This is instead feasible in nearby active galaxies, which are ideal laboratories to explore outflow structure and properties, as well as the effects of AGN on their host galaxies. In this proceeding we present preliminary results from the MAGNUM survey, which comprises nearby Seyfert galaxies observed with the integral field spectrograph VLT/MUSE. We focus on two sources, NGC 1365 and NGC 4945, that exhibit double conical outflows extending on distances >1 kpc. We disentangle the dominant contributions to ionization of the various gas components observed in the central $\sim$5.3 kpc of NGC 1365. An attempt to infer outflow 3D structure in NGC 4945 is made via simple kinematic modeling, suggesting a hollow cone geometry.

We investigate the morphology of the [CII] emission in a sample of "normal" star-forming galaxies at $5<z<7.2$ in relation to their UV (rest-frame) counterpart. We use new ALMA observations of galaxies at $z\sim6-7$ as well as a careful re-analysis of archival ALMA data. In total 29 galaxies were analysed, 21 of which are detected in [CII]. For several of the latter the [CII] emission breaks into multiple components. Only a fraction of these [CII] components, if any, is associated with the primary UV systems, while the bulk of the [CII] emission is associated either with fainter UV components, or not associated with any UV counterpart at the current limits. By taking into account the presence of all these components, we find that the L$_{\rm [CII]}$-SFR relation at early epochs is fully consistent with the local relation, but it has a dispersion of 0.48$\pm$0.07 dex, which is about two times larger than observed locally. We also find that the deviation from the local L$_{\rm [CII]}$-SFR relation has a weak anti-correlation with the EW(Ly$\alpha$). The morphological analysis also reveals that [CII] emission is generally much more extended than the UV emission. As a consequence, these primordial galaxies are characterised by a [CII] surface brightness generally much lower than expected from the local ${\rm \Sigma_{[CII}-\Sigma_{SFR}}$ relation. These properties are likely a consequence of a combination of different effects, namely: gas metallicity, [CII] emission from obscured star-forming regions, strong variations of the ionisation parameter, and circumgalactic gas in accretion or ejected by these primeval galaxies.

Himiko is one of the most luminous Ly{\alpha} emitters at z = 6.595. It has three star forming clumps detected in the rest-frame UV, with a total SFR = 20 M$_\odot$/yr. We report the ALMA detection of the [CII]158$\mu$m line emission in this galaxy with a significance of 9$\sigma$. The total [CII] luminosity (L[CII]= (1.2$\pm$0.2)$\times$10$^{8}$ L$_{\odot}$) is fully consistent with the local L[CII]-SFR relation. The ALMA high-angular resolution reveals that the [CII] emission is made of two distinct components. The brightest [CII] clump is extended over 4 kpc and is located on the peak of the Ly{\alpha} nebula, which is spatially offset by 1 kpc relative to the brightest UV clump. The second [CII] component is spatially unresolved (size $<$2 kpc) and coincident with one of the three UV clumps. While the latter component is consitent with the local L[CII]-SFR relation, the other components are scattered above and below the local relation. We shortly discuss the possible origin of the [CII] components and their relation with the star forming clumps traced by the UV emission.

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 new ALMA observations of the [OIII]88$\mu$m line and high angular resolution observations of the [CII]158$\mu$m line in a normal star forming galaxy at z$=$7.1. Previous [CII] observations of this galaxy had detected [CII] emission consistent with the Ly$\alpha$ redshift but spatially slightly offset relative to the optical (UV-rest frame) emission. The new [CII] observations reveal that the [CII] emission is partly clumpy and partly diffuse on scales larger than about 1kpc. [OIII] emission is also detected at high significance, offset relative to the optical counterpart in the same direction as the [CII] clumps, but mostly not overlapping with the bulk of the [CII] emission. The offset between different emission components (optical/UV and different far-IR tracers) is similar to what observed in much more powerful starbursts at high redshift. We show that the [OIII] emitting clump cannot be explained in terms of diffuse gas excited by the UV radiation emitted by the optical galaxy, but it requires excitation by in-situ (slightly dust obscured) star formation, at a rate of about 7 M$_{\odot}$/yr. Within 20 kpc from the optical galaxy the ALMA data reveal two additional [OIII] emitting systems, which must be star forming companions. We discuss that the complex properties revealed by ALMA in the z$\sim$7.1 galaxy are consistent with expectations by recent models and cosmological simulations, in which differential dust extinction, differential excitation and different metal enrichment levels, associated with different subsystems assembling a galaxy, are responsible for the different appearance of the system when observed with different tracers.

Recent observations have revealed massive galactic molecular outflows that may have physical conditions (high gas densities) required to form stars. Indeed, several recent models predict that such massive galactic outflows may ignite star formation within the outflow itself. This star-formation mode, in which stars form with high radial velocities, could contribute to the morphological evolution of galaxies, to the evolution in size and velocity dispersion of the spheroidal component of galaxies, and would contribute to the population of high-velocity stars, which could even escape the galaxy. Such star formation could provide in-situ chemical enrichment of the circumgalactic and intergalactic medium (through supernova explosions of young stars on large orbits), and some models also predict that it may contribute substantially to the global star formation rate observed in distant galaxies. Although there exists observational evidence for star formation triggered by outflows or jets into their host galaxy, as a consequence of gas compression, evidence for star formation occurring within galactic outflows is still missing. Here we report new spectroscopic observations that unambiguously reveal star formation occurring in a galactic outflow at a redshift of 0.0448. The inferred star formation rate in the outflow is larger than 15 Msun/yr. Star formation may also be occurring in other galactic outflows, but may have been missed by previous observations owing to the lack of adequate diagnostics.

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 report ALMA observations of the most massive (star forming) galaxy in the redshift range 3<z<4 within the whole GOODS-S field. We detect a large elongated structure of molecular gas around the massive primeval galaxy, traced by the CO(4-3) emission, and extended over 40 kpc. We infer a mass of the large gaseous structure of Mgas~2-6x10^11 Msun. About 60% of this mass is not directly associated with either the central galaxy or its two lower mass satellites. The CO extended structure is also detected in continuum thermal emission. The kinematics of the molecular gas shows the presence of different components, which cannot be ascribed to simple rotation. Furthermore, on even larger scales, we detect nine additional CO systems within a radius of 250 kpc from the massive galaxy and mostly distributed in the same direction as the CO elongated structure found in the central 40 kpc. The stacked images of these CO systems show detections in the thermal continuum and in the X-rays, suggesting that these systems are forming stars at a rate of 30-120 Msun/yr. We suggest that the extended gas structure, combined with its kinematic properties, and the gas rich star forming systems detected on larger scales, are tracing the inner and densest regions of large scale accreting streams, feeding the central massive galaxy. These results corroborate models of galaxy formation, in which accreting streams are clumpy and undergo some star formation (hence enriching the streams with metals) even before accreting onto the central galaxy.

We present Adaptive Optics assisted near-IR integral field spectroscopic observations of a luminous quasar at $z = 2.4$, previously observed as the first known example at high redshift of large scale quasar-driven outflow quenching star formation in its host galaxy. The nuclear spectrum shows broad and blueshifted H$\beta$ in absorption, which is tracing outflowing gas with high densities ($>10^8$ - $10^9$ cm$^{-3}$) and velocities in excess of 10,000 km s$^{-1}$. The properties of the outflowing clouds (covering factor, density, column density and inferred location) indicate that they likely originate from the Broad Line Region. The energetics of such nuclear regions are consistent with that observed in the large scale outflow, supporting models in which quasar driven outflows originate from the nuclear region and are energy conserving. We note that the asymmetric profile of both the H$\beta$ and H$\alpha$ emission lines is likely due to absorption by the dense outflowing gas along the line of sight. This outflow-induced asymmetry has implications on the estimation of the black hole mass using virial estimators, and warns about such effects for several other quasars characterized by similar line asymmetries. More generally, our findings may suggest a broader revision of the decomposition and interpretation of quasar spectral features, in order to take into account the presence of potential broad blueshifted Balmer absorption lines. Our high spatial resolution data also reveals redshifted, dynamically colder nebular emission lines, likely tracing an inflowing stream.

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 results on [CII]158$\mu$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$\alpha$ emission line in their spectra. The Ly$\alpha$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$\alpha$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$\alpha$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$\alpha$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

AGN driven outflows are invoked in numerical simulations to reproduce several observed properties of local galaxies. The z > 1 epoch is of particular interest as it was during this time that the volume averaged star formation and the accretion rate of black holes were maximum. Radiatively driven outflows are therefore believed to be common during this epoch. We aim to trace and characterize outflows in AGN hosts with high mass accretion rates at z > 1 using integral field spectroscopy. We obtain spatially-resolved kinematics of the [OIII]5007 line in two targets which reveal the morphology and spatial extension of the outflows. We present J and H+K band SINFONI observations of 5 AGNs at 1.2 < z < 2.2. To maximize the chance of observing radiatively driven outflows, our sample was pre-selected based on peculiar values of the Eddington ratio and the hydrogen column density of the surrounding interstellar medium. We observe high velocity (~600-1900 km/s) and kiloparsec scale extended ionized outflows in at least 3 of our targets, using [OIII]5007 line kinematics tracing the AGN narrow line region. We estimate the total mass of the outflow, the mass outflow rate, and the kinetic power of the outflows based on theoretical models and report on the uncertainties associated with them. We find mass outflow rates of ~1-10 M_sun/yr for the sample presented in this paper. Based on the high star formation rates of the host galaxies, the observed outflow kinetic power and the expected power due to the AGN, we infer that both star formation and AGN radiation could be the dominant source for the outflows. The outflow models suffer from large uncertainties, hence we call for further detailed observations for an accurate determination of the outflow properties to confirm the exact source of these outflows.

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.

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

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$\mu$Jy, to derive differential number counts down to 60$\mu$Jy and 100$\mu$Jy at $\lambda=$1.3 mm and $\lambda=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$\mu$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{\lambda = 1.1~mm} = 130~\mu$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{\lambda = 1.3~mm} = 60~\mu$Jy. At the faintest flux limits, i.e. 30$\mu$Jy and 40$\mu$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

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.

We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5-15 Msun/yr these systems are much more representative of the high-z galaxy population than other systems targeted in the past by millimeter observations. For the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected in any of the three galaxies, with 3sigma upper limits significantly lower than the [CII] emission observed in lower reshift galaxies. These results suggest that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. These findings are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.

We performed a kinematical analysis of the [CII] line emission of the BR 1202-0725 system at z~4,7 using ALMA observations. The most prominent sources of this system are a quasar and a submillimeter galaxy, separated by a projected distance of about 24 kpc and characterized by very high SFR, higher than 1000 Msun/yr. However, the ALMA observations reveal that these galaxies apparently have undisturbed rotating disks, which is at variance with the commonly accepted scenario in which strong star formation activity is induced by a major merger. We also detected faint components which, after spectral deblending, were spatially resolved from the main QSO and SMG emissions. The relative velocities and positions of these components are compatible with orbital motions within the gravitational potentials generated by the QSO host galaxy and the SMG, suggesting that they are smaller galaxies in interaction or gas clouds in accretion flows of tidal streams. We did not find any clear spectral evidence for outflows caused by AGN or stellar feedback. This suggests that the high star formation rates might be induced by interactions or minor mergers with these companions, which do not affect the large-scale kinematics of the disks, however. Our kinematical analysis also indicates that the QSO and the SMG have similar Mdyn, mostly in the form of molecular gas, and that the QSO host galaxy and the SMG are seen close to face-on with slightly different disk inclinations: the QSO host galaxy is seen almost face-on (i~15), while the SMG is seen at higher inclinations (i~25). Finally, the ratio between the black hole mass of the QSO, obtained from XShooter spectroscopy, and the Mdyn of the host galaxy is similar to value found in very massive local galaxies, suggesting that the evolution of black hole galaxy relations is probably better studied with dynamical than with stellar host galaxy masses.