Issue |
A&A
Volume 505, Number 2, October II 2009
|
|
---|---|---|
Page(s) | 629 - 639 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200810784 | |
Published online | 22 July 2009 |
Hydrogenation reactions in interstellar CO ice analogues
A combined experimental/theoretical approach
G. W. Fuchs - H. M. Cuppen - S. Ioppolo - C. Romanzin - S. E. Bisschop - S. Andersson - E. F. van Dishoeck - H. Linnartz
Raymond & Beverly Sackler Laboratory for Astrophysics, Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
Received 11 August 2008 / Accepted 11 June 2009
Abstract
Context. Hydrogenation reactions of CO in inter- and circumstellar ices are regarded as an important starting point in the formation of more complex species. Previous laboratory measurements by two groups of the hydrogenation of CO ices provided controversial results about the formation rate of methanol.
Aims. Our aim is to resolve this controversy by an independent investigation of the reaction scheme for a range of H-atom fluxes and different ice temperatures and thicknesses. To fully understand the laboratory data, the results are interpreted theoretically by means of continuous-time, random-walk Monte Carlo simulations.
Methods. Reaction rates are determined by using a state-of-the-art ultra high vacuum experimental setup to bombard an interstellar CO ice analog with H atoms at room temperature. The reaction of CO + H into H2CO and subsequently CH3OH is monitored by a Fourier transform infrared spectrometer in a reflection absorption mode. In addition, after each completed measurement, a temperature programmed desorption experiment is performed to identify the produced species according to their mass spectra and to determine their abundance. Different H-atom fluxes, morphologies, and ice thicknesses are tested. The experimental results are interpreted using Monte Carlo simulations. This technique takes into account the layered structure of CO ice.
Results. The formation of both formaldehyde and methanol via CO hydrogenation is confirmed at low temperature (
T = 12-20 K). We confirm that the discrepancy between the two Japanese studies is caused mainly by a difference in the applied hydrogen atom flux, as proposed by Hidaka and coworkers. The production rate of formaldehyde is found to decrease and the penetration column to increase with temperature. Temperature-dependent reaction barriers and diffusion rates are inferred using a Monte Carlo physical chemical model. The model is extended to interstellar conditions to compare with observational H2CO/CH3OH data.
Key words: astrochemistry - infrared: ISM - ISM: atoms - ISM: molecules - methods: laboratory
1 Introduction
An increasing number of experimental and theoretical studies have focused on the characterisation of solid state astrochemical processes. These studies were triggered by the recognition that many simple and more complex molecules in the interstellar medium are most likely to have formed on the surfaces of dust grains. Astronomical observations, detailed laboratory studies, and progress in UHV surface techniques allow experimental verification of the initial surface reaction schemes, as introduced by Tielens, Hagen, and Charnley (Tielens & Charnley 1997; Tielens & Hagen 1982). Recently the formation of water was demonstrated in hydrogenation schemes starting from solid molecular oxygen (Ioppolo et al. 2008; Miyauchi et al. 2008) and the formation of ethanol from acetaldehyde (Bisschop et al. 2007a). The first solid-state astrochemical laboratory studies focused on the formation of formaldehyde and methanol by H-atom bombardment of CO ice. Methanol is observed abundantly in interstellar ices and is considered to be a resource for the formation of more complex molecules through surface reactions and after evaporation in the gas phase (Charnley et al. 1992). The hydrogenation scheme for the solid state formation of methanol was proposed to be![]() |
(1) |
Laboratory studies of H-atom bombardment of CO ice were performed independently by two groups (Watanabe & Kouchi 2002; Hiraoka et al. 2002). Hiraoka et al. (2002) observed only formaldehyde formation, whereas Watanabe & Kouchi (2002) also found efficient methanol production. In a series of papers, these conflicting results were discussed (Hiraoka et al. 2002; Watanabe et al. 2004,2003) and the prevailing discrepancy between results was proposed to be a consequence of different experimental conditions, most noticeable the adopted H-atom flux (Hidaka et al. 2004). Understanding the solid-state formation route to methanol became even more important following an experimental finding that the gas-phase formation route via ion-neutral reactions is less efficient than previously estimated and cannot explain the observed interstellar abundances (Geppert et al. 2005; Garrod et al. 2006).
![]() |
Figure 1: Schematic representation of the experimental setup. CO ice is deposited through the dosing line and the products are monitored by means of infrared spectroscopy and quadrupole mass spectroscopy (QMS). |
Open with DEXTER |
The present paper strongly supports the flux argument given by Hidaka et al. (2004). Furthermore, we present a systematic study of the physical dependencies involved in the CO-ice hydrogenation to place previous work in a context that allows an extension of solid state astrochemical processes to more complex species. Special emphasis is placed on the flux and temperature dependence of the formation rate. An analysis of the spectral changes of CO ice during hydrogenation is included to provide insight into the structure of the reactive layer. Furthermore, Monte Carlo simulations are presented that allow us to interpret the experimental results in greater detail and to vary parameters that are difficult to study independently by experiment. We conclude with a simulation of H2CO/CH3OH formation under interstellar conditions, in particular for low H-atom fluxes. The outcome is compared with astronomical observations.
2 Experimental procedure
The experiments are performed under UHV conditions. The room
temperature base pressure of the vacuum system is better than
mbar. Figure 1 shows a schematic
representation of the setup (see Ioppolo et al. (2008) for additional
information). Amorphous CO ices ranging from a few to several
monolayers are grown on a gold coated copper substrate that is
located at the centre of the main vacuum chamber and mounted on the
tip of a cold finger of a 10 K He-cryostat. The temperature of the
ice can be controlled between 12.0 K and 300 K with 0.5 K relative
precision between experiments. The absolute accuracy is better than 2 K. During deposition, the layer thickness is monitored by
simultaneous recording of reflection absorption infrared (RAIR)
spectra. To exclude the effects of potential pollution, ices are
grown using CO, 13CO, or C18O isotopologues.
During the experiment the ice layers are exposed to a hydrogen atom
beam. The atoms are produced by a well characterised commercial
thermal-cracking source (Tschersich 2000; Tschersich & von Bonin 1998) that
provides H-atom fluxes on the sample surface of between 1012and 1014 atoms cm-2 s-1. By comparison, the Hiraoka
group used fluxes below 1013 atoms cm-2 s-1 and the
Watanabe group worked in the 10
14-1015 atoms cm-2 s-1 regime. The hot (2000 K) hydrogen atoms are
cooled to room temperature by surface collisions in a nose-like
shaped quartz pipe between the atomic source and the ice sample. In
this way, hot hydrogen atoms cannot affect the ice directly. H-atom
recombination in this connecting pipe results in a lower final flux.
Details about the flux determination are given in Appendix A. The
absolute fluxes are estimated to be within a factor of two and the
relative fluxes to within 50%.
The relatively high temperature of the incident atoms of 300 K does not affect the process; previous experiments with colder H atoms did not show any substantial temperature dependence because the atoms are immediately thermalized on the surface (Watanabe & Kouchi 2002). It is argued that the surface is covered with a thin layer of hydrogen molecules under these conditions. These molecules are either formed on the surface or originate from the partially dissociated beam. Since the incoming atoms have to penetrate this cold H2 layer, they are thermally adjusted to the surface temperature once they come in contact with the CO molecules.
Information about the reaction products is obtained using two complementary techniques. During the H-atom bombardment, reactants and products are monitored by recording RAIR spectra. The RAIR spectra are recorded using a Fourier transform infrared spectrometer with 1 and 4 cm-1 resolution and covering the spectral region in which CO (2143 (s) cm-1), formaldehyde (1732 (s), 1479, and 2812 (m), and 1246, 1175, 2991, 2880, and 2812 (mw) cm-1) and methanol (1035 (s) and 1125 (w) cm-1) exhibit strong (s), medium (m), or weak (w) absorption bands. The intensity of spectral features is directly related to the density in the ice. The products are monitored by mass spectrometry using temperature programmed desorption (TPD) once a hydrogenation experiment is completed.
3 Experimental results
3.1 A sample experiment
To illustrate the experimental method, we start by discussing a
sample experiment in which a CO ice of
molecules cm-2 is bombarded with H atoms with a flux of
cm-2 s-1 for three hours at a surface
temperature of 12.0 K. This corresponds to a fluence of
cm-2. Figure 2 shows the
RAIR difference spectrum (
)
after these three hours
of exposure (after - before). The CO, the H2CO and the CH3OH
spectral signatures are indicated with respect to the spectrum
recorded before the H-atom bombardment started. The CO appears as a
negative band, indicating its use-up, and the other bands are
positive, indicating the formation of H2CO and CH3OH. Neither
the intermediate species, HCO and H3CO, nor more complex species
are observed.
The
column density NX (molecules cm-2) of
species X in the ice is calculated using
where




![]() |
(3) |
for a set of different experiments. In addition, the results discussed in the present paper are all in a regime where the proportionality relation (Teolis et al. 2007) still holds (<

![]() |
Figure 2:
RAIR difference spectrum of a CO ice at 12.0 K exposed to
|
Open with DEXTER |
The CO band shape can change when molecules other than CO are formed. Figure 4 shows the 2143 cm-1 IR peak before and after the H-atom exposure. A clear decrease in the peak height can be observed caused by the use-up of CO during the experiment, as expected. However, an additional peak appears at 2135 cm-1 (see inset Fig. 4), which is due to a CH3OH-CO ice interaction. Transmission IR spectra of a CH3OH:CO mixture show a band at 2136 cm-1(Palumbo & Strazzulla 1993; Bisschop 2007). When the methanol bands grow, the band at 2135 cm-1 also increases. Figure 5 shows how the peak position of CO shifts with the methanol content in the reflection spectra. The RAIR spectra on which this graph is based, are taken for ice layers that are formed by co-deposition of CO and CH3OH of known ratio. The CO stretching mode in H2O:CO and NH3:CO mixtures shows a similar behaviour (Sandford et al. 1988; Bouwman et al. 2007). As for both H2O and NH3, CH3OH is able to form hydrogen bonds and these hydrogen bonds most likely cause the redshift of the CO band. By comparing the position of the peak in Fig. 4 at 2135 cm-1 to the data in Fig. 5, we can constrain the methanol fraction in the top layers, assuming that the formed CH3OH:CO mixture has the same spectral behaviour as the deposited mixtures. The observed data after three hours correspond to a CH3OH:CO mixture of at least 90%. This means that the top layer of the ice is completely converted into H2CO and CH3OH and that no or very little additional mixing with CO occurs. For the H2CO and CH3OH band, no spectral changes are observed during the experiments.
![]() |
Figure 3: The decrease in integrated absorbance of CO and CH3OH (1035 cm-1) following desorption at a constant temperature of 29 and 135 K, respectively. The arrows indicate the transition points from the multi- to sub-monolayer regime. |
Open with DEXTER |
![]() |
Figure 4: Spectral change of the CO 2143 cm-1 RAIR band before and after H-atom bombardment. The inset shows the corresponding difference spectrum. |
Open with DEXTER |
To quantify the use-up of CO and the formation of new products, we have to assume that the apparent absorption band strength is constant during an experiment, i.e., independent of the ice composition. Bouwman et al. (2007) found that the band strength of the 2143 cm-1 CO feature is indeed unaffected within the experimental error by the water content in H2O:CO-ice mixtures up to 4:1. The band strength is expected to behave similarly for a CO:CH3OH-mixture. Furthermore, if the band strength is strongly affected by the ice composition, the total ice thickness determined using a constant band strength should vary in time, whereas the real thickness is constant. Since this does not occur, we estimate that the change in band strength caused by the change in ice composition is negligible and well within our margins of error.
Figure 6a shows the time evolution of the integrated CO, H2CO, and CH3OH signals in symbols. It shows how the amount of CO decreases as the abundance of H2CO increases for four different temperatures. After bombardment with
H atoms cm-2, the formation of methanol kicks off at the expense of the growth of the H2CO abundance. Similar trends of abundance evolution as a function of fluence are reported by Watanabe et al. (2006). This indicates that the fluence is determined with relatively high accuracy since in both experiments different atomic sources (Tschersich vs. microwave induced plasma) and different calibration methods are used.
![]() |
Figure 5: CO RAIR band position as a function of CH3OH content in a CO:CH3OH mixed ice obtained by codeposition experiments. |
Open with DEXTER |
![]() |
Figure 6:
Time evolution of the surface abundance (in molecules cm-2) of CO, H2CO and CH3OH during H-atom bombardment of CO ice with a H-atom flux of
|
Open with DEXTER |
3.2 Flux dependence
As mentioned in the introduction, the apparent discrepancy between the results by Hiraoka et al. (2002) and Watanabe & Kouchi (2002) has been attributed to a difference in the H-atom flux used in the respective experiments. The setup in our laboratory is able to cover the entire flux range from 1012 to 1014 cm-2 s-1. For high flux, both formaldehyde and methanol are formed as can be seen in Figs. 2 and 6 and in the corresponding work of Watanabe & Kouchi (2002).
A difference spectrum of a similar experiment but with a much lower flux of 1012 cm-2 s-1 is plotted in Fig. 7. The exposure time here is four hours to reach a higher total fluence of
cm-2, which is still significantly less than the sample experiment shown in Fig. 2. Note that the vertical scales in Figs. 2 and 7 cover the same range. For longer exposures, surface contamination would become a problem, but methanol features would eventually become detectable. As Fig. 7 clearly shows, a smaller fraction of CO is transformed into H2CO and the sensitivity of the RAIR spectrometer is not sufficiently high to confirm the formation of CH3OH in these circumstances. TPD, however, is more sensitive as a diagnostic tool, although harder to use for a quantitative or time resolved analysis. Figure 8 plots several TPD spectra. It shows a small methanol desorption peak around 150 K. We checked experimentally that the carrier of this peak is indeed formed in the ice during the hydrogen exposure and that the observed CH3OH is not a contaminant in the UHV chamber. This is a strong indication that the formation mechanism of formaldehyde and methanol does not fundamentally change with varying flux. The H2O desorption at 20-30 K originates in frozen background water on the surrounding parts of the cryohead.
The arrows in Fig. 6a indicate the corresponding fluences for the low and high flux experiments shown in Figs. 7 and 2, respectively. From this, it is immediately apparent that only a limited amount of methanol can be formed under low flux circumstances. Note that Hiraoka et al. (2002) probably used an even lower fluence since their exposure time was four times shorter than in our experiment. In addition, they used a slightly lower temperature of 10 K.
![]() |
Figure 7:
RAIR spectrum of a CO ice at 12.0 K exposed to a fluence of
|
Open with DEXTER |
![]() |
Figure 8: The TPD spectra corresponding to Fig. 7. |
Open with DEXTER |
3.3 Thickness dependence
The effect of the initial layer thickness on the formation yield of
H2CO and CH3OH is investigated by repeating the sample
experiment for different CO layer thicknesses. Figure 9 shows the absolute reaction yield after a fluence of
H atoms cm-2 as a function of the layer
thickness. In all cases, a steady state value for H2CO is reached for this
fluence. The figure clearly shows that for CO layers
thicker than
molecules cm-2, the absolute yield
is layer thickness independent and the results are reproducible within
the measurement error. The combined H2CO and CH3OH yield of
molecules cm-2 is lower than the
molecules cm-2 penetration column. From these
experiments, we conclude that the penetration column of the H atoms
into the CO ice is at most
molecules cm-2 at 12.0 K. This corresponds to 4 monolayers (ML) of solid
(bulk) CO molecules. At least half of the CO molecules in the active
layer is converted into H2CO and CH3OH. The
determination of the penetration column by this experiment is only an upper limit because of
the low thickness resolution in Fig. 9. It however
agrees with the previous estimate of nearly 100% conversion.
![]() |
Figure 9:
The absolute reaction yield of H2CO and CH3OH after a fluence of
|
Open with DEXTER |
3.4 Temperature dependence
Several experiments for different surface temperatures have been performed. The initial layer thickness and flux values are comparable to the values used in the sample experiment. Figures 6b-d show the results for hydrogenation experiments at 13.5, 15.0, and 16.5 K, respectively. These clearly indicate the very different evolution of CO, H2CO, and CH3OH abundance with temperature. Table 1 gives the initial formation rate of formaldehyde (slope at t=0) and the final H2CO and CH3OH yields, and also indicates whether or not a steady state is reached. The table shows that at early times, the formation rate of H2CO is much lower for higher temperatures than for 12.0 K. We will return to this point later. The final yield of CH3OH is however larger at 13.5 and 15.0 K. For T>15 K, the production rate of H2CO is simply so low that a steady state is not reached. Minimal amounts of formed methanol were also detected in experiments at 18.0 and 20.0 K, but since some CO desorption and redeposition occurs at these temperatures, they are not presented here for a quantitative discussion.
The appearance of the extra CO band at 2135 cm-1 indicates that for temperatures between 12.0 and 15.0 K a nearly pure methanol layer is formed. We expect a similar behaviour for formaldehyde. This means that the active CO layer involved in the reactions can be determined directly from the steady state yield of H2CO and CH3OH. This active layer increases with temperature indicating that the penetration column of H atoms into CO ice increases with temperature as one would expect. The CO molecules in the ice are more mobile at higher temperatures making it easier for H atoms to penetrate the CO ice, since the ice becomes less rigid. Note that the absolute temperature calibration in the set-up of Watanabe and ours appears to differ by 1-2 K (comparing Fig. 3 in Watanabe et al. 2006 and Fig. 6 here), but the observed trends are identical.
Table 1: The production rate and yield of H2CO and the yield of CH3OH.
4 Monte Carlo simulations
4.1 The method
To infer the underlying mechanisms leading to the formation of methanol, a detailed physical-chemical model is required. The present section discusses an approach based on the continuous time, random-walk Monte Carlo simulation technique. This method differs from previous studies based on rate equations and enables the study of surface processes to be performed in more detail. In addition, it provides a clearer understanding about what occurs physically on the surface. In contrast to an analysis using rate laws, the Monte Carlo method determines the H surface abundance by taking into account the layered structure of the ice, the H-atom flux, diffusion, reaction and desorption. This allows an extension of the results to conditions with much lower fluxes such as in the interstellar medium (ISM). For a detailed description of the method and program, we refer to Cuppen & Herbst (2007).
During a simulation, a sequence of processes - hopping, desorption,
deposition, and reaction - is performed, where this sequence is chosen
by means of a random number generator in combination with the rates
for the different processes. First, an initial ice layer is created by
deposition of CO on a surface. The resulting surface roughness of this
layer depends on temperature and flux. For the experimental
conditions simulated here, the CO ice is compact with a maximum
height difference across the surface of only 2-3 monolayers.
Hydrogen atoms and hydrogen molecules are subsequently deposited,
following their relative abundance in the H-atom flux, at an angle
perpendicular to the surface to mimic the experimental conditions.
They move, react, and desorb according to rates of a form similar to
that used in gas-grain models
![]() |
(4) |
where Ex is the activation energy for process X, and A is the pre-exponential factor for which a constant number of

This expression ensures microscopic reversibility between the different types of sites. The parameter


Diffusion into the ice is also considered. Minimum energy path
calculations suggest that CO and H can swap position enabling an H
atom to penetrate into the CO ice (see Appendix B). The barrier for
this process strongly depends on the layer in which the H atom is
situated. In the simulations, the barrier for this event is (
350 +
2(z1 + z2)) K for an H atom to swap between layer z1 and
z2. This compares to a hopping barrier of
K and a desorption energy of
K (see next section).
Hiraoka et al. (1998) found that hydrogen atoms can relatively easily
diffuse through the CO ice. Moreover, the current experiments show
that hydrogen atoms can penetrate into a maximum of four monolayers
for 12.0 K. Hydrogen atoms are also allowed to swap with
formaldehyde and methanol, but here the initial barrier is chosen to
be higher (450 and 500 K) since these species are heavier and are
more strongly bound within the ice matrix.
4.2 The CO ice layer
Although the experimental CO layers are probably amorphous
(Kouchi 1990), crystalline layers are used in the Monte Carlo
simulations discussed here. In this way, a lattice-gas Monte Carlo
method can be used, which enables far longer simulation times than
in off-lattice methods. We expect the crystalline assumption to be
reasonable since the local structure of the CO layers is probably
close to crystallinity. The energy released during deposition may
help the molecules to rearrange slightly, leading to
micro-crystalline domains. The -CO structure
(Vegard 1930) is used with layers in the (110) orientation.
The dominant faces of a CO crystal will have this crystallographic
orientation. The CO surface consists of alternating carbon and
oxygen terminated bi-layers. In the bulk configuration, each CO molecule has 14 nearest neighbours:
five in layers below, five in layers above, and four in the same
layer. The additive energy contribution of these neighbours is 2Efor the layers below and E for the neighbours in the same layer or
of lower z, where z is the depth in layers with respect to the
top layer. The different treatment for sites below the particle is
to add a contribution for longer range interactions from the ice
layer. E is chosen to be 32 K for atomic hydrogen, and 63 K for
CO. This leads to a binding energy of
K for H on top of a flat CO ice layer, and
K and
K for H
and CO, respectively, embedded in a CO layer. These values agree
very well with binding energies obtained by calculations using
accurate H-CO and CO-CO potentials of 320, 440, and 850 K,
respectively (see Appendix B).
4.3 Comparison to the experiment
The solid lines in Fig. 6 represent the results from the Monte Carlo calculations. The exact mechanisms included in these simulations are discussed in more detail in the following sections. The resulting time evolution series are in very good agreement for 12.0 K. The agreement for 13.5, 15.0, and 16.5 K is far less good, probably because of missing mechanisms that promote the penetration into the ice. In the current simulations, only swapping of species is included. Because of the thermal motion of the CO molecules, ``real'' penetration in which the H atoms penetrate in the CO matrix may also be possible. The shape of the curves is reproduced and only the H2CO abundance levels off at too low yields.
The main parameters varied to fit the experimental data are the reaction barriers and the diffusion rates. The best-fit model barriers are summarised in Table 2. Since the intermediate species HCO and H3CO are not experimentally detected, the barriers for hydrogenation of these species are significantly lower than for the other two reactions, presumably even zero. The HCO and H3CO abundances stay below detectable levels in the simulations. The reaction barriers for H + CO and H + H2CO are temperature dependent and increase with temperature. Our values are in good absolute agreement with the barriers found by Awad et al. (2005), who also found a similar temperature behaviour. Their values were obtained using a rate equation analysis for T=10, 15, and 20 K using the data from Watanabe et al. (2006). The temperature dependence suggests that there is a clear tunnelling component for the reaction at low temperature. The two barriers for forming H2CO and CH3OH show different temperature dependencies. The formation of methanol becomes relatively more important at higher temperature. Note that the Monte Carlo method automatically treats a reaction in competition with desorption and hopping. This is in contrast to gas-grain codes, where it must be included explicitly. To describe the chemical processes properly, one should introduce this competition into the gas-grain model.
Table 2: Reaction rates and barriers for CO + H and H2CO + H for different temperatures.
The errors in the energy barriers reflect the errors due to the uncertainties in the sticking probability, H-atom flux, diffusion, and exact structure of the CO ice.Molecular hydrogen is formed on the surface with efficiencies ranging from 3% (T =16.5 K) to 70% (T =12.0 K). However, because of the large excess energy of the formation reaction, the majority of the formed H2 molecules leaves the surface, and the H2 surface abundance is predominantly determined by impinging H2 molecules.
4.4 Effect of diffusion
Since the diffusion rates are uncertain, this section discusses the effect of diffusion in more detail. Minimum energy path calculations of the diffusion of a single hydrogen atom on a CO (110) surface (see Appendix B) results in energy barriers ranging from 70 to 170 K (




![]() |
Figure 10:
Monte Carlo simulations of the time evolution of the surface abundance of CO, H2CO and CH3OH during H-atom bombardment of CO ice at 12.0 K ( top) and 15.0 K ( bottom). Reaction barriers for H + CO and H + H2CO can be found in Table 2. The diffusion is varied via the parameter |
Open with DEXTER |
4.5 Effect of H2 molecules on the hydrogenation
All simulations include the deposition of both H atoms and H2 molecules, which results from the undissociated H2 molecules in the H-beam. If the H2 molecules are excluded from the simulations, the formation of H2CO and CH3OH is affected in only a limited number of cases of both fast diffusion and high temperature. The presence of H2 appears to have two principle effects: it limits the penetration into the ice and decelerates the H atoms, since they move through a ``sea'' of H2. The first effect reduces the production rate, whereas the second depends on the reaction barrier.
The experimental results at temperatures higher than 12.0 K show non-first-order behaviour at early times (exponential decay of CO). The H2CO production rate increases until 30 and 50 min of exposure for T = 13.5 and 15.0 K, respectively. After this time, the H2CO and CH3OH follow the expected first order behaviour. None of the simulations in Fig. 10 show this trend. The only mechanism capable of describing this phenomenon is an increasing effective H-atom flux with time. This increase in effective flux can be cased by an increase in the sticking of atomic hydrogen to the surface. Since the incoming H atoms are relatively warm, they need to dissipate this extra energy into the surface to stick. Because CO is relatively heavy compared to the H atoms, this energy dissipation will be inefficient, and most of the H atoms will scatter back into the gas phase. Once the surface abundance of the much lighter H2 molecules increases, the sticking of the H atoms to the surface will increase as well. We assume a 1% sticking for H atoms and H2 molecules on a bare CO surface and a 65% sticking of H atoms on a surface that is fully covered with H2. The sticking probability is further assumed to grow linearly with the H2 coverage. The H2 surface abundance reaches a constant level of 0.39 ML after a few minutes for T is 12.0 K. This results in a sticking of H atoms to the surface of 26%. For higher temperatures, it is noticeably longer time before to a steady state is reached, explaining the non-linear behaviour at early times and inferring a lower final sticking probability. The solid lines in Fig. 6 include this mechanism.
As mentioned earlier, Watanabe & Kouchi (2002) concluded that the temperature of the beam has little effect on the hydrogenation process, which seems to contradict our H2 argument. However, their experiments were carried out at 10 K, where the surfaces are covered with hydrogen atoms early on in the experiment because of the enhanced sticking at low temperatures. They further reported an unknown flux difference between the cold and warm beam, which makes quantifying the sticking probability using these experiments not possible. In conclusion, the temperature of the beam can affect the effective flux of H atoms landing on the surfaces, but it does not introduce additional energetic effects that influence the crossing of the barrier.
5 CO hydrogenation under interstellar conditions
Based on the fitting results in the previous section, the Monte Carlo routine can now be used to simulate CO hydrogenation reactions under interstellar conditions. An important ingredient is the H-atom density in the cloud. As in our laboratory beam, the gas in dense clouds consists of a mix of H and H2. Under steady-state conditions, the balance of the rates of H2 formation on grains and H2 destruction by cosmic rays infers an H-atom density of around 1 cm-3 (Hollenbach & Salpeter 1971). This H-atom number density is independent of the total density because both the formation and destruction rates scale with density. Before steady-state is reached, however, the H-atom density may be higher because the timescale for H to H2 conversion is long (


The starting configuration for the simulations is a layer of pure CO ice. This is believed to be representative for the top layers of the grain mantles at the centre of a high-density collapsing cloud. Here, the ice layer is observed to consist of predominantly CO ice as the result of ``catastrophic'' CO freeze-out (Pontoppidan 2006; Pontoppidan et al. 2008). More heterogeneous ice layers are formed at lower densities where CO and H2O are mixed, or towards the centre of proto-stellar envelops or proto-planetary disks where the dust has been heated and CO has desorbed from the top layers.
![]() |
Figure 11: Monte Carlo simulations of CO-ice hydrogenation at 12.0 ( top) and 16.5 K ( bottom). A constant atomic hydrogen gas phase density of 10 cm-3 and a gas temperature of 20 K is assumed. Thick lines represent interstellar conditions, thin lines are the scaled experimental simulations. The results are shown as the change in column density compared with t=0 yr. |
Open with DEXTER |
The bottom panel in Fig. 11 shows similar trends for 16.5 K.
Again the onset of H2CO and CH3OH formation is at much lower
fluences than in the experiment. At the end of the simulation,
nearly all H2CO has been converted into CH3OH. This is in
contrast to the 12.0 K simulations, where a constant non-zero amount
of H2CO remains after
yr. The crossover point from
H2CO-rich to CH3OH-rich ice occurs at slightly later times at 16.5 K compared to 12.0 K. This can clearly be seen in
Fig. 12, which plots the H2CO/CH3OH ratio for
both temperatures. At early times, this ratio is similar for 12.0
and 16.5 K. At t > 103 yr, the ratio starts to level off at 12.0 K, while it still decreases rapidly at 16.5 K. The noise in the
curve for 16.5 K below
yr is due to the low
abundances of H2CO and CH3OH.
In space, the H2CO/CH3OH ice ratio has been determined
directly for only three high-mass young stellar objects (YSOs): W 33A, NGC 7538 IRS9, and AFGL 70009S, with inferred ratios ranging
from 0.09 to 0.51 (Keane et al. 2001; Gibb et al. 2004). The laboratory
curves for the H2CO and CH3OH production show that H2CO is
more or equally abundant during most of our experiments. Thus,
values as low as 0.09-0.51 cannot easily be reproduced in the
experiments. However, the Monte Carlo simulations for interstellar
conditions have a crossover from H2CO-rich to CH3OH-rich ice
at significantly earlier times than the experimental curves and a
H2CO/CH3OH ratio of 0.51 is obtained after
yr
at
K. Grains at higher temperatures will have this
crossover at even earlier times and for grains with
K, a H2CO/CH3OH ratio of even 0.09 is obtained
after
yr. Thus, the observed ratios are in agreement
with the models discussed above for chemical timescales >
yr, which is consistent with the estimated ages of these
high-mass protostars of a few
104-105 yr (Hoare et al. 2007).
CH3OH ice has also been detected toward low-mass YSOs with
abundances ranging from <1% to more than 25% of the H2O ice
abundance (Pontoppidan et al. 2003; Boogert et al. 2008). An interesting
example is the Class 0 protostar Serpens SMM 4, for which
a particularly high CH3OH abundance of 28% with respect to
H2O ice was deduced for the outer envelope
(Pontoppidan et al. 2004). The upper limit to the H2CO-ice
abundance implies a H2CO/CH3OH ratio <0.18, implying an
age >
yr at 16.5 K. This is consistent with the
estimated timescale for heavy freeze out in low-mass YSOs of
yr, including both the pre-stellar and proto-stellar
phases (Jørgensen et al. 2005).
![]() |
Figure 12: The H2CO/CH3OH ratio as a function of time obtained from the Monte Carlo simulations of CO hydrogenation at 12.0 and 16.5 K under ISM conditions (see Fig. 11). The grey box indicates Spitzer ice observations, the black box gas phase observations. |
Open with DEXTER |

In contrast, since the CH3OH ice abundance with respect to that
of H2O is known to vary by more than an order of magnitude, both
local conditions and timescales appear to play a role. Note,
however, that for CH3OH abundances as large as 25% (columns as
large as 1018 cm-2), the CH3OH layer is approximately 25 ML thick (
ML), much more than can be produced
from just the upper 4 ML of the CO ice. Thus, conversion of CO into
CH3OH ice must in these cases occur simultaneously with the
freeze-out and building up of the CO layer. Pure CO ice can also
easily desorb as soon as the protostar heats up. This complicates
the use of CH3OH/CO ice as an evolutionary probe. A proper model
of interstellar CH3OH ice formation should therefore include the
effects of the variations in CO-ice abundances and dust temperatures
in the pre- and protostellar phases, and account for the differences
in timescales for CH3OH-ice formation compared with those of CO
adsorption and desorption. This paper provides the necessary
molecular data to compile such a model.
6 Conclusion
The present paper shows that the formation of methanol by successive hydrogenation of CO and H2CO is efficient under various laboratory conditions covering




Energy barriers for the H + CO and H2CO + H reactions are obtained by fitting Monte Carlo simulation results to the experimental data. Using these barriers, the methanol production is simulated for interstellar conditions. The obtained H2CO and CH3OH abundances do not scale directly with fluence because of the different relative importance of H2 production and CO hydrogenation in space compared with the laboratory, as can be clearly seen by comparing the thick and thin lines in Fig. 11. However, laboratory experiments are required to derive the necessary rates that serve as input to the Monte Carlo program. The obtained H2CO/CH3OH ratios for the interstellar simulations are in closer agreement with observational limits than a direct translation of the experimental observations.
Monte Carlo simulations of the hydrogenation process show that the presence of H2 has three effects: it promotes the sticking of the warm H atoms, it limits the penetration into the ice, and it slows down the diffusion of H atoms. The first effect will be negligible under interstellar conditions since the incoming H atoms will be cold already and the sticking probability will therefore be high regardless of the substrate. The latter two effects will be important and are similar to the conditions in the laboratory for also a high H2 abundance.
The experiments show that the hydrogenation process is thickness
independent for layers thicker than
cm-2,
and that the active layer, which contains only a limited amount of
CO after a steady state is reached, becomes slightly thicker with
temperature. For temperatures higher than 15.0 K, a clear drop in
the production rate of methanol is observed. This is probably due to
two effects: the desorption of H atoms becoming important and the
sticking of H atoms being reduced because of the low H2 surface
abundance. Both effects cause the H surface abundance to drop
substantially at those temperatures and therefore reduce the
probability of hydrogenation reactions occurring
in the laboratory. Simulations of CO hydrogenation in space show a strong temperature dependence of the H2CO/CH3OH ratio over several orders of magnitude. The CH3OH abundance changes with time, temperature, and fluence.
Acknowledgements
Part of this work was supported by the Netherlands Research School for Astronomy, NOVA, and Netherlands Organization for Scientific Research (NWO) through a VENI grant. We thank Stephan Schlemmer and Helen Fraser for their contribution during the first construction phase and Gijsbert Verdoes, Martijn Witlox and Ewie de Kuyper from the Fijn Mechanische Dienst for their support. Ayman Al-Halabi, Lou Allamandola, Eric Herbst, Xander Tielens, Klaus Pontoppidan, and Zainab Awad have contributed to this work through long and inspiring discussions.
Appendix A: Absolute and relative H-atom flux determination
A.1 Absolute flux determination
The (accuracy of the) absolute value of the H-atom flux at the ice surface is obtained by estimating lower and upper limits in two independent ways. We exemplify here the H-atom flux determination for the case of our standard values with an H2 pressure in the chamber of
The lower limit to the absolute flux is directly available from the
experimental results presented in Ioppolo et al. (2008). That paper
discusses the H2O2 and H2O production from H-atom
bombardment of O2-ice in time using the same setup and settings.
During the first hour, H2O2 and H2O are produced with an
almost constant production rate of
molecules cm-2 s-1. Since both molecules contain two hydrogen atoms,
this means that the H-atom flux should be at least twice this value.
Assuming a conservative sticking probability of hydrogen atoms at 300 K to O2-ice between 12 K and 28 K of at most 50%, we
determine a lower limit to the flux of
cm-2 s-1.
The determination of the upper limit to the H-atom flux is more elaborate and involves several steps. Figure 1 shows that the hydrogen atoms travel from the source through the atomic-line chamber to a quartz pipe, where the atoms are collisionally cooled and then travel through the main chamber onto the substrate. The final H-atom flux is then determined by
![]() |
(A.1) |
where

Our specific hydrogen source, used
in the experiments described here, has been tested prior to delivery at the Forschungszentrum in Jülich where the flux, solid angle, and dissociation rate have been measured for a wide range of H2 pressures and filament temperatures. The set-up used for these calibration experiments is described in Tschersich & von Bonin (1998). These measurements confirmed that there is little variation between individual instruments, since nearly identical rates were obtained by Tschersich & von Bonin (1998) and Tschersich (2000) and later by Tschersich et al. (2008) for different H-atom sources of the same type. From the flux and dissociation rate measured in Jülich,
can be obtained as well as k1 using the solid angle information. In our example case,
H-atoms s-1 leave the H-atom source and 44% of these atoms enter the quartz pipe, which is located at a distance of 1.5 cm.
The pipe is designed so that the atoms cannot reach the substrate directly and that the number of hydrogen recombinations is kept to a minimum. This is achieved by using a short pipe with a high diameter/length ratio and choosing quartz, which is known to have a low recombination efficiency. Following Walraven & Silvera (1982), a theoretical estimate of the number of recombinations in the pipe can be determined, considering the specific shape and material. This reduces the H-atom flux by another 27%. The pipe ends in close proximity to the cryogenic surface. The use of a pipe instead of a pinhole or a slit results in a focused H-atom beam for which the flux can be determined to relatively low uncertainty. From geometric considerations, a minimum solid angle can be estimated. This will suffice, since our aim is to obtain an upper limit to the flux. The H-atom beam covers A=4.9 cm2 of the substrate that is located 3 cm behind the quartz pipe. This spot falls completely on the surface and k2 can readily be assumed to be unity.
Finally, the pressure drop between the source and the main chamber can be determined in two ways: by a calculation using the conductance of the pipe and the pumping speed, and by measuring the pressures in both chambers using undissociated beams. Both results are in reasonable agreement, leading to
.
Our upper limit to the flux is now
![]() |
= | ![]() |
|
= | ![]() |
(A.2) |
Deviations from this upper limit are expected to be due to a lower k1 value, because of misalignments between the source and the entrance of the quartz pipe, an underestimation of the solid angle of the exiting beam from the quartz pipe (lower k1 and higher A), more recombinations in the pipe or backscattering of atoms from the quartz pipe, to the chamber of the H-atom source.
The value of flux adopted in the present paper is the resultant intermediate H-atom flux of
cm-2 s-1, which is within a factor of 2 of the upper and lower limits. It should be noted that this is a conservative error, since the true lower and higher flux limits are likely to be higher and lower, respectively.
A.2 Relative flux determination
The accuracy in the relative flux is particularly important to the conclusion presented in this paper, more so than the absolute value. For this we use the CO-hydrogenation data obtained from the experiments. Figure A.1 shows the CO, H2CO and CH3OH evolution as a function of fluence for three different fluxes. The fluences are calculated using the flux determination as described above. The three curves clearly overlap, which means that the accuracy of the relative fluxes is well within our error bars. We conclude that the accuracy in the relative flux is substantially higher than the accuracy of the absolute flux, well below 50%. One of the main conclusions of the paper, that the discrepancy between the two Japanese groups is due to a difference in flux, as envisaged by Hidaka et al. (2004), is therefore robust.Finally, reproducing the same experiments on different days over the course of several months showed that reproducibility over periods from day-to-day to months is excellent, to within a few percent.
![]() |
Figure A.1:
Time evolution of the surface abundance (in molecules cm-2) of CO, H2CO and CH3OH during H-atom bombardment of CO ice at 12.0 K with three different H fluxes of
|
Open with DEXTER |
Appendix B: Binding energy calculations
To calculate binding energies and barriers to diffusion, recently-developed CO-CO and H-CO potentials are used. Takahashi and van Hemert (in prep.) have fitted high level electronic structure (coupled cluster) calculations of the CO-CO dimer to an analytic potential consisting of partial charges on the atoms and the centres of mass of the CO molecules, atom-based Lennard-Jones type interactions, and Morse potentials for the intramolecular C-O interaction. In the work by Andersson et al. (in prep.), a potential for the interaction between a hydrogen atom and CO has been calculated by fitting damped dispersion and exponential repulsion potentials to coupled cluster calculations.
Using the CO-CO potential, a CO (110) surface has been created consisting of
528 CO molecules in 11 monolayers in a cell with dimensions
Å in the surface plane. By applying periodic boundary
conditions, an infinite surface
is created. Binding energies have been calculated by performing energy
minimisations for H atoms at different sites on top
of and inside the CO surface and by comparing to the energy when the hydrogen
is far away from the surface. In the same manner,
the binding energy for a CO molecule in the top layer has been calculated.
In all instances, the top 3 monolayers of the ice have been allowed to relax.
To calculate the energy barriers to diffusion both on and into the surface, the nudged elastic band (NEB) method (Jónsson et al. 1998) is used initially to map out the minimum energy path (MEP) connecting two potential minima. To fine-tune the barrier height, the Lanczos method is used to optimise the saddle point of the potential energy (Olsen et al. 2004).
References
- Awad, Z., Chigai, T., Kimura, Y., Shalabiea, O. M., & Yamamoto, T. 2005, ApJ, 626, 262 [NASA ADS] [CrossRef] (In the text)
- Bisschop, S. E. 2007, Ph.D. Thesis, Leiden Observatory, Leiden University
- Bisschop, S. E., Fuchs, G. W., van Dishoeck, E. F., & Linnartz, H. 2007a, A&A, 474, 1061 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Bisschop, S. E., Jørgensen, J. K., van Dishoeck, E. F., & de Wachter, E. B. M. 2007b, A&A, 465, 913 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Boogert, A. C. A., Pontoppidan, K. M., Knez, C., et al. 2008, ApJ, 678, 985 [NASA ADS] [CrossRef]
- Bouwman, J., Ludwig, W., Awad, Z., et al. 2007, A&A, 476, 995 [NASA ADS] [CrossRef] [EDP Sciences]
- Charnley, S. B., Tielens, A. G. G. M., & Millar, T. J. 1992, ApJ, 399, L71 [NASA ADS] [CrossRef] (In the text)
- Cuppen, H. M., & Herbst, E. 2007, ApJ, 668, 294 [NASA ADS] [CrossRef] (In the text)
- Garrod, R., Park, I. H., Caselli, P., & Herbst, E. 2006, Disc. Faraday Soc., 133, 5
- Geppert, W. D., Thomas, R. D., Ehlerding, A., et al. 2005, J. Phys.: Conf. Ser., 4, 26 [NASA ADS] [CrossRef]
- Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A., & Tielens, A. G. G. M. 2004, ApJS, 151, 35 [NASA ADS] [CrossRef]
- Goldsmith, P. F., Li, D., & Krco, M. 2007, ApJ, 654, 273 [NASA ADS] [CrossRef] (In the text)
- Hidaka, H., Watanabe, N., Shiraki, T., Nagaoka, A., & Kouchi, A. 2004, ApJ, 614, 1124 [NASA ADS] [CrossRef] (In the text)
- Hiraoka, K., Miyagoshi, T., Takayama, T., Yamamoto, K., & Kihara, Y. 1998, ApJ, 498, 710 [NASA ADS] [CrossRef] (In the text)
- Hiraoka, K., Sato, T., Sato, S., et al. 2002, ApJ, 577, 265 [NASA ADS] [CrossRef]
- Hoare, M. G., Kurtz, S. E., Lizano, S., Keto, E., & Hofner, P. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 181 (In the text)
- Hollenbach, D., & Salpeter, E. E. 1971, ApJ, 163, 155 [NASA ADS] [CrossRef] (In the text)
- Ioppolo, S., Cuppen, H. M., Romanzin, C., van Dishoeck, E. F., & Linnartz, H. 2008, ApJ, 686, 1474 [NASA ADS] [CrossRef]
- Jónsson, H., Mills, G., & Jacobsen, K. W. 1998, Nudged Elastic Band Method for Finding Minimum Energy Paths of Transitions, ed. B. J. Berne, G. Ciccotti, & D. F. Coker (Singapore: World Scientific), 385 (In the text)
- Jørgensen, J. K., Schöier, F. L., & van Dishoeck, E. F. 2005, A&A, 435, 177 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Katz, N., Furman, I., Biham, O., Pirronello, V., & Vidali, G. 1999, ApJ, 522, 305 [NASA ADS] [CrossRef] (In the text)
- Keane, J. V., Tielens, A. G. G. M., Boogert, A. C. A., Schutte, W. A., & Whittet, D. C. B. 2001, A&A, 376, 254 [NASA ADS] [CrossRef] [EDP Sciences]
- Kouchi, A. 1990, J. Cryst. Growth, 99, 1220 [NASA ADS] [CrossRef] (In the text)
- Miyauchi, N., Hidaka, H., Chigai, T., et al. 2008, Chem. Phys. Lett., 456, 27 [NASA ADS] [CrossRef]
- Nagaoka, A., Watanabe, N., & Kouchi, A. 2005, ApJS, 624, L29 [NASA ADS] [CrossRef]
- Olsen, R. A., Kroes, G. J., Henkelman, G., Arnaldsson, A., & Jónsson, H. 2004, J. Chem. Phys., 121, 9776 [NASA ADS] [CrossRef] (In the text)
- Palumbo, M. E., & Strazzulla, G. 1993, A&A, 269, 568 [NASA ADS]
- Pontoppidan, K. M. 2006, A&A, 453, L47 [NASA ADS] [CrossRef] [EDP Sciences]
- Pontoppidan, K. M., Dartois, E., van Dishoeck, E. F., Thi, W.-F., & d'Hendecourt, L. 2003, A&A, 404, L17 [NASA ADS] [CrossRef] [EDP Sciences]
- Pontoppidan, K. M., van Dishoeck, E. F., & Dartois, E. 2004, A&A, 426, 925 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Pontoppidan, K. M., Boogert, A. C. A., Fraser, H. J., et al. 2008, ApJ, 678, 1005 [NASA ADS] [CrossRef]
- Sandford, S. A., Allamandola, L. J., Tielens, A. G. G. M., & Valero, G. J. 1988, ApJ, 329, 498 [NASA ADS] [CrossRef]
- Teolis, B. D., Loeffler, M. J., Raut, U., Famá, M., & Baragiola, R. A. 2007, Icarus, 190, 274 [NASA ADS] [CrossRef] (In the text)
- Tielens, A. G. G. M., & Charnley, S. B. 1997, Orig. Life Evol. Biosph., 27, 23 [NASA ADS] [CrossRef]
- Tielens, A. G. G. M., & Hagen, W. 1982, A&A, 114, 245 [NASA ADS]
- Tschersich, K. G. 2000, J. Appl. Phys., 87, 2565 [NASA ADS] [CrossRef]
- Tschersich, K. G., & von Bonin, V. 1998, J. Appl. Phys., 84, 4065 [NASA ADS] [CrossRef]
- Tschersich, K. G., Fleischhauer, J. P., & Schuler, H. 2008, J. Appl. Phys., 104, 034908 [NASA ADS] [CrossRef] (In the text)
- Vegard, I. 1930, Z. Physik, 61, 185 [NASA ADS] [CrossRef] (In the text)
- Walraven, J. T. M., & Silvera, I. F. 1982, Rev. Sci. Instr., 53, 1167 [NASA ADS] [CrossRef] (In the text)
- Watanabe, N., & Kouchi, A. 2002, ApJ, 571, L173 [NASA ADS] [CrossRef]
- Watanabe, N., Shiraki, T., & Kouchi, A. 2003, ApJ, 588, L121 [NASA ADS] [CrossRef]
- Watanabe, N., Nagaoka, A., Shiraki, T., & Kouchi, A. 2004, ApJ, 616, 638 [NASA ADS] [CrossRef]
- Watanabe, N., Hidaka, H., & Kouchi, A. 2006, in Astrochemistry - From Laboratory Studies to Astronomical Observations, ed. R. I. Kaiser, P. Bernath, Y. Osamura, S. Petrie, & A. M. Mebel, AIP Conf. Ser., 855, 122
All Tables
Table 1: The production rate and yield of H2CO and the yield of CH3OH.
Table 2: Reaction rates and barriers for CO + H and H2CO + H for different temperatures.
All Figures
![]() |
Figure 1: Schematic representation of the experimental setup. CO ice is deposited through the dosing line and the products are monitored by means of infrared spectroscopy and quadrupole mass spectroscopy (QMS). |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
RAIR difference spectrum of a CO ice at 12.0 K exposed to
|
Open with DEXTER | |
In the text |
![]() |
Figure 3: The decrease in integrated absorbance of CO and CH3OH (1035 cm-1) following desorption at a constant temperature of 29 and 135 K, respectively. The arrows indicate the transition points from the multi- to sub-monolayer regime. |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Spectral change of the CO 2143 cm-1 RAIR band before and after H-atom bombardment. The inset shows the corresponding difference spectrum. |
Open with DEXTER | |
In the text |
![]() |
Figure 5: CO RAIR band position as a function of CH3OH content in a CO:CH3OH mixed ice obtained by codeposition experiments. |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Time evolution of the surface abundance (in molecules cm-2) of CO, H2CO and CH3OH during H-atom bombardment of CO ice with a H-atom flux of
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
RAIR spectrum of a CO ice at 12.0 K exposed to a fluence of
|
Open with DEXTER | |
In the text |
![]() |
Figure 8: The TPD spectra corresponding to Fig. 7. |
Open with DEXTER | |
In the text |
![]() |
Figure 9:
The absolute reaction yield of H2CO and CH3OH after a fluence of
|
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Monte Carlo simulations of the time evolution of the surface abundance of CO, H2CO and CH3OH during H-atom bombardment of CO ice at 12.0 K ( top) and 15.0 K ( bottom). Reaction barriers for H + CO and H + H2CO can be found in Table 2. The diffusion is varied via the parameter |
Open with DEXTER | |
In the text |
![]() |
Figure 11: Monte Carlo simulations of CO-ice hydrogenation at 12.0 ( top) and 16.5 K ( bottom). A constant atomic hydrogen gas phase density of 10 cm-3 and a gas temperature of 20 K is assumed. Thick lines represent interstellar conditions, thin lines are the scaled experimental simulations. The results are shown as the change in column density compared with t=0 yr. |
Open with DEXTER | |
In the text |
![]() |
Figure 12: The H2CO/CH3OH ratio as a function of time obtained from the Monte Carlo simulations of CO hydrogenation at 12.0 and 16.5 K under ISM conditions (see Fig. 11). The grey box indicates Spitzer ice observations, the black box gas phase observations. |
Open with DEXTER | |
In the text |
![]() |
Figure A.1:
Time evolution of the surface abundance (in molecules cm-2) of CO, H2CO and CH3OH during H-atom bombardment of CO ice at 12.0 K with three different H fluxes of
|
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.