Search

In recent years, the continued detection of complex organic molecules of prebiotic interest has refueled the interest on a panspermic origin of life. The prebiotic molecule glyceraldehyde is proposed to be formed from (Z)-1,2-ethenediol, a molecule recently detected towards the G+0.693-0.027 molecular cloud at the galactic center. In this work, we computationally simulate the formation of (Z)-1,2-ethenediol from vinyl alcohol on the surface of amorphous solid water in a two-step synthesis involving a OH addition and a H abstraction reaction. In total, we considered all reaction possibilities of the 1,1 and 1,2-OH addition to vinyl alcohol followed by H-abstraction or H-addition reactions on the resulting radicals. The combination of these reactions is capable of explaining the formation of (Z)-1,2-ethenediol provided a suprathermal diffusion of OH. We also conclude that our proposed formation pathway is not selective and also yields other abstraction and addition products. Key in our findings is the connection between the adsorption modes of the reactants and intermediates and the stereoselectivity of the reactions., Comment: This version is the unedited accepted version in the Astrophysical Journal

Methane (CH4) is one of the major components of the icy mantle of cosmic dust prevalent in cold, dense regions of interstellar media, playing an important role in the synthesis of complex organic molecules and prebiotic molecules. Solid CH4 is considered to be formed via the successive hydrogenation of C atoms accreting onto dust: C + 4H -> CH4. However, most astrochemical models assume this reaction on the ice mantles of dust to be barrierless and efficient, without considering the states of adsorption. Recently, we found that C atoms exist in either the physisorbed or chemisorbed state on compact amorphous solid water, which is analogous to an interstellar ice mantle. These distinct adsorption states considerably affect the hydrogenation reactivity of the C atom. Herein, we elucidate the reactivities of physisorbed and chemisorbed C atoms with H atoms via sequential deposition and co-deposition processes. The results indicate that only physisorbed C atoms can produce CH4 on ice. Combining this finding with a previous estimate for the fraction of physisorbed C atoms on ice, we determined the upper limit for the conversion of C atoms into CH4 to be 30%., Comment: 24 pages, 4 figures, 1 table; accepted for publication in ApJ

Chemical processes taking place on ice-grain mantles are pivotal to the complex chemistry of interstellar environments. In this study, we conducted a comprehensive analysis of the catalytic effects of an amorphous solid water (ASW) surface on the reaction between ammonia (NH$_3$) and formaldehyde (H$_2$CO) to form aminomethanol (NH$_2$CH$_2$OH) using density functional theory. We identified potential catalytic sites based on the binding energy distribution of NH$_3$ and H$_2$CO reactants, on a set-of-clusters surface model composed of 22 water molecules and found a total of 14 reaction paths. Our results indicate that the catalytic sites can be categorized into four groups, depending on the interactions of the carbonyl oxygen and the amino group with the ice surface in the reactant complex. A detailed analysis of the reaction mechanism using Intrinsic Reaction Coordinate and reaction force analysis revealed three distinct chemical events for this reaction: formation of the C--N bond, breaking of the N--H bond, and formation of the O--H hydroxyl bond. Depending on the type of catalytic site, these events can occur within a single, concerted, albeit asynchronous, step, or can be isolated in a step-wise mechanism, with the lowest overall transition state energy observed at 1.3 kcal mol$^{-1}$. A key requirement for the low-energy mechanism is the presence of a pair of dangling OH bonds on the surface, found at 5\% of the potential catalytic sites on an ASW porous surface.

We investigate the role of carbon monoxide ice in the chemical evolution of prestellar cores using astrochemical rate equation models. We constrain the ratios of the binding energies on CO ice and H$_{2}$O ice for a series of adsorbates deemed important in diffusive chemistry on H$_{2}$O ices. We later include these ratios in our chemical reaction network model, where the binding and diffusion energies of icy species vary as a function of the surface composition. When the surface coverage of CO increases, the model shows an enhancement of O-bearing complex organic molecules, especially those formed from the intermediate products of CO hydrogenation (e.g. HCO) and CH$_{3}$/CH$_{2}$. Because the binding energy of CH$_{3}$/CH$_{2}$ is in the right range, its diffusion rate increases significantly with CO coverage. At $T>$14 K and with less influence, enhanced diffusion of HCO also contributes to the increase of the abundances of COM. We find, however, that chemistry is not always enhanced on CO ice and that the temperature and cosmic ray ionization rate of each astronomical object is crucial for this particular chemistry, revealing a highly non-trivial behavior that needs to be addressed on a per-case basis. Our results are highly relevant in the context of interstellar ice observations with JWST., Comment: Accepted for publication in Astronomy and Astrophysics

We report the first detection in the interstellar medium of $N$-cyanomethanimine (H$_2$CNCN), the stable dimer of HCN of highest energy, and the most complex organic molecule identified in space containing the prebiotically relevant NCN backbone. We have identified a plethora of $a$-type rotational transitions with 3 $\leq J_\text{up} \leq$ 11 and $K_\text{a} \leq$ 2 that belong to this species towards the Galactic Center G+0.693-0.027 molecular cloud, the only interstellar source showing the three cyanomethanimine isomers (including the $Z$- and $E$- isomers of $C$-cyanomethanimine, HNCHCN). We have derived a total column density for H$_2$CNCN of (2.9$\, \pm \,$0.1)$\times$10$^{12}$ cm$^{-2}$, which translates into a total molecular abundance with respect to H$_2$ of (2.1$\, \pm \,$0.3)$\times$10$^{-11}$. We have also revisited the previous detection of $E$- and $Z$-HNCHCN, and found a total $C/N$-cyanomethanimine abundance ratio of 31.8$\, \pm \,$1.8 and a $Z/E$-HNCHCN ratio of 4.5$\, \pm \,$0.2. While the latter can be explained on the basis of thermodynamic equilibrium, chemical kinetics are more likely responsible for the observed $C/N$-cyanomethanimine abundance ratio, where the gas-phase reaction between methanimine (CH$_2$NH) and the cyanogen radical (CN) arises as the primary formation route., Comment: 22 pages, 8 figures, 3 tables, 2 appendix - Accepted for publication in The Astrophysical Journal

Many interstellar complex organic molecules (COMs) are believed to be produced on the surfaces of icy grains at low temperatures. Atomic carbon is considered responsible for the skeletal evolution processes, such as C-C bond formation, via insertion or addition reactions. Before reactions, C atoms must diffuse on the surface to encounter reaction partners; therefore, information on their diffusion process is critically important for evaluating the role of C atoms in the formation of COMs. In situ detection of C atoms on ice was achieved by a combination of photostimulated desorption and resonance enhanced multiphoton ionization methods. We found that C atoms weakly bound to the ice surface diffused approximately above 30 K and produced C2 molecules. The activation energy for C-atom surface diffusion was experimentally determined to be 88 meV (1,020 K), indicating that the diffusive reaction of C atoms is activated at approximately 22 K on interstellar ice. The facile diffusion of C at T > 22 K atoms on interstellar ice opens a previously overlooked chemical regime where the increase in complexity of COMs as driven by C atoms. Carbon addition chemistry can be an alternative source of chemical complexity in translucent clouds and protoplanetary disks with crucial implications in our current understanding on the origin and evolution of organic chemistry in our Universe., Comment: 33 pages (main + SI), 14 figures, 1 table

CO2 is one of the dominant components of the interstellar ice. Recent observations show CO2 exists more abundantly in polar (H2O-dominated) ice than in apolar (H2O-poor) ice. CO2 ice formation is primarily attributed to the reaction between CO and OH, which has a barrier. Highly accurate quantum chemical calculations were employed to analyze the stationary points of the potential energy surfaces of the title reaction in the gas phase on a H2O and CO clusters. Microcanonical transition state theory was used as a diagnostic tool for the efficiency of the reaction under ISM conditions. We simulate the kinetics of ice chemistry, considering different scenarios involving non-thermal processes and energy dissipation. The CO + OH reaction proceeds through the remarkably stable intermediate HOCO radical. On the H2O cluster, the formation of this intermediate is efficient, but the subsequent reaction leading to CO2 formation is not. Conversely, HOCO formation on the CO cluster is inefficient without external energy input. Thus, CO2 ice cannot be formed by the title reaction alone either on the H2O cluster or CO cluster. In the polar ice, CO2 ice formation is possible via CO + OH -> HOCO, followed by HOCO + H ->CO2 + H2, as demonstrated by abundant experimental literature. In apolar ice, CO2 formation is less efficient because HOCO formation requires external energy. Our finding is consistent with the JWST observations. Further experimental work is encouraged using low-temperature OH radicals., Comment: Accepted for publication in Astronomy and Astrophysics

Energy redistribution after a chemical reaction is one of the few mechanisms to explain the diffusion and desorption of molecules which require more energy than the thermal energy available in quiescent molecular clouds (10 K). This energy distribution can be important in phosphorous hydrides, elusive yet fundamental molecules for interstellar prebiotic chemistry. We studied the reaction dynamics of the \ce{P + H -> PH} reaction on amorphous solid water, a reaction of astrophysical interest, using \emph{ab-initio} molecular dynamics with atomic forces evaluated by a neural network interatomic potential. We found that the exact nature of the initial phosphorous binding sites is less relevant for the energy dissipation process because the nascent PH molecule rapidly migrates to sites with higher binding energy after the reaction. Non-thermal diffusion and desorption-after-reaction were observed and occurred early in the dynamics, essentially decoupled from the dissipation of the chemical reaction energy. From an extensive sampling of reactions on sites, we constrained the average dissipated reaction energy within the simulation time (50 ps) to be between 50 and 70 %. Most importantly, the fraction of translational energy acquired by the formed molecule was found to be mostly between 1 and 5 %. Including these values, specifically for the test cases of 2% and 5% of translational energy conversion, in astrochemical models, reveals very low gas-phase abundances of PH$_{x}$ molecules and reflects that considering binding energy distributions is paramount for correctly merging microscopic and macroscopic modelling of non-thermal surface astrochemical processes. Finally, we found that PD molecules dissipate more of the reaction energy. This effect can be relevant for the deuterium fractionation and preferential distillation of molecules in the interstellar medium., Comment: Accepted for publication in Astronomy and Astrophysics

The abundance ratios between isomers of a COM observed in the ISM provides valuable information about the chemistry and physics of the gas and eventually, the history of molecular clouds. In this context, the origin of an abundance of c-HCOOH acid of only 6% the isomer c-HCOOH abundance in cold cores, remains unknown. Herein, we explain the presence of c-HCOOH in dark molecular clouds through the destruction and back formation of c-HCOOH and t-HCOOH in a cyclic process that involves HCOOH and highly abundant molecules such as HCO+ and NH3. We use high-level ab initio methods to compute the potential energy profiles for the cyclic destruction/formation routes of c-HCOOH and t-HCOOH. Global rate constants and branching ratios were calculated based on the transition state theory and the master equation formalism under the typical conditions of the ISM. The destruction of HCOOH by reaction with HCO+ in the gas phase leads to three isomers of the cation HC(OH)2+. The most abundant cation can react in a second step with other abundant molecules of the ISM like NH3 to form back c-HCOOH and t-HCOOH. This mechanism explains the formation of c-HCOOH in dark molecular clouds. Considering this mechanism, the fraction of c-HCOOH with respect t-HCOOH is 25.7%. To explain the 6% reported by the observations we propose that further destruction mechanisms of the cations of HCOOH should be taken into account. The sequential acid-base (SAB) mechanism proposed in this work involves fast processes with very abundant molecules in the ISM. Thus, HCOOH very likely suffers our proposed transformations in the conditions of dark molecular clouds. This is a new approach in the framework of the isomerism of organic molecules in the ISM which has the potential to try to explain the ratio between isomers of organic molecules detected in the ISM.

Recent interstellar detections include a significant number of molecules containing vinyl (C2H3) and ethyl (C2H5) groups in their structure. For several of these molecules, there is not a clear experimental or theoretical evidence that support their formation from simpler precursors. We carried out a systematic search of viable reactions starting from closed shell hydrocarbons containing two carbon atoms (ethane, C2H6; ethylene, C2H4; and acetylene, C2H2) with the goal of determining viable chemical routes for the formation of vinyl and ethyl molecules on top of interstellar dust grains. Our results show that both H and OH radicals are key in converting acetylene and ethylene into more complex radicals that are susceptible to continue reacting and forming interstellar complex organic molecules. The relevant reactions, for example OH additions, present rate constants above 10$^{1}$ s$^{-1}$ that are likely competitive with OH diffusion on grains. Similarly, H atom addition to acetylene and ethylene is a very fast process with rate constants above 10$^{4}$ s$^{-1}$ in all cases, and greatly enhanced by quantum tunneling. Hydrogen abstraction reactions are less relevant, but may play a role in specific cases involving the OH radical. Reactions with other radicals NH2, CH3 are likely to have a much lesser impact in the chemistry of ethyl and vinyl bearing molecules., Comment: Accepted for publication in Astronomy and Astrophysics. Supporting structures available (see text). Numerical values of the rate constants are provided as supplementary data

The isomerism of molecules in the interstellar medium and the mechanisms behind it are essential questions in the chemistry of organic molecules in space. In particular, for the simple formic and thioformic acids, the low temperatures found in molecular clouds indicate that cis-trans isomerization in the gas-phase must be impeded. Reactions happening on top of interstellar dust grains may explain the isomer interconversion at low temperatures. We studied the isomerization processes of formic and thioformic acid susceptible to happen on the surface of interstellar dust grains and initiated by H abstraction reactions. Similarly, deuterium enrichment of the acids can occur by the same mechanism. Our objective is to shed light on both topics to increase our understanding of key precursors of organic molecules in space., Comment: Accepted for publication in Astronomy and Astrophysics

The diffusion of atoms and radicals on interstellar dust grains is a fundamental ingredient for predicting accurate molecular abundances in astronomical environments. Quantitative values of diffusivity and diffusion barriers usually rely heavily on empirical rules. In this paper, we compute the diffusion coefficients of adsorbed nitrogen atoms by combining machine-learned interatomic potentials, metadynamics, and kinetic Monte Carlo simulations. With this approach, we obtain a diffusion coefficient of nitrogen atoms on the surface of amorphous solid water of merely $(3.5 \pm 1.1)10^{-34}$cm$^2$s$^{-1}$ at 10 K for a bare ice surface. Thus, we find that nitrogen, as a paradigmatic case for light and weakly bound adsorbates, is unable to diffuse on bare amorphous solid water at 10 K. Surface coverage has a strong effect on the diffusion coefficient by modulating its value over 9--12 orders of magnitude at 10 K and enables diffusion for specific conditions. In addition, we have found that atom tunneling has a negligible effect. Average diffusion barriers of the potential energy surface (2.56 kJ mol$^{-1}$) differ strongly from the effective diffusion barrier obtained from the diffusion coefficient for a bare surface (6.06 kJ mol$^{-1}$) and are, thus, inappropriate for diffusion modeling. Our findings suggest that the thermal diffusion of N on water ice is a process that is highly dependent on the physical conditions of the ice., Comment: Accepted in MNRAS

We report new computational and experimental evidence of an efficient and astrochemically relevant formation route to formaldehyde (H$_2$CO). This simplest carbonylic compound is central to the formation of complex organics in cold interstellar clouds, and is generally regarded to be formed by the hydrogenation of solid-state carbon monoxide. We demonstrate H$_2$CO formation via the reaction of carbon atoms with amorphous solid water. Crucial to our proposed mechanism is a concerted proton transfer catalyzed by the water hydrogen bonding network. Consequently, the reactions $^3$C + H$_2$O -> $^3$HCOH and $^1$HCOH -> $^1$H$_2$CO can take place with low or without barriers, contrary to the high-barrier traditional internal hydrogen migration. These low barriers or absence thereof explain the very small kinetic isotope effect in our experiments when comparing the formation of H$_2$CO to D$_2$CO. Our results reconcile the disagreement found in the literature on the reaction route: C + H$_2$O -> H$_2$CO., Comment: Accepted for publication in JPCL

With the presence of evermore complex S-bearing molecules being detected lately, studies on their chemical formation routes need to keep up the pace to rationalize observations, suggest new candidates for detection, and provide input for chemical evolution models. In this paper, we theoretically characterize the hydrogenation channels of OCS on top of amorphous solid water as an interstellar dust grain analog in molecular clouds. Our results show that the significant reaction outcome is trans-HC(O)SH, a recently detected prebiotic molecule toward G+0.693. The reaction is diastereoselective, explaining the seemingly absence of the cis isomer in the astronomical observations. We found that the reaction proceeds through a highly localized radical intermediate (cis-OCSH), which could be essential in the formation of other sulfur-bearing complex organic molecules due to its slow isomerization dynamics on top of amorphous solid water., Comment: Accepted in ApJ

The recent wave of detections of interstellar aromatic molecules has sparked interest in the chemical behavior of aromatic molecules under astrophysical conditions. In most cases, these detections have been made through chemically related molecules, called proxies, that implicitly indicate the presence of a parent molecule. In this study, we present the results of the theoretical evaluation of the hydrogenation reactions of different aromatic molecules (benzene, pyridine, pyrrole, furan, thiophene, silabenzene, and phosphorine). The viability of these reactions allows us to evaluate the resilience of these molecules to the most important reducing agent in the interstellar medium, the hydrogen atom (H). All significant reactions are exothermic and most of them present activation barriers, which are, in several cases, overcome by quantum tunneling. Instanton reaction rate constants are provided between 50 K and 500 K. For the most efficiently formed radicals, a second hydrogenation step has been studied. We propose that hydrogenated derivatives of furan, pyrrole, and specially 2,3-dihydropyrrole, 2,5-dihydropyrrole, 2,3-dihydrofuran, and 2,5-dihydrofuran are promising candidates for future interstellar detections.

Context. The diffusion of volatile species on amorphous solid water ice affects the chemistry on dust grains in the interstellar medium as well as the trapping of gases enriching planetary atmospheres or present in cometary material. Aims. The aim of the work is to provide diffusion coefficients of CH$_4$ on amorphous solid water (ASW), and to understand how they are affected by the ASW structure. Methods. Ice mixtures of H$_2$O and CH$_4$ were grown in different conditions and the sublimation of CH$_4$ was monitored via infrared spectroscopy or via the mass loss of a cryogenic quartz crystal microbalance. Diffusion coefficients were obtained from the experimental data assuming the systems obey Fick's law of diffusion. Monte Carlo simulations modeled the different amorphous solid water ice structures investigated and were used to reproduce and interpret the experimental results. Results. Diffusion coefficients of methane on amorphous solid water have been measured to be between 10$^{-12}$ and 10$^{-13}$ cm$^2$ s$^{-1}$ for temperatures ranging between 42 K and 60 K. We showed that diffusion can differ by one order of magnitude depending on the morphology of amorphous solid water. The porosity within water ice, and the network created by pore coalescence, enhance the diffusion of species within the pores.The diffusion rates derived experimentally cannot be used in our Monte Carlo simulations to reproduce the measurements. Conclusions. We conclude that Fick's law can be used to describe diffusion at the macroscopic scale, while Monte Carlo simulations describe the microscopic scale where trapping of species in the ices (and their movement) is considered., Comment: 17 pages, 15 figures

Dynamics of adsorption and desorption of (4S)-N on amorphous solid water are analyzed using molecular dynamics simulations. The underlying potential energy surface was provided by machine-learned interatomic potentials. Binding energies confirm the latest available theoretical and experimental results. The nitrogen sticking coefficient is close to unity at dust temperatures of 10 K but decreases at higher temperatures. We estimate a desorption time scale of 1 {\mu}s at 28 K. The estimated time scale allows chemical processes mediated by diffusion to happen before desorption, even at higher temperatures. We found that the energy dissipation process after a sticking event happens on the picosecond timescale at dust temperatures of 10 K, even for high energies of the incoming adsorbate. Our approach allows the simulation of large systems for reasonable time scales at an affordable computational cost and ab-initio accuracy. Moreover, it is generally applicable for the study of adsorption dynamics of interstellar radicals on dust surfaces.

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., Comment: The Astrophysical Journal, accepted: 17 pages, 9 figures, 7 tables