skip to main content

Title: Red quasars blow out molecular gas from galaxies during the peak of cosmic star formation

Recent studies have suggested that red quasars are a phase in quasar evolution when feedback from black hole accretion evacuates obscuring gas from the nucleus of the host galaxy. Here, we report a direct link between dust-reddening and molecular outflows in quasars at z ∼ 2.5. By examining the dynamics of warm molecular gas in the inner region of galaxies, we find evidence for outflows with velocities 500–1000 km s−1 and time-scales of ≈0.1 Myr that are due to ongoing quasar energy output. We infer outflows only in systems where quasar radiation pressure on dust in the vicinity of the black hole is sufficiently large to expel their obscuring gas column densities. This result is in agreement with theoretical models that predict radiative feedback regulates gas in the nuclear regions of galaxies and is a major driving mechanism of galactic-scale outflows of cold gas. Our findings suggest that radiative quasar feedback ejects star-forming gas from within nascent stellar bulges at velocities comparable to those seen on larger scales, and that molecules survive in outflows even from the most luminous quasars.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 3377-3391
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We present a comprehensive study of the molecular gas properties of 17 Type 2 quasars at z < 0.2 from the Quasar Feedback Survey (L$_{\rm [O~{\small III}]}$ > 1042.1 $\rm ergs^{-1}$), selected by their high [O iii] luminosities and displaying a large diversity of radio jet properties, but dominated by LIRG-like galaxies. With these data, we are able to investigate the impact of AGN and AGN feedback mechanisms on the global molecular interstellar medium. Using Atacama Pathfinder EXperiment and ALMA ACA observations, we measure the total molecular gas content using the CO(1-0) emission and homogeneously sample the carbon monoxide (CO) spectral line energy distributions, observing CO transitions (Jup  = 1, 2, 3, 6, 7). We observe high r21 ratios (r21  = L’CO(2-1)/L’CO(1-0)) with a median r21  = 1.06, similar to local (U)LIRGs (with r21 ∼ 1) and higher than normal star-forming galaxies (with r21 ∼ 0.65). Despite the high r21 values, for the seven targets with the required data, we find low excitation in CO(6-5) & CO(7-6) (r61 and r62 < 0.6 in all but one target), unlike high-redshift quasars in the literature, which are far more luminous and show higher line ratios. The ionized gas traced by [O iii] exhibits systematically higher velocities than the molecular gas traced by CO. We conclude that any effects of quasar feedback (e.g. via outflows and radio jets) do not have a significant instantaneous impact on the global molecular gas content and excitation and we suggest that it only occurs on more localized scales.

    more » « less

    Molecular gas flows are analysed in 14 cluster galaxies (BCGs) centred in cooling hot atmospheres. The BCGs contain $10^{9}\!-\!10^{11}~\rm M_\odot$ of molecular gas, much of which is being moved by radio jets and lobes. The molecular flows and radio jet powers are compared to molecular outflows in 45 active galaxies within z < 0.2. We seek to understand the relative efficacy of radio, quasar, and starburst feedback over a range of active galaxy types. Molecular flows powered by radio feedback in BCGs are ∼10–1000 times larger in extent compared to contemporary galaxies hosting quasar nuclei and starbursts. Radio feedback yields lower flow velocities but higher momenta compared to quasar nuclei, as the molecular gas flows in BCGs are usually ∼10–100 times more massive. The product of the molecular gas mass and lifting altitude divided by the AGN or starburst power – a parameter referred to as the lifting factor – exceeds starbursts and quasar nuclei by 2–3 orders of magnitude, respectively. When active, radio feedback is generally more effective at lifting gas in galaxies compared to quasars and starburst winds. The kinetic energy flux of molecular clouds generally lies below and often substantially below a few per cent of the driving power. We find tentatively that star formation is suppressed in BCGs relative to other active galaxies, perhaps because these systems rarely form molecular discs that are more impervious to feedback and are better able to promote star formation.

    more » « less

    Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) galaxies, including normal star-forming galaxies, quasar hosts, and dusty star-forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star-formation histories and galactic outflow properties of galaxies in massive haloes ($10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$) at z ≳ 5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ∼0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium. We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems.

    more » « less
  4. Abstract

    We use medium-resolution Keck/Echellette Spectrograph and Imager spectroscopy of bright quasars to study cool gas traced by Caiiλλ3934, 3969 and Naiλλ5891, 5897 absorption in the interstellar/circumgalactic media of 21 foreground star-forming galaxies at redshifts 0.03 <z< 0.20 with stellar masses 7.4 ≤ logM*/M≤ 10.6. The quasar–galaxy pairs were drawn from a unique sample of Sloan Digital Sky Survey quasar spectra with intervening nebular emission, and thus have exceptionally close impact parameters (R< 13 kpc). The strength of this line emission implies that the galaxies’ star formation rates (SFRs) span a broad range, with several lying well above the star-forming sequence. We use Voigt profile modeling to derive column densities and component velocities for each absorber, finding that column densitiesN(Caii) > 1012.5cm−2(N(Nai) > 1012.0cm−2) occur with an incidencefC(Caii) = 0.63+0.10−0.11(fC(Nai) = 0.57+0.10−0.11). We find no evidence for a dependence offCor the rest-frame equivalent widthsWr(CaiiK) orWr(Nai5891) onRorM*. Instead,Wr(CaiiK) is correlated with local SFR at >3σsignificance, suggesting that Caiitraces star formation-driven outflows. While most of the absorbers have velocities within ±50 km s−1of the host redshift, their velocity widths (characterized by Δv90) are universally 30–177 km s−1larger than that implied by tilted-ring modeling of the velocities of interstellar material. These kinematics must trace galactic fountain flows and demonstrate that they persist atR> 5 kpc. Finally, we assess the relationship between dust reddening andWr(CaiiK) (Wr(Nai5891)), finding that 33% (24%) of the absorbers are inconsistent with the best-fit Milky WayE(B−V)-Wrrelations at >3σsignificance.

    more » « less
  5. Water vapor (H2O) is one of the brightest molecular emitters after carbon monoxide (CO) in galaxies with high infrared (IR) luminosity, allowing us to investigate the warm and dense phase of the interstellar medium (ISM) where star formation occurs. However, due to the complexity of its radiative spectrum, H2O is not frequently exploited as an ISM tracer in distant galaxies. Therefore, H2O studies of the warm and dense gas at high-zremain largely unexplored. In this work, we present observations conducted with the Northern Extended Millimeter Array (NOEMA) toward threez > 6 IR-bright quasarsJ2310+1855,J1148+5251, andJ0439+1634targeted in their multiple para- and ortho-H2O transitions (312 − 303, 111 − 000, 220 − 211, and 422 − 413), as well as their far-IR (FIR) dust continuum. By combining our data with previous measurements from the literature, we estimated the dust masses and temperatures, continuum optical depths, IR luminosities, and star formation rates (SFR) from the FIR continuum. We modeled the H2O lines using the MOLPOP-CEP radiative transfer code, finding that water vapor lines in our quasar host galaxies are primarily excited in the warm, dense (with a gas kinetic temperature and density ofTkin = 50 K,nH2 ∼ 104.5 − 105 cm−3) molecular medium with a water vapor column density ofNH2O ∼ 2 × 1017 − 3 × 1018 cm−3. High-JH2O lines are mainly radiatively pumped by the intense optically-thin far-IR radiation field associated with a warm dust component at temperatures ofTdust ∼ 80 − 190 K that account for < 5 − 10% of the total dust mass. In the case of J2310+1855, our analysis points to a relatively high value of the continuum optical depth at 100 μm (τ100 ∼ 1). Our results are in agreement with expectations based on the H2O spectral line energy distribution of local and high-zultra-luminous IR galaxies and active galactic nuclei (AGN). The analysis of the Boltzmann diagrams highlights the interplay between collisions and IR pumping in populating the high H2O energy levels and it allows us to directly compare the excitation conditions in the targeted quasar host galaxies. In addition, the observations enable us to sample the high-luminosity part of the H2O–total-IR (TIR) luminosity relations (LH2O − LTIR). Overall, our results point to supralinear trends that suggest H2O–TIR relations are likely driven by IR pumping, rather than the mere co-spatiality between the FIR continuum- and line-emitting regions. The observedLH2O/LTIRratios in ourz > 6 quasars do not show any strong deviations with respect to those measured in star-forming galaxies and AGN at lower redshifts. This supports the notion that H2O can be likely used to trace the star formation activity buried deep within the dense molecular clouds.

    more » « less