-
Pathways to Interstellar Amides via Carbamoyl (NH2CO) Isomers by Radical-Neutral Reactions on Ice Grain Mantles
Authors:
Gabriela Silva-Vera,
Giulia M. Bovolenta,
Namrata Rani,
Sebastian Vera,
Stefan Vogt-Geisse
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
In-Depth Exploration of Catalytic Sites on Amorphous Solid Water: I. The Astrosynthesis of Aminomethanol
Authors:
Giulia M. Bovolenta,
Gabriela Silva-Vera,
Stefano Bovino,
German Molpeceres,
Johannes Kästner,
Stefan Vogt-Geisse
Abstract:
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 potent…
▽ More
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.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
The Astrochemistry Low-energy Electron Cross-Section (ALeCS) database I. Semi-empirical electron-impact ionization cross-section calculations and ionization rates
Authors:
Brandt A. L. Gaches,
Tommaso Grassi,
Stefan Vogt-Geisse,
Giulia M. Bovolenta,
Claire Vallance,
David Heathcote,
Marco Padovani,
Stefano Bovino,
Prasanta Gorai
Abstract:
(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-ener…
▽ More
(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.
△ Less
Submitted 13 February, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
Binding Energy Evaluation Platform: A database of quantum chemical binding energy distributions for the astrochemical community
Authors:
Giulia M. Bovolenta,
Stefan Vogt-Geisse,
Stefano Bovino,
Tommaso Grassi
Abstract:
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)…
▽ More
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.
△ Less
Submitted 26 July, 2022;
originally announced July 2022.
-
High level ab initio binding energy distribution of molecules on interstellar ices: Hydrogen fluoride
Authors:
Giulia Bovolenta,
Stefano Bovino,
Esteban Vöhringer-Martinez,
David A. Saez,
Tommaso Grassi,
Stefan Vogt-Geisse
Abstract:
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 me…
▽ More
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.
△ Less
Submitted 18 October, 2020;
originally announced October 2020.
-
A novel framework to study the impact of binding energy distributions on the chemistry of dust grains
Authors:
T. Grassi,
S. Bovino,
P. Caselli,
G. Bovolenta,
S. Vogt-Geisse,
B. Ercolano
Abstract:
The evaporation of molecules from dust grains is crucial to understand some key aspects of the star- and the planet-formation processes. During the warm-up phase the presence of young protostellar objects induces molecules to evaporate from the dust surface into the gas phase, enhancing its chemical complexity. Similarly, in circumstellar disks, the position of the so-called snow-lines is determin…
▽ More
The evaporation of molecules from dust grains is crucial to understand some key aspects of the star- and the planet-formation processes. During the warm-up phase the presence of young protostellar objects induces molecules to evaporate from the dust surface into the gas phase, enhancing its chemical complexity. Similarly, in circumstellar disks, the position of the so-called snow-lines is determined by evaporation, with important consequences for the formation of planets. The amount of molecules that are desorbed depends on the interaction between the species and the grain surface, which is controlled by the binding energy. Recent theoretical and experimental works point towards a distribution of values for this parameter instead of the single value often employed in astrochemical models.We present here a new "multi-binding energy" framework, to assess the effects that a distribution of binding energies has on the amount of species bound to the grains. We find that the efficiency of the surface chemistry is significantly influenced by this process with crucial consequences on the theoretical estimates of the desorbed species.
△ Less
Submitted 8 September, 2020;
originally announced September 2020.
-
Determination of pulsation periods and other parameters of 2875 stars classified as MIRA in the All Sky Automated Survey (ASAS)
Authors:
N. Vogt,
A. Contreras-Quijada,
I. Fuentes-Morales,
S. Vogt-Geisse,
C. Arcos,
C. Abarca,
C. Agurto-Gangas,
M. Caviedes,
H. DaSilva,
J. Flores,
V. Gotta,
F. Peñaloza,
K. Rojas,
J. I. Villaseñor
Abstract:
We have developed an interactive PYTHON code and derived crucial ephemeris data of 99.4% of all stars classified as 'Mira' in the ASAS data base, referring to pulsation periods, mean maximum magnitudes and, whenever possible, the amplitudes among others. We present a statistical comparison between our results and those given by the AAVSO International Variable Star Index (VSX), as well as those de…
▽ More
We have developed an interactive PYTHON code and derived crucial ephemeris data of 99.4% of all stars classified as 'Mira' in the ASAS data base, referring to pulsation periods, mean maximum magnitudes and, whenever possible, the amplitudes among others. We present a statistical comparison between our results and those given by the AAVSO International Variable Star Index (VSX), as well as those determined with the machine learning automatic procedure of Richards et al. 2012. Our periods are in good agreement with those of the VSX in more than 95% of the stars. However, when comparing our periods with those of Richards et al, the coincidence rate is only 76% and most of the remaining cases refer to aliases. We conclude that automatic codes require still more refinements in order to provide reliable period values. Period distributions of the target stars show three local maxima around 215, 275 and 330 d, apparently of universal validity, their relative strength seems to depend on galactic longitude. Our visual amplitude distribution turns out to be bimodal, however 1/3 of the targets have rather small amplitudes (A $<$ 2.5$^{m}$) and could refer to semi-regular variables (SR). We estimate that about 20% of our targets belong to the SR class. We also provide a list of 63 candidates for period variations and a sample of 35 multiperiodic stars which seem to confirm the universal validity of typical sequences in the double period and in the Petersen diagrams
△ Less
Submitted 16 September, 2016;
originally announced September 2016.