Skip to Main content Skip to Navigation
Journal articles

Chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core-collapse calculations

Abstract : 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.

The multi-dimensional multi-species equilibrium abundance table and a copy of the code are only available at the CDS via anonymous ftp to ( or via

Document type :
Journal articles
Complete list of metadata
Contributor : Nathalie POTHIER Connect in order to contact the contributor
Submitted on : Thursday, June 30, 2022 - 5:52:14 PM
Last modification on : Friday, August 5, 2022 - 3:37:49 PM


Publisher files allowed on an open archive



P. Marchand, J. Masson, G. Chabrier, P. Hennebelle, B. Commerçon, et al.. Chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core-collapse calculations. Astronomy and Astrophysics - A&A, EDP Sciences, 2016, 592, ⟨10.1051/0004-6361/201526780⟩. ⟨insu-03710577⟩



Record views


Files downloads