Search

Explaining the formation pathways of amides on ice-grain mantels is crucial to understanding the prebiotic chemistry in an interstellar medium. In this computational study, we explore different radical-neutral formation pathways for some of the observed amides (formamide, acetamide, urea, and N-methylformamide) via intermediate carbamoyl (NH2CO) radical precursors and their isomers. We assess the relative energy of four NH2CO isomers in the gas phase and evaluate their binding energy on small water clusters to discern the affinity that the isomers present to an ice model. We consider three possible reaction pathways for the formation of the carbamoyl radicals, namely, the OH + HCN, CN + H2O, and NH2 + CO reaction channels. We computed the binding energy distribution for the HCN and CH3CN precursors on an ice model consisting of a set of clusters of 22 water molecules each to serve as a starting point for the reactivity study on the ice surface. The computations revealed that the lowest barrier to the formation of an NH2CO isomer corresponds to the NH2 + CO reaction (12.6 kJ/mol). The OH + HCN reaction pathway results in the exothermic formation of the N-radical form of carbamoyl HN(C=O)H with a reaction barrier of 26.7 kJ/mol. We found that the CN + H2O reaction displays a high energy barrier of 70.6 kJ/mol. Finally, we also probed the direct formation of the acetamide radical precursor via the OH + CH3CN reaction and found that the most probable outcome on interstellar ices is the H-abstraction reaction to yield CH2CN and H2O. Based on these results, we believe that including alternative reaction pathways, leading to the formation of amides via the N-radical form of carbamoyl, would provide an improvement in the prediction of the amide abundances in astrochemical models, especially regarding the chemistry of star-forming regions., Comment: Journal is ACS Earth Space Chem, 2024. ASAP (As Soon As Publishable)

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.

(Abridged) Electron-molecule interaction is a fundamental process in radiation-driven chemistry in space, from the interstellar medium to comets. Therefore, knowledge of interaction cross-sections is key. While there has been a plethora of studies of total ionization cross-sections, data is often spread over many sources, or not public or readily available. We introduce the Astrochemistry Low-energy Electron Cross-Section (ALeCS) database, a public database for electron interaction cross-sections and ionization rates for molecules of astrochemical interest. In this work, we present the first data release comprising total ionization cross-sections and ionization rates for over 200 neutral molecules. We include optimized geometries and molecular orbital energies at various levels of theory, and for a subset of the molecules, the ionization potentials. We compute total ionization cross-sections using the binary-encounter Bethe model and screening-corrected additivity rule, and ionization rates and reaction network coefficients for molecular cloud environments for $>$200 neutral molecules ranging from diatomics to complex organics. We demonstrate that our binary-encounter Bethe cross-sections agree well with experimental data. We show that the ionization rates scale roughly linearly with the number of constituent atoms in the molecule. We introduce and describe the public ALeCS database. For the initial release, we include total ionization cross-sections for $>$200 neutral molecules and several cations and anions calculated with different levels of quantum chemistry theory, the chemical reaction rates for the ionization, and network files in the formats of the two most popular astrochemical networks, the KIDA and UMIST. The database will be continuously updated for more molecules and interactions., Comment: Accepted to A&A. Public database online at https://github.com/AstroBrandt/ALeCS

The quality of astrochemical models is highly dependent on reliable binding energy (BE) values that consider the morphological and energetic variety of binding sites on the surface of ice-grain mantles. Here, we present the Binding Energy Evaluation Platform (BEEP) and database that, using quantum chemical methods, produces full BE distributions of molecules bound to an amorphous solid water (ASW) surface model. BEEP is highly automatized and allows to sample binding sites on set of water clusters and to compute accurate BEs. Using our protocol, we computed 21 BE distributions of interstellar molecules and radicals on an amorphized set of 15-18 water clusters of 22 molecules each. The distributions contain between 225 and 250 unique binding sites. We apply a Gaussian fit and report the mean and standard deviation for each distribution. We compare with existing experimental results and find that the low and high coverage experimental BEs coincide well with the high BE tail and mean value of our distributions, respectively. Previously reported single BE theoretical values are broadly in line with ours, even though in some cases significant differences can be appreciated. We show how the use of different BE values impact a typical problem in astrophysics, such as the computation of snow lines in protoplanetary discs. BEEP will be publicly released so that the database can be expanded to other molecules or ice-models in a community effort., Comment: Accepted on ApJS

The knowledge of the binding energy of molecules on astrophysically relevant ices can help to obtain an estimate of the desorption rate, i.e. the molecules residence time on the surface. This represents an important parameter for astrochemical models, crucial to determine the chemical fate of complex organic molecules formed on dust grains and observed in the densest regions of the interstellar medium. In this work, we propose a new robust procedure to study the interaction of atoms and molecules with interstellar ices, based on \textit{ab initio} molecular dynamics and density functional theory, validated by high-level \textit{ab initio} methods at a CCSD(T)/CBS level. We have applied this procedure to hydrogen fluoride (HF), a promising tracer of the molecular content of galaxies. In total we found 13 unique equilibrium structures of HF binding to small water clusters of up to 4 molecules, with binding energies ranging from 1208 to 7162 K. We computed a 22-molecules model of amorphous solid water (ASW) surface using \textit{ab initio} molecular dynamics simulations and carried out a systematic analysis of the binding sites of HF, in terms of binding modes and binding energies. Considering 10 different water clusters, we found a binding energy distribution with an average value of $5313\pm74$ K, and a dispersion of $921\pm115$ K. Finally, the effect of the electrostatic field of the 22 water molecules on the binding energies was investigated incrementally by symmetry adapted perturbation theory, in order to gauge the effect of the water environment. The results indicate that the extent of the electrostatic interaction of HF with ASW depends strongly on the properties of the binding site. We expect that this work will provide a solid foundation for a systematic development of a binding energy distribution database of molecules on interstellar surfaces., Comment: 12 pages, 9 figures, 6 tables, accepted to Mol. Astrophys