MNRAS 000,1{17 (2016) Preprint 24 March 2016 Compiled using MNRAS LATEX style file v3.0 Evidence for the Thermal Sunyaev-Zel'dovich Effect Associated with Quasar Feedback Devin Crichton1?, Megan B. Gralla1;2, Kirsten Hall1, Tobias A. Marriage1, Nadia L. Zakamska1, Nick Battaglia3, J. Richard Bond4, Mark J. Devlin5, J. Colin Hill6, Matt Hilton7, Adam D. Hincks8, Kevin M. Huffenberger9, John P. Hughes10, Arthur Kosowsky11, Kavilan Moodley7, Michael D. Niemack12, Lyman A. Page13, Bruce Partridge14, Jonathan L. Sievers15, Crist´obal Sif´on16, Suzanne T. Staggs13, Marco P. Viero17 and Edward J. Wollack18 1 Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 3 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA 4 Canadian Institute for Theoretical Astrophysics, University of Toronto, Toronto, ON, M5S 3H8, Canada 5 Department of Physics and Astronomy, University of Pennsylvania, 209 South 33rd Street, Philadelphia, PA 19104 USA 6 Department of Astronomy, Pupin Hall, Columbia University, New York, NY 10027, USA 7 Astrophysics and Cosmology Research Unit, School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Durban 4041, South Africa 8 Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada 9 Department of Physics, Florida State University, PO Box 3064350, Tallahassee, FL 32306-4350, USA 10 Department of Physics and Astronomy, Rutgers University, 136 Frelinghuysen Road, Piscataway, NJ 08854-8019, USA 11 Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA, 15260, USA 12 Department of Physics, Cornell University, Ithaca, NY 14853, USA 13 Joseph Henry Laboratories of Physics, Jadwin Hall, Princeton University, Princeton, NJ, 08544, USA 14 Department of Physics and Astronomy, Haverford College, Haverford, PA, 19041, USA 15 Astrophysics and Cosmology Research Unit, School of Chemistry and Physics, University of KwaZulu-Natal, Durban 4041, South Africa 16 Leiden Observatory, Leiden University, PO Box 9513, NL-2300 RA Leiden, Netherlands 17 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA 18 NASA/Goddard Space Flight Center, Greenbelt, MD, 20771, USA Accepted 2016 February 10. Received 2016 January 29; in original form 2015 October 20 ABSTRACT Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalogue, spanning redshifts 0.5{3.5, we derive the mean millimetre and far-infrared quasar spectral energy distributions (SEDs) via a stacking analysis of Atacama Cos- mology Telescope and Herschel-Spectral and Photometric Imaging REceiver data. We constrain the form of the far-infrared emission and find 3σ{4σ evidence for the ther- mal Sunyaev-Zel'dovich (SZ) effect, characteristic of a hot ionized gas component with thermal energy (6:2 ± 1:7) × 1060 erg. This amount of thermal energy is greater than expected assuming only hot gas in virial equilibrium with the dark matter haloes of 12 −1 arXiv:1510.05656v3 [astro-ph.GA] 23 Mar 2016 (1 − 5) × 10 h M that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this ± −1 energy. Our measurements are consistent with quasars depositing up to (14:5 3:3) τ8 per cent of their radiative energy into their circumgalactic environment if their typical 8 13 −1 period of quasar activity is τ8 × 10 yr. For high quasar host masses, ∼ 10 h M , this percentage will be reduced. Furthermore, the uncertainty on this percentage is only statistical and additional systematic uncertainties enter at the 40 per cent level. The SEDs are dust dominated in all bands and we consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3σ{4σ level for most models. Key words: galaxies: active { galaxies: intergalactic medium { quasars: general c 2016 The Authors ? E-mail: [email protected] 2 D. Crichton et al. 1 INTRODUCTION evidence for large-scale impact of these outflows has long been elusive. Only the last few years have brought about Feedback from accreting black holes has become a key ele- a flood of observations of powerful, galaxy-wide outflows ment in modelling galaxy evolution (Tabor & Binney 1993; launched by radio-quiet quasars (Arav et al. 2008; Moe et al. Silk & Rees 1998; Springel et al. 2005). Quasar feedback 2009; Alexander et al. 2010; Greene et al. 2011; Cano-D´ıaz is routinely invoked in galaxy formation models to quench et al. 2012; Liu et al. 2013a,b; Hainline et al. 2013; Rupke star formation and to explain the steep decline of the bright & Veilleux 2013; Veilleux et al. 2013; Harrison et al. 2014; end of the luminosity function of galaxies (Thoul & Wein- Hainline et al. 2014; Zakamska & Greene 2014; Brusa et al. berg 1995; Croton et al. 2006) and to reheat the intracluster 2015a; Carniani et al. 2015; Harrison et al. 2015; Perna et al. medium (e.g., Rawlings & Jarvis 2004; Scannapieco & Oh 2015; Brusa et al. 2015b). 2004). The discovery that the masses of supermassive black Quasar-driven winds are likely to be inhomogeneous, holes in inactive galaxies strongly correlate with the velocity with different phases of the wind medium observable in dispersions and masses of their hosts' stellar bulges (Magor- different domains of the electromagnetic spectrum. The rian et al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt simplest models for these winds involve hot, low density, 2000) also suggests that the energy output of the black hole volume-filling plasma (Faucher-Gigu`ere & Quataert 2012; in its quasar phase must be somehow coupled to the gas Zubovas & King 2012), with higher density clumps, shells from which the stars form (Hopkins et al. 2006). Finally, or filaments. The warmer of those (T ' 104K) produce opti- the similarity of the history of star formation activity and cal emission lines, and cooler (and denser) clumps may exist that of supermassive black hole accretion in the Universe in neutral or even molecular form. Thus, the physical con- (Boyle & Terlevich 1998; Hopkins et al. 2008) suggests that ditions in quasar-driven winds are qualitatively similar to these two phenomena physically affect each other or result those in starburst-driven winds, where the different phases from the same physical process. Both star formation activity are directly observable and spatially resolved in some nearby and black hole accretion activity peak at z > 2 and rapidly galaxies (Heckman et al. 1990; Veilleux et al. 1994). It is at decline afterward; thus it is at these high redshifts that the present not known which of the phases of the quasar winds physical connection between supermassive black holes and carries most of the mass, momentum and energy; obtain- their hosts was established. ing these measurements is critical for understanding the full How can supermassive black holes so profoundly af- impact of quasars on galaxy formation. fect their large-scale galactic and even intergalactic environ- The hot volume-filling component is particularly elu- ments? Even a small fraction of the binding energy of mate- sive because of its extremely low density, which implies that rial accreted by the supermassive black hole is in principle any emission associated with it would be very weak (Greene sufficient to liberate the galaxy-scale gas from the galaxy et al. 2014). The Sunyaev-Zeldovich (SZ) effect (Sunyaev & potential, but coupling this energy output to the surround- Zeldovich 1970) offers a unique opportunity to detect and ing gas is no trivial task. However, there are circumstances characterize the bulk of quasar-driven winds. The SZ ef- in which accretion energy clearly has had an impact on its fect manifests itself as a spectral distortion of the cosmic environment. Radio jets in massive elliptical galaxies and microwave background (CMB) that occurs when CMB pho- brightest cluster galaxies deposit energy into the hot gas en- tons inverse-Compton scatter off intervening hot ionized gas. velope (McNamara & Nulsen 2007; Fabian 2012). Likewise, Most importantly, the magnitude of the total integrated SZ powerful radio jets entrain warm gas and carry significant effect is proportional to the total thermal energy of the hot amounts of material out of their host galaxies, especially ionized gas and thus allows us to measure the thermalized at high redshifts (van Breugel et al. 1986; Tadhunter 1991; energetic output from the AGN, which is not detectable by Villar-Mart´ın et al. 1999; Nesvadba et al. 2006, 2008; Fu & other methods. Furthermore, the surface brightness of the Stockton 2009). While a narrow collimated outflow might SZ effect is not subject to dimming with redshift, facilitat- seem inefficient at impacting and removing the surround- ing the study of high-z systems of low-density with ionized ing matter, in practice the power of the jet is quickly and gas that is too faint to be studied in emission with current efficiently thermalized if it propagates into high density ma- facilities. terial (Begelman & Cioffi 1989). Fortunately, measurements of the SZ effect have re- Radio-quiet quasars (those that do not seem to have cently reached a maturity that enables investigations of powerful radio jets) can also have a strong effect on AGN feedback. Although the SZ effect was theorized over 40 their large-scale environment. Radiation pressure, especially years ago as a CMB spectral distortion by hot gas around through absorption in bound-bound transitions (Murray galaxy clusters, only within the past 5 years have mea- et al. 1995; Proga et al. 2000), can accelerate material in the surements of the SZ effect begun to reach their potential immediate vicinity of the black hole (. 1000 Schwarzschild with detections of hundreds of clusters in millimetre-wave radii) to velocities .