The reaction between NH3+ and CH3COOH: a possible process for the formation of glycine precursors in the interstellar medium

L. Largo; P. Redondo; V. M. Rayón; A. Largo; C. Barrientos

doi:10.1051/0004-6361/201014057

Home

All issues

Volume 516 (June-July 2010)

A&A, 516 (2010) A79

Full HTML

Free Access

Issue		A&A Volume 516, June-July 2010


Article Number		A79
Number of page(s)		8
Section		Atomic, molecular, and nuclear data
DOI		https://doi.org/10.1051/0004-6361/201014057
Published online		19 July 2010

A&A 516, A79 (2010)

The reaction between NH₃⁺ and CH₃COOH: a possible process for the formation of glycine precursors in the interstellar medium

L. Largo - P. Redondo - V. M. Rayón - A. Largo - C. Barrientos

Computational Chemistry Group, Universidad de Valladolid, Departamento de Química Física Facultad de Ciencias, 47005 Valladolid, Spain

Received 13 January 2010 / Accepted 22 March 2010

Abstract
Context. The formation of glycine is strongly relevant to our understanding of the interstellar medium and is most accuretely studied computationally.
Aims. We carry out a theoretical study of the reactions between the radical cation of ammonia and CH₃COOH/CH₂COOH as possible processes leading to glycine derivatives.
Methods. The gas-phase reactions were theoretically studied using ab initio methods. We employed the second-order Moller-Plesset level in conjunction with the cc-pVTZ basis set. In addition, the electronic energies were refined by means of single-point calculations at the CCSD(T) level on the MP2/cc-pVTZ geometries with the aug-cc-pVTZ basis set.
Results. We report accurate potential energy surfaces for the reactions considered in this work. The different intermediate species as well as the most relevant transition states for these reactions are characterized.
Conclusions. Formation of protonated glycine from the reaction of NH₃⁺ with acetic acid is an exothermic (-9.15 kcal/mol at CCSD(T) level) barrier free process. However, the results obtained indicate that the hydrogen-transfer process forming NH₄⁺ and CH₂COOH is clearly the dominating path, in agreement with the experimental evidence. The formation of ionized glycine from the reaction of product CH₂COOH with NH₃⁺ is a quasi-isoenergetic (2.03 kcal/mol at CCSD(T) level) barrier free process that leads to a highly stable intermediate: protonated glycine.

Key words: astrochemistry - astrobiology - molecular processes - molecular data - methods: numerical

1 Introduction

There has been an increasing fascination with the complexity of interstellar chemistry and the existence of interstellar molecules of biological interest. In particular, the formation of biomolecules in space constitutes an open question about the origin of life. It has been argued that interstellar amino acid formation could precede the syntheses of more complex systems in space (Wincel et al. 2000). Astrochemical models, laboratory experiments, and theoretical studies have aided in understanding the formation of different organic molecules in interstellar clouds, the regions where stars and planetary systems could be formed. If organic matter were to have been formed in the interstellar regions by gas phase reactions or reactions catalyzed by dust grains, interstellar bodies might have deposited this matter on Earth in the pre-life era which would indicate that these bodies play an essential role in the formation of life. Among the biologically important molecules the amino acids, as key building blocks of proteins and peptides, deserve special attention. The presence of amino acids in the interstellar medium would be of crucial importance to understanding the possible pathways in the origin of life in the Universe. Nevertheless, after more than two decades of continuous searches in space no convincing detection of even the simplest amino acid, glycine, has been reported (Brown et al. 1979; Hollis et al. 1980; Snyder et al. 1983; Combes et al. 1996; Ceccarelli et al. 2000; Kuan et al. 2003; Hollis et al. 2003a,b; Snyder et al. 2005; Guélin et al. 2008).

The observational difficulties associated with the detection of amino acids in line-rich, compact, hot molecular cores can be summarized in two points:

amino acids have limited photostability (Ehrenfreund et al. 2001). Amino acids in the gas phase will probably be destroyed during the lifetime of a typical interstellar cloud, even when exposed to ultraviolet photons of relatively low energy. However, this does not preclude the amino acids in meteorites having formed in the interstellar medium; it only requires that amino acids be shielded from UV in protected environments, such as dense molecular clouds and hot cores (star-forming regions);
amino acids constitute relatively large molecules and their rotational spectrum contains weak lines due to large molecular partition functions and contamination of the target transitions by emission from many other molecules in the hot molecular cores and their accompanying extended molecular envelopes (Combes et al. 1996).

Different reaction pathways leading to glycine have been proposed involving relatively abundant interstellar organic molecules. Hoyle & Wickramasinghe (1976) suggested the formation of glycine from the reaction of methanimine (HNCH₂) and formic acid (HCOOH) since these compounds had been observed in dense molecular clouds of comparable abundance. In this context, Basiuk & Bogillo (2000), Basiuk (2001), and Basiuk & Kobayashi (2002) carried out a theoretical study of amino-acid-precursor formation in the interstellar medium considering different pathways that include the participation of methanimine and methanimine-related cations, HNCH⁺ and H₂NCH₂⁺. Their studies demonstrated that although the reactions considered are exothermic they are not feasible under the interstellar conditions because of the high energy transition states and the possibility of competitive formation of side products.

Based on the Miller experiment (Miller 1953), involving the synthesis of amino acids under conditions that simulated primitive Earth's atmosphere from simple organic molecules, Maeda & Ono (2004) considered a two-step synthetic route of glycine with no activation barriers. The route is composed of an addition reaction of ammonia with singlet methylene yielding ammonium ylide (NH₃ + CH₂ $\rightarrow$ H₂CNH₃) followed by carboxilation of the ylide with carbon dioxide (H₂CNH₃ + CO₂ $\rightarrow$ H₂NCH₂COOH).

Other authors (Bernstein et al. 2002) have reported a laboratory demonstration that amino acids can be naturally formed from ultraviolet photolysis of icy interstellar grains containing small amounts of materials such as methanol, ammonia, and hydrogen cyanide. The study suggests that at least some meteoritic amino acids are the result of interstellar photochemistry rather than formation in liquid water on an early Solar System body. In addition, different theoretical works (Woon 2001, 2002a,b) have been devoted to the study of the viability of various pathways in forming amino acids in astrophysical ices.

The possibility of synthesizing amino acids from simple interstellar species by means of gas-phase ion-molecule reactions initiated by cosmic and ultraviolet radiation was discussed by Herbst (2001) in his work about the chemistry of interstellar space.

In this context, different synthetic means of producing glycine and alanine have been proposed (Blagojevic et al. 2003; Snow et al. 2007; Jackson et al. 2005a,b) involving the joining of amine and carboxyl functional groups. In particular, theoretical and experimental works (Blagojevic et al. 2003; Snow et al. 2007) have dealt with the study of gas-phase synthesis of amino acids from smaller molecules found in the space. From these studies, the formation of both ionized/protonated glycine and $\beta$ -alanine in the interstellar medium was inferred to occur by means of the reactions of ionized/protonated hydroxylamine, NH₂,₃OH⁺ with acetic and propanoic acids, respectively. However, we note that NH₂,₃OH⁺ has not yet been detected in the interstellar medium.

To seek a viable formation pathway for producing amino acids from interstellar molecules, Jackson et al. (2005a,b) carried out a selected ion flow tube (SIFT) study of two types of ion-molecule reactions: a) reactions between a sequence of ions such as N⁺, N₂⁺, HCOOH₂⁺, CH₃COOH₂⁺ $\dots$ with neutral amine species, CH₃NH₂ and CH₃CH₂NH₂ (Jackson et al.2005a,b); and b) reactions of positively charged ammonia and amine ion fragments, NH₂⁺, NH₃⁺, and HCNH⁺ with neutral species, HCOOH, CH₃COOH, and CH₃OCHO (Jackson et al. 2005b).

We performed a preliminary theoretical study of possible ion-molecule reactions leading to precursors of interstellar glycine (Largo et al. 2004), including predictions of their reaction enthalpies at different levels of theory. One of our main conclusions was that the reactions between both NH₂⁺ and NH₃⁺ with acetic acid were promising because the formation of glycine precursors (ionized and protonated glycine, respectively) were exothermic. For these exothermic processes, and to ascertain whether energetic barriers are present, a full exploration of the potential energy surface was required. To complete our early study, we carried out a theoretical study of the reaction of NH₂⁺ with acetic acid considering both the singlet and triplet potential energy surfaces (Largo et al. 2008). Our work suggests that the formation of ionized glycine from the route considered is a feasible process under interstellar conditions (exothermic and barrier free). However, we found two more favourable channels, namely proton transfer and the formation of NH₃ + CH₂COOH⁺.

The reaction between the radical cation of ammonia and acetic acid was previously studied by selected ion flow tube experiments (Jackson et al. 2005b). The authors concluded that NH₃⁺ reacts with CH₃COOH via H atom abstraction giving NH₄⁺ + CH₂COOH, this product being more thermodynamically stable than the CH₃COO radical. In this work, we consider whether the radicals generated by the NH₃⁺ abstraction reaction, CH₂COOH, play a role in forming interstellar amino acids. On the other hand, Blagojevic et al. (2003) argued that the reaction of NH₃⁺ with CH₃COOH does not produce ionized glycine but instead the acetohydroximic acid cation, CH₃CONHOH⁺.

Following on from our previous studies (Largo et al. 2004, 2008), we focused on the reaction of the radical cation of ammonia with acetic acid as a possible gas-phase reaction that could produce interstellar protonated glycine. In addition, to address previous experimental suggestions (Jackson et al. 2005b), we considered the reaction between NH₃⁺ and the CH₂COOH radical. We present estimations of the geometries of the different species involved in both reactions and reaction enthalpies for the possible products. We also carried out a detailed analysis of the corresponding potential energy surfaces. We note that the reactants considered in this study have been identified in the interstellar medium. Acetic acid was detected in one core source within the giant molecular cloud complex Sagitarius B2 (Mehringher et al. 1997) and hot cores are rich in ammonia (Ehrenfreund 2000). Theoretical studies in the field of the interstellar chemistry are especially suitable because the conditions of the molecular clouds are appropriate for the application of theoretical methods.

2 Computational methods

The geometries and harmonic vibrational frequencies of the minima and transitions states involved in the reactions of NH₃⁺ with CH₃COOH or CH₂COOH were computed at the second-order Moller-Plesset level in conjunction with the cc-pVTZ (correlation-consistent polarized valence triple-zeta) basis set developed by Dunning (1989). The nature of each optimized structure on the corresponding potential surfaces was checked by calculating vibrational frequencies; these computations also allow us to estimate zero-point vibrational energy (ZPVE) corrections.

Electronic energies were refined by means of single-point calculations at the CCSD(T) level (coupled cluster single and double excitation model augmented with a non-iterative treatment of triple excitations) (Raghavachari et al. 1989) on the MP2/cc-pVTZ geometries with the aug-cc-pVTZ basis set, which also includes diffuse functions. In this case, ZPVE corrections were included at the MP2/cc-pVTZ level.

The intrinsic reaction coordinate (IRC) technique (Gonzalez & Schelegel 1989, 1990) was used to check the connections between transition-state structures and adjacent minima.

All calculations were performed with the Gaussian-98 program package (Frisch at al. 1998).

3 Results and discussion

3.1 Reaction of NH₃⁺ with CH₃COOH

Given the multiplicity of the reactants, NH₃⁺ (²A $_{2}^{\prime \prime }$ ), CH₃COOH (¹A $^{\prime }$ ), the reaction takes place on the doublet potential energy surface; we therefore characterized the possible intermediate species as well as the relevant transition states in this potential energy surface. For systems involving polyatomic molecules of the present level, building an explicit potential energy surface, is difficult because of a large number of possible competitive processes, which are difficult to account for. Consequently, we limited our study to the search for the most important stationary points and estimates of energetic feasibility for some selected reactions.

$\begin{figure} \par\includegraphics[width=4.8cm,angle=-90]{14057fg1a.ps} \includegraphics[width=6.1cm,angle=-90]{14057fg1b.ps} \end{figure}$	Figure 1: MP2/cc-pVTZ optimized geometries for the different (NO₂C₂H₇)⁺ minima. Distances are given in angstroms and angles in degrees.
Open with DEXTER

Figure 1 depicts the optimized geometrical parameters, at the MP2/cc-pVTZ level of theory, of the (NO₂C₂H₇)⁺ isomers that are relevant to the reaction of NH₃⁺ with CH₃COOH. The relative energies, including zero-point vibrational energy differences (MP2/cc-pVTZ level), at both MP2/cc-pVTZ//MP2/cc-pVTZ and CCSD(T)/aug-cc-pVTZ//MP2/cc-pVTZ levels, are collected in Table 1. All the intermediates shown in this table present real vibrational frequencies and are therefore true minima on the potential surface.

From the examination of the different structures shown in Fig. 1, it can be inferred that isomer I1 is obtained from the direct approach of the hydrogen atom of the NH₃ group to the carbonylic oxygen of acetic acid. On the other hand, isomer I2 results from the interaction of the hydrogen atom of NH₃⁺ with the carbonylic oxygen of CH₃COOH and abstraction of the hydrogen atom of the methyl group to the nitrogen atom. Isomer I5 is produced by the approach of one hydrogen atom of NH₃⁺ to the hydroxylic oxygen of acetic acid. Isomers I6 and I7 correspond to the interaction between the nitrogen atom of NH₃⁺ and the carbonylic oxygen of CH₃COOH. Finally, isomers I3 and I4 are produced by the approach of the nitrogen atom of the NH₃⁺ ion to the carbon atom of methyl group of acetic acid and abstraction of hydrogen atom of the methyl group to carbonylic oxygen.

Analyzing the optimized geometries from a quantitative viewpoint, we can see that, apart from isomer I2, all intermediates depicted in Fig. 1 have similar C-C bond distances (in the range of 1.486 Å to 1.493 Å), which are typical of carbon-carbon single bonds. The C-C bond distance in isomer I2, 1.445 Å, suggests a certain degree of double character.

Table 1: Relative energies (kcal/mol) for the relevant intermediates.

In term of energetics, Table 1 shows that the two levels of theory employed predict similar stability order for the different intermediates. The CCSD(T)/aug-cc-pVTZ level tends to predict slightly lower relative energies than the MP2/cc-pVTZ level. At both MP2/cc-pVTZ and CCSD(T)/aug-cc-pVTZ levels of theory, the global minimum is predicted to be clearly isomer I2. The stability order of the different isomers is (>means more stable than)

$\begin{eqnarray*}{\rm I2 > I1 > I3 > I4 > I6 > I5 > I7}.$

In principle, acetic acid may react with the radical cation of ammonia to produce protonated glycine, although other different products could be formed:

$\rm NH_{3}^{+} + CH_{3}COOH$ $\longrightarrow$

$\left\{ \begin{array}{lr} \rm NH_{3}CH_{2}COOH^{+} (c1) + H & [A] \\ \rm NH_... ...{3}COOH^{+} & [I]\\ \rm NH_{4} + CH_{2}COOH^{+} & [J].\\ \end{array}\right.$

As we previously mentioned, because of the typical conditions reigning in the interstellar medium (basically low density and low temperature), only reactions that are exothermic and barrier free may play a significant role in interstellar chemistry. Accordingly, we first will calculate the reaction energies for the different processes considered.

Table 2: Relative energies (kcal/mol) for the possible products.

In Table 2, we collected the relative energies with respect to reactants, at MP2/cc-pVTZ//MP2/cc-pVTZ and CCSD(T)/aug-cc-pVTZ//MP2/cc-pVTZ levels of theory, of the possible products that can be formed in the reaction between NH₃⁺ and CH₃COOH. Zero-point vibrational energy differences were included (MP2/cc-pVTZ level). In some cases, two possible conformations (referred to as a quasi cis conformer) and t- (a quasi trans conformer) are considered. The optimized geometries for reactants and products are given, as Online Material, in Fig. A.1. Examination of the energetics showed in Table 2 allows us to infer that the formation of protonated glycine, through channels A and B, leading, respectively, to the lowest conformation c1 (consisting of an intramolecular hydrogen bond between a hydrogen atom of NH₃ and the carbonylic oxygen) and to the second most stable conformer c2 (where the hydrogen bond is connected to the less basic hydroxilic oxygen) of protonated glycine are exothermic processes. However, the most favorable process, from thermodynamical arguments, is the formation of ammonium cation and CH₂COOH (channel C) in consonance with the SIFT experiments (Jackson et al. 2005b). In this case, we emphasize that the reaction may also form the CH₃COO radical via extraction of the carboxylic H, although this was reported to be less thermodynamically and structurally stable (Schalley et al. 1998; Butkovskaya et al. 2004). Our results for this process are in agreement with this assertion. The relative energy with respect to the reactants for the formation of NH₄⁺ + CH₃COO is ${\Delta} E = - 10.91$ kcal/mol at CCSD(T) level and -6.95 kcal/mol at MP2 level. On the other hand, the reaction between NH₃⁺ + CH₃COOH to produce ionized glycine (Channel D) is a slightly exothermic process (quasi-isoenergetic). Following the comment of Blagojevic et al. (2003), we also considered the formation of acetohydroximic acid cation CH₃CONHOH⁺ and CH₃CONH₂OH⁺ produced by inserting NH₃⁺ into the C-OH bond and the subsequent elimination of H₂ or H, respectively (channels E and H). Nevertheless, these reactions clearly seem endothermic at both levels of theory ( ${\Delta} E=31.81$ kcal/mol and 48.62 kcal/mol respectively at CCSD(T) level). Charge transfer products and the remaining processes studied are endothermic.

$\begin{figure} \par\includegraphics[width=4.8cm,angle=-90]{14057fg2a.ps} \includegraphics[width=5.7cm,angle=-90]{14057fg2b.ps} \vspace*{2.5mm} \end{figure}$	Figure 2: MP2/cc-pVTZ optimized geometries for the relevant transition states involved in the reaction of NH₃⁺ (²A $_{2}^{\prime \prime }$ ) with CH₃COOH (¹A $^{\prime }$ ). Distances are given in angstroms and angles in degrees.
Open with DEXTER

The exothermicity is a necessary, but not a sufficient argument for a reaction to be plausible under interstellar conditions. In addition to this energetic requirement, interstellar reactions should be barrier free. We therefore carried out a full exploration of the potential surface, to analyze the feasibility of the reactions studied as a source of precursors of glycine in the interstellar medium.

$\begin{figure} \par\includegraphics[width=9cm,clip]{14057fg3.ps} \end{figure}$	Figure 3: Energy profile, in kcal/mol, for the reaction of NH₃⁺ (²A $_{2}^{\prime \prime }$ ) with CH₃COOH (¹A $^{\prime }$ ) at the CCSD(T)/aug-cc-pVTZ//MP2/cc-pVTZ and MP2/cc-pVTZ//MP2/cc-pVTZ (in parentheses) levels (including zero-point vibrational energy differences).
Open with DEXTER

Figure 2 shows the geometries of the more relevant transition states, whereas the energy profile for the NH₃⁺ + CH₃COOH reaction is represented in Fig. 3. As already mentioned, the MP2 and CCSD(T) energetics are in good agreement. Reaction between NH₃⁺ and the acetic acid starts with the barrierless approximation of the radical cation of ammonia to the carbonyl oxygen atom in CH₃COOH giving rise to a rather stable intermediate I1 ( ${\Delta} E = -28.8$ kcal/mol at CCSD(T) level)

$\begin{eqnarray*}{\rm NH_{3}^{+} + CH_{3}COOH \longrightarrow I1}. \end{eqnarray*}$

The exothermic formation of intermediate I1 produces an energy reservoir that is used as the reaction proceeds toward the products. After forming the intermediate I1, different processes can take place. Intermediate I1 can proceed to the most stable conformer of protonated glycine through an exothermic and barrier free path, a.1: