The formation of massive stars is a long standing problem. Although a number of theories of massive star formation exist, ideas appear to converge to a disk-mediated accretion scenario. Here we present radiative hydrodynamic simulations of a star accreting mass via a disk embedded in a torus. We use a Monte Carlo based radiation hydrodynamics code to investigate the impact that ionizing radiation has on the torus. Ionized regions in the torus midplane are found to be either gravitationally trapped or in pressure driven expansion depending on whether or not the size of the ionized region exceeds a critical radius. Trapped Hii regions in the torus plane allow accretion to progress, while expanding Hii regions disrupt the accretion torus preventing the central star from aggregating more mass, thereby setting the star's final mass. We obtain constraints for the luminosities and torus densities that lead to both scenarios.

We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionization is treated self-consistently. This leads to a recurring collapse of the HII region over time. We derive a semi-analytic model to explain this instability, and test it using spatially converged 1D radiation hydrodynamical simulations. We discuss the implications of the 1D instability on 3D radiation hydrodynamics simulations of supersonic accreting flows.

This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flows, we find gravitationally trapped HII regions with accretion through the ionisation front onto the source. We found that, for a constant ionising luminosity, fluctuations in the ionisation front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionising luminosity increases with time as the source accretes mass. The maximum radius of the recurring HII region increases on the accretion timescale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionised gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionisation front, the accretion rate onto the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of HII regions include a continuation of accretion of material onto forming stars for a period after the star starts to emit ionising radiation, and an extension of the lifetime of ultracompact HII regions.