Search

While observations of molecular gas at cosmic noon and beyond have focused on the gas within galaxies (i.e., the interstellar medium; ISM), it is also crucial to study the molecular gas reservoirs surrounding each galaxy (i.e., in the circumgalactic medium; CGM). Recent observations of galaxies and quasars hosts at high redshift (z>2) have revealed evidence for cold gaseous halos of scale r_CGM~10kpc, with one discovery of a molecular halo with r_CGM~200kpc and a molecular gas mass one order of magnitude larger than the ISM of the central galaxy. As a follow-up, we present deep ACA and ALMA observations of CO(3-2) from this source and two other quasar host galaxies at z~2.2. While we find evidence for CO emission on scales of r~10kpc, we do not find evidence for molecular gas on scales larger than r>20 kpc. Therefore, our deep data do not confirm the existence of massive molecular halos on scales of ~100 kpc for these X-ray selected quasars. As an interesting by-product of our deep observations, we obtain the tentative detection of a negative continuum signal on scales larger than r>200kpc, which might be tracing the Sunyaev-Zeldovich effect associated with the halo heated by the active galactic nucleus (AGN). If confirmed with deeper data, this could be direct evidence of the preventive AGN feedback process expected by cosmological simulations., 17 pages, 12 figures. Accepted for publication in MNRAS

The final phase of the reionization process can be probed by rest-frame UV absorption spectra of quasars at z>6, shedding light on the properties of the diffuse intergalactic medium within the first Gyr of the Universe. The ESO Large Programme "XQR-30: the ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" dedicated ~250 hours of observations at the VLT to create a homogeneous and high-quality sample of spectra of 30 luminous quasars at z~6, covering the rest wavelength range from the Lyman limit to beyond the MgII emission. Twelve quasar spectra of similar quality from the XSHOOTER archive were added to form the enlarged XQR-30 sample, corresponding to a total of ~350 hours of on-source exposure time. The median effective resolving power of the 42 spectra is R~11400 and 9800 in the VIS and NIR arm, respectively. The signal-to-noise ratio per 10 km/s pixel ranges from ~11 to 114 at $\lambda \simeq 1285$ \AA rest frame, with a median value of ~29. We describe the observations, data reduction and analysis of the spectra, together with some first results based on the E-XQR-30 sample. New photometry in the H and K bands are provided for the XQR-30 quasars, together with composite spectra whose characteristics reflect the large absolute magnitudes of the sample. The composite and the reduced spectra are released to the community through a public repository, and will enable a range of studies addressing outstanding questions regarding the first Gyr of the Universe., Comment: 21 pages, 10 figures. Final version accepted by MNRAS

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

ALMA observations have revealed that many high redshift galaxies are surrounded by extended (10-15 kpc) [CII]-emitting halos which are not predicted by even the most advanced zoom-in simulations. Using a semi-analytical model, in a previous work we suggested that such halos are produced by starburst-driven, catastrophically cooling outflows. Here, we further improve the model and compare its predictions with data from 7 star-forming ($10\lesssim \rm SFR/ M_\odot \rm yr^{-1}, Comment: Accepted for publication in MNRAS. 11 pages, 6 figures

Large-scale outflows are believed to be an important mechanism in the evolution of galaxies. We can determine the impact of these outflows by studying either current galaxy outflows and their effect in the galaxy or by studying the effect of past outflows on the gas surrounding the galaxy. In this work, we examine the CO(7−6), [C i] ($^{3} \rm P_{1} \rightarrow {\rm ^3 P}_{0}$), H2O 211–202, and dust continuum emission of 15 extremely red quasars at z ∼ 2.3 using ALMA. By investigating the radial surface brightness profiles of both the individual sources and the stacked emission, we detect extended cold gas and dust emission on scales of ∼14 kpc in CO(7−6), [C i](2−1), and dust continuum. This is the first time that the presence of a large amount of molecular gas was detected on large, circumgalactic medium scales around quasar host galaxies using [C i] extended emission. We estimate the dust and molecular gas mass of these haloes to be 107.6 and 1010.6 M⊙, indicating significant dust and molecular gas reservoirs around these extreme quasars. By estimating the time-scale at which this gas can reach these distances by molecular gas outflows (7–32 Myr), we conclude that these haloes are a relic of past AGN or starburst activity, rather than an effect of the current episode of extreme quasar activity.

At redshift z > 5, the far-infrared (FIR) continuum spectra of main-sequence galaxies are sparsely sampled, often with a single data point. The dust temperature Td,SED, thus has to be assumed in the FIR continuum fitting. This introduces large uncertainties regarding the derived dust mass (Md), FIR luminosity, and obscured fraction of the star formation rate. These are crucial quantities to quantify the effect of dust obscuration in high-z galaxies. To overcome observation limitations, we introduce a new method that combines dust continuum information with the overlying [C $\scriptstyle \rm II$] 158 µm line emission. By breaking the Md–Td,SED degeneracy, with our method, we can reliably constrain the dust temperature with a single observation at 158 µm. This method can be applied to all Atacama Large Millimeter Array (ALMA) and NOEMA [C $\scriptstyle \rm II$] observations, and exploited in ALMA Large Programs such as ALPINE and REBELS targeting [C $\scriptstyle \rm II$] emitters at high-z. We also provide a physical interpretation of the empirical relation recently found between molecular gas mass and [C $\scriptstyle \rm II$] luminosity. We derive an analogous relation linking the total gas surface density and [C $\scriptstyle \rm II$] surface brightness. By combining the two, we predict the cosmic evolution of the surface density ratio $\Sigma _{\rm H_2} / \Sigma _{\rm gas}$. We find that $\Sigma _{\rm H_2} / \Sigma _{\rm gas}$ slowly increases with redshift, which is compatible with current observations at 0 < z < 4.

We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. In dusty (${M_{\rm d}\gtrsim 3\times 10^7 M_{\rm \odot}}$) galaxies, $\gtrsim 50-90 \%$ of the UV radiation is obscured by dust inhomogeneities on scales ${\gtrsim 100}$ pc. In runs with AGN, a clumpy, warm ($\approx 250$ K) dust component co-exists with a colder ($\approx 60$ K) and more diffuse one, heated by stars. Warm dust provides up to ${50 \%}$ of the total IR luminosity, but only $\lesssim 0.1 \%$ of the total mass content. The AGN boosts the MIR flux by ${10-100 \times}$ with respect to star forming galaxies, without significantly affecting the FIR. Our simulations successfully reproduce the observed SED of bright (${M_{\rm UV}\sim -26}$) ${z\sim 6}$ quasars, and show that these objects are part of complex, dust-rich merging systems, containing multiple sources (accreting BHs and/or star forming galaxies) in agreement with recent HST and ALMA observations. Our results show that the proposed ORIGINS missions will be able to investigate the MIR properties of dusty star forming galaxies and to obtain good quality spectra of bright quasars at $z\sim 6$. Finally, the MIR-to-FIR flux ratio of faint (${M_{\rm UV}\sim -24}$) AGN is ${>10\times}$ higher than for normal star forming galaxies. This implies that combined JWST/ORIGINS/ALMA observations will be crucial to identify faint and/or dust-obscured AGN in the distant Universe., 22 pages, 13 figures, accepted for publication on MNRAS

We investigate the spatially-resolved morphology of galaxies in the early Universe. We consider a typical redshift z = 6 Lyman Break galaxy, "Althaea" from the SERRA hydrodynamical simulations. We create mock rest-frame ultraviolet, optical, and far-infrared observations, and perform a two-dimensional morphological analysis to de-blend the galaxy disk from substructures (merging satellites or star-forming regions). We find that the [CII]158um emitting region has an effective radius 1.5 - 2.5 times larger than the optical one, consistent with recent observations. This [CII] halo in our simulated galaxy arises as the joint effect of stellar outflows and carbon photoionization by the galaxy UV field, rather than from the emission of unresolved nearby satellites. At the typical angular resolution of current observations (> 0.15") only merging satellites can be detected; detection of star-forming regions requires resolutions of < 0.05". The [CII]-detected satellite has a 2.5 kpc projected distance from the galaxy disk, whereas the star-forming regions are embedded in the disk itself (distance < 1 kpc). This suggests that multi-component systems reported in the literature, which have separations > 2 kpc, are merging satellites, rather than galactic substructures. Finally, the star-forming regions found in our mock maps follow the local L[CII] - SFR_UV relation of galaxy disks, although sampling the low-luminosity, low-SFR tail of the distribution. We show that future JWST observations, bridging UV and [CII] datasets, will be exceptionally suited to characterize galaxy substructures thanks to their exquisite spatial resolution and sensitivity to both low-metallicity and dust-obscured regions that are bright at infrared wavelengths., Accepted for publication in MNRAS; 17 pages (plus appendix), 7 figures, 4 tables

Aims. We investigate the ionizing properties of the pair of bright Lyα emitting galaxies BDF521 and BDF2195 at z = 7.012 in order to constrain their contribution to the formation of the Bremer Deep Field (BDF) ‘reionized bubble’ in which they have been shown to reside. Methods. We obtain constraints on four UV emission lines (the CIVλ1548 doublet, HeIIλ1640, the OIII]λ1660 doublet, and the CIII]λ1909 doublet) from deep VLT X-shooter observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and active galactic nucleus (AGN) emission. We use this spectroscopic information, together with the photometry available in the field, to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. Results. We do not detect any significant emission at the expected position of the UV lines, with 3σ upper limits of equivalent width (EW) ≲2–7 Å rest-frame. We find that the two objects have a lower CIII] emission than expected on the basis of the correlation between the Lyα and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at ∼2 − 3σ significance, and only models with a ≲40% AGN contribution are compatible with the observations. The two objects are found to be relatively young (∼20–30 Myrs) and metal-poor (≲0.3 Z⊙), with stellar masses of a few 109 M⊙. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log(ξion*/Hz erg−1) = 25.02–25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in z > 7 galaxies with strong CIII] and/or optical line emission (≃25.6–25.7). Conclusions. The two BDF emitters show no evidence of higher-than-average ionizing capabilities and are not capable of reionizing their surroundings by their own means, under realistic assumptions of the escape fraction of ionizing photons. Therefore, a dominant contribution to the formation of the reionized bubble must have been provided by fainter companion galaxies. The capabilities of the James Webb Space Telescope will be needed for spectroscopic confirmation of these objects.

Active galactic nuclei (AGN) feedback has a major impact onto the supermassive black-hole (SMBH) growth, the properties of the host galaxies, and their cosmic evolution. We investigate the effects of different kinetic feedback prescriptions on the observable properties of AGN and their host galaxies at $z>6$ in a suite of zoom-in cosmological simulations. We find that kinetic feedback decreases the column density of the interstellar medium (ISM) in the host galaxy by up to a factor of $\approx10$, especially when the SMBHs reach high accretion rates ($\approx10-30\,\mathrm{M_\odot\,yr^{-1}}$). In particular, kinetic feedback is required to extend the ISM size to $>1$ kpc and match the observed sizes of the gas reservoirs in $z>6$ AGN host galaxies. Moreover, it produces unobscured lines of sight along which the AGN can be detected in the rest-frame UV band with magnitudes consistent with observed values of $z>6$ AGN. The assumed geometry of the outflow plays an important role in shaping the observed properties of high-redshift AGN. We find that a biconical geometry is favored over a spherical one to reproduce the observed properties, but it overestimates the number of multiple AGN systems detectable in X-ray observations. This result suggests that simplistic BH seeding recipes widely employed in cosmological simulations produce too many X-ray detectable multiple AGN at $z=6-7$, thus soliciting the adoption of more physically motivated seeding prescriptions., Comment: Accepted for publication on MNRAS, 18 pages

Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, although direct observational evidence is still scarce and debated. Here we present Early Release Science JWST NIRSpec IFU observations of the z=1.59 prototypical obscured Active Galactic Nucleus (AGN) XID2028: This target represents a unique test case for studying quasar feedback at the peak epoch of AGN-galaxy co-evolution because extensive multi-wavelength coverage is available and a massive and extended outflow is detected in the ionised and molecular components. With the unprecedented sensitivity and spatial resolution of the JWST, the NIRSpec dataset reveals a wealth of structures in the ionised gas kinematics and morphology that were previously hidden in the seeing-limited ground-based data. In particular, we find evidence of an interaction between the interstellar medium of the galaxy and the quasar-driven outflow and radio jet that produces an expanding bubble from which the fast and extended wind detected in previous observations emerges. The new observations confirm the complex interplay between the AGN jet, wind and the interstellar medium of the host galaxy, highlighting the role of low-luminosity radio jets in AGN feedback. They also clearly show the new window that NIRSpec opens for detailed studies of feedback at high redshift., 12 pages, 11 figures, accepted for publication in A&A

Surveys with James Webb Space Telescope (JWST) have discovered candidate galaxies in the first 400 Myr of cosmic time. Preliminary indications have suggested these candidate galaxies may be more massive and abundant than previously thought. However, without confirmed distances, their inferred properties remain uncertain. Here we identify four galaxies located in the JWST Advanced Deep Extragalactic Survey (JADES) Near-Infrared Camera (NIRCam) imaging with photometric redshifts z~10-13. These galaxies include the first redshift z>12 systems discovered with distances spectroscopically confirmed by JWST in a companion paper. Using stellar population modelling, we find the galaxies typically contain a hundred million solar masses in stars, in stellar populations that are less than one hundred million years old. The moderate star formation rates and compact sizes suggest elevated star formation rate surface densities, a key indicator of their formation pathways. Taken together, these measurements show that the first galaxies contributing to cosmic reionisation formed rapidly and with intense internal radiation fields., Author version of manuscript, please visit Nature Astronomy for the version published 04 April 2023

Recent ALMA large surveys unveiled the presence of significant dust continuum emission in star-forming galaxies at $z>4$. Unfortunately, such large programs -- i.e. ALPINE ($z\sim 5$) and REBELS ($z \sim 7$) -- only provide us with a single Far-Infrared (FIR) continuum data point for their individual targets. Therefore, high-$z$ galaxies FIR spectral energy densities (SEDs) remain mostly unconstrained, hinging on an assumption for their dust temperature ($T_{\rm d}$) in the SED fitting procedure. This introduces uncertainties in the inferred dust masses ($M_{\rm d }$), infrared luminosities ($L_{\rm IR}$), and obscured Star Formation Rate (SFR) fraction at $z > 4$. In this work we use a method that allows us to constrain $T_{\rm d}$ with a single band measurement by combining the $158\ \mathrm{μm}$ continuum information with the overlying [CII] emission line. We analyse the $21$ [CII] and FIR continuum detected $z\sim 5$ galaxies in ALPINE, finding a range of $T_{\rm d}=25-60\ \mathrm{K}$ and $M_{\rm d} = 0.6-25.1\ \times 10^{7}\ \mathrm{M_{\odot}}$. Given the measured stellar masses of ALPINE galaxies, the inferred dust yields are around $M_{\rm d}/M_{\star} = (0.2-8) \times 10^{-3}$, consistent with theoretical dust-production constraints. We find that $8$ out of $21$ ALPINE galaxies have $L_{\rm IR} \geq 10^{12}\ \mathrm{L_{\odot}}$, comparable to UltraLuminous IR Galaxies (ULIRGs). Relying on ultraviolet-to-optical SED fitting, the SFR was underestimated by up to $2$ orders of magnitude in $4$ of these $8$ ULIRGs-like galaxies. We conclude that these $4$ peculiar sources should be characterised by a two-phase interstellar medium structure with "spatially-segregated" FIR and ultraviolet emitting regions., Accepted for publication in MNRAS. 13 pages, 4 figures

We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 202 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $\kappa_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios can be reproduced by a part of the SERRA galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$)., Comment: 23 pages, 14 figures, accepted by MNRAS; part of the data used for this study is available at the website http://cosmology.sns.it/data_access.html

The existence of multiple active galactic nuclei (AGN) at small projected distances on the sky is due to either the presence of multiple, in-spiraling SMBHs, or to gravitational lensing of a single AGN. Both phenomena allow us to address important astrophysical and cosmological questions. However, few kpc-separation multiple AGN are currently known. Recently, the newly-developed Gaia Multi peak (GMP) method provided numerous new candidate members of these populations. We present spatially resolved, integral-field spectroscopy of a sample of four GMP-selected multiple AGNs candidates. In all of these systems, we detect two or more components with sub-arcsec separations. We find that two of the systems are dual AGNs, one is either an intrinsic triple or a lensed dual AGN, while the last system is a chance AGN/star alignment. Our observations double the number of confirmed multiple AGNs at projected separations below 7 kpc at z > 0.5, present the first detection of a possible triple AGN in a single galaxy at z > 0.5, and successfully test the GMP method as a novel technique to discover previously unknown multiple AGNs., 6 pages, 2 figures

The Near-Infrared Spectrograph (NIRSpec) is one of the four focal plane instruments on the James Webb Space Telescope. In this paper, we summarize the in-orbit performance of NIRSpec, as derived from data collected during its commissioning campaign and the first few months of nominal science operations. More specifically, we discuss the performance of some critical hardware components such as the two NIRSpec Hawaii-2RG (H2RG) detectors, wheel mechanisms, and the micro-shutter array. We also summarize the accuracy of the two target acquisition procedures used to accurately place science targets into the slit apertures, discuss the current status of the spectro-photometric and wavelength calibration of NIRSpec spectra, and provide the as measured sensitivity in all NIRSpec science modes. Finally, we point out a few important considerations for the preparation of NIRSpec science programs., accepted by PASP for special issue on JWST in-orbit performance

Baryonic feedback is expected to play a key role in regulating the star formation of low-mass galaxies by producing galaxy-scale winds associated with mass-loading factors of β ∼ 1 − 50. We test this prediction using a sample of 19 nearby systems with stellar masses of 107 M⋆/M⊙ 10, mostly lying above the main sequence of star-forming galaxies. We used MUSE at VLT optical integral field spectroscopy to study the warm ionised gas kinematics of these galaxies via a detailed modelling of their Hα emission line. The ionised gas is characterised by irregular velocity fields, indicating the presence of non-circular motions of a few tens of km s−1 within galaxy discs, but with intrinsic velocity dispersion of 40 − 60 km s−1 that are only marginally larger than those measured in main-sequence galaxies. Galactic winds, defined as gas at velocities larger than the galaxy escape speed, encompass only a few percent of the observed fluxes. Mass outflow rates and loading factors are strongly dependent on M⋆, the star formation rate (SFR), SFR surface density, and specific SFR (sSFR). For M⋆ of 108 M⊙ we find β ≃ 0.02, which is more than two orders of magnitude smaller than the values predicted by theoretical models of galaxy evolution. In our galaxy sample, baryonic feedback stimulates a gentle gas cycle rather than causing a large-scale blow-out.

We investigate the attenuation law in $z\sim 6$ quasars by combining cosmological zoom-in hydrodynamical simulations of quasar host galaxies, with multi-frequency radiative transfer calculations. We consider several dust models differing in terms of grain size distributions, dust mass and chemical composition, and compare the resulting synthetic Spectral Energy Distributions (SEDs) with data from bright, early quasars. We show that only dust models with grain size distributions in which small grains ($a < 0.1~\mu$m, corresponding to $\approx 60\%$ of the total dust mass) are selectively removed from the dusty medium provide a good fit to the data. Removal can occur if small grains are efficiently destroyed in quasar environments and/or early dust production preferentially results in large grains. Attenuation curves for these models are close to flat, and consistent with recent data; they correspond to an effective dust-to-metal ratio $f_d \simeq 0.38$, i.e. close to the Milky Way value., STFC ERC

We study the impact of deviations from the Kennicutt-Schmidt relation (quantified by the `burstiness' parameter $\kappa_s$), gas metallicity ($Z$), and density ($n$) on the observed [OIII]88$\mu$m/[CII]158$\mu$m surface brightness ratios ($\Sigma_{[OIII]}/\Sigma_{[CII]}$) in nine galaxies at $z\approx6-9$. We first discuss possible biases in the measured $\Sigma_{[OIII]}/\Sigma_{[CII]}$ ratios by comparing the data with zoom-in cosmological simulations, and then use a Markov Chain Monte Carlo algorithm to derive the best fit values of ($\kappa_s, Z, n$). We find that (i) the strongest dependence of $\Sigma_{[OIII]}/\Sigma_{[CII]}$ is on $\kappa_s$; (ii) high ratios identify starburst galaxies with short gas depletion times ($t_{dep}=6-49\,\rm Myr$); (iii) a secondary dependence on density is found, with $\Sigma_{[OIII]}/\Sigma_{[CII]}$ anticorrelating with $n$ as a result of the lower [OIII] critical density, (iv) the ratio only weakly depends on $Z$. The nine galaxies are significantly enriched (Z=0.2-0.5 $Z_\odot$), and dense ($n=10^{1-3} {\rm cm}^{-3}$). This lends further support to the starburst scenario in which a rapid enrichment of the interstellar medium is expected., Comment: accepted for publication in MNRAS, 12 pages, 8 figures

We present recent results from our MAGNUM survey of nearby active galactic nuclei (AGN), which exploits observations from the optical/near-IR integral field spectrograph MUSE at VLT. We detect strongly enhanced line widths in emission line maps of four galaxies perpendicularly to their low-power jets and AGN ionisation cones, indicative of turbulent/outflowing material. The observation of a similar phenomenon in other works suggests that it originates from an interaction mechanism between the jet and the galaxy disc through which it propagates.

We study the structure of spatially resolved, line-of-sight velocity dispersion for galaxies in the Epoch of Reionization (EoR) traced by [CII] $158\mu\rm{m}$ line emission. Our laboratory is a simulated prototypical Lyman-break galaxy, "Freesia", part of the SERRA suite. The analysis encompasses the redshift range 6 < z < 8, when Freesia is in a very active assembling phase. We build velocity dispersion maps for three dynamically distinct evolutionary stages (Spiral Disk at z=7.4, Merger at z=8.0, and Disturbed Disk at z=6.5) using [CII] hyperspectral data cubes. We find that, at a high spatial resolution of 0.005" ($\simeq 30 pc$), the luminosity-weighted average velocity dispersion is $\sigma_{\rm{CII}}$~23-38 km/s with the highest value belonging to the highly-structured Disturbed Disk stage. Low resolution observations tend to overestimate $\sigma_{\rm CII}$ values due to beam smearing effects that depend on the specific galaxy structure. For an angular resolution of 0.02" (0.1"), the average velocity dispersion is 16-34% (52-115%) larger than the actual one. The [CII] emitting gas in Freesia has a Toomre parameter $\mathcal{Q}$~0.2 and a rotational-to-dispersion ratio of $v_{\rm c}/\sigma$~ 7 similar to that observed in z=2-3 galaxies. The primary energy source for the velocity dispersion is due to gravitational processes, such as merging/accretion events; energy input from stellar feedback is generally subdominant (< 10%). Finally, we find that the resolved $\sigma_{\rm{CII}} - {\Sigma}_{\rm SFR}$ relation is relatively flat for $0.02, Comment: Accepted for publication in MNRAS; 16 pages (plus appendix), 8 figures, 3 tables

Context. Sources at the brightest end of the quasi-stellar object (QSO) luminosity function, during the peak epoch in the history of star formation and black hole accretion (z ∼ 2-4, often referred to as "Cosmic noon") are privileged sites to study the cycle of feeding & feedback processes in massive galaxies. Aims. We aim to perform the first systematic study of cold gas properties in the most luminous QSOs, by characterising their host-galaxies and environment. These targets exhibit indeed widespread evidence of outflows at nuclear and galactic scales. Methods. We analyse ALMA, NOEMA and JVLA observations of the far-infrared continuum, CO and [CII] emission lines in eight QSOs (bolometric luminosity LBol ≳ 3 × 1047 erg s-1) from the WISE-SDSS selected hyper-luminous (WISSH) QSOs sample at z ∼ 2.4-4.7. Results. We report a 100% emission line detection rate and a 80% detection rate in continuum emission, and we find CO emission to be consistent with the steepest CO ladders observed so far. Sub-millimetre data reveal presence of (one or more) bright companion galaxies around ∼80% of WISSH QSOs, at projected distances of ∼6-130 kpc. We observe a variety of sizes for the molecular gas reservoirs (∼1.7-10 kpc), mostly associated with rotating disks with disturbed kinematics. WISSH QSOs typically show lower CO luminosity and higher star formation efficiency than infrared matched, z ∼ 0-3 main-sequence galaxies, implying that, given the observed SFR ∼170-1100 M⊙ yr-1, molecular gas is converted into stars in ≲ 50 Myr. Most targets show extreme dynamical to black-hole mass ratios Mdyn/MBH ∼ 3-10, two orders of magnitude smaller than local relations. The molecular gas fraction in the host-galaxies of WISSH is lower by a factor of ∼10-100 than in star forming galaxies with similar M∗. Conclusions. Our analysis reveals that hyper-luminous QSOs at Cosmic noon undergo an intense growth phase of both the central super-massive black hole and of the host-galaxy. These systems pinpoint the high-density sites where giant galaxies assemble, where we show that mergers play a major role in the build-up of the final host-galaxy mass. We suggest that the observed low molecular gas fraction and short depletion timescale are due to AGN feedback, whose presence is indicated by fast AGN-driven ionised outflows in all our targets. © ESO 2020., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00417.S, ADS/JAO.ALMA#2019.1.01070.S, ADS/JAO.ALMA#2015.1.01602.S, ADS/JAO.ALMA#2016.1.00718.S, ADS/JAO. ALMA#2016.1.01515.S, ADS/JAO.ALMA#2018.1.01806.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work is based on observations carried out under project number S17BW and W17DT with the IRAM NOEMA Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). We are grateful to the anonymous referee for useful feedback which helped us to improve the paper. We thank Dr. Y. Wang for providing us with SUBARU images of J0209-0005 and Dr. P. Santini for useful discussion about SFE in high-z galaxies. MBi, MBr, CF, FF, AM, EP acknowledge support from PRIN MIUR project "Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper", contract #2017PH3WAT. MBi, CF and FF acknowledge support from INAF under PRIN SKA/CTA FORECaST and PRIN MAINSTREAM 2018 "Black hole winds and the baryon cycle". EP, LZ and MBi acknowledge financial support under ASI-INAF contract 2017-14-H.0. GB acknowledges financial support under the INTEGRAL ASI-INAF agreement 2019-35-HH.0 MPT acknowledges financial support from the State Agency for Research of the Spanish MCIU through the Center of Excellence Severo Ochoa award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709) and through grant PGC2018-098915-B-C21 (MCI/AEI/FEDER, UE). RM acknowledges ERC Advanced Grant 695671 "QUENCH" and support by the Science and Technology Facilities Council (STFC). GVe acknowledges support from CONICYT Basal-CATA AFB-170002 and FONDECYT Postdoctorado 3200802 grants. CC acknowledges support from the Royal Society. MP is supported by the Programa Atraccion de Talento de la Comunidad de Madrid via grant 2018-T2/TIC-11715. GVi acknowledges financial support from Premiale 2015 MITiC (PI B. Garilli).

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-$��$ relation, (b) alleviates the problem of the uncomfortably large dust masses deduced from observations of some EoR galaxies., 14 pages, 6 figures, accepted for publication in MNRAS

ALMA observations have revealed that [CII] 158μm line emission in high-z galaxies is ~2-3× 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[OIII]/L[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′′, implying that L[CII] is underestimated by a factor ≈2 in data at low (6 galaxies lies, on average, slightly below the local L[CII]−SFR relation (Δz=6−9=−0.07±0.3), but within the intrinsic dispersion of the relation. SBD correction also yields L[OIII]/L[CII], STFC ERC

Aims. 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 AGNs at redshift z∼2 with high (∼2 kpc) spatial resolution, and to correlate them with the properties of their host galaxy and central black hole. The main aims of this paper are: (a) to derive reliable estimates for the masses of the black holes and accretion rates for the Type-1 AGNs in this survey; and (b) to characterise the properties of the AGN-driven winds in the broad line region (BLR). Methods. We analysed rest-frame optical and UV spectra of 21 Type-1 AGNs. We used Hα, Hβ, and MgII line profiles to estimate the masses of the black holes. We used the blueshift of the CIV line profile to trace the presence of winds in the BLR. Results. We find that the Hα and Hβ line widths are strongly correlated, as is the line continuum luminosity at 5100 Å with Hα line luminosity, resulting in a well-defined correlation between black hole masses estimated from Hα and Hβ. Using these lines, we estimate that the black hole masses for our objects are in the range Log (MBH/M·) = 8.4-10.8 and are accreting at λEdd = 0.04-1.3. Furthermore, we confirm the well-known finding that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blueshifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR for which we estimated velocities up to ∼4700 km s-1. We confirm the strong dependence of the BLR wind velocity on the UV-to-X-ray continuum slope, the bolometric luminosity, and Eddington ratio. We infer BLR mass outflow rates in the range 0.005-3 M· yr-1, revealing a correlation with the bolometric luminosity consistent with that observed for ionised winds in the narrow line region (NLR), and X-ray winds detected in local AGNs, and kinetic power ∼10-7-10-4 × LBol. The coupling efficiencies predicted by AGN-feedback models are much higher than the values reported for the BLR winds in the SUPER sample; although it should be noted that only a fraction of the energy injected by the AGN into the surrounding medium is expected to become kinetic power in the outflow. Finally, we find an anti-correlation between the equivalent width of the [OIII] line and the CIV velocity shift, and a positive correlation between this latter parameter and [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGNs at z∼2, support a scenario where BLR winds are connected to galaxy-scale detected outflows, and are therefore capable of affecting the gas in the NLR located at kiloparsec scale distances., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

We present ALMA observations of a merging system atz ∼ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII]158 μm and far-infrared continuum data with multi-wavelength ancillary data, we find that the system is composed of two massive (M⋆ ≳ 1010 M⊙) star-forming galaxies experiencing a major merger (stellar mass ratiormass ≳ 0.9) at close spatial (∼13 kpc; projected) and velocity (Δv < 300 km s−1) 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. Additionally, ALMA reveals the presence of [CII] emission arising from a circumgalactic gas structure, extending up to a diameter-scale of ∼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 characterised 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-zgalaxies, and thus play a key role in the galaxy baryon cycle at early epochs.

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%., Accepted for publication in A&A, 34 pages, 14 figures and 5 tables

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., Comment: Accepted for publication in MNRAS. 12 pages, 6 Figures. Comments welcome

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., Comment: 14 pages, 12 figures, accepted for publication in MNRAS

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., Accepted 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, Comment: 28 pages, 8 figures, 5 tables. Accepted for publication in A&A

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��$ 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., 13 pages, 9 figures, 3 tables; MNRAS journal