Simulating the Formation of Massive Protostars: I. Radiative Feedback and Accretion Disks

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities.