In this work, we present the results of our investigation into the chemistry of Z- and E-cyanomethanimine (HNCHCN), both of which are possible precursors to the nucleobase adenine. Ab initio quantum chemical calculations for a number of reactions with atomic hydrogen were carried out. We find that the reaction H + Z/E-HNCHCN leading both to H-addition as well as H$_2$-abstraction proceed via similar short-range barriers with bimolecular rate coefficients on the order of $\sim10^{-17}$ cm$^{3}$ s$^{-1}$. These results were then incorporated into astrochemical models and used in simulations of the giant molecular cloud G+0.693. The calculated abundances obtained from these models were compared with previous observational data and found to be in good agreement, with a predicted [Z/E] ratio of $\sim3$ - somewhat smaller than the previously derived value of $6.1\pm2.4$. We find that the [Z/E] ratio in our simulations is due mostly to ion-molecule destruction rates driven by the different permanent dipoles of the two conformers. Based on these results, we propose a general rule-of-thumb for estimating the abundances of isomers in interstellar environments.

We investigated the behavior of H$_2$, main constituent of the gas phase in dense clouds, after collision with amorphous solid water (ASW) surfaces, one of the most abundant chemical species of interstellar ices. We developed a general framework to study the adsorption dynamics of light species on interstellar ices. We provide binding energies and their distribution, sticking probabilities for incident energies between 1 meV and 60 meV, and thermal sticking coefficients between 10 and 300 K for surface temperatures from 10 to 110 K. We found that the sticking probability depends strongly on the adsorbate kinetic energy and the surface temperature, but hardly on the angle of incidence. We observed finite sticking probabilities above the thermal desorption temperature. Adsorption and thermal desorption should be considered as separate events with separate time scales. Laboratory results for these species have shown a gap in the trends attributed to the differently employed experimental techniques. Our results complement observations and extend them, increasing the range of gas temperatures under consideration. We plan to employ our method to study a variety of adsorbates, including radical and charged species.

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.

This paper investigates the impact of the increased reaction rate constant due to tunneling effects on planet-forming disks. Our aim is to quantify the astrophysical implications of atom tunneling for simple molecules that are frequently used to infer disk structure information or to define the initial conditions for planet (atmosphere) formation. We explain the tunneling effect on reaction rates by using a scholarly example in comparison to previous UMIST2012 rate constants. In a chemical network with 1299 reactions, we identify all reactions that could show atom tunneling. We devise a simple formulation of reaction rate constants that overestimates tunneling and screen a standard T Tauri disk model for changes in species abundances. For those reactions found to be relevant, we find values of the most recent literature for the rate constants including tunneling and compare the resulting disk chemistry to the standard disk models. The rate constants in the UMIST2012 database in many cases already capture tunneling effects implicitly. A rigorous screening procedure identified three neutral-neutral reactions where atom tunneling could change simple molecule abundances. However, by adopting recent values of the rate constants of these reactions and due to the layered structure of planet-forming disks, the effects are limited to a small region where cold neutral-neutral chemistry dominates. Abundances of water close to the midplane snowline can increase by a factor of two at most compared to previous results with UMIST2012 rates. Observables from the disk surface, such as high excitation (> 500 K) water line fluxes, decrease by 60% at most when tunneling effects are explicitly excluded. On the other hand, disk midplane quantities relevant for planet formation such as the C-to-O ratio and also the ice-to-rock ratio are clearly affected by these gas-phase tunneling effects.

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.

To date, two isomers of H$_2$C$_3$O have been detected, namely, propynal (HCCCHO) and cylclopropenone (c-H$_2$C$_3$O). A third, propadienone (CH$_2$CCO), has thus far eluded observers despite the fact that it is the lowest in energy of the three. This previously noted result is in contradiction of the minimum energy principle, which posits that the abundances of isomers in interstellar environments can be predicted based on their relative stabilities - and suggests, rather, the importance of kinetic over thermodynamic effects in explaining the role of such species. Here, we report results of \textit{ab initio} quantum chemical calculations of the reaction between H and (a) HC$_3$O, (b) H$_2$C$_3$O (both propynal and propadienone), and (c) CH$_2$CHCO. We have found that, among all possible reactions between atomic hydrogen and either propadienone or propynal, only the destruction of propadienone is barrierless and exothermic. That this destruction pathway is indeed behind the non-detection of CH$_2$CCO is further suggested by our finding that the product of this process, the radical CH$_2$CHCO, can subsequently react barrierlessly with H to form propenal (CH$_2$CHCHO) which has, in fact, been detected in regions where the other two H$_2$C$_3$O isomers are observed. Thus, these results not only shed light on a previously unresolved astrochemical mystery, but also further highlight the importance of kinetics in understanding the abundances of interstellar molecules.

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.

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.

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.

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.

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.

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.

Investigating how formamide forms in the interstellar medium is a hot topic in astrochemistry, which can contribute to our understanding of the origin of life on Earth. We have constructed a QM/MM model to simulate the hydrogenation of isocyanic acid on amorphous solid water surfaces to form formamide. The binding energy of HNCO on the ASW surface varies significantly between different binding sites, we found values between $\sim$0 and 100 kJ mol$^{-1}$. The barrier for the hydrogenation reaction is almost independent of the binding energy, though. We calculated tunneling rate constants of H + HNCO $\rightarrow$ NH$_2$CO at temperatures down to 103 K combining QM/MM with instanton theory. Tunneling dominates the reaction at such low temperatures. The tunneling reaction is hardly accelerated by the amorphous solid water surface compared to the gas phase for this system, even though the activation energy of the surface reaction is lower than the one of the gas-phase reaction. Both the height and width of the barrier affect the tunneling rate in practice. Strong kinetic isotope effects were observed by comparing to rate constants of D + HNCO $\rightarrow$ NHDCO. At 103 K we found a KIE of 231 on the surface and 146 in the gas phase. Furthermore, we investigated the gas-phase reaction NH$_2$ + H$_2$CO $\rightarrow$ NH$_2$CHO + H and found it unlikely to occur at cryogenic temperatures. The data of our tunneling rate constants are expected to significantly influence astrochemical models.

We calculated reaction rate constants including atom tunneling of the reaction of dihydrogen with the hydroxy radical down to a temperature of 50 K. Instanton theory and canonical variational theory with microcanonical optimized multidimensional tunneling (CVT/$\mu$OMT) were applied using a fitted potential energy surface [J. Chem. Phys. 138, 154301 (2013)]. All possible protium/deuterium isotopologues were considered. Atom tunneling increases at about 250 K (200 K for deuterium transfer). Even at 50 K the rate constants of all isotopologues remain in the interval $ 4 \cdot 10^{-20}$ to $4 \cdot 10^{-17}$ cm$^3$ s$^{-1}$ , demonstrating that even deuterated versions of the title reaction are possibly relevant to astrochemical processes in molecular clouds. The transferred hydrogen atom dominates the kinetic isotope effect at all temperatures.