LETTER Oy Ntttsfravne Td,Teuiest Ftoky of University the 2 Study, Advanced for Institutes Tokyo Ate Bershady Matthew Ai .Wake A
Total Page:16
File Type:pdf, Size:1020Kb
LETTER Suppressing star formation in quiescent galaxies with supermassive black hole winds Edmond Cheung1 , Kevin Bundy1, Michele Cappellari2, Sébastien Peirani1,3, Wiphu Rujopakarn1,4, Kyle Westfall5, Renbin Yan6, Matthew Bershady7, Jenny E. Greene8, Timothy M. Heckman9, Niv Drory10, David R. Law11, Karen L. Masters5, Daniel Thomas5, David A. Wake7,12, Anne-Marie Weijmans13, Kate Rubin14, Francesco Belfiore15,16, Benedetta Vulcani1, Yan-mei Chen17, Kai Zhang6, Joseph D. Gelfand18,19, Dmitry Bizyaev20,21, A. Roman-Lopes22, & Donald P. Schneider23,24 Quiescent galaxies with little or no ongoing star formation domi- To illuminate the salient features of this class, we focus on a pro- 10 nate the galaxy population above M∗ 2 10 M⊙, where their totypical example, nicknamed “Akira” (Fig. 1). The SDSS imaging numbers have increased by a factor of∼ ×25 since z 21–4. Once shows Akira to be an unremarkable spheroidal galaxy of moderate ∼ ∼ star formation is initially shut down, perhaps during the quasar stellar mass (log M∗/M⊙ = 10.78) that is interacting with a low-mass phase of rapid accretion onto a supermassive black hole5–7, an un- companion (nicknamed “Tetsuo”) at a projected separation of 32 kpc known mechanism must remove or heat subsequently accreted gas (67′′); they are not classified as members of a larger galaxy group≈ 18 and from stellar mass loss8 or mergers that would otherwise cool to the properties of both galaxies are listed in Table 1. Spectral energy dis- form stars9,10. Energy output from a black hole accreting at a low tribution (SED) fitting indicates that Akira is nearly dormant, with al- rate has been proposed11–13, but observational evidence for this in most no detection of ongoing star formation19. Resolved spectroscopy, the form of expanding hot gas shells is indirect and limited to radio however, reveals intriguing and complex patterns among spectral trac- galaxies at the centers of clusters14,15, which are too rare to explain ers of gas in Akira that point to a much more active internal state. With the vast majority of the quiescent population16. Here we report ionized gas emission detected across the entire galaxy, the map of Hα bisymmetric emission features co-aligned with strong ionized gas EW (which measures the line flux relative to the stellar continuum; velocity gradients from which we infer the presence of centrally- Fig. 1c) reveals a prominent and somewhat twisted bisymmetric pat- driven winds in typical quiescent galaxies that host low-luminosity tern with a position angle (PA) of 46◦. The projected velocity gra- ∼ −1 −1 active nuclei. These galaxies are surprisingly common, accounting dient ranges from vionized gas = −225 km s to vionized gas = 200 km s 10 ◦ for as much as 10% of the population at M∗ 2 10 M⊙. In a along the kinematic major axis, which is at a PA of 26 (Fig. 1h). We prototypical example, we calculate that the energy∼ × input from the observe high ionized gas velocity dispersions across∼ the galaxy with galaxy’s low-level active nucleus is capable of driving the observed interesting internal structure and maxima that reach σionized gas 200 −1 −1 ∼ wind, which contains sufficient mechanical energy to heat ambient, km s (Fig. 1i) and W80 500 km s (see Methods) perpendicular to cooler gas (also detected) and thereby suppress star formation. the major kinematic axis.∼ Meanwhile, stellar motions reveal a mini- − Using optical imaging spectroscopy from the Sloan Digital Sky mal gradient ( 30 km s 1; Fig. 1f) that follows the PA of the galaxy’s ± ◦ Survey-IV Mapping Nearby Galaxies at Apache Point Observatory17 elliptical isophotes of 53 (contours in Fig. 1c). We also detect a ∼ (SDSS-IV MaNGA) program, we define a new class of quiescent spatially offset enhancement in Na D absorption (Fig. 1d) that is coin- galaxies (required to have red rest-frame colors, NUV − r > 5) that is cident with excess dust in our derived extinction map (see Methods). characterized by the presence of narrow bisymmetric patterns in equiv- Measurements of the Na D line center trace a separate and distinct ve- alent width (EW) maps of strong emission lines, such as Hα and [O III]. Our selection employs multiband imaging to exclude galaxies with dust WI 53706, USA lanes and other disk signatures. The observed enhanced emission fea- 8Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA 9 tures are oriented randomly with respect to the optical surface bright- Center for Astrophysical Sciences, Department of Physics & Astronomy, The Johns Hopkins Uni- arXiv:1605.07626v1 [astro-ph.GA] 24 May 2016 versity, Baltimore, Maryland 21218 ness morphology, but roughly align with strong, systematic velocity 10McDonald Observatory, Department of Astronomy, University of Texas at Austin, 1 University Sta- gradients as traced by the ionized gas emission lines. The gas velocity tion, Austin, TX 78712-0259, USA 11Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA fields in these galaxies are decoupled from their stellar motions. These 12Department of Physical Sciences, The Open University, Milton Keynes, MK7 6AA, UK galaxies are surprisingly common among the quiescent population, ac- 13School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife KY16 9SS, UK counting for 10% of quiescent galaxies with log M∗/M⊙ 10.3. ∼ ∼ 14Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 15Cavendish Laboratory, University of Cambridge, 19 J.J. Thomson Ave, CB3 0HE Cambridge, UK 16University of Cambridge, Kavli Institute for Cosmology, CB3 0HE Cambridge, UK 1Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of 17Department of Astronomy, Nanjing University, Nanjing 210093, China Tokyo Institutes for Advanced Study, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 18NYU Abu Dhabi, P.O. Box 129188, Abu Dhabi, UAE 2Sub-department of Astrophysics, Department of Physics, University of Oxford, Denys Wilkinson 19Center for Cosmology and Particle Physics, New York University, Meyer Hall of Physics, 4 Wash- Building, Keble Road, Oxford OX1 3RH, UK ington Place, New York, NY 10003, USA 3Institut d’Astrophysique de Paris (UMR 7095: CNRS and UPMC), 98 bis Bd Arago F-75014 Paris, 20Apache Point Observatory and New Mexico State University, P.O. Box 59, Sunspot, NM, 88349- France 0059, USA 4Department of Physics, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, 21Sternberg Astronomical Institute, Moscow State University, Moscow Pathumwan, Bangkok 10330, Thailand 22Departamento de Física y Astronomía, Facultad de Ciencias, Universidad de La Serena, Cister- 5Institute for Cosmology and Gravitation, University of Portsmouth, Dennis Sciama Building, Burn- nas 1200, La Serena, Chile aby Road, Portsmouth PO1 3FX 23Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, 6Department of Physics and Astronomy, University of Kentucky, 505 Rose Street, Lexington, KY PA 16802 40506-0055, USA 24Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 7Department of Astronomy, University of Wisconsin-Madison, 475 North Charter Street, Madison, 16802 1 41 −1 locity gradient field across the offset absorption (Fig. 1e) that ranges the AGN’s mechanical output (Pmech = 8.1 10 erg s ) is sufficient −1 −1 ˙ × 39 −1 from approximately vNa D = −80 km s to vNa D = 60 km s . to supply the wind’s kinetic power (Ewind 10 erg s ; see Methods). ∼ These observations indicate the presence of multiple gas compo- Moreover, the wind can inject sufficient energy, coupled to the ambient nents with different temperatures and velocity structures. We inter- gas through the turbulent dynamics observed (Fig. 1i and Fig. 3a-c), to ˙ 39 −1 pret the ionized gas velocity field as resulting from a centrally driven balance cooling in both the ionized and cool gas (Egas 10 erg s ). 8∼ (volume-filling) wind with a wide opening angle. The projected flux Indeed, the amount of cool Na D gas (Mcool gas 10 M⊙) implies a −2 −1 ∼ distribution of this ionized component largely follows the stellar sur- star formation rate of SFR 1 10 M⊙ yr , which is much higher 19 ∼ × −5 −1 face brightness, suggesting that its primary ionization source is the than the estimated SFRAkira = 7 10 M⊙ yr that leverages well- × local radiation field from evolved stars20–22. The bisymmetric EW fea- detected WISE photometry. The picture that emerges is one in which tures represent enhanced emission due to shocks or over-densities along cool gas inflow in Akira, triggered by the minor merger, has initiated a the wind’s central axis. A distinct and cooler gas component is indi- relatively low-power AGN-driven wind that is nonetheless able to heat cated by the Na D absorption. Because it is spatially confined with the surrounding gas through turbulence and shocks and thereby prevent its own velocity structure, this cooler foreground material is likely to any substantial star formation. be within 1–2 Re of Akira. Simulations of galaxy mergers constrained As with Akira, the other galaxies in this class show little or no by the data (see Methods) suggest that the cool component is part of ongoing star formation, and the majority harbor similarly weak radio a tidal stream and is arcing towards the observer from the far West point sources (according to followup Jansky VLA observations) that (blueshifted; Fig. 1e) before plunging back towards Akira’s center (red- would be classified as “jet mode,” “kinetic mode,” or “radio mode” shifted; Fig. 1e). AGN7,15. With similar levels of fast-moving ionized gas oriented along Previous work has noted similar objects20,23–25 but has typically at- enhanced ionized emission, we conclude that AGN-driven winds are tributed their gaseous dynamics and unusual emission line features to present in these systems as well and represent an important heating accreted, rotating disks26.