MNRAS 449, 1057–1075 (2015) doi:10.1093/mnras/stv324
Modelling the supernova-driven ISM in different environments
A. Gatto,1‹ S. Walch,1,2 M.-M. Mac Low,3 T. Naab,1 P. Girichidis,1 S. C. O. Glover,4 R. Wunsch,¨ 5 R. S. Klessen,4,6,7 P. C. Clark,4,8 C. Baczynski,3,4 T. Peters,1,9 J. P. Ostriker,10,11 J. C. Iba´nez-Mej˜ ´ıa3,4 and S. Haid2 1Max-Planck-Institut fur¨ Astrophysik, Karl-Schwarzschild-Strasse 1, D-85748 Garching, Germany 2Physikalisches Institut, Universitat¨ Koln,¨ Zulpicher¨ Strasse 77, D-50937 Koln,¨ Germany 3Department of Astrophysics, American Museum of Natural History, 79th Street at Central Park West, New York, NY 10024, USA 4Institut fur¨ Theoretische Astrophysik, Zentrum fur¨ Astronomie, Universitat¨ Heidelberg, Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany 5Astronomicky´ Ustav,´ Akademie ved˘ Cesk˘ e´ Republiky, Bocn˘ ´ı II 1401, CZ-14131 Praha, Czech Republic 6Department of Astronomy and Astrophysics, University of California, 1156 High Street, Santa Cruz, CA 95064, USA 7Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, SLAC Nat. Acc. Lab., Menlo Park, CA 94025, USA 8School of Physics and Astronomy, Cardiff University, 5 The Parade, Cardiff CF24 3AA, UK 9Institut fur¨ Computergestutzte¨ Wissenschaften, Universitat¨ Zurich,¨ Winterthurerstrasse 190, CH-8057 Zurich,¨ Switzerland 10Department of Astrophysical Sciences, Princeton University, 4 Ivy Ln, Princeton, NJ 08544, USA 11Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY 10027, USA
Accepted 2015 January 28. Received 2015 January 21; in original form 2014 October 31
ABSTRACT We use hydrodynamical simulations in a (256 pc)3 periodic box to model the impact of supernova (SN) explosions on the multiphase interstellar medium (ISM) for initial densities n = 0.5–30 cm−3 and SN rates 1–720 Myr−1. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas (T 106 K) filling >90 per cent of the volume. This gas reaches high 4 7 −3 pressures (10 < P/kB < 10 Kcm ) due to the combination of SN explosions in the hot, low-density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass (>50 per cent), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ( 1 per cent). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity −1 dispersion in H I remains 10 km s in all cases. For peak driving, the velocity dispersion in Hα can be as high as 70 km s−1 due to the contribution from young, embedded SN remnants. Key words: methods: numerical – ISM: evolution – ISM: kinematics and dynamics – ISM: structure – ISM: supernova remnants – galaxies: evolution.
motions in the ISM gas (Norman & Ferrara 1996; Joung & Mac 1 INTRODUCTION Low 2006; de Avillez & Breitschwerdt 2007; Gent et al. 2013), Supernova (SN) explosions are an important component for shaping although their relative importance with respect to other driving the interstellar medium (ISM). In particular, they produce its hottest mechanisms has yet to be fully determined (see e.g. Mestel & phase, with temperature T 106 K (Cox & Smith 1974; McKee & Spitzer 1956; Sellwood & Balbus 1999; Kritsuk & Norman 2002; Ostriker 1977). In the Milky Way, this hot phase fills from 20 to Wada, Meurer & Norman 2002; Tamburro et al. 2009; Klessen & 80 per cent of the volume with increasing height above the disc Hennebelle 2010). (Ferriere` 2001; Kalberla & Dedes 2008). SNe also drive turbulent Supersonic turbulent motions in warm and cold gas, observed in emission from H I, have been found in many local and extragalactic environments, with a typical one-dimensional velocity dispersion of −1 E-mail: [email protected] σ 10 km s (e.g. Heiles & Troland 2003; Petric & Rupen 2007;