Track the frontier, from condensed matter to cosmology.
Search papers, labs, and topics across Lattice.
The authors simulated an M82-like starburst galaxy across 4 to 32 pc resolutions to resolve the persistent numerical divergence of sub-grid star formation and feedback models in hydrodynamic galaxy simulations. By anchoring sink-particle accretion to invariant ~20 pc giant molecular cloud scales and coupling supernova energy directly to the remnant's physical cooling radius, the prescriptions isolate physical processes from grid cell size. This framework successfully caps cross-resolution variance in total stellar mass and kinetic energy at ~20%, while reducing relative errors in cool-phase mass outflow rates to under 40% across an 8× resolution range.
Anchoring sub-grid star formation and supernova feedback to invariant physical scales rather than numerical cell sizes cuts cross-resolution variance in simulated stellar mass down to ~20% across an 8× grid refinement.
The predictive power of hydrodynamic simulations of galaxies is often hampered by the resolution dependence of sub-grid models for star formation (SF) and stellar feedback. We present three-dimensional hydrodynamic simulations of an M82-like low-mass starburst galaxy at resolutions of 4, 8, 16, 32 pc and develop optimized sub-grid prescriptions to mitigate these numerical inconsistencies. For star formation, we adopt a $\sim20$ pc physical-scale control volume, corresponding to the typical size of giant molecular clouds (GMCs), for gas accretion onto sink particles, reducing its dependence on numerical resolution. We further incorporate an empirical, resolution- and density-dependent star formation scheme that accounts for variations in gas surface density across spatial scales, as suggested by observations. Furthermore, we refine the supernova feedback coupling by anchoring the injection volume to the physical cooling radius of the supernova remnant ($R_{\rm cool}$) and utilizing a resolution-dependent thermal-to-kinetic energy partition. These improvements substantially enhance convergence across 4 to 32 pc: variations in the total stellar mass formed and kinetic energy are reduced to $\sim20\%$, while relative errors in the cool-phase mass outflow rate and integrated X-ray luminosity decrease to below $\sim40\%$. The warm and hot outflow phases remain less well converged, with variations of $\sim50\%$ and $\sim75\%$, respectively. The upgraded model also improves the convergence of normalized wind properties, with relative errors below $\sim50\%$, $\sim80\%$, and $\sim20\%$ for the time-averaged mass loading factor, energy loading factor, and X-ray emission efficiency, respectively. Our results provide a practical framework for achieving more robust simulations of starburst galaxies.