Simulations of protostellar collapse using multigroup radiation hydrodynamics

Neil Vaytet (1) , Gilles Chabrier (2, 1) , Edouard Audit (3, 4) , Benoît Commerçon (5) , Jacques Masson (1) , Jason Ferguson (6) , Franck Delahaye (5)

1 ENS de Lyon - École normale supérieure de Lyon
2 EXETER - School of Physics
3 AIM (UMR7158 / UMR_E_9005 / UM_112) - Astrophysique Interprétation Modélisation
4 MDLS - Maison de la Simulation
5 LERMA - Laboratoire d'Etude du Rayonnement et de la Matière en Astrophysique
6 Wichita State University

Context. Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite, and the protostar is born. This sequence of events entails many physical processes, of which radiative transfer is of paramount importance. Simulated collapsing cores without radiative transfer rapidly become thermally supported before reaching high enough temperatures and densities, preventing the formation of stars.Aims. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency, which can potentially be overlooked by grey models and lead to significantly different results. In this paper, we follow up on a previous paper on the collapse and formation of Larson’s first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson’s second core.Methods. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1, and 10 M⊙ cloud cores.Results. We find that, for an identical central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass, and temperature of the second cores are independent of initial cloud mass, size, and temperature.Conclusions. Our simulations support the idea of a standard (universal) initial second core size of ~3 × 10-3 AU and mass ~1.4 × 10-3 M⊙. The grey approximation for radiative transfer appears to perform well in one-dimensional simulations of protostellar collapse, most probably because of the high optical thickness of the majority of the protostar-envelope system. A simple estimate of the characteristic timescale of the second core suggests that the effects of using multigroup radiative transfer may be more important in the long-term evolution of the protostar.

Keywords: stars: formation; methods: numerical; hydrodynamics; radiative transfer; stars: protostars

Domains: Astrophysics [astro-ph]

Neil Vaytet, Gilles Chabrier, Edouard Audit, Benoît Commerçon, Jacques Masson, et al.. Simulations of protostellar collapse using multigroup radiation hydrodynamics. Astronomy and Astrophysics - A&A, 2013, 557, pp.A90. ⟨10.1051/0004-6361/201321423⟩. ⟨cea-01135422⟩