Simulating Radiative Feedback and the Formation of Massive Stars

Abstract

This thesis is a study of massive star formation: the environments in which they form and the effect that their radiation feedback has on their environments. We present high-performance supercomputer simulations of massive star formation inside molecular cloud clumps and cores. First, we present a novel radiative transfer code that hybridizes two previous approaches to radiative transfer (raytracing and flux-limited diffusion) and implements it in a Cartesian grid-based code with adaptive mesh refinement, representing the first of such implementations. This hybrid radiative transfer code allows for more accurate calculations of the radiation pressure and irradiated gas temperature that are the hallmark of massive star formation and which threaten to limit the mass which stars can ultimately obtain. Next, we apply this hybrid radiative transfer code in simulations of massive protostellar cores. We simulate their gravitational collapse and the formation of a massive protostar surrounded by a Keplerian accretion disk. These disks become gravitationally unstable, increasing the accretion rate onto the star, but do not fragment to form additional stars. We demonstrate that massive stars accrete material predominantly through their circumstellar disks, and via radiation pressure drive large outflow bubbles that appear stable to classic fluid instabilities. Finally, we present simulations of the larger context of star formation: turbulent, magnetised, filamentary cloud clumps. We study the magnetic field geometry and accretion flows. We find that in clouds where the turbulent and magnetic energies are approximately equal, the gravitational energy must dominate the kinetic energy for there to be a coherent magnetic field structure. Star cluster formation takes place inside the primary filament and the photoionisation feedback from a single massive star drives the creation of a bubble of hot, ionised gas that ultimately engulfs the star cluster and destroys the filament.