Sensitivity analyses of dense cloud chemical models*
V. Wakelam1,2, E. Herbst3, J. Le Bourlot4, F. Hersant1,2, F. Selsis1,2 and S. Guilloteau1,2
Université de Bordeaux, Observatoire Aquitain des Sciences de l'Univers, BP 89, 33271 Floirac Cedex, France e-mail: firstname.lastname@example.org
2 CNRS, UMR 5804, Laboratoire d'Astrophysique de Bordeaux, BP 89, 33271 Floirac Cedex, France
3 Departments of Physics, Astronomy, and Chemistry, The Ohio State University, Columbus, OH 43210, USA
4 Observatoire de Paris, LUTH and Université Paris-Diderot, Place J. Janssen, 92190 Meudon, France
Accepted: 2 April 2010
Context. Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters.
Aims. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances.
Methods. We used a Monte Carlo method to randomly vary individual parameters and groups of parameters within realistic ranges. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively.
Results. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4 × 105 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature.
Conclusions. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints.
Key words: astrochemistry / molecular processes / ISM: abundances / ISM: molecules / ISM: individual objects: L134N / ISM: individual objects: TMC-1 (CP)
The osu.03.2008 network is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (22.214.171.124) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/517/A21
© ESO, 2010