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

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ömgren 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 ($\sim 8000$ 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.

19 pages, 14 figures; MNRAS published

References in corpus (19)

• Array Programming with NumPy
• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• Dust Masses, PAH Abundances, and Starlight Intensities in the SINGS Galaxy Sample
• Fast and inefficient star formation due to short-lived molecular clouds and rapid feedback
• Before the first supernova: combined effects of HII regions and winds on molecular clouds
• A radiative transfer scheme for cosmological reionization based on a local Eddington tensor
• A Detailed Study of Feedback from a Massive Star
• Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies
• The GRIFFIN project -- Formation of star clusters with individual massive stars in a simulated dwarf galaxy starburst
• LyC escape from SPHINX galaxies in the Epoch of Reionization
• Numerical Simulations of Turbulent, Molecular Clouds Regulated by Radiation Feedback Forces I: Star Formation Rate and Efficiency
• Disruption of giant molecular clouds and formation of bound star clusters under the influence of momentum stellar feedback
• An implementation of radiative transfer in the cosmological simulation code GADGET
• H-alpha emission in local galaxies: star formation, time variability and the diffuse ionized gas
• Momentum feedback from marginally-resolved HII regions in isolated disc galaxies
• Structure and Expansion Law of HII Regions in structured Molecular Clouds
• Theoretically Modelling Photoionized Regions with Fractal Geometry in Three Dimensions
• Modelling of ionising feedback with Smoothed Particle Hydrodynamics and Monte Carlo Radiative Transfer on a Voronoi grid
• How D-type HII region expansion depends on numerical resolution