A grid of one-dimensional low-mass star formation collapse models

Context. Numerical simulations of star formation are becoming ever more sophisticated, incorporating new physical processes in increasingly realistic set-ups. These models are being compared to the latest observations through state-of-the-art synthetic renderings that trace the different chemical species present in the protostellar systems. The chemical evolution of the interstellar and protostellar matter is very topical, with more and more chemical databases and reaction solvers available online to the community. Aims. The current study was developed to provide a database of relatively simple numerical simulations of protostellar collapse as a template library for observations of cores and very young protostars, and for researchers who wish to test their chemical modelling under dynamic astrophysical conditions. It was also designed to identify statistical trends that may appear when running many models of the formation of low-mass stars by varying the initial conditions. Methods. A large set of 143 calculations of the gravitational collapse of an isolated sphere of gas with uniform temperature and a Bonnor-Ebert-like density profile was undertaken using a 1D fully implicit Lagrangian radiation hydrodynamics code. The parameter space covered initial masses from 0.2 to 8 M ⊙ , temperatures of 5−30 K, and radii 3000 ≤ R 0 ≤ 30 000 AU. Results. A spread due to differing initial conditions and optical depths, was found in the thermal evolutionary tracks of the runs. Within less than an order of magnitude, all first and second Larson cores had masses and radii essentially independent of the initial conditions. Radial profiles of the gas density, velocity, and temperature were found to vary much more outside of the first core than inside. The time elapsed between the formation of the first and second cores was found to strongly depend on the first core mass accretion rate, and no first core in our grid of models lived for longer than 2000 years before the onset of the second collapse. Conclusions. The end product of a protostellar cloud collapse, the second Larson core, is at birth a canonical object with a mass and radius of about 3 M J and 8 R J , independent of its initial conditions. The evolution sequence which brings the gas to stellar densities can, however, proceed in a variety of scenarios, on different timescales or along different isentropes, but each story line can largely be predicted by the initial conditions. All the data from the simulations are publicly available.

[1]  T. Henning,et al.  Synthetic observations of first hydrostatic cores in collapsing low-mass dense cores. I. Spectral energy distributions and evolutionary sequence , 2012, 1207.0656.

[2]  Tucson,et al.  The AGORA High-resolution Galaxy Simulations Comparison Project. III. Cosmological Zoom-in Simulation of a Milky Way–mass Halo , 2013, The Astrophysical Journal.

[3]  T. Henning,et al.  Rosseland and Planck mean opacities for protoplanetary discs , 2003, astro-ph/0308344.

[4]  R. Klessen,et al.  Disc formation in turbulent massive cores: circumventing the magnetic braking catastrophe , 2012, 1201.5302.

[5]  P. Hennebelle,et al.  FRAGMENTATION OF MASSIVE DENSE CORES DOWN TO ≲ 1000 AU: RELATION BETWEEN FRAGMENTATION AND DENSITY STRUCTURE , 2014, 1401.8292.

[6]  P. Hennebelle,et al.  Ambipolar diffusion in low-mass star formation - I. General comparison with the ideal magnetohydrodynamic case , 2015, 1509.05630.

[7]  C. Loore,et al.  The Formation of Stars , 1992 .

[8]  T. Haugboelle,et al.  A detailed framework to incorporate dust in hydrodynamical simulations , 2016, 1606.01229.

[9]  P. Padoan,et al.  The Stellar Initial Mass Function from Turbulent Fragmentation , 2000, astro-ph/0011465.

[10]  N. Peretto,et al.  A census of dense cores in the Aquila cloud complex: SPIRE/PACS observations from the Herschel Gould Belt survey , 2015, 1507.05926.

[11]  S. Okuzumi,et al.  RADIATION MAGNETOHYDRODYNAMIC SIMULATIONS OF PROTOSTELLAR COLLAPSE: PROTOSTELLAR CORE FORMATION , 2012, 1206.3567.

[12]  R. Teyssier,et al.  Protostellar collapse: radiative and magnetic feedbacks on small-scale fragmentation , 2009, Astronomy and Astrophysics.

[13]  A. Frank,et al.  Triggered star formation and its consequences , 2014, 1404.5924.

[14]  R. Klessen,et al.  MODELING JET AND OUTFLOW FEEDBACK DURING STAR CLUSTER FORMATION , 2014, 1406.3625.

[15]  C. J. Zeippen,et al.  Updated opacities from the Opacity Project , 2004, astro-ph/0410744.

[16]  J. Jørgensen,et al.  Strong irradiation of protostellar cores in Corona Australis , 2012, 1209.4817.

[17]  Hirohiko Masunaga,et al.  A Radiation Hydrodynamic Model for Protostellar Collapse. I. The First Collapse , 1998 .

[18]  J. Ferguson,et al.  Simulations of protostellar collapse using multigroup radiation hydrodynamics , 2013, Astronomy & Astrophysics.

[19]  F. Gianturco,et al.  KROME - a package to embed chemistry in astrophysical simulations , 2013, 1311.1070.

[20]  Shu-ichiro Inutsuka,et al.  A Radiation Hydrodynamic Model for Protostellar Collapse. II. The Second Collapse and the Birth of a Protostar , 2000 .

[21]  E. Ostriker,et al.  Theory of Star Formation , 2007, 0707.3514.

[22]  C. Hunter The collapse of unstable isothermal spheres. , 1977 .

[23]  E. Salpeter,et al.  Primordial stellar evolution - The protostar phase , 1986 .

[24]  P. Padoan,et al.  INFALL-DRIVEN PROTOSTELLAR ACCRETION AND THE SOLUTION TO THE LUMINOSITY PROBLEM , 2014, 1407.1445.

[25]  S. Inutsuka,et al.  FIRST DIRECT SIMULATION OF BROWN DWARF FORMATION IN A COMPACT CLOUD CORE , 2009, 0907.3255.

[26]  T. Millar,et al.  The UMIST database for astrochemistry 2012 , 2012, 1212.6362.

[27]  S. Frimann,et al.  Large-scale numerical simulations of star formation put to the test: Comparing synthetic images and actual observations for statistical samples of protostars , 2015, 1510.07827.

[28]  S. Frimann,et al.  Protostellar accretion traced with chemistry: Comparing synthetic C18O maps of embedded protostars to real observations , 2015, 1512.00416.

[29]  K. Dobashi,et al.  STAR FORMATION IN TURBULENT MOLECULAR CLOUDS WITH COLLIDING FLOW , 2014, 1412.5741.

[30]  B. Dubroca,et al.  Etude théorique et numérique d'une hiérarchie de modèles aux moments pour le transfert radiatif , 1999 .

[31]  Christoph Federrath,et al.  Inefficient star formation through turbulence, magnetic fields and feedback , 2015, 1504.03690.

[32]  R. Klein,et al.  RADIATION-HYDRODYNAMIC SIMULATIONS OF THE FORMATION OF ORION-LIKE STAR CLUSTERS. II. THE INITIAL MASS FUNCTION FROM WINDS, TURBULENCE, AND RADIATION , 2012, 1203.2620.

[33]  Shu-ichiro Inutsuka,et al.  Does “τ≈1” Terminate the Isothermal Evolution of Collapsing Clouds? , 1999 .

[34]  Gilles Chabrier,et al.  An Equation of State for Low-Mass Stars and Giant Planets , 1995 .

[35]  P. Hopkins,et al.  The fundamentally different dynamics of dust and gas in molecular clouds , 2015, 1510.02477.

[36]  David R. Alexander,et al.  Low-Temperature Opacities , 2005, astro-ph/0502045.

[37]  M. Penston Dynamics of Self-Gravitating Gaseous Spheres—III: Analytical Results in the Free-fall of Isothermal Cases , 1969 .

[38]  Richard B. Larson,et al.  Numerical Calculations of the Dynamics of a Collapsing Proto-Star , 1969 .

[39]  N. Vaytet,et al.  Multigroup radiation hydrodynamics with flux-limited diffusion and adaptive mesh refinement , 2015, 1504.01894.

[40]  J. Ferreira Braking down an accreting protostar: disc-locking, disc winds, stellar winds, X-winds and Magnetospheric Ejecta , 2013, 1310.4203.

[41]  C. D. Levermore,et al.  Relating Eddington factors to flux limiters , 1984 .

[42]  S. Okuzumi,et al.  RADIATION MAGNETOHYDRODYNAMIC SIMULATIONS OF PROTOSTELLAR COLLAPSE: NONIDEAL MAGNETOHYDRODYNAMIC EFFECTS AND EARLY FORMATION OF CIRCUMSTELLAR DISKS , 2015, 1501.04102.

[43]  Pierre Charrier,et al.  An HLLC Scheme to Solve The M1 Model of Radiative Transfer in Two Space Dimensions , 2007, J. Sci. Comput..

[44]  P. Foster,et al.  Gravitational collapse of an isothermal sphere , 1992 .

[45]  F. Shu Self-similar collapse of isothermal spheres and star formation. , 1977 .

[46]  H. Shibai,et al.  MILLIMETER-WAVE POLARIZATION OF PROTOPLANETARY DISKS DUE TO DUST SCATTERING , 2015, 1504.04812.

[47]  P. Hennebelle,et al.  COLLAPSE OF MASSIVE MAGNETIZED DENSE CORES USING RADIATION MAGNETOHYDRODYNAMICS: EARLY FRAGMENTATION INHIBITION , 2011, 1110.2955.

[48]  G. Chabrier,et al.  The influence of frequency-dependent radiative transfer on the structures of radiative shocks , 2013, 1303.1385.

[49]  Moustafa T. Chahine,et al.  Foundations of Radiation Hydrodynamics (Dimitri Mihalas and Barbara Weibel Mihalas) , 1987 .

[50]  P. K. Leung,et al.  Self-similar fragmentation regulated by magnetic fields in a region forming massive stars , 2015, Nature.

[51]  Daniel J. Price,et al.  Can non-ideal magnetohydrodynamics solve the magnetic braking catastrophe? , 2015, 1512.01597.

[52]  U. Michigan,et al.  Chemical tracers of episodic accretion in low-mass protostars , 2015, 1503.04951.

[53]  Á. Sánchez-Monge,et al.  Revealing the dynamics of Class 0 protostellar discs with ALMA , 2016, 1601.02384.

[54]  U. Exeter,et al.  Simulations of protostellar collapse using multigroup radiation hydrodynamics I. The first collapse , 2012, 1205.5143.

[55]  V. Wakelam,et al.  CHEMICAL AND PHYSICAL CHARACTERIZATION OF COLLAPSING LOW-MASS PRESTELLAR DENSE CORES , 2016, 1603.02529.

[56]  G. Nyman,et al.  THE 2014 KIDA NETWORK FOR INTERSTELLAR CHEMISTRY , 2015, 1503.01594.

[57]  W. Bonnar,et al.  Boyle's Law and gravitational instability , 1956 .

[58]  Daniel J. Price,et al.  Collapse of a molecular cloud core to stellar densities: stellar-core and outflow formation in radiation magnetohydrodynamic simulations , 2013, 1310.1092.