arXiv preprint ·

Toward Resolution Convergence for Star Formation and Galactic Winds in Galaxy Simulations

Xue-Fu Li, Weishan Zhu, Tian-Rui Wang, Huan Chen and Long-Long Feng

Published
Source
arXiv preprint, posted before peer review
Original
Read it at arXiv
Topics
Astrophysics of GalaxiesSolar & Stellar Astrophysics

In brief

Boosting Galaxy Simulations' Accuracy

This paper introduces refined sub-grid models and optimized simulations for low-mass starburst galaxies, significantly improving convergence across various resolutions. The authors present a practical framework that enhances the reliability of simulations, making them more robust for future research.

What they found: The study demonstrates that by optimizing sub-grid models, the authors achieved significant improvements in the convergence of galaxy simulations. Variations in total stellar mass formed and kinetic energy were reduced to around 20%, while errors in cool-phase mass outflow rates dropped below 40%. Additionally, normalized wind properties showed substantial enhancements, with relative errors now below 50%.

AI-written summary. Check anything you rely on against the original paper.

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 $\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.

The authors’ abstract, as arXiv lists it. arXiv’s listings, abstracts included, are free for anyone to reuse.

Read it on ReadSearch

Each opens after you create an account, and brings you back to this paper.

Already a reader? Sign in

All new papers in Astrophysics of Galaxies

Get the few papers worth your time in Astrophysics of Galaxies every morning.

ReadSearch reads the new papers in your field each day and brings you the few worth reading, each with its key finding in a sentence.

Papers belong to their authors and publishers; every page here links to the original. All research topics