Unveiling the chemical kinetics of aminomethanol (NH2CH2OH): insights into O. H and O2 photo-oxidation reactions and formamide dominance

Aminomethanol is released into the atmosphere through various sources, including biomass burning. In this study, we have expounded the chemical kinetics of aminomethanol in the reaction pathways initiated by the hydroxyl radical ( O˙ H) with the aid of ab initio//density functional theory (DFT) i.e., coupled-cluster theory (CCSD(T))//hybrid-DFT (M06-2X/6-311++G (3df, 3pd). We have explored various possible directions of the O˙ H radical on aminomethanol, as well as the formation of distinct pre-reactive complexes. Our computational findings reveal that the H transfer necessitates activation energies ranging from 4.1 to 6.5 kcal/mol from the –CH2 group, 3.5–6.5 kcal/mol from the –NH2 group and 7–9.3 kcal/mol from the –OH group of three rotational conformers. The H transfer from –CH2, –NH2 and –OH exhibits an estimated total rate constant (k OH) of approximately 1.97 × 10−11 cm3 molecule−1 s−1 at 300 K. The branching fraction analysis indicates a pronounced dominance of C-centered NH2 C˙ HOH radicals with a favorability of 77%, surpassing the N-centered N˙ HCH2OH (20%) and O-centered NH2CH2 O˙ (3%) radicals. Moreover, our investigation delves into the oxidation of the prominently favored carbon-centered NH2 C˙ HOH radical through its interaction with atmospheric oxygen molecules. Intriguingly, our findings reveal that formamide (NH2CHO) emerges as the predominant product in the NH2 C˙ HOH + 3O2 reaction, eclipsing alternative outcomes such as amino formic acid (NH2COOH) and formimidic acid (HN = C(H)-OH). At atmospheric conditions pertinent to the troposphere, the branching fraction value for the formation of formamide is about 99%, coupled with a rate constant of 5.5 × 10−12 cm3 molecule−1 s−1. Finally, we have scrutinized the detrimental impact of formamide on the atmosphere. Interaction of formamide with atmospheric hydroxyl radicals could give rise to the production of potentially perilous compounds such as HNCO. Further, unreacted N˙ HCH2OH radicals may initiate the formation of carcinogenic nitrosamines when reacting with trace N-oxides (namely, NO and NO2). This, in turn, escalates the environmental risk factors.


GRAPHICAL ABSTRACT 1 Introduction
Amines, especially those bearing the -NH 2 functional group (RNH 2 ), find extensive applications in the chemical industry, including uses in solvents, catalysis, surfactants, pharmaceuticals, adhesives, dyes and pigments, etc. (Ge et al., 2011a).For instance, carbon capture and storage (CCS) technology utilizes amine solvents to separate the CO 2 emission from the acid gas treatment of natural gas and fossil fuel power plants (Dai et al., 2012;Yamada, 2021).The rapid development of CCS technology consequently leads to a significant increase in alkanolamines emissions into the atmosphere due to their high vapor pressure (Kapteina et al., 2005).Additionally, the amines are also released into the atmosphere by various other sources such as thorough industrial emission, wastewater treatment, animal husbandry and automotive activities, emissions from ocean organisms, biomass combustion and the degradation of proteins, etc. (Ge et al., 2011b).These amines undergo degradation and participate in conversion reactions both in gas and aqueous phase (clouds, raindrops, fog) in the atmosphere.As a result, they have detrimental effects on air quality, leading to issues like acid rain, urban smog and tropospheric ozone.Also, the amines and alkanolamines have the potential to contribute to greenhouse gas formation, notably nitrous oxide (N 2 O) and the production of highly carcinogenic compounds and nitrosamines ( Schade and Crutzen, 1995;Ge et al., 2011a;Ge et al., 2011b).Hence, elucidating the transformation, mechanisms and kinetics of amines is of potential importance.
The Strecker synthesis, involving a series of chemical reactions, has received great scholarly attention as a potential technique to synthesize chiral α-amino acids using ammonia (NH 3 ), hydrogen cyanide (HCN) and aldehyde (RCHO) (Strecker, 1850;Gröger, 2003; Nájera and Sansano, 2007;Wang et al., 2011).However, aminomethanol (NH 2 CH 2 OH), a pivotal intermediate for the synthesis of the simplest amino acid namely, glycine, has not been observed in the laboratories (Nielsen et al., 1979;Schutte et al., 1993).This is may be due to the decomposition of aminomethanol to methanimine (CH 2 NH) and water (H 2 O).Previous theoretical investigations have indicated that aminomethanol is kinetically stable in the gas phase with a substantial barrier of 230 kJ/mol toward dehydration to methanimine (Feldmann et al., 2005).Our recent study also supports that the formation of methanimine and water from aminomethanol is negligibly small and forbidden by an energy barrier of 234 kJ/mol (Ali, 2019).Bossa et al., have observed aminomethanol at low temperatures through the thermal reaction between ammonia and formaldehyde (Bossa et al., 2009).They have also pointed out that aminomethanol may exist in hot corinos in gas phase.The recent experimental study on aminomethanol in astrophysical-like conditions also validates that aminomethanol could be generated from amines on ice grains in proto-stellar cores or protoplanetary disks (Singh et al., 2022).These studies strongly corroborate the existence of aminomethanol in the atmosphere.
In our recent study (Ali, 2019), the calculated rate constant of CH 2 O + NH 3 has suggested that aminomethanol could potentially form at higher temperatures rather than under atmospheric conditions, as illustrated in Figure 1.Additionally, the significant lifetime (~4 days) of aminomethanol could facilitate the initiation of various chemical reactions with other atmospheric species (Ali, 2019).Nevertheless, the atmospheric chemistry of aminomethanol remains largely unexplored to date, primarily due to the considerable challenges associated with experimental synthesis.Therefore, it is crucial to delve into the reaction mechanisms and kinetics of aminomethanol to evaluate the potential formation of various compounds, including toxicants or carcinogenic byproducts.However, there has been a lack of experimental studies specifically examining the reaction kinetics of the photo-oxidation of aminomethanol.In light of this gap, theoretical approaches, coupled with state-of-the-art computational methodologies and advanced statistical rate theories, are essential for exploring the intricate atmospheric chemistry of aminomethanol, particularly in extreme environmental conditions.
The fate of aminomethanol, in general, is greatly influenced by its gas-phase oxidation by various oxidizing agents.Hydroxyl radicals ( _ OH) are widely present in the atmosphere and the reactions initiated by them are of significant importance among with by various other potential oxidizing agents.It is noteworthy to mention that the higher-order analogue of aminomethanol, namely, monoethanolamine (MEA) (NH 2 CH 2 CH 2 OH), is commercially available and the corresponding gas-phase species are easily formed.Numerous experimental studies have been conducted to characterize MEA and to determine its rate constants (k OH ) with the _ OH radical, along with the resulting products.For example, Murphy and coworkers have evaluated the rate constant, k OH of MEA by reacting with the _ OH radical at room temperature.They determined the rate constant to be 7.02 ± 0.46 × 10 −11 cm 3 molecule −1 s -1 (Borduas et al., 2013).Similarly, Onel et al., reported the k OH of MEA reacting with _ OH as 7.61 ± 0.76 × 10 −11 , while Karl et al., obtained a value of 9.2 ± 1.1 × 10 −11 cm 3 molecule −1 s -1 (Karl et al., 2012;Onel et al., 2012).These studies demonstrate a rapid reaction between MEA and the _ OH radical, resulting in formamide and formaldehyde as major products, along with the potential formation of hazardous HNCO species.These experimental findings suggest that the initial reaction of MEA with the _ OH, followed by MEA radicals reacting with O 2 , plays a crucial role in determining the atmospheric fate of MEA.In a parallel fashion, the abstraction of a hydrogen atom from NH 2 CH 2 OH by the _ OH radical leads to the formation of three different aminomethanol radicals, including a carbon-centered NH 2 _ CHOH radical, a nitrogen-centered _ NHCH 2 OH radical and an oxygen-centered NH 2 CH 2 _ O radical as well as H 2 O, as depicted in Figure 1.It is evident from previous studies on analogous electronic systems that carbon-centered radicals are more prevalent than nitrogenand oxygen-centered radicals (Xie et al., 2014;Franco et al., 2021).These carbon-centered radicals are subsequently targeted by atmospheric oxygen, leading to the formation of _ OH radical initiated products of aminomethanol namely, formamide (NH 2 CHO), formimidic acid (HN = C(H)-OH) and amino formic acid (NH 2 COOH) in the atmosphere.
Among these product species, formamide is simplest and naturally occurring amide molecule.It is a potential precursor for the synthesis of a broad range of organic molecules (Saladino et al., 2012a;Saitta and Saija, 2014), which are vital for life as it poses all the important elements including carbon C), hydrogen H), oxygen O) and nitrogen N) except the heavy elements like sulphur S) and phosphorus P) (Saladino et al., 2005;2007;2012b).Additionally, formamide has an amide functional (-N-C (=O)-) group, which is essential for the formation of the chains of amino acids in order to build up proteins.It promotes the synthesis of four nucleobases namely, adenine, guanine, cytosine and uracil of ribonucleic acid (RNA) (Saladino et al., 2005;Ferus et al., 2015).It also serves as the pre-cursor for carboxylic acids, amino acids and sugars (Botta et al., 2018).These points clearly demonstrating the potential role of formamide in biochemical processes, that are essential for life's sustenance and evolution.
Also, formamide is a highly ubiquitous molecule in the Universe.It was first detected in space dates back to 1971 when Rubin et al., utilized the 140-foot telescope at the National Radio Astronomy Observatory (NRAO) to scrutinize the Sagittarius B2 region (Sgr B2), situated proximate to the Galactic center (Rubin et al., 1971).Their observations distinctly revealed the three hyperfine components (ΔF = 0) of the 2 1,1 -2 1,2 rotational transition of NH 2 CHO, resonating at 4.62 GHz (6.5 cm).Sgr B2, renowned as the most prodigious star formation region within our Milky Way and marked the maiden identification of an interstellar compound harboring carbon C), hydrogen H), oxygen O), and nitrogen N).Subsequent detections of formamide have been made in diverse astronomical locales, including the vicinity of the young stellar entity W33A (Schutte et al., 1999), within cometary bodies such as C/1995 O1 (Bockelée-Morvan et al., 2000) and C/1996 B2 (Lis et al., 1997) and within the broader expanse of the interstellar medium (Solomon, 1973).
In the current investigation, our focus lies in exploring the intricate mechanistic aspects and branching ratios corresponding to the pivotal stages of the _ OH-initiated oxidation reactions of aminomethanol.The primary objectives of this study are described as follows, 1) unveiling the intricate details about the reaction pathways, energetics and kinetics involved in the interaction of aminomethanol with _ OH radicals and identifying the potential aminomethanol-derived radicals, 2) exploring the reaction pathways, energetics and kinetics involved in the further reactions of aminomethanol-derived radicals with molecular oxygen (O 2 ). 3) analyzing the atmospheric implications arising from the exclusive product molecules resulting from the reactions of aminomethanol with OH and O 2 and 4) conducting a comparative analysis of the pre-reactive complexes, transition states, intermediates, other reactive species and rate constants derived from this study with those reported in previous investigations involving the atmospheric reactions of analogous electronic systems to validate our methodology.
2 Theoretical methodology 2.1 Electronic structure calculations All gas-phase ab initio/density functional theory (DFT) calculations were performed using the Gaussian 16 suite of programs (Frisch et al., 2019).The stationary points of all molecular systems, including reactants, products, intermediates and transition states on the potential energy surface (PES) were optimized using the Minnesota 2006 exchange correlation functional such as M06-2X, in conjunction with the Pople's splitvalence 6-311++G (3df, 3pd) basis set (Frisch et al., 1984;Zhao and Truhlar, 2008).The long-range van der Waals interactions between the reactive species were accounted using Grimme's empirical dispersion (GD3) corrections (Grimme et al., 2010).Previous studies have demonstrated that the current level of theory exhibits a reliable performance in addressing noncovalent interactions between gaseous molecules and in locating the transition states of atmospheric and combustion reactions (Ali et al., 2016;2018;2022;Ali, 2019;2020;Dash and Ali, 2022;2023;Ali and Balaganesh, 2023).Tight convergence criteria were applied during the wave function optimization of the reactive species, complexes, products and transition states.Unscaled vibrational frequencies at the same level of theory (M06-2X/6-311++G (3df, 3pd)) were utilized to compute zero-point energy (ZPE) corrections to the total energies of all molecular systems, to characterize the stationary points on the PES and for rate-constant calculations.Vibrational frequency analysis confirmed all positive frequencies for the reactants, complexes, intermediates and products, while a single imaginary frequency was observed for the transition states.
Additionally, single-point energy calculations were conducted at a higher-level of theory on the molecular structures optimized at a lower-level of theory to ensure an accurate description of the energetic parameters.Specifically, the CCSD(T)/6-311++G (3df, 3pd) level of theory was utilized to estimate the single-point energies of the gas-phase molecular geometries, which were initially optimized at the M06-2X/6-311++G (3df, 3pd) level of theory (Raghavachari et al., 1989).The basis set superposition error (BSSE) calculations were also performed using the counterpoise (CP) corrected method (Boys and Bernardi, 1970;Simon et al., 1996).The <Ŝ 2 > eigenvalues were monitored to evaluate the spin contamination for the wavefunction of the open-shell radicals.The T1-diagnostic values obtained at the CCSD(T)/6-311++G (3df, 3pd) level of theory were analyzed to validate the single-reference method and were found to be within the acceptable range (i.e., ≤0.02) for all important species ( Lee and Taylor, 1989).Overall, the combination of CCSD(T)//M06-2X functionals has been employed in numerous research studies, providing a reasonably accurate description of the thermochemistry and chemical kinetics of many atmospheric reactions (Ali et al., 2016;Ali et al., 2018;Ali, 2019;Ali, 2020;Ali et al., 2022;Dash and Ali, 2022;Ali and Balaganesh, 2023; Dash and Ali, 2023).

Chemical kinetic analysis
Comprehensive chemical kinetic calculations for the _ OH initiated oxidation reaction of aminomethanol (AM) using MultiWell suite of codes (Barker, 2001;Barker, 2009;Barker et al., 2023) were performed.This oxidation reaction involves a fast pre-equilibrium between the reactive species such as H 2 NCH 2 OH + _ OH and the pre-reactive complex [H 2 NCH 2 OH/OH] • stabilized by the van der Waals forces followed by a second step leading to the respective products as follows, where, k 1 and k -1 are the forward and reverse rate constant for the first bimolecular reaction and the k 2 is the rate constant for the second unimolecular reaction.The kinetic rate constants for these bimolecular (k, in cm 3 molecule −1 s -1 ) and unimolecular (k uni , s -1 ) reactions in the high-pressure limit defined by transition state theory are represented as follows, Assuming that the pre-reactive complex was in equilibrium with the reactants and was at a steady state, then the overall rate constants is expressed as; This kinetic model is reasonably correct at the high-pressure limit, where the pre-reactive complex can be stabilized by collisions with other atmospherics species.This approach was widely used in literature for the water-assisted reaction and the predicted rate coefficients are reasonably good agreement with the experimental values (Ali, 2019;Ali et al., 2022, Ali andBalaganesh 2023).
The different parameters of Eqs The temperature-and pressure-dependent microscopic rate constants k(E) have also been computed for the O 2 addition reaction to the aminomethanol radicals generated in Eq. 2. This was accomplished using the Rice−Ramsperger−Kassel−Marcus (RRKM)/master equation (ME) theory, implemented in the MultiWell suite of programs.The MultiWell code facilitates the computation of non-steady-state effects including unimolecular decomposition processes, isomerization, collision energy transfer and chemical activation for the complex rateconstant calculations.To perform these calculations, molecular and energetic parameters such as vibrational frequencies, moments of inertia and reaction barriers are required as input data.Using this data, the MultiWell suite computes sum and density-of-states, followed by the evaluation of microscopic rateconstant k(E).The RRKM/ME microscopic rate-constant k(E) is defined as follows, The details of each term can be found in the Supplementary Material.Temperature and pressure-dependent rate constants and branching ratios of the products were evaluated by incorporating N 2 gas as the bath gas.The collisional energy transfer process was addressed using the conventional temperature-dependent exponential-down model with a <ΔE > down parameter (which represents the average energy loss per the collision of the active compound with the bath gas molecule), with an approximate value of ~200*(T/300) 0.85 cm −1 (Goldsmith et al., 2012).Lennard-Jones (L-J) parameters were employed to account for the frequency of collisions between the active compound and the bath gas (N 2 ) collider.The L-J parameters for N 2 gas, specifically σ (N 2 ) = 3.74 Å and ε/k B (N 2 ) = 82 K, were sourced from the literature, while the same parameters for all wells were adopted from our previous study ( Dash and Ali, 2022).
For the barrierless reactions i.e., _ AM +O 2 →AM-O _ O and product complexes to individual product molecules Inverse Laplace Transform (ILT) method was incorporated to determine the rate-coefficients (Robertson et al., 1995).Consistent with studies on numerous analogous reactions, this method has proven effective, with the Arrhenius's activation energy equating to the critical energy of the reaction (E 0 ) (Firaha et al., 2018).Additionally, statistical rate theories that neglect non-statistical effects were applied, including slow intramolecular vibrational energy redistribution (IVR) as discussed in previous study (Mazarei and Barker, 2022).
The pressure-dependent total rate constants k bimol (T, M) for aminomethanol radical ( _ AM) + O 2 have been computed using, where, Γ represents the quantum mechanical tunneling corrections, f _ AM+O2 is the branching fraction ( f ) of the chemical reaction returning to the respective reactive species and k uni ∞ is the rate constant at the high-pressure limit.The tunneling was incorporated for the chemical activation distribution in all our chemical kinetic calculations using the keywords "CHEMACT" and "TUN".
Finally, the calculated rate constants were fitted at the highpressure limit (k ∞ ) in the temperature range of 200-400 K to the modified Arrhenius expression, which is as follows, where A is the pre-exponential factor, T is the temperature, n is the temperature exponent and E a is the activation energy.The coordinates of equilibrium geometries, vibrational harmonic frequencies and rational constants of all important species involved in the _ OH initiated oxidation reaction of AM are listed in the Supplementary Material.

Rotational conformers of aminomethanol
The oxidation reaction between the aminomethanol (AM) and _ OH radical occur through the abstraction of H-atom from three different H-bearing functional groups (i.e., -CH 2 , -NH 2 and -OH) of AM.Generally, the H atom of -CH 2 , -NH 2 and -OH groups along with the presence of lone pair of electrons on the N and O atoms facilitate the hydrogen bonding (H-bonding) between the AM and the _ OH radical.Interestingly, the -NH 2 and -OH groups of AM can freely rotate around the single bonds to form the inter molecular H-bonding with the _ OH radical based on its attacking direction.It leads to different rotational conformations for the AM-_ OH radical complex.In addition, formation of such intermolecular H-bonding in the pre-reactive and the transitions state (TS) structures has a great effect on the energetics of the reaction.Hence, we have comprehensively explored the different reaction pathways and the corresponding energy barriers using different rotational conformations of AM.It is evident from previous studies that the AM exhibits four different rotational conformations, which we named as AM0, AM1, AM2 and AM3 as presented in Figure 2 (Hays et al., 2013).It shows that the AM0 is energetically more stable and can be consider as ground state conformation of AM.The geometrical features of ground state conformation of AM are good in agreement with that of previous study (Hays et al., 2013).Other rotational conformations such as AM1, AM2 and AM3 exhibit ~0.2, 0.8 and 4 kcal/mol higher energy compared to the AM0, respectively.The trend in the conformational stabilities of AM is in accordance with the study of Weaver and co-workers (Hays et al., 2013).Hence, the AM0 is used as a reference to provide the further intricate details about the mechanistic pathways of the AM oxidation reaction with the _ OH radical.

Potential energy surface of aminomethanol + _ OH
The zero-point energy (ZPE) corrected potential energy surface (PES) for the H-abstraction reaction of the most stable conformation of aminomethanol (AM0) by the _ OH radical is depicted in Figure 3.In this reaction, the _ OH radical abstracts the H-atom from three different H-bearing moieties (i.e., -CH 2 , -NH 2 and -OH) of AM0, leading to the formation of three distinct For example, the bond angle between the H atom of the -CH 2 group and _ OH is found to be 94.6 °(i.e., ;HC-H••• _ OH = 94.6 °) in the TS-CH (see Figure 4).The presence of two inequivalent H-atoms on -NH 2 group form two different TSs.The TSs associated with the back and front H-atoms are designated as TS-NH a and TS-NH b , respectively.The bond angles are found to be ; Different rotational conformations of aminomethanol (AM).The values in the parenthesis indicate the relative energies in kcal/mol.To account for the effect of other rotational conformers on the oxidation reaction, the energies of PRCs, TSs and PORCs have been computed for the abstraction of H atom from -CH 2 , -NH 2 and -OH groups of other rotational conformations of AM (i.e., AM1 and AM2).The complete reaction profile for the H-abstraction reaction of AM1 along with the energies is presented Supplementary Figure S1.The H abstraction from the -CH 2 group of AM1 proceeds via three TSs such as TS-CH a TS-CH b and TS-CH c .Among these TSs, TS-CH a is linked to the PRC1 while TS-CH b and TS-CH c are linked to the PRC2 (see Supplementary Figure S2).However, all these-CH TSs are linked to the same product complex and lead to the formation of a single NH 2 _ CHOH radical conformation.On the other hand, the abstraction of H from -NH and -OH groups of AM1 occurs via the TSs namely, TS-NH and TS-OH, which are linked to the pre-reactive complex, PRC1.These TSs of AM1 subsequently form the N-and O-centered radicals similar to that of AM0+ _ OH radical reaction.Similarly, we evaluated the energetics of all important species of AM2+ _ OH oxidation reaction and the reaction profile is shown in Supplementary Figure S3.The optimized geometries of the reactive species, intermediates, TSs, post-reactive complexes and C-, N-, O-centered radicals of AM2+ _ OH reactions are presented in Supplementary Figure S4.
Finally, the energetics and the barrier height values of H abstraction reactions from AM0, AM1 and AM2 by _ OH radical have been compared and presented in Supplementary Table S1.It shows that the energies of PRCs of AM0, AM1 and AM2 are in the range of −4.6 to −5.8 kcal/mol.These values are in good agreement with the PRC energies of similar electronic systems.For example, Franco et al., have investigated the abstraction of H atom from the different conformers of methanediol (Franco et al., 2021).They showed that, the PRC between OHCH 2 OH and _ OH is stabilized by −5.22 kcal/mol with respect to the reactants in the most stable conformation of methanediol.-NH abstraction of alkylamines around 4.4-9.5 kcal/mol (Tian et al., 2009;Onel et al., 2013).Our computed results show that the barrier height values for the abstraction of H from -CH 2 , -NH 2 and -OH groups are in the range of 4.12-6.50,3.76-5.50and 6.5-9.25 kcal/ mol, respectively.These results also corroborate that the barrier heights are good in agreement with the previous studies.This analysis not only provides more confidence but also substantiate our choice of ab initio and DFT methods.
It is worth to mention here that, the TS energies for the major H-abstraction pathways i.e., from the -CH 2 and -NH 2 channels of NH 2 CH 2 OH + _ OH radical reaction are found to be around 1-2 kcal/ mol.Previous studies show that the _ OH radical initiated H-abstraction reactions from different reactants are favorable even with the slight positive TS energies.For example, Nguyen et al., have studied the H-abstraction from the CH 3 OH using ab initio/RRKM methods (Nguyen et al., 2019).They reported the TS energies for H-abstraction around 2.3 and 0.5 kcal/mol.Further, they satisfactorily reproduced the experimental rate constants using these TS energies.Baidya et al., have explored the H-abstraction reaction of CHF 2 CH 2 OH with _ OH radical to unveil the atmospheric implications of chlorofluorocarbons (Baidya et al., 2018).All the TSs reported in this study exhibit positive energies in the range of 0.8-3.1 kcal/mol.However, they reported that their calculated k OH value is in good agreement with the experimental value.Rahbar et al., have investigated the kinetics and mechanism of the _ OH radical initiated atmospheric oxidation of catechol over the temperature range 200-400 K (Rahbar et al., 2021).However, catechol shows high positive TS values (1.27-8.8kcal/mol) for the H-abstraction reaction.The authors pointed out that the computed rate constants are close to that of experimental value.In light of these analogous studies, we believe that the slight positive energies (~1-2 kcal/mol) for the TSs are acceptable within the troposphere, which is the focal point of interest in our present investigation.
The rate constant (k OH ) and the branching fraction values have been computed for the hydrogen abstraction reactions at 300 K using all possible conformations.The k OH values for the abstraction of H from different channels (-CH 2 , -NH 2 and -OH) of AM0 are plotted in Figure 5A and listed in Supplementary Table S2.The sum of the rate constant values from different channels of AM0 (i.e., k OH -AM0) are presented in Table 1.This data indicates that the rate constants for H abstraction from different channels gradually decrease within the temperature range 200-400 K.The k OH -AM0 value for _ OH radical reaction with AM0 rotational conformation is found to be 1.40 × 10 −12 cm 3 molecule −1 s −1 at 300 K. Similarly, the k OH for each H abstraction channel of AM1 and AM2 has been computed and plotted in Figures 5B, C and the numerical values are presented in Supplementary Tables S3,  S4.The sum of the rate constant values from different channels of AM1 and AM2 (i.e., k OH -AM1 and k OH -AM1) has been found to be 4.40 × 10 −12 and 1.39 × 10 −11 cm 3 molecule −1 s −1 , respectively at 300 K (see Table 1).This results in an overall rate constant for the hydrogen abstraction from the all AM conformations to be around 1.97 × 10 −11 cm 3 molecule −1 s −1 (see Figure 5D; Table 1).The computed total k OH value of AM is closely aligns with that of CH 3 NH 2 (1.97 × 10 −11 cm 3 molecule −1 s −1 ) (Onel et al., 2013), CH 3 NHCH 3 (6.27× 10 −11 cm 3 molecule −1 s −1 ) (Onel et al., 2013), CH 3 CH 2 NH 2 (2.50 × 10 −11 cm 3 molecule −1 s −1 ) (Onel et al., 2013), NH 2 CH 2 CH 2 OH (7.27 × 10 −11 cm 3 molecule −1 s −1 ) (Xie et al., 2014), CH 3 CH 2 OH (9.06 × 10 −11 cm 3 molecule −1 s −1 ) at 300 K.The k OH values are also computed using basis set super position error (BSSE) corrections.It is evident from Supplementary Table S5 that the BSSE corrections does not altered the rate constants appreciably.Overall, a negative temperature dependence in the k OH of AM has been observed over the temperature range 200-400 K as can be seen from Figure 5D, similar to that of amines and alcohol reactions.The energetic values coupled with the computed rate constant results clearly align with those of previous studies, which further corroborating the reliability of the computational methodology and the scheme adopted in this study.Additionally, the branching fraction values for the abstraction of H from -CH 2 , -NH 2 and  Frontiers in Chemistry frontiersin.orgaround 82%, 17% and 1% for the abstraction of H from -CH 2 , -NH 2 and -OH groups of monoethanolamine (MEA) (Xie et al., 2014).The computed branching fraction values of AM are in good agreement with that of MEA (Xie et al., 2014).The branching fraction values clearly suggest a strong preference for the C-centered H-abstraction over N-and O-centered H-abstraction.Interestingly, the computed branching fraction values also strongly adhere to the Evans−Polanyi relationship bond enthalpy of O−H > N−H > CH relationship (Evans and Polanyi, 1938).
The bimolecular reactions between the aminomethanol (NH 2 CH 2 OH) + _ OH radical and the NH 2 _ CHOH + O 2 would also be possible in their excited electronic states.It is evident from previous studies that, the computationally characterized potential energy surface (PES) and the corresponding rate constants of (for example, ethyl alcohol (H 3 C-CH 2 -OH) (Xu et al., 2019), monoethanolamine (NH 2 CH 2 CH 2 OH) (Xie et al., 2014), methanediol (OH-CH 2 -OH) (Franco et al., 2021)) similar electronic systems agrees well with that of experimental studies.However, these studies not included any excited state photochemical interactions to obtain the rate constants.Further, Al-Hashimi and co-workers reported that the interactions between the O-Anisidine and _ OH radical occurs in the ground state rather than in the excited state (Abdel-Rahman et al., 2021).Priya et al., have explored the abstraction reaction mechanism of _ OH radical with 2methoxyphenol ( Priya and Lakshmipathi, 2017).They showed that the H-abstraction reaction can occur in ground state than the in the excited state.Similar to these studies, we believe that the interactions between the aminomethanol and the gaseous species could be corresponding to ground state rather than to excited state.
In recent decades, significant advances in atmospheric chemistry have spurred the development of new theoretical approaches for exploring intricate details of ground-state chemical reactions and their underlying mechanisms.Nevertheless, an equivalent synergy between theory and experimentation remains absent in the realm of atmospheric photochemistry involving electronically excited states.The modeling of molecular photochemistry necessitates a meticulous consideration of non-adiabatic effects, specifically, the coupling between electronic states and molecular motion.This presents formidable challenges, as it contradicts several conventional approximations in theoretical chemistry.Notably, non-adiabatic effects challenge the venerable Born-Oppenheimer approximation, while classical treatments of nuclear dynamics may prove inadequate and non-equilibrium phenomena can challenge established reaction rate theories.A plethora of methodologies have emerged to address these challenges, including MCTDH (Multi Configuration Time Dependent Hartree) (Manthe et al., 1992;Beck, 2000;Meyer, 2012), trajectory surface hopping (TSH) (Tully and Preston, 1971;Hammes-Schiffer and Tully, 1994;Fabiano et al., 2008) and ab initio multiple spawning (AIMS) (Ben-Nun et al., 2000;Ben-Nun and Martínez, 2002;Yang and Martínez, 2011).However, the application of these techniques to investigate atmospheric photochemistry encounters various challenges.These encompass complexities in simulating spectroscopically relevant choices and atmospheric modeling interests, the intricate electronic architecture of multichromophoric volatile organic compounds (VOCs), the diverse excited-state dynamics triggered by solar irradiation, the protracted dynamics of excited VOCs, the prominence of intersystem crossings or collisional processes and the modulating influence of aqueous environments such as those found in atmospheric aerosols and clouds.Given that non-adiabatic dynamic approaches demand considerable computational resources beyond our current capacity and each method carries its own limitations, we envision these studies as prospects for future exploration.

Potential energy surface of NH
The aminomethanol radicals produced in the initial reactions of aminomethanol + _ OH are highly reactive and undergoes subsequent reactions with atmospheric O 2. The NH 2 _ CHOH has been considered for further reactions with atmospheric oxygen due to its large branching fraction.The zero-point energy (ZPE) corrected potential energy surface of NH 2 _ CHOH + O 2 reaction has been computed and shown in Figure 6, while the optimized geometries of intermediates, transition state structures and the product complexes (PCs) are depicted in Figure 7.The triplet oxygen molecule ( 3 O 2 ) reacts with NH 2 _ CHOH radical and barrierlessly added to the C-site to form a peroxy radical (IM-0A) intermediate, as shown in Figure 6.However, the attacking direction of O 2 molecule on NH 2 _ CHOH and the feasibility for the rotation of -NH 2 and -OH groups around the CO and CN single bonds lead to the formation of different peroxy radical intermediate rotational conformers.Similar to aminomethanol, the IM-0A also exhibits other rotational conformations, for example, IM-0B and IM-0C as presented in Supplementary Figure S5.Previous studies have shown that the rotational conformations of the peroxy radical have a negligible effect on the formation of the end products (Xie et al., 2014).Hence, we have chosen the highly stable NH 2 C(O _ O)HOH radical (IM-0A) for further studies and the reactions between other rotational conformations and O 2 have been excluded.The IM-0A strongly stabilized and located well below the reactants with a relative energy of −39 kcal/mol with respect to the reactants.The transfer of H occurred from −CH, -NH 2 and -OH groups to the O-site within the peroxy radical intermediate (IM-0A).We consider the H transfer followed by the breaking of CO or OO bonds occurs through two consecutive steps as an indirect mechanism while the same occurs in a single step as a direct mechanism.Overall, five indirect (via TS-1A, TS-3A, TS-5A, TS-7A and TS-8A) and three direct reaction paths (via TS-10A, TS-11A and TS-12A) have been studied for the NH 2 _ CHOH + O 2 reaction as shown in Figure 7.In detail, the transfer of H atom from -CH and -OH groups to the O-site of peroxy radical (IM-0A) is favored through the transition states TS-1A and TS-8A.On the other hand, the presence of two different H atoms in -NH 2 group and different attacking directions of 3 O 2 lead to three different reaction pathways, which proceed through the TSs, TS-3A, TS-5A and TS-7A for the H abstraction from -NH 2 group.Among all these TSs, the TS-1A and TS-3A are situated above the reactants with an energy of 1.5 and 0.5 kcal/mol followed by TS-5A, TS-7A and TS-8A (−1.3, −3.7 and −10.0 kcal/mol).It indicates that the H-transfer from -CH group to O-site of IM-0A requires a high energy of ~40 kcal/mol to form the intermediate IM-1A.The transfer of two different H atoms from -NH 2 group respectively requires 40 (via TS-3A) and 35 (via TS-7A) kcal/mol to form the intermediates, IM-2A and IM-3A.Further, different spatial arrangement of O 2 also facilitates a distinct transition state, TS-5A for the H transfer from -NH 2 group and it leads to the previous intermediate, IM-2A.The reaction proceed through TS-5A shows a barrier height more than 35 kcal/ mol.The H transfer process from that of -OH group requires a relatively low energy of 29 kcal/mol via the TS-8A to form IM-4A.The barrier heights for these H-transfer reactions strongly comply with that of H transfer from the similar functional groups.The intermediates IM-1A to IM-4A are stabilized by -29, −20.2, −20.6 and −16 kcal/mol with respect to the reactants, respectively.The breaking of O-O and C-O bond in IM-1A and IM-4A through the transition states TS-2A and TS-9A form the product complexes, PC-1A and PC-4A, respectively.While, PC-2A and PC-3A can be formed by the breaking of C-O bond of IM-2A via the TSs, TS-4A1 and TS-4A2.The IM-3A also leads to the same product complexes, PC-2A and PC-4A via the TSs, TS-6A1 and TS-6A2.Figure 6 shows that, the IM-1A is almost barrierlessly (0.3 kcal/mol) dissociated into the product complex PC-1A, which is composed of hydrogen bonded amino formic acid (H 2 NCOOH) and _ OH radical (see Figure 7).The PC-1A is a highly stable product complex with an energy of -101 kcal/mol, which eventually separates into amino formic acid (NH 2 COOH) and _ OH radical.While, the IM-2A, IM-3A and IM-4A exhibit a barrier height of ~13, 16 and 13 kcal/mol for the conversion into PC-2A, PC-3A and PC-4A, which are the hydrogen bond stabilized product complexes of formimidic acid (HN = C(H)-OH), formamide (H 2 N-CHO) and _ O 2 H radical.The formimidic acid and formamide product complexes (PC-2A, PC-3A and PC-4A) located above than that of amino formic acid (PC-1A) with an energy of −32, −36 and −49 kcal/mol, respectively (see Figure 6).On the other hand, the reactions proceed via the transition states, TS-10A, TS-11A and TS-12A directly form the product complexes.In detail, these TSs initiates the simultaneous transfer of H atom from -NH 2 and -OH groups to the O-site of peroxy radical and the breaking of C-O bond to form the post-reactive product complexes, PC-2A, PC-3A and PC-4A, respectively.Similar to TS-3A and TS-7A, the TS-10A and TS-11A involve the transfer of inequivalent H atoms of -NH 2 to the O-site.The energy barriers have been found to be 18, 17 and 4.3 kcal/mol for the transfer of H atom from -NH 2 (via TS-10A, TS-11A) and -OH (via TS-12A) groups to the O-site of peroxy radical.Overall, it is evident from Figure 6 that, the reaction pathway proceeding through the transition state, TS-12A to form P-4A (NH 2 CHO) is energetically more favorable followed by that proceeding through TS-10A and TS-11A to form P-2A and P-3A (NHCHOH) by the simultaneous H transfer and CO bond breaking process rather than reaction pathways that proceeding through TS-1A, TS-3A, TS-5A, TS-7A and TS-8A.
Overall, the computed results suggest that the reaction pathways, which proceed through the TS-1A, TS-3A, TS-5A, TS-7A and TS-8A from the intermediate IM0-A exhibit high energy barriers.These energy barriers are in the range of (30-40 kcal/mol) as shown in Figure 7. Hence, these reaction pathways are excluded from the further studies of reaction kinetics and branching ratio analysis due to the high energy barriers.Relatively low energy barrier pathways which occurred via the transition states TS-10A, TS-11A and TS-12A were adopted for further kinetics and branching ratio calculations.Computed rate constants for the NH 2 _ CHOH + 3 O 2 reaction at 1 atm pressure for the formation of formamide over the temperature range 200-400 K.
Frontiers in Chemistry frontiersin.org12 Nulakani and Ali 10.3389/fchem.2024.1407355The rate constants and the branching ratio values of 3 O 2 + NH 2 _ CHOH reaction have been computed for the most favorable reaction pathways (via TS-10A, TS-11A and TS-12A), which form the product compounds P-2A, P-3A and P-4A.We have treated all these reaction pathways in a master equation to evaluate the temperature and pressure-dependent rate constants.The computed rate constants at the atmospheric conditions relevant to the troposphere, i.e., at 1 atm pressure and 298 K for the overall reaction is 5.5 × 10 −12 cm 3 molecule −1 s −1 .The rate constant of O 2 + NH 2 _ CHOH was compared with those of isoelectronic systems, namely, CH 3 _ CHOH + O 2 (Silva et al., 2009) and OH _ CHOH + O 2 (Ali and Balaganesh, 2023) reactions.We have found that the overall rate constant of O 2 + NH 2 _ CHOH decreases with the temperature similar to that of OH _ CHOH + O 2 as shown in Figure 8.The total rate constant of NH 2 _ CHOH + O 2 reaction is mainly contributed from the rate constant of P-4A reaction rather than that of P-2A and P-3A (see Supplementary Figure S6).Further, branching ratio analysis suggests that 3 O 2 + NH 2 _ CHOH reaction predominantly leads to the formation of P-4A (NH 2 CHO) with the branching fraction of ~99% over the temperature range 200-400 K, which is akin to the previous analysis on similar reaction system (Ali et al., 2021).The branching fraction values for the other products are negligibly small and independent of temperature and pressure.
We have also computed the ZPE corrected PES of nitrogen centered radical, namely, _ NHCH 2 OH + O 2 reaction and presented in Supplementary Figure S8.The CH-direct, CH-indirect and OHindirect reaction pathways shows the high positive energies for the TSs around 8.9 (TS5), 27.6 (TS1) and 19.5 (TS2) kcal/mol and hindered by the strong positive energy barriers, which are around 10.5, 29.2 and 22.9 kcal/mol, respectively.The large positive energies for the TSs are not encouraging and the corresponding reaction pathways will likely not be traversed in the upper atmosphere.

Molecular dynamics simulations
We have conducted Born-Oppenheimer molecular dynamics (BOMD) simulations to investigate the formation of formamide at 300 K through the reactions involving the reactive species ( _ OH + AM and NH 2 _ CHOH radical + O 2 ) using the CP2K code (VandeVondele et al., 2005).The computational methodology details are provided in the Supplementary Material.The snapshots captured at different time intervals from the trajectory of the 3ps BOMD simulation for the H abstraction from each channel of AM are shown in Figure 9.These snapshots clearly demonstrates that the abstraction of H from -CH 2 , -NH 2 and -OH groups during the _ OH + AM reaction follows the similar paths as discussed in Section 3.2.For example, the initial geometry for the H abstraction by the _ OH radical from -NH 2 group is depicted in Figure 9 at 0 femtoseconds (fs).During the BOMD simulation, a pre-reactive complex like geometry formed between the AM and _ OH radical at 62 fs.Subsequently, we have observed the transfer of H from -NH 2 to _ OH radical occurring through a transition state-like geometry at 230 fs.Further, this forms a _ NHCH 2 OH + H 2 O complex at 250 fs.Similar mechanisms are also observed in the cases of H abstraction from -CH 2 and -OH groups of AM as illustrated in Figure 9. Additionally, BOMD simulations have

Atmospheric implications
Our previous investigation (Ali, 2019) has demonstrated the occurrence of a high-temperature reaction between ammonia (NH₃) and formaldehyde (H 2 CO) during biomass burning, resulting in the formation of aminomethanol (NH₂CH₂OH).This prompts a keen interest in elucidating the oxidation reaction mechanism and kinetics of aminomethanol, leading to the formation of various compounds, including potentially hazardous or carcinogenic byproducts.The possible atmospheric decomposition pathway resulting from _ OH radical initiated oxidation of NH₂CH₂OH is portrayed in Figure 10.We have computed the atmospheric lifetime (τ) of NH₂CH₂OH due to its interactions with _ OH radicals using the following the equation, where, k denotes the rate constant of NH₂CH₂OH + _ OH radical reaction (5.33 × 10 −11 cm 3 molecule −1 s −1 ) at 230 K, relevant to an altitude of 12 km.

The averaged concentration of _
OH radicals in the upper troposphere is denoted by [ _ OH] and has been taken as 1.0 × 10 6 molecule cm -3 akin to that of earlier studies ( Ali and Balaganesh, 2023).Computed results reveal that, the aminomethanol has a lifetime of 5 h in the presence of _ OH radicals and produce the NH 2 _ CHOH radicals as the major product when compared to _ NHCH 2 OH and NH 2 CH 2 _ O radicals.Under tropospheric conditions, the major radical product, i.e., NH 2 _ CHOH undergoes further reaction with molecular oxygen ( 3 O 2 ) to yield formamide (NH 2 CHO).We have employed k of 7.80 × 10 −12 cm 3 molecule −1 s −1 at 225 K to calculate the atmospheric lifetime of the NH 2 _ CHOH radicals.Our results show that, the lifetime of NH 2 _ CHOH is approximately 13 microseconds (μs), indicating a rapid formation of formamide under atmospheric conditions.The current study also suggests that the formations of amino formic acid (NH 2 COOH) and formimidic acid (NHCHOH) are unfavorable.It is important to mention here that the formation of formamide involves the reaction of aminomethanol with _ OH and O 2 .Analysis of NH 2 CH 2 OH+ _ OH reactions indicates that the transition states (TSs) tend to exhibit a slight positive value for the H-abstraction reactions.Consequently, these reactions pose challenges under astrochemical conditions, where temperatures typically remain below 100 K. Nevertheless, the reactions involving the carbon-centred NH 2 _ CHOH radical with O 2 may be viable in such conditions, as the corresponding aminomethanol-based peroxy radical intermediate, NH 2 CH(O _ O)OH required less than 5 kcal/mol for the formamide formation.Therefore, investigating the comprehensive reaction pathways for the formamide formation from aminomethanol under astrochemical conditions holds significant promise for future studies in this field.Generalized reaction mechanism for the formation of aminomethanol from ammonia and formaldehyde and its subsequent decomposition reaction pathways to potentially hazardous isocyanic acid.Barnes et al., have reviewed the mechanistic details and atmospheric chemistry of amides (Barnes et al., 2010).They showed that the formamide could further react with _ OH radicals and produce C-and N-centered formamide radicals as shown in Figure 10.However, Zhu and co-workers have revealed that the C-centered formamide radicals have exclusively been formed and rule out the possibility for the formation of N-centered radicals using the combination of experimental and theoretical studies (Bunkan et al., 2016).They have also found that the isocyanic acid (HNCO), which is known to be a potentially hazardous compound for the human health, is only the product formed during the NH 2 CHO + OH and O 2 reaction.
On the other hand, the branching fraction (20%) of N-centered radicals indicates the formation of _ NHCH 2 OH in significant quantities.It suggests that the N-centered radicals are also prone to further reactions with atmospheric gaseous compounds.However, N-centered radicals usually react slowly with the atmospheric oxygen and potentially leads to the formation of carcinogenic nitrosamines or nitramines through the bi-molecular reactions with the other trace compounds of the atmosphere (NO and NO 2 ).
Overall, investigating the branching fraction and ensuing rate constants associated with the hydrogen abstraction from -CH 2 , -NH 2 , and -OH groups of aminomethanol by hydroxyl radicals offers initial insights into the predominant formation of formamide.Additionally, this study provides additional insights into the formation of N-centered radicals in significant quantities.Nevertheless, conducting further investigations on the reactions of O 2 , NO and NO 2 with N-centered radicals would undoubtedly contribute to significantly understand their impact on the atmosphere.

Conclusion
In this study, we have comprehensively explored the mechanistic details of aminomethanol oxidation reaction initiated by atmospheric _ OH/O 2 using CCSD(T)/M06-2X level of theory.The NH 2 CH 2 OH + _ OH reaction predominantly yields the C-centered NH 2 _ CHOH radicals when compared to the N-and O-centered _ NHCH 2 OH and NH 2 CH _ O radicals.Our findings reveal that, the total rate constants (k OH ) for the H transfer from aminomethanol is approximately 1.97 × 10 −11 cm 3 molecule −1 s −1 at 300 K.This k OH value closely matches with that of CH 3 NH 2 (1.97 × 10 −11 cm 3 molecule −1 s −1 ), CH 3 NHCH 3 (6.27× 10 −11 cm 3 molecule −1 s −1 ), CH 3 CH 2 NH 2 (2.50 × 10 −11 cm 3 molecule −1 s −1 ), NH 2 CH 2 CH 2 OH (7.27 × 10 −11 cm 3 molecule −1 s −1 ), indicating the accuracy of our findings.Moreover, the computed atmospheric lifetime of NH 2 _ CHOH radicals is remarkably short, estimated at 13 μs, suggesting rapid reactions with atmospheric oxygen.We observed that the reaction between NH 2 _ CHOH radicals and O 2 exclusively forms formamide (NH 2 CHO) as the product, with a branching fraction of approximately 99% and a rate constant of 5.5 × 10 −12 cm 3 molecule −1 s −1 under the tropospheric conditions.The rate constant of NH 2 _ CHOH + O 2 reaction aligns well with that of HOCHOH + O 2 and CH 3 CHOH + O 2 reactions.Our Born-Oppenheimer molecular dynamics (BOMD) simulations also substantiate the formation of formamide as the prime product.Finally, we have addressed the environmental implications of exclusively formed formamide.The formamide can leads to formation of potentially hazardous compounds like HNCO on further reactions with atmospheric hydroxyl radicals.The unreacted aminomethanol _ NHCH 2 OH radicals may form the carcinogenic nitrosamines on reacting with trace N-oxides (viz., NO and NO 2 ).Consequently, it leads to increase the environmental risk factors.
AM0 radicals, namely, the carbon-centered NH 2 _ CHOH radical, nitrogen-centered _ NHCH 2 OH radical and oxygen-centered NH 2 CH 2 _ O radical.These radicals are obtained through an energetically favorable pre-reactive complex (PRC), followed by the transition states (TSs) for the transfer of the H-atom.The optimized geometries of AM0, _ OH radical, PRC, TSs, postreactive complexes (PORCs) and the final radical species of _ OH initiated AM oxidation reaction are presented in Figure 4.In detail, the _ OH radical attacks the AM0 and forms a PRC, which is stabilized by the formation of H-bonding between the H atom of _ OH radical and the O atom of AM0 as shown in Figure 4.The PRC is stabilized with respect to the reactants by −5.26 kcal/mol.Subsequently, the H-bonding stabilized _ OH radical rearranges in the reactive space of AM0 in a way to form a bond angle between the H-atom of -CH 2 , -NH 2 and -OH groups and the _ OH radical close to that of H 2 O molecule, to viably eliminate as a H 2 O molecule through a transition state.

FIGURE 3
FIGURE 3The zero-point energy corrected potential energy surface of AM0 + _ OH radical.

FIGURE 4
FIGURE 4 Optimized geometries of aminomethanol (AM0), PRC, TSs, PORCs and the radicals of _ OH initiated AM0 oxidation reaction.All the values in the parenthesis are in kcal/mol.

FIGURE 5
FIGURE 5 Calculated rate constants (k OH ) for the abstraction of H from -CH 2 , -NH 2 and -OH groups of different rotational conformations of aminomethanol.(A) AM0, (B) AM1 and (C) AM2.(D) The total rate constants from each channel of AM0, AM1 and AM2.

FIGURE 6
FIGURE 6The zero-point energy (ZPE) corrected potential energy surface of AM radical (NH 2 _ CHOH) + 3 O 2 radical.All the energy values are in kcal/mol.

FIGURE 9
FIGURE 9 Snapshots taken at different time intervals during the BOMD simulations between _ OH radical and AM0.Abstraction of H from (A) -CH 2 , (B) -NH 2 and (C) -OH groups by _ O H radical.

FIGURE 10
FIGURE 10 3, 6 were breakdown and elaborately discussed the specifics of each component in Supplementary Material to prevent redundancy from earlier research.The k and k 2 for the other plausible oxidation reactions of _ OH initiated oxidation reaction of AM such as for H 2 NCH 2 OH+ _ OH→H _ NCH 2 OH + H 2 O and H 2 NCH 2 OH+ _ OH→H 2 NCH 2 _ OH + H 2 O were also computed.

TABLE 1
Calculated temperature dependent rate constants (k OH , cm 3 molecule −1 s −1 ) for the _ O H initiated oxidation reaction of aminomethanol using the rotational conformations, AM0, AM1 and AM2.k OH -Total is the total rate constant using all the conformations.