Showing 1–50 of 130 results for author: Feruglio, C

SUPER VIII. Fast and Furious at $z\sim2$: obscured type-2 active nuclei host faster ionised winds than type-1 systems

Authors: G. Tozzi, G. Cresci, M. Perna, V. Mainieri, F. Mannucci, A. Marconi, D. Kakkad, A. Marasco, M. Brusa, E. Bertola, M. Bischetti, S. Carniani, C. Cicone, C. Circosta, F. Fiore, C. Feruglio, C. M. Harrison, I. Lamperti, H. Netzer, E. Piconcelli, A. Puglisi, J. Scholtz, G. Vietri, C. Vignali, G. Zamorani

Abstract: 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 exami… ▽ More 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]$λ$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. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 21 pages, 13 figures. Key figure is 8. Accepted for publication in A&A

The resolved star formation law in NGC 7469 from JWST, ALMA and VLA

Authors: MariaVittoria Zanchettin, Marcella Massardi, Francesco Salvestrini, Manuela Bischetti, Chiara Feruglio, Fabrizio Fiore, Andrea Lapi

Abstract: We investigate the star formation process within the central 3.3 kpc region of the nearby luminous infrared Seyfert NGC 7469, probing scales ranging from 88 to 330 pc. We combine JWST/MIRI imaging with the F770W filter, with CO(2-1) and the underlying 1.3 mm dust continuum data from ALMA, along with VLA radio continuum observations at 22 GHz. NGC 7469 hosts a starburst ring which dominates the ove… ▽ More We investigate the star formation process within the central 3.3 kpc region of the nearby luminous infrared Seyfert NGC 7469, probing scales ranging from 88 to 330 pc. We combine JWST/MIRI imaging with the F770W filter, with CO(2-1) and the underlying 1.3 mm dust continuum data from ALMA, along with VLA radio continuum observations at 22 GHz. NGC 7469 hosts a starburst ring which dominates the overall star formation activity. We estimate a global star formation rate SFR $\sim 11.5$ $\rm M_{\odot}~yr^{-1}$ from the radio at 22 GHz, and a cold molecular gas mass M(H2) $\sim$ 6.4 $\times$ $\rm 10^9 M_{\odot}$ from the CO(2-1) emission. We find that the 1.3 mm map shows a morphology remarkably similar to those traced by the 22 GHz and the 7.7 $\rm μm$ polycyclic aromatic hydrocarbon (PAH) emission observed with JWST. The three tracers reproduce the morphology of the starburst ring with good agreement. We further investigate the correlations between the PAHs, the star formation rate and the cold molecular gas. We find a stronger correlation of the PAHs with the star formation than with the CO, with steeper correlations within the starburst ring ($n > 2$) than in the outer region ($n < 1$). We derive the correlation between the star formation rate and the cold molecular gas mass surface densities, the Kennicutt-Schmidt star formation law. Comparisons with other galaxy populations, including starburst galaxies and active galactic nuclei, highlighted that NGC 7469 exhibits an intermediate behavior to the Kennicutt-Schmidt relations found for these galaxy populations. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ

Multi-phase black-hole feedback and a bright [CII] halo in a Lo-BAL quasar at $z\sim6.6$

Authors: Manuela Bischetti, Hyunseop Choi, Fabrizio Fiore, Chiara Feruglio, Stefano Carniani, Valentina D'Odorico, Eduardo Bañados, Huanqing Chen, Roberto Decarli, Simona Gallerani, Julie Hlavacek-Larrondo, Samuel Lai, Karen M. Leighly, Chiara Mazzucchelli, Laurence Perreault-Levasseur, Roberta Tripodi, Fabian Walter, Feige Wang, Jinyi Yang, Maria Vittoria Zanchettin, Yongda Zhu

Abstract: Although the mass growth of supermassive black holes during the Epoch of Reionisation is expected to play a role in shaping the concurrent growth of their host-galaxies, observational evidence of feedback at z$\gtrsim$6 is still sparse. We perform the first multi-scale and multi-phase characterisation of black-hole driven outflows in the $z\sim6.6$ quasar J0923+0402 and assess how these winds impa… ▽ More Although the mass growth of supermassive black holes during the Epoch of Reionisation is expected to play a role in shaping the concurrent growth of their host-galaxies, observational evidence of feedback at z$\gtrsim$6 is still sparse. We perform the first multi-scale and multi-phase characterisation of black-hole driven outflows in the $z\sim6.6$ quasar J0923+0402 and assess how these winds impact the cold gas reservoir. We employ the SimBAL spectral synthesis to fit broad absorption line (BAL) features and find a powerful ionized outflow on $\lesssim210$ pc scale, with a kinetic power $\sim2-100$\% of the quasar luminosity. ALMA observations of [CII] emission allow us to study the morphology and kinematics of the cold gas. We detect high-velocity [CII] emission, likely associated with a cold neutral outflow at $\sim0.5-2$ kpc scale in the host-galaxy, and a bright extended [CII] halo with a size of $\sim15$ kpc. For the first time at such an early epoch, we accurately constrain the outflow energetics in both the ionized and the atomic neutral gas phases. We find such energetics to be consistent with expectations for an efficient feedback mechanism, and both ejective and preventative feedback modes are likely at play. The scales and energetics of the ionized and atomic outflows suggest that they might be associated with different quasar accretion episodes. The results of this work indicate that strong black hole feedback is occurring in quasars at $z\gtrsim6$ and is likely responsible for shaping the properties of the cold gas reservoir up to circum-galactic scales. △ Less

Submitted 16 May, 2024; v1 submitted 18 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ

MUSE view of PDS 456: kpc-scale wind, extended ionized gas and close environment

Authors: A. Travascio, E. Piconcelli, M. Bischetti, G. Cresci, C. Feruglio, M. Perna, G. Vietri, S. Carniani, S. Cantalupo, C. Cicone, M. Ginolfi, G. Venturi, K. Zubovas, A. Bongiorno, M. Brusa, A. Luminari, V. Mainieri, A. Marconi, N. Menci, E. Nardini, A. Pensabene, C. Ramos Almeida, F. Tombesi, C. Vignali, L. Zappacosta , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 15 pages, 13 figures accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A250 (2024)

HYPERION. Coevolution of supermassive black holes and galaxies at $z>6$ and the build-up of massive galaxies

Authors: R. Tripodi, C. Feruglio, F. Fiore, L. Zappacosta, E. Piconcelli, M. Bischetti, A. Bongiorno, S. Carniani, F. Civano, C. -C. Chen, S. Cristiani, G. Cupani, F. Di Mascia, V. D'Odorico, X. Fan, A. Ferrara, S. Gallerani, M. Ginolfi, R. Maiolino, V. Mainieri, A. Marconi, I. Saccheo, F. Salvestrini, A. Tortosa, R. Valiante

Abstract: 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… ▽ More 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 ($τ_{\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). △ Less

Submitted 28 June, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 29 pages; 6 tables; 18 figures. Accepted by A&A. A section about SF efficiency has been added compared to the previous version

A nebular origin for the persistent radio emission of fast radio bursts

Authors: Gabriele Bruni, Luigi Piro, Yuan-Pei Yang, Salvatore Quai, Bing Zhang, Eliana Palazzi, Luciano Nicastro, Chiara Feruglio, Roberta Tripodi, Brendan O'Connor, Angela Gardini, Sandra Savaglio, Andrea Rossi, A. M. Nicuesa Guelbenzu, Rosita Paladino

Abstract: Fast radio bursts (FRBs) are millisecond-duration, bright ($\sim$Jy) extragalactic bursts, whose production mechanism is still unclear. Recently, two repeating FRBs were found to have a physically associated persistent radio source of non-thermal origin. These two FRBs have unusually large Faraday rotation measure values likely tracing a dense magneto-ionic medium, consistent with synchrotron radi… ▽ More Fast radio bursts (FRBs) are millisecond-duration, bright ($\sim$Jy) extragalactic bursts, whose production mechanism is still unclear. Recently, two repeating FRBs were found to have a physically associated persistent radio source of non-thermal origin. These two FRBs have unusually large Faraday rotation measure values likely tracing a dense magneto-ionic medium, consistent with synchrotron radiation originating from a nebula surrounding the FRB source. Recent theoretical arguments predict that, if the observed Faraday rotation measure mostly arises from the persistent radio source region, there should be a simple relation between the luminosity of the latter and the first. We report here the detection of a third, less luminous persistent radio source associated with the repeating FRB source FRB20201124A at a distance of 413 Mpc, significantly expanding the predicted relation into the low luminosity - low Faraday rotation measure regime ($<$1000 rad m-2). At lower values of the Faraday rotation measure, the expected radio luminosity falls below the limit of detection threshold for present-day radio telescopes. These findings support the idea that the persistent radio sources observed so far are generated by a nebula in the FRB environment, and that FRBs with low Faraday rotation measure may not show a persistent radio source because of a weaker magneto-ionic medium. This is generally consistent with models invoking a young magnetar as the central engine of the FRB, where the surrounding ionized nebula - or the interacting shock in a binary system - powers the persistent radio source. △ Less

Submitted 10 June, 2024; v1 submitted 23 December, 2023; originally announced December 2023.

Comments: Accepted for publication on Nature

Boost recall in QSO selection from highly imbalanced photometric datasets

Authors: Giorgio Calderone, Francesco Guarneri, Matteo Porru, Stefano Cristiani, Andrea Grazian, Luciano Nicastro, Manuela Bischetti, Konstantina Boutsia, Guido Cupani, Valentina D'Odorico, Chiara Feruglio, Fabio Fontanot

Abstract: Context. The identification of bright QSOs is of great importance to probe the intergalactic medium and address open questions in cosmology. Several approaches have been adopted to find such sources in currently available photometric surveys, including machine learning methods. However, the rarity of bright QSOs at high redshifts compared to contaminating sources (such as stars and galaxies) makes… ▽ More Context. The identification of bright QSOs is of great importance to probe the intergalactic medium and address open questions in cosmology. Several approaches have been adopted to find such sources in currently available photometric surveys, including machine learning methods. However, the rarity of bright QSOs at high redshifts compared to contaminating sources (such as stars and galaxies) makes the selection of reliable candidates a difficult task, especially when high completeness is required. Aims. We present a novel technique to boost recall (i.e., completeness within the considered sample) in the selection of QSOs from photometric datasets dominated by stars, galaxies, and low-z QSOs (imbalanced datasets). Methods. Our method operates by iteratively removing sources whose probability of belonging to a noninteresting class exceeds a user-defined threshold, until the remaining dataset contains mainly high-z QSOs. Any existing machine learning method can be used as underlying classifier, provided it allows for a classification probability to be estimated. We applied the method to a dataset obtained by cross-matching PanSTARRS1, Gaia, and WISE, and identified the high-z QSO candidates using both our method and its direct multi-label counterpart. Results. We ran several tests by randomly choosing the training and test datasets, and achieved significant improvements in recall which increased from 50% to 85% for QSOs with z>2.5, and from 70% to 90% for QSOs with z>3. Also, we identified a sample of 3098 new QSO candidates on a sample of 2.6x10^6 sources with no known classification. We obtained follow-up spectroscopy for 121 candidates, confirming 107 new QSOs with z>2.5. Finally, a comparison of our candidates with those selected by an independent method shows that the two samples overlap by more than 90% and that both methods are capable of achieving a high level of completeness. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: Accepted for publication on A&A

Galaxy Formation and Symbiotic Evolution with the Inter-Galactic Medium in the Age of ELT-ANDES

Authors: Valentina D'Odorico, James S. Bolton, Lise Christensen, Annalisa De Cia, Erik Zackrisson, Aron Kordt, Luca Izzo, Jiangtao Li, Roberto Maiolino, Alessandro Marconi, Philipp Richter, Andrea Saccardi, Stefania Salvadori, Irene Vanni, Chiara Feruglio, Michele Fumagalli, Johan P. U. Fynbo, Pasquier Noterdaeme, Polychronis Papaderos, Celine Peroux, Aprajita Verma, Paolo Di Marcantonio, Livia Origlia, Alessio Zanutta

Abstract: High-resolution absorption spectroscopy toward bright background sources has had a paramount role in understanding early galaxy formation, the evolution of the intergalactic medium and the reionisation of the Universe. However, these studies are now approaching the boundaries of what can be achieved at ground-based 8-10m class telescopes. The identification of primeval systems at the highest redsh… ▽ More High-resolution absorption spectroscopy toward bright background sources has had a paramount role in understanding early galaxy formation, the evolution of the intergalactic medium and the reionisation of the Universe. However, these studies are now approaching the boundaries of what can be achieved at ground-based 8-10m class telescopes. The identification of primeval systems at the highest redshifts, within the reionisation epoch and even into the dark ages, and of the products of the first generation of stars and the chemical enrichment of the early Universe, requires observing very faint targets with a signal-to-noise ratio high enough to detect very faint spectral signatures. In this paper, we describe the giant leap forward that will be enabled by ANDES, the high-resolution spectrograph for the ELT, in these key science fields, together with a brief, non-exhaustive overview of other extragalactic research topics that will be pursued by this instrument, and its synergistic use with other facilities that will become available in the early 2030s. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 40 pages, 7 figures; submitted to Experimental Astronomy on behalf of the ANDES Science Team

NGC 2992: The interplay between the multiphase disk, wind and radio bubbles

Authors: Maria Vittoria Zanchettin, Chiara Feruglio, Marcella Massardi, Andrea Lapi, Manuela Bischetti, Sebastiano Cantalupo, Fabrizio Fiore, Angela Bongiorno, Angela Malizia, Andrea Marinucci, Manuela Molina, Enrico Piconcelli, Francesco Tombesi, Andrea Travascio, Giulia Tozzi, Roberta Tripodi

Abstract: We present an analysis of the gas kinematics in NGC 2992, based on VLT/MUSE, ALMA and VLA data, aimed at characterising the disk, the wind and their interplay in the cold molecular and warm ionised phases. CO(2-1) and H$\rm α~$ arise from a multiphase disk with inclination 80 deg and radii 1.5 and 1.8 kpc, respectively. We find that the velocity dispersion of the cold molecular phase is consistent… ▽ More We present an analysis of the gas kinematics in NGC 2992, based on VLT/MUSE, ALMA and VLA data, aimed at characterising the disk, the wind and their interplay in the cold molecular and warm ionised phases. CO(2-1) and H$\rm α~$ arise from a multiphase disk with inclination 80 deg and radii 1.5 and 1.8 kpc, respectively. We find that the velocity dispersion of the cold molecular phase is consistent with that of star forming galaxies at the same redshift, except in the inner 600 pc region, and in the region between the cone walls and the disk. This suggests that a disk-wind interaction locally boosts the gas turbulence. We detect a clumpy ionised wind distributed in two wide opening angle ionisation cones reaching scales of 7 kpc. The [O III] wind expands with velocity exceeding -1000 km/s in the inner 600 pc, a factor of 5 larger than the previously reported wind velocity. Based on spatially resolved electron density and ionisation parameter maps, we infer an ionised outflow mass of $M_{\rm of,ion} = (3.2 \pm 0.3) \times \, 10^7 \, M_{\odot}$, and a total ionised outflow rate of $\dot M_{\rm of,ion}=13.5\pm1$ \sfr. We detected clumps of cold molecular gas located above and below the disk reaching maximum projected distances and velocities of 1.7 kpc and 200 km/s, respectively. On these scales, the wind is multiphase, with a fast ionised component and a slower molecular one, and a total mass of $M_{\rm of, ion+mol}= 5.8 \times 10^7 \, M_{\odot}$, of which the molecular component carries the bulk of the mass. The dusty molecular outflowing clumps and the turbulent ionised gas are located at the edges of the radio bubbles, suggesting that the bubbles interact with the surrounding medium through shocks. We detect a dust reservoir co-spatial with the molecular disk, with a cold dust mass $M_{\rm dust} = (4.04 \pm 0.03) \times \, 10^{6} \, M_{\odot}$. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: 19 pages, 17 figures, 6 tables; Accepted by A&A

Journal ref: A&A 679, A88 (2023)

Crossing the Rubicon of Reionization with z~5 QSOs

Authors: A. Grazian, K. Boutsia, E. Giallongo, S. Cristiani, F. Fontanot, M. Bischetti, A. Bongiorno, G. Calderone, G. Cupani, V. D'Odorico, C. Feruglio, F. Fiore, F. Guarneri, M. Porru, I. Saccheo

Abstract: One of the key open questions in Cosmology is the nature of the sources that completed the cosmological hydrogen Reionization at z~5.2. High-z primeval galaxies have been long considered the main drivers for Reionization, with a minor role played by high-z AGN. However, in order to confirm this scenario, it is fundamental to measure the photo-ionization rate produced by active SMBHs close to the e… ▽ More One of the key open questions in Cosmology is the nature of the sources that completed the cosmological hydrogen Reionization at z~5.2. High-z primeval galaxies have been long considered the main drivers for Reionization, with a minor role played by high-z AGN. However, in order to confirm this scenario, it is fundamental to measure the photo-ionization rate produced by active SMBHs close to the epoch of Reionization. Given the pivotal role played by spectroscopically complete observations of high-z QSOs, in this paper we present the first results of the RUBICON (Reionizing the Universe with BrIght COsmological Nuclei) survey. It consists of a color selected sample of bona-fide z~5 QSO candidates from the Hyper Suprime-Cam Subaru Strategic Survey. Our QSO candidates have been validated both by photometric redshifts based on SED fitting and by spectroscopic redshifts, confirming that they lie at 4.5<z_spec<5.2. A relatively large space density of QSOs (Phi~1.4x10^-8 cMpc^-3) is thus confirmed at z~5 and M1450~-27, consistent with a pure density evolution of the AGN luminosity function from z=4 to z=5, with a mild density evolution rate of 0.25 dex. This indicates that AGN could play a non-negligible role in the cosmic Reionization. The Rubicon of Reionization has been crossed. △ Less

Submitted 3 August, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ, 25 pages, 11 figures, 3 tables, figures 1, 3, 5, and 6 resized

Star formation efficiency and AGN feedback in narrow-line Seyfert 1 galaxies with fast X-ray nuclear winds

Authors: Quentin Salomé, Yair Krongold, Anna Lia Longinotti, Manuela Bischetti, Santiago García-Burillo, Olga Vega, Miguel Sánchez-Portal, Chiara Feruglio, María Jesús Jiménez-Donaire, Maria Vittoria Zanchettin

Abstract: We present the first systematic study of the molecular gas and star formation efficiency in a sample of ten narrow-line Seyfert 1 galaxies selected to have X-ray Ultra Fast Outflows and, therefore, to potentially show AGN feedback effects. CO observations were obtained with the IRAM 30m telescope in six galaxies and from the literature for four galaxies. We derived the stellar mass, star formation… ▽ More We present the first systematic study of the molecular gas and star formation efficiency in a sample of ten narrow-line Seyfert 1 galaxies selected to have X-ray Ultra Fast Outflows and, therefore, to potentially show AGN feedback effects. CO observations were obtained with the IRAM 30m telescope in six galaxies and from the literature for four galaxies. We derived the stellar mass, star formation rate, AGN and FIR dust luminosities by fitting the multi-band spectral energy distributions with the CIGALE code. Most of the galaxies in our sample lie above the main sequence (MS) and the molecular depletion time is one to two orders of magnitude shorter than the one typically measured in local star-forming galaxies. Moreover, we found a promising correlation between the star formation efficiency and the Eddington ratio, as well as a tentative correlation with the AGN luminosity. The role played by the AGN activity in the regulation of star formation within the host galaxies of our sample remains uncertain (little or no effect? positive feedback?). Nevertheless, we can conclude that quenching by the AGN activity is minor and that star formation will likely stop in a short time due to gas exhaustion by the current starburst episode. △ Less

Submitted 24 July, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: Published in MNRAS, Volume 524, Issue 2, Pages 3130-3145

Radio WISSH: tuning on the most luminous quasars in the Universe

Authors: Gabriele Bruni, Javier Moldón, Enrico Piconcelli, Francesca Panessa, Miguel Pérez-Torres, Manuela Bischetti, Chiara Feruglio, Giustina Vietri, Cristian Vignali, Luca Zappacosta, Ivano Saccheo

Abstract: In the past years, the results obtained by the WISSH quasar project provided a novel general picture on the distinctive multi-band properties of hyper-luminous ($L_{bol}>10^{47}$ erg/s) quasars at high redshift (z$\sim$2-4), unveiling interesting relations among active galactic nuclei, winds and interstellar medium, in these powerful sources at cosmic noon. Since 2022, we are performing a systemat… ▽ More In the past years, the results obtained by the WISSH quasar project provided a novel general picture on the distinctive multi-band properties of hyper-luminous ($L_{bol}>10^{47}$ erg/s) quasars at high redshift (z$\sim$2-4), unveiling interesting relations among active galactic nuclei, winds and interstellar medium, in these powerful sources at cosmic noon. Since 2022, we are performing a systematic and statistically-significant VLA study of the radio properties of WISSH. We carried out high-resolution VLA observations aiming at: 1) identifying young radio source from the broad-band spectral shape of these objects; 2) sample an unexplored high redshift/high luminosity regime, tracking possible evolutionary effects on the radio-loud/radio-quiet dichotomy; 3) quantifying orientation effects on the observed winds/outflows properties. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: To be published on "Proceedings of the International Astronomical Union", IAU symposium #378: "Black hole winds at all scales"

XQR-30: Black Hole Masses and Accretion Rates of 42 z>6 Quasars

Authors: C. Mazzucchelli, M. Bischetti, V. D'Odorico, C. Feruglio, J. -T. Schindler, M. Onoue, E. Bañados, G. D. Becker, F. Bian, S. Carniani, R. Decarli, A. -C. Eilers, E. P. Farina, S. Gallerani, S. Lai, R. A. Meyer, S. Rojas-Ruiz, S. Satyavolu, B. P. Venemans, F. Wang, J. Yang, Y. Zhu

Abstract: We present bolometric luminosities, black hole masses and Eddington ratios for 42 luminous quasars at z>6 using high signal-to-noise ratio VLT/X-Shooter spectra, acquired in the enlarged ESO Large Programme XQR-30. In particular, we derive bolometric luminosities from the rest-frame 3000 A, luminosities using a bolometric correction from the literature, and the black hole masses by modelling the s… ▽ More We present bolometric luminosities, black hole masses and Eddington ratios for 42 luminous quasars at z>6 using high signal-to-noise ratio VLT/X-Shooter spectra, acquired in the enlarged ESO Large Programme XQR-30. In particular, we derive bolometric luminosities from the rest-frame 3000 A, luminosities using a bolometric correction from the literature, and the black hole masses by modelling the spectral regions around the CIV 1549A and the MgII 2798A emission lines, with scaling relations calibrated in the local universe. We find that the black hole masses derived from both emission lines are in the same range, and the scatter of the measurements agrees with expectations from the scaling relations. The MgII-derived masses are between ~(0.8-12) x 10^9 Msun, and the derived Eddington ratios are within ~0.13-1.73, with a mean (median) of 0.84 (0.72). By comparing the total sample of quasars at z>5.8, from this work and from the literature, to a bolometric luminosity distribution-matched sample at z~1.5, we find that quasars at high redshift host slightly less massive black holes which accrete slightly more rapidly than at lower-z, with a difference in the mean Eddington ratios of the two samples of ~0.27, in agreement with recent literature work. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 9 pages; 5 figures; accepted for publication in A&A

Journal ref: A&A 676, A71 (2023)

HYPERION. Interacting companion and outflow in the most luminous $z>6$ quasar

Authors: R. Tripodi, J. Scholtz, R. Maiolino, S. Fujimoto, S. Carniani, J. D. Silverman, C. Feruglio, M. Ginolfi, L. Zappacosta, T. Costa, G. C. Jones, E. Piconcelli, M. Bischetti, F. Fiore

Abstract: We present ALMA deep observations of the [CII] 158 $μ$m emission line and the continuum at 253 GHz and 99 GHz towards SDSS J0100+2802 at $z\simeq 6.3$, the most luminous quasi-stellar object (QSO) at z$>$6. J0100+2802 belongs to the HYPERION sample of luminous QSOs at $z\sim 6-7.5$. The observations (at 2.2" resolution in Band 3 and 0.9" resolution in Band 6) are optimized to detect extended emiss… ▽ More We present ALMA deep observations of the [CII] 158 $μ$m emission line and the continuum at 253 GHz and 99 GHz towards SDSS J0100+2802 at $z\simeq 6.3$, the most luminous quasi-stellar object (QSO) at z$>$6. J0100+2802 belongs to the HYPERION sample of luminous QSOs at $z\sim 6-7.5$. The observations (at 2.2" resolution in Band 3 and 0.9" resolution in Band 6) are optimized to detect extended emission around the QSO. We detect an interacting, tidally disrupted companion both in [CII], peaking at $z\sim 6.332$, and in continuum, stretching on scales up to 20 kpc from the quasar, with a knotty morphology. The higher velocity dispersion in the direction of the companion emission and the complex morphology of tidally stretched galaxy suggest a possible ongoing or future merger. For the newly-detected companion we derive a range of dust mass, $M_{\rm dust}=(0.3-2.6)\times 10^7\ \rm M_\odot$, and of star formation rate, SFR$=[35-344]\ \rm M_\odot$. This shows that both the QSO and its companion are gas rich and that a major merger may be at the origin of the boosted star formation. We also detect a broad blueshifted component in the [CII] spectrum, that we interpret as a gaseous outflow for which we estimate a mass outflow rate in the range $\dot{M}_{\rm out}=(118-269)\ \rm M_\odot\ yr^{-1}$. J0100+2802 was recently found to reside in a strong overdensity, however this close companion remained undetected by both previous higher resolution ALMA observations and by JWST-NIRCAM imaging. Our results highlight the importance of deep medium-resolution ALMA observations for the study of QSOs and their environment at the Epoch of Reionization. △ Less

Submitted 9 January, 2024; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 11 pages, 7 figures, 3 Tables. Accepted by A&A. In press

XQR-30: the ultimate XSHOOTER quasar sample at the reionization epoch

Authors: Valentina D'Odorico, E. Banados, G. D. Becker, M. Bischetti, S. E. I. Bosman, G. Cupani, R. Davies, E. P. Farina, A. Ferrara, C. Feruglio, C. Mazzucchelli, E. Ryan-Weber, J. -T. Schindler, A. Sodini, B. P. Venemans, F. Walter, H. Chen, S. Lai, Y. Zhu, F. Bian, S. Campo, S. Carniani, S. Cristiani, F. Davies, R. Decarli , et al. (24 additional authors not shown)

Abstract: 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-quali… ▽ More 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 $λ\simeq 1285$ Å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. △ Less

Submitted 19 June, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 21 pages, 10 figures. Final version accepted by MNRAS

Journal ref: MNRAS, Volume 523, Issue 1, pp.1399-1420 (2023)

The lively accretion disc in NGC 2992. III. Tentative evidence of rapid Ultra Fast Outflow variability

Authors: Alfredo Luminari, Andrea Marinucci, Stefano Bianchi, Barbara de Marco, Chiara Feruglio, Giorgio Matt, Riccardo Middei, Emanuele Nardini, Enrico Piconcelli, Simonetta Puccetti, Francesco Tombesi

Abstract: We report on the 2019 XMM-Newton+NuSTAR monitoring campaign of the Seyfert galaxy NGC 2992, observed at one of its highest flux levels in the X-rays. The time-averaged spectra of the two XMM-Newton orbits show Ultra Fast Outflows (UFOs) absorbing structures above 9 keV with $> 3 σ$ significance. A detailed investigation of the temporal evolution on a $\sim$5 ks time scale reveals UFO absorption li… ▽ More We report on the 2019 XMM-Newton+NuSTAR monitoring campaign of the Seyfert galaxy NGC 2992, observed at one of its highest flux levels in the X-rays. The time-averaged spectra of the two XMM-Newton orbits show Ultra Fast Outflows (UFOs) absorbing structures above 9 keV with $> 3 σ$ significance. A detailed investigation of the temporal evolution on a $\sim$5 ks time scale reveals UFO absorption lines at a confidence level $>$95% (2$σ$) in 8 out of 50 XMM-Newton segments, estimated via Monte Carlo simulations. We observe a wind variability corresponding to a length scale of 5 Schwarzschild radii $r_S$. Adopting the novel Wind in the Ionised Nuclear Environment (WINE) model, we estimate the outflowing gas velocity ($v=0.21-0.45 c$), column density ($N_H=4-8\cdot 10^{24} cm^{-2}$) and ionisation state ($\log(ξ_0/erg\ cm\ s^{-1})=3.7-4.7$), taking into account geometrical and special relativity corrections. These parameters lead to instantaneous mass outflow rates $\dot{M}_{out}\simeq0.3-0.8 M_{\odot} yr^{-1}$, with associated outflow momentum rates $\dot{p}_{out}\simeq 20-90 L_{Bol}/c$ and kinetic energy rates $\dot{E}_K \simeq 2-25 L_{Bol}$. We estimate a wind duty cycle $\approx$ 12% and a total mechanical power $\approx$ 2 times the AGN bolometric luminosity, suggesting the wind may drive significant feedback effects between the AGN and the host galaxy. Notably, we also provide an estimate for the wind launching radius and density $\approx 5 r_S, 10^{11} {cm}^{-3}$, respectively. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Astrophysical Journal. 21 pages, 11 figures

HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). A new regime for the X-ray nuclear properties of the first quasars

Authors: L. Zappacosta, E. Piconcelli, F. Fiore, I. Saccheo, R. Valiante, C. Vignali, F. Vito, M. Volonteri, M. Bischetti, A. Comastri, C. Done, M. Elvis, E. Giallongo, F. La Franca, G. Lanzuisi, M. Laurenti, G. Miniutti, A. Bongiorno, M. Brusa, F. Civano, S. Carniani, V. D'Odorico, C. Feruglio, S. Gallerani, R. Gilli , et al. (18 additional authors not shown)

Abstract: The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous Q… ▽ More The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($Γ\sim2.4\pm0.1$). Absorption is not required. The average $Γ$ is inconsistent at $\geq4σ$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $Γ=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe. △ Less

Submitted 19 July, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

Journal ref: A&A 678, A201 (2023)

First constraints of dense molecular gas at z~7.5 from the quasar Pōniuā'ena

Authors: Chiara Feruglio, Umberto Maio, Roberta Tripodi, Jan Martin Winters, Luca Zappacosta, Manuela Bischetti, Francesca Civano, Stefano Carniani, Valentina D'Odorico, Fabrizio Fiore, Simona Gallerani, Michele Ginolfi, Roberto Maiolino, Enrico Piconcelli, Rosa Valiante, Maria Vittoria Zanchettin

Abstract: We report the detection of CO(6-5) and CO(7-6) and their underlying continua from the host galaxy of quasar J100758.264+211529.207 (Pōniuā'ena) at z=7.5419, obtained with the NOrthern Extended Millimeter Array (NOEMA). Pōniuā'ena belongs to the HYPerluminous quasars at the Epoch of ReionizatION (HYPERION) sample of 17 $z>6$ quasars selected to be powered by supermassive black holes (SMBH) which ex… ▽ More We report the detection of CO(6-5) and CO(7-6) and their underlying continua from the host galaxy of quasar J100758.264+211529.207 (Pōniuā'ena) at z=7.5419, obtained with the NOrthern Extended Millimeter Array (NOEMA). Pōniuā'ena belongs to the HYPerluminous quasars at the Epoch of ReionizatION (HYPERION) sample of 17 $z>6$ quasars selected to be powered by supermassive black holes (SMBH) which experienced the fastest mass growth in the first Gyr of the Universe. The one reported here is the highest-redshift measurement of the cold and dense molecular gas to date. The host galaxy is unresolved and the line luminosity implies a molecular reservoir of $\rm M(H_2)=(2.2\pm0.2)\times 10^{10}$ $\rm M_\odot$, assuming a CO spectral line energy distribution typical of high-redshift quasars and a conversion factor $α=0.8$ $\rm M_{\odot} (K\,km \, s^{-1} \,pc^{2})^{-1} $. We model the cold dust spectral energy distribution (SED) to derive a dust mass of M$_{\rm dust} =(2.1\pm 0.7)\times 10^8$ $\rm M_\odot$, and thus a gas to dust ratio $\sim100$. Both the gas and dust mass are not dissimilar from the reservoir found for luminous quasars at $z\sim6$. We use the CO detection to derive an estimate of the cosmic mass density of $\rm H_2$, $Ω_{H_2} \simeq 1.31 \times 10^{-5}$. This value is in line with the general trend suggested by literature estimates at $ z < 7 $ and agrees fairly well with the latest theoretical expectations of non-equilibrium molecular-chemistry cosmological simulations of cold gas at early times. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ Letters

Outflows in the Gaseous Discs of Active Galaxies and their impact on Black Hole Scaling Relations

Authors: N. Menci, F. Fiore, F. Shankar, L. Zanisi, C. Feruglio

Abstract: To tackle the still unsolved and fundamental problem of the role of Active Galactic Nuclei (AGN) feedback in shaping galaxies, in this work we implement a new physical treatment of AGN-driven winds into our semi-analytic model of galaxy formation. To each galaxy in our model, we associate solutions for the outflow expansion and the mass outflow rates in different directions, depending on the AGN l… ▽ More To tackle the still unsolved and fundamental problem of the role of Active Galactic Nuclei (AGN) feedback in shaping galaxies, in this work we implement a new physical treatment of AGN-driven winds into our semi-analytic model of galaxy formation. To each galaxy in our model, we associate solutions for the outflow expansion and the mass outflow rates in different directions, depending on the AGN luminosity, on the circular velocity of the host halo, and on gas content of the considered galaxy. To each galaxy we also assign an effective radius derived from energy conservation during merger events, and a stellar velocity dispersion self-consistently computed via Jeans modelling. We derive all the main scaling relations between Black hole (BH) mass and total/bulge stellar mass, velocity dispersion, host halo dark matter mass, and star formation efficiency. We find that our improved AGN feedback mostly controls the dispersion around the relations but plays a subdominant role in shaping slopes and/or normalizations of the scaling relations. Including possible limited-resolution selection biases in the model provides better agreement with the available data. The model does not point to any more fundamental galactic property linked to BH mass, with velocity dispersion playing a similar role with respect to stellar mass, in tension with present data. In line with other independent studies carried out on comprehensive semi-analytic and hydrodynamic galaxy-BH evolution models, our current results signal either an inadequacy of present cosmological models of galaxy formation in fully reproducing the local scaling relations, in terms of both shape and residuals, and/or point to an incompleteness issue affecting the local sample of dynamically-measured BHs. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 21 pages, accepted for publications in Astronomy & Astrophysics

Accurate dust temperature and star formation rate in the most luminous $z>6$ quasar in the HYPerluminous quasars at the Epoch of ReionizatION (HYPERION) sample

Authors: Roberta Tripodi, Chiara Feruglio, Francisca Kemper, Francesca Civano, Tiago Costa, Martin Elvis, Manuela Bischetti, Stefano Carniani, Fabio Di Mascia, Valentina D'Odorico, Fabrizio Fiore, Simona Gallerani, Michele Ginolfi, Roberto Maiolino, Enrico Piconcelli, Rosa Valiante, Luca Zappacosta

Abstract: We present ALMA Band 9 continuum observation of the ultraluminous quasi-stellar object (QSO) SDSS J0100+2802, providing a $\sim 10σ$ detection at $\sim 670$ GHz. SDSS J0100+2802 is the brightest QSO with the most massive super massive black hole (SMBH) known at $z>6$, and we study its dust spectral energy distribution in order to determine the dust properties and the star formation rate (SFR) of i… ▽ More We present ALMA Band 9 continuum observation of the ultraluminous quasi-stellar object (QSO) SDSS J0100+2802, providing a $\sim 10σ$ detection at $\sim 670$ GHz. SDSS J0100+2802 is the brightest QSO with the most massive super massive black hole (SMBH) known at $z>6$, and we study its dust spectral energy distribution in order to determine the dust properties and the star formation rate (SFR) of its host-galaxy. We obtain the most accurate estimate so far of the temperature, mass and emissivity index of the dust, having $T_{\rm dust}=48.4\pm2.3$ K, $M_{\rm dust}=(2.29\pm0.83)\times 10^7$ M$_\odot$, $β=2.63\pm 0.23$. This allows us to measure the SFR with the smallest statistical error for this QSO, SFR$=265\pm 32\ \rm M_\odot yr^{-1}$. Our results enable us to evaluate the relative growth of the SMBH and host galaxy of J0100+2802, finding that the SMBH is dominating the process of BH-galaxy growth in this QSO at $z=6.327$, when the Universe was $865$ Myr old. Such unprecedented constraints on the host galaxy SFR and dust temperature can only be obtained through high frequency observations, and highlight the importance of ALMA Band 9 to obtain a robust overview of the build-up of the first quasars' host galaxies at $z>6$. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 10 pages, 4 figures, 1 table. Accepted for publication in ApJL

NOEMA spatially resolved view of the multi-phase outflow in IRAS17020+4544: a shocked wind in action?

Authors: Anna Lia Longinotti, Quentin Salomé, Chiara Feruglio, Yair Krongold, Santiago García-Burillo, Marcello Giroletti, Francesca Panessa, Carlo Stanghellini, Olga Vega, Victor Manuel Patiño-Álvarez, Vahram Chavushyan, Mauricio Elías-Chavez, Aitor Robleto-Orús

Abstract: The Narrow Line Seyfert 1 Galaxy IRAS17020+4544 is one of the few AGN where a galaxy-scale energy-conserving outflow was revealed. This paper reports on NOEMA observations addressed to constrain the spatial scale of the CO emission in outflow. The molecular outflowing gas is resolved in five components tracing approaching and receding gas, all located at a distance of 2-3~kpc on the West and East… ▽ More The Narrow Line Seyfert 1 Galaxy IRAS17020+4544 is one of the few AGN where a galaxy-scale energy-conserving outflow was revealed. This paper reports on NOEMA observations addressed to constrain the spatial scale of the CO emission in outflow. The molecular outflowing gas is resolved in five components tracing approaching and receding gas, all located at a distance of 2-3~kpc on the West and East side of the active nucleus. This high velocity gas (up to v_out=~1900 km/s) is not coincident with the rotation pattern of the CO gas in the host galaxy disk. The estimated mass outflow rate shows that with a global mass output of $\dot{M}_{H_2}$=~139$\pm$20$~M_\odot$~yr$^{-1}$, this powerful galaxy-scale outflow is consistent with the wind conserving its energy, and with a momentum rate boost of a factor of ~30 compared to the momentum rate of the nuclear X-ray wind. Preliminary results from ancillary X-ray (Chandra) and radio images (e-MERLIN) are reported. While the nature of the radio source is not conclusive, the Chandra image may tentatively trace extended emission, as expected by an expanding bubble of hot X-ray gas. The outcome of the NOEMA analysis and of past and ongoing publications dedicated to the description of the outflow multi-band phenomenology in IRAS17020+4544 concur to provide compelling reasons to postulate that an outflow shocking with the galaxy interstellar medium is driving the multi-phase wind in this peculiar AGN. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: Accepted for publication on MNRAS, 17 pages, 7 figures

X-Ray Detection of the Galaxy's Missing Baryons in the Circum-Galactic Medium of L$^*$ Galaxies

Authors: Fabrizio Nicastro, Yair Krongold, Taotao Fang, Filippo Fraternali, Smita Mathur, Stefano Bianchi, Alessandra De Rosa, Enrico Piconcelli, Luca Zappacosta, Manuela Bischetti, Chiara Feruglio, Anjali Gupta, Zheng Zhou

Abstract: The amount of baryons hosted in the disks of galaxies is lower than expected based on the mass of their dark-matter halos and the fraction of baryon-to-total matter in the universe, giving rise to the so called galaxy missing-baryon problem. The presence of cool circum-galactic matter gravitationally bound to its galaxy's halo up to distances of at least ten times the size of the galaxy's disk, mi… ▽ More The amount of baryons hosted in the disks of galaxies is lower than expected based on the mass of their dark-matter halos and the fraction of baryon-to-total matter in the universe, giving rise to the so called galaxy missing-baryon problem. The presence of cool circum-galactic matter gravitationally bound to its galaxy's halo up to distances of at least ten times the size of the galaxy's disk, mitigates the problem but is far from being sufficient for its solution. It has instead been suggested, that the galaxy missing baryons may hide in a much hotter gaseous phase of the circum-galactic medium, possibly near the halo virial temperature and co-existing with the cool phase. Here we exploit the best available X-ray spectra of known cool circum-galactic absorbers of L$^*$ galaxies to report the first direct high-statistical-significance (best estimates ranging from $4.2-5.6σ$, depending on fitting methodology)} detection of associated OVII absorption in the stacked XMM and Chandra spectra of three quasars. We show that these absorbers trace hot medium in the X-ray halo of these systems, at logT(in k)$\simeq 5.8-6.3$ K (comprising the halo virial temperature T$_{vir} \simeq 10^6$ K). We estimate masses of the X-ray halo within 1 virial radius within the interval M$_{hot-CGM}\simeq (1-1.7)\times 10^{11} (Z/0.3 Z_{\odot})^{-1}$ M$_{\odot}$. For these systems, this corresponds to galaxy missing baryon fractions in the range $ξ_b = M_{hot-CGM}/M_{missing}\simeq (0.7-1.2) (Z/0.3 Z_{\odot})^{-1}$, thus potentially closing the galaxy baryon census in typical L$^*$ galaxies. Our measurements contribute significantly to the solution of the long-standing galaxy missing baryon problem and to the understanding of the continuous cycle of baryons in-and-out of galaxies throughout the life of the universe. △ Less

Submitted 1 August, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: 26 pages, 11 figures, accepted for publication in the ApJL

The fraction and kinematics of broad absorption line quasars across cosmic time

Authors: Manuela Bischetti, Fabrizio Fiore, Chiara Feruglio, Valentina D'Odorico, Nahum Arav, Tiago Costa, Kastytis Zubovas, George Becker, Sarah E. I. Bosman, Guido Cupani, Rebecca Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Andrea Ferrara, Massimo Gaspari, Chiara Mazzucchelli, Masafusa Onoue, Enrico Piconcelli, Maria-Vittoria Zanchettin, Yongda Zhu

Abstract: Luminous quasars are powerful targets to investigate the role of feedback from supermassive black-holes (BHs) in regulating the growth phases of BHs themselves and of their host galaxies, up to the highest redshifts. Here we investigate the cosmic evolution of the occurrence and kinematics of BH-driven outflows, as traced by broad absorption line (BAL) features, due to the C IV ionic transition. W… ▽ More Luminous quasars are powerful targets to investigate the role of feedback from supermassive black-holes (BHs) in regulating the growth phases of BHs themselves and of their host galaxies, up to the highest redshifts. Here we investigate the cosmic evolution of the occurrence and kinematics of BH-driven outflows, as traced by broad absorption line (BAL) features, due to the C IV ionic transition. We exploit a sample of 1935 quasars quasars at $z=2.1-6.6$ with bolometric luminosity log($L_{\rm bol}/$erg s$^{-1})\gtrsim46.5$, drawn from the Sloan Digital Sky Survey and from the X-shooter legacy survey of Quasars at Reionisation (XQR-30). We consider rest-frame optical bright quasars to minimise observational biases due to quasar selection criteria. We apply a homogeneous BAL identification analysis, based on employing composite template spectra to estimate the quasar intrinsic emission. We find a BAL quasar fraction close to 20\% at $z\sim2-4$, while it increases to almost 50\% at $z\sim6$. The velocity and width of the BAL features also increase at $z\gtrsim4.5$. We exclude that the redshift evolution of the BAL properties is due to differences in terms of quasar luminosity and accretion rate. These results suggest significant BH feedback occurring in the 1 Gyr old Universe, likely affecting the growth of BHs and, possibly, of their host galaxies, as supported by models of early BH and galaxy evolution. △ Less

Submitted 18 April, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: Accepted for publication in ApJ

Dynamical signature of a stellar bulge in a quasar host galaxy at $z\simeq 6$

Authors: Roberta Tripodi, Federico Lelli, Chiara Feruglio, Fabrizio Fiore, Fabio Fontanot, Manuela Bischetti, Roberto Maiolino

Abstract: We present a dynamical analysis of a quasar-host galaxy at $z\simeq 6$ (SDSS J2310+1855) using a high-resolution ALMA observation of the [CII] emission line. The observed rotation curve was fitted with mass models that considered the gravitational contribution of a thick gas disc, a thick star-forming stellar disc, and a central mass concentration, which is likely due to a combination of a spheroi… ▽ More We present a dynamical analysis of a quasar-host galaxy at $z\simeq 6$ (SDSS J2310+1855) using a high-resolution ALMA observation of the [CII] emission line. The observed rotation curve was fitted with mass models that considered the gravitational contribution of a thick gas disc, a thick star-forming stellar disc, and a central mass concentration, which is likely due to a combination of a spheroidal component (i.e. a stellar bulge) and a supermassive black hole (SMBH). The SMBH mass of $5\times 10^9\ \rm M_{\odot}$, previously measured using the CIV and MgII emission lines, is not sufficient to explain the high velocities in the central regions. Our dynamical model suggests the presence of a stellar bulge with a mass of $\rm M_{bulge}\sim 10^{10}\ \rm M_{\odot}$ in this object, when the Universe was younger than 1 Gyr. To finally be located on the local $M_{\rm SMBH}-M_{\rm bulge}$ relation, the bulge mass should increase by a factor of $\sim$40 from $z=6$ to 0, while the SMBH mass should grow by a factor of 4 at most. This points towards asynchronous galaxy-BH co-evolution. Imaging with the JWST will allow us to validate this scenario. △ Less

Submitted 23 January, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 8 pages, 4 figures, 3 tables. Accepted by A&A

Journal ref: A&A 671, A44 (2023)

Dusty winds clear JWST super-early galaxies

Authors: Fabrizio Fiore, Andrea Ferrara, Manuela Bischetti, Chiara Feruglio, Andrea Travascio

Abstract: The JWST discovery of a number of super-early (redshift $z>10$), blue galaxies requires these systems to be essentially dust-free in spite of their large stellar masses. A possible explanation is that dust is evacuated by radiatively-driven outflows. We test this hypothesis by deriving the Eddington ratio $λ_E=L_{\rm bol}/L_{E}$, where $L_{\rm bol}$ is the bolometric luminosity produced by star-fo… ▽ More The JWST discovery of a number of super-early (redshift $z>10$), blue galaxies requires these systems to be essentially dust-free in spite of their large stellar masses. A possible explanation is that dust is evacuated by radiatively-driven outflows. We test this hypothesis by deriving the Eddington ratio $λ_E=L_{\rm bol}/L_{E}$, where $L_{\rm bol}$ is the bolometric luminosity produced by star-formation and possible black hole accretion, for 134 galaxies at $6.5< z <16$. We find a strong anti-correlation between $λ_E$ and dust UV optical depth, $τ_{1500} \propto λ_E^{-0.63}$; also, $λ_E$ increases with redshift. We confirm that galaxies exceeding a specific star formation rate ${\rm sSFR} > 13\, \rm Gyr^{-1}$ develop powerful outflows clearing the galaxy from its dust. This result is supported by ALMA dust continuum non-detections in three super-early systems. △ Less

Submitted 25 January, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: ApJ Letters in press

The WISSH quasars project XI. The mean Spectral Energy Distribution and Bolometric Corrections of the most luminous quasars

Authors: Ivano Saccheo, Angela Bongiorno, Enrico Piconcelli, Manuela Bischetti, Gabriele Bruni, Giovanni Cresci, Chiara Feruglio, Fabrizio Fiore, Andrea Grazian, Alfredo Luminari, Elisabeta Lusso, Vincenzo Mainieri, Roberto Maiolino, Alessandro Marconi, Federica Ricci, Francesco Tombesi, Andrea Travascio, Giustina Vietri, Cristian Vignali, Luca Zappacosta, Fabio La Franca

Abstract: Hyper-luminous Quasi-Stellar Objects (QSOs) represent the ideal laboratory to investigate Active Galactic Nuclei (AGN) feedback mechanism since their formidable energy release causes powerful winds at all scales and thus the maximum feedback is expected. We aim at deriving the mean Spectral Energy Distribution (SED) of a sample of 85 WISE-SDSS Selected Hyper-luminous (WISSH) quasars. Since the S… ▽ More Hyper-luminous Quasi-Stellar Objects (QSOs) represent the ideal laboratory to investigate Active Galactic Nuclei (AGN) feedback mechanism since their formidable energy release causes powerful winds at all scales and thus the maximum feedback is expected. We aim at deriving the mean Spectral Energy Distribution (SED) of a sample of 85 WISE-SDSS Selected Hyper-luminous (WISSH) quasars. Since the SED provides a direct way to investigate the AGN structure, our goal is to understand if quasars at the bright end of the luminosity function have peculiar properties compared to the bulk of the population. We built a mean intrinsic SED after correcting for the dust extinction, absorption and emission lines and intergalactic medium absorption. We also derived bolometric, IR band and monochromatic luminosities together with bolometric corrections at lambda = 5100 A and 3 micron. We define a new relation for the 3 micron bolometric correction. We find that the mean SED of hyper-luminous WISSH QSOs is different from that of less luminous sources, i.e. a relatively lower X-ray emission and a near and mid IR excess which can be explained assuming a larger dust contribution. WISSH QSOs have stronger emission from both warm and very hot dust, the latter being responsible for shifting the typical dip of the AGN SED from 1.3 to 1.1 micron. We also derived the mean SEDs of two sub-samples created according to the presence of Broad Absorption Lines and equivalent width of CIV line. We confirm that BALs are X-ray weak and that they have a reddened UV-optical continuum. We also find that BALs tend to have stronger emission from the hot dust component. This analysis suggests that hyper-luminous QSOs have a peculiar SED compared to less luminous objects. It is therefore critical to use SED templates constructed exclusively from very bright quasars samples when dealing with particularly luminous sources. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A. 20 pages, 15 figures

Black hole-galaxy co-evolution and the role of feedback

Authors: Pedro R. Capelo, Chiara Feruglio, Ryan C. Hickox, Francesco Tombesi

Abstract: Active galactic nuclei (AGN) are accreting supermassive black holes co-evolving with their host galaxies through a complex interplay of feeding and feedback. In this chapter, we first discuss AGN fuelling in galaxies, both in interacting and isolated systems, focusing on the role that instabilities have on the angular momentum budget of the gas. We then review observations and models of feedback t… ▽ More Active galactic nuclei (AGN) are accreting supermassive black holes co-evolving with their host galaxies through a complex interplay of feeding and feedback. In this chapter, we first discuss AGN fuelling in galaxies, both in interacting and isolated systems, focusing on the role that instabilities have on the angular momentum budget of the gas. We then review observations and models of feedback through AGN-driven winds from nuclear, sub-pc scales out to galactic and circumgalactic medium scales. We continue with an overview of surveys and statistical properties of the AGN population, before concluding with a discussion on the prospects for the future facilities, focusing in particular on Athena. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: Invited Chapter of Section "Active Galactic Nuclei in X- and Gamma-rays" (Section Editors: A. De Rosa and C. Vignali) of the "Handbook of X-ray and Gamma-ray Astrophysics" (Editors: C. Bambi and A. Santangelo), Springer Nature

The 700 ks Chandra Spiderweb Field II: Evidence for inverse-Compton and thermal diffuse emission in the Spiderweb galaxy

Authors: P. Tozzi, R. Gilli, A. Liu, S. Borgani, M. Lepore, L. Di Mascolo, A. Saro, L. Pentericci, C. Carilli, G. Miley, T. Mroczkowski, M. Pannella, E. Rasia, P. Rosati, C. S. Anderson, A. Calabro', E. Churazov, H. Dannerbauer, C. Feruglio, F. Fiore, R. Gobat, S. Jin, M. Nonino, C. Norman, H. J. A. Rottgering

Abstract: We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spider… ▽ More We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spiderweb galaxy hosts a mildly absorbed quasar, with modest yet significant variability on a timescale of ~1 year. We find that the emission in the jet regions is well described by a power law with Gamma~2-2.5, and it is consistent with IC upscattering of the CMB photons by the relativistic electrons. We also find a roughly symmetric, diffuse emission within a radius of ~100 kpc. This emission is consistent with thermal bremsstrahlung from a hot ICM with a temperature of kT=2.0_{-0.4}^{+0.7} keV, and a metallicity of Z<1.6Z_sun. The average electron density within 100 kpc is n_e=(1.51+-0.24+-0.14)E-2 cm^{-3}, corresponding to an upper limit for the total ICM mass of <=(1.76+-0.30+-0.17)E+12 M_sun (where error bars are 1 sigma statistical and systematic, respectively). If we apply hydrostatic equilibrium to the ICM, we measure a total gravitational mass M(<100 kpc)=(1.5^{+0.5}_{-0.3})E+13 M_sun and, extrapolating at larger radii, we estimate a total mass M_{500}=(3.2^{+1.1}_{-0.6})E+13 M_sun within a radius of r_{500}=(220+-30) kpc. We conclude that the Spiderweb protocluster shows significant diffuse emission within a radius of 12 arcsec, whose major contribution is provided by IC scattering associated with the radio jets. Outside the jet regions, we also identified thermal emission within a radius of ~100 kpc, revealing the presence of hot, diffuse baryons that may represent the embryonic virialized halo of the forming cluster. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 19 pages, 16 figures, accepted in Astronomy & Astrophysics

The black hole and host galaxy growth in an isolated $z\sim 6$ QSO observed with ALMA

Authors: R. Tripodi, C. Feruglio, F. Fiore, M. Bischetti, V. D'Odorico, S. Carniani, S. Cristiani, S. Gallerani, R. Maiolino, A. Marconi, A. Pallottini, E. Piconcelli, L. Vallini, T. Zana

Abstract: The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass (… ▽ More The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass ($M_{\rm dyn}$) plane, we compare the star formation rate (SFR), derived from the accurate estimate of the dust temperature and the dust mass ($T_{\rm dust}, M_{\rm dust}$), with the BH accretion rate. To this aim, we analysed a deep, $900$ pc resolution ALMA observation of the sub-mm continuum, [CII] and H$_2$O of the $z\sim 6$ QSO J2310+1855, enabling a detailed study of dust properties and cold gas kinematics. We performed an accurate SED analysis obtaining a dust temperature of $T_{\rm dust} = 71$ K and a dust mass of $M_{\rm dust}= 4.4 \times 10^8\ \rm M_{\odot}$. The implied AGN-corrected SFR is $1240 \ \rm M_{\odot}yr^{-1}$, a factor of 2 smaller than previously reported for this QSO. We derived the best estimate of the dynamical mass $M_{\rm dyn} = 5.2\times 10^{10}\ \rm M_{\odot}$ within $r = 1.7$ kpc, based on a dynamical model of the system. We found that ${\rm SFR}/M_{\rm dyn}>\dot M_{\rm BH}/M_{\rm BH}$, suggesting that AGN feedback might be efficiently acting to slow down the SMBH accretion, while the stellar mass assembly is still vigorously taking place in the host galaxy. In addition, we were also able to detect high-velocity emission on the red and blue sides of the [CII] emission line, that traces a gaseous outflow, and for the first time, we mapped a spatially-resolved water vapour disk through the H$_2$O v=0 $3_{(2,2)}-3_{(1,3)}$ emission line detected at $ν_{\rm obs} = 274.074$ GHz, whose kinematic properties and size are broadly consistent with those of the [CII] disk. △ Less

Submitted 30 July, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures, 4 tables. Accepted in A&A

Journal ref: A&A 665, A107 (2022)

The WISSH quasars project X. Discovery of a multi-component and highly-variable UV ultra-fast outflow in a z=3.6 quasar

Authors: G. Vietri, T. Misawa, E. Piconcelli, P. Franzetti, A. Luminari, A. Travascio, M. Bischetti, S. Bisogni, A. Bongiorno, G. Bruni, C. Feruglio, A. Giunta, F. Nicastro, I. Saccheo, V. Testa, F. Tombesi, C. Vignali, L. Zappacosta, F. Fiore

Abstract: We report on the variability of a multi-component broad absorption line (BAL) system observed in the hyper-luminous quasar J1538+0855 at z=3.6. Observations from SDSS, VLT, LBT and Subaru telescopes taken at five different epochs, spanning 17 yr in the observed frame, are presented. We detect three (A, B, C) CIV variable troughs exhibiting extreme velocities ($\sim$40,000-54,000 km s$^{-1}$) simil… ▽ More We report on the variability of a multi-component broad absorption line (BAL) system observed in the hyper-luminous quasar J1538+0855 at z=3.6. Observations from SDSS, VLT, LBT and Subaru telescopes taken at five different epochs, spanning 17 yr in the observed frame, are presented. We detect three (A, B, C) CIV variable troughs exhibiting extreme velocities ($\sim$40,000-54,000 km s$^{-1}$) similar to the ultra-fast outflows (UFOs) typically observed in the X-ray spectra. The A component of the BAL UFO ($\rm v_{ufo}$ $\sim$0.17 c) shows strength variations, while B ($\rm v_{ufo}$ $\sim$0.15 c) and C ($\rm v_{ufo}$ $\sim$0.13 c) components show changes both in shape and strength, appearing and disappearing at different epochs. In addition, during the last observation on June 2021 the entire BAL system disappears. The variability trends observed during the first two epochs (1.30 yr rest-frame) in the CIV, SiIV, OVI and NV absorption spectral regions are the same for B and C troughs, while the A component of the BAL varies independently. This suggests a change in the ionization state of the absorbing gas for B and C components and tangential motion for the A component, as causes of this temporal behavior. Accordingly, it is possible to provide an upper limit for distance of the gas responsible for the A component of $R\rm_{out}^{A}$$\le$58 pc, and in turn, a kinetic power of $\dot{E}\rm_{K,ufo}$ $\le$5.2 $\times$ 10$^{44}$ erg s$\rm^{-1}$. We also obtain $R\rm_{out}^{B,C}$ $\le$2.7 kpc for B and C components, which implies an upper limit estimation of $\dot{E}\rm_{K,ufo}$ $\le$2.1$\times$10$^{46}$ erg s$\rm^{-1}$ and $\dot{E}\rm_{K,ufo}$ $\le$1.4$\times$10$^{46}$ erg s$\rm^{-1}$, respectively. Future spectral monitoring with high-resolution instruments is mandatory to accurately constrain physical properties of the BAL UFO discovered in the UV spectrum of J1538+0855. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 668, A87 (2022)

Widespread, strong outflows in XQR-30 quasars at the Reionisation epoch

Authors: M. Bischetti, C. Feruglio, V. D'Odorico, N. Arav, E. Bañados, G. Becker, S. E. I. Bosman, S. Carniani, S. Cristiani, G. Cupani, R. Davies, A. C. Eilers, E. P. Farina, A. Ferrara, R. Maiolino, C. Mazzucchelli, A. Mesinger, R. Meyer, M. Onoue, E. Piconcelli, E. Ryan-Weber, J-T. Schindler, F. Wang, J. Yang, Y. Zhu , et al. (1 additional authors not shown)

Abstract: Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dy… ▽ More Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dynamical mass correlation, thus following the black hole-dominance growth path. It is unknown what slowed down the black hole growth, leading towards the symbiotic growth observed in the local Universe, and when this process started, although black hole feedback is a likely driver. This deadlock is due to the lack of large, homogeneous samples of high-redshift quasars with high-quality, broad-band spectroscopic information. Here we report results from a VLT/X-shooter survey of 30 quasars at redshift 5.8$\le$z$\le$6.6 (XQR-30). About 50% of their spectra reveal broad blue-shifted absorption line (BAL) throughs, tracing powerful ionised winds. The BAL fraction in z$\gtrsim$6 quasars is 2-3 times higher than in quasars at z~2-4.5. XQR-30 BAL quasars exhibit extreme outflow velocities, up to 17% of the light speed, rarely observed at lower redshift. These outflows inject large amounts of energy into the galaxy interstellar medium, which can contrast nuclear gas accretion, slowing down the black-hole growth. The star-formation rate in high-z quasar hosts is generally $>$100 M$_\odot$/yr, so these galaxies are growing at a fast rate. The BAL phase may then mark the beginning of significant feedback, acting first on black hole growth and possibly later on galaxy growth. The red optical colors of BAL quasars at z$\gtrsim$6 indeed suggest that these systems are dusty and may be caught during an initial quenching phase of obscured accretion. △ Less

Submitted 17 May, 2022; v1 submitted 29 April, 2022; originally announced May 2022.

Comments: This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature, and is available online at https://doi.org/10.1038/s41586-022-04608-1

Journal ref: Nature, 605, 244-247 (2022)

The 700 ks Chandra Spiderweb Field I: evidence for widespread nuclear activity in the Protocluster

Authors: P. Tozzi, L. Pentericci, R. Gilli, M. Pannella, F. Fiore, G. Miley, M. Nonino, H. J. A. Rottgering, V. Strazzullo, C. S. Anderson, S. Borgani, A. Calabro', C. Carilli, H. Dannerbauer, L. Di Mascolo, C. Feruglio, R. Gobat, S. Jin, A. Liu, T. Mroczkowski, C. Norman, E. Rasia, P. Rosati, A. Saro

Abstract: (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and sub… ▽ More (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and submm bands to identify X-ray sources belonging to the protocluster. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11<z<2.20, and 1 source with photometric redshift consistent with this range. Our X-ray spectral analysis shows that their intrinsic spectral slope is consistent with an average <Γ>~1.84+-0.04. The best-fit intrinsic absorption for 5 protocluster X-ray members is N_H>10^{23} cm^{-2}, while other 6 have upper limits of the order of fewX10^{22} cm^{-2}. Two sources can only be fitted with very flat Γ<=1, and are therefore considered Compton-thick candidates. Their 0.5-10 keV rest frame luminosities are larger than 2X10^{43} erg/s, significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of AGN in the volume associated to the Spiderweb protocluster in the range 10^{43}<L_X<10^{44.5} erg/s, is at least 10 times higher than that in the field at the same redshift and significantly flatter. The X-ray AGN fraction is measured to be (25.5+-4.5)% in the stellar mass range log(M*/M_sun)>10.5, corresponding to an enhancement of 6.0^{+9.0}_{-3.0} with respect to the COSMOS field at comparable redshifts and stellar mass range. We conclude that the galaxy population in the Spiderweb Protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on Mpc scales. △ Less

Submitted 4 April, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: Astronomy & Astrophysics. Minor changes in this version. Version accepted in A&A

Journal ref: A&A 662, A54 (2022)

The diverse cold molecular gas contents, morphologies, and kinematics of type-2 quasars as seen by ALMA

Authors: C. Ramos Almeida, M. Bischetti, S. Garcia-Burillo, A. Alonso-Herrero, A. Audibert, C. Cicone, C. Feruglio, C. N. Tadhunter, J. C. S. Pierce, M. Pereira-Santaella, P. S. Bessiere

Abstract: We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, w… ▽ More We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, whereas in the spirals it is only 5-12%. J1010+0612 and J1430+1339 show double-peaked CO morphologies which do not have optical counterparts. Based on our analysis of the ionized and molecular kinematics and mm continuum emission, these CO morphologies are most likely produced by AGN feedback in the form of outflows, jets, and/or shocks. The CO kinematics of the QSO2s are dominated by rotation but also reveal noncircular motions. According to our analysis of the kinematics, these noncircular motions correspond to molecular outflows mostly coplanar with the CO discs in four of the QSO2s, and either to a coplanar inflow or vertical outflow in the case of J1010+0612. These outflows represent 0.2-0.7% of the QSO2s' total molecular gas mass and have maximum velocities of 200-350 km/s, radii from 0.4 to 1.3 kpc, and outflow rates of 8-16 Msun/yr. These properties are intermediate between those of the mild molecular outflows measured for Seyferts, and the fast and energetic outflows of ULIRGs. This suggests that it is not only AGN luminosity that drives massive molecular outflows. Other factors such as jet power, coupling between winds, jets, and/or ionized outflows and the CO discs, and amount or geometry of dense gas in the nuclear regions might be also relevant. Thus, although we do not find evidence for a significant impact of quasar feedback on the total molecular gas reservoirs and SFRs, it appears to be modifying the distribution of cold molecular gas in the central kpc of the galaxies. △ Less

Submitted 15 March, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 29 pages including 18 figures, 7 tables and appendices. A&A accepted version including minor edits

Journal ref: A&A 658, A155 (2022)

SUPER V. ALMA continuum observations of z~2 AGN and the elusive evidence of outflows influencing star formation

Authors: I. Lamperti, C. M. Harrison, V. Mainieri, D. Kakkad, M. Perna, C. Circosta, J. Scholtz, S. Carniani, C. Cicone, D. M. Alexander, M. Bischetti, G. Calistro Rivera, C. -C. Chen, G. Cresci, C. Feruglio, F. Fiore, F. Mannucci, A. Marconi, L. N. Martínez-Ramírez, H. Netzer, E. Piconcelli, A. Puglisi, D. J. Rosario, M. Schramm, G. Vietri , et al. (2 additional authors not shown)

Abstract: We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$α$ 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… ▽ More We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$α$ 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$α$ 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$α$ 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. △ Less

Submitted 3 August, 2023; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 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

Journal ref: A&A 654, A90 (2021)

Capturing dual AGN activity and kiloparsec-scale outflows in IRAS 20210+1121

Authors: Francesco Gabriele Saturni, Giustina Vietri, Enrico Piconcelli, Christian Vignali, Manuela Bischetti, Angela Bongiorno, Sara Cazzoli, Chiara Feruglio, Fabrizio Fiore, Bernd Husemann, Cristina Ramos Almeida

Abstract: The most accepted scenario for the evolution of massive galaxies across cosmic time predicts a regulation based on the interplay between AGN feedback, which injects large amounts of energy in the host environment, and galaxy mergers, being able to trigger massive star formation events and accretion onto the supermassive black holes. Interacting systems hosting AGN are useful laboratories to get ke… ▽ More The most accepted scenario for the evolution of massive galaxies across cosmic time predicts a regulation based on the interplay between AGN feedback, which injects large amounts of energy in the host environment, and galaxy mergers, being able to trigger massive star formation events and accretion onto the supermassive black holes. Interacting systems hosting AGN are useful laboratories to get key insights into both phenomena. In this context, we present the analysis of the optical spectral properties of IRAS 20210+1121 (I20210), a merging system at $z = 0.056$. According to X-ray data, this object comprises two interacting galaxies, each hosting an obscured AGN. The optical spectra confirm the presence of AGN features in both galaxies. In particular, we are able to provide a Seyfert classification for I20210 North. The spectrum of I20120 South shows broad blueshifted components associated with the most intense emission lines that indicate the presence of an ionized outflow, for which we derive a maximum velocity of $\sim$2000 km s$^{-1}$, an extension of $\sim$2 kpc and a mass rate of $\sim$0.6 M$_\odot$ yr$^{-1}$. We also report the existence of an ionized nebular component with $v \sim 1000$ km s$^{-1}$ at $\sim$6.5 kpc Southwards of I20210 South, that can be interpreted as disrupted gas ejected from the host galaxy by the action of the outflow. I20120 therefore exhibits a double obscured AGN, with one of them showing evidence of ongoing events for AGN-powered outflows. Future spatially-resolved spectroscopy will allow to accurately map on the gas kinematics in this AGN pair and evaluate the impact of the outflow on both the interstellar medium and galaxy environment. △ Less

Submitted 9 August, 2021; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 16 pages, 9 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 654, A154 (2021)

The IBISCO survey: I. Multiphase discs and winds in the Seyfert galaxy Markarian 509

Authors: M. V. Zanchettin, C. Feruglio, M. Bischetti, A. Malizia, M. Molina, A. Bongiorno, M. Dadina, C. Gruppioni, E. Piconcelli, F. Tombesi, A. Travascio, F. Fiore

Abstract: We present the analysis of the ALMA CO(2-1) emission line and the underlying 1.2 mm continuum of Mrk509 with spatial resolution of 270 pc. This local Seyfert 1.5 galaxy, optically classified as a spheroid, is known to host a ionised disc, a starburst ring, and ionised gas winds on both nuclear and galactic scales. From CO(2-1) we estimate a molecular gas mass… ▽ More We present the analysis of the ALMA CO(2-1) emission line and the underlying 1.2 mm continuum of Mrk509 with spatial resolution of 270 pc. This local Seyfert 1.5 galaxy, optically classified as a spheroid, is known to host a ionised disc, a starburst ring, and ionised gas winds on both nuclear and galactic scales. From CO(2-1) we estimate a molecular gas mass $M_{H_2}=1.7\times 10^9\, \rm M_{\odot}$, located within a disc of size 5.2 kpc, with $M_{dyn}$=(2.0$\pm$1.1) $\times$ $10^{10}\, \rm M_{\odot}$ inclined at $44\pm10$ deg. The molecular gas fraction within the disc is $μ_{gas}=5\%$. The gas kinematics in the nuclear region within r=700 pc suggests the presence of a warped nuclear disc. Both the presence of a molecular disc with ongoing star-formation in a starburst ring, and the signatures of a minor merger, are in agreement with the scenario where galaxy mergers produce gas destabilization, feeding both star-formation and AGN activity. The spatially-resolved Toomre Q-parameter across the molecular disc is in the range $Q_{gas}=0.5-10$, and shows that the disc is marginally unstable across the starburst ring, and stable at nucleus and in a lopsided ring-like structure located inside of the starburst ring. We find complex molecular gas kinematics and significant kinematics perturbations at two locations, one within 300 pc from the nucleus, and one 1.4 kpc away close to the region with high $Q_{gas}$, that we interpret as molecular winds. The total molecular outflow rate is in the range 6.4-17.0 $\rm M_\odot/yr$. The molecular wind total kinetic energy is consistent with a multiphase momentum-conserving wind driven by the AGN with $\dot{P}_{of}/\dot{P}_{rad}$ in the range 0.06-0.5. The spatial overlap of the inner molecular wind with the ionised wind, and their similar velocity suggest a cooling sequence within a multiphase AGN driven wind. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 13 pages, 11 figures, 3 tables Accepted by A&A

Journal ref: A&A 655, A25 (2021)

SUPER IV. CO(J=3-2) properties of active galactic nucleus hosts at cosmic noon revealed by ALMA

Authors: C. Circosta, V. Mainieri, I. Lamperti, P. Padovani, M. Bischetti, C. M. Harrison, D. Kakkad, A. Zanella, G. Vietri, G. Lanzuisi, M. Salvato, M. Brusa, S. Carniani, C. Cicone, G. Cresci, C. Feruglio, B. Husemann, F. Mannucci, A. Marconi, M. Perna, E. Piconcelli, A. Puglisi, A. Saintonge, M. Schramm, C. Vignali , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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) △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: Key figures: 2, 3, 4. 22 pages, 9 figures (5 in appendices), 3 tables. Accepted for publication in A&A

Journal ref: A&A 646, A96 (2021)

Evidence of galaxy interaction in the Narrow-line Seyfert 1 galaxy IRAS17020+4544 seen by NOEMA

Authors: Q. Salomé, A. L. Longinotti, Y. Krongold, C. Feruglio, V. Chavushyan, O. Vega, S. García-Burillo, A. Fuente, A. Olguín-Iglesias, V. M. Patiño-Álvarez, I. Puerari, A. Robleto-Orús

Abstract: The narrow-line Seyfert 1 galaxy IRAS17020+4544 is one of the few sources where both an X-ray ultra-fast outflow and a molecular outflow were observed to be consistent with energy conservation. However, IRAS17020+4544 is less massive and has a much more modest active galactic nucleus (AGN) luminosity than the other examples. Using recent CO(1-0) observations with the NOrthern Extended Millimeter A… ▽ More The narrow-line Seyfert 1 galaxy IRAS17020+4544 is one of the few sources where both an X-ray ultra-fast outflow and a molecular outflow were observed to be consistent with energy conservation. However, IRAS17020+4544 is less massive and has a much more modest active galactic nucleus (AGN) luminosity than the other examples. Using recent CO(1-0) observations with the NOrthern Extended Millimeter Array (NOEMA), we characterised the molecular gas content of the host galaxy for the first time. We found that the molecular gas is distributed into an apparent central disc of 1.1x10^9 Msun, and a northern extension located up to 8 kpc from the centre with a molecular gas mass M_H2~10^8 Msun. The molecular gas mass and the CO dynamics in the northern extension reveal that IRAS 17020+4544 is not a standard spiral galaxy, instead it is interacting with a dwarf object corresponding to the northern extension. This interaction possibly triggers the high accretion rate onto the super massive black hole. Within the main galaxy, which hosts the AGN, a simple analytical model predicts that the molecular gas may lie in a ring, with less molecular gas in the nuclear region. Such distribution may be the result of the AGN activity which removes or photodissociates the molecular gas in the nuclear region (AGN feedback). Finally, we have detected a molecular outflow of mass M_H2=(0.7-1.2)x10^7 Msun in projection at the location of the northern galaxy, with a similar velocity to that of the massive outflow reported in previous millimeter data obtained by the Large Millimeter Telescope. △ Less

Submitted 6 January, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: Published in MNRAS, Volume 501, Issue 1, Pages 219-228

SUPER III. Broad Line Region properties of AGN at z$\sim$2

Authors: G. Vietri, V. Mainieri, D. Kakkad, H. Netzer, M. Perna, C. Circosta, C. M. Harrison, L. Zappacosta, B. Husemann, P. Padovani, M. Bischetti, A. Bongiorno, M. Brusa, S. Carniani, C. Cicone, A. Comastri, G. Cresci, C. Feruglio, F. Fiore, G. Lanzuisi, F. Mannucci, A. Marconi, E. Piconcelli, A. Puglisi, M. Salvato , et al. (5 additional authors not shown)

Abstract: 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… ▽ More 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$α$ and H$β$ 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λ_{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λ_{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. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 644, A175 (2020)

The WISSH QSOs project IX. Cold gas content and environment of luminous QSOs at z~2.4-4.7

Authors: M. Bischetti, C. Feruglio, E. Piconcelli, F. Duras, M. Pèrez-Torres, R. Herrero, G. Venturi, S. Carniani, G. Bruni, I. Gavignaud, V. Testa, A. Bongiorno, M. Brusa, C. Circosta, G. Cresci, V. D'Odorico, R. Maiolino, A. Marconi, M. Mingozzi, C. Pappalardo, M. Perna, E. Traianou, A. Travascio, G. Vietri, L. Zappacosta , et al. (1 additional authors not shown)

Abstract: Sources at the brightest end of QSO luminosity function during the peak epoch of star formation and black hole accretion (z~2-4, i.e. Cosmic noon) are privileged sites to study the feeding & feedback cycle of massive galaxies. We perform the first systematic study of cold gas properties in the most luminous QSOs, by characterising their host-galaxies and environment. We analyse ALMA, NOEMA and JVL… ▽ More Sources at the brightest end of QSO luminosity function during the peak epoch of star formation and black hole accretion (z~2-4, i.e. Cosmic noon) are privileged sites to study the feeding & feedback cycle of massive galaxies. We perform the first systematic study of cold gas properties in the most luminous QSOs, by characterising their host-galaxies and environment. We analyse ALMA, NOEMA and JVLA observations of FIR continuum, CO and [CII] emission lines in eight QSOs ($L_{\rm Bol}>3\times10^{47}$ erg/s) from the WISSH sample at z~2.4-4.7. 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-mm 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), associated with rotating disks with disturbed kinematics. WISSH QSOs typically show lower CO luminosity and higher star formation efficiency than FIR matched, z~0-3 main-sequence galaxies, implying that, given the observed SFR ~170-1100 $M_\odot$/yr, molecular gas is converted into stars on <50 Myr. Most targets show extreme dynamical to black-hole mass ratios $M_{\rm dyn}/M_{\rm BH}\sim3-10$, two orders of magnitude smaller than local relations. The molecular gas fraction in WISSH hosts is lower by a factor of ~10-100 than in star forming galaxies with similar $M_*$. WISSH QSOs undergo an intense growth phase of both the central SMBH and host-galaxy. They pinpoint high-density sites where giant galaxies assemble and mergers play a major role in the build-up of the final host-galaxy mass. The observed low molecular gas fraction and short depletion timescale are likely due to AGN feedback, as traced by fast AGN-driven ionised outflows in all our targets. △ Less

Submitted 13 October, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 645, A33 (2021)

SUPER-II: Spatially resolved ionized gas kinematics and scaling relations in z~2 AGN host galaxies

Authors: D. Kakkad, V. Mainieri, G. Vietri, S. Carniani, C. M. Harrison, M. Perna, J. Scholtz, C. Circosta, G. Cresci, B. Husemann, M. Bischetti, C. Feruglio, F. Fiore, A. Marconi, P. Padovani, M. Brusa, C. Cicone, A. Comastri, G. Lanzuisi, F. Mannucci, N. Menci, H. Netzer, E. Piconcelli, A. Puglisi, M. Salvato , et al. (5 additional authors not shown)

Abstract: 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}}$ =… ▽ More 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%. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: Accepted for publication in A&A, 34 pages, 14 figures and 5 tables

Journal ref: A&A 642, A147 (2020)

The Rise of Active Galactic Nuclei in the GAlaxy Evolution and Assembly semi-analytic model

Authors: Fabio Fontanot, Gabriella De Lucia, Michaela Hirschmann, Lizhi Xie, Pierluigi Monaco, Nicola Menci, Fabrizio Fiore, Chiara Feruglio, Stefano Cristiani, Francesco Shankar

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

Submitted 12 June, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 19 pages, 12 figures, 1 table, replaced with MNRAS accepted version

The WISSH quasars project VII. The impact of extreme radiative field in the accretion disk and X-ray corona interplay

Authors: L. Zappacosta, E. Piconcelli, M. Giustini, G. Vietri, F. Duras, G. Miniutti, M. Bischetti, A. Bongiorno, M. Brusa, M. Chiaberge, A. Comastri, C. Feruglio, A. Luminari, A. Marconi, C. Ricci, C. Vignali, F. Fiore

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

Submitted 17 March, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 10 pages (including Appendix), 4 figures. Astronomy & Astrophysics Letters in press. Title numbering changed to match the published version

Journal ref: A&A 635, L5 (2020)

The WISSH quasars project VII. Outflows and metals in the circumgalactic medium around the hyper-luminous z~3.6 quasar J1538+08

Authors: A. Travascio, L. Zappacosta, S. Cantalupo, E. Piconcelli, F. Arrigoni Battaia, M. Ginolfi, M. Bischetti, G. Vietri, A. Bongiorno, V. D'Odorico, F. Duras, C. Feruglio, C. Vignali, F. Fiore

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

Submitted 20 January, 2020; originally announced January 2020.

Comments: 14 pages, 16 figures, 2 tables. Accepted for publication in A&A Journal

Enhanced UV radiation and dense clumps in Mrk231's molecular outflow

Authors: Claudia Cicone, Roberto Maiolino, Susanne Aalto, Sebastien Muller, Chiara Feruglio

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

Submitted 25 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A, 14 pages, 7 figures, 4 tables

Journal ref: A&A 633, A163 (2020)

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

Authors: D. Lutz, E. Sturm, A. Janssen, S. Veilleux, S. Aalto, C. Cicone, A. Contursi, R. I. Davies, C. Feruglio, J. Fischer, S. Garcia-Burillo, R. Genzel, E. González-Alfonso, J. Gracía-Carpio R. Herrera-Camus, R. Maiolino, A. Schruba, T. Shimizu, A. Sternberg, L. J. Tacconi, A. Weiß

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

Submitted 13 November, 2019; originally announced November 2019.

Comments: 37 pages. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 633, A134 (2020)

The WISSH quasars project VI. Fraction and properties of BAL quasars in the hyper-luminosity regime

Authors: G. Bruni, E. Piconcelli, T. Misawa, L. Zappacosta, F. G. Saturni, G. Vietri, C. Vignali, A. Bongiorno, F. Duras, C. Feruglio, F. Tombesi, F. Fiore

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

Submitted 26 August, 2019; originally announced August 2019.

Comments: Accepted for publication on A&A

Constraints on gamma-ray and neutrino emission from NGC 1068 with the MAGIC telescopes

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (150 additional authors not shown)

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

Submitted 12 August, 2019; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 11 pages, 4 figures, accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal 883 (2019) 135

Outflows in the Discs of Active Galaxies

Authors: N. Menci, F. Fiore, C. Feruglio, A. Lamastra, F. Shankar, E. Piconcelli, E. Giallongo, A. Grazian

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

Submitted 16 April, 2019; originally announced April 2019.

Comments: 34 pages, accepted for publication in The Astrophysical Journal

Multiphase gas flows in the nearby Seyfert galaxy ESO428-G14

Authors: C. Feruglio, G. Fabbiano, M. Bischetti, M. Elvis, A. Travascio, F. Fiore

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

Submitted 2 April, 2019; originally announced April 2019.

Comments: Submitted to ApJ