The dissociation of glycolate — astrochemical and prebiotic relevance †

On the basis of mass spectrometric experiments and quantum chemical calculations, including detailed kinetic and dynamics calculations, we report the unimolecular dissociation of an isolated glycolate anion. The dominating processes are: loss of formaldehyde; loss of carbon monoxide; loss of carbon dioxide; and loss of a hydrogen molecule, with the latter having the lowest energetic threshold. At higher energies, CO loss is the dominating reaction. The loss of CO may be followed by a second CO loss, leading to the H H2O complex in close mechanistic relationship to the Nibbering reaction. The results provide valuable insights into possible mechanisms for interstellar and prebiotic formation of glycolate via the reverse of the unimolecular dissociation reactions. We propose that the addition of the complex of OH and CO to CH2O is the most feasible route to gas phase synthesis of glycolate, since all species are abundant in interstellar space.


Introduction
Among the large number of organic molecules that are plausibly involved in the first prebiotic reactions, glycolic acid (HOCH 2 COOH), being the smallest a-hydroxy carboxylic acid is of special interest.This interest results from the fact that the molecule not only is involved in present biological processes, in particular since the conversion of glycolic acid leads to the simplest a-amino acid, glycine, but also was very likely present on the early Earth.
The hypothesis that animate matter originates from inanimate matter through a spontaneous and gradual increase in molecular complexity was first formulated in 1924 by Alexander Oparin, 1 and is referred to as chemical evolution.The underlying so-called chemical continuity principle states that there is a gradual increase in complexity during this part of evolution.The process started already when atoms combined to form small inorganic and organic molecules that in turn formed larger organic molecules.Eventually, macromolecules formed and molecular complexes were built and became involved in metabolic networks-which in turn lead to cellular life.After the formation of the first unicellular living organisms, one may say that chemical evolution turned into biological evolution.
Prebiotic organic compounds on the early Earth can be divided into two groups depending on whether they were of terrestrial or extraterrestrial origin. 2 Terrestrial sources include endogenous organic synthesis driven by different energy sources (UV radiation, electric discharge, hydrothermal energy) and have been subject to laboratory simulations.Extraterrestrial sources of organic material include cosmic dust, meteorites and comets, and are subject to direct examination because these objects continue to deliver organic molecules to Earth.
Glycolic acid (alongside with formic and lactic acid, glycine and alanine) is one of the dominant products in electric discharge syntheses in strongly reducing atmospheres (mixtures of CH 4 , NH 3 or N 2 , and H 2 O, with or without H 2 ). 3 In mildly reducing and non-reducing atmospheres, the yields of hydroxy acids (as well as amino acids) are low. 3If such conditions were prevailing, one should consider extraterrestrial objects as more likely primary sources of organic compounds.Glycolic acid has been detected in the Murchison meteorite. 4,5It is one of the most studied meteorites due to its large mass (>100 kg), the fact that it was an observed fall, and that it belongs to a group of meteorites called carbonaceous chondrites, which are rich in organic compounds.The formation of complex organic molecules in interstellar space has been considered to occur mainly in the gas phase by bimolecular radical or cationic reactions, since reactions without significant activation energies are essential for efficient synthesis under the extremely low pressures and temperatures encountered in typical molecular clouds. 6Interestingly, the existence of glycolic acid or glycine in the much-studied giant molecular cloud Sgr B(N) has not yet been confirmed, 7,8 although related compounds have been detected including acetaldehyde, 9 acetic acid, 10 glycolaldehyde, 11 and methyl formate. 12ore recently, it has been realized that gas-grain processes also should be considered to explain the pre-comet or pre-meteor formation of organic molecules. 13In this respect it is relevant that laboratory studies show that glycolic acid is efficiently formed by UV radiation of organic ices. 14,15In addition to gas-grain and radical and cationic bimolecular reactions, it is also necessary to include anionic reactions.Traditionally, there has been little interest in this topic, 16 but recent observations of organic anions 17 in giant molecular clouds have challenged this view.For this reason it is also pertinent to look for plausible anionic routes to key prebiotic molecules.In the present study, we have investigated the deprotonated form of glycolic acid, the glycolate ion, in order to elucidate the energetics and mechanisms of its unimolecular dissociation.
Direct investigation of the gas phase reactions that may lead to a given target molecule is complicated, even for a small molecule such as a glycolate anion, since the number of potential combination of reactants leading to this molecule is large.Studying the reverse reaction may in many cases be a more rewarding strategy-an approach that is analogous to the retrosynthetic analysis 18 used in modern organic synthesis; consequently, the spontaneous dissociation of the molecule of interest can be studied to obtain insight into the reverse synthetic pathways to the molecule.0][21] Very recently, a partial computational study on some possible dissociation mechanisms has also been published. 22espite that these studies contain useful details and provide some chemical insight, none of them is comprehensive, and the mass spectra presented in the two experimental publications are widely different.The purpose of our contribution is therefore to clarify the unimolecular chemistry, the associated mechanistic features, and the energetics in order to provide a full and consistent picture.Based on this, our second goal is to identify likely reactions for the interstellar synthesis of glycolic acid.

Low-energy collisional activation mass spectrometry
The experiments were conducted using a three-sector mass spectrometer with quadrupole/hexapole/time-of-flight (QHT) geometry (QTOF 2, Micromass/Waters, Manchester, U.K.) equipped with an electrospray ionization (ESI) source operated in the negative ion mode.Glycolic acid solution was purchased from Sigma-Aldrich, and no further purification was done.Solutions were prepared by diluting the glycolic acid with 1 : 1 (by weight) water-methanol down to a concentration of 50 mg ml À1 .The solution was injected at a rate of 0.01 ml min À1 into the electrospray ionization (ESI) source via a syringe pump.Collisionally induced dissociation (CID) experiments were performed at two different pressures for mass selected deprotonated glycolic acid (m/z = 75).Argon was used as the collision gas at nominal pressures set to 4.8 Â 10 À3 mbar and 6.4 Â 10 À4 mbar for the high and low pressure experiment, respectively.The gas inlet was controlled by a leak valve.
In addition, energy-resolved CID mass spectra were obtained by varying the central electrical potential difference between the quadrupole and hexapole units.The correspondence between the potential readout and the actual value has been calibrated, and is known to be accurate within 0.1 V.By plotting and analyzing the breakdown curves-the relative abundance of each fragment ion as a function of the center-of-mass energy E CM (i.e., the fraction of E lab that is available for inelastic scattering, E CM = E lab [M/(M + m)], where M is the mass of the neutral collision gas molecule and m is the mass of the parent ion) was obtained.Attempts to estimate the threshold energies, E 0 , for the most abundant fragment ions were done by a deconvolution procedure of the energy-resolved, collision-induced, dissociation cross-sections. 23,24antum-chemical calculations Gaussian 09, the most recent version of the Gaussian program suite, was used for most quantum-chemical calculations. 25nitial geometry optimizations were performed using B3LYP/ 6-31G(2df,p).The results of frequency calculations (N.B. specifically monitoring the number of imaginary frequencies) were used to identify the nature of the stationary points (minimum or first-order saddle points).Intrinsic reaction path calculations (IRC) were employed for all transition state geometries to verify that they actually connect the presupposed minima.To obtain more accurate estimates of the thermochemical quantities we took advantage of the G4 (Gaussian-4 theory) compound method, 26 for which geometry optimization is performed at a moderate level of theory, B3LYP/6-31G(2df,p).Then the equilibrium structures obtained are subject to a sequence of single point energy calculations; CCSD(T) calculations are performed with a moderate-sized basis set as are MP4 calculations with a relatively large basis set.Finally, the results of the calculations are combined using an extrapolation scheme, also including zero point vibrational energy corrections, to approximate the energies of more expensive calculations; final energies are estimated to be accurate within AE10 kJ mol À1 .
Multireference quantum chemical methods are often essential in addressing questions of energetics and even geometries of molecular species with multiconfigurational structures.Methylene (CH 2 ) and closely related species, such as the hydroxymethyl anion and the radical of present relevance, have a long and controversial history in quantum chemistry, 27 and are generally only reliably described using multireference methods.Second-order Generalized van Vleck perturbation theory (GVVPT2) 28 is a particularly robust variant of multireference perturbation theory 29 that can be used to describe whole potential energy surfaces, including those of excited electronic states.GVVPT2 has been used in this study to provide additional insight into the electronic structures of the hydroxymethyl anion and the neutral radical.Although GVVPT2 also supports incomplete model spaces, complete model spaces of the CASSCF variety were used in this study.Specifically, a 10 electron, 9 orbital CASSCF was used to describe the anion; this generated an active space of 2699 configuration state functions This journal is c the Owner Societies 2013 Phys.Chem.Chem.Phys., 2013, 15, 16615--16625 16617 (CSFs) in 1 A 0 symmetry.The augmented cc-pVTZ single particle basis set (aug-cc-pVTZ) was used for all calculations. 30All single and double excitations related to any one of the reference CSFs were included in the description of the external space; this resulted in a space spanned by 875, 039, 238 spin-and spaceadapted CSFs.Geometry optimizations were performed using an energy-based algorithm to a precision of 2 Â 10 À4 a.u. in coordinates and 4 Â 10 À8 a.u. in energy.GVVPT2 and CASSCF calculations were performed using a local electronic structure package, referred to as UNDMOL. 31elected CCSD(T)/aug-cc-pV(TQ)Z//CCSD(T)/aug-cc-pVTZ calculations, where aug-cc-pV(TQ)Z denotes basis set extrapolation after the extrapolation scheme of Helgaker et al., 32 were performed using Molpro 2012.1, [33][34][35][36][37][38] to corroborate the G4 and GVVPT2 predictions.
Product distribution of the most important first generation dissociations as a function of effective temperature after the collision cell was calculated using a master equation model with a simplified potential energy surface, as shown in the ESI.† Microcanonical rate coefficients for the loose transition states were calculated using inverse Laplace transform of Langevin capture rates; for the tight transition states, RRKM theory with rovibrational data from the quantum chemical calculations was used.In the simulations, 1 mbar of Ar was used as bath gas, but the calculated product distributions showed only minimal pressure dependence.The product distribution calculations were performed using MESMER 1.0. 39he H 2 loss results from a 1,2-elimination mechanism, as reported by Baker and Gabryelski on the basis of the mass shifts observed in their isotopic labelling experiments. 19The fact that H 2 loss is only evident in our CID mass spectra recorded at relatively high pressure and low collision energy is indicative of a process requiring multiple low energy activating collisions, i.e. slow heating of the ions.This is also consistent with the fact that the peak resulting from H 2 loss is the most abundant one in the mass spectra of Baker and Gabryelski, while it is not seen in the spectra reported by Bialecki et al. 20 In the former case, the experiments were conducted using an ion trap under slow heating conditions, while in the latter case the experiments were done using a triple quadrupole mass spectrometer under near single collision conditions at relatively high collision energy (E lab = 20 eV, Ar collision gas).In addition to H 2 loss, CO loss (giving rise to m/z 47) is the only other major fragmentation that can be inferred from the CID spectra from the ion trap experiment, indicating also a relatively low energy threshold for this process.In the spectra obtained under higher energy single collision conditions, both seen in Fig. 1a and in agreement with the observation by Bialecki et al., the peak at m/z 45, hardly seen in the spectra of Baker and Gabryelski, dominates.Bialecki et al. attribute this to the hydroxycarbonyl ion, as discussed in more detail by Sheldon and Bowie, 40 a signature peak in CID mass spectra of a-hydroxy carboxylic acids, [41][42][43] as carefully demonstrated by isotopic labelling and quantum chemical calculations. 20In the case of the glycolate ion, this is most likely due to loss of formaldehyde, although sequential CO + H 2 loss cannot be ruled out.

Mass spectra
The remaining product ions (m/z 19 and 17) have the elemental compositions H 3 O À and OH À , respectivelyalthough the corresponding fragmentation reactions are not necessarily trivial.
The peak at m/z 31, corresponding to CO 2 loss, is seen in both spectra.We note that this peak is neither seen in the ion trap spectra of Baker and Gabryelski nor the triple quadrupole spectra of Bialecki et al.
The peak at m/z 19 is seen in the spectrum presented by Bialecki et al., but was not commented on by them. 20The mass range of the ion trap used in the experiments by Baker and Gabryelski was not set wide enough to include the low mass region.Interestingly, a peak at m/z 35, probably due to OH À Á H 2 O, was observed but not commented on.The H 3 O À ion (m/z 19) differs in mass from the precursor ion by 56 Da, and is most likely generated by the consecutive losses of two CO molecules from glycolate.The existence of H 3 O À was first reported in 1982 by Paulson and Henchman, 44 when it was observed as a product ion in the reaction OH À ÁH 2 O + H 2 -H 3 O À + H 2 O in a tandem mass spectrometry experiment.In 1983, Kleingeld and Nibbering generated H 3 O À in a FT-ICR mass spectrometer by the reaction between OH À and formaldehyde, resulting from dehydrogenation.The results of the experiments in the presence of D 2 O provided a further understanding of both the structure of the H 3 O À ion and the mechanism of its formation.These experiments showed that the hydrogen atoms in H 3 O À are not equivalent, and the ion structure is best described as a water-solvated hydride ion, H À ÁH 2 O. 45 The formation reaction was summarized as follows 45 and proceeds through formation of a collision complex [OH À Á H 2 CO], which is followed by proton transfer.The complex formed upon the proton transfer [H 2 OÁHCO À ] is not likely to dissociate into H 2 O + HCO À since water is more acidic than formaldehyde.The intermediates of this reaction have the formula CH Further evidence about the kinetics and energetics of the dissociation of glycolate was obtained from energy resolved CID (Fig. 2).Inspection shows that the onset of m/z 47 is slightly lower than for m/z 45, but m/z 47 quickly becomes more intense.Both ions have onset in the region E CM = 2.5-3 eV.Attempts to accurately determine the energetic thresholds by applying the program L-CID, developed by Chen and co-workers, 24 failed due to the complexity of the dissociation reactions (N.B.L-CID is most reliable for one or two fragmentation reactions).

Isomers
The glycolate anion may exist in several isomeric forms.Starting from glycolic acid, deprotonation may occur from three different sites: each of the two non-equivalent hydroxyl groups and the alpha carbon.For each deprotonated form there exist a number of different conformers.The ten most stable isomers, according to our computational modelling, are summarized in Fig. 3. Starting geometries were generated by systematically rotating around the rotatable bonds in a molecule by taking advantage of this functionality provided by OpenBabel. 47

Isomerization
According to the results of our quantum chemical survey, the majority of the dissociation reactions occur from the immediate precursors 4 and 5.However, the most stable isomer, 1, corresponds to deprotonation of the carboxylic group.For that reason, the dissociation processes are induced by isomerization leading from 1 to 4 or 5 involving successive intramolecular proton transfer and bond rotations.The energy diagram describing the relevant isomerization reactions is shown in Fig. 4. Proton transfer from the hydroxymethyl group to the carboxylate group of isomer 1 gives rise to intermediate isomer 7. Continuing, the subsequent steps are: C-C bond rotation giving isomer 5; and, in order to reach 4, a final rotation around the C-O bond is  required.These details are in accordance with previous studies. 48,49ss of H 2 CO and CO 2 Loss of formaldehyde (H 2 CO) results in (C, O 2 , H) À (m/z 45).As mentioned above, the corresponding peak in the mass spectrum is diagnostic of a-hydroxy carboxylic acids and the product ion has been proposed to be the hydroxycarbonyl anion.20 On the basis of our calculations, we present a mechanism (see Fig. 5) that in broad terms is in agreement with earlier mechanistic proposals.20,43 It involves isomer 5 as the first intermediate, which rearranges into the ion-neutral complex [OCOH À , CH 2 O] (26).The direct dissociation of this complex leads to the hydroxycarbonyl anion (11) plus formaldehyde, for which the reverse process is barrierless.The calculations predict an energetic threshold for OCOH À + CH 2 O of E crit 265 kJ mol À1 .
In contrast to larger a-hydroxy acids, glycolate is unique by lacking b-hydrogens, thereby making the alternative elimination of formic acid from the complex more difficult. 43More importantly, however, is the fact that formic acid is much more acidic than formaldehyde, thereby effectively blocking this pathway.
The hydroxycarbonyl ion (COOH À ) and formate (HCOO À ) are the two isomeric anions resulting from deprotonation of formic acid.When formic acid reacts with OH À in the gas phase, both ions are produced, but in a ratio of 11 : 1 in favour of the formate ion. 50Our quantum chemical calculations predict that the formate ion is 145 kJ mol À1 more stable than the hydroxycarbonyl ion.The rearrangement of the hydroxycarbonyl ion to formate, COOH À -HCOO À , is, in spite of this, unfavourable under low energy CID conditions, since it requires overcoming an additional barrier of 122 kJ mol À1 (Fig. 6).
It should also be mentioned that glycolate isomers 4, 5, 6 and 7 are of nearly the same relative energy (150 kJ mol À1 ) and are also potential direct precursors for formaldehyde plus hydroxycarbonyl anions.Isomer 4 is the least preferable since, during its dissociation, the hydroxycarbonyl anion is not formed in its most stable form (with dihedral angle j(O-C-O-H) = 01) but instead in a form with j(O-C-O-H) = 1801, which is 5 kJ mol À1 higher in energy (at the G4 level of theory).The dissociation of the remaining isomers 5, 6, and 7, on the other hand, gives the most stable conformer of the hydroxycarbonyl anion, 11 directly (at 265 kJ mol À1 ).
Loss of carbon dioxide, giving either À CH 2 OH or CH 3 O À at m/z 31, is evident from the mass spectra presented in Fig. 1; and it appears from Fig. 2 that this process has a critical energy similar to the other major fragmentations.We have investigated alternative reaction pathways by which the CO 2 loss may occur.Starting from the most stable isomer of glycolate, 1, we first investigated the direct dissociation leading to À CH 2 OH (29) + CO 2 at 309 kJ mol À1 (Fig. 7), which at first sight appears  to be feasible.According to the calculations, the reverse reaction is barrierless.To our knowledge, the hydroxymethyl anion has not been observed, and its existence has been seriously questioned by Bowie on the basis of MP2 calculations. 51For this reason we conducted high level calculations to investigate the electronic structure of this species in greatest possible detail, including the use of multiconfigurational methods.B3LYP/6-31G(2df,p), B3LYP/ 6-311++G(d,p) and GVVPT2/aug-cc-pVTZ predict it to have a potential energy minimum, with all methods giving very similar geometries (bond lengths within 0.01 Å).The vertical electron detachment energies (CCSD(T)/aug-cc-pVTZ//GVVPT2/aug-cc-pVTZ and GVVPT2/aug-cc-pVTZ) are 66 and 48 kJ mol À1 , respectively.However, it turns out that the structure of the CH 2 OH radical is very different from that of the anion.As a result, the corresponding adiabatic detachment energy is À4 and À17 kJ mol À1 , meaning that the anion is inherently unstable towards spontaneous electron loss.In our experiments, in which we register product ions several milliseconds after fragmentation, it seems unlikely that any ions of the CH 2 OH À structure will survive the flight all the way to the ToF mass analyzer.We also investigated the lowest-lying triplet electronic state of the À CH 2 OH anion, which is 55 kJ mol À1 and 51 kJ mol À1 higher in energy than the singlet at CCSD(T)/ aug-cc-pVTZ//GVVPT2/aug-cc-pVTZ and GVVPT2/aug-cc-pVTZ levels of theory, respectively.
According to our G4 computational modelling, the methoxy anion, CH 3 O À (31), is 132 kJ mol À1 lower in potential energy than the isomeric (putative, at best short-lived) hydroxymethyl anion (29), but a considerable isomerization barrier of 94 kJ mol À1 separates them (Fig. 7).In order to reach the products CH 3 O À + CO 2 but avoid this barrier, rearrangement mechanisms requiring hydrogen rearrangements prior to dissociation were considered.During the dissociation leading to formaldehyde loss, we had already noted the intermediacy of a key [OCOH À , CH 2 O] complex, which may exist in several forms of similar energy (26-28, Fig. 8).It turns out that this complex may dissociate after intracomplex hydride transfer: ( These features are illustrated in Fig. 8, and the transition state geometry is shown in Fig. 9.According to our calculations, this rearrangement mechanism, having an intracomplex hydride transfer as the limiting step, has a significant barrier of E crit = 300 kJ mol À1 , leading to the products at 177 kJ mol À1 (confirmed by IRC calculations).However, the nature of the transition state led us to consider an alternative outcome, namely that the same transition state geometry (TSQ) could serve as a link to the even more stable products, i.e.
This would require a situation where the transient product complex [CH 3 O À ,CO 2 ] would be sufficiently long-lived and that the relative geometric requirements would allow for back-transfer  of the hydride to the CO 2 unit, but this time to the central carbon atom giving formate as the ultimate ionic product.The interesting hydroxycarbonyl-formate dichotomy has previously been discussed by Sheldon and Bowie, 40 and the present results provide a catalytic mechanism for the isomerization.The formate-formaldehyde product pair has a potential energy of 120 kJ mol À1 , and we now see that this also would correspond to an alternative mechanism for the formaldehyde loss already discussed.In order to provide better insight into this interesting possibility, we performed Born Oppenheimer molecular dynamics calculations.A total of 150 trajectories were calculated for which we varied the initial conditions by picking 50 samples of a Boltzmann ensemble for each of three temperatures (Table 1).
It turned out that the majority of trajectories (82 out of 150) produced the methoxy anion-carbon dioxide product pairs, but a total of 48 trajectories led to the formation of the formate anion-formaldehyde product pair.The statistically small sample neither allows for definitive conclusions about product distribution nor the effect of temperature.A priori it would appear to be reasonable that higher internal energies at TSQ would discriminate against back-transfer of the hydride and thereby this formaldehyde loss mechanism.Separate calculations show that the hydride transfer from CH 3 O À to CO 2 is barrierless.

Loss of H 2
At high pressure and low collision energy we observe a peak m/z 73, which corresponds to an ion resulting from the loss of a hydrogen molecule.According to our computations, the two hydrogen atoms are released from adjacent atoms, a 1,2-elimination (cf.Fig. 5).This finding is in agreement with experiments on isotope-labeled glycolate. 19The low energetic threshold for this reaction of 220 kJ mol À1 is in good agreement with the appearance of the m/z 73 ion in CID mass spectra obtained by slow heating of the ions (see above).Furthermore, the unfavourably tight transition state structure TSD also explains the disappearance of this reaction channel at slightly higher energies when processes with looser transition states emerge, as evident from the CID experiments conducted at low pressure/high energy; see further discussion of the reaction kinetics in the Discussion and Conclusion section.It should also be mentioned that the alternative 1,4-elimination of H 2 from 9 has a considerably higher barrier at 310 kJ mol À1 .
The product of the hydrogen molecule elimination -glyoxylate (13) -could in principle fragment further by eliminating carbon monoxide, thereby forming the formate ion and contributing to the m/z 45 peak.However, our calculations indicate that the barrier to carbon monoxide elimination from the glyoxylate anion is prohibitively high, being 18 kJ mol À1 higher in energy than the barrier for the formaldehyde elimination from the glycolate anion (see Fig. 5) as well as requiring an extremely tight transition state geometry (TSE).For this reason we do not consider that this consecutive H 2 + CO loss mechanism contributes to the m/z 45 signal to any significant degree, in contrast to a previous proposition. 53

Loss of CO
The peak at m/z 47 in the CID mass spectra corresponds to the neutral loss of CO.On the basis of experiments with  isotopically labelled glycolic acid, Baker and Gabryelski found that the major source of CO is from the carboxylate group. 19heir reaction mechanism involves nucleophilic attack of a carboxylate oxygen on the a-carbon followed by C-C bond cleavage.The intermediate ion formed in that process finally eliminates CO.
On the basis of our computational investigation, we suggest a mechanism that has these essential features but is somewhat more complex.The first step is proton transfer accompanied by conformeric rotation, which leads from the most stable glycolate isomer 1 to isomer 5 (Fig. 3).After that, nucleophilic attack of the hydroxyl group at the a-carbon takes place.This process (via TSF at 247 kJ mol À1 ) results in irreversible OH insertion into the C-C bond, which eventually leads to the elimination of CO (see Fig. 6).
We also investigated the possibility of further fragmentation of the product hydroxymethanolate ion ( 18) by hydrogen molecule elimination as an additional source of the m/z 45 product ion.We find that the barrier for vicinal elimination of a hydrogen molecule leading to the formation of formate anion (15) is more than 50 kJ mol À1 higher (Fig. 10) than the barrier for formaldehyde elimination giving COOH À (11) (Fig. 6).The barrier for geminal elimination leading to the formation of hydroxycarbonyl anion is even higher.Together with the results of the calculations on the consecutive loss of a hydrogen molecule followed by the carbon monoxide elimination mentioned in the previous section, we consider formaldehyde elimination as the major source of the m/z 45 product ion and that the product ion has the hydroxycarbonyl ion (11)  geometry rather than that of the formate ion (15).

Second CO loss
We propose that the m/z 47 product ion, hydroxymethanolate, may dissociate further and give rise to the peak at m/z 19.The loss of 28 mass units corresponds to a second consecutive CO elimination; the elemental composition of the product ion is H 3 O À .
The mechanism for this process, based on the results of our computations, consists of 3 steps.First, the elongation of the C-O bond to the hydroxyl group leads to the transfer of the hydroxyl group towards the hydrogen atoms of the hydroxymethyl group.In concert with this, a proton is transferred from the emerging formaldehyde moiety to the emerging hydroxide.Second, the resulting intermediate (20) is then subject to conformeric rotation around one of the O-H bonds of the water molecule.The rearrangement in the last of the three steps can be understood as hydride transfer from the formyl ion to the water molecule, which results in elimination of CO and formation of H 3 O À (23) (Fig. 11).These reaction steps are essentially identical to those occurring in the direct reaction between OH À and CH 2 O, the Nibbering reaction, eqn (1).

Dissociation of H 3 O À
Paulson and Henchman reported the first observation of the H 3 O À ion and also estimated the ion binding energies with respect to two different dissociation reactions 44 -the one which corresponds to H À + H 2 O products was found to be 72 AE 5 kJ mol À1 , and the other which corresponds to OH À + H 2 products was found to be 29 AE 5 kJ mol À1 .Later it was shown by de Lange and Nibbering 46 that upon collisional activation the dissociation follows only the lowest energy pathway and OH À is the only product ion of H 3 O À collision-induced dissociation.A number of ab initio calculations by Chałasins ´ki, Kendall, and Simons, 54 by Ortiz, 55 and by Xantheas and Dunning, 56 and also more recent measurements by Miller et al. 57 are in a good agreement with the results quoted above.
On the basis of the experimental observations and quantum chemical computational modelling, we present a consistent model for the unimolecular decomposition of glycolate anions (m/z 75).The lowest energy processes and estimated minimum energy requirements are: While H 2 elimination (eqn ( 4)) has the lowest energetic threshold, the more favourable transition states for the loss of CO and in particular for the loss of CH 2 O make the latter two processes dominate at the energies probed in our experiments.In order to check the consistency of our reaction model against experimental observations, we simulated the reaction kinetics using the quantum chemical data (without modifications) at various pressures.As evident from Fig. 12, the kinetic simulations reproduce the experimental results well by showing that H 2 loss dominates at low internal energies and CO 2 becomes the most important reaction at high internal energies.In the broad intermediate energy range CO loss and CH 2 O loss predominate.
In a homogenous gas at low pressure and temperature, typical conditions found in the interstellar medium, only barrierless bimolecular reactions will be of interest in accounting for the synthesis of larger molecules from smaller.Herbst 16 considered the presence of anionic molecules in the interstellar medium based on radiative attachment of electrons, and estimated that their maximum abundance would be of the order of 1% of the neutral species.Noticeable cross sections for the production of H À , and somewhat lower for O À and OH À , have been measured for electron collisions with water in the gas phase. 582 O + e À -H À + OH, (10)   This leads us to propose the following series of reactions in an atmosphere of CO, H 2 O, CH 2 O, etc., based on the reverse of CO loss observed for glycolate: Although this scenario is reasonable, the significant barriers of reaction 15 (Fig. 10) has the consequence that the overall reaction is inefficient at the low temperatures typical of giant molecular clouds.Alternatively, we consider a route corresponding to the reverse of the reaction described in Fig. 5, i.e.

Conclusions
The unimolecular dissociation characteristics of the glycolate ion have been probed by mass spectrometry, and simulated by kinetic master equation model calculations in conjunction with a detailed survey of the relevant parts of the potential energy surfaces using high-level quantum chemical methods.Good consistency between the model calculations and the experimental observations is reported.The most significant reactions are loss of CO and loss of CH 2 O. Loss of CO 2 and H 2 , respectively, is also reported but are of minor importance.Interestingly, it was found that there is direct competition between the methoxy anion/carbon dioxide and the formate anion/formaldehyde product pairs in the sense that both product pairs are accessible from the same transition state so that the product distribution is determined by the detailed dynamics in the post-TS region of the potential energy surface.In addition, the CO loss may be followed by a second CO loss at Fig. 12 Calculated branching ratios for the first generation fragmentation reactions calculated using a master equation model based on a simplified potential energy surface shown in the ESI.† The glycolate anion was allowed to dissociate directly to give H 2 COH À plus CO 2 or to isomerize to 5, which could in turn dissociate to 11 and 12, 13 and 14 or isomerize to 17 and then dissociate to higher energies leading to H 3 O À and subsequently to OH À .The reverse of the unimolecular dissociation reactions with regard to possible synthetic pathways to glyoxylic acid in the gas phase via anionic mechanisms is discussed briefly.While such reaction pathways from simpler molecules appear likely under ambient conditions, it is concluded that under the conditions typically found in giant molecular clouds formation of glycolate by these reactions is not likely.

Fig. 1
Fig. 1 shows two CID mass spectra of deprotonated glycolic acid, recorded under different experimental conditions.Five fragment peaks are detected and the relative intensities depend on the gas pressure and the collision energy: (elemental composition of the neutral loss indicated)-at m/z 73 (-H 2 ), m/z 47 (-CO), m/z 45 (-COH 2 ), m/z 31 (-CO 2 ), m/z 19 (-C 2 O 2 ) and m/z 17 (-C 2 O 2 H 2 ).The H 2 loss results from a 1,2-elimination mechanism, as reported by Baker and Gabryelski on the basis of the mass shifts observed in their isotopic labelling experiments.19The fact that H 2 loss is only evident in our CID mass spectra recorded at relatively high pressure and low collision energy is indicative of a process requiring multiple low energy activating collisions, i.e. slow heating of the ions.This is also consistent with the fact that the peak resulting from H 2 loss is the most abundant one in the mass spectra of Baker and Gabryelski, while it is not seen in the spectra reported by Bialecki et al.20 In the former case, the experiments were conducted using an ion trap under slow heating conditions, while in the latter case the experiments were done using a triple quadrupole mass spectrometer under near single collision conditions at relatively high collision energy (E lab = 20 eV, Ar collision gas).In addition to H 2 loss, CO loss (giving rise to m/z 47) is the only other major fragmentation that can be inferred from the CID spectra from the ion trap experiment, indicating also a relatively low energy threshold for this process.In the spectra obtained under higher
3 O 2 À , which is identical to the m/z 45 ion formed by CO loss from the glycolate anion; a fact that both hints at the origin of the m/z 19 seen in our experiments and the mode of formation.It was shown by de Lange and Nibbering that collision-induced dissociation of H 3 O À ion results exclusively in formation of OH À by elimination of H 2 , 46 thereby also connecting the OH À (m/z 17) observed in our experiments to the same reaction chain.

Fig. 2
Fig.2Energy-resolved CID cross section (arbitrary intensity units) for the various dissociation routes.The experimental data were obtained under near single collision conditions with p = 6.4 Â 10 À5 mbar.

Fig. 3
Fig.3The lowest energy forms of the glycolate anion, with relative energies in kJ mol À1 given in italics.

Fig. 4
Fig. 4 Potential energy diagram for the isomerization reactions prior to the dissociation processes.

Fig. 5
Fig. 5 Potential energy diagram illustrating formaldehyde and hydrogen molecule loss from glycolate anions.

Fig. 6
Fig. 6 Potential energy diagram of the interconversion of formate and hydroxycarbonyl ions.

Fig. 7
Fig. 7 Potential energy diagram of the direct loss of CO 2 from 1 and the isomerization to CH 3 O À .

Fig. 8 .Fig. 9
Fig. 8 Potential energy diagram of the process leading to both CH 3 O À and HCO 2 À .

Fig. 11
Fig. 11 Potential energy diagram illustrating carbon monoxide loss from the hydroxymethanolate anion.

Table 1
Number of different reactive trajectories out of 50 for each temperature c the Owner Societies 2013