Single and double addition of oxygen atoms to propyne on surfaces at low temperatures

Experiments designed to simulate the low temperature surface chemistry occurring in interstellar clouds provide clear evidence of a reaction between oxygen atoms and propyne ice. The reactants are dosed onto a surface held at a ﬁ xed temperature between 14 and 100 K. After the dosing period, temperature programmed desorption (TPD), coupled with time-of-ﬂ ight mass spectrometry, are used to identify two reaction products with molecular formulae C 3 H 4 O and C 3 H 4 O 2 . These products result from the addition of a single oxygen atom, or two oxygen atoms, to a propyne reactant. A simple model has been used to extract kinetic data from the measured yield of the single-addition (C 3 H 4 O) product at surface temperatures from 30 – 100 K. This modelling reveals that the barrier of the solid-state reaction between propyne and a single oxygen atom (160 (cid:1) 10 K) is an order of magnitude less than that reported for the gas-phase reaction. In addition, estimates for the desorption energy of propyne and reaction rate coe ﬃ cient, as a function of temperature, are determined for the single addition process from the modelling. The yield of the single addition product falls as the surface temperature decreases from 50 K to 30K, but rises again as the surface temperature falls below 30 K. This increase in the rate of reaction at low surface temperatures is indicative of an alternative, perhaps barrierless, pathway to the single addition product which is only important at low surface temperatures. The kinetic model has been further developed to characterize the double addition reaction, which appears to involve the addition of a second oxygen atom to C 3 H 4 O. This modelling indicates that this second addition is a barrierless process. The kinetic parameters we extract from our experiments indicate that the reaction between atomic oxygen and propyne could occur under on interstellar dust grains on an astrophysical time scale.


Introduction
The elemental composition of the known universe comprises almost exclusively light atoms ($99.9%hydrogen and helium).However, to date, over 160 different molecules have been detected in the interstellar medium. 1 The vast majority of these interstellar molecules contain hydrogen.The distribution of these molecules is far from uniform across the interstellar medium (ISM), where dense interstellar clouds are observed to harbour relatively large densities and varieties of complex molecules that have been proposed as precursors to biologically relevant species. 2,3In these interstellar clouds, molecular lifetimes are extended with respect to other more diffuse parts of the ISM.Specically, the relative opacity of these dense interstellar clouds shields the cloud's centre from high energy photons, allowing relatively fragile molecules to survive for extended periods.
In dense interstellar clouds, molecules are detected in abundances that cannot be completely accounted for by known gas-phase kinetics. 4Hence, it is now widely accepted that there is a contribution to these molecular abundances from reactions on and within the molecular ices on the surfaces of interstellar dust.This interstellar dust comprises predominantly silicate or carbonaceous particles, with a typical diameter of 100 nm, 5 making up typically 1% of the mass of an interstellar cloud. 6n general, bimolecular reactions on surfaces are thought to follow one of two prototypical reaction pathways.The rst is the Langmuir-Hinshelwood (LH) mechanism where both reactive species are initially adsorbed and thermalised on the surface. 7Diffusive processes then allow the reactants to encounter one another and react.The second mechanism, the Eley-Rideal (ER) pathway, involves a thermalised surface molecule undergoing direct reaction with an incident, and potentially energetic, gas-phase species. 7To correctly model the contribution of these heterogeneous reactions to the abundances of interstellar molecules, kinetic data is required.The measurement of such kinetic data is the objective of our experimental efforts. 8,9ne relatively abundant interstellar molecule, that is considered to play an important role in astrochemical processes, is the small unsaturated hydrocarbon, propyne (CH 3 CCH). 10,11Interstellar propyne was rst identied in the Milky Way in the giant molecular cloud Sagittarius B2. 11,12 Propyne has also been observed in cold cloud cores and lukewarm corinos in the Milky Way. 13 Propyne has also been detected across the disc of the distant Messier 82 (M82) galaxy. 14In M82, the propyne-to-methanol ratio ([CH 3 CCH]/[CH 3 OH] > 8) is higher than for comparable starburst galaxies such as NGC 253 ([CH 3 CCH]/[CH 3 OH] z 1). 15The propyne component of a prototypical starburst galaxy such M82 and NGC 253 resides primarily in interstellar clouds, where thermal dust grain chemistry is potentially important.Here we distinguish thermal surface chemistry, where reactions occur on the surface without the input of additional energy, with activated chemistry, where reactions are initiated and molecules subsequently processed, by external agents such as cosmic rays. 16o our knowledge, despite propyne being repeatedly detected in the ISM, the reactivity of propyne on cold surfaces, such as those of dust grains, has not yet been investigated.8][19][20] Given their comparable abundances, the surface reactivity of propyne appears to be overdue an investigation.Such an investigation is the target of this study.
In addition to its detection in interstellar clouds, propyne has also been observed in planetary environments in our Solar System.Specically, propyne has been detected in stratospheres of Jupiter and Saturn, displaying column densities of 1.5 AE 4Â10 14 molecules cm À2 and 2Â10 15 molecules cm À2 respectively. 21,22nterstellar ice composition is considered to be important when discussing the abundance of propyne in the Saturnian and Jovian atmospheres.4][25][26][27][28] Propyne has also recently been identied in the atmosphere of Uranus 29 and the atmosphere of Saturn's moon Titan, where it appears enriched at northern latitudes. 30,31In all cases propyne is identied by its infrared absorption band (n 9 ) centred at 15.8 mm.
Oxygen is the third most abundant element in the ISM, aer hydrogen and helium. 324][35] It has been proposed that interstellar grains could act as a sink for oxygen atoms, but it appears that the necessary depletion cannot be generated simply by the incorporation of oxygen atoms into the structures of the dust grains themselves. 34One possible additional sink for oxygen atoms is perhaps their reaction with the organic component of the molecular ice mantles that coat the dust grains in some interstellar clouds.The reaction we investigate in this paper, the addition of oxygen atoms to propyne on a cold surface, could be considered representative of oxygen atom depletion by reaction with organic ices.
Oxygen atom diffusion has recently been studied on an amorphous water ice surface. 36In this work Minissale et al 36 report that oxygen atom diffusion exhibits a quantum-classical transition at a surface temperature of 20 K.That is, above 20 K a classical motion involving barrier 'hopping' between adsorption sites is the dominant O atom migration mechanism, but at surface temperatures below 20 K quantum tunnelling of the O atoms becomes dominant.
8][39][40][41][42][43][44] In the gas phase the reaction produces CO and the CH 3 CH diradical which form from a vibrationally excited methyl ketene. 45At high pressures, or in a condensed medium where intermediates can be collisionally stabilized, the adduct that results from O atom addition is expected to simply relax to a more stable isomer such as methyl ketene or propenal. 38,46n this paper we present a study of the reaction of propyne with atomic oxygen on a cold surface.These experiments generate kinetic data for this reaction under astrophysically relevant conditions.These studies build on our previous investigations of the reaction of oxygen atoms with alkenes 9 and sulphur containing molecules. 8

Experimental procedure
The experimental apparatus employed in this study, which has been described in detail before, 8,9 is designed to probe the reactions of molecules with atomic species on molecular ices deposited on a graphite surface held at a specic temperature.The surface temperatures investigated are pertinent to heterogeneous reactivity on dust surfaces in the ISM.The products of any surface reactions are detected by a temperature-programmed desorption (TPD) methodology.In the TPD experiment, either an electron beam or a laser is used to ionize the molecules formed on the surface when they are desorbed as the surface is heated.The ions resulting from this ionization of the desorbed molecules are detected and iden-tied by a time-of-ight mass spectrometer (TOFMS). 47he experimental apparatus consists of two vacuum chambers.The 'source' chamber has a base pressure of approximately 10 À7 Torr.This chamber houses a microwave discharge cell used to generate O atoms from O 2 .Previous experiments have determined an O 2 dissociation efficiency of approximately 20% from this source, 9 so O 2 is still the majority species in our "O atom beam".
In investigating the low-temperature reactions induced on surfaces by our O atom source one must be aware, as discussed before, 9 that ozone can be formed by reactions of O atoms with O 2 on the surface.However, the oxygen atom addition reactions we observe proceed efficiently at surface temperatures well above 30 K, where ozone is not present on the surface. 9This observation indicates that the oxygenated products we observe are not the result of the reaction between propyne and ozone.Supporting this analysis is the observation that ozone only reacts with unsaturated hydrocarbons at low temperatures following irradiation. 48e also note that no reaction is observed if no microwaves are applied to the O atoms source and simply O 2 is dosed onto the surface.
Gas from the O atom source, undergoes signicant differential pumping in the source chamber before being piped into the 'target' chamber.A second deposition line attached to the source chamber allows the dosing of stable molecules, in this case propyne, onto the cold surface in addition to the oxygen atoms.Again,  the propyne dosing line undergoes signicant differential pumping in the source chamber before being directed into the target chamber.This differential pumping arrangement serves to allow sufficient pressure in the microwave source for stable operation, whilst dosing the target surface at an acceptably low rate.The UHV target chamber has a base pressure of approximately 10 À10 Torr when the dosing gases are switched off, but has a pressure of approximately 10 À8 Torr during the dosing process.
The PTFE tubes which direct the reactants from the source chamber into the target chamber terminate close above a highly oriented pyrolytic graphite (HOPG) substrate.This substrate can be cooled to close to 10 K and heated to above 500 K. Propyne and O 2 /O $ are co-deposited on the HOPG substrate at a xed surface temperature (in the range from 14 K to 100 K).The uxes of the propyne and O atoms are 1.0 Â 10 15 cm À2 s À1 and 1.3 Â 10 14 cm À2 s À1 respectively. 9Aer one hour, the delivery lines are evacuated and the substrate is allowed to cool from the xed deposition temperature to 14.0 AE 0.5 K.A current is then passed through a tantalum strip heater to slowly raise the surface temperature to $200 K.As the substrate temperature increases, the molecules on the surface desorb at their specic sublimation temperatures and enter the source region of the TOFMS which is located in front of the sample.It should be noted that the dosing uxes employed in our experiments rapidly result in a multilayer propyne/O 2 /O ice.Thus, formally the reactivity we observe should be considered representative of this surface.However, since the propyne molecules are physisorbed at the surface, we can assume that their electronic structure is largely unperturbed upon adsorption, and thus the reactivity we observe should be broadly representative of physisorbed propyne, irrespective of the precise nature of the surface on which the molecule is adsorbed.
Molecules desorbed from the surface are ionized in the source region of the TOFMS by a pulsed beam of 200 eV electrons running with a duty cycle of 32 ms.The electron pulses have a duration of approximately 1 ms.Aer each pulse of electrons crosses the source, a voltage is applied to the repeller plate of the TOFMS to accelerate any ions formed in the spectrometer's source region towards the detector for analysis.The detector is a pair of multichannel plates (MCPs), the output from which is amplied and discriminated before being passed to a timeto-digital convertor (TDC).The TDC is triggered by the same pulse generator that controls the electron gun and the repeller plate voltage.Times from the TDC are recorded as a histogram of time-of-ight versus counts, i.e. a mass spectrum.
During the TPD phase, as the temperature of the surface is slowly increased, a new mass spectrum is generated every second.At the end of each TPD experiment, these mass spectra are combined to generate a two-dimensional data set of ion intensity at each mass-to-charge (m/z) ratio against the substrate temperature during the heating process.This dataset can then be used to generate, for example, a mass spectrum for a particular temperature range or to see how the intensity of a particular ion (at a given m/z) varies as the surface temperature increases.In this work we recorded these TPD datasets aer dosing propyne and O atoms onto the surface at a range of surface temperatures between 14 K and 100 K, performing four separate dosing/TPD measurements at each surface temperature investigated.

Results
Whenever propyne and O atoms are allowed to react at surface temperatures below 100 K, we observe signals at mass-to-charge ratios (m/z) of 56 and 72 in the TOF mass spectra recorded during the TPD phase (Figure 1).These two signals possess different TPD proles, with the signal at m/z ¼ 56 appearing at lower desorption temperatures.Thus, it seems evident that these two mass spectral signals correspond to two different products.These two signals in the mass To determine the yield of each of these products, the total integrated mass spectral signals at m/z ¼ 56 and m/z¼72 in each TPD measurement is determined.This procedure involves selecting the temperature range over which each product desorbs and summing the mass spectra in this temperature range for one TPD measurement to form an integrated mass spectrum.A background mass spectrum, recorded in an identical experiment carried out with the microwave discharge off, is then subtracted from this integrated mass spectrum.This background correction procedure conrms the product signals are a result of the oxygen atoms generated in the discharge.The intensities of the resulting mass spectral peaks are determined from this "corrected" integrated spectrum by peak tting.This integration procedure is repeated for each separate experiment, and the results averaged for each surface temperature at which the dosing was carried out.The above integration procedure gives the relative ion signals at m/z ¼ 56 and m/z ¼ 72 in terms of ion counts at the MCP detector.To compare with the kinetic model, with which we interpret our results (see below), we require the product yields in terms of the number of molecules formed on the surface.To transform our mass spectral ion counts to these absolute units requires calibration of the mass spectral intensities.To carry out this transformation, as reported before, 8,9 we measure the mass spectrum resulting from a TPD dataset recorded following the adsorption of a known dose of propyne when the surface is held at 12 K.The integrated propyne mass spectral intensity in such a TPD dataset allows us to calculate the proportionality constant, the detection efficiency, between the dose of molecules and the mass spectral signal, assuming a sticking probability of 1 in the calibration experiment.Such a sticking coefficient is likely to be an excellent characterization of the interaction of propyne with the 12 K surface.The detection efficiency we measure for propyne must now be adjusted to represent the detection of the two product molecules.Since the ionization cross sections at 200 eV are not known for propyne and the two products (C 3 H 4 O and C 3 H 4 O 2 ), we must assume they are the same for all three species.Since we integrate just the parent ion signal in the mass spectrum, we must also adjust the detection efficiency to reect the parent to fragment ion ratio in the mass spectrum of propyne and in the mass spectra of the product species.For the single addition product we used the parent to fragment ion ratio in a propanal mass spectrum measured in our apparatus.For the double addition product we estimated the fragment to parent ratio in the standard mass spectrum of methyl glyoxal from the NIST reference database. 49The above procedure allows us to estimate the yield of the product ions on an absolute scale from our TPD spectra, and report the yield as a function of the dosing temperature.
We initially report the total product yield, the summed yields of the single and double addition products, as a function of dosing surface temperature (Figure 2).These experimental results show that below 30 K the yield of C 3 H 4 O gently increases with decreasing temperature.Above a surface temperature of 30 K the product yield increases to a maximum at 50 K and then decreases to zero at a surface temperature of 100 K.We also report the individual yields for the single addition product and the double addition product (Figure 3). Figure 3 shows that the same general trend with surface temperature is exhibited by the yields of both of the individual addition products as was described for the total product yield (Figure 2).We also note the yields for the single and double addition products are very similar below a surface temperature of 40 K. Above a surface temperature of 40K the yield of the single addition product is much larger than that of the double addition product (Figure 3).

Data reduction
To extract estimates of kinetic parameters from our experimental product yields, we t the experimental data with a simple kinetic model.This kinetic model has been described before in the literature, 8,9 but has also been extended here to model sequential addition of two oxygen atoms to propyne.
The kinetic model involves the two prototypical mechanisms by which the reaction of O atoms with propyne can proceed on the surface.Firstly, the LH mechanism where both reactants are adsorbed and thermalised on the substrate: Secondly, the ER mechanism, where an adsorbed species undergoes reaction with a gas-phase partner: In our experiments, as we have described before, 8,9 the molecular species should have a signicantly larger adsorption energy than the O atom.Thus, to simplify our model we only consider the form of ER reactivity where the O atom is the gas-phase reactant and the propyne molecule is adsorbed on the surface.This is because when both reactants are adsorbed on the surface, the rate for LH reactivity should dominate.Given the arguments above, the rates for the LH and ER mechanisms are therefore:

(ads)] [O(ads)]
(3) where r i is the rate of reaction and k i is the rate coefficient.In these equations i denotes the LH or ER mechanism, [C 3 H 4 (ads)] and [O(ads)] are the surface concentrations (in molecules cm À2 ) of the reactants absorbed on the surface and F O is the deposition uence of oxygen atoms.
At the low surface temperatures employed in this study, and given the inertness of the surface, it is safe to assume that the reactants are physisorbed and, as physisorption is dominant, there is no limit on the number of accessible adsorption sites.That is, once the rst layer of molecular ice, comprising reactants and products, is deposited a second layer will then readily form on top of the rst layer, eventually building up a multi-layer ice.In our model, we assume only the uppermost monolayer is accessible to the incident reactants, that is, we assume the reactants radicals cannot penetrate into the ice.Thus, the maximum surface coverage of each species is constrained to 10 15 molecules cm À2 in our model, 8,9 the standard value for an accessible monolayer.
The temperature variation of the rate coefficients for the two surface reactions can be described by an Arrhenius expression: where i again denotes the LH or ER mechanism.In equation 5, A i is a pre-exponential factor and E i is the activation energy of the pathway.To simplify the modelling and reduce the number of free parameters, the activation energy is assumed to be the same for both the LH and ER mechanisms.
Since both LH and ER reactions can occur simultaneously, the overall rate r of formation of the single addition product is given by: As we see in Eq. 6, the reaction rate is dependent upon the surface concentration of propyne and oxygen atoms.For propyne the three factors that control the surface concentration are given in Equation 7: In equation 7, F propyne is the ux of propyne onto the substrate and r Des,propyne is the rate of propyne desorption from the surface.The surface concentration of propyne can be found by integrating Equation 7with respect to deposition time t.As described before, the propyne ux, F propyne , is estimated experimentally from the propyne ow rate and the pressure/pumping characteristics of the vacuum chamber. 33,35The depletion of propyne due to the reaction r is calculated using equation 6.Finally, the desorption rate of propyne, r Des,propyne , can be evaluated using a second Arrhenius equation at each surface temperature: where A Des,propyne is the pre-exponential factor for the desorption of propyne and E Des,propyne is the desorption energy for propyne.As discussed before, since the desorption of the multi-layer ice in our experiment just reveals another layer of the ice for reactions, only desorption in the monolayer regime will affect the surface concentrations and hence rst-order desorption kinetics are used in Equation ( 9). 33,35he concentration of oxygen atoms is evaluated analogously to the concentration of propyne, allowing for the dissociation efficiency of the source. 33,35esorption of the products of the reactions are not considered in our model.Such a simplication is justiable as the products from the surface reaction between propyne and oxygen are heavier and more polar than the reactants and so should have larger desorption energies.This conclusion is conrmed by the TPD proles we observe for the products.
As discussed above, the different TPD proles of the two products, m/z ¼ 56 and m/z ¼ 72, conrm that the single addition signal, m/z ¼56, is from a separate product to that responsible for the double addition signal at m/z ¼ 72.That is, the signal at m/z ¼ 56 is clearly not a mass spectral fragment of the species at m/z ¼ 72.Given that the reaction conditions rule out a concerted three-body reaction, it seems clear that the reaction that forms the double addition product occurs in two steps.Firstly, a single oxygen atom adds to propyne to form a C 3 H 4 O species.Then, in the second step, a second oxygen atom adds to the C 3 H 4 O to form the double addition product (C 3 H 4 O 2 ).Hence, all the double addition product we detect was at some point during the reaction a single addition product.Thus, to extract the kinetic parameters for the initial single addition step, we t our kinetic model to the summed experimental yields of the single and double addition products.
Some preliminary considerations need to be addressed before tting the kinetic model to the experimental data (Figure 2) in order to extract the kinetic parameters.The experimental yield curves we determined in previous investigations of the reactions of O displayed unambiguous evidence of the parallel operation of both the LH and ER mechanisms. 9Specically, a clear peak in the yield was observed at surface temperatures just below the desorption temperature of the O atoms.At higher surface temperatures another peak in the yield was observed due to efficient operation of the ER pathway.The experimental yield prole for O + propyne (Figure 2) does not show this distinct double peaked structure.However it is apparent that the LH mechanism is still a signicant contributor to the yield of C 3 H 4 O, since we see a sharp drop in the yield at surface temperatures where the O atoms no longer stick efficiently to the surface, above approximately 50K.Our kinetic model does not give a good t to the experimental data using the LH mechanism alone since above 50 K, where the O atoms no longer stick efficiently, there is still a signicant yield of C 3 H 4 O.Therefore, although there is no observed double peaked structure in the experimental yield, both the LH and ER processes appear to be contributing to the product yield.Thus we have used both the LH and ER reactions in our kinetic model to t the experimental data.
To model the yield of C 3 H 4 O at a given surface temperature, we numerically integrate Equation 6for a time period equal to the experimental dosing period.The total number of product molecules formed can then be compared with the experimental data.To achieve a t with the experimental data, we can vary the preexponential factors and the activation energy for the chemical reaction and also the desorption energies of the reactants.Again to constrain the number of free parameters, the pre-exponential factors for the desorption of the two reactants from the surface are kept xed.Usually such pre-exponential factors (A Des,O and A Des,propyne ) are taken as the vibrational frequency of the adsorbate-surface bond.In previous work we employed a value for A Des,O of 3.10Â10 12 s À1 and we employ this value again here. 9,50To the best of our knowledge there are no available surface vibrational frequencies for propyne, so for A Des,propyne we employ the value of 2.33 x 10 12 s À1 , which is that calculated for an acetylene-graphite bond. 9,51e rst attempted to t the kinetic model to the experimental data, constraining the pre-exponential factors for the LH and ER mechanisms to be equal (A LH ¼ A ER ) and also constraining the activation energies (E LH ¼ E ER ) for both mechanisms to be equal.Under these constraints the t to the experimental data is poor.From this poor t of the highly constrained model we concluded that we needed a larger rate coefficient for the ER mechanism than for the LH mechanism.If the activation energy for surface diffusion of the O atoms is not signicant, we would expect the activation energies for the LH and ER mechanisms to be equal.Thus, to allow a larger ER rate coefficient in the model we permit A ER to take a different value to A LH .This additional degree of freedom allows us to achieve a much better, although not perfect, t to the experimental data (Figure 2) as discussed below.In attempting this tting it is clear that the barrier for the LH Open Access Article.Published on 21 January 2014.Downloaded on 10/29/2023 12:36:09 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.mechanism is well constrained by the form of the data to 160 AE 10 K. Similarly, the tting shows the propyne binding energy is also well constrained by the experimental data; we comment further on the O atom binding energy below.The desorption energies for propyne and oxygen are found to be 20.8AE 0.3 kJ mol À1 and 14.0 AE 1.2 kJ mol À1 respectively.
The parameters for this satisfactory t are reported in Table 1 and the t is shown graphically in Figure 2. As is clear in Figure 2, the reported t still does not satisfactorily reproduce experimental yield at T ¼ 70 K.Possible physical explanations for this poor t at a surface temperature of 70 K include the presence of more than one oxygen atom binding site, as discussed further below.
As mentioned earlier, we have also extended our model to attempt to extract kinetic parameters for the double addition reaction.This extension involves the addition of two further rate equations for the formation of the double addition product from the single addition product: Again, we integrate these equations to derive the yield of double addition product at the different surface temperatures, allowing the single addition product to form with the kinetic parameters listed in Table 1.Under these constraints we nd that to t the yield of the double addition product (Figure 3) requires a negative activation energy (À60 AE 20 K).

Discussion
The oxygen atom desorption energy which gives the best t to the experimental data in our kinetic model is 15.3 kJ mol À1 .The interaction energy of an oxygen atom and pyrene (representative of a bridge site in graphite) is calculated to be 11.6 kJ mol À1 . 50As discussed above, in our experiments the graphite substrate is saturated aer a matter of seconds during dosing.This means that although our substrate is graphite, the majority of the oxygen atoms are interacting with a propyne-oxygen matrix.Previous experimental work has shown that the desorption energy of an oxygen atom from an ethene-oxygen matrix and from a propene-oxygen matrix is, similarly to the oxygen-graphite system, about 12 kJ mol À1 . 9The oxygen atom desorption energy which best ts our data is markedly larger than these previous values for the binding of oxygen atoms to small organic molecules.If we use a value of 12.9 kJ mol À1 in our simulation (Figure 4), we This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
satisfactorily reproduce the rise in the yield of the product over the surface temperature range from 30 K to 50 K, but a larger desorption energy, or a component of the oxygen atoms possessing a larger desorption energy, is required to t the data point at a surface temperature of 60 K (Figure 2).From this analysis it appears that the experimental data might be best represented with at least two desorption energies for the oxygen atoms.Given the consequent increase in the number of free parameters, we are reluctant to use such a desorption energy distribution for the oxygen atoms as we feel we risk over-interpreting our experimental data.We simply report that the oxygen atom desorption energy is not that well dened by the form of our experimental data, and is best represented by a value of 14 AE 2 kJ mol À1 , this value perhaps representing a distribution of desorption energies.The propyne desorption energy extracted using our kinetic model is 20.8 AE 0.3 kJ mol À1 .To our knowledge there are no literature values for the desorption energy of propyne from interstellar ice analogues or a graphite surface.The value of this desorption energy agrees nicely with the observed trend that an alkyne's desorption energy is smaller than that of its alkene equivalent. 51Specically, the work of Rubes et al. 51 shows that ethene's desorption energy from a graphite (0001) surface is 17.36 AE 0.0361 kJ mol À1 whilst acetylene's desorption energy from the same surface is 14.67 AE 0.0273 kJ mol À1 .In agreement with this trend, previous work has shown that the desorption energy of propene from propene ice is 21.4 AE 0.3 kJ mol À1 , 9 a value that is slightly larger than the equivalent value we extract here for the desorption of propyne.
Our modelling yields a value for the activation energy of single O atom addition to propyne of 160 AE 10 K, a value larger than the comparable reaction of O atoms with propene (145 AE 10 K). 9 Such an increase in the activation energy is in accord with chemical intuition, as the reactive C-C bond in propyne is stronger than that in propene.Such a trend is also supported by the fact we observe no reaction between O atoms and acetylene at a range of surface temperatures between 14 K and 100 K. Previous work has shown that the activation energy (for the LH mechanism) for the reaction between ethene and O atoms is 190 AE 45 K. 9 The absence of any signals in our experiments for the addition of O atoms to acetylene indicates the barrier for this process is signicantly above 200 K.Such a conclusion is in accord with the trend of a rise in the activation energy, for the reaction with O atoms, on changing the reactant from the alkene to the corresponding alkyne.More generally, our results show that, under our surface conditions, the methyl substituted organic reactant (propene and propyne) is more reactive with oxygen atoms that the corresponding smaller molecule (ethene and acetylene).A similar trend is observed in combustion chemistry 43 where the activation energy for the reaction of O atoms with propyne is smaller than that for acetylene.
The pre-exponential factors we extract from the t of the kinetic model to the experimental data are 9.8 AE 0.2 x 10 À17 cm 1 molecule À1 s À1 and 2.45 AE 0.4 x 10 À16 cm 1 molecule À1 s À1 for the LH and ER reactions respectively.These pre-exponential factors are approximately a factor of 10 smaller than we report for the reaction between oxygen atoms and alkenes. 9This trend in pre-exponential factors is also present in the gas phase, where the pre-exponential factors for the reaction of O atoms with alkynes are smaller than those for the reaction with alkenes. 38,52As encompassed by these pre-exponential factors, and discussed above, to t our experimental data (Figure 2), the rate coefficient for the ER reaction must be larger than the rate coefficient for the LH process.One reason for this difference might be the fact that in the ER process the O atoms are not thermalized with the surface.Indeed, in our experiment the O atoms will possess a kinetic energy distribution representative of 300 K.One way to test such a hypothesis experimentally would be to cool the incident O atom beam, but the methodology for such cooling is not well established experimentally.A greater rate constant for the ER mechanism, in comparison with the LH mechanism, has also been reported for the reactions between O atoms and O 2 on cold surfaces and again assigned to non-thermal O atoms arriving at the surface. 36s discussed above, at surface temperatures below 30 K the yield of the single and double addition products increases as the surface temperature decreases.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Our simple kinetic model does not account for the increase in the product yield at temperatures below 30 K. A similar reactive pathway, with an increasing yield at low temperatures, has been observed before for reactions of O atoms with CS 2 . 8eactivity at such surface temperatures is particularly pertinent to the ISM.Since the rate of this low temperature channel rises with falling surface temperature, it could be that this low temperature process is effectively barrierless.Since this low temperature reactivity is only present between surface temperatures of 14 K and 30 K, performing any sort of kinetic t for our data is not statistically meaningful.Hence we simply report that the rate coefficient at 14 K for the formation of the single addition product is (6 AE 2)Â10 À19 cm 1 molecule À1 s À1 .We note that this low temperature rate is at least comparable to that observed for the reaction of O atoms with CS 2 , 8 and that tunneling processes have been recently reported for O atoms at low surface temperatures. 36ur key nding for the addition of a second oxygen atom to the C 3 H 4 O single addition product is that this is a "barrierless" process.Such a conclusion is not unexpected since the initial addition of an oxygen atom to propyne is likely to initially result in a radical species.The subsequent reaction of this radical with a further oxygen atom might be expected to exhibit a negative activation energy. 53s shown in Figure 3, the yields of the single (C 3 H 4 O) and double (C 3 H 4 O 2 ) addition products are similar at surface temperatures below 40 K. Above 40 K there is a dramatic increase in the yield of the single addition product with respect to the double addition product.The observed increase in the yield of the single addition product at surface temperatures above 40 K perhaps implies that the reactive intermediate resulting from single addition isomerizes above 40 K to yield a more stable isomer of C 3 H 4 O, an isomer that is less susceptible to attack by a second oxygen atom.
Our previous work has shown that the surface reaction between ethene (C 2 H 4 ) and O atoms results in the formation of ethylene oxide (Structure I, Figure 5). 9imilarly it is reasonable to propose that the rst intermediate formed by the reaction of O atoms with propyne is an oxirene (II, Figure 5).We would expect this oxirene to be reactive, but potentially to be stabilized on a low temperature surface.Such an interpretation agrees nicely with our deduction above that at surface temperatures below 40 K we have a reactive single addition product on the surface, but above a surface temperature of 40 K this reactive single addition product isomerizes to a less reactive species.There are many (10) possible isomerisation products for the reactive oxirene, but several of these would also be highly reactive.However, one stable molecule which can be formed from the isomerisation of the oxirene is propenal (III, Figure 5), the isomerisation occurring via a pathway involving a carbene. 54It has been shown that the rearrangement of an oxirene to a ketocarbene is barrierless. 55Indeed, the rst investigation of the gas-phase reactivity of O atoms with propyne concluded that, if energy could be efficiently lost from the initial addition product, the ketocarbene (VI Figure 5) formed from the oxirene will relax to methyl ketene (VII Figure 5) or to propenal (III Figure 5). 38he structure of the double addition product is more speculative.An obvious candidate for the initial product of the second addition is the bicyclic structure IV (Figure 5).This strained intermediate might be expected to be very reactive and will probably rapidly rearrange.There are a number of possible rearranged structures for this primary bicyclic adduct, with one stable rearrangement product being methyl glyoxal (V Figure 5).Our experiments show that on interstellar dust grain analogues, under UHV conditions and at low temperature, oxygen atoms can readily add to the carboncarbon triple bond of a propyne molecule.The reaction is detectable at surface temperatures of 14 K and above, and is most efficient at 50 K.We estimate our experiments to be the equivalent of approximately 10 5 -10 6 years of exposure to O atoms in an interstellar cloud. 9,33,35It is therefore possible that this pathway is active in interstellar clouds.In addition, given the observed kinetics, the yield from the reaction of propyne with oxygen atoms may also be signicant when the cloud begins to warm. 56ur experiments reveal a marked difference in the heterogeneous reactivity of acetylene and propyne on low temperature surfaces, a difference in reactivity that may be pertinent in the ISM.Indeed, we note that the reactions of molecular ices with incident oxygen atoms are not well represented in recent descriptions of gasgrain reaction networks. 57s discussed in the Introduction, propyne has been observed in a variety of environments in the ISM, and in these environments the surface chemistry we report in this article may be relevant.We also previously discussed the need for an O atom sink in interstellar dust clouds, to account for the observed gas-phase depletion.We postulate that the stable molecule resulting from the single addition of O atoms to propyne on a cold surface is propenal (III Figure 5).This reaction, along with other O atom reactions, may contribute to that O atom sink.Propenal (C 3 H 4 O) has been observed in the ISM towards the star forming region Sagittarius B2(N). 58Hollis et al. postulate that the formation of interstellar propenal is due to hydrogen addition to propynal. 58However, the work reported in this paper suggests that, in addition, the surface chemistry of propyne may also be a source of propenal.Clearly modelling work is required to support this suggestion.Encouragingly, propyne has also been detected in Sagittarius B2 59 in accord with the idea that propenal may be formed from the heterogeneous oxidation of propyne by oxygen atoms.
To the best of our knowledge no molecule with the empirical formula of the double addition product (C 3 H 4 O 2 ) has been observed in the interstellar medium.However, given the observed facility of the double addition process, our experimental data would suggest that at least one isomer of C 3 H 4 O 2 , perhaps methylglyoxal, could be present during the warm-up phase of an interstellar cloud where propyne ice has been thermally processed by oxygen atoms.

Conclusions
This paper reports the rst laboratory investigation of the heterogeneous reaction of propyne with oxygen atoms at temperatures relevant to the ISM.Our data shows this reaction, which initially forms an adduct due to single addition (C 3 H 4 O) proceeds efficiently at surface temperatures below 100 K.As the surface temperature decreases from 100 K the yield of the reaction increases to reach a maximum at 50 K and then falls to a minimum at 30K.Modelling the yield of the reaction at surface temperatures above 30 K allows the extraction of an activation energy for the reaction of 160 AE 10 K. We also detect a product corresponding to the addition of two oxygen atoms to propyne.This double-addition product appears to form in a sequential reaction of the single addition product with another oxygen atom.The addition of this second oxygen atom to the single addition product appears to be barrierless.We propose structures for the intermediates and products of the O + propyne reaction which are consistent with our observations and the known behaviour of these organic compounds.We postulate that the addition of oxygen atoms to propyne to form propenal could act as a source of this molecule in the ISM.

Fig. 1
Fig. 1 Sections of representative mass spectra recorded during the TPD phase showing peaks for m/z ¼ 56 and m/z ¼ 72, corresponding to the single and double addition products.See text for details.

Fig. 2
Fig. 2 Summed yield of the single and double addition product, formed following the codeposition of propyne (C 3 H 4 ) and O atoms, as a function of surface temperature.Squares: experimental data; solid line: model; dash LH mechanism; dot ER mechanism.The error bars associated with the experimental results represent two standard deviations from four repeats at each surface temperature.The kinetic model used to derive the fit shown employs the paramenters listed inTable 1 and an O atom desorption energy of 15.3 kJ mol À1 .

Fig. 3
Fig. 3 Experimental yield of the single (dashed line) and double (dotted line) addition products following the co-deposition of propyne (C 3 H 4 ) and O atoms.The error bars associated with the experimental results represent two standard deviations from four experiments at each surface temperature.The lines linking the points serve only to guide the eye.

Fig. 4
Fig. 4 Summed yield of the single and double addition product, formed following the codeposition of propyne (C 3 H 4 ) and O atoms, as a function of surface temperature.Squares: experimental data; solid line: model; dash LH mechanism; dot ER mechanism.The error bars associated with the experimental results represent two standard deviations from four repeats at each surface temperature.The kinetic model used to derive the fits shown employs the paramenters listed in Table1and an O atom desorption energy of 12.9 kJ mol À1 .

Fig. 5
Fig. 5 Organic structures relevant to the discussion of the reaction of O atoms with propyne.

Table 1
Kinetic parameters characterizing the single addition of O atoms to propyne, to form the single addition product C 3 H 4 O, as a function of surface temperature.As discussed in the text, these parameters have been extracted by fitting a kinetic model to the experimental data we record for this surface reaction(E i /R) /K i ¼LH or ER 10 16 A LH / cm 2 molecule À1 s À1 10 16 A ER / cm 2 molecule À1 s À1 Open Access Article.Published on 21 January 2014.Downloaded on 10/29/2023 12:36:09 AM.