Formation of interstellar complex organic molecules is currently thought to be dominated by the barrierless coupling between radicals on the interstellar icy grain surfaces. Previous standard DFT results on the reactivity between CH$_3$ and HCO on amorphous water surfaces, showed that formation of CH$_4$ + CO by H transfer from HCO to CH$_3$ assisted by water molecules of the ice was the dominant channel. However, the adopted description of the electronic structure of the biradical (i.e., CH$_3$/HCO) system was inadequate (without the broken-symmetry (BS) approach). In this work, we revisit the original results by means of BS-DFT both in gas phase and with one water molecule simulating the role of the ice. Results indicate that adoption of BS-DFT is mandatory to describe properly biradical systems. In the presence of the single water molecule, the water-assisted H transfer exhibits a high energy barrier. In contrast, CH$_3$CHO formation is found to be barrierless. However, direct H transfer from HCO to CH$_3$ to give CO and CH$_4$ presents a very low energy barrier, hence being a potential competitive channel to the radical coupling and indicating, moreover, that the physical insights ofthe original work remain valid.

Glycoaldehyde, ethylene glycol, and methyl formate are complex organic molecules that have been observed in dark molecular clouds. Because there is no efficient gas-phase route to produce these species, it is expected that a low-temperature surface route existst that does not require energetic processing. CO hydrogenation experiments at low temperatures showed that this is indeed the case. Glyoxal can form through recombination of two HCO radicals and is then further hydrogenated. Here we aim to constrain the methyl formate, glycolaldehyde, and ethylene glycol formation on the surface of interstellar dust grains through this cold and dark formation route. We also probe the dependence of the grain mantle composition on the initial gas-phase composition and the dust temperature. A full CO hydrogenation reaction network was built based on quantum chemical calculations for the rate constants and branching ratios. This network was used in combination with a microscopic kinetic Monte Carlo simulation to simulate ice chemistry, taking into account all positional information. After benchmarking the model against CO-hydrogenation experiments, simulations under molecular cloud conditions were performed. COMs are formed in all interstellar conditions we studied, even at temperatures as low as 8 K. This is because the HCO + HCO reaction can occur when HCO radicals are formed close to each other and do not require to diffuse. Relatively low abundances of methyl formate are formed. The final COM abundances depend more on the H-to-CO ratio and less on temperature. Only above 16 K, where CO build-up is less efficient, does temperature start to play a role. Molecular hydrogen is predominantly formed through abstraction reactions on the surface. Our simulations are in agreement with observed COM ratios for mantles that have been formed at low temperatures.

In this work, we reexamine sulfur chemistry occurring on and in the ice mantles of interstellar dust grains, and report the effects of two new modifications to standard astrochemical models; namely, (a) the incorporation of cosmic ray-driven radiation chemistry and (b) the assumption of fast, non-diffusive reactions for key radicals in the bulk. Results from our models of dense molecular clouds show that these changes can have a profound influence on the abundances of sulfur-bearing species in ice mantles, including a reduction in the abundance of solid-phase H$_2$S and HS, and a significant increase in the abundances of OCS, SO$_2$, as well as pure allotropes of sulfur, especially S$_8$. These pure-sulfur species - though nearly impossible to observe directly - have long been speculated to be potential sulfur reservoirs and our results represent possibly the most accurate estimates yet of their abundances in the dense ISM. Moreover, the results of these updated models are found to be in good agreement with available observational data. Finally, we examine the implications of our findings with regard to the as-yet-unknown sulfur reservoir thought to exist in dense interstellar environments.

Complex organic molecules (COMs) have been detected in the gas-phase in cold and lightless molecular cores. Recent solid-state laboratory experiments have provided strong evidence that COMs can be formed on icy grains through 'non-energetic' processes. In this contribution, we show that propanal and 1-propanol can be formed in this way at the low temperature of 10 K. Propanal has already been detected in space. 1-propanol is an astrobiologically relevant molecule, as it is a primary alcohol, and has not been astronomically detected. Propanal is the major product formed in the C2H2 + CO + H experiment, and 1-propanol is detected in the subsequent propanal + H experiment. The results are published in Qasim et al. (2019c). ALMA observations towards IRAS 16293-2422B are discussed and provide a 1-propanol:propanal upper limit of < 0.35 - 0.55, which are complemented by computationally-derived activation barriers in addition to the performed laboratory experiments.

1-propanol (CH3CH2CH2OH) is a three carbon-bearing representative of primary linear alcohols that may have its origin in the cold dark cores in interstellar space. To test this, we investigated in the laboratory whether 1-propanol ice can be formed along pathways possibly relevant to the prestellar core phase. We aim to show in a two-step approach that 1-propanol can be formed through reaction steps that are expected to take place during the heavy CO freeze-out stage by adding C2H2 into the CO + H hydrogenation network via the formation of propanal (CH3CH2CHO) as an intermediate and its subsequent hydrogenation. Temperature programmed desorption-quadrupole mass spectrometry (TPD-QMS) is used to identify the newly formed propanal and 1-propanol. Reflection absorption infrared spectroscopy (RAIRS) is used as a complementary diagnostic tool. The mechanisms that can contribute to the formation of solid-state propanal and 1-propanol, as well as other organic compounds, during the heavy CO freeze-out stage are constrained by both laboratory experiments and theoretical calculations. Here it is shown that recombination of HCO radicals, formed upon CO hydrogenation, with radicals formed upon C2H2 processing - H2CCH and H3CCH2 - offers possible reaction pathways to solid-state propanal and 1-propanol formation. This extends the already important role of the CO hydrogenation chain in the formation of larger COMs (complex organic molecules). The results are used to compare with ALMA observations. The resulting 1-propanol:propanal ratio concludes an upper limit of < 0:35-0:55, which is complemented by computationally-derived activation barriers in addition to the experimental results.

Formic acid (HCOOH) and carbon dioxide (CO2) are simple species that have been detected in the interstellar medium. The solid-state formation pathways of these species under experimental conditions relevant to prestellar cores are primarily based off of weak infrared transitions of the HOCO complex and usually pertain to the H2O-rich ice phase, and therefore more experimental data are desired. In this article, we present a new and additional solid-state reaction pathway that can form HCOOH and CO2 ice at 10 K 'non-energetically' in the laboratory under conditions related to the "heavy" CO freeze-out stage in dense interstellar clouds, i.e., by the hydrogenation of an H2CO:O2 ice mixture. This pathway is used to piece together the HCOOH and CO2 formation routes when H2CO or CO reacts with H and OH radicals. Temperature programmed desorption - quadrupole mass spectrometry (TPD-QMS) is used to confirm the formation and pathways of newly synthesized ice species as well as to provide information on relative molecular abundances. Reflection absorption infrared spectroscopy (RAIRS) is additionally employed to characterize reaction products and determine relative molecular abundances. We find that for the conditions investigated in conjunction with theoretical results from the literature, H+HOCO and HCO+OH lead to the formation of HCOOH ice in our experiments. Which reaction is more dominant can be determined if the H+HOCO branching ratio is more constrained by computational simulations, as the HCOOH:CO2 abundance ratio is experimentally measured to be around 1.8:1. H+HOCO is more likely than OH+CO (without HOCO formation) to form CO2. Isotope experiments presented here further validate that H+HOCO is the dominant route for HCOOH ice formation in a CO-rich CO:O2 ice mixture that is hydrogenated. These data will help in the search and positive identification of HCOOH ice in prestellar cores.

A number of recent experimental studies have shown that solid-state complex organic molecules (COMs) can form under conditions that are relevant to the CO freeze-out stage in dense clouds. In this work, we show that alcohols can be formed well before the CO freeze-out stage (i.e., in the H2O-rich ice phase). This joint experimental and computational investigation shows that the isomers, n- and i-propanol (H3CCH2CH2OH and H3CCHOHCH3) and n- and i-propenol (H3CCHCHOH and H3CCOHCH2), can be formed in radical-addition reactions starting from propyne (H3CCCH) + OH at the low temperature of 10 K, where H3CCCH is one of the simplest representatives of stable carbon chains already identified in the interstellar medium (ISM). The resulting average abundance ratio of 1:1 for n-propanol:i-propanol is aligned with the conclusions from the computational work that the geometric orientation of strongly interacting species is influential to the extent of which 'mechanism' is participating, and that an assortment of geometries leads to an averaged-out effect. Three isomers of propanediol are also tentatively identified in the experiments. It is also shown that propene and propane (H3CCHCH2 and H3CCH2CH3) are formed from the hydrogenation of H3CCCH. Computationally-derived activation barriers give additional insight into what types of reactions and mechanisms are more likely to occur in the laboratory and in the ISM. Our findings not only suggest that the alcohols studied here share common chemical pathways and therefore can show up simultaneously in astronomical surveys, but also that their extended counterparts that derive from polyynes containing H3C(CC)nH structures may exist in the ISM. Such larger species, like fatty alcohols, are the possible constituents of simple lipids that primitive cell membranes on the early Earth are thought to be partially composed of.

The radicals HCO and CH$_3$ on carbon monoxide ice surfaces were simulated using density functional theory. Their binding energy on amorphous CO ice shows broad distributions, with approximative average values of 500 K for HCO and 200 K for CH$_3$. If they are located on the surface close to each other (3 to 4 \AA), molecular dynamics calculations based on density functional theory show that they can form acetaldehyde (CH$_3$CHO) or CH$_4$ + CO in barrier-less reactions, depending on the initial orientation of the molecules with respect to each other. In some orientations, no spontaneous reactions were found, the products remained bound to the surface. Sufficient configurational sampling, inclusion of the vibrational zero point energy, and a thorough benchmark of the applied electronic structure method are important to predict reliable binding energies for such weakly interacting systems. From these results it is clear that complex organic molecules, like acetaldehyde, can be formed by recombination reactions of radicals on CO surfaces.

In molecular clouds at temperatures as low as 10 K, all species except hydrogen and helium should be locked in the heterogeneous ice on dust grain surfaces. Nevertheless, astronomical observations have detected over 150 different species in the gas phase in these clouds. The mechanism by which molecules are released from the dust surface below thermal desorption temperatures to be detectable in the gas phase is crucial for understanding the chemical evolution in such cold clouds. Chemical desorption, caused by the excess energy of an exothermic reaction, was first proposed as a key molecular release mechanism almost 50 years ago. Chemical desorption can, in principle, take place at any temperature, even below the thermal desorption temperature. Therefore, astrochemical net- work models commonly include this process. Although there have been a few previous experimental efforts, no infrared measurement of the surface (which has a strong advantage to quantify chemical desorption) has been performed. Here, we report the first infrared in situ measurement of chemical desorption during the reactions H + H2S -> HS + H2 (reaction 1) and HS + H -> H2S (reaction 2), which are key to interstellar sulphur chemistry. The present study clearly demonstrates that chemical desorption is a more efficient process for releasing H2S into the gas phase than was previously believed. The obtained effective cross-section for chemical desorption indicates that the chemical desorption rate exceeds the photodesorption rate in typical interstellar environments.

Hydrogen addition and abstraction reactions play an important role as surface reactions in the buildup of complex organic molecules in the dense interstellar medium. Addition reactions allow unsaturated bonds to be fully hydrogenated, while abstraction reactions recreate radicals that may undergo radical-radical recombination reactions. Previous experimental work has indicated that double and triple C--C bonds are easily hydrogenated, but aldehyde -C=O bonds are not. Here, we investigate a total of 29 reactions of the hydrogen atom with propynal, propargyl alcohol, propenal, allyl alcohol, and propanal by means of quantum chemical methods to quantify the reaction rate constants involved. First of all, our results are in good agreement with and can explain the observed experimental findings. The hydrogen addition to the aldehyde group, either on the C or O side, is indeed slow for all molecules considered. Abstraction of the H atom of the aldehyde group, on the other hand, is among the faster reactions. Furthermore, hydrogen addition to C--C double bonds is generally faster than to triple bonds. In both cases, addition on the terminal carbon atom that is not connected to other functional groups is easiest. Finally, we wish to stress that it is not possible to predict rate constants based solely on the type of reaction: the specific functional groups attached to a backbone play a crucial role and can lead to a spread of several orders of magnitude in the rate constant.

The 29 reactions linking carbon monosulfide (CS) to methyl mercaptan (CH3SH) via ten intermediate radicals and molecules have been characterized with relevance to surface chemistry in cold interstellar ices. More intermediate species than previously considered are found likely to be present in these ices, such as trans- and cis-HCSH. Both activation and reaction energies have been calculated, along with low-temperature (T > 45~K) rate constants for the radical-neutral reactions. For barrierless radical-radical reactions on the other hand, branching ratios have been determined. The combination of these two sets of information provides, for the first time, quantitative information on the full H + CS reaction network. Early on in this network, that is, early on in the lifetime of an interstellar cloud, HCS is the main radical, while later on this becomes first CH2SH and finally CH3S.

Radical recombination has been proposed to lead to the formation of complex organic molecules (COMs) in CO-rich ices in the early stages of star formation. These COMs can then undergo hydrogen addition and abstraction reactions leading to a higher or lower degree of saturation. Here, we have studied 14 hydrogen transfer reactions for the molecules glyoxal, glycoaldehyde, ethylene glycol, and methylformate and an additional three reactions where \ce{CH_nO} fragments are involved. Over-the-barrier reactions are possible only if tunneling is invoked in the description at low temperature. Therefore the rate constants for the studied reactions are calculated using instanton theory that takes quantum effects into account inherently. The reactions were characterized in the gas phase, but this is expected to yield meaningful results for CO-rich ices due to the minimal alteration of reaction landscapes by the CO molecules. We found that rate constants should not be extrapolated based on the height of the barrier alone, since the shape of the barrier plays an increasingly larger role at decreasing temperature. It is neither possible to predict rate constants based only on considering the type of reaction, the specific reactants and functional groups play a crucial role. Within a single molecule, though, hydrogen abstraction from an aldehyde group seems to be always faster than hydrogen addition to the same carbon atom. Reactions that involve heavy-atom tunneling, e.g., breaking or forming a C--C or C--O bond, have rate constants that are much lower than those where H transfer is involved.

The hydrogen abstraction reaction between H and H$_2$S, yielding HS and H$_2$ as products, has been studied within the framework of interstellar surface chemistry. High-temperature rate constants up to 2000 K are calculated in the gas phase and are in agreement with previously reported values. Subsequently low-temperature rate constants down to 55 K are presented for the first time that are of interest to astrochemistry, i.e., covering both bimolecular and unimolecular reaction mechanisms. For this, a so-called implicit surface model is used. Strictly speaking, this is a structural gas-phase model in which the restriction of the rotation in the solid state is taken into account. The calculated kinetic isotope effects are explained in terms of difference in activation and delocalization. All rate constants are calculated at the UCCSD(T)-F12/cc-VTZ-F12 level of theory. Finally, we show that the energetics of the reaction is only affected to a small extent by the presence of H$_2$O or H$_2$S molecular clusters that simulate an ice surface, calculated at the MPWB1K/def2-TZVP level of theory.

On the surface of icy dust grains in the dense regions of the interstellar medium a rich chemistry can take place. Due to the low temperature, reactions that proceed via a barrier can only take place through tunneling. The reaction H + H$_2$O$_2$ $\rightarrow$ H$_2$O + OH is such a case with a gas-phase barrier of $\sim$26.5 kJ/mol. Still the reaction is known to be involved in water formation on interstellar grains. Here, we investigate the influence of a water ice surface and of bulk ice on the reaction rate constant. Rate constants are calculated using instanton theory down to 74 K. The ice is taken into account via multiscale modeling, describing the reactants and the direct surrounding at the quantum mechanical level with density functional theory (DFT), while the rest of the ice is modeled on the molecular mechanical level with a force field. We find that H$_2$O$_2$ binding energies cannot be captured by a single value, but rather depend on the number of hydrogen bonds with surface molecules. In highly amorphous surroundings the binding site can block the routes of attack and impede the reaction. Furthermore, the activation energies do not correlate with the binding energies of the same sites. The unimolecular rate constants related to the Langmuir-Hinshelwood mechanism increase as the activation energy decreases. Thus, we provide a lower limit for the rate constant and argue that rate constants can have values up to two order of magnitude larger than this limit.

OH radicals play a key role as an intermediate in the water formation chemistry of the interstellar medium. For example the reaction of OH radicals with H$_2$ molecules is among the final steps in the astrochemical reaction network starting from O, O$_2$, and O$_3$. Experimentally it was shown that even at 10 K this reaction occurs on ice surfaces. As the reaction has a high activation energy only atom tunneling can explain such experimental findings. In this study we calculated reaction rate constants for the title reaction on a water-ice I$_h$ surface. To our knowledge, low-temperature rate constants on a surface are not available in the literature. All surface calculations were done using a QM/MM framework (BHLYP/TIP3P) after a thorough benchmark of different density functionals and basis sets to highly accurate correlation methods. Reaction rate constants are obtained using instanton theory which takes atom tunneling into account inherently, with constants down to 110 K for the Eley-Rideal mechanism and down to 60 K for the Langmuir-Hinshelwood mechanism. We found that the reaction is nearly temperature independent below 80 K. We give kinetic isotope effects for all possible deuteration patterns for both reaction mechanisms. For the implementation in astrochemical networks, we also give fit parameters to a modified Arrhenius equation. Finally, several different binding sites and binding energies of OH radicals on the I$_h$ surface are discussed and the corresponding rate constants are compared to the gas-phase case.

Solid state astrochemical reaction pathways have the potential to link the formation of small nitrogen-bearing species, like NH3 and HNCO, and prebiotic molecules, specifically amino acids. To date, the chemical origin of such small nitrogen containing species is still not well understood, despite the fact that ammonia is an abundant constituent of interstellar ices toward young stellar objects and quiescent molecular clouds. This is mainly because of the lack of dedicated laboratory studies. The aim of the present work is to experimentally investigate the formation routes of NH3 and HNCO through non-energetic surface reactions in interstellar ice analogues under fully controlled laboratory conditions and at astrochemically relevant temperatures. This study focuses on the formation of NH3 and HNCO in CO-rich (non-polar) interstellar ices that simulate the CO freeze-out stage in dark interstellar cloud regions, well before thermal and energetic processing start to become relevant. We demonstrate and discuss the surface formation of solid HNCO through the interaction of CO molecules with NH radicals - one of the intermediates in the formation of solid NH3 upon sequential hydrogenation of N atoms. The importance of HNCO for astrobiology is discussed.

Hydroxylamine (NH2OH) is one of the potential precursors of complex pre-biotic species in space. Here we present a detailed experimental study of hydroxylamine formation through nitric oxide (NO) surface hydrogenation for astronomically relevant conditions. The aim of this work is to investigate hydroxylamine formation efficiencies in polar (water-rich) and non-polar (carbon monoxide-rich) interstellar ice analogues. A complex reaction network involving both final (N2O, NH2OH) and intermediate (HNO, NH2O, etc.) products is discussed. The main conclusion is that hydroxylamine formation takes place via a fast and barrierless mechanism and it is found to be even more abundantly formed in a water-rich environment at lower temperatures. In parallel, we experimentally verify the non-formation of hydroxylamine upon UV photolysis of NO ice at cryogenic temperatures as well as the non-detection of NC- and NCO-bond bearing species after UV processing of NO in carbon monoxide-rich ices. Our results are implemented into an astrochemical reaction model, which shows that NH2OH is abundant in the solid phase under dark molecular cloud conditions. Once NH2OH desorbs from the ice grains, it becomes available to form more complex species (e.g., glycine and beta-alanine) in gas phase reaction schemes.

This study focuses on the formation of two molecules of astrobiological importance - glycolaldehyde (HC(O)CH2OH) and ethylene glycol (H2C(OH)CH2OH) - by surface hydrogenation of CO molecules. Our experiments aim at simulating the CO freeze-out stage in interstellar dark cloud regions, well before thermal and energetic processing become dominant. It is shown that along with the formation of H2CO and CH3OH - two well established products of CO hydrogenation - also molecules with more than one carbon atom form. The key step in this process is believed to be the recombination of two HCO radicals followed by the formation of a C-C bond. The experimentally established reaction pathways are implemented into a continuous-time random-walk Monte Carlo model, previously used to model the formation of CH3OH on astrochemical time-scales, to study their impact on the solid-state abundances in dense interstellar clouds of glycolaldehyde and ethylene glycol.

We present a combined experimental and theoretical study focussing on the quantum tunneling of atoms in the reaction between CH4 and OH. The importance of this reaction pathway is derived by investigating isotope substituted analogs. Quantitative reaction rates needed for astrochemical models at low temperature are currently unavailable both in the solid state and in the gas phase. Here, we study tunneling effects upon hydrogen abstraction in CH4 + OH by focusing on two reactions: CH4 + OD -> CH3 + HDO and CD4 + OH -> CD3 + HDO. The experimental study shows that the solid-state reaction rate R(CH4 + OD) is higher than R(CD4 + OH) at 15 K. Experimental results are accompanied by calculations of the corresponding unimolecular and bimolecular reaction rate constants using instanton theory taking into account surface effects. From the work presented here, the unimolecular reactions are particularly interesting as these provide insight into reactions following a Langmuir-Hinshelwood process. The resulting ratio of the rate constants shows that the H abstraction (k(CH4 + OD)) is approximately ten times faster than D-abstraction (k(CD4 + OH)) at 65 K. We conclude that tunneling is involved at low temperatures in the abstraction reactions studied here. The unimolecular rate constants can be used by the modeling community as a first approach to describe OH-mediated abstraction reactions in the solid phase. For this reason we provide fits of our calculated rate constants that allow the inclusion of these reactions in models in a straightforward fashion.

The final step of the water formation network on interstellar grain surfaces starting from the H + O$_2$ route is the reaction between H and H$_2$O$_2$. This reaction is known to have a high activation energy and therefore at low temperatures it can only proceed via tunneling. To date, however, no rate constants are available at temperatures below 200 K. In this work, we use instanton theory to compute rate constants for the title reaction with and without isotopic substitutions down to temperatures of 50 K. The calculations are based on density functional theory, with additional benchmarks for the activation energy using unrestricted single-reference and multireference coupled-cluster single-point energies. Gas-phase bimolecular rate constants are calculated and compared with available experimental data not only for H + H$_2$O$_2$ $\rightarrow$ H$_2$O + OH, but also for H + H$_2$O$_2$ $\rightarrow$ H$_2$ + HO$_2$. We find a branching ratio where the title reaction is favored by at least two orders of magnitude at 114 K. In the interstellar medium this reaction predominantly occurs on water surfaces, which increases the probability that the two reactants meet. To mimic this one, two, or three spectator H2O molecules are added to the system. Eley-Rideal bimolecular and Langmuir-Hinshelwood unimolecular rate constants are presented here. The kinetic isotope effects for the various cases are compared to experimental data as well as to expressions commonly used in astrochemical models. Both the rectangular barrier and the Eckart approximations lead to errors of about an order of magnitude. Finally, fits of the rate constants are provided as input for astrochemical models.

Sets of systematic laboratory experiments are presented -- combining Ultra High Vacuum cryogenic and plasma-line deposition techniques -- that allow us to compare H/D isotopic effects in the reaction of H2O (D2O) ice with the hydroxyl radical OD (OH). The latter is known to play a key role as intermediate species in the solid-state formation of water on icy grains in space. The main finding of our work is that the reaction H2O + OD --> OH + HDO occurs and that this may affect the HDO/H2O abundances in space. The opposite reaction D2O + OH --> OD + HDO is much less effective, and also given the lower D2O abundances in space not expected to be of astronomical relevance. The experimental results are extended to the other four possible reactions between hydroxyl and water isotopes and are subsequently used as input for Kinetic Monte Carlo simulations. This way we interpret our findings in an astronomical context, qualitatively testing the influence of the reaction rates.

We have investigated the thermally induced proton/deuteron exchange in mixed amorphous H$_2$O:D$_2$O ices by monitoring the change in intensity of characteristic vibrational bending modes of H$_2$O, HDO, and D$_2$O with time and as function of temperature. The experiments have been performed using an ultra-high vacuum setup equipped with an infrared spectrometer that is used to investigate the spectral evolution of homogeneously mixed ice upon co-deposition in thin films, for temperatures in the 90 to 140 K domain. With this non-energetic detection method we find a significantly lower activation energy for H/D exchange -- $3840 \pm 125$ K -- than previously reported. Very likely this is due to the amorphous nature of the interstellar ice analogues involved. This provides reactive timescales ($\tau<10^4$ years at $T$ $>70$ K) fast enough for the process to be important in interstellar environments. Consequently, an astronomical detection of D$_2$O will be even more challenging because of its potential to react with H$_2$O to form HDO. Furthermore, additional experiments, along with previous studies, show that proton/deuteron swapping also occurs in ice mixtures of water with other hydrogen bonded molecules, in particular on the OH and NH moieties. We conclude that H/D exchange in ices is a more general process that should be incorporated into ice models that are applied to protoplanetary disks or to simulate the warming up of cometary ices in their passage of the perihelion, to examine the extent of its influence on the final deuteron over hydrogen ratio.

The formation of interstellar water has been commonly accepted to occur on the surfaces of icy dust grains in dark molecular clouds at low temperatures (10-20 K), involving hydrogenation reactions of oxygen allotropes. As a result of the large abundances of molecular hydrogen and atomic oxygen in these regions, the reaction H2 + O has been proposed to contribute significantly to the formation of water as well. However, gas phase experiments and calculations, as well as solid-phase experimental work contradict this hypothesis. Here, we use precisely executed temperature programmed desorption (TPD) experiments in an ultra-high vacuum setup combined with kinetic Monte Carlo simulations to establish an upper limit of the water production starting from H2 and O. These reactants are brought together in a matrix of CO2 in a series of (control) experiments at different temperatures and with different isotopological compositions. The amount of water detected with the quadrupole mass spectrometer upon TPD is found to originate mainly from contamination in the chamber itself. However, if water is produced in small quantities on the surface through H2 + O, this can only be explained by a combined classical and tunneled reaction mechanism. An absolutely conservative upper limit for the reaction rate is derived with a microscopic kinetic Monte Carlo model that converts the upper limit into a maximal possible reaction rate. Incorporating this rate into simulations run for astrochemically relevant parameters, shows that the upper limit to the contribution of the reaction H2 + O in OH, and hence water formation, is 11% in dense interstellar clouds. Our combined experimental and theoretical results indicate however, that this contribution is likely to be much lower.