DEPENDENCE OF INTERSTELLAR TURBULENT PRESSURE ON SUPERNOVA RATE

Feedback from massive stars is one of the least understood aspects of galaxy formation. We perform a suite of vertically stratified local interstellar medium (ISM) models in which supernova (SN) rates and vertical gas column densities are systematically varied based on the Schmidt-Kennicutt law. Our simulations have a sufficiently high spatial resolution (1.95 pc) to follow the hydrodynamic interactions among multiple SNe that structure the interstellar medium. At a given SN rate, we find that the mean mass-weighted sound speed and velocity dispersion decrease as the inverse square root of gas density. The sum of thermal and turbulent pressures is nearly constant in the midplane, so the effective equation of state is isobaric. In contrast, across our four models having SN rates that range from 1 to 512 times the Galactic SN rate, the mass-weighted velocity dispersion remains in the range 4-6 km s(-1). Hence, gas averaged over similar to 100 pc regions follows P proportional to rho(alpha) with alpha approximate to 1, indicating that the effective equation of state on this scale is close to isothermal. Simulated Hi emission lines have widths of 10-18 km s(-1), comparable to observed values. In our highest SN rate model, superbubble blowouts occur, and the turbulent pressure on large scales is greater than or similar to 4 times higher than the thermal pressure. We find a tight correlation between the thermal and turbulent pressures averaged over similar to 100 pc regions in the midplane of each model, as well as across the four ISM models. We construct a subgrid model for turbulent pressure based on analytic arguments and explicitly calibrate it against our stratified ISM simulations. The subgrid model provides a simple yet physically motivated way to include SN feedback in cosmological simulations.