Search

The energy released by AGN has the potential to heat or remove the gas of the ISM, thus likely impacting the cold molecular gas reservoir of host galaxies at first, with star formation following on longer timescales. Previous works on high-z galaxies have yielded conflicting results, possibly due to selection biases and other systematics. To provide a reliable benchmark for galaxy evolution models at cosmic noon (z=1-3), two surveys were conceived: SUPER and KASHz, both targeting unbiased X-ray-selected AGN at z>1 that span a wide bolometric luminosity range. In this paper, we assess the effects of AGN feedback on the molecular gas content of host galaxies in a statistically robust, uniformly selected, coherently analyzed sample of AGN at z=1-2.6, drawn from the KASHz and SUPER surveys. By using ALMA data in combination with dedicated SED modeling, we retrieve CO and FIR luminosity as well as $M_*$ of SUPER and KASHz AGN. We selected non-active galaxies from PHIBBS, ASPECS and multiple ALMA/NOEMA surveys of sub-mm galaxies. By matching the samples in z, $M_*$ and $L_{FIR}$, we compared the properties of AGN and non-active galaxies within a Bayesian framework. We find that AGN hosts at given $L_{FIR}$ are on average CO depleted compared to non-active galaxies, confirming what was previously found in the SUPER survey. Moreover, the molecular gas fraction distributions of AGN and non-active galaxies are statistically different, with that of of AGN being skewed to lower values. Our results indicate that AGN can indeed reduce the total cold molecular gas reservoir of their host galaxies. Lastly, by comparing our results with predictions from three cosmological simulations (TNG, Eagle and Simba) filtered to match the observed properties, we confirm already known discrepancies and highlight new ones between observations and simulations.[Abridged], Comment: 26 pages (17 in the main body of the paper and 9 in the appendix), 14 figures (11 in the main body of the paper and 3 in the appendix), 8 tables (2 in the main body of the paper and 6 in the appendix). Accepted for publication in A&A, in press

We present spatially resolved VLT/SINFONI spectroscopy with adaptive optics of type-2 active galactic nuclei (AGN) from the SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback (SUPER), which targeted X-ray bright ($L_{2-10 keV}\gtrsim10^{42}$ erg s$^{-1}$) AGN at Cosmic Noon ($z\sim2$). Our analysis of the rest-frame optical spectra unveils ionised outflows in all seven examined targets, as traced via [OIII]$\lambda$5007 line emission, moving at $v\gtrsim600$ km s$^{-1}$. In six objects these outflows are clearly spatially resolved and extend on 2-4 kpc scales, whereas marginally resolved in the remaining one. Interestingly, these SUPER type-2 AGN are all heavily obscured sources ($N_{H}\gtrsim10^{23}$ cm$^{-2}$) and host faster ionised outflows than their type-1 counterparts within the same range of bolometric luminosity ($L_{bol} \sim 10^{44.8-46.5}$ erg s$^{-1}$). SUPER has hence provided observational evidence that the type-1/type-2 dichotomy at $z\sim2$ might not be driven simply by projection effects, but might reflect two distinct obscuring life stages of active galaxies, as predicted by evolutionary models. Within this picture, SUPER type-2 AGN might be undergoing the 'blow-out' phase, where the large amount of obscuring material efficiently accelerates large-scale outflows via radiation pressure on dust, eventually unveiling the central active nucleus and signal the start of the bright, unobscured type-1 AGN phase. Moreover, the overall population of ionised outflows detected in SUPER has velocities comparable with the escape speed of their dark matter halos, and in general high enough to reach 30-50 kpc distances from the centre. These outflows are hence likely to sweep away the gas (at least) out of the baryonic disk and/or to heat the host gas reservoir, thus reducing and possibly quenching star formation., Comment: 21 pages, 13 figures. Key figure is 8. Accepted for publication in A&A

PDS 456 is the most luminous RQQ at z<0.3 and can be regarded as a local counterpart of the powerful QSOs shining at Cosmic Noon. It hosts a strong nuclear X-ray ultra-fast outflow, and a massive and clumpy CO(3-2) molecular outflow extending up to 5 kpc from the nucleus. We analyzed the first MUSE WFM and AO-NFM optical integral field spectroscopic observations of PDS456. The AO-NFM observations provide an unprecedented spatial resolution, reaching up to 280 pc. Our findings reveal a complex circumgalactic medium around PDS 456, extending up to a maximum projected size of ~46 kpc. This includes a reservoir of gas with a mass of ~1e7-1e8 Modot, along with eight companion galaxies, and a multi-phase outflow. WFM and NFM MUSE data reveal an outflow on a large scale (~12 kpc from the quasar) in [OIII], and on smaller scales (within 3 kpc) with higher resolution (about 280 pc) in Halpha, respectively. The [OIII] outflow mass rate is 2.3 +/- 0.2 Modot/yr which is significantly lower than those typically found in other luminous quasars. Remarkably, the Ha outflow shows a similar scale, morphology, and kinematics to the CO(3-2) molecular outflow, with the latter dominating in terms of kinetic energy and mass outflow rate by two and one orders of magnitude, respectively. Our results therefore indicate that mergers, powerful AGN activity, and feedback through AGN-driven winds will collectively contribute to shaping the host galaxy evolution of PDS 456, and likely, that of similar objects at the brightest end of the AGN luminosity function across all redshifts. Moreover, the finding that the momentum boost of the total outflow deviates from the expected energy-conserving expansion for large-scale outflows highlights the need of novel AGN-driven outflow models to comprehensively interpret these phenomena., Comment: 15 pages, 13 figures accepted for publication in Astronomy & Astrophysics

Observations at UV and optical wavelengths have revealed that galaxies at z~1-4 host star-forming regions, dubbed "clumps", which are believed to form due to the fragmentation of gravitationally unstable, gas-rich disks. However, the detection of the parent molecular clouds that give birth to such clumps is still possible only in a minority of galaxies, mostly at z~1. We investigated the [CII] and dust morphology of a z~3.4 lensed galaxy hosting four clumps detected in the UV continuum. We aimed to observe the [CII] emission of individual clumps that, unlike the UV, is not affected by dust extinction, to probe their nature and cold gas content. We conducted ALMA observations probing scales down to ~300 pc and detected three [CII] clumps. One (dubbed "NE") coincides with the brightest UV clump, while the other two ("SW" and "C") are not detected in the UV continuum. We do not detect the dust continuum. We converted the [CII] luminosity of individual clumps into molecular gas mass and found Mmol~10^8 Msun. By complementing it with the star formation rate (SFR) estimate from the UV continuum, we estimated the gas depletion time (tdep) of clumps and investigated their location in the Schmidt-Kennicutt plane. While the NE clump has a short tdep=0.16 Gyr, comparable with high-redshift starbursts, the SW and C clumps instead have longer tdep>0.65 Gyr and are likely probing the initial phases of star formation. The lack of dust continuum detection is consistent with the blue UV continuum slope estimated for this galaxy (beta~-2.5) and it indicates that dust inhomogeneities do not significantly affect the detection of UV clumps in this target. We pushed the observation of the cold gas content of individual clumps up to z~3.4 and showed that the [C II] line emission is a promising tracer of molecular clouds at high redshift, allowing the detection of clumps with a large range of depletion times., Comment: Accepted for publication in A&A

We used low- to high-frequency ALMA observations to investigate the cold gas and dust in ten QSOs at $z\gtrsim 6$. Our analysis of the CO(6-5) and CO(7-6) emission lines in the selected QSOs provided insights into their molecular gas masses, which average around $10^{10}\ \rm M_\odot$, consistent with typical values for high-redshift QSOs. Proprietary and archival ALMA observations in bands 8 and 9 enabled precise constraints on the dust properties and star formation rate (SFR) of four QSOs in our sample for the first time. The examination of the redshift distribution of dust temperatures revealed a general trend of increasing $T_{\rm dust}$ with redshift, which agrees with theoretical expectations. We computed a mean cold dust spectral energy distribution considering all ten QSOs. This offers a comprehensive view of the dust properties of high-$z$ QSOs. The QSOs marked by a more intense growth of the supermassive black hole (HYPERION QSOs) showed lower dust masses and higher gas-to-dust ratios on average, but their $\rm H_2$ gas reservoirs are consistent with those of other QSOs at the same redshift. The observed high SFR in our sample yields high SF efficiencies and thus very short gas depletion timescales ($\tau_{\rm dep}\sim 10^{-2}$ Gyr). Beyond supporting the paradigm that high-$z$ QSOs reside in highly star-forming galaxies, our findings portrayed an interesting evolutionary path at $z>6$. Our study suggests that they are undergoing rapid galaxy growth that might be regulated by strong outflows. Their inferred evolutionary path shows a convergence toward the massive end of the local relation, which supports the idea that they are candidate progenitors of local massive galaxies. The observed pathway involves intense BH growth followed by substantial galaxy growth, in contrast with a symbiotic growth scenario. The abstract has been shortened (full version in the article)., Comment: 29 pages; 6 tables; 18 figures. Accepted by A&A. A section about SF efficiency has been added compared to the previous version

We present a study of molecular gas, traced via CO (3-2) from ALMA data, of four z< 0.2, `radio quiet', type 2 quasars (log [L(bol)/(erg/s)] = 45.3 - 46.2; log [L(1.4 GHz)/(W/Hz)] = 23.7 - 24.3). Targets were selected to have extended radio lobes (>= 10 kpc), and compact, moderate-power jets (1 - 10 kpc; log [Pjet/(erg/s)]= 43.2 - 43.7). All targets show evidence of central molecular outflows, or injected turbulence, within the gas disks (traced via high-velocity wing components in CO emission-line profiles). The inferred velocities (Vout = 250 - 440 km/s) and spatial scales (0.6 - 1.6 kpc), are consistent with those of other samples of luminous low-redshift AGN. In two targets, we observe extended molecular gas structures beyond the central disks, containing 9 - 53 % of the total molecular gas mass. These structures tend to be elongated, extending from the core, and wrap-around (or along) the radio lobes. Their properties are similar to the molecular gas filaments observed around radio lobes of, mostly `radio loud', Brightest Cluster Galaxies. They have: projected distances of 5 - 13 kpc; bulk velocities of 100 - 340 km/s; velocity dispersion of 30 - 130 km/s; inferred mass outflow rates of 4 - 20 Msolar/yr; and estimated kinetic powers of log [Ekin/(erg/s)]= 40.3 - 41.7. Our observations are consistent with simulations that suggest moderate-power jets can have a direct (but modest) impact on molecular gas on small scales, through direct jet-cloud interactions. Then, on larger scales, jet-cocoons can push gas aside. Both processes could contribute to the long-term regulation of star formation., Comment: Accepted for publication in MNRAS; references updated; typos corrected

We present a comprehensive study of the molecular gas properties of 17 Type 2 quasars at $z <$ 0.2 from the Quasar Feedback Survey (L$_{[OIII]}$ > $10^{42.1}$ $\rm ergs^{-1}$), selected by their high [OIII] luminosities and displaying a large diversity of radio jet properties, but dominated by LIRG-like galaxies. With these data, we are able to investigate the impact of AGN and AGN feedback mechanisms on the global molecular interstellar medium. Using APEX and ALMA ACA observations, we measure the total molecular gas content using the CO(1-0) emission and homogeneously sample the CO spectral line energy distributions (SLEDs), observing CO transitions (J$_{up}$ = 1, 2, 3, 6, 7). We observe high $r_{21}$ ratios (r$_{21}$ = L'$_{CO(2-1)}$/L'$_{CO(1-0)}$) with a median $r_{21}$ = 1.06, similar to local (U)LIRGs (with $r_{21}$ $\sim$ 1) and higher than normal star-forming galaxies (with r$_{21}$ $\sim$ 0.65). Despite the high $r_{21}$ values, for the 7 targets with the required data we find low excitation in CO(6-5) & CO(7-6) ($r_{61}$ and $r_{62}$ < 0.6 in all but one target), unlike high redshift quasars in the literature, which are far more luminous and show higher line ratios. The ionised gas traced by [OIII] exhibit systematically higher velocities than the molecular gas traced by CO. We conclude that any effects of quasar feedback (e.g. via outflows and radio jets) do not have a significant instantaneous impact on the global molecular gas content and excitation and we suggest that it only occurs on more localised scales., Comment: 32 pages (20 in the main body of the paper and 12 in the appendix), 28 figures (10 in main body of paper and 18 in appendix) Accepted for publication in MNRAS. Data available at https://doi.org/10.25405/data.ncl.24312502

The [OIII] 5007 A emission line is the most common tracer of warm, ionized outflows in active galactic nuclei across cosmic time. JWST newly allows us to use mid-infrared spectral features at both high spatial and spectral resolution to probe these same winds. Here we present a comparison of ground-based, seeing-limited [OIII] and space-based, diffraction-limited [SIV] 10.51 micron maps of the powerful, kiloparsec-scale outflow in the Type 1 red quasar SDSS J110648.32+480712.3. The JWST data are from the Mid-InfraRed Instrument (MIRI). There is a close match in resolution between the datasets (0."6), in ionization potential of the O$^{+2}$ and S$^{+3}$ ions (35 eV), and in line sensitivity (1e-17 to 2e-17 erg/s/cm$^2$/arcsec$^2$). The [OIII] and [SIV] line shapes match in velocity and linewidth over much of the 20 kpc outflowing nebula, and [SIV] is the brightest line in the rest-frame 3.5-19.5 micron range, demonstrating its usefulness as a mid-IR probe of quasar outflows. [OIII] is nevertheless intriniscally brighter and provides better contrast with the point-source continuum, which is strong in the mid-IR. There is a strong anticorrelation of [OIII]/[SIV] with average velocity, which is consistent with a scenario of differential obscuration between the approaching (blueshifted) and receding (redshifted) sides of the flow. The dust in the wind may also obscure the central quasar, consistent with models that attribute red quasar extinction to dusty winds., Comment: resubmitted to match published version; minor changes

We present the hot molecular and warm ionised gas kinematics for 33 nearby ($0.001\lesssim z\lesssim0.056$) X-ray selected active galaxies using the H$_2 2.1218 \mu$m and Br$\gamma$ emission lines observed in the K-band with the Gemini Near-Infrared Field Spectrograph (NIFS). The observations cover the inner 0.04$-$2 kpc of each AGN at spatial resolutions of 4$-$250 pc with a velocity resolution of $\sigma_{\rm inst}\approx$20 ${\rm km s^{-1}}$. We find that 31 objects (94 per cent) present a kinematically disturbed region (KDR) seen in ionised gas, while such regions are observed in hot molecular gas for 25 galaxies (76 per cent). We interpret the KDR as being due to outflows with masses of 10$^2-$10$^7$ M$_\odot$ and 10$^0-$10$^4$ M$_\odot$ for the ionised and hot molecular gas, respectively. The ranges of mass-outflow rates ($\dot{M}_{\rm out}$) and kinetic power ($\dot{E}_{\rm K}$) of the outflows are 10$^{-3}-$10$^{1}$ M$_\odot$yr$^{-1}$ and $\sim$10$^{37}$$-$10$^{43}$ erg s$^{-1}$ for the ionised gas outflows, and 10$^{-5}$$-$10$^{-2}$ M$_\odot$ yr$^{-1}$ and 10$^{35}$$-$10$^{39}$ erg s$^{-1}$ for the hot molecular gas outflows. The median coupling efficiency in our sample is $\dot{E}_{K}/L_{\rm bol}\approx1.8\times10^{-3}$ and the estimated momentum fluxes of the outflows suggest they are produced by radiation-pressure in low-density environment, with possible contribution from shocks., Comment: 37 pages, published in MNRAS - A few typos in the text and in the label of Fg 1 were corrected in this version

We present spatially resolved H$\alpha$ properties of 21 type 1 AGN host galaxies at z$\sim$2 derived from the SUPER survey. These targets were observed with the adaptive optics capabilities of the SINFONI spectrograph, a near-infrared integral field spectrograph, that provided a median spatial resolution of 0.3 arcsec ($\sim$2 kpc). We model the H$\alpha$ emission line profile in each pixel to investigate whether it traces gas in the narrow line region or if it is associated with star formation. To do this, we first investigate the presence of resolved H$\alpha$ emission by removing the contribution of the AGN PSF. We find extended H$\alpha$ emission in sixteen out of the 21 type 1 AGN host galaxies (76%). Based on the BPT diagnostics, optical line flux ratios and the line widths (FWHM), we show that the H$\alpha$ emission in five galaxies is ionised by the AGN (30%), in four galaxies by star formation (25%) and for the rest (45%), the ionisation source is unconstrained. Two galaxies show extended H$\alpha$ FWHM $>$600 km/s, which is interpreted as a part of an AGN-driven outflow. Morphological and kinematic maps of H$\alpha$ emission in targets with sufficient signal-to-noise ratio suggest the presence of rotationally supported disks in six galaxies and possible presence of companions in four galaxies. In two galaxies, we find an anti-correlation between the locations of extended H$\alpha$ emission and [OIII]-based ionised outflows, indicating possible negative feedback at play. However, in the majority of galaxies, we do not find evidence of outflows impacting H$\alpha$ based star formation., Comment: 16 pages, 11 figures, 2 tables, accepted for publication in MNRAS; Appendix = 2 tables and 3 figures

The phase transition between galaxies and quasars is often identified with the rare population of hyper-luminous, hot dust-obscured galaxies. Galaxy formation models predict these systems to grow via mergers, that can deliver large amounts of gas toward their centers, induce intense bursts of star formation and feed their supermassive black holes. Here we report the detection of 24 galaxies emitting Lyman-alpha emission on projected physical scales of about 400 kpc around the hyper-luminous hot dust-obscured galaxy W0410-0913, at redshift z = 3.631, using Very Large Telescope observations. While this indicates that W0410-0913 evolves in a very dense environment, we do not find clear signs of mergers that could sustain its growth. Data suggest that if mergers occurred, as models expect, these would involve less massive satellites, with only a moderate impact on the internal interstellar medium of W0410-0913, which is sustained by a rotationally-supported fast-rotating molecular disk, as Atacama Large Millimeter Array observations suggest., Comment: 26 pages, 12 figures, 1 table. Authors' version. Published in Nature Communications on 05 August 2022

We present results from a combined radio polarization and emission line study of five type 2 quasars at $z<0.2$ with the Karl G. Jansky Very Large Array (VLA) B-array at 5 GHz and Hubble Space Telescope (HST) [O III] observations. These five sources are known to exhibit close association between radio structures and ionized gas morphology and kinematics. Four sources (J0945+1737, J1000+1242, J1356+1026 and J1430+1339) show polarization in the current data. J1010+1413 is the unpolarized source in our sample. We detect $0.5-1\%$ fractional polarization in the radio cores and a high fractional polarization ($10-30\%$) in the lobes of these sources. The morphological, spectral and polarization properties suggest a jet origin for radio emission in J0945+1737, J1000+1242, J1010+1413 and J1430+1339 whereas the current data cannot fully discern the origin of radio emission (jet or wind) in J1356+1026. An anti-correlation between various polarized knots in the radio and [O III] emission is observed in our sources, similar to that observed in some radio-loud AGN in the literature. This suggests that the radio emission is likely to be depolarized by the emission-line gas. By modeling the depolarization effects, we estimate the size of the emission-line gas clouds to be $\sim(2.8\pm1.7)\times10^{-5}$ parsec and the amount of thermal material mixed with the synchrotron plasma to be $\sim(1.01\pm0.08)\times10^{6}$M$_\odot$ in the lobe of J0945+1737 (which exhibits the most prominent polarization signature in its lobe). The current work demonstrates that the interplay of jets/winds and emission-line gas is most likely responsible for the nature of radio outflows in radio-quiet AGN., Comment: 17 pages, 8 figures; v2: matches published version

We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas properties. The jets are seen to drive galaxy-wide bi-conical turbulent outflows, reaching W80 = 1000-1300 km/s, in the ionised phase (traced via optical emission-lines), which also have increased electron densities compared to the quiescent gas. The turbulent gas is driven perpendicular to the jet axis and is escaping along the galaxy minor axis, reaching 7.5 kpc on both sides. Traced via CO(3-2) emission, the turbulent material in molecular gas phase is one-third as spatially extended and has 3 times lower velocity-dispersion as compared to ionised gas. The jets are seen to be strongly interacting with the interstellar medium (ISM) through enhanced ionised emission and disturbed/depleted molecular gas at the jet termini. We see further evidence for jet-induced feedback through significantly higher stellar velocity-dispersion aligned, and co-spatial with, the jet axis (<5 deg). We discuss possible negative and positive feedback scenarios arising due to the interaction of the low-power jets with the ISM in the context of recent jet-ISM interaction simulations, which qualitatively agree with our observations. We discuss how jet-induced feedback could be an important feedback mechanism even in bolometrically luminous 'radio-quiet' quasars., Comment: Accepted for publication in MNRAS. Main manuscript has 21 pages with 8 figures. Supplementary material is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

We study the gas distribution and kinematics of the inner kpc of six moderately luminous (43.43$\leq$log$L_{\rm bol}\leq$44.83) nearby ($0.004 \leq z \leq 0.014$) Seyfert galaxies observed with the Near-infrared Integral Field Spectrograph (NIFS) in the J(1.25$\mu$m) and K(2.2$\mu$m) bands. We analyse the most intense emission lines detected on these spectral wavebands: [Fe II] 1.2570 $\mu$m and Pa$\beta$, which trace the ionised gas in the partially and fully ionised regions, and H$_2$2.1218$\mu$m, that traces the hot ($\sim 2000$ K) molecular gas. The dominant kinematic component is rotation in the disc of the galaxies, except for the ionised gas in NGC 5899 which shows only weak signatures of a disc component. We find ionised gas outflow in four galaxies, while signatures of H$_2$ outflows are seen in three galaxies. The ionised gas outflows display velocities of a few hundred km\,s$^{-1}$, and their mass-outflow rates are in the range $0.005 - 12.49$M$_{\odot}$yr$^{-1}$ . Their kinetic powers correspond to $0.005-0.7$ per cent of the AGN bolometric luminosities. Besides rotation and outflows signatures in some cases, the H 2 kinematics reveals also inflows in three galaxies. The inflow velocities are $50-80$km\,s$^{-1}$ and the mass inflow rates are in the range $1 - 9 \times 10^{-4}$M$_{\odot}$yr$^{-1}$ for hot molecular gas. These inflows might be only the hot skin of the total inflowing gas, which is expected to be dominated by colder gas. The mass inflow rates are lower than the current accretion rates to the AGN, and the ionised outflows are apparently disturbing the gas in the inner kpc., Comment: Accepted for publication in November 26, 2021

Giant star-forming regions (clumps) are widespread features of galaxies at $z \approx 1-4$. Theory predicts that they can play a crucial role in galaxy evolution if they survive to stellar feedback for > 50 Myr. Numerical simulations show that clumps' survival depends on the stellar feedback recipes that are adopted. Up to date, observational constraints on both clumps' outflows strength and gas removal timescale are still uncertain. In this context, we study a line-emitting galaxy at redshift $z \simeq 3.4$ lensed by the foreground galaxy cluster Abell 2895. Four compact clumps with sizes $\lesssim$ 280 pc and representative of the low-mass end of clumps' mass distribution (stellar masses $\lesssim 2\times10^8\ {\rm M}_\odot$) dominate the galaxy morphology. The clumps are likely forming stars in a starbursting mode and have a young stellar population ($\sim$ 10 Myr). The properties of the Lyman-$\alpha$ (Ly$\alpha$) emission and nebular far-ultraviolet absorption lines indicate the presence of ejected material with global outflowing velocities of $\sim$ 200-300 km/s. Assuming that the detected outflows are the consequence of star formation feedback, we infer an average mass loading factor ($\eta$) for the clumps of $\sim$ 1.8 - 2.4 consistent with results obtained from hydro-dynamical simulations of clumpy galaxies that assume relatively strong stellar feedback. Assuming no gas inflows (semi-closed box model), the estimates of $\eta$ suggest that the timescale over which the outflows expel the molecular gas reservoir ($\simeq 7\times 10^8\ \text{M}_\odot$) of the four detected low-mass clumps is $\lesssim$ 50 Myr., Comment: 19 pages, 11 figures + Appendix. Published in MNRAS

We present the discovery of copious molecular gas in the halo of cid_346, a $z=2.2$ quasar studied as part of the SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER). New Atacama Compact Array (ACA) CO(3-2) observations detect a much higher flux (by a factor of $14\pm5$) than measured on kiloparsec scales ($r\lesssim8$ kpc) using previous snapshot Atacama Large Millimeter/submillimeter Array data. Such additional CO(3-2) emission traces a structure that extends out to $r\sim200$ kpc in projected size, as inferred through direct imaging and confirmed by an analysis of the uv visibilities. This is the most extended molecular circumgalactic medium (CGM) reservoir that has ever been mapped. It shows complex kinematics, with an overall broad line profile (FWHM $= 1000$ km/s) that is skewed towards redshifted velocities up to at least $v\sim1000$ km/s. Using the optically thin assumption, we estimate a strict lower limit for the total molecular CGM mass observed by ACA of $M_{mol}^{CGM}>10^{10}~M_{\odot}$. There is however room for up to $M^{CGM}_{mol}\sim 1.7\times 10^{12}$ $M_{\odot}$, once optically thick CO emission with $\alpha_{\rm CO}=3.6$ $\rm M_{\odot}~(K~km~s^{-1}~pc^2)^{-1}$ and $L^{\prime}_{CO(3-2)}/L^{\prime}_{CO(1-0)}=0.5$ are assumed. Since cid_346 hosts quasar-driven ionised outflows and since there is no evidence of merging companions or an overdensity, we suggest that outflows may have played a crucial rule in seeding metal-enriched, dense gas on halo scales. However, the origin of such an extended molecular CGM remains unclear., Comment: Accepted for publication in A&A Letters. 9 pages and 9 figures including Appendix. Version updated after language editing

We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$\alpha$ emission and [OIII] outflows, with matched-resolution observations of the rest-frame FIR emission. We present high-resolution ALMA Band 7 observations of eight X-ray selected AGN at z~2 from the SUPER sample, targeting the rest-frame ~260 um continuum at ~2 kpc (0.2'') resolution. We detected 6 out of 8 targets with S/N>10 in the ALMA maps, with continuum flux densities F = 0.27-2.58 mJy and FIR half-light radii Re = 0.8-2.1 kpc. The FIR Re of our sample are comparable to other AGN and star-forming galaxies at a similar redshift from the literature. However, we find that the mean FIR size in X-ray AGN (Re = 1.16+/- 0.11 kpc) is slightly smaller than in non-AGN (Re = 1.69+/-0.13 kpc). From SED fitting, we find that the main contribution to the 260 um flux density is dust heated by star formation, with < 4% contribution from AGN-heated dust and < 1% from synchrotron emission. The majority of our sample show different morphologies for the FIR (mostly due to reprocessed stellar emission) and the ionised gas emission (H$\alpha$ and [OIII], mostly due to AGN emission). This could be due to the different locations of dust and ionised gas, the different sources of the emission (stars and AGN), or the effect of dust obscuration. We are unable to identify any residual H$\alpha$ emission, above that dominated by AGN, that could be attributed to star formation. Under the assumption that the FIR emission is a reliable tracer of obscured star formation, we find that the obscured star formation activity in these AGN host galaxies is not clearly affected by the ionised outflows. However, we cannot rule out that star formation suppression is happening on smaller spatial scales than the ones we probe with our observations (< 2 kpc) or on different timescales., Comment: 32 pages, 13 figures. Accepted for publication in A&A. The ALMA and SINFONI maps presented in this paper are available at: https://doi.org/10.25405/data.ncl.23694765

We present high-resolution ($\sim$2.4\,kpc) ALMA band 7 observations (rest-frame $\lambda \sim 250\mu$m) of three powerful z$\sim$2.5 quasars ($L_{\rm bol}=10^{47.3}$-$10^{47.5}$ ergs s$^{-1}$). These targets have previously been reported as showing evidence for suppressed star formation based on cavities in the narrow H$\alpha$ emission at the location of outflows traced with [O~{\sc iii}] emission. Here we combine the ALMA observations with a re-analysis of the VLT/SINFONI data to map the rest-frame far-infrared emission, H$\alpha$ emission, and [O~{\sc iii}] emission. In all targets we observe high velocity [O~{\sc iii}] gas (i.e., W80$\sim$1000--2000\,km\,s$^{-1}$) across the whole galaxy. We do not identify any H$\alpha$ emission that is free from contamination from AGN-related processes; however, based on SED analyses, we show that the ALMA data contains a significant dust-obscured star formation component in two out of the three systems. This dust emission is found to be extended over $\approx$1.5--5.5\,kpc in the nuclear regions, overlaps with the previously reported H$\alpha$ cavities and is co-spatial with the peak in surface brightness of the [O~{\sc iii}] outflows. In summary, within the resolution and sensitivity limits of the data, we do not see any evidence for a instantaneous shut down of in-situ star formation caused directly by the outflows. However, similar to the conclusions of previous studies and based on our measured star formation rates, we do not rule out that the global host galaxy star formation could be suppressed on longer timescales by the cumulative effect of quasar episodes during the growth of these massive black holes., Comment: Accepted for publication in MNRAS

We present the first results from the Quasar Feedback Survey, a sample of 42 z<0.2, [O III] luminous AGN (L[O III]>10^42.1 ergs/s) with moderate radio luminosities (i.e. L(1.4GHz)>10^23.4 W/Hz; median L(1.4GHz)=5.9x10^23 W/Hz). Using high spatial resolution (~0.3-1 arcsec), 1.5-6 GHz radio images from the Very Large Array, we find that 67 percent of the sample have spatially extended radio features, on ~1-60 kpc scales. The radio sizes and morphologies suggest that these may be lower radio luminosity versions of compact, radio-loud AGN. By combining the radio-to-infrared excess parameter, spectral index, radio morphology and brightness temperature, we find radio emission in at least 57 percent of the sample that is associated with AGN-related processes (e.g. jets, quasar-driven winds or coronal emission). This is despite only 9.5-21 percent being classified as radio-loud using traditional criteria. The origin of the radio emission in the remainder of the sample is unclear. We find that both the established anti-correlation between radio size and the width of the [O III] line, and the known trend for the most [O III] luminous AGN to be associated with spatially-extended radio emission, also hold for our sample of moderate radio luminosity quasars. These observations add to the growing evidence of a connection between the radio emission and ionised gas in quasar host galaxies. This work lays the foundation for deeper investigations into the drivers and impact of feedback in this unique sample., Comment: Accepted for publication in MNRAS. Data products from this paper and the survey pilot papers are available through our website: https://blogs.ncl.ac.uk/quasarfeedbacksurvey/. The extensive supplementary material (containing additional figures and information on individual targets) is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

Outflows driven by active galactic nuclei (AGN) are expected to have a significant impact on the host galaxy evolution, but it is still debated how they are accelerated and propagate on galaxy-wide scales. This work addresses these questions by studying the link between X-ray, nuclear ultra-fast outflows (UFOs) and extended ionised outflows, for the first time in two quasars close to the peak of AGN activity ($z\sim2$), where AGN feedback is expected to be more effective. As targets, we selected two multiple-lensed quasars at $z\sim1.5$, HS 0810+2554 and SDSS J1353+1138, known to host UFOs and observed with the near-IR integral field spectrometer SINFONI at the VLT. We performed a kinematical analysis of the [O III]$\lambda$5007 optical emission line, in order to trace the presence of ionised outflows. We detected spatially resolved ionised outflows in both galaxies, extended more than 8 kpc and moving up to $v>2000$ km/s. We derived mass outflow rates of $\sim$12 M$_{sun}$/yr and $\sim$2 M$_{sun}$/yr for HS 0810+2554 and SDSS J1353+1138. Comparing with the co-hosted UFO energetics, the ionised outflow energetics in HS 0810+2554 is broadly consistent with a momentum-driven regime of wind propagation, while in SDSS J1353+1138 it differs by a factor of $\sim$100 from theoretical predictions, requiring either a massive molecular outflow or a high variability of the AGN activity to account for such a discrepancy. By additionally considering our results with those from the small sample of well-studied objects (all local but one), with both UFO and extended (ionised/atomic/molecular) outflow detections, we found that in 10 out of 12 galaxies the large-scale outflow energetics is consistent with the theoretical predictions of either a momentum- or an energy-driven scenario. This suggests that such models explain relatively well the acceleration mechanism of AGN-driven winds on large scales., Comment: 19 pages, 9 figures. Accepted for publication in A&A

Feedback from AGN is thought to be key in shaping the life cycle of their host galaxies by regulating star-formation activity. Therefore, to understand the impact of AGN on star formation, it is essential to trace the molecular gas out of which stars form. In this paper we present the first systematic study of the CO properties of AGN hosts at z~2 for a sample of 27 X-ray selected AGN spanning two orders of magnitude in AGN bolometric luminosity (Lbol= 10^44.7-10^46.9 erg/s) by using ALMA Band 3 observations of the CO(3-2) transition (~1" angular resolution). To search for evidence of AGN feedback on the CO properties of the host galaxies, we compared our AGN with a sample of inactive (i.e., non-AGN) galaxies from the PHIBSS survey with similar redshift, stellar masses, and SFRs. We used the same CO transition as a consistent proxy for the gas mass for the two samples in order to avoid systematics involved when assuming conversion factors. By adopting a Bayesian approach to take upper limits into account, we analyzed CO luminosities as a function of stellar masses and SFRs, as well as the ratio LCO(3-2)/M* (proxy for the gas fraction). The two samples show statistically consistent trends in the LCO(3-2)-Lfir and LCO(3-2)-M* planes. However, there are indications that AGN feature lower CO(3-2) luminosities (0.4-0.7 dex) than inactive galaxies at the 2-3sigma level when we focus on the subset of parameters where the results are better constrained and on the distribution of the mean LCO(3-2)/M*. Therefore, even by conservatively assuming the same excitation factor r31, we would find lower molecular gas masses in AGN, and assuming higher r31 would exacerbate this difference. We interpret our result as a hint of the potential effect of AGN activity (e.g., radiation and outflows), which may be able to heat, excite, dissociate, and/or deplete the gas reservoir of the host galaxies. (abridged), Comment: Key figures: 2, 3, 4. 22 pages, 9 figures (5 in appendices), 3 tables. Accepted for publication in A&A

The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) was designed to conduct a blind search for AGN-driven outflows on X-ray selected AGN at redshift z$\sim$2 with high ($\sim$2 kpc) spatial resolution, and correlate them to the properties of the host galaxy and central black hole. The main aims of this paper are: a) to derive reliable estimates for the BH mass and accretion rates for the Type-1 AGN in this survey; b) to characterize the properties of the AGN driven winds in the BLR. We analyzed rest-frame optical and UV spectra of 21 Type-1 AGN. We found that the BH masses estimated from H$\alpha$ and H$\beta$ lines are in agreement. We estimate BH masses in the range Log(M$\rm_{BH}/M_{\odot}$)=8.4-10.8 and Eddington ratios $\rm\lambda_{Edd}$ =0.04-1.3. We confirm that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blue-shifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR, with velocities up to $\sim$4700 km/s. We confirm the strong dependence of the BLR wind velocity with the UV-to-Xray continuum slope, L$\rm_{Bol}$ and $\rm\lambda_{Edd}$. We inferred BLR mass outflow rates in the range 0.005-3 M$_{\odot}$/yr, showing a correlation with the bolometric luminosity consistent with that observed for ionized winds in the NLR and X-ray winds detected in local AGN, and kinetic power $\sim$10$^{[-7:-4]}\times$ L$\rm_{Bol}$. Finally, we found an anti-correlation between the equivalent width of the [OIII] line with respect to the CIV shift, and a positive correlation with [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGN at z$\sim$2, support a scenario where BLR winds are connected to galaxy scale detected outflows, and are capable of affecting the gas in the NLR located at kpc scale., Comment: Accepted for publication in A&A

The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) aims at tracing and characterizing ionized gas outflows and their impact on star formation in a statistical sample of X-ray selected Active Galactic Nuclei (AGN) at z$\sim$2. We present the first SINFONI results for a sample of 21 Type-1 AGN spanning a wide range in bolometric luminosity (log $\mathrm{L_{bol}}$ = 45.4-47.9 erg/s). The main aims of this paper are determining the extension of the ionized gas, characterizing the occurrence of AGN-driven outflows, and linking the properties of such outflows with those of the AGN. We use Adaptive Optics-assisted SINFONI observations to trace ionized gas in the extended narrow line region using the [OIII]5007 line. We classify a target as hosting an outflow if its non-parametric velocity of the [OIII] line, $\mathrm{w_{80}}$, is larger than 600 km/s. We study the presence of extended emission using dedicated point-spread function (PSF) observations, after modelling the PSF from the Balmer lines originating from the Broad Line Region. We detect outflows in all the Type-1 AGN sample based on the $\mathrm{w_{80}}$ value from the integrated spectrum, which is in the range 650-2700 km/s. There is a clear positive correlation between $\mathrm{w_{80}}$ and the AGN bolometric luminosity (99% correlation probability), but a weaker correlation with the black hole mass (80% correlation probability). A comparison of the PSF and the [OIII] radial profile shows that the [OIII] emission is spatially resolved for $\sim$35% of the Type-1 sample and the outflows show an extension up to $\sim$6 kpc. The relation between maximum velocity and the bolometric luminosity is consistent with model predictions for shocks from an AGN driven outflow. The escape fraction of the outflowing gas increase with the AGN luminosity, although for most galaxies, this fraction is less than 10%., Comment: Accepted for publication in A&A, 34 pages, 14 figures and 5 tables

We use a sample of powerful z~0.1 type 2 quasars ('obscured'; log[L(AGN)/erg/s]>~45), which host kiloparsec-scale ionized outflows and jets, to identify possible signatures of AGN feedback on the total molecular gas reservoirs of their host galaxies. Specifically, we present Atacama Pathfinder EXperiment (APEX) observations of the CO(2-1) transition for nine sources and the CO(6-5) for a subset of three. We find that the majority of our sample reside in starburst galaxies (average specific star formation rates of 1.7/Gyr), with the seven CO-detected quasars also having large molecular gas reservoirs (average Mgas = 1.3x10^10Msun), even though we had no pre-selection on the star formation or molecular gas properties. Despite the presence of quasars and outflows, we find that the molecular gas fractions (Mgas/Mstar = 0.1-1.2) and depletion times (Mgas/SFR = 0.16-0.95Gyr) are consistent with those expected for the overall galaxy population with matched stellar masses and specific star formation rates. Furthermore, for at least two of the three targets with the required measurements, the CO(6-5)/CO(2-1) emission-line ratios are consistent with star formation dominating the CO excitation over this range of transitions. The targets in our study represent a gas-rich phase of galaxy evolution with simultaneously high levels of star formation and nuclear activity; furthermore, the jets and outflows do not have an immediate appreciable impact on the global molecular gas reservoirs., Comment: Accepted for publication in MNRAS. Supplementary material is available for download under "Ancillary files"

We present X-ray spectra of 185 bright sources detected in the XMM-Newton deep survey of the Chandra Deep Field South, combining the three EPIC cameras. The 2-10 keV flux limit of the sample is 2e-15 erg/s/cm2. The sources are distributed over a redshift range of z=0.1-3.8. Eleven new X-ray redshift measurements are included. A spectral analysis was performed using a simple model to obtain absorbing column densities, rest-frame 2-10 keV luminosities and Fe K line properties of 180 sources at z>0.4. Obscured AGN are found to be more abundant toward higher redshifts. Using the XMM-Newton data alone, seven Compton-thick AGN candidates are identified, which makes the Compton-thick AGN fraction to be ~4%. An exploratory spectral inspection method with two rest-frame X-ray colours and an Fe line strength indicator is introduced and tested against the results from spectral fitting. This method works reasonably well to characterise a spectral shape and can be useful for a pre-selection of Compton-thick AGN candidates. We found six objects exhibiting broad Fe K lines out of 21 unobscured AGN of best data quality, implying a detection rate of ~30%. Five redshift spikes, each of which has more than six sources, are identified in the redshift distribution of the X-ray sources. Contrary to the overall trend, the sources at the two higher-redshift spikes at z=1.61 and z=2.57 have puzzlingly low obscuration., Comment: Accepted for publication in Astronomy & Astrophysics. Table 1 and Fig. 2 are available in electronic form the CDS. Typos corrected and some language editing and references added

As part of our KMOS AGN Survey at High-redshift (KASHz), we present spatially-resolved VLT/KMOS and VLT/SINFONI spectroscopic data and ALMA 870$\mu$m continuum imaging of eight $z$=1.4--2.6 moderate AGN ($L_{\rm 2-10 \rm kev}$ = $10^{42} - 10^{45}$ ergs s$^{-1}$). We map [OIII], H$\alpha$ and rest-frame FIR emission to search for any spatial anti-correlation between ionised outflows (traced by the [OIII] line) and star formation (SF; traced by H$\alpha$ and FIR), that has previously been claimed for some high-z AGN and used as evidence for negative and/or positive AGN feedback. Firstly, we conclude that H$\alpha$ is unreliable to map SF inside our AGN host galaxies based on: (i) SF rates inferred from attenuation-corrected H$\alpha$ can lie below those inferred from FIR; (ii) the FIR continuum is more compact than the H$\alpha$ emission by a factor of $\sim 2$ on average; (iii) in half of our sample, we observe significant spatial offsets between the FIR and H$\alpha$ emission, with an average offset of $1.4\pm0.6$ kpc. Secondly, for the five targets with outflows we find no evidence for a spatial anti-correlation between outflows and SF using either H$\alpha$ or FIR as a tracer. This holds for our re-analysis of a famous $z$=1.6 X-ray AGN (`XID 2028') where positive and negative feedback has been previously claimed. Based on our results, any impact on SF by ionised outflows must be subtle, either occurring on scales below our resolution, or on long timescales., Comment: Accepted for publication by MNRAS, 24 pages, 12 Figures

HATLAS J084933.4+021443 was identified as a dusty starburst via its rest-frame far-IR emission. Multi-frequency imaging and spectroscopy revealed a cluster of four dusty galaxies at $z = 2.41$, covering 80 kpc. Here, we use ALMA to confirm a more distant, fifth protocluster member, and present X-ray and rest-frame optical imaging spectroscopy of the brightest, an unlensed hyperluminous IR galaxy (HyLIRG). The data reveal broad H \alpha and bright [N II] lines, and bright X-ray emission, characteristics that betray a Type-1 active galactic nucleus (AGN), strengthening evidence that AGN are ubiquitous amongst HyLIRGs. The accreting black hole is super massive, $M_{\rm bh}\approx 2\times 10^9$ M$_\odot$, with little intrinsic absorption, $N_{\rm H}\approx 5\times 10^{21}$ cm$^{-2}$. The X-ray properties suggest the accretion luminosity rivals that of the starburst, yet it is not obvious where this emerges in its panchromatic spectral energy distribution (SED). We outline three scenarios that could give rise to the observed characteristics, and how we might distinguish between them. In the first, we see the AGN through the host galaxy because of the cavity it excavates. In the others, the AGN is not co-spatial with the starburst, having been ejected via asymmetric gravitational radiation, or having evolved towards the naked quasar phase in an unseen companion., Comment: 10 pages; 5 figures; in press at MNRAS

We present ALMA Band~7 observations at $850\mu$m of 20 luminous ($\log\, L_{\rm bol}>46.9$ [erg s$^{-1}$]) unobscured quasars at $z\sim2$. We detect continuum emission for 19/20 quasars. After subtracting an AGN contribution, we measure the total far-IR luminosity for 18 quasars, assuming a modified blackbody model, and attribute the emission as indicative of the star formation rate (SFR). Our sample can be characterized with a log-normal SFR distribution having a mean of 140 $M_\odot$ yr$^{-1}$ and a dispersion of 0.5 dex. Based on an inference of their stellar masses, the SFRs are similar, in both the mean and dispersion, with star-forming main-sequence galaxies at the equivalent epoch. Thus, there is no evidence for a systematic enhancement or suppression (i.e., regulation or quenching) of star formation in the hosts of the most luminous quasars at $z\sim2$. These results are consistent with the Magneticum cosmological simulation, while in disagreement with a widely recognized phenomenological model that predicts higher SFRs than observed here based on the high bolometric luminosities of this sample. Furthermore, there is only a weak relation between SFR and accretion rate onto their supermassive black holes both for average and individual measurements. We interpret these results as indicative of star formation and quasar accretion being fed from the available gas reservoir(s) in their host with a disconnect due to their different physical sizes, temporal scales, and means of gas processing., Comment: 20 pages, 12 figures, re-submitted to MNRAS

We select a sample of 90 obscured (type2) AGN with 1.45

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., Comment: Accepted in A&A

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

We present 1-7 GHz high-resolution radio imaging (VLA and e-MERLIN) and spatially-resolved ionized gas kinematics for ten z<0.2 type~2 `obscured' quasars (log [L(AGN)/(erg/s)]>~45) with moderate radio luminosities (log [L(1.4GHz)/(W/Hz)]=23.3-24.4). These targets were selected to have known ionized outflows based on broad [OIII] emission-line components (FWHM~800-1800 km/s). Although `radio-quiet' and not `radio AGN' by many traditional criteria, we show that for nine of the targets, star formation likely accounts for <~10 per cent of the radio emission. We find that ~80-90 per cent of these nine targets exhibit extended radio structures on 1-25 kpc scales. The quasars' radio morphologies, spectral indices and position on the radio size-luminosity relationship reveals that these sources are consistent with being low power compact radio galaxies. Therefore, we favour radio jets as dominating the radio emission in the majority of these quasars. The radio jets we observe are associated with morphologically and kinematically distinct features in the ionized gas, such as increased turbulence and outflowing bubbles, revealing jet-gas interaction on galactic scales. Importantly, such conclusions could not have been drawn from current low-resolution radio surveys such as FIRST. Our observations support a scenario where compact radio jets, with modest radio luminosities, are a crucial feedback mechanism for massive galaxies during a quasar phase., Comment: Accepted for publication in MNRAS. Main manuscript is 22 pages with 12 figures. Supplementary material is available for download under "Ancillary files" or by downloading the source file listed under "Other formats"

We present a multi-wavelength study of seven AGN at spectroscopic redshift >2.5 in the 7 Ms Chandra Deep Field South, selected to have good FIR/sub-mm detections. Our aim is to investigate the possibility that the obscuration observed in the X-rays can be produced by the interstellar medium (ISM) of the host galaxy. Based on the 7 Ms Chandra spectra, we measured obscuring column densities N$_{H, X}$ in excess of 7x10$^{22}$ cm$^{-2}$ and intrinsic X-ray luminosities L$_{X}$>10$^{44}$ erg s$^{-1}$ for our targets, as well as equivalent widths for the Fe K$\alpha$ emission line EW>0.5-1 keV. We built the UV-to-FIR spectral energy distributions by using broad-band photometry from CANDELS and Herschel catalogs. By means of an SED decomposition technique, we derived stellar masses (M$_{*}$~10$^{11}$ Msun), IR luminosities (L$_{IR}$>10$^{12}$ Lsun), star formation rates (SFR~190-1680 Msun yr$^{-1}$) and AGN bolometric luminosities (L$_{bol}$~10$^{46}$ erg s$^{-1}$) for our sample. We used an empirically-calibrated relation between gas masses and FIR/sub-mm luminosities and derived M$_{gas}$~0.8-5.4x10$^{10}$ Msun. High-resolution (0.3-0.7'') ALMA data (when available, CANDELS data otherwise) were used to estimate the galaxy size and hence the volume enclosing most of the ISM under simple geometrical assumptions. These measurements were then combined to derive the column density associated with the ISM of the host, on the order of N$_{H, ISM}$~10$^{23-24}$ cm$^{-2}$. The comparison between the ISM column densities and those measured from the X-ray spectral analysis shows that they are similar. This suggests that, at least at high redshift, significant absorption on kpc scales by the dense ISM in the host likely adds to or substitutes that produced by circumnuclear gas on pc scales (i.e., the torus of unified models). The lack of unobscured AGN among our ISM-rich targets supports this scenario., Comment: 15 pages, 3 figures. Accepted for publication in A&A

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., Comment: 28 pages, 19 figures. Forthcoming article in A&A

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

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., Comment: Accepted for publication on Astronomy and Astrophysics. 22 pages, 14 figures

Although studying outflows in the host galaxies of AGN have become the forefront of extra-galactic astronomy in recent years, estimating the energy associated with these outflows have been a major challenge. Determination of the energy associated with an outflow often involves an assumption of uniform density in the NLR, which span a wide range in literature leading to large systematic uncertainties in energy estimation. In this paper, we present electron density maps for a sample of outflowing and non-outflowing Seyfert galaxies at z<0.02 drawn from the S7 survey and understand the origin and values of the observed density structures to reduce the systematic uncertainties in outflow energy estimation. We use the ratio of the [SII]6716,6731 emission lines to derive spatially resolved electron densities (<50-2000 cm$^{-3}$). Using optical IFU observations, we are able to measure densities across the central 2-5 kpc of the selected AGN host galaxies. We compare the density maps with the positions of the HII regions derived from the narrow H$\alpha$ component, ionization maps from [OIII], and spatially resolved BPT diagrams, to infer the origin of the observed density structures. We also use the electron density maps to construct density profiles as a function of distance from the central AGN. We find a spatial correlation between the sites of high star formation and high electron density for targets without an active ionized outflow. The non-outflowing targets also show an exponential drop in the electron density as a function of distance from the center, with a mean exponential index of ~0.15. The correlation between the star forming sites and electron density ceases for targets with an outflow. The density within the outflowing medium is not uniform and shows both low and high density sites, most likely due to the presence of shocks and highly turbulent medium., Comment: Accepted for publication in A&A, 16 pages, 10 figures, 2 tables

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged], Comment: Submitted to the Astrophysical Journal Supplement Series. 64 pages, 36 figures. The reduced data sets will be made available once the paper is accepted for publication. A version with full resolution Figures is available at http://www.mpe.mpg.de/~forster/FS18_AOsurvey_ApJSsubm.html

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

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., Comment: 15 pages, 11 figures (including Appendix). Version matches A&A published version (after proof)

We present a comprehensive study of the Na I $\lambda$5890, 5895 (Na I D) resonant lines in the Sloan Digital Sky Survey (SDSS, DR7) spectroscopic sample to look for neutral gas outflows in the local galaxies. Individual galaxy spectra are stacked in bins of M${\star}$ and SFR to investigate the dependence of galactic wind occurrence and velocity as a function of the galaxy position in the SFR-$M{\star}$ plane. In massive galaxies at the high SFR tail we find evidence of a significant blue-shifted Na I D absorption, which we interpret as evidence of neutral outflowing gas. The occurrence of the blue-shifted absorption is observed at the same significance for purely SF galaxies, AGN and composite systems at fixed SFR. In all classes of objects the blue-shift is the largest and the Na I D equivalent width the smallest for face-on galaxies while the absorption feature is at the systemic velocity for edge-on systems. This indicates that the neutral outflow is mostly perpendicular or biconical with respect to the galactic disk. We also compare the kinematics of the neutral gas with the ionized gas phase as traced by the [OIII]$\lambda$5007, H$\alpha$, [NII]$\lambda6548$ and [NII]$\lambda6584$ emission lines. Differently for the neutral gas phase, all the emission lines show evidence of perturbed kinematics only in galaxies with a significant level of nuclear activity and, they are independent from the disk inclination. In conclusion, we find that, in the local Universe, galactic winds show two faces which are related to two different ejection mechanisms, namely the neutral outflowing gas phase related to the SF activity along the galaxy disk and the ionized phase related to the AGN feedback. In both the neutral and ionized gas phases, the observed wind velocities suggest that the outflowing gas remains bound to the galaxy with no definitive effect on the gas reservoir., Comment: Accepted to A&A, 13 pages, 9 figures

Active Galactic Nuclei (AGN) are energetic astrophysical sources powered by accretion onto supermassive black holes in galaxies, and present unique observational signatures that cover the full electromagnetic spectrum over more than twenty orders of magnitude in frequency. The rich phenomenology of AGN has resulted in a large number of different "flavours" in the literature that now comprise a complex and confusing AGN "zoo". It is increasingly clear that these classifications are only partially related to intrinsic differences between AGN, and primarily reflect variations in a relatively small number of astrophysical parameters as well the method by which each class of AGN is selected. Taken together, observations in different electromagnetic bands as well as variations over time provide complementary windows on the physics of different sub-structures in the AGN. In this review, we present an overview of AGN multi-wavelength properties with the aim of painting their "big picture" through observations in each electromagnetic band from radio to gamma-rays as well as AGN variability. We address what we can learn from each observational method, the impact of selection effects, the physics behind the emission at each wavelength, and the potential for future studies. To conclude we use these observations to piece together the basic architecture of AGN, discuss our current understanding of unification models, and highlight some open questions that present opportunities for future observational and theoretical progress., Comment: Accepted for publication in Astronomy & Astrophysics Review, 56 pages, 25 figures

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., Comment: 14 pages, 13 figures, 4 tables, accepted for publication in A&A

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

Similarly to the cosmic star formation history, the black hole accretion rate density of the Universe peaked at 199% 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., Comment: 12 pages, 5 figures, 3 tables, accepted for publication in MNRAS

We present rest-frame optical spectra from the FMOS-COSMOS survey of twelve $z \sim 1.6$ \textit{Herschel} starburst galaxies, with Star Formation Rate (SFR) elevated by $\times$8, on average, above the star-forming Main Sequence (MS). Comparing the H$\alpha$ to IR luminosity ratio and the Balmer Decrement we find that the optically-thin regions of the sources contain on average only $\sim 10$ percent of the total SFR whereas $\sim90$ percent comes from an extremely obscured component which is revealed only by far-IR observations and is optically-thick even in H$\alpha$. We measure the [NII]$_{6583}$/H$\alpha$ ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity-SFR anticorrelation observed at fixed stellar mass for the MS population. The [SII]$_{6732}$/[SII]$_{6717}$ ratio from the average spectrum indicates an electron density $n_{\rm e} \sim 1,100\ \mathrm{cm}^{-3}$, larger than what estimated for MS galaxies but only at the 1.5$\sigma$ level. Our results provide supporting evidence that high-$z$ MS outliers are the analogous of local ULIRGs, and are consistent with a major merger origin for the starburst event., Comment: 6 pages, 4 figures, Accepted for publication in ApJ Letters

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., Comment: A&A in press

We analyze a complete spectroscopic sample of galaxies ($\sim$600,000 ) drawn from Sloan Digital Sky Survey (SDSS, DR7) to look for evidence of galactic winds in the local Universe. We focus on the shape of the [OIII]$\lambda$5007 emission line as a tracer of ionizing gas outflows. We stack our spectra in a fine grid of star formation rate (SFR) and stellar mass to analyze the dependence of winds on the position of galaxies in the SFR versus mass diagram. We do not find any significant evidence of broad and shifted [OIII]$\lambda$5007 emission line which we interpret as no evidence of outflowing ionized gas in the global population. We have also classified these galaxies as star-forming or AGN dominated according to their position in the standard BPT diagram. We show how the average [OIII]$\lambda$5007 profile changes as function of nature of the dominant ionizing source. We find that in the star-forming dominated source the oxygen line is symmetric and governed by the gravitational potential well. The AGN or composite AGN$\setminus$star-formation activity objects, in contrast, display a prominent and asymmetric profile that can be well described by a broad gaussian component that is blue-shifted from a narrow symmetric core. In particular, we find that the blue wings of the average [OIII]$\lambda$5007 profiles are increasingly prominent in the LINERs and Seyfert galaxies. We conclude that, in the low-redshift Universe, "pure" star-formation activity does not seem capable of driving ionized-gas outflows, while, the presence of optically selected AGN seems to play a primary role to drive such winds. We discuss the implications of these results for the role of the quenching mechanism in the present day Universe., Comment: Re-submitted to A&A after addressing referee's comments,14 pages, 17 figures, 1 table