Vertical Equilibrium, Energetics, and Star Formation Rates in Magnetized Galactic Disks Regulated by Momentum Feedback from Supernovae

TL;DR

This study uses 3D magnetohydrodynamic simulations to explore how magnetic fields influence the equilibrium state, star formation rates, and pressure support in galactic disks, revealing magnetic effects can reduce star formation by about 25%.

Contribution

It extends previous hydrodynamic models to include magnetic fields, demonstrating that magnetic energy saturates similarly to kinetic and thermal energies, and quantifies magnetic feedback effects on star formation.

Findings

Magnetic fields reach saturation quickly, similar to other energies.

Magnetic energy proportions are consistent across a wide range of initial conditions.

Magnetic feedback reduces the predicted star formation rate by approximately 25%.

Abstract

Recent hydrodynamic (HD) simulations have shown that galactic disks evolve to reach well-defined statistical equilibrium states. The star formation rate (SFR) self-regulates until energy injection by star formation feedback balances dissipation and cooling in the interstellar medium (ISM), and provides vertical pressure support to balance gravity. In this paper, we extend our previous models to allow for a range of initial magnetic field strengths and configurations, utilizing three-dimensional, magnetohydrodynamic (MHD) simulations. We show that a quasi-steady equilibrium state is established as rapidly for MHD as for HD models unless the initial magnetic field is very strong or very weak, which requires more time to reach saturation. Remarkably, models with initial magnetic energy varying by two orders of magnitude approach the same asymptotic state. In the fully saturated state of the fiducial model, the integrated energy proportions E_kin:E_th:E_mag,t:E_mag,o are 0.35:0.39:0.15:0.11, while the proportions of midplane support P_turb:P_th:\Pi_mag,t:\Pi_mag,o are 0.49:0.18:0.18:0.15. Vertical profiles of total effective pressure satisfy vertical dynamical equilibrium with the total gas weight at all heights. We measure the "feedback yields" \eta_c=P_c/\Sigma_SFR (in suitable units) for each pressure component, finding that \eta_turb~4 and \eta_th~1 are the same for MHD as in previous HD simulations, and \eta_mag,t~1. These yields can be used to predict the equilibrium SFR for a local region in a galaxy based on its observed gas and stellar surface densities and velocity dispersions. As the ISM weight (or dynamical equilibrium pressure) is fixed, an increase in $\eta$ from turbulent magnetic fields reduces the predicted \Sigma_SFR by ~25% relative to the HD case.