In the interstellar medium, six molecules have been conclusively detected in the solid state in interstellar ices, and a few dozen have been hypothesized and modeled to be present in the solid state ...as well. The icy mantles covering micrometer-sized dust grains are, in fact, thought to be at the core of complex molecule formation as a consequence of the local high density of molecules that are simultaneously adsorbed. From a structural perspective, the icy mantle is considered to be layered, with an amorphous water-rich inner layer surrounding the dust grain, covered by an amorphous CO-rich outer layer. Moreover, recent studies have suggested that the CO-rich layer might be crystalline and possibly even be segregated as a single crystal atop the ice mantle. If so, there are far-reaching consequences for the formation of more complex organic molecules, such as methanol and sugars, that use CO as a backbone. Validation of these claims requires further investigation, in particular on acquiring atomistic insight into surface processes, such as adsorption, diffusion, and reactivity on CO ices. Here, we present the first detailed computational study toward treating the weak interaction of (pure) CO ices. We provide a benchmark of the performance of various density functional theory methods in treating the binding of pure CO ices. Furthermore, we perform an atomistic and in-depth study of the binding energy of CO on amorphous and crystalline CO ices using a pair-potential-based force field. We find that CO adsorption is represented by a large distribution of binding energies (200–1600 K) on amorphous CO, including a significant amount of weak binding sites (<350 K). Increasing both the cluster size and the number of neighbors increases the mean of the observed binding energy distribution. Finally, we find that CO binding energies are dominated by dispersion and, as such, exchange-correlation functionals need to include a treatment of dispersion to accurately simulate surface processes on CO ices. In particular, we find the ωB97M-V functional to be a strong candidate for such simulations.
Icy dust grains in space act as catalytic surfaces onto which complex molecules form. These molecules are synthesized through exothermic reactions from precursor radicals and, mostly, hydrogen atom ...additions. Among the resulting products are species of biological relevance, such as hydroxylamine-NH sub(2)OH-a precursor molecule in the formation of amino acids. In this Letter, laboratory experiments are described that demonstrate NH sub(2)OH formation in interstellar ice analogs for astronomically relevant temperatures via successive hydrogenation reactions of solid nitric oxide (NO). Inclusion of the experimental results in an astrochemical gas-grain model proves the importance of a solid-state NO + H reaction channel as a starting point for prebiotic species in dark interstellar clouds and adds a new perspective to the way molecules of biological importance may form in space.
Icy dust grains in space act as catalytic surfaces onto which complex molecules form. These molecules are synthesized through exothermic reactions from precursor radicals and, mostly, hydrogen atom ...additions. Among the resulting products are species of biological relevance, such as hydroxylamine-NH{sub 2}OH-a precursor molecule in the formation of amino acids. In this Letter, laboratory experiments are described that demonstrate NH{sub 2}OH formation in interstellar ice analogs for astronomically relevant temperatures via successive hydrogenation reactions of solid nitric oxide (NO). Inclusion of the experimental results in an astrochemical gas-grain model proves the importance of a solid-state NO+H reaction channel as a starting point for prebiotic species in dark interstellar clouds and adds a new perspective to the way molecules of biological importance may form in space.
The impact of including the reactions of C and CH with molecular hydrogen in a gas-grain network is assessed via a sensitivity analysis. To this end, we vary 3 parameters, namely, the efficiency for ...the reaction \ce{C + H2 -> CH2}, and the cosmic ray ionisation rate, with the third parameter being the final density of the collapsing dark cloud. A grid of 12 models is run to investigate the effect of all parameters on the final molecular abundances of the chemical network. We find that including reactions with molecular hydrogen alters the hydrogen economy of the network; since some species are hydrogenated by molecular hydrogen, atomic hydrogen is freed up. The abundances of simple molecules produced from hydrogenation, such as \ce{CH4}, \ce{CH3OH} and \ce{NH3}, increase, and at the same time, more complex species such as glycine and its precursors see a significant decrease in their final abundances. We find that the precursors of glycine are being preferentially hydrogenated, and therefore glycine itself is produced less efficiently.
The chemical network governing interstellar sulfur has been the topic of unrelenting discussion for the past decades due to the conspicuous discrepancy between its expected and observed abundances in ...different interstellar environments. More recently, the astronomical detections of CH3CH2SH and CH2CS highlighted the importance of interstellar formation routes for sulfur-bearing organic molecules with two carbon atoms. In this work, we perform a laboratory investigation of the solid-state chemistry resulting from the interaction between C2H2 molecules and SH radicals -- both thought to be present in interstellar icy mantles -- at 10 K. Reflection absorption infrared spectroscopy and quadrupole mass spectrometry combined with temperature-programmed desorption experiments are employed as analytical techniques. We confirm that SH radicals can kick-start a sulfur reaction network under interstellar cloud conditions and identify at least six sulfurated products: CH3CH2SH, CH2CHSH, HSCH2CH2SH, H2S2, and tentatively CH3CHS and CH2CS. Complementarily, we utilize computational calculations to pinpoint the reaction routes that play a role in the chemical network behind our experimental results. The main sulfur-bearing organic molecule formed under our experimental conditions is CH3CH2SH and its formation yield increases with the ratios of H to other reactants. It serves as a sink to the sulfur budget within the network, being formed at the expense of the other unsaturated products. The astrophysical implications of the chemical network proposed here are discussed.
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.