Volume 592, August 2016
|Number of page(s)||14|
|Section||Stellar structure and evolution|
|Published online||04 July 2016|
Chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core-collapse calculations⋆
1 École normale supérieure de Lyon, CRAL, UMR CNRS 5574, 69364 Lyon Cedex 07, France
2 School of Physics, University of Exeter, Exeter, EX4 4QL, UK
3 Laboratoire AIM, Paris-Saclay, CEA/IRFU/SAp-CNRS-Université Paris Diderot, 91191 Gif-sur-Yvette Cedex, France
4 Centre for Star and Planet Formation, Niels Bohr Institute and Natural History Museum of Denmark, University of Copenhagen, Øster Voldgade 5–7, 1350 Copenhagen K, Denmark
Received: 18 June 2015
Accepted: 19 April 2016
We develop a detailed chemical network relevant to calculate the conditions that are characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of potassium, sodium, and hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to nH = 1012 cm-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.
Key words: magnetohydrodynamics (MHD) / ISM: molecules / stars: formation
The multi-dimensional multi-species equilibrium abundance table and a copy of the code are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (188.8.131.52) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/592/A18
© ESO, 2016
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.