The birth and early evolution of a low mass protostar

Abstract: Understanding the collapse of dense molecular cloud cores to stellar densities and the subsequent evolution of the protostar is of importance to model the feedback effects such an object has on its surrounding environment, as well as describing the conditions with which it enters the stellar evolutionary track. This process is fundamentally multi-scale and necessitates the use of robust numerical simulations. We aim to model the birth and early evolution of a low-mass protostar. We also seek to describe the interior structure of the protostar and the radiative behavior of its accretion shock front. We carry out a high resolution numerical simulation of the collapse of a gravitationally unstable $1$ $\mathrm{M_{\odot}}$ dense molecular cloud core to stellar densities using three-dimensional radiation hydrodynamics under the gray flux-limited diffusion approximation. We follow the initial isothermal phase, the first adiabatic contraction, the second gravitational collapse triggered by the dissociation of $\mathrm{H}_{2}$ molecules, and $\approx 247$ days of the subsequent main accretion phase. We find that the sub-critical radiative behavior of the protostar's shock front causes it to swell as it accretes matter. We also find that the protostar is turbulent from the moment of its inception despite its radiative stability. This turbulence causes significant entropy mixing inside the protostar, which regulates the swelling. Furthermore, we find that the protostar is not fully ionized at birth, but the relative amount of ionized material within it increases as it accretes matter from its surroundings. Finally, we report in the appendix the results of the first 3D calculations involving a frequency-dependent treatment of radiative transfer, which has not produced any major differences with its gray counterpart.