Contents

A&A 466, 1197-1204 (2007)
DOI: 10.1051/0004-6361:20064981

The UMIST database for astrochemistry 2006[*]

J. Woodall1 - M. Agúndez2 - A. J. Markwick-Kemper1,3 - T. J. Millar1,4


1 - Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Sackville Street Building, Manchester M60 1QD, UK
2 - Departamento de Astrofisica Molecular e Infrarroja, Instituto de Estructura de la Materia, CSIC, Serrano 121 28006 Madrid, Spain
3 - Department of Astronomy, University of Virginia, PO Box 3818, Charlottesville, VA 22903, USA
4 - Astrophysics Research Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, UK

Received 7 February 2006 / Accepted 30 January 2007

Abstract
Aims. We present a new version of the UMIST Database for Astrochemistry, the fourth such version to be released to the public. The current version contains some 4573 binary gas-phase reactions, an increase of 10% from the previous (1999) version, among 420 species, of which 23 are new to the database.
Methods. Major updates have been made to ion-neutral reactions, neutral-neutral reactions, particularly at low temperature, and dissociative recombination reactions. We have included for the first time the interstellar chemistry of fluorine. In addition to the usual database, we have also released a reaction set in which the effects of dipole-enhanced ion-neutral rate coefficients are included.
Results. These two reactions sets have been used in a dark cloud model and the results of these models are presented and discussed briefly. The database and associated software are available on the World Wide Web at www.udfa.net.

Key words: molecular data - molecular processes - ISM: molecules - astrochemistry

   
1 Introduction

The use of astrochemical models to interpret observations of interstellar and circumstellar molecular line observations is widespread; recent years have seen the application of such models to protoplanetary disks (Aikawa & Herbst 1999; Aikawa et al. 2002; Markwick et al. 2002; Ilgner et al. 2004; Semenov et al. 2005), protoplanetary nebulae (Cernicharo 2004; Woods et al. 2005), hot molecular cores (Viti et al. 2004; Nomura & Millar 2004; Wakelam et al. 2004), the early universe (Flower & Pineau des Forêts 2003; Lepp 2003) as well as the more traditional applications to diffuse and dark interstellar clouds (Stantcheva & Herbst 2004; Lee et al. 2004), interstellar shocks (Bergin et al. 1999; Flower et al. 2003), YSOs (Rodgers & Charnley 2003; Doty et al. 2004; Stäuber et al. 2005), and circumstellar envelopes (Brown & Millar 2003; Willacy 2004). Many of these papers consider grain chemistry either implicitly or explicitly and several have been aimed at exploring deuterium fractionation (Roberts et al. 2003, 2004) and, in particular, the possible role that grains may play in the creation of the high fractionation ratios observed in many multi-deuterated species.

The role of surface chemistry in the synthesis of interstellar molecules is a matter of much debate presently and, given the difficulties in finding an appropriate numerical description of the surface kinetics, most arguments to support the influence of grain reactions start from the failure of gas-phase chemistry to reproduce abundances of particular molecules or classes of molecule. In this regard it is therefore important to have as complete a system of gas-phase reactions as possible to describe molecular synthesis.

In this paper we discuss the latest release of the UMIST Database for Astrochemistry which has undergone a detailed re-examination over the past two years. Although the number of reactions has increased by only 10% we have critically re-assessed the rate coefficients of all reactions and have incorporated the latest laboratory data. Indeed in this release, there is experimental data available on just over one-third of all reactions, an increase of about 6%. In the following section we present our conventions, Sect. 3 discusses the updates we have made to the data while in Sect. 4 we present the results of a dark cloud model.

   
2 Species and related data

He, C, N, O, F, Na, Mg, Si, P, S, Cl, and Fe, are listed in order of mass in Table 1.

The permanent electric dipole moments of the neutral molecules are listed in Table 2.

Table 3, available in full only in electronic form, contains all of the reactions and their associated rate coefficients.


   
Table 3: List of the first ten reactions in the database. The full database is available in electronic form at www.udfa.net.

Index
Type Reactants Products Alpha Beta Gamma Source $T_{\rm l}$ $T_{\rm u}$ Acc. Ref.

1
NN H CH   C H2     $1.31\times 10^{-10}$ 0.00 80.0 C 300 2000 B NIST
2 NN H CH2   CH H2     $6.64\times 10^{-11}$ 0.00 0.0 L 300 2500 A NIST
3 NN H NH   N H2     $1.73\times 10^{-11}$ 0.50 2400.0 L 80 300 C  
4 NN H CH3   CH2 H2     $1.00\times 10^{-10}$ 0.00 7600.0 L 300 2500 A NIST
5 NN H NH2   NH H2     $5.25\times 10^{-12}$ 0.79 2200.0 L 73 300 C  
6 NN H NH2   NH H2     $1.05\times 10^{-10}$ 0.00 4450.0 M 1100 3000 A NIST
7 NN H CH4   CH3 H2     $5.94\times 10^{-13}$ 3.00 4045.0 L 300 2500 A NIST
8 NN H OH   O H2     $6.99\times 10^{-14}$ 2.80 1950.0 L 300 2500 A NIST
9 NN H NH3   NH2 H2     $7.80\times 10^{-13}$ 2.40 4990.0 M 200 2500 C NIST
10 NN H H2O   OH H2     $1.59\times 10^{-11}$ 1.20 9610.0 L 250 3000 A NIST

Our reaction format is:
I, type, R1, R2, R3, P1, P2, P3, P4, $\alpha $, $\beta $, $\gamma $, source, $T_{\rm l}$, $T_{\rm u}$, accuracy, reference code where I is the reaction number; type the reaction type (see Table 7); R1, R2 and R3 are reactants (to allow for the inclusion of three-body reactions); P1, P2, P3, P4 are products; $T_{\rm l}$ ($T_{\rm u}$) is the lower (upper) temperature limit at which the rate coefficient holds; source indicates the source of the data ("M'' measured, "E'' estimated, "C'' calculated, and "L'' literature search, with the latter most often used to note rate coefficients which are derived using several experimental values). The accuracy is described by a letter - A, B, C, D, E - where the errors are <25%, <50%, within a factor of 2, within an order of magnitude, and highly uncertain, respectively. Thenumerical values of the reference code is taken from the Anicich (2004) database for ion-neutral reactions, while the other codes are listed in Table 4.

As before, some reactions are duplicated when their temperature dependence cannot be approximated by one Arrhenius-type formula. Examples include the important reaction, $\rm O + OH \longrightarrow
O_2 + H$ which is best fit by a formula which includes a (negative) activation energy barrier above the experimental lower limit of 158 K and, at low temperatures, a formula which has no activation energy barrier and, at most, only a weak temperature dependence (Smith et al. 2004; Sims 2006; Carty et al. 2006). This reaction is also an example of one which has a negative activation energy, usually the result of fitting a rate coefficient over a particular temperature range. Indiscriminate use of such rate coefficients at low temperatures, indeed outside their range of validity, can lead to serious errors in any model calculation. For those reactions for which more than one temperature range is included, we have been careful to ensure that the temperature ranges do not overlap. A more pressing concern is how to evaluate a rate coefficient outside of its given temperature range. In general, we recommend that one uses the entry in the database closest to the temperature of interest. In most cases, this will give a reasonable value for the rate coefficient. However, care needs to be taken when the rate coefficient has a negative activation energy; in such cases one should check to see whether the rate coefficient remains physical. For example, the reaction between C2H and C2H2 has a rate coefficient best fit with a negative activation energy of 25 K over the range 143-3000 K. At 10 K, it is probably still reasonable to use this fit to determine the rate coefficient. On the other hand, the reaction between O and NO2, which has a negative activation energy of 120 K over the range 200-2500 K, has an unreasonably large rate coefficient at 10 K. In our experience, the rate coefficients of the following reactions should be set to zero in a model at 10 K: reaction numbers 6, 93, 140, 171, 266, 282, 312, 316, 351, 353, 377, 419, 431, 443, 446, 448, 488, 493, 501-2, 520, 533, 540, 660, 666-7, 949, 1684, 1731, 2939, 2946, 3533-4, 3536-7, 3552-3, 3555-6, 3558-9, 3581-6, 4007-8, 4013-4, 4016, 4079, 4109, 4112, 4115, 4138, and 4558. Some of these are duplicated reactions.

In addition, we also make available a second reaction file, called Rate06-dipole, which includes the enhancement of ion-neutral rate coefficients for cases in which the neutral has a large, permanent electric dipole moment. Such enhancements, due to "dipole locking'' effects (Troe 1987; Herbst & Leung 1986b), result in rate coefficients which have a T-1/2 dependence at low temperatures. In our reaction file we include this power-law behaviour in all reactions for which (i) the neutral has a dipole moment in excess of 0.9 Debye; (ii) the reaction does not already have a temperature dependence; and (iii) the reaction does not have a measured rate coefficient at low temperatures. In general, the use of the dipole-enhanced rate file is to be preferred at the low temperatures of dark interstellar clouds. In particular, the rate coefficients are more accurate for exothermic proton-transfer reactions and many other ion-neutral reactions which proceed with unit efficiency. It is less clear that rate coefficients are enhanced when the reaction is slow at room temperature. One of the main advantages in using the dipole-enhanced rate coefficients is that they give a much better fit to the observations of protonated ions detected in the cold interstellar clouds, e.g. HCS+, HCNH+, HC3NH+, and are, as a result, to be preferred for calculations at 10 K. We discuss the influence of such coefficients on calculated abundances in Sect. 4.

   
2.1 Calculation of the rates from $\alpha $, $\beta $ and $\gamma $

For two-body reactions, the rate coefficient is given by:

\begin{displaymath}k = \alpha \; (T/300)^{\beta} \; \exp~(-\gamma /T) \qquad \mbox{cm$^3$ s$^{-1}$ }
\end{displaymath} (1)

where T is the gas temperature.

For direct cosmic-ray ionisation (R2 = CRP):

\begin{displaymath}k = \alpha \qquad \mbox{s$^{-1}$ }
\end{displaymath} (2)

whereas for cosmic-ray-induced photoreactions (R2 = CRPHOT):

\begin{displaymath}k = \alpha \; (T/300)^{\beta} \; \gamma /(1-\omega) \qquad \mbox{s$^{-1}$ }
\end{displaymath} (3)

where $\alpha $ is the cosmic-ray ionisation rate, $\gamma $ is the probability per cosmic-ray ionisation that the appropriate photoreaction takes place, and $\omega$ is the dust grain albedo in the far ultraviolet (typically 0.6 at 150 nm). The cosmic-ray ionisation rates listed here are normalized to a total ionization rate, $\zeta_0$, for H2 of $1.36\times 10^{-17}$ ionisations per sec. Rates for both cosmic-ray ionisation and cosmic-ray-induced photoreactions can be scaled to other choices of the ionization rate, $\zeta$, by multiplying by the appropriate rate coefficients by $\zeta/\zeta_0$.

For interstellar photoreactions (R2 = PHOTON), the rate is derived as:

\begin{displaymath}k = \alpha \; \exp~(-\gamma A_V) \ \mbox{s$^{-1}$ }
\end{displaymath} (4)

where $\alpha $ represents the rate in the unshielded interstellar ultraviolet radiation field, $A_{\rm V}$ is the extinction at visible wavelengths caused by interstellar dust, and $\gamma $ is the parameter used to take into account the increased extinction of dust at ultraviolet wavelengths.

   
3 Updates since RATE99


   
Table 5: New species and alterations to chemical formulae. For the latter, we specify the old representation followed by the new.

New species
F F+ HF+ HF H2F+ CF+ C3H6+ C3H7+
NH2CN NH2CNH+ N2O+ HN2O+ HCOOH+ SiF+ NCCN NCCNH+
H3C3N+ C4H7+ CH3CS+ C4N NCCNCH3+ C6H6+ C10  

Altered chemical formulae
C3H4 by CH3CCH H4C2N+ by CH3CNH+ H5C2O+ by CH3CHOH+ C2H6OH+ by CH3OCH4+
C2H7O+ by C2H5OH2+ H2C3N+ by HC3NH+ H3C3N by CH2CHCN C3H6OH+ by C2H6COH+
C5H4 by CH3C4H H3C4N+ by CH3C3N+ H3C4N by CH3C3N H4C4N+ by CH3C3NH+
C7H4 by CH3C6H H3C6N by CH3C5N H4C6N+ by CH3C5NH+ H3C8N by CH3C7N
H4C8N+ by CH3C7NH+            

Including 32 ternary reactions, the database has increased in size from 4113 in RATE99 to 4605 entries in RATE06, covering 4564 distinct chemical reactions as opposed to 4082 in RATE99. Table 5 lists newly included species and changes to nomenclature from those adopted in RATE99. There are 577 reactions new to RATE06, and of the reactions common to both versions of the database, around 450 of them have had their rate coefficient updated. Of these, the 156 cosmic-ray photon reaction rates were updated as described below, and for the remainder, the changes are either due to new data being available, or a reappraisal of the existing data. The reaction rate coefficients changing the most are largely neutral-neutral reactions whose low temperature behaviour has been reexamined in the literature, such as reactions between atomic oxygen and polyynes - rate coefficients for these reactions have typically increased from 10-25 in RATE99 to 10-11 in RATE06. There are only 10 reactions whose 10 K rate coefficient has increased by more than a factor of 100. There are 7 reactions whose 10 K rate coefficient has decreased by more than a factor of 100. These reactions and their data are shown in Table 6. On the database web site www.udfa.net, the changes to each rate coefficient between RATE95, RATE99 and RATE06 can be examined in more detail on a reaction by reaction basis.


 

 
Table 6: The reactions whose 10 K rate coefficient has changed the most between the current and previous releases of the database. In most cases this is due to new data for the reaction being available. The reference codes follow our usual convention.
  RATE06 RATE99
Reaction $\alpha $ $\beta $ $\gamma $ Reference $\alpha $ $\beta $ $\gamma $ Reference

Biggest rate coefficient decrease

C2H4 + C3H $_3^+\rightarrow$ C5H5+ + H2
1.10E-19 0.00 0.0 0001 1.10e-09 0.00 0.0 HL89
C2H3 + O $_2\rightarrow$ H2CO + HCO 4.62E-12 0.00 -171.0 NIST 2.73e-11 -1.39 510.0 NIST
OH + CO $\rightarrow$ CO2 + H 2.81E-13 0.00 176.0 SH04 1.17e-13 0.95 -74.0 NIST
H+ + H $\rightarrow$ H2+ + $h\nu$ 5.13E-19 1.85 0.0 N/A 1.00e-17 0.00 0.0 PH80
OH + OH $\rightarrow$ H2O + O 1.65E-12 1.14 50.0 SH04 3.87e-13 1.69 -469.0 NIST
H + CH $\rightarrow$ C + H2 1.31E-10 0.00 80.0 NIST 2.70e-11 0.38 0.0 HG93
CH2+ + $h\nu \rightarrow$ CH+ + H 4.60E-11 0.00 1.8 VD06 1.70e-09 0.00 1.7  

Biggest rate coefficient increase

O + NS $\rightarrow$ SO + N
1.00E-11 0.00 0.0 SH04 1.73e-11 0.50 4000.0 MI84
O + C6H $\rightarrow$ C5H + CO 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
O + C7H $\rightarrow$ C6H + CO 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
O + C8H $\rightarrow$ C7H + CO 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
O + C9H $\rightarrow$ C8H + CO 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
O + C3H $\rightarrow$ CO + C2H 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
O + C4H $\rightarrow$ C3H + CO 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
O + C5H $\rightarrow$ C4H + CO 1.70E-11 0.00 0.0 SH04 1.00e-10 0.00 250.0 HL89
N + C4H $_2^+\rightarrow $ HC4N+ + H 9.50E-12 0.00 0.0 9902 1.00e-18 0.00 0.0  
CN + O $_2\rightarrow$ NO + CO 5.01E-12 -0.46 -8.0 NIST 5.30e-13 0.00 0.0  


The table of reactions is ordered in broad categories as listed in Table 7. Within each category, ordering is done by mass such that species R1 is the lightest of the collision partners, and, for a given R1, R2 is listed in order of increasing mass. Figure 1 shows the changes in the source of the data and accuracy of the rate coefficients between RATE99 and RATE06. As expected, the fraction of measured reactions and the accuracy of the data has improved with time. In addition to the one-third or so reactions which have measured rate coefficients, many of the photorates are, of course, based on experimental cross-sections. We describe each category in turn, noting important changes to the data where appropriate.


  \begin{figure}
\par\mbox{\includegraphics[width=7cm]{4981f1a.ps}\hspace{3mm}
\includegraphics[width=7cm]{4981f1b.ps} }
\par
\end{figure} Figure 1: The rate coefficient data in RATE06 (dark) is more accurate, and based on more laboratory measurements than RATE99 (light). The panels show a comparison of the numbers of reactions versus their data source ( left) and quoted accuracy ( right), between the two releases.


 

 
Table 7: The types of chemical reactions present and their count in the database, together with the number for which experimental data is available. The first column gives the abbreviation used for the types.
  Reaction type Count Measured

NN
Neutral-Neutral 549 275
IN Ion-Neutral 2387 956
CE Charge Exchange 552 201
II Atomic Ion-Ion Neutralisation 31 0
DR Dissociative Recombination 486 95
RR Radiative Recombination 25 0
AD Associative Detachment 47 1
RA Radiative Association 91 17
PH Photoprocess 216 0
CP Cosmic-Ray Proton (CRP) 11 0
CR Cosmic-Ray Photon (CRPHOT) 156 0
CD Collisional Dissociation 16 2
CI Chemi-ionisation 1 0
IM Ion-Molecular Ion Neutralisation 5 0
CL Collider 32 11


3.1 Neutral-neutral reactions

Since the time of our last release we have increased the number of neutral-neutral reactions by just over 100, removing around 10 reactions in the process, mostly because they are endothermic or because new laboratory measurements have indicated different products or negligible rate coefficients. Our new reactions are the result of an extensive search of the NIST Chemistry Webbook (http://webbook.nist.gov/chemistry) as well as more recent literature (the NIST Webbook has not included new data since around 2001) for reactions likely to be of interest in interstellar, circumstellar and shock chemistry. Many of the rate coefficients have been taken from the series of papers on evaluated neutral-neutral rate coefficients published by the IUPAC Subcommittee on Gas Kinetic Data Evaluation for Atmospheric Chemistry (much of its data is available on the web site http://www.iupac-kinetic.ch.cam.ac.uk), while we have made extensive use of the paper by Smith et al. (2004) which discusses the low-temperature behaviour of neutral-neutral reactions.

3.2 Cation-neutral reactions

Almost 2950 reactions in the database fall into the category of binary reactions between a cation and a neutral species, with around 550 of these involving charge exchange rather than bond breaking. We have reviewed these reactions thoroughly, mostly using the recent - and final - release of the Anicich (2004) database which includes "all'' ion-neutral reactions measured from 1936 to the end of 2003. Around 1160 of the ion-neutral reactions listed in our database have measured rate coefficients, mostly at room temperature. Again, where low-temperature information is available - and more is available for this category than for neutral-neutral reactions - we have included this in the database. For exothermic reactions involving species with small electric dipole moments, we expect that the rate coefficients will be temperature-independent. For reactions involving species with a large electric dipole moment, we expect a dependence on the order of T-1/2 and we adopt this in our "dipole-enhanced'' rate file, see Sect. 2.

3.3 Dissociative recombination reactions

This category of reaction now involves some 500 product channels of which close to 100 have been measured. This is an area which is significantly different from previous releases due to the large amount of high-quality data which have resulted from storage ring measurements, many of which involve recombination at low collisional energies. The ability of these experiments to determine product channels has given rise to some unexpected results which are challenging traditional models of gas-phase interstellar chemistry. They include the fact that the dissociative recombination of N2H+ breaks the N-N bond in two-thirds of the collisions (Geppert et al. 2004), and that of CH3OH2+ produces methanol, CH3OH, in only 3-4% of recombinations (Geppert et al. 2006).

3.4 Radiative association

The database includes 91 radiative associations of which 17 have rate coefficients determined with the aid of experimental data, often three-body data. New information on low-temperature radiative association processes has become available due to the development of ion-trap systems which allow reactants to thermalise at low temperature. Of particular importance here is the association of CH3+ and water to form protonated methanol, CH3OH2+, for which Herbst (1985) calculated a rate coefficient of $5.5\times10^{-12}$(T/300 K)-1.7 cm3 s-1, based on experimental three-body measurements. More recently, Luca et al. (2002) have reported a direct measurement of the rate coefficient below 80 K and find it to be $2.0\times10^{-12}$ cm3 s-1, independent of temperature. This low value, some 1000 times smaller than the Herbst value at 10 K, makes it difficult to form CH3OH at interstellar temperatures. In spite of the fact that this result has not, as yet, been confirmed by the publication of a refereed article, we adopt it in the new version of database, and the model results presented in Sect. 4 use it.

3.5 Photoprocesses

Some additional photoprocesses have been incorporated in to the database and a total of 216 product channels are now included. In part, these are due to the inclusion of photodissociation of the larger, organic molecules (Herbst & Leung 1986). Additional data for a few species have been provided by van Dishoeck (private communication).

3.6 Cosmic ray photon reactions

Doty noticed that the probabilities, Pi, adopted in RATE99 and in the OSU databases are too large by a factor of two, reproducing as they do the original Pi listed by Sternberg et al. (1987). Thus, the probabilities in the older implementations at UMIST and OSU should be reduced by a factor of two. This error occurred because in these databases the cosmic-ray ionisation rate is per H2 molecule - as opposed to per H-atom in Sternberg et al. (1987) - while the dust extinction, implicit in the determination of the probabilities, is proportional to the total density, i.e. 2n(H2), in dense clouds. The probabilities of the 156 cosmic-ray photo-induced reactions listed here have been corrected for this error.

3.7 Fluorine chemistry

Following the detection of HF and CF+ in interstellar clouds (Neufeld et al. 1997; Neufeld et al. 2006) we have included a small number of reactions describing the chemistry of fluorine. This is based on the work of Neufeld et al. (2005) but we have, in some cases, used alternative rate coefficients where these have been better determined experimentally. For example, we adopt the experimental rate coefficient of F+ + H2 from Anicich (2004), while for F + H2O, we adopt the critically-evaluated rate coefficient of $1.4\times 10^{-11}$ cm3 s-1 (Atkinson et al. 2005), ten times less than that adopted by Neufeld et al. (2005). Finally, we have used the measured rate coefficient for the dissociative recombination of the observed ion CF+ (Novotny et al. 2005). The measured rate coefficient of $5.2\times 10^{-8}$(T/300)-0.8 cm3 s-1 is smaller than that adopted by Neufeld et al. (2005), $2.0\times 10^{-7}$(T/300)-0.5 cm3 s-1, although at 10 K the measured rate coefficient is smaller only by about 30%.
  \begin{figure}
\par\includegraphics[angle=-90,width=5.8cm,clip]{4981f2a.ps}\hspace*{9mm}
\includegraphics[angle=-90,width=5.8cm,clip]{4981f2b.ps}\end{figure} Figure 2: The time evolution of the major sulphur-bearing species with RATE99 ( left) and with RATE06 ( right).

   
4 Results

Before discussing model results with RATE06, we compare in Fig. 2, the model results from RATE99 and RATE06 for a subset of species, namely the sulphur species. Here we have chosen parameters appropriate for a cold, dense cloud, namely T = 10 K, $n(\rm H_2) = 10^4$ cm-3, $A_{\rm V} = 15$ mag, and with initial elemental abundances, typical of the low-metal abundances often adopted for TMC-1, given in Table 8. This figure shows that while some changes in abundances can be seen, the overall tendency is for the abundances to evolve in a similar manner and to similar abundances. Indeed, most of the changes which do occur in these abundances are the result of changes to the oxygen chemistry, in particular the O and OH abundances. Since OH is a very reactive radical at low temperatures, it can play an important role in determining the abundances of simple oxides while O atoms can destroy hydrocarbons very efficiently. As noted by Smith et al. (2004) many of the reactions which cause the greatest sensitivity in model results involve atomic oxygen.


   
Table 8: Initial elemental abundances relative to total hydrogen nuclei.

Element
Abundance Element Abundance Element Abundance Element Abundance

He
0.14 C+ $7.30\times 10^{-5}$ N $2.14\times 10^{-5}$ O $1.76\times 10^{-4}$
F $2.00\times 10^{-8}$ Na+ $3.00\times 10^{-9}$ Mg+ $3.00\times 10^{-9}$ Si+ $3.00\times 10^{-9}$
P+ $3.00\times 10^{-9}$ S+ $2.00\times 10^{-8}$ Cl+ $3.00\times 10^{-9}$ Fe+ $3.00\times 10^{-9}$

Table 9 shows steady state output for both the standard database and the dipole-enhanced database, while Figs. 3 and 4 show the time-dependent evolution of the major oxygen and sulfur-bearing species both with and without the inclusion of dipole-enhanced rate coefficients. Because of the relatively low level of ionisation in dark clouds, radicals are most rapidly destroyed by fast neutral-neutral reactions if available, whereas stable neutrals are destroyed through reaction with cations such as H3+, He+ and HCO+. As can be seen from the table and figures, the inclusion of dipole-enhanced rate coefficients reduces the abundances of the relevant neutrals, e.g. H2O, HNC, HCN, SO2 and H2S, by factors of 3-5 whereas radicals such as OH and NO are relatively unaffected. CH3CN and the cyanopolyynes, HC2n+1N, are among the most affected species because of their limited chemistry; they decrease by an order of magnitude when dipole-enhanced rate coefficients are used.


  \begin{figure}
\par\includegraphics[angle=-90,width=5.8cm,clip]{4981f3a.ps}\hspace*{9mm}
\includegraphics[angle=-90,width=5.8cm,clip]{4981f3b.ps}\end{figure} Figure 3: The time evolution of the major oxygen-bearing species without ( left) and with ( right) dipole-enhanced rate coefficients in Rate06.


  \begin{figure}
\par\includegraphics[angle=-90,width=5.8cm,clip]{4981f4a.ps}\hspace*{9mm}
\includegraphics[angle=-90,width=5.8cm,clip]{4981f4b.ps}\end{figure} Figure 4: The time evolution of the major sulfur-bearing species without ( left) and with ( right) dipole-enhanced rate coefficients in Rate06.


  \begin{figure}
\par\includegraphics[angle=-90,width=5.8cm,clip]{4981f5a.ps}\hspace*{9mm}
\includegraphics[angle=-90,width=5.8cm,clip]{4981f5b.ps}\par\end{figure} Figure 5: The time evolution of the major sulphur-bearing species with ( left) the rate coefficient for the O + OH reaction taken from Smith et al. (2004) and with ( right) the value taken from Carty et al. (2006). Both model calculations use dipole-enhanced rate coefficients.

The low temperature dependence of the rate coefficient for the O + OH reaction plays a critical role in determining the abundances of many species in which OH plays a synthetic role. As an example of the effect of the sensitivity of calculated abundances to this rate coefficient, we plot in Fig. 5 the time-dependent abundances of several sulphur-bearing species for two cases. One in which we use a rate coefficient of $7.15\times 10^{-11}$(T/300 K)-0.25 cm3 s-1 (Smith et al. 2004) - $1.67\times 10^{-10}$ cm3 s-1 at 10 K - the other with a rate coefficient of $3.50\times 10^{-11}$ cm3 s-1 measured by Carty et al. (2006) over the range 39-142 K and adopted here for 10 K. The lower rate coefficient in the latter case leads to a larger fractional abundance of OH and hence to significantly larger fractional abundances of SO and SO2. In the non-dipole enhanced calculation, the use of the Carty et al. (2006) rate coefficient leads to abundances of SO and SO2 comparable to S (Fig. 4).


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[angle=-90,clip]{4981f6.ps}}\end{figure} Figure 6: The time evolution of the methanol abundance for various choices of key rate coefficients. The curves calculated using the latest RATE06 rate coefficients are labelled "RATE06'' and "RATE06-dipole'' for the non-dipole and dipole versions respectively. The curve labelled "RATE99'' uses the rate coefficients for the radiative association of CH3+ and H2O and the dissociative recombination of CH3OH2+ from that database.

Finally, we plot the methanol abundance as a function of time for several choices of rate coefficients in Fig 6. The curve labelled "RATE99'' uses the data from RATE99, in particular it assumes that 50% of the dissociative recombinations of protonated methanol lead to methanol, and that the rate coefficient for the radiative association of CH3+ and H2O is that determined by Herbst (1985); in contrast, the curve labelled "RATE06'' uses the latest experimental data for the branching ratios of the dissociative recombination of protonated methanol from Geppert et al. (2006), and the low experimental rate for the radiative association of CH3+ and H2O as measured by Luca et al. (2002). The "RATE99'' curve gives values of the methanol abundance close to those observed in dark clouds, but the abundance falls by over an order of magnitude when the new experimental dissociative recombination rate coefficients are introduced. The inclusion of dipole-enhanced rate coefficients (curve "RATE06-dipole'') also has an effect since the water, which has a large dipole moment, has a lower abundance in such a case. For the most extreme case, the methanol fractional abundance is never above 10-10 and at steady-state is about 1000 times less than that predicted by RATE99. This is the most significant astrochemical difference between RATE99 and RATE06.

Acknowledgements
Several individuals, including S. Doty, A. Sternberg and E. F. van Dishoeck, made useful comments on earlier versions of the database. We are grateful to PPARC for supporting the work of J.W. and T.J.M. through the UMIST Rolling Grant and to the EU through contract number MRTN-CT-2004-512302. M.A. acknowledges the Spanish MEC for support under project AYA 2003-02785-E and grant AP2003-4619.

References

 

  
5 Online Material


 

 
Table 1: List of species contained in the database ordered by mass.

e-

H H- H+ H2+ H2 H3+ He+
He HeH+ C+ C C- CH CH+ CH2
CH2+ N+ N NH+ NH CH3 CH3+ NH2+
O+ CH4+ NH2 O- CH4 O NH3 OH-
OH+ NH3+ CH5+ OH NH4+ H2O H2O+ F
H3O+ F+ HF HF+ H2F+ Na+ Na C2+
Mg+ Mg C2 C2H C2H+ C2H2 C2H2+ CN
CN- CN+ HNC HCN C2H3 C2H3+ HCN+ H2CN
CO+ Si+ N2+ Si N2 C2H4 H2NC+ HCNH+
C2H4+ CO SiH+ C2H5 HCO HN2+ SiH HCO+
HOC+ C2H5+ CH2NH H2CO+ SiH2 NO CH3CH3 SiH2+
CH4N+ NO+ CH2NH2+ CH3CH3+ H2CO HOCH SiH3 SiH3+
H3CO+ C2H7+ P P+ HNO HNO+ CF+ S+
SiH4 O2+ SiH4+ S CH3OH+ S- PH H2NO+
PH+ CH3OH O2 O2H PH2+ O2H+ SiH5+ CH3OH2+
HS+ PH2 HS H2S+ PH3+ H2O2 H2S H3S+
Cl Cl+ C3+ C3 HCl HCl+ C3H H2Cl+
C3H+ C3H2 C3H2+ C2N+ H2CCC CNC+ C2N H2C3H+
C2NH+ C3H3 C3H3+ C2O+ CH2CN CH2CN+ CH3CCH+ C2O
CH3CCH SiC+ SiC C3H5+ CH3CN CH3CN+ HC2O+ HCSi
HCSi+ CH2CO+ OCN+ SiCH2+ NH2CN C3H6+ CH2CO SiN+
SiCH2 OCN CH3CNH+ SiN C3H7+ SiCH3 CP+ HNSi+
HNSi SiCH3+ CH3CO+ CP HNCO+ NH2CNH+ N2O+ SiO+
CS+ CO2+ SiCH4+ N2O HCP+ CH3CHO+ SiNH2+ CS
HCP CH3CHO CO2 SiO HCS HCS+ SiOH+ HN2O+
PN PN+ CH3CHOH+ HCO2+ PCH2+ HPN+ HCOOH+ C2H5OH+
CH3OCH3 PCH3+ C2H5OH NS+ H2CS+ CH2PH NO2+ H2SiO
H2CS NO2 NS H2SiO+ HCOOH CH3OCH3+ HNS+ H3SiO+
CCl PO SiF+ PNH2+ PCH4+ CCl+ PO+ H3CS+
CH3OCH4+ C2H5OH2+ HCOOH2+ HPO+ C4+ C4 HPO PNH3+
SO+ SO C4H+ H2PO+ C4H H2CCl+ HSO+ C3N+
C4H2 C3N C4H2+ C4H3 ClO HC3N+ ClO+ HC3N
C4H3+ C2N2+ HC3NH+ NCCN SiC2 C4H4+ SiC2+ C3O+
C3O C4H5+ HC3O+ CH2CHCN SiC2H+ NCCNH+ H3C3N+ SiC2H
SiNC+ C3H2O+ H4C3N+ SiC2H2+ SiNC SiC2H2 CCP+ SiC2H3+
H3C3O+ C4H7+ SiNCH+ CCP C2S Fe+ Fe C2S+
HC2P HC2P+ HC2S+ PC2H2+ PC2H3+ CH3COCH3 CH3COCH3+ CH3COCH4+
CH3CS+ PC2H4+ COOCH4+ C5 OCS+ C5+ HCOOCH3 SiS+
OCS SiS SiO2 HOCS+ C5H HSiS+ C5H+ HSiO2+
H5C2O2+ C4N C5H2+ C4N+ C5H2 C5H3+ HC4N+ SO2
H2C4N+ C5H4+ S2 SiC3+ CH3C4H SO2+ SiC3 S2+
SiC3H C5H5+ HSO2+ CH3C3N HS2+ CH3C3N+ HS2 SiC3H+
H2S2+ CH3C3NH+ SiC3H2+ H2S2 C3P NCCNCH3+ H3S2+ PC3H+
C3S+ C3S HC3S+ C6+ C6 C6H C6H+ C5N+
C6H2 C5N C6H2+ HC5N HC5N+ C6H3+ C6H4+ H2C5N+
SiC4+ SiC4 SiC4H+ C6H5+ H3C5N+ C6H6+ C6H6 C6H7+
C4P C4P+ PC4H+ C4S+ C4S HC4S+ C7+ C7
C7H+ C7H C7H2 C7H2+ C7H3+ CH3C6H C7H4+ CH3C5N
C7H5+ CH3C5NH+ C8 C8+ C8H+ C8H C7N+ C7N
C8H2 C8H2+ HC7N+ C8H3+ HC7N C8H4+ H2C7N+ H3C7N+
C8H5+ C9+ C9 C9H+ C9H C9H2 C9H2+ C9H3+
C9H4+ C9H5+ CH3C7N CH3C7NH+ C10+ C10 C9N C9N+
HC9N+ HC9N H2C9N+ H3C9N+        



   
Table 2: Permanent electric dipole moments in Debye of the neutral molecules. References are given for those cases in which the values of the dipole moments have been updated since RATE99.

Species
$\mu_{\rm D}$ Species $\mu_{\rm D}$ Species $\mu_{\rm D}$ Species $\mu_{\rm D}$ Species $\mu_{\rm D}$

H2
0 CH 1.46 CH2 0.57 NH 1.3 CH3 0
NH2 1.83 CH4 0 NH3 1.47 OH 1.66 H2O 1.85
HF 1.83 C2 0 C2H 0.8 C2H2 0 CN 1.45
HNC 2.7 HCN 2.98 C2H3 1.5 H2CN 2.54 N2 0
C2H4 0 CO 0.11 C2H5 0 HCO 1.0 SiH 0.12
CH2NH 2.02 SiH2 0.18 NO 0.15 C2H6 0 H2CO 2.33
HOCH 0 SiH3 0 HNO 1.67 SiH4 0 PH 0.64
CH3OH 1.7 O2 0 O2H 2.09 PH2 0.0 HS 0.76
H2O2 1.57 H2S 0.97 C3 0 HCl 1.08 l-C3H 3.1
c-C3H2 3.4 H2CCC 4.1 C2N 0.6 C3H3 4.0 CH2CN 3.99
C2O 1.3 CH3CCH 0.78 SiC 1.7 CH3CN 3.92 HCSi 0.07
NH2CN 0 CH2CO 1.42 SiCH2 ... OCN 0.64 SiN 2.3
SiCH3 0.65 HNSi 0.16 CP 0.86 N2O 0.16 CS 1.96
HCP 0.3 CH3CHO 2.69 CO2 0 SiO 3.1 HCS 1.02
PN 2.75 CH3OCH3 1.3 C2H5OH 1.44 CH2PH 0 H2SiO 3.82
H2CS 1.65 NO2 0.32 NS 1.81 HCOOH 1.41 CCl 1.22
PO 1.88 C4 0 HPO 2.33 SO 1.55 C4H 0.9
C4H2 0 C3N 2.2 C4H3 0 ClO 1.24 HC3N 3.6
NCCN 0 SiC2 2.39 C3O 2.39 CH2CHCN 3.89 SiC2H 1.4
SiNC 2.03 SiC2H2 2.5 CCP 3.5 C2S 2.8 HC2P ...
C2H6CO 2.8 C5 0 HCOOCH3 1.77 OCS 0.71 SiS 1.73
SiO2 0.5 C5H 4.3 C4N 0.14 C5H2 5.9 SO2 1.63
S2 0 CH3C4H 1.21 SiC3 4.2 SiC3H 2.0 ? CH3C3N 4.91
HS2 1.4 H2S2 1.2 C3P 2.77 C3S 3.7 C6 0
C6H 5.0 C6H2 6.2 C5N 3.38 HC5N 4.33 SiC4 6.3
C6H6 0 C4P 4.18 C4S 3.8 C7 0 C7H 4.5
C7H2 7.7 CH3C6H 1.5 CH3C5N 5.75 C8 0 C8H 5.0
C7N 3.0 C8H2 8.2 HC7N 4.62 C9 0 C9H 4.7
C9H2 9.7 CH3C7N 5.47 C10 0 C9N 3.3 HC9N 4.84
Notes: c-C3H = 2.4D; HCCCCH = 0D; H2C3H(propargyl) = 0.14D; H2CCCC = 4.1D from Oswald & Botschwina (1995); HCSi from Smith et al. (2001); HCS from Ochsenfeld et al. (1999); HS2 from Zhuo et al. (1994); CCP from El-Yazal et al. (1997); CCl from Largo et al. (2001); H2SiO from Ma & Schaeffer (1994); CH2CN from Gutsev & Adamowicz (1995); SiCH3 from Kaiser & Osamura (2005); C3P from del Rio et al. (1996); C4P from Yu et al. (2005); C4S from Pascoli & Lavendry (1998); C5N from Botschwina (1996); CnH2, n = 5-9, from Maluendes & McLean (1992).


   
Table 4: List of the database reaction reference codes and their correspondance with the bibliographical references. NOTE: Any digit-only code found in the database is the same code as that used for referencing reactions in Anicich (2004).

Code
References Code References         Code Reference

AA83
Alge et al. (1983) AB93 Abouelaziz et al. (1993)         AM06 Andreazza et al. (2006)
AS95 Andreazza et al. (1995) AS84 Adams et al. (1984)         AS97 Andreazza & Singh (1997)
AS04 Asvany et al. (2004) BA84 Barlow (1984)         BH95 Bettens & Herbst (1995)
BS97 Brownsword et al. (1997) BV06 Barinovs & van Hemert (2006)         CA01 Canosa et al. (2001)
CG06 Carty et al. (2006) CH94 Clary et al. (1994)         DA97 Decker & Adams (1997)
DD90 Dalgarno et al. (1990) EH04 Ehlerding et al. (2004)         FA80 Field et al. (1980)
FM01 Fikri et al. (2001) FO02 Fockenberg & Preses (2002)         GA04 Geppert et al. (2004b)
GA06 Geppert et al. (2006) GB04 Geppert et al. (2004c)         GC04 Geppert et al. (2004d)
GE05 Montaigne et al. (2005) GH04 Geppert et al. (2004a)         GL87 Gredel et al. (1987)
GL89 Gredel et al. (1989) GT04 Geppert et al. (2004g)         GW04 Geppert et al. (2004f)
HD89 Herbst et al. (1989a) HE85 Herbst (1985)         HE87 Herbst (1987)
HE06 Heber et al. (2006) HH93 Haider & Husain (1993)         HL86 Herbst & Leung (1986a)
HL89 Herbst & Leung (1989) HL90 Herbst & Leung (1990)         HM89 Herbst et al. (1989b)
HT00 Herbst et al. (2000) IS89 Smith (1989)         JE00 Jensen et al. (2000)
KA02 Kalhori et al. (2002) KD93 Kimura & Dalgarno (1993)         KS96 Kaiser et al. (1996)
LG88 Leen & Graff (1988) LL98 Larson et al. (1998)         LU02 Luca et al. (2002)
MB73 McGregor & Berry 1973 MB87 Millar et al. (1987)         MC04 McCall et al. (2004)
MH90 Millar & Herbst (1990) MH91 Millar et al. (1991)         MI84 Mitchell (1984)
MI90 Mitchell (1990) MI91 Millar (1991)         MM93 Maluendes et al. (1993)
MS99 McEwan et al. (1999) NA99 Nahar (1999)         NIST NIST Chemical Kinetics Database
NM90 Nesbitt et al. (1990) NM05 Novotny et al. (2005)         NP97 Nahar & Pradhan (1997)
NW05 Neufeld et al. (2005) PD89 Petuchowski et al. (1989)         PE98 Peterson et al. (1998)
PH80 Prasad & Huntress (1980) PM03 Petrie et al. (2003)         RA92 Rawlings (1992)
RI99 Rim & Hershberger (1999) RJ91 Roberge et al. (1991)         RO00 Rosen et al. (2000)
SA00 Singh & Andreazza (2000) SD98 Stancil & Dalgarno (1998)         SE98 Semaniak et al. (1998)
SE01 Semaniak et al. (2001) SH04 Smith et al. (2004)         SI99 Singh et al. (1999)
SI06 Sims (2006) SL98 Stancil et al. (1998)         SS94 Smith et al. (1994)
SS99 Stancil et al. (1999) ST00 Stancil et al. (2000)         SY92 Suzuki et al. (1992)
TE96 Talbi et al. (1996) VA99 Vikor et al. (1999)         VD87 van Dishoeck (1987)
VD88 van Dishoeck (1988) VD06 van Dishoeck et al. (2006)         VP83 Viggiano & Paulson (1983)
WB88 Wlodek et al. (1988) WG04 Geppert et al. (2004e)         ZS98 Zygelman et al. (1998)

The NIST Chemical Database is available on-line at http://kinetics.nist.gov/index.php


 

 
Table 9: Steady state fractional abundances with respect to H2. "Rate06'' refers to the standard database while "Dipole'' refers to the database with dipole-enhanced rate coefficients.

Species
Rate06 Dipole Species Rate06 Dipole Species Rate06 Dipole

H
2.21E-04 2.22E-04 H- 9.37E-13 9.43E-13 H+ 1.26E-10 1.24E-10
H2+ 5.68E-13 5.68E-13 H2 1.00E+00 1.00E+00 H3+ 3.19E-09 3.15E-09
He+ 4.97E-10 4.96E-10 He 2.80E-01 2.80E-01 HeH+ 1.38E-14 1.38E-14
C+ 1.37E-09 1.34E-09 C 6.21E-10 6.15E-10 C- 1.71E-19 1.69E-19
CH 4.38E-11 4.43E-11 CH+ 1.75E-17 1.79E-17 CH2 1.38E-11 1.34E-11
CH2+ 7.08E-16 6.97E-16 N+ 3.37E-11 3.34E-11 N 1.89E-06 1.93E-06
NH+ 4.57E-15 4.55E-15 NH 5.58E-09 5.34E-09 CH3 5.06E-11 4.79E-11
CH3+ 2.28E-12 2.24E-12 NH2+ 2.25E-14 2.34E-14 O+ 6.08E-14 6.02E-14
CH4+ 6.66E-16 6.53E-16 NH2 7.83E-10 7.91E-10 O- 3.73E-18 3.85E-18
CH4 6.58E-08 6.53E-08 O 4.91E-05 5.07E-05 NH3 1.93E-07 4.16E-08
OH- 2.32E-14 4.94E-15 OH+ 2.35E-13 2.37E-13 NH3+ 5.27E-12 5.48E-12
CH5+ 6.33E-12 6.18E-12 OH 1.46E-07 1.40E-07 NH4+ 5.35E-11 6.00E-11
H2O 1.99E-06 3.95E-07 H2O+ 4.64E-13 4.79E-13 F 1.59E-08 3.47E-08
H3O+ 4.35E-09 4.52E-09 F+ 1.31E-16 1.59E-16 HF 2.41E-08 5.33E-09
HF+ 3.84E-17 4.64E-17 H2F+ 6.86E-12 8.23E-12 Na+ 5.72E-09 5.72E-09
Na 2.78E-10 2.79E-10 C2+ 1.10E-18 1.24E-18 Mg+ 5.47E-09 5.46E-09
Mg 5.29E-10 5.39E-10 C2 5.47E-12 6.31E-12 C2H 9.79E-11 9.29E-11
C2H+ 2.41E-18 2.53E-18 C2H2 1.38E-09 1.36E-09 C2H2+ 3.64E-15 3.56E-15
CN 1.04E-11 1.12E-11 CN- 7.35E-17 1.78E-17 CN+ 2.98E-17 5.76E-17
HNC 1.02E-08 2.40E-09 HCN 5.89E-09 1.43E-09 C2H3 8.68E-14 8.27E-14
C2H3+ 4.94E-13 4.81E-13 HCN+ 6.94E-17 9.98E-17 H2CN 4.68E-14 4.47E-14
CO+ 2.88E-14 2.89E-14 Si+ 8.54E-12 2.76E-11 N2+ 2.59E-15 2.61E-15
Si 5.04E-12 2.59E-11 N2 1.62E-05 1.64E-05 C2H4 1.43E-11 1.49E-11
H2NC+ 1.15E-11 1.33E-11 HCNH+ 1.28E-11 1.67E-11 C2H4+ 3.99E-14 3.90E-14
CO 1.46E-04 1.46E-04 SiH+ 3.57E-16 1.82E-15 C2H5 2.68E-16 1.05E-16
HCO 2.46E-12 2.32E-12 HN2+ 2.81E-10 2.79E-10 SiH 1.17E-15 3.97E-15
HCO+ 9.09E-09 8.88E-09 HOC+ 1.03E-13 1.04E-13 C2H5+ 3.06E-15 3.09E-15
CH2NH 9.86E-11 4.68E-12 H2CO+ 3.54E-14 4.90E-14 SiH2 1.21E-17 5.94E-17
NO 8.23E-06 8.07E-06 C2H6 4.15E-17 8.58E-18 SiH2+ 1.91E-16 6.13E-16
CH4N+ 8.58E-15 1.01E-14 NO+ 3.45E-10 3.36E-10 CH2NH2+ 2.19E-14 5.77E-15
C2H6+ 5.21E-20 1.98E-20 H2CO 6.56E-11 5.94E-11 HOCH 0.00E+00 0.00E+00
SiH3 2.26E-21 1.10E-20 SiH3+ 3.97E-18 2.01E-17 H3CO+ 5.34E-14 1.93E-13
C2H7+ 7.45E-22 1.50E-22 P 2.20E-09 4.53E-09 P+ 2.46E-13 4.97E-13
HNO 2.23E-09 2.18E-09 HNO+ 3.35E-10 3.22E-10 CF+ 8.12E-12 9.66E-12
S+ 2.61E-11 5.26E-11 SiH4 7.72E-19 3.99E-18 O2+ 3.82E-10 3.79E-10
SiH4+ 1.60E-28 8.10E-28 S 1.27E-08 2.62E-08 CH3OH+ 6.04E-17 2.35E-17
S- 2.43E-17 4.98E-17 PH 4.94E-12 9.80E-12 H2NO+ 2.37E-13 1.24E-12
PH+ 2.47E-13 5.00E-13 CH3OH 1.69E-12 4.05E-13 O2 7.32E-05 7.32E-05
O2H 6.20E-10 8.42E-10 PH2+ 3.38E-15 6.75E-15 O2H+ 1.39E-16 1.39E-16
SiH5+ 5.26E-22 2.64E-21 CH3OH2+ 3.19E-16 1.91E-16 HS+ 1.47E-12 2.96E-12
PH2 5.07E-14 9.98E-14 HS 1.63E-12 3.33E-12 H2S+ 1.39E-14 2.81E-14
PH3+ 1.13E-14 2.29E-14 H2O2 0.00E+00 0.00E+00 H2S 6.67E-11 3.11E-11
H3S+ 8.66E-14 1.88E-13 Cl 3.88E-09 5.40E-09 Cl+ 1.11E-17 1.59E-17
C3+ 8.72E-19 1.42E-18 C3 8.26E-12 9.75E-12 HCl 2.12E-09 6.03E-10
HCl+ 1.81E-17 2.54E-17 C3H 2.35E-12 1.78E-12 H2Cl+ 2.76E-13 4.06E-13
C3H+ 4.52E-17 4.65E-17 C3H2 1.71E-10 5.81E-11 C3H2+ 3.19E-15 5.57E-15
C2N+ 6.97E-13 8.84E-13 H2CCC 1.63E-12 1.30E-12 CNC+ 3.87E-13 5.04E-13
C2N 7.30E-12 6.94E-12 H2C3H+ 7.80E-14 8.51E-14 C2NH+ 8.78E-15 9.30E-15
C3H3 8.40E-14 6.82E-14 C3H3+ 5.12E-14 8.30E-14 C2O+ 1.95E-18 1.01E-17
CH2CN 9.18E-11 5.21E-12 CH2CN+ 5.34E-15 1.63E-15 C3H4+ 2.76E-18 9.73E-18
C2O 7.88E-14 7.59E-14 CH3CCH 1.01E-13 9.78E-14 SiC+ 5.15E-22 9.85E-21
SiC 7.39E-16 3.19E-15 C3H5+ 1.43E-16 1.39E-16 CH3CN 1.39E-11 8.27E-13
CH3CN+ 2.04E-18 6.37E-19 HC2O+ 2.05E-14 2.07E-14 HCSi 4.22E-16 3.88E-15
HCSi+ 8.95E-22 1.14E-20 CH2CO+ 4.46E-15 3.52E-15 OCN+ 5.91E-15 5.88E-15
SiCH2+ 1.22E-17 1.25E-16 NH2CN 2.02E-10 4.75E-11 C3H6+ 1.69E-26 3.45E-27
CH2CO 1.09E-10 2.19E-11 SiN+ 1.01E-15 1.08E-14 SiCH2 3.98E-15 1.66E-14
OCN 1.14E-08 1.22E-08 CH3CNH+ 5.10E-14 1.54E-14 SiN 3.88E-13 1.43E-12
C3H7+ 1.45E-25 2.95E-26 SiCH3 9.81E-16 3.43E-15 CP+ 2.30E-23 5.06E-23
HNSi+ 7.34E-21 6.78E-20 HNSi 2.56E-11 9.94E-11 SiCH3+ 1.63E-18 6.71E-18
CH3CO+ 5.63E-14 5.85E-14 CP 1.48E-14 2.89E-14 HNCO+ 6.22E-19 1.31E-18
NH2CNH+ 6.11E-14 1.41E-14 N2O+ 1.36E-13 1.26E-13 SiO+ 1.02E-17 8.35E-17
CS+ 1.02E-17 4.51E-17 CO2+ 1.32E-17 2.67E-17 SiCH4+ 1.25E-18 4.35E-18
N2O 2.44E-09 2.51E-09 HCP+ 1.51E-21 2.98E-21 CH3CHO+ 1.58E-15 1.17E-15
SiNH2+ 3.15E-14 1.22E-13 CS 1.69E-10 8.19E-11 HCP 1.53E-12 3.09E-12
CH3CHO 2.43E-11 4.77E-12 CO2 1.47E-07 2.98E-07 SiO 5.83E-09 5.71E-09
HCS 5.41E-15 1.07E-14 HCS+ 3.57E-14 1.26E-13 SiOH+ 2.23E-12 1.21E-11
HN2O+ 4.95E-14 5.02E-14 PN 3.69E-09 1.27E-09 PN+ 1.00E-14 1.85E-14
CH3CHOH+ 6.13E-14 6.35E-14 HCO2+ 4.18E-12 1.58E-11 PCH2+ 7.80E-16 1.57E-15
HPN+ 1.01E-12 1.90E-12 HCOOH+ 1.61E-25 7.47E-27 C2H5OH+ 2.47E-17 1.68E-17
CH3OCH3 2.43E-17 1.18E-18 PCH3+ 1.40E-16 2.82E-16 C2H5OH 2.51E-13 6.14E-14
NS+ 3.56E-14 8.21E-14 H2CS+ 1.65E-16 3.57E-16 CH2PH 1.44E-16 2.79E-16
NO2+ 1.24E-16 1.42E-16 H2SiO 2.11E-17 1.98E-17 H2CS 2.80E-11 1.15E-11
NO2 5.87E-09 5.83E-09 NS 1.27E-12 2.42E-12 H2SiO+ 6.41E-19 3.35E-18
HCOOH 3.26E-08 1.51E-09 CH3OCH3+ 1.78E-21 3.41E-22 HNS+ 3.96E-16 4.06E-15
H3SiO+ 2.54E-20 1.29E-19 CCl 1.71E-15 2.54E-15 PO 9.78E-11 1.80E-10
SiF+ 2.90E-14 1.13E-13 PNH2+ 3.34E-16 8.01E-16 PCH4+ 9.37E-18 1.88E-17
CCl+ 5.22E-14 7.98E-14 PO+ 8.84E-13 1.85E-12 H3CS+ 9.91E-15 2.08E-14
CH3OCH4+ 3.59E-20 8.90E-21 C2H5OH2+ 4.64E-16 5.38E-16 HCOOH2+ 1.93E-11 4.78E-12
HPO+ 1.24E-13 3.67E-13 C4+ 9.24E-19 1.84E-18 C4 1.57E-15 2.35E-15
HPO 1.33E-13 2.61E-13 PNH3+ 4.43E-16 1.06E-15 SO+ 3.53E-12 6.40E-12
SO 1.06E-08 9.78E-09 C4H+ 3.17E-18 6.07E-18 H2PO+ 1.72E-15 4.11E-15
C4H 7.51E-13 1.16E-12 H2CCl+ 1.38E-16 2.11E-16 HSO+ 3.36E-12 1.67E-11
C3N+ 1.09E-21 1.20E-21 C4H2 1.20E-11 1.53E-11 C3N 4.66E-14 3.78E-14
C4H2+ 6.34E-15 1.31E-14 C4H3 1.67E-13 2.50E-13 ClO 0.00E+00 0.00E+00
HC3N+ 7.47E-19 2.26E-18 ClO+ 0.00E+00 0.00E+00 HC3N 1.16E-12 1.41E-13
C4H3+ 3.58E-15 5.41E-15 C2N2+ 3.22E-23 8.19E-24 HC3NH+ 3.01E-16 2.44E-16
NCCN 6.13E-13 1.54E-13 SiC2 1.73E-15 6.72E-15 C4H4+ 5.72E-17 8.56E-17
SiC2+ 5.77E-22 1.84E-21 C3O+ 4.91E-20 1.44E-19 C3O 6.76E-11 3.68E-11
C4H5+ 1.95E-19 2.32E-19 HC3O+ 3.28E-14 9.70E-14 CH2CHCN 4.80E-17 9.17E-18
SiC2H+ 5.62E-17 1.95E-16 NCCNH+ 1.05E-16 2.60E-17 H3C3N+ 3.33E-21 3.29E-21
SiC2H 2.38E-14 1.54E-14 SiNC+ 2.39E-15 9.57E-15 C3H2O+ 8.64E-18 8.45E-18
H4C3N+ 3.42E-20 3.45E-20 SiC2H2+ 1.37E-17 4.84E-17 SiNC 4.34E-17 3.22E-17
SiC2H2 1.68E-14 1.11E-14 CCP+ 1.72E-17 3.43E-17 SiC2H3+ 1.56E-17 5.45E-17
H3C3O+ 1.17E-14 1.14E-14 C4H7+ 7.39E-22 7.55E-22 SiNCH+ 7.50E-19 7.37E-19
CCP 1.26E-15 2.21E-15 C2S 3.53E-12 1.82E-12 Fe+ 5.66E-09 5.65E-09
Fe 3.43E-10 3.48E-10 C2S+ 5.46E-21 1.29E-20 HC2P 7.97E-17 1.57E-16
HC2P+ 7.39E-18 1.74E-17 HC2S+ 1.73E-15 4.28E-15 PC2H2+ 7.47E-18 1.50E-17
PC2H3+ 1.97E-19 4.02E-19 C2H6CO 9.08E-17 3.40E-18 C2H6CO+ 1.12E-20 1.93E-21
C2H6COH+ 2.48E-19 4.88E-20 CH3CS+ 1.46E-18 3.07E-18 PC2H4+ 8.40E-22 1.76E-21
COOCH4+ 6.08E-25 1.60E-24 C5 4.60E-13 7.30E-13 OCS+ 6.13E-15 1.29E-14
C5+ 2.20E-21 1.45E-20 HCOOCH3 6.39E-21 3.84E-21 SiS+ 3.82E-18 2.57E-17
OCS 2.32E-10 4.84E-10 SiS 3.17E-17 3.91E-17 SiO2 1.28E-10 1.23E-10
HOCS+ 4.27E-14 8.76E-14 C5H 1.01E-14 8.88E-15 HSiS+ 7.12E-20 4.71E-19
C5H+ 2.12E-16 3.25E-16 HSiO2+ 5.19E-14 4.88E-14 H5C2O2+ 1.24E-23 3.83E-23
C4N 5.07E-14 7.69E-14 C5H2+ 7.73E-18 2.95E-17 C4N+ 1.01E-19 6.58E-20
C5H2 1.70E-13 1.82E-14 C5H3+ 2.94E-17 1.79E-17 HC4N+ 3.06E-22 4.81E-22
SO2 1.62E-08 3.37E-09 H2C4N+ 3.06E-22 4.81E-22 C5H4+ 7.44E-19 1.45E-18
S2 1.86E-15 7.66E-15 SiC3+ 8.86E-19 6.91E-18 C5H4 4.81E-15 9.04E-15
SO2+ 3.03E-16 3.41E-16 SiC3 3.86E-16 2.21E-15 S2+ 1.54E-16 6.71E-16
SiC3H 2.59E-17 6.41E-17 C5H5+ 2.79E-18 5.16E-18 HSO2+ 9.59E-13 1.08E-12
CH3C3N 1.40E-15 2.06E-15 HS2+ 1.24E-18 6.16E-18 CH3C3N+ 3.06E-22 4.81E-22
HS2 5.01E-14 4.27E-14 SiC3H+ 2.52E-17 1.55E-16 H2S2+ 4.55E-17 2.09E-16
CH3C3NH+ 2.60E-18 3.77E-18 SiC3H2+ 2.44E-20 3.23E-19 H2S2 9.78E-15 8.92E-15
C3P 3.76E-18 1.35E-17 NCCNCH3+ 9.21E-23 2.28E-23 H3S2+ 1.14E-17 5.70E-17
PC3H+ 1.17E-19 4.35E-19 C3S+ 4.84E-18 2.26E-17 C3S 2.87E-13 1.98E-13
HC3S+ 1.80E-16 6.51E-16 C6+ 1.78E-20 2.53E-20 C6 1.29E-16 1.71E-16
C6H 1.23E-15 1.51E-15 C6H+ 8.43E-19 1.04E-18 C5N+ 5.34E-23 8.36E-23
C6H2 3.56E-14 1.09E-14 C5N 2.11E-16 2.04E-16 C6H2+ 5.33E-18 7.41E-18
HC5N 1.78E-15 4.51E-16 HC5N+ 9.50E-23 1.96E-22 C6H3+ 3.35E-18 5.59E-18
C6H4+ 1.62E-19 3.69E-19 H2C5N+ 4.53E-18 7.20E-18 SiC4+ 1.07E-19 2.90E-18
SiC4 2.01E-21 1.74E-20 SiC4H+ 1.31E-22 1.24E-21 C6H5+ 8.91E-24 3.19E-23
H3C5N+ 3.47E-21 6.25E-21 C6H6+ 1.77E-21 6.23E-21 C6H6 1.30E-17 4.70E-17
C6H7+ 4.76E-21 1.71E-20 C4P 8.70E-18 3.21E-17 C4P+ 1.80E-23 3.58E-22
PC4H+ 2.91E-19 1.29E-18 C4S+ 4.04E-19 5.76E-18 C4S 6.15E-15 3.96E-15
HC4S+ 6.23E-18 2.02E-17 C7+ 3.22E-22 5.37E-22 C7 4.17E-15 7.11E-15
C7H+ 9.26E-19 1.54E-18 C7H 3.06E-17 2.27E-17 C7H2 2.74E-16 6.87E-17
C7H2+ 6.85E-21 2.91E-20 C7H3+ 5.38E-20 1.05E-19 C7H4 3.03E-17 4.89E-17
C7H4+ 5.22E-21 1.36E-20 CH3C5N 2.63E-19 1.26E-20 C7H5+ 2.57E-20 4.07E-20
CH3C5NH+ 2.82E-22 7.11E-23 C8 6.37E-19 1.08E-18 C8+ 1.29E-21 2.18E-21
C8H+ 6.55E-24 1.09E-23 C8H 3.98E-18 5.11E-18 C7N+ 2.33E-25 3.95E-25
C7N 6.77E-19 5.97E-19 C8H2 5.29E-17 5.32E-18 C8H2+ 3.12E-20 5.47E-20
HC7N+ 7.49E-23 1.43E-22 C8H3+ 5.21E-21 3.88E-21 HC7N 7.89E-18 1.76E-18
C8H4+ 1.72E-22 2.27E-22 H2C7N+ 4.42E-21 6.68E-21 H3C7N+ 3.82E-22 1.40E-22
C8H5+ 3.49E-22 6.72E-22 C9+ 2.40E-24 2.53E-24 C9 6.47E-18 5.63E-18
C9H+ 1.43E-21 1.19E-21 C9H 4.96E-20 2.37E-20 C9H2 5.06E-19 1.04E-19
C9H2+ 1.44E-23 4.67E-23 C9H3+ 1.19E-22 1.70E-22 C9H4+ 1.26E-23 2.02E-23
C9H5+ 2.87E-24 2.23E-24 CH3C7N 1.14E-21 4.79E-23 CH3C7NH+ 1.26E-24 2.79E-25
C10+ 5.25E-23 4.91E-23 C10 1.25E-24 5.74E-25 C9N 1.06E-21 5.77E-22
C9N+ 3.61E-28 3.06E-28 HC9N+ 1.16E-25 1.17E-25 HC9N 1.01E-20 1.45E-21
H2C9N+ 6.40E-24 5.78E-24 H3C9N+ 5.54E-25 3.03E-25 e- 3.70E-08 3.70E-08




Copyright ESO 2007