The Effect of Radiation and Supernovae Feedback on LyC Escape in Local Star-forming Galaxies

Feedback is widely recognized as an essential condition for Lyman continuum (LyC) escape in star-forming galaxies. However, the mechanisms by which galactic outflows clear neutral gas and dust remain unclear. In this paper, we model the Mg II 2796\r{A}, 2804\r{A} absorption + emission lines in 29 galaxies taken from the Low-z LyC Survey (LzLCS) to investigate the impact of (radiation + mechanical) feedback on LyC escape. Using constraints on Mg$^+$ and photoionization models, we map the outflows' neutral hydrogen content and predict $f_{esc}^{LyC}$ with a multiphase wind model. We measure mass, momentum, and energy loading factors for the neutral winds, which carry up to 10% of the momentum and 1% of the energy in SFR-based deposition rates. We use SED template fitting to determine the relative ages of stellar populations, allowing us to identify radiation feedback dominant systems. We then examine feedback related properties (stellar age, loading factors, etc.) under conditions that optimize feedback efficiency, specifically high star formation rate surface density and compact UV half-light radii. Our findings indicate that the strongest leakers are radiation feedback dominant, lack Mg II outflows, but have extended broad components in higher ionization lines like [O III] 5007\r{A}, as observed by Amor\'in et al. (2024). In contrast, galaxies experiencing supernovae feedback typically exhibit weaker $f_{esc}^{LyC}$ and show evidence of outflows in both Mg II and higher ionization lines. We attribute these findings to rapid or "catastrophic" cooling in the radiation-dominant systems, which, given the low metallicities in our sample, are likely experiencing delayed supernovae.