Simulating ionization feedback from young massive stars: impact of numerical resolution

TL;DR

This study investigates how numerical resolution affects the simulation of ionization feedback from young massive stars, emphasizing the importance of resolving the Strömgren sphere and using appropriate time steps for accurate modeling.

Contribution

The paper introduces new numerical methods and guidelines to improve the accuracy of ionization feedback simulations under limited resolution conditions.

Findings

Resolution of 10-100 elements is needed to accurately model HII regions.

Insufficient spatial resolution affects ionization fraction and feedback strength.

Temporal resolution critically influences the simulation of ionization feedback effects.

Abstract

Modelling galaxy formation in hydrodynamic simulations has increasingly adopted various radiative transfer methods to account for photoionization feedback from young massive stars. However, the evolution of HII regions around stars begins in dense star-forming clouds and spans large dynamical ranges in both space and time, posing severe challenges for numerical simulations in terms of both spatial and temporal resolution that depends strongly on gas density ($\propto n^{-1}$). In this work, we perform a series of idealized HII region simulations using the moving-mesh radiation-hydrodynamic code Arepo-RT to study the effects of numerical resolution. The simulated results match the analytical solutions and the ionization feedback converges only if the Str\"omgren sphere is resolved by at least $10$--$100$ resolution elements and the size of each time integration step is smaller than $0.1$ times the recombination timescale. Insufficient spatial resolution leads to reduced ionization fraction but enhanced ionized gas mass and momentum feedback from the HII regions, as well as degrading the multi-phase interstellar medium into a diffuse, partially ionized, warm ($\sim8000$ K) gas. On the other hand, insufficient temporal resolution strongly suppresses the effects of ionizing feedback. This is because longer timesteps are not able to resolve the rapid variation of the thermochemistry properties of the gas cells around massive stars, especially when the photon injection and thermochemistry are performed with different cadences. Finally, we provide novel numerical implementations to overcome the above issues when strict resolution requirements are not achievable in practice.