2021MNRAS.501.4136A

Gas metallicity (Z) and the related dust-to-gas ratio (f_d) can influence the growth of H II regions via metal line cooling and ultraviolet (UV) absorption. We model these effects in star-forming regions containing massive stars. We compute stellar feedback from photoionization and radiation pressure (RP) using Monte Carlo radiative transfer coupled with hydrodynamics, including stellar and diffuse radiation fields. We follow a 10⁵ M_☉ turbulent cloud with Z/Z_☉ = 2, 1, 0.5, and 0.1, and f_\textrmd = 0.01 Z/Z_☉ with a cluster-sink particle method for star formation. The models evolve for at least 1.5 Myr under feedback. Lower Z results in higher temperatures and therefore larger H II regions. For Z >= Z_☉, RP (P_rad) can dominate locally over the gas pressure (P_gas) in the inner half-parsec around sink particles. Globally, the ratio of P_rad/P_gas is around 1 (2 Z_☉), 0.3 (Z_☉), 0.1 (0.5 Z_☉), and 0.03 (0.1 Z_☉). In the solar model, excluding RP results in an ionized volume several times smaller than the fiducial model with both mechanisms. Excluding RP and UV attenuation by dust results in a larger ionized volume than the fiducial case. That is, UV absorption hinders growth more than RP helps it. The radial expansion velocity of ionized gas reaches +15 km s^–1 outwards, while neutral gas has inward velocities for most of the runtime, except for 0.1 Z_☉ that exceeds +4 km s^–1. Z and f_d do not significantly alter the star formation efficiency, rate, or cluster half-mass radius, with the exception of 0.1 Z_☉ due to the earlier expulsion of neutral gas.