Toward Resolution Convergence for Star Formation and Galactic Winds in Galaxy Simulations
Abstract
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 ∼20 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 (Rcool) 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 ∼20%, while relative errors in the cool-phase mass outflow rate and integrated X-ray luminosity decrease to below ∼40%. The warm and hot outflow phases remain less well converged, with variations of ∼50% and ∼75%, respectively. The upgraded model also improves the convergence of normalized wind properties, with relative errors below ∼50%, ∼80%, and ∼20% 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.
Related
- Galactic Morpholography: Linking present-day galaxy compactness to assembly history and relic formation in the COLIBRE simulations
- Comment on "Fermi-LAT Galactic Center Excess Morphology of Dark Matter in Simulations of the Milky Way Galaxy"
- From short-lived to long-lived clouds: Impact of star formation models on giant molecular cloud evolution in simulations of an NGC 300-like galaxy
- Comparing Infrared & X-ray-selected Active Galactic Nuclei host galaxy star formation rates to the galaxy main sequence
- Comparing infrared and X-ray-selected active galactic nuclei host galaxy star formation rates to the galaxy main sequence
BibTeX
@misc{li2026toward,
title = {{Toward Resolution Convergence for Star Formation and Galactic Winds in Galaxy Simulations}},
author = {Xue-Fu Li and Weishan Zhu and Tian-Rui Wang and Huan Chen and Long-Long Feng},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.34336},
eprint = {2609.34336},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.34336},
}
RIS
TY - UNPB
TI - Toward Resolution Convergence for Star Formation and Galactic Winds in Galaxy Simulations
AU - Xue-Fu Li
AU - Weishan Zhu
AU - Tian-Rui Wang
AU - Huan Chen
AU - Long-Long Feng
PY - 2026
DA - 2026/09/28
PB - arXiv
DO - 10.48550/arxiv.2609.34336
UR - https://arxiv.org/abs/2609.34336
AB - 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 ∼20 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_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 ∼20%, while relative errors in the cool-phase mass outflow rate and integrated X-ray luminosity decrease to below ∼40%. The warm and hot outflow phases remain less well converged, with variations of ∼50% and ∼75%, respectively. The upgraded model also improves the convergence of normalized wind properties, with relative errors below ∼50%, ∼80%, and ∼20% 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.
ER -