Free Access

This article has an erratum: [https://doi.org/10.1051/0004-6361/201118013e]


Issue
A&A
Volume 538, February 2012
Article Number A91
Number of page(s) 15
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201118013
Published online 10 February 2012

© ESO, 2012

1. Introduction

Oxygen is the most abundant “metal” element in the cosmos (Savage & Sembach 1996; Asplund et al. 2009). In the cold dense interstellar clouds, gas-phase chemical models predict that oxygen mainly resides in CO and O2 molecules (Herbst & Leung 1989; Millar & Herbst 1990; Wakelam et al. 2006). However, although CO is ubiquitously distributed in the interstellar medium, O2 is not. The latter is only detected very recently in ρ Oph A at a low abundance (relative to molecular hydrogen) of 5 × 10-8 (Larsson et al. 2007), and in Orion at an abundance of (0.3–7)  ×  10-6 (Goldsmith et al. 2011). On the other hand, the observed water (gas or ice) abundance can be as high as 10-4 (van Dishoeck 2004). Thus it seems that water, instead of O2, is a main reservoir of oxygen in addition to CO. When only gas phase chemistry is included, the H2O abundance can be on the order 10-7 at most (see, for example, Bergin et al. 2000; Roberts & Herbst 2002) for typical dark cloud conditions. That O2 is overproduced and H2O is underproduced in gas phase chemistry suggests that adsorption onto the grain surfaces and the reactions on the surfaces may play important roles.

On the grain surface, H2O can form through successive additions of hydrogen atoms to an oxygen atom: \begin{eqnarray} & &{\rm H} + {\rm O} \rightarrow {\rm OH}, \label{eq:HO} \\ & &{\rm H} + {\rm OH} \rightarrow {\rm H}_2{\rm O}, \label{eq:HOH} \end{eqnarray}both of which are barrierless (Allen & Robinson 1977). It can also form via hydrogen addition to molecular oxygen: H+O2O2H,H+O2HH2O2,H+H2O2H2O+OH.\begin{eqnarray} \label{eq:HO2} &&{\rm H} + {\rm O}_2 \rightarrow {\rm O}_2{\rm H} , \\ \label{eq:HHO2} &&{\rm H} + {\rm O}_2{\rm H} \rightarrow {\rm H}_2{\rm O}_2 , \\ \label{eq:HH2O2} &&{\rm H} + {\rm H}_2{\rm O}_2 \rightarrow {\rm H}_2{\rm O} + {\rm OH} . \end{eqnarray}Reaction (3) was assumed to have an activation barrier of 1200 K in Tielens & Hagen (1982). However, based on experimental results, Cuppen et al. (2010) have recently concluded that it is barrierless. Other possible formation pathways of water include the reaction between H2 and OH and the route with O3 as an intermediate (Tielens & Hagen 1982).

In the second route described above (Eqs. (3)–(5)), hydrogen peroxide (HOOH, also written as H2O2, which is adopted in this paper) appears as an intermediate product. Thus if this route is indeed important, a significant amount of H2O2 might form on the grain, and its gas phase counterpart could also be detectable if effective desorption mechanisms exist.

In the current mainstream gas phase reaction networks for astrochemistry, H2O2 is not efficiently formed in the gas phase. For example, in the 2009 version of the OSU network1, the only two reactions leading to the formation of H2O2 are H2+O2HH2O2+H,OH+OHH2O2+.\begin{eqnarray*} & &{\rm H}_2 + {\rm O}_2{\rm H} \rightarrow {\rm H}_2{\rm O}_2 + {\rm H}, \\ &&{\rm OH} + {\rm OH} \rightarrow {\rm H}_2{\rm O}_2 + h\nu. \end{eqnarray*}The first one has a large activation barrier of 104 K, rendering it inactive at low temperatures. H2O2 is mainly consumed by H2O2+OH+OH,OH+H2O2H2O+O2H,\begin{eqnarray*} & &{\rm H}_2{\rm O}_2 + h\nu \rightarrow {\rm OH} + {\rm OH}, \\ & &{\rm OH} + {\rm H}_2{\rm O}_2 \rightarrow {\rm H}_2{\rm O} + {\rm O}_2{\rm H}, \end{eqnarray*}the first of which is dissociation by cosmic-ray induced radiation. Other destruction channels by reacting with H and O are ineffective due to large activation barriers. At a temperature of 10 K and an H2 density of 104 cm-3, the steady-state abundance of H2O2 can be approximated by X(H2O2)103X2(OH)5×10-12.\begin{equation*} X({\rm H}_2{\rm O}_2) \simeq 10^3 X^2({\rm OH}) \simeq 5\times10^{-12}. \end{equation*}At a higher density, the abundance of H2O2 will be even less because OH is less abundant in this case. With the UMIST RATE06 network the abundance of H2O2 is essentially zero (Woodall et al. 2007). Thus if a substantial amount of H2O2 can be detected in the interstellar medium, then it must have been synthesized on the dust grains, rather than in the gas phase. This would also provide information and constraints on the formation route of H2O.

It has indeed recently been detected (for the first time) in the ρ Oph A cloud by Bergman et al. (2011b), at an abundance of  ~10-10, which is well above what would be predicted by the gas phase chemistry, indicating that chemical processes on the grains are responsible for this detection. Why this molecule has not been detected in the past seems to be a puzzle and will be discussed later.

In the present work we aim at modeling the gas phase abundance of H2O2 at a physical condition relevant to ρ Oph A, to demonstrate whether the grain chemistry is able to explain its observed abundance. The model is also required to give consistent abundances for other species detected earlier in this region. Also, ice and gas-phase abundance predictions for previously undetected species are made.

The remaining part of this paper is organized as follows. In Sect. 2 we describe the chemical model used in this work. In Sect. 3 we present the results of our modeling. The conclusions are in Sect. 4. Appendix A contains an explanation to a spike-like feature in the evolution curves of some species. The surface reaction network we use is listed in Appendix B, and the enthalpies of the surface species, which are needed in the chemical desorption mechanism (see Sect. 2) are listed in Appendix C.

2. Chemical model

For the gas phase chemistry, we use a subset of the UMIST RATE06 network2 (Woodall et al. 2007). Species containing Fe, Na, Mg, and Cl are excluded. In total 284 gas phase species and 3075 gas phase reactions are included. The cosmic-ray ionization rate is taken to be the canonical value of 1.36 × 10-17 s-1 (Woodall et al. 2007).

The surface chemical network is a combination of a selection of the reactions in Allen & Robinson (1977), Tielens & Hagen (1982), and Hasegawa et al. (1992), with the rates of a few reactions updated according to the recent experimental results and/or theoretical calculations. In total 56 surface species and 151 surface reactions are included (see Appendix B).

The binding energies of the surface species are either taken from Hasegawa & Herbst (1993) or estimated from the value of a similar species as in Garrod et al. (2008). These values are applicable to bare grains, i.e., grains without an ice mantle. A grain will be covered by ice (typically water) as adsorption and reaction proceeds, thus these values are not always appropriate. Ideally, they should be varied according to the real-time composition of the grain. The water ice mantle mainly affects the binding energies of species with a hydrogen bond, such as OH and H2O. The effect of this should be minor for our purpose, because most of the reactions involving a species with a hydrogen bond are primarily mediated by another reaction partner that does not have a hydrogen bond.

The barriers against surface diffusion are taken to be a fixed fraction of the binding energies. A range of values have been used for this fraction in the past, from 0.3 (Hasegawa et al. 1992) through 0.5 (Garrod et al. 2008) to 0.77 (Ruffle & Herbst 2000). We use a value of 0.77, based on the analysis of Katz et al. (1999). Because our model is mainly for a relatively high temperature (~20 K) in comparison with most of the previous models (predominantly for a temperature of  ~10 K), a low diffusion barrier for the surface species would lead to an unrealistic ice mantle composition. The effect of changing this parameter is discussed later. We allow H and H2 on a dust grain to migrate through quantum tunneling or thermal hopping, depending on which is faster. All the heavier species are only allowed to move by thermal hopping. The quantum tunneling and thermal hopping rates are calculated using the formulation of Hasegawa et al. (1992). For calculating the quantum tunneling rates, we use a barrier width of 1 Å. The exact value of this width depends on the composition and structure of the surface, which has not been fully quantified.

About the activation barrier of reaction (3), as there is a big discrepancy between the value adopted in the past and the value proposed recently based on experiments (Cuppen et al. 2010), we adopt an intermediate value of 600 K. However, the effect of varying this parameter is also tested during the modeling, and will be discussed later. Reaction (5) has a barrier of 1400 K in Tielens & Hagen (1982), while in Cuppen et al. (2010) this reaction is barrierless. We choose to adopt the latter result in this case, because too high a barrier for it would result in building up too much H2O2 on the grain surface.

Surface reactions with an activation barrier are allowed to proceed thermally or through quantum tunneling, depending on whichever is faster. The formula used to calculate the rates is also the same as in Hasegawa et al. (1992), and the reaction barriers are assumed to have a width of 1 Å, although different values are possible (Garrod & Pauly 2011).

The reaction rate of a two-body surface reaction

A+BC+···$$ {\rm A+B\rightarrow C+\cdots} $$is

[kdiff(A)+kdiff(B)]N(A)N(B)/NS,$$ [k_{\rm diff}({\rm A}) + k_{\rm diff}({\rm B})]N({\rm A})N({\rm B)}/N_{\rm S}, $$if A ≠ B. Here NS is the number of reaction sites of a grain, kdiff(A) and kdiff(B) are the diffusion rates of A and B, and N(A) and N(B) are the number of species A and B on a single grain. If A = B, then the reaction rate should be

kdiff(A)N(A)(N(A)1)/NS.$$ k_{\rm diff}(A)N({\rm A})(N({\rm A})-1)/N_{\rm S}. $$With a number density of reaction sites being 1015 cm-2, a dust grain with radius 0.1 μm has an NS of about 106.

As we are mainly concerned with the gas phase abundances of several species, their desorption mechanism must be treated carefully, especially if they are mainly produced on the grains. Besides the normal thermal desorption, species can also get evaporated episodically when a cosmic-ray hits a grain. This is treated in the same manner as in Hasegawa & Herbst (1993). Furthermore, the nonthermal desorption mechanism via exothermic surface reactions (for brevity we call it “chemical desorption”) proposed by Garrod et al. (2006; see also Watson & Salpeter 1972; Garrod et al. 2007; Cazaux et al. 2010) is also included. Here the products of the exothermic reactions on the grain have a probability of being directly ejected into the gas phase. The rate of such a desorption mechanism depends on the exoergicity of the reaction, as well as on the desorption energy of the products. A parameter characterizing the efficiency of this mechanism (the “a” parameter in Garrod et al. 2007) is introduced, which we take to be 0.1. The yield of chemical desorption is directly proportional to this “a” parameter, although it is not well constrained. The value we adopt here gives a good match to the observational results. See Sect. 3.6 for further discussion. The exoergicities of these reactions are estimated from the enthalpies of the reactants and products in the same manner as in Allen & Robinson (1977, their Eqs. (3) and (4)), and the enthalpies of the species involved in these reactions are taken from Binnewies & Milke (1999), from the NIST chemistry web book3, or from some other sporadic sources (see Appendix C).

However, the desorption mechanisms described above are not always sufficient to provide enough gas phase abundances for some species, especially at late times. Even if a large amount of a species is produced on the grain and released into the gas phase at early times, it would later be accreted back to the grain surface. If at this later time its production is no longer active (due to the exhaustion of the precursor species), its gas phase abundance cannot be maintained. Dust sputtering (Tielens et al. 1994) and photodesorption (Öberg et al. 2007) might help release them to the gas phase, both of which should not be very important in a quiet cold dark cloud. Another possible mechanism is that cosmic-ray induced radiation can dissociate the species on the grain, and when the fragments recombine, the products can possibly be ejected into the gas phase directly because of the energy release of the reaction, as described before. We implement this mechanism in the same way as in Ruffle & Herbst (2001a) and Garrod & Pauly (2011; see also Cuppen & Herbst 2007), namely, the cosmic-ray induced photodissociation rates for the surface species are taken to be the same as in the gas phase. Several dissociation branches from Garrod et al. (2008, their Table 1) are included.

In our model the numbers of all the species on a single grain are solved with the hybrid moment equation (HME) approach (Du & Parise 2011). It has been shown in Du & Parise (2011) that the rate equation method can be inaccurate in some cases, and the HME approach provides a major improvement over the rate equation method. Since the HME approach is relatively new, in several cases we also benchmarked our HME results with the exact Monte Carlo method (similar to the one of Vasyunin et al. 2009) as in Du & Parise (2011), and the agreement is satisfactory. It is impractical to run all the models with the Monte Carlo method because the run time would be too long.

At present the layered structure of the grain mantle is not incorporated into our model. Such a structure might be more realistic, and it is important for retaining some of the ice species. However, it is also possible that the interstellar dust grains may have an amorphous structure, which renders the layered structure an inaccurate description. On the other hand, particles landing on a grain are able to penetrate the interior by several to tens of layers, as demonstrated by experiments (see, e.g., Ioppolo et al. 2010), thus although a model neglecting the layered structure is not accurate, one that deactivates all the layers below the outermost surface does not reflect the whole reality either. In fact, for the H-addition reactions, whether the layered structure is taken into account only plays a minor role in determining the reaction rates in the accretion limit (i.e. when the accretion and evaporation processes are much slower than the reactions; Garrod 2008), because the species involved in these reactions never build up a full layer.

3. Results and discussions

3.1. Modeling ρ Oph A

We ran the chemical model for the physical parameters that are appropriate for ρ Oph A, where H2O2 has first been detected by Bergman et al. (2011b). In Fig. 1 we show the abundances of several species as a function of time. In this model, a temperature of 21 K and a hydrogen density of 6 × 105 cm-3 have been assumed, which are determined for ρ Oph A observationally by Bergman et al. (2011a,b). The dust temperature and gas temperature are assumed to be the same. A fixed value of 15 for the visual extinction AV has been adopted. We assume a canonical grain size of 0.1 μm and a reaction site density of 1015 cm-2. The dust-to-gas mass ratio is set to 0.01, and the dust grain material is assumed to have a mass density of 2 g cm-3. We also assume that the ratio between the diffusion barrier and the binding energy is 0.77, which is at the higher end of the values used in the past, to give a reasonable ice composition. The initial condition is atomic except for H2. The elemental abundances are the same as in Garrod & Pauly (2011).

The observed abundance of gas phase H2O2 is  ~10-10 relative to H2 (Bergman et al. 2011b). This value is best matched at a time of  ~6 × 105 yr. In the early time, before about 2 × 105 years, the gas phase H2O2 abundance can be as high as  ~5 × 10-7. At late times, the H2O2 abundance decreases to a very low value, due to the exhaustion of O2 on the grain and a full conversion of H2O2 into H2O.

thumbnail Fig. 1

The time evolution of the abundances of selected species. A prefix “g” means a grain surface species (right panel), while a name without a prefix means a gas phase species (left panel). The dotted lines are the observed abundances (relative to H2) of O2 (blue), H2CO (cyan), CH3OH (magenta), and H2O2 (black), respectively, in the ρ Oph A source.

H2O2 is mainly formed through reaction (4) on the grain, followed by immediate desorption of the product into the gas phase caused by the reaction heat. About 7% of the produced H2O2 is released this way. Its gas phase abundance is determined by the adsorption and chemical desorption processes. The dissociation of gas phase H2O2 by cosmic-ray-induced photons is unimportant in consuming it, compared to adsorption.

Ioppolo et al. (2008) modeled the abundance of H2O2 ice briefly, giving a value of 10-14–10-10 relative to molecular hydrogen, depending on which energy barriers of several relevant reactions were used. Our major goal is to model the gas phase H2O2 abundance, rather than the H2O2 ice. In our model results, the gas phase abundance of H2O2 is much higher than its ice counterpart. The H2O2 ice does not have a significant abundance in the later stage, being well below the upper limit (5.2% with respect to H2O ice) given by Boudin et al. (1998), because it is constantly transformed into H2O by reacting with the accreted H atoms. However, by irradiating thin water ice film with low energy ions, Gomis et al. (2004) found that it is possible to obtain an H2O2 to H2O ratio in the solid phase up to a few percent. This direct processing of the grain mantle by cosmic rays is not included in our model. The production of H2O2 inside water ice in an O2-rich environment triggered by UV radiation (Shi et al. 2011) should also be of little importance here. We notice in Fig. 1 (right panel) that in the early stage (before  ~5 × 104 yr) the H2O2 ice can achieve a rather high abundance,  ~5 × 10-6 relative to H nucleus or  ~10% relative to H2O. During this early period the water formation on the grain mainly proceeds through reaction (5) with H2O2 as an intermediate, which is responsible for about half of the final water ice repository on the grain mantle. In the later stage, reaction (2) takes over. In the results of Ioppolo et al. (2008, their Fig. 4), we do not see a similar feature (i.e. a high abundance of H2O2 ice in the early stage). In our current model the layered structure of the grain mantle is not taken into account. It is possible that, if such a structure is considered, the inner layers with a relatively high H2O2 content might be maintained, which would give a value of a few percent for the H2O2-to-H2O ratio in the solid phase.

Methanol (CH3OH) and formaldehyde (H2CO) are also detected in the ρ Oph A SM1 core (Bergman et al. 2011a), at an abundance of  ~2 × 10-9 and  ~5 × 10-9, respectively. Their abundances are also reproduced very well at a time of  ~6 × 105 yr in our model. At early times, both CH3OH and H2CO have a high abundance. Their abundances also have a peak in the period between 2 × 105 yr to 107 yr. In our current network, CH3OH is mainly formed on the grains, and mainly through the addition of H atom to CH2OH, while the latter is mainly produced from the reaction between C and OH to form HOC followed by successive H additions. Thus the abundance of CH3OH decreases at very late times due to the depletion of atomic C (which is mainly in CH4 ice in the late stage). The normal formation channel through successive hydrogenation of CO is important at around 5 × 105–2 × 106 yr. The gas phase H2CO mainly forms in the gas phase in the early stage (<105 yr), and mainly through the reaction CH3 + O → H2CO + H. Later it is mainly formed through successive hydrogenation of CO on the grain surface followed by chemical desorption. The abundance of methanol and formaldehyde ice relative to water ice can be as high as  ~20% at their peaks at a time of  ~2 × 106 yr, but falls to a very low value in the late times. The late-time abundances are consistent with the upper limit derived for quiescent environment and low-mass young stellar objects in Gibb et al. (2004). However, in Pontoppidan et al. (2004) a much higher abundance of CH3OH ice is observed along the line-of-sight of SVS4 (a dense cluster of pre-main sequence stars), which is close to a class 0 protostar, and this is consistent with the peak abundances in our model.

From Fig. 1 (left panel) it can be seen that the abundance of gaseous O2 at an intermediate time of 6 × 105 year is  ~6  ×  10-9, which is within one order of magnitude of the observed abundances of 5  ×  10-8 for O2 (Larsson et al. 2007). The late-time abundance of O2 drops to a very low value, while its observed abundance is best matched at a time of  ~2  ×  105 yr. We notice that during the period  ~(0.6–2)  ×  105 yr, the abundance of O2 ice has a prominent bump, reaching a peak abundance of  ~10-5 relative to H2. At the same time, the gas phase O2 also reaches an abundance of  ~ (1–5)  ×  10-7. These values can be compared with the recent detection of gas phase O2 at an abundance of (0.3–7)  ×  10-6 in Orion by Herschel (Goldsmith et al. 2011). Warm-up of the dust grain at this stage may release a large amount of O2 molecule into the gas phase.

As a precursor of H2O2, O2H mainly forms from the reaction between O and OH on the grain, which does not have a barrier according to Hasegawa et al. (1992). The ratio between the gas phase O2H and H2O2 is almost constant throughout the evolution, being approximately 3. Thus in our current network, the gas phase O2H also has a remarkable abundance, which might be detectable in the future.

Except at very early times, the grain mantle is mainly composed of water ice. The abundances of CO and CO2 are comparable at an intermediate time of (0.3–1)  ×  106 yr, being about 40–60% of water ice. This is in rough agreement with the ice composition for intermediate-mass YSOs (Gibb et al. 2004) (see also Öberg et al. 2011), and is also in line with the suggestion of An et al. (2011) that CO2 ice is mixed with CH3OH ice (the latter is about 10% of the former in our model). Water ice mainly forms from reaction (5) in the early time, and from reaction (2) in the late time. The gas phase formation route of water only plays a minor role. In contrast, the CO ice mantle is mainly from accretion of CO molecules formed in the gas phase. For the CO2 ice, it is mainly accreted from its gas phase counterpart in the early time, and its abundance is increased further through the reaction OH + CO → CO2 + H in the late stage. At late times (>3  ×  106 yr), most of the carbon resides in the form of CH4 ice, with the latter about half of the water ice. However, as far as we know, such a high abundance of CH4 ice (see also Garrod & Pauly 2011) has not been observed in the interstellar medium.

The gas phase CO is heavily depleted, with an abundance  ~10-6 relative to the H nucleus or  ~1% relative to its ice counterpart, at an intermediate time (105–106 yr). Its abundance is mainly determined by the balance between the adsorption and cosmic-ray induced evaporation processes. At late times the CO ice abundance drops to a very low value. This is because CO is continuously hydrogenated into H2CO or CH3OH, and the dissociation of CH3OH by cosmic-ray induced photons produces CH3, which quickly becomes CH4 by hydrogenation. Taking the layered structure of grain mantle into account can retain a fair amount of CO ice.

The abundance of gas phase H2O in the intermediate to late time is on the order 10-8. At these times, its grain surface production route (2) followed by chemical desorption and the gas phase production route H3O +  + e −  → H2O + H plays a similarly important role. H3O+ itself is mainly formed from successive protonation of atomic oxygen at this stage.

The hydroxyl radical (OH) has an abundance  ~10-8–10-9 in the gas phase at the intermediate-to-late times. These values are comparable to the observed abundance of  ~(0.5–1)  ×  10-8 in the envelope around the high-mass star-forming region W3 IRS 5 obtained with HIFI onboard Herschel (Wampfler et al. 2011). The physical condition in this source is different from ρ Oph A. However, since OH is readily produced and recycled in the gas phase by the reactions H3O +  + e −  → OH + 2H and H3++OHH2O++H2\hbox{${\rm H_3^++OH\rightarrow H_2O^++H_2}$}, grain processes will not play a dominant role in determining its abundance, especially when the temperature is not too low. However, Goicoechea et al. (2006) detected a much higher abundance (0.5–1)  ×  10-6 of OH in the Orion KL outflows, which seems to require other formation pathways (e.g. shock destruction of H2O ice).

A relatively high abundance (~10-7 and  ~10-5) of gas phase and grain surface ozone (O3) is also obtained for t ≲ 105 yr. However, these values should be treated with caution, because the chemistry of O3 is very incomplete in our current model: neither the OSU gas phase network nor the UMIST RATE06 network includes it. Possible gas phase destruction pathways involving atomic O and S (according to the NIST chemistry web book; see footnote 3) may lower its gas phase abundance significantly, but its ice mantle abundance should not be severely affected.

The abundance of atomic hydrogen in the gas phase is quite high at late times as seen from Fig. 1, which seems to be at odds with the usual results. In many gas phase chemical models, the abundance of atomic hydrogen in the gas phase is determined by the balance between its adsorption onto the dust grains and the dissociation of H2 molecules by cosmic rays. The adsorption process is thought to have a sticking coefficient close to unity, and evaporation is assumed not to occur (which is appropriate at low temperatures). In such a framework it can be found that the gas phase atomic hydrogen will always have a fixed density on the order of 1 cm-3. However, in our case with a temperature of 21 K, evaporation is very fast and cannot be neglected. Furthermore, the dissociation of ice mantle by cosmic-ray induced photons generates atomic hydrogen, which enhances its gas phase abundance significantly in the late stage.

We note that, at a time of about 5  ×  104 years, the abundances of several species change very quickly, and the abundances of some other species show a spike-like feature. At first sight this may resemble an erroneous behavior caused by the differential equation solver, so we ran the model with the same parameters using a Monte Carlo code (also used for benchmark purpose in Du & Parise 2011), which is free of such problems, and it turns out that these features are genuine. A semi-quantitative explanation of this feature is in Appendix A.

3.2. Chemical age versus dynamical time scale

As noted before, the time of the closest agreement between our modeling results and the observational results of Bergman et al. (2011a,b) is at  ~6  ×  105 year. Interestingly, this time scale is quite close to those derived in André et al. (2007, their Table 7). For example, the evolution time scale for ρ Oph A as estimated to be three times the free-fall time is (0.5–2)  ×  105 years, which is close to the statistically estimated age, while the collisional time scale of 5.5  ×  105 years and the cross time scale of 8  ×  105 years are also within the same order of magnitude.

It is important to define the starting point when talking about age. In the chemical evolution model described above, the whole system starts to evolve from a simple initial state where all the elements except for hydrogen are in atomic form, and the grains are bare. How relevant is such an initial condition when we talk about the dynamical evolution of a cloud condensation? At a density of  ~103 cm-3, a temperature of  ~20 K, and a visual extinction of 2, which are typical of diffuse or translucent molecular clouds, a chemical model starting from an atomic initial condition (except for H2) reaches steady state for most of the species in several 103 years. Because such a time scale is much shorter than the dynamical time scale of a cloud, it should not make much difference for a time-dependent model (i.e. one in which temperature and density etc. vary with time) to use either atomic or molecular initial conditions. For a time-independent cloud model (such as the one we are using) with constant physical conditions, adopting an atomic initial condition and a high density is equivalent to assuming that the cloud was compressed from the diffuse interstellar medium very quickly. Converging flows in the interstellar medium might play such a role, although this seems unlikely in ρ Oph A because only a very small velocity gradient was observed (André et al. 2007).

However, the time of best match and the predicted abundances may also be very dependent on the value adopted for some of the modeling parameters that are not very well constrained. It is necessary to see how the results would be changed if these parameters are varied.

3.3. Effects of changing the energy barrier of the surface reaction H + O2  →  HO2

The activation barrier of reaction (3) was set to 1200 K in Tielens & Hagen (1982), which was taken from the theoretical calculation by Melius & Blint (1979) for the gas phase case. Cuppen et al. (2010) conclude that this reaction has a negligible barrier. We choose to use an intermediate value for the energy barrier of reaction (3), namely, 600 K for the modeling. Here we study how the uncertainties in this parameter would affect the abundances of several species of interest.

In Fig. 2 the abundances of several species at the best-match time (6  ×  105 yr) are plotted as a function of the activation energy barrier of reaction (3). The temperature and density are fixed at 21 K and 6  ×  105 cm-3, and the ratio between the diffusion energy barrier and the binding energy is set to 0.77.

The abundance of gas phase O2 is not affected by changing the energy barrier of reaction (3) because its abundance is mainly determined by the gas phase production process O + OH → O2 + H and adsorption. However, the abundance of grain surface O2 increases by five orders of magnitude as the barrier changes from 0 to 1200 K, because reaction (3) is one of the main reactions for the consumption of grain O2. However, even with a barrier of 1200 K for reaction (3), the abundance of O2 on the grain at a late stage is too low to be detected. Although the rate coefficient of reaction (3) is reduced by six orders of magnitude when its barrier changes from 0 to 1200 K, the abundance of O2H on the grain does not change significantly. The reason is that when reaction (3) becomes slower, more O2 will build up on the grain, compensating for the effect of increasing the barrier of reaction (3). The abundance of gas phase H2O2 does not change much either, as its abundance is mainly determined by accretion onto the grain, and production by hydrogenation of O2H followed by partial chemical desorption.

thumbnail Fig. 2

Dependence of the abundances of several species at the best-match time (6  ×  105 yr) on the reaction barrier of reaction (3), with density fixed to 6  ×  105 cm-3 and temperature fixed to 21 K.

3.4. Effects of changing the diffusion energy barriers

The diffusion energy barrier of a species on the grain determines how fast it can migrate on the grain, so it basically determines the pace of the grain chemistry. Usually it is set to a fixed fraction (here we denote it by η) of the binding energy of each species. The latter determines how fast a species evaporates into the gas phase. However, this parameter might depend on the material and morphology of the dust grain, as well as on the property of the adsorbate itself, so it is very uncertain. Values in the range of 0.3–0.77 have been adopted in the literature. We mainly used 0.77 for our modeling. Here we investigate how different values of η would affect the abundances of several species of interest.

In Fig. 3 we plot the abundances of several species at the time of best match (6  ×  105 yr) as a function of the ratio between the diffusion barrier and binding energy (η). Changing this parameter has a strong effect on the abundances of some species. For example, the abundance of CO2 ice is reduced with a higher η. This is because at the time of concern it mainly forms through the reaction OH + CO → CO2 + H, which requires the migration of two relatively heavy species. With a low diffusion energy and a moderate temperature (~20 K), OH and CO can thermally hop quite fast, leading to a high abundance of CO2, which is not observed. However, if a lower temperature (~10 K) is adopted, the problem becomes the opposite: the mobilities of OH and CO are so low that it is difficult for them to meet each other to form enough CO2 ice, and some intricate mechanisms (e.g. three body reaction) have to be introduced to account for this (Garrod & Pauly 2011).

The abundance of H2O ice increases as η increases. At this stage it is mainly formed by hydrogenation of OH. As the mobility of atomic hydrogen on the grain is not greatly affected by the value of the diffusion energy barrier because it is allowed to migrate through quantum tunneling, the reaction rate of H + OH → H2O is not significantly affected by η; however, a larger η leaves more OH available for water because a lower amount of it is consumed in forming CO2. This is also the reason for a higher CO (and species dependent on it such as H2CO and CH3OH) ice abundance when η is larger. The abundances of gas phase H2O2 and O2H also increase for larger η, albeit only mildly.

thumbnail Fig. 3

Dependence of the abundances of several species at the best-match time (6  ×  105 yr) on the ratio between the diffusion and desorption energy, with density fixed to 6  ×  105 cm-3 and temperature fixed to 21 K.

3.5. Dependence on the temperature and density

The physical conditions (temperature and density) of ρ Oph A are determined by non-LTE radiative transfer modeling (Bergman et al. 2011a), which is usually subject to uncertainties from many aspects, such as the excitation condition, source geometry, beam filling factor, etc. In this section we study how the uncertainties in the temperature and density of the system would affect the abundances of several species in our model. In Figs. 4 and 6 we plot the abundances of several species at the time of best match (6  ×  105 yr) as a function of temperature and density.

Apparently, temperature plays a much more drastic role than density. This is intuitively easy to understand because temperature enters the calculation of rates exponentially for the surface reactions. The general trend is that when the temperature is either too low or too high, the grain surface chemistry tends to be inactive or unimportant. In the former case the mobilities of species other than atomic hydrogen (which migrates through quantum tunneling in our present model) are low, while in the latter case the surface abundances of many species are low due to elevated evaporation rates.

thumbnail Fig. 4

Dependence of the abundances of several species at the time of best match (6  ×  105 yr) on temperature, with density fixed to 6  ×  105 cm-3. Left panel: gas phase species; right panel: surface species.

As can be seen in Fig. 4, the abundance of CO ice starts to decrease at a temperature of around 20 K. This value can be estimated as the temperature at which the gas phase and grain surface abundance of CO are equal (see also Tielens & Hagen 1982), when only the adsorption and evaporation processes are taken into account (see also Hollenbach et al. 2009): Tevap=ED/ln[νnHRG1πr218kTgas/πm]ED/60+ln[(105cm-3nH)(20KTgas)1/2(m28au)1/2],\begin{eqnarray} T_{\rm evap} &=& E_{\rm D}/\ln\left[\frac{\nu}{n_{\rm H}\;{R_{\rm G}}}\frac{1}{\pi r^2}\frac{1}{\sqrt{8kT_{\rm gas}/\pi m}}\right] \nonumber \\ &&\simeq E_{\rm D}/\left\{60+ \ln\left[ \left(\frac{10^5\ {\rm cm}^{-3}}{n_{\rm H}}\right) \left(\frac{20\ {\rm K}}{T_{\rm gas}}\right)^{1/2} \left(\frac{m}{28\ {\rm au}}\right)^{1/2} \right]\right\}, \end{eqnarray}(6)where ED is the evaporation energy barrier of a species on the grain surface, ν the vibrational frequency of a species on the grain, RG the dust-to-gas number ratio, r the grain radius, while m is the molecular mass of the species being considered. A typical value of 1012 s-1 for ν, a dust grain radius of 0.1 μm, and a dust-to-gas ratio of 2.8  ×  10-12 have been adopted in deriving the number 60. Since the logarithmic part in this equation is usually small for a typical gas density and temperature, the evaporation temperature can be approximated simply by Tevap ≃ ED/60. For CO, a canonical value of ED is 1210 K (Allen & Robinson 1977), which gives an evaporation temperature of 20 K, while for water, an ED of 1860 K (Hasegawa & Herbst 1993) on bare graphite grains gives an evaporation temperature of  ~30 K. In Garrod & Herbst (2006) a much higher desorption energy of 5700 K for water is used (appropriate for water ice), which gives a evaporation temperature of  ~95 K, close to the observed evaporation temperature of water in envelopes surrounding protostars (Maret et al. 2002). The evaporation temperature of CH4 ice is close to that of CO ice. The evaporation time scale at the evaporation temperature can be estimated to be roughly ν-1exp(ED/Tevap) ≃ 10-12e60  s ≃ 3.6 × 106 yr. For higher temperature, the evaporation will be much faster.

thumbnail Fig. 5

The abundances of gas phase H2O2 and CH3OH as a function of time for different temperatures. The density is fixed to 6  ×  105 cm-3. The horizontal dotted lines mark the observed abundances of the corresponding species.

thumbnail Fig. 6

Dependence of the abundances of several species at the time of best match (6  ×  105 yr) on nH, with temperature fixed to 21 K. Left panel: gas phase species; right panel: surface species.

The abundance of CO2 ice initially increases with increasing temperature. This is because it mainly forms from the reaction CO + OH → CO2 + H with a barrier of 80 K, and an increase in temperature greatly enhances the mobility of the reacting species, as well as the probability of overcoming the reaction barrier. But when temperature increases more, the abundance of CO ice becomes so low that the CO2 abundance also drops. The evaporation temperature of CO2 is about 40 K, so evaporation is not responsible for the decline in CO2 ice abundance seen in Fig. 4. A similar trend is seen in other species, such as H2CO ice, CH3OH ice, as well as gas phase H2O2, CH3OH, H2CO, and H2O. On the other hand, for species efficiently produced in the gas phase, e.g. O2, its abundance increases with temperature owing to the faster evaporation at a higher temperature.

The abundances of gas phase H2O2, CH3OH, CO, and O2 have a sensitive dependence on temperature at around 21–24 K as seen from Fig. 4. For example, changing the temperature from 20 K to 22 K increases the abundance of H2O2 at a given time (6  ×  105 yr) by about one order of magnitude. This rather small change of 2 K in temperature is normally below the accuracy of the kinetic temperature as determined from observational data. The dependence of the evolution curves of H2O2 and CH3OH on temperature can be seen more clearly in Fig. 5. Changing the temperature not only shifts the evolution curves horizontally, it also changes their shapes significantly. We can see that, although for CH3OH it is possible to match the observed abundance at multiple stages, for H2O2 the best match is only possible at 3  ×  105–106 yr (if we ignore the match in the very early stage).

Regarding the density dependence, Fig. 6 shows that the abundances of the gas phase species at 6  ×  105 yr generally decrease with increasing density, because the accretion of molecules onto the dust grains is faster for a higher density. This does not necessarily mean that the abundances of the surface species always increase with density. For example, the abundances of CH4 ice and CH3OH ice decrease with higher density, while the abundance of H2CO ice has the opposite trend. One important factor is the time when we look at the system. In Fig. 7 we plot the abundances of H2CO ice and CH3OH ice as a function of time for different densities, while the temperature is fixed at 21 K. It can be seen that at a given time, the abundances of H2CO ice or CH3OH ice can either increase or decrease when the density is increased. We note that the evolution curves of these species have a quasi-oscillatory feature. For the same species, the evolution curves have a similar shape for different densities, except that with a lower density the evolution is slower and thus the curves are shifted toward the right.

3.6. Discussions and limits of the model

It can be seen from Figs. 1, 6, and 5 that H2O2 can be produced with a rather high abundance at certain evolutionary stages or with certain physical parameters, which is frequently higher than the currently observed value (Bergman et al. 2011b), especially in the early-to-intermediate times (~103–105 yr). One natural question to ask is why H2O2 has not been commonly detected in the interstellar medium. One possibility is that its spectral lines have been overlooked in the past. Another simple explanation would be that it is overproduced in our model, because in the past the chemistry of H2O2 (as well as O2H and O3, etc.) may not have been studied in detail in the astrochemistry context (especially the destruction reactions in the gas phase), so that its destruction routes might be incomplete. For example, there is no reaction in which H2O2 is destroyed by reacting with H+3\hbox{$_3^+$} in the current mainstream gas phase chemical networks. As a rough estimate, if we take the rate of destruction by such a reaction to be the same as the rate of the reaction H3++H2OH3O++H2\hbox{${\rm H_3^+ + H_2O\rightarrow H_3O^+ + H_2}$}, then the abundance of gas phase H2O2 can be reduced by about one order of magnitude. Further theoretical/experimental studies of the H2O2 chemistry would thus be very helpful, given that it has been detected recently and it plays a potentially important role in the grain chemistry of water. A third possibility is that the rarity of H2O2 might be an age effect. From Fig. 5 we notice that the abundance of H2O2 is only very high in a relatively early stage (before  ~5  ×  105 yr). If for certain reasons most of these cloud cores that are being observationally studied are older than this (due to some selection effects), then the H2O2 abundance in these objects would be too low to detect. Basically, at least three physical parameters are relevant, namely age, density, and temperature. A probability distribution of these three parameters of the cloud cores would help to give the detection probability of H2O2 (and any other molecules). Thus ρ Oph A may be considered special in the sense that it has a relatively high density (~106 cm-3) and temperature (20–30 K), while most dark clouds with a high density ( ≳ 104 cm-3) tend to be very cold ( ≲ 15 K) (Bergin & Tafalla 2007). An inhomogeneous physical condition would make the situation more complex, which may require a self-consistent dynamical-chemical model. However, a thorough study of these possibilities has to be left to future work.

thumbnail Fig. 7

The abundances of H2CO and CH3OH ice as a function of time for different densities. The temperature is fixed to 21 K.

On the other hand, although CH3OH has been studied quite extensively in the past, we notice that the gas phase reactions associated with it contain some important (although not decisive) differences between the OSU09 network and the UMIST RATE06 network. For example, the reaction between CH and CH3OH to form CH3 and H2CO has a rate 2.49  ×  10-10(T/300)-1.93 in the UMIST RATE06 network, but it does not exist in the OSU09 network. One possible problem is that the temperature of cold interstellar medium (at most several tens of Kelvins) is out of the indicated valid range for many reactions in RATE06, and it is not clear how to extrapolate these reaction rates correctly, although we have closely followed the instructions in Woodall et al. (2007). In our modeling we have been using the RATE06 network.

The energy barriers for the hydrogenation of CO and H2CO on the grain are both taken to be 2500 K. Woon (2002) calculated the barrier heights of these two reactions, giving a value of  ~2740 K and 3100 K in the case three water molecules are present, with zero-point energy corrections added. If these values are adopted, then the observed abundances of H2CO and CH3OH can only be reproduced within one order of magnitude at best. However, Goumans (2011) gives a lower barrier height (~2200 K) for the hydrogenation of H2CO, which would give better agreement with the observational results than if Woon (2002) were used in our model.

The chemical desorption is very important for the abundances of the gas phase H2O2 and CH3OH. However, the efficiency of this mechanism (the “a” parameter in Garrod et al. 2007) is uncertain. Garrod & Pauly (2011) adopt a low value of 0.01 for it, to avoid overproduction of some gas phase species (Garrod et al. 2007). We use a value of 0.1 in our study, and the abundances of H2CO and CH3OH are not overproduced, except possibly in the early stages of the evolution. We note that the temperature and density of major concern in our study is around 20 K and 6  ×  105 cm-3, while in Garrod et al. (2007) the temperature is set to 10 K and density to 2  ×  104 cm-3. As a test we ran our model with the latter physical condition, and in this case CH3OH and H2CO are indeed overproduced by about one order of magnitude. Thus it seems that the efficiency of chemical desorption depends on the temperature of the dust grain, in the sense that at higher temperature the probability that the product of a surface exothermic reaction gets ejected to the gas phase is higher. However, a detailed study of this possibility is beyond the scope of the present paper.

Regarding the formation of water ice, Bergin et al. (1998) propose an interesting mechanism in which water is first formed in the high-temperature shocked gas, and then gets adsorbed onto the dust grains in the post-shock phase. This mechanism may act as an alternative or supplement to the grain chemistry route. It is not our aim here to discuss to what extent this mechanism contributes to the water ice budget. However, we point out that even with a temperature 1000–2000 K, the amount of H2O2 produced in a pure gas phase chemistry (using the UMIST RATE06 network) is still much lower than the detected level.

4. Conclusions

  • With a gas-grain chemical model that properly takes thedesorption of grain surface species by the heat released bychemical reactions into account, we reproduced the observedabundances of H2O2 in ρ Oph A at a time of  ~6  ×  105 yr. The solid phase H2O2 abundance is very low at this stage. However, a H2O2 to H2O ratio of a few percent might be obtained in the solid phase if the layered structure of grain mantle is considered.

  • The abundances of other species such as H2CO, CH3OH, and O2 detected in the same object can also be reasonably reproduced at a time of  ~6  ×  105 yr. Such a time scale is consistent with the evolution time scale estimated through dynamical considerations.

  • The O2H radical is a precursor of H2O2 on the dust grain, and we predict that it has a gas-phase abundance with the same order-of-magnitude of H2O2 and should thus be detectable. Observational searches for it are under way.

  • For physical conditions relevant to ρ Oph A, water is mainly in solid form, being the dominant grain mantle material. Its gas phase abundance is only on the order of 10-8 according to our model.

  • We note that the abundance of gas-phase H2O2 in our model results can be much higher than the current observed level for a range of physical conditions. This may suggest that its gas-phase destroying channels are incomplete. Due to the potentially important role played by H2O2 in the formation of water, its reaction network needs to be studied more thoroughly in the future. Other uncertainties in our modeling include the ratio between the diffusion energy barrier to the binding energy of a species on the grain surface, the activation energy barriers of certain key reactions, and the efficiency of the chemical desorption mechanism. In the present work we mainly make use of their canonical values or of values that give good matches to the observational results, and we also vary them to see the effects on the resulting abundances, which are significant in many cases.

Online material

Appendix A: An explanation of the spike-like features in the evolution curves

In Fig. 1 we note that at a time of  ~ 5  ×  104 yr a spike-like feature appears in the evolution curves of some species (e.g. gas phase H2O2 and CH3OH), while the abundances of some other species (e.g. H2O2 and O2 on the grain) change very rapidly at the same time. These features might appear to be caused by faults in the program for solving the set of differential equations. However, after solving the same problem with the Monte Carlo method, which is immune to such numerical instabilities, we find that these features are still present, indicating that they are genuine. How could a smooth ordinary differential equation system generate such an almost-singular feature? In Fig. A.1 we make a zoom-in of Fig. 1 (with several species added and several species removed). It can be seen that, although the time scale of the spike-like feature is relatively short, the evolution is always smooth (except the discreteness caused by the finite sampling of the curve). Then what determines the appearance of such a feature?

thumbnail Fig. A.1

A zoom-in plot of Fig. 1 (with several species removed and several added), focusing on the spike-like feature (though nothing looks like a spike anymore). The abundance of gH has been multiplied by a factor of 1012 (to make the plot look more compact).

Because atomic hydrogen is central to the surface chemistry, we first look at the most important reactions governing its abundance on the grain. The main reactions consuming atomic hydrogen on the grains are (in the following a species name preceded by a “g” means a species on the grain, otherwise it is in the gas phase), gH+gO2H0k1{gOH+gOH,gH2O2,OH+OH,H2O2,gH+gH2O20k2{gH2O+gOH,H2O+OH,\appendix \setcounter{section}{1} \begin{eqnarray} {\rm gH} &+& {\rm gO_2H} \xrightarrow{k_1} \left\{ \begin{aligned} & {\rm gOH + gOH,} \\ & {\rm gH_2O_2,}\\ & {\rm OH + OH,} \\ & {\rm H_2O_2,} \end{aligned} \right. \\ {\rm gH} &+& {\rm gH_2O_2} \xrightarrow{k_2} \left\{ \begin{aligned} & {\rm gH_2O + gOH},\\ & {\rm H_2O + OH}, \end{aligned} \right. \end{eqnarray}while others such as the reactions with gO3 and gCH2OH are relatively unimportant in the early times. It is mainly produced by H0k3gH,gOH+gCO0k4gCO2+gH,gH2O0k5gOH+gH.\appendix \setcounter{section}{1} \begin{eqnarray} && {\rm H \xrightarrow{k_3} gH,} \\ &&{\rm gOH + gCO \xrightarrow{k_4} gCO_2 + gH,} \\ &&{\rm gH_2O \xrightarrow{k_5} gOH + gH.} \end{eqnarray}As the abundance of atomic H in the gas phase and the abundance of gCO and gH2O change relatively smoothly, they can be viewed as constant on a short time scale. The main reactions for the consumption and production of gOH are gH+gH2O20k9gH2O+gOH,gH2O0k10gH+gOH.\appendix \setcounter{section}{1} \begin{eqnarray} && {\rm gO + gOH} \xrightarrow{k_6} \left\{ \begin{aligned} &{\rm gO_2H,} \\ &{\rm O_2H,} \\ \end{aligned} \right. \label{eq:gOgOH} \\ &&{\rm gCO + gOH} \xrightarrow{k_7} \left\{ \begin{aligned} &{\rm gCO_2 + gH,} \\ &{\rm CO_2 + H,} \\ \end{aligned} \right. \label{eq:gCOgOH} \\ && {\rm gH + gO_2H} \xrightarrow{k_8} {\rm gOH + gOH,} \label{eq:gHgO2H} \\ &&{\rm gH + gH_2O_2} \xrightarrow{k_9} {\rm gH_2O + gOH}, \\ &&{\rm gH_2O} \xrightarrow{k_{10}} {\rm gH+gOH}. \end{eqnarray}From the above reaction list we may write the evolution equation of gH and gOH as4tgH=k1gHgO2Hk2gHgH2O2\appendix \setcounter{section}{1} \begin{eqnarray} \partial_t{\rm gH} &= & -k_1 {\rm gH gO_2H} - k_2 {\rm gH gH_2O_2} \nonumber \\ & & +k_3 {\rm H} + k_4 {\rm gOH gCO} + k_5 {\rm gH_2O}, \label{eq:gHEq} \end{eqnarray}(A.11)tgOH=k6gOgOHk7gCOgOH\appendix \setcounter{section}{1} \begin{eqnarray} \partial_t{\rm gOH}& = & -k_6 {\rm gO gOH} - k_7 {\rm gCO gOH} \nonumber \\ & & +2k_8 {\rm gH gO_2H} + k_9 {\rm gH gH_2O_2} + k_{10}{\rm gH_2O}, \label{eq:gOHEq} \end{eqnarray}(A.12)or in a more succinct form tgH=κ1gH+κ2gOH+b1,\appendix \setcounter{section}{1} \begin{eqnarray} &&\partial_t{\rm gH} = \kappa_1 {\rm gH} + \kappa_2 {\rm gOH} + b_1, \nonumber \\ &&\partial_t{\rm gOH} = \kappa_3 {\rm gH} + \kappa_4 {\rm gOH} + b_2, \label{eq:gHgOH} \end{eqnarray}(A.13)where κ1=k1gO2Hk2gH2O2,κ2=k4gCO,κ3=2k8gO2H+k9gH2O2,κ4=k6gOk7gCO,b1=k3H+k5gH2O,b2=k10gH2O.\appendix \setcounter{section}{1} \begin{eqnarray*} \kappa_1 &=& -k_1{\rm gO_2H} - k_2{\rm gH_2O_2},\\ \kappa_2 &=& k_4{\rm gCO},\\ \kappa_3 &=& 2k_8{\rm gO_2H}+k_9{\rm gH_2O_2},\\ \kappa_4 &=& -k_6{\rm gO}-k_7{\rm gCO},\\ b_1 &=& k_3{\rm H} + k_5{\rm gH_2O},\\ b_2 &=& k_{10}{\rm gH_2O}. \end{eqnarray*}If we view κ1κ4, as well as b1 and b2, as constants (of course they are not), then Eq. (A.13) can be solved exactly; the solution contains an exponential part and a constant part. The amplitude of the exponential part will be inversely proportional to the determinant of the coefficient matrix det(κ)=().\appendix \setcounter{section}{1} \begin{eqnarray*} {\rm det}(\kappa) = \left( \begin{array}{cc} \kappa_1 & \kappa_2\\ \kappa_3 & \kappa_4 \end{array} \right). \end{eqnarray*}Since κ1κ4 are not really constant in our problem, we expect that when they become such that det(κ) is close to zero, a spike-like or jump-like behavior would appear. Namely, we require κ1κ4κ2κ3=(k1gO2H+k2gH2O2)(k6gO+k7gCO)(k4gCO)(2k8gO2H+k9gH2O2)gO2H/gH2O2+11.2gO2H/gH2O2+1×(2.5×108gO/gCO+1)\appendix \setcounter{section}{1} \begin{eqnarray} \frac{\kappa_1\kappa_4}{\kappa_2\kappa_3} &=&\frac{(k_1{\rm gO_2H} + k_2{\rm gH_2O_2})(k_6{\rm gO}+k_7{\rm gCO})} {(k_4{\rm gCO})(2k_8{\rm gO_2H}+k_9{\rm gH_2O_2})} \nonumber \\ &\simeq&\frac{{\rm gO_2H/gH_2O_2 + 1}}{1.2{\rm gO_2H/gH_2O_2}+1} \times \left(2.5\times10^8{\rm gO/gCO + 1}\right) \nonumber \\ &\simeq& 1. \label{eq:spikecondition} \end{eqnarray}(A.14)The actual value of the parameters have been inserted in the second line of the above equation. These parameters depend on the physical conditions.

To satisfy this condition (at least approximately), gO/gCO should be very small. In Fig. A.2 the ratio (κ1κ4)/(κ2κ3) and the value of 2.5  ×  108 gO/gCO are plotted as a function of time. The abundances of O, gO, gO2, gH, and gOH are also plotted for reference (not to scale). It can be seen that the (κ1κ4)/(κ2κ3) ratio does decrease, and it approaches a value of unity before the time of the spike/jump feature, and the gO/gCO ratio drops to a very low value monotonically.

thumbnail Fig. A.2

The ratio (κ1κ4)/(κ2κ3) and the value of 2.5  ×  108 × gO/gCO as a function of time. They are to scale. The abundances of several other species are plotted only for reference (not to scale).

The above mathematical argument can also be understood intuitively. When the gO abundance is so low that reaction (A.6) can be neglected, gOH is only destroyed by reaction (A.7). Each time a gOH radical is consumed, one gH is created (if we neglect the desorption process), and this gH will quickly react with gH2O2 or gO2H to create one or two gOHs. There will then be a net gain in the gOH abundance, leading to its fast growth, and the gH abundance will increase accordingly. Thus we see that reaction (A.8) is crucial in that it produced two gOH radicals by consuming only one gH.

Therefore for such a spike-like feature to occur, the abundance of gO must decrease to a low value such that the reaction between gOH and gO becomes unimportant for consuming

gOH. Namely, we require  [ gO ]  ≲  [ gCO ]   k7/k6 ≃ 5 × 10-13. The abundances of atomic oxygen in the gas phase and on the grain surface are related by O/gO ≃ kevap(gO)/kad(O) ≃ 3 × 106, so equivalently we require  [ O ]  < 1.5 × 10-6 (i.e., about a factor of 200 less than the initial O abundance) at the time of the spike-like feature. Atomic oxygen is mainly consumed on the grain surface by reacting with another gO to form gO2 (for t ≲ 103 yr) or by reacting with gOH to form gO2H. Only the latter is relevant here. As the abundance of gOH does not change much before the spike-like feature, the time scale for the consumption of atomic oxygen can be estimated to be ln200 × (O/gO)/(k6gOH) ≃ 105 yr, which is on the same order of magnitude as the time of occurrence of the spike-like feature.

The time scale for the endurance of the spike-like feature itself can be estimated to be the time scale for the exhaustion of gO2H or gH2O2 (so that Eqs. (A.11) and (A.12) do not hold anymore) by reacting with gH, which is 1/(k1gH)103yr,\begin{equation*} 1/(k_1{\rm gH})\simeq 10^3\ {\rm yr}, \end{equation*}where the gH population (the average number of atomic H on a single grain) is taken to be a median value (10-14) during the rapidly varying period. That equation (A.14) seems to hold after this period (see Fig. A.2) does not mean that gH will keep increasing rapidly, simply because the premise of our argument, namely Eqs. (A.11) and (A.12) are not a good description of the evolution of gH and gOH anymore.

As the abundance of H on the grain increases, almost all the O3 on the grain are converted into O2 and OH. O2 molecules on the grain are then consumed by the slower reaction gH + gO2 → gO2H, with a time scale  ~ 105 yr. This explains the prominent peak in the evolution curve of gO2 (see Fig. 1).

Would the spike-like feature have any practical significance (especially observationally)? Ideally, such a short-time feature may be used to constrain the age of a dense cloud, by distinguishing the abundances of certain species between their early-time values (before the spike-like feature) and late-time values (after the spike-like feature). However, due to its dependence on the reaction network being used, which usually contains a lot of uncertainties and is subject to change when new experiments are carried out, the question whether this feature is really relevant for the study of interstellar medium can only be answered by future investigations.

Appendix B: The surface reaction network used in this work

Table B.1

Surface network used in this work.

Appendix C: The enthalpies of the surface species

Table C.1

Enthalpies of the surface species considered in this work that are used to calculate the exoergicities of the surface reactions.


4

Here for brevity we use the name of a species to denote its average population on a single grain; for example, if gCO = 100, it mains on average there are 100 CO molecules on a single grain.

Acknowledgments

We thank the anonymous referee for very detailed and constructive comments. We also thank A. G. G. M. Tielens for interesting discussions. F. Du and B. Parise are financially supported by the Deutsche Forschungsgemeinschaft Emmy Noether program under grant PA1692/1-1.

References

  1. Allen, M., & Robinson, G. W. 1977, ApJ, 212, 396 [NASA ADS] [CrossRef] [Google Scholar]
  2. An, D., Ramírez, S. V., Sellgren, K., et al. 2011, ApJ, 736, 133 [NASA ADS] [CrossRef] [Google Scholar]
  3. André, P., Belloche, A., Motte, F., & Peretto, N. 2007, A&A, 472, 519 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  5. Atkinson, R., Baulch, D. L., Cox, R. A., et al. 2004, Atmos. Chem. & Phys., 4, 1461 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bergin, E. A., & Tafalla, M. 2007, ARA&A, 45, 339 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bergin, E. A., Neufeld, D. A., & Melnick, G. J. 1998, ApJ, 499, 777 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  8. Bergin, E. A., Melnick, G. J., Stauffer, J. R., et al. 2000, ApJ, 539, L129 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bergman, P., Parise, B., Liseau, R., & Larsson, B. 2011a, A&A, 527, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bergman, P., Parise, B., Liseau, R., et al. 2011b, A&A, 531, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Binnewies, M., & Milke, E. 1999, Thermochemical data of elements and compounds (Wiley-VCH) [Google Scholar]
  12. Boudin, N., Schutte, W. A., & Greenberg, J. M. 1998, A&A, 331, 749 [NASA ADS] [Google Scholar]
  13. Cazaux, S., Cobut, V., Marseille, M., Spaans, M., & Caselli, P. 2010, A&A, 522, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Cuppen, H., & Herbst, E. 2007, ApJ, 668, 294 [Google Scholar]
  15. Cuppen, H. M., Ioppolo, S., Romanzin, C., & Linnartz, H. 2010, Phys. Chem. Chem. Phys. (Incorporating Faraday Transactions), 12, 12077 [NASA ADS] [CrossRef] [Google Scholar]
  16. Du, F., & Parise, B. 2011, A&A, 530, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Fuchs, G. W., Cuppen, H. M., Ioppolo, S., et al. 2009, A&A, 505, 629 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Garrod, R. 2008, A&A, 491, 239 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Garrod, R. T., & Herbst, E. 2006, A&A, 457, 927 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Garrod, R. T., & Pauly, T. 2011, ApJ, 735, 15 [NASA ADS] [CrossRef] [Google Scholar]
  21. Garrod, R., Park, I. H., Caselli, P., & Herbst, E. 2006, Faraday Discussions, 133, 51 [Google Scholar]
  22. Garrod, R., Wakelam, V., & Herbst, E. 2007, A&A, 467, 1103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Garrod, R. T., Weaver, S. L. W., & Herbst, E. 2008, ApJ, 682, 283 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A., & Tielens, A. G. G. M. 2004, ApJS, 151, 35 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  25. Goicoechea, J. R., Cernicharo, J., Lerate, M. R., et al. 2006, ApJ, 641, L49 [NASA ADS] [CrossRef] [Google Scholar]
  26. Goldsmith, P. F., Liseau, R., Bell, T. A., et al. 2011, ApJ, 737, 96 [Google Scholar]
  27. Gomis, O., Leto, G., & Strazzulla, G. 2004, A&A, 420, 405 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Goumans, T. P. M. 2011, MNRAS, 413, 2615 [NASA ADS] [CrossRef] [Google Scholar]
  29. Goumans, T. P. M., Uppal, M. A., & Brown, W. A. 2008, MNRAS, 384, 1158 [NASA ADS] [CrossRef] [Google Scholar]
  30. Hasegawa, T. I., & Herbst, E. 1993, MNRAS, 261, 83 [NASA ADS] [CrossRef] [Google Scholar]
  31. Hasegawa, T., Herbst, E., & Leung, C. M. 1992, ApJS, 82, 167 [NASA ADS] [CrossRef] [Google Scholar]
  32. Herbst, E., & Leung, C. M. 1989, ApJS, 69, 271 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Hollenbach, D., Kaufman, M. J., Bergin, E. A., & Melnick, G. J. 2009, ApJ, 690, 1497 [CrossRef] [Google Scholar]
  34. Ioppolo, S., Cuppen, H. M., Romanzin, C., van Dishoeck, E. F., & Linnartz, H. 2008, ApJ, 686, 1474 [NASA ADS] [CrossRef] [Google Scholar]
  35. Ioppolo, S., Cuppen, H. M., Romanzin, C., van Dishoeck, E. F., & Linnartz, H. 2010, Phys. Chem. Chem. Phys. (Incorporating Faraday Transactions), 12, 12065 [NASA ADS] [CrossRef] [Google Scholar]
  36. Kaiser, R. I., Ochsenfeld, C., Head-Gordon, M., & Lee, Y. T. 1998, Science, 279, 1181 [NASA ADS] [CrossRef] [Google Scholar]
  37. Katz, N., Furman, I., Biham, O., Pirronello, V., & Vidali, G. 1999, ApJ, 522, 305 [NASA ADS] [CrossRef] [Google Scholar]
  38. Larsson, B., Liseau, R., Pagani, L., et al. 2007, A&A, 466, 999 [Google Scholar]
  39. Maret, S., Ceccarelli, C., Caux, E., Tielens, A. G. G. M., & Castets, A. 2002, A&A, 395, 573 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Melius, C. F., & Blint, R. J. 1979, Chem. Phys. Lett., 64, 183 [NASA ADS] [CrossRef] [Google Scholar]
  41. Millar, T. J., & Herbst, E. 1990, MNRAS, 242, 92 [NASA ADS] [Google Scholar]
  42. Nagy, B., Csontos, J., Kallay, M., & Tasi, G. 2010, J. Phys. Chem. A, 114, 13213 [CrossRef] [Google Scholar]
  43. Öberg, K. I., Fuchs, G. W., Awad, Z., et al. 2007, ApJ, 662, L23 [NASA ADS] [CrossRef] [Google Scholar]
  44. Öberg, K. I., Boogert, A. C. A., Pontoppidan, K. M., et al. 2011, ApJ, 740, 109 [NASA ADS] [CrossRef] [Google Scholar]
  45. Pontoppidan, K. M., van Dishoeck, E. F., & Dartois, E. 2004, A&A, 426, 925 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Roberts, H., & Herbst, E. 2002, A&A, 395, 233 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Ruffle, D. P., & Herbst, E. 2000, MNRAS, 319, 837 [NASA ADS] [CrossRef] [Google Scholar]
  48. Ruffle, D. P., & Herbst, E. 2001a, MNRAS, 322, 770 [NASA ADS] [CrossRef] [Google Scholar]
  49. Ruffle, D. P., & Herbst, E. 2001b, MNRAS, 324, 1054 [NASA ADS] [CrossRef] [Google Scholar]
  50. Savage, B. D., & Sembach, K. R. 1996, ARA&A, 34, 279 [NASA ADS] [CrossRef] [Google Scholar]
  51. Shi, J., Raut, U., Kim, J.-H., Loeffler, M., & Baragiola, R. A. 2011, ApJ, 738, L3 [NASA ADS] [CrossRef] [Google Scholar]
  52. Tielens, A., & Hagen, W. 1982, A&A, 114, 245 [NASA ADS] [Google Scholar]
  53. Tielens, A. G. G. M., McKee, C. F., Seab, C. G., & Hollenbach, D. J. 1994, ApJ, 431, 321 [NASA ADS] [CrossRef] [Google Scholar]
  54. van Dishoeck, E. F. 2004, ARA&A, 42, 119 [NASA ADS] [CrossRef] [Google Scholar]
  55. Vandooren, J., Sarkisov, O. M., Balakhnin, V. P., & van Tiggelen, P. J. 1991, Chem. Phys. Lett., 184, 294 [NASA ADS] [CrossRef] [Google Scholar]
  56. Vasyunin, A., Semenov, D., Wiebe, D., & Henning, T. 2009, ApJ, 691, 1459 [NASA ADS] [CrossRef] [Google Scholar]
  57. Wakelam, V., Herbst, E., & Selsis, F. 2006, A&A, 451, 551 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Wampfler, S. F., Bruderer, S., Kristensen, L. E., et al. 2011, A&A, 531, L16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Watson, W. D., & Salpeter, E. E. 1972, ApJ, 174, 321 [NASA ADS] [CrossRef] [Google Scholar]
  60. Woodall, J., Agúndez, M., Markwick-Kemper, A., & Millar, T. 2007, A&A, 466, 1197 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Woon, D. E. 2002, ApJ, 569, 541 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table B.1

Surface network used in this work.

Table C.1

Enthalpies of the surface species considered in this work that are used to calculate the exoergicities of the surface reactions.

All Figures

thumbnail Fig. 1

The time evolution of the abundances of selected species. A prefix “g” means a grain surface species (right panel), while a name without a prefix means a gas phase species (left panel). The dotted lines are the observed abundances (relative to H2) of O2 (blue), H2CO (cyan), CH3OH (magenta), and H2O2 (black), respectively, in the ρ Oph A source.

In the text
thumbnail Fig. 2

Dependence of the abundances of several species at the best-match time (6  ×  105 yr) on the reaction barrier of reaction (3), with density fixed to 6  ×  105 cm-3 and temperature fixed to 21 K.

In the text
thumbnail Fig. 3

Dependence of the abundances of several species at the best-match time (6  ×  105 yr) on the ratio between the diffusion and desorption energy, with density fixed to 6  ×  105 cm-3 and temperature fixed to 21 K.

In the text
thumbnail Fig. 4

Dependence of the abundances of several species at the time of best match (6  ×  105 yr) on temperature, with density fixed to 6  ×  105 cm-3. Left panel: gas phase species; right panel: surface species.

In the text
thumbnail Fig. 5

The abundances of gas phase H2O2 and CH3OH as a function of time for different temperatures. The density is fixed to 6  ×  105 cm-3. The horizontal dotted lines mark the observed abundances of the corresponding species.

In the text
thumbnail Fig. 6

Dependence of the abundances of several species at the time of best match (6  ×  105 yr) on nH, with temperature fixed to 21 K. Left panel: gas phase species; right panel: surface species.

In the text
thumbnail Fig. 7

The abundances of H2CO and CH3OH ice as a function of time for different densities. The temperature is fixed to 21 K.

In the text
thumbnail Fig. A.1

A zoom-in plot of Fig. 1 (with several species removed and several added), focusing on the spike-like feature (though nothing looks like a spike anymore). The abundance of gH has been multiplied by a factor of 1012 (to make the plot look more compact).

In the text
thumbnail Fig. A.2

The ratio (κ1κ4)/(κ2κ3) and the value of 2.5  ×  108 × gO/gCO as a function of time. They are to scale. The abundances of several other species are plotted only for reference (not to scale).

In the text

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.