Reaction of atomic hydrogen with formic acid †

We study the reaction of atomic hydrogen with formic acid and characterize the radical products using IR spectroscopy in a Kr matrix and quantum chemical calculations. The reaction first leads to the formation of an intermediate radical trans -H 2 COOH, which converts to the more stable radical trans – cis -HC(OH) 2 via hydrogen atom tunneling on a timescale of hours at 4.3 K. These open-shell species are observed for the first time as well as a reaction between atomic hydrogen and formic acid. The structural assignment is aided by extensive deuteration experiments and ab initio calculations at the UMP2 and UCCSD(T) levels of theory. The simplest geminal diol radical trans – cis -HC(OH) 2 identified in the present work as the final product of the reaction should be very reactive, and further reaction channels are of particular interest. These reactions and species may constitute new channels for the initiation and propagation of more complex organic species in the interstellar clouds.


Introduction
Reactions of atomic hydrogen with various molecules are fundamental in combustion, atmospheric, and interstellar chemistry as well as in biology.2][3][4][5][6][7] Reactions of atomic hydrogen with organic compounds, especially with hydrocarbons, have been extensively studied; [1][2][3] however, relatively little information is available on reactions with organic acids.A very early report on the hydrogen abstraction reaction in acetic acid (H + CH 3 COOH -CH 2 COOH + H 2 ) in aqueous solution seems to be the only experimental example of a reaction between an organic acid and atomic hydrogen. 8ormic acid (FA) is the simplest organic acid that is an important intermediate in chemical synthesis and significant in atmospheric and interstellar chemistry as well as in human metabolism.Reactions of FA with OH and Cl radicals have been studied, which are important sources of the hydrocarboxyl radical (HOCO) in the atmosphere. 9The reaction of Cl with FA was found to proceed predominantly via abstraction of the hydrogen atom from the carbon atom to form HOCO rather than the abstraction of the hydroxyl hydrogen to form HCO 2 . 10,11FA is a model system of conformational isomerism resulting from different orientations of the OH group (trans and cis), 12 which is typical for carboxylic acids, including amino acids.3][14] Hydrogen atom tunneling has also been recently reported for the simplest aromatic carboxylic acid, benzoic acid. 15he conformation-dependent reaction between FA and atomic oxygen was observed in Kr and Xe matrices. 16The FA + O reaction leads to peroxyformic acid (HCOOOH) for the ground-state trans-FA conformer, and it results in the hydrogen bonded FAÁ Á ÁO complex for the higher-energy cis-FA conformer. 16To the best of our knowledge, no direct experimental and theoretical data are available for the reaction of FA with atomic hydrogen to date.It should be emphasized that the present work investigates the reaction of a neutral hydrogen atom with FA forming radicals.This reaction is different from FA protonation, which is a known process leading to closed-shell species.
Reactions of atomic hydrogen with organic compounds can produce radicals; for example, the hydrogen atom abstraction and molecular hydrogen formation are observed for ethane whereas the hydrogen atom addition to the unsaturated carbon atom occurs for ethylene resulting in the formation of a vibrationally excited C 2 H 5 radical. 1,2Hydrogen atom transfer plays a significant role in many organic and biological reactions.In particular, the radical isomerization via intramolecular hydrogen atom transfer in species with the C-H, O-H, and CQO bonds often occurs in peptide and protein radicals as a result of radiation or oxidative damage. 17,18his journal is © the Owner Societies 2014 It is probable that reaction of atomic hydrogen with FA produces radicals via hydrogen atom addition, which may then isomerize via intramolecular hydrogen transfer.
Matrix-isolation IR spectroscopy is a valuable method to study reactive intermediates obtained in reactions of atomic hydrogen.In our group, the H + C 2 H 2 -C 2 H 3 , 19 H + HCN -H 2 CN, 20 H + HNCO -H 2 NCO, 21 and H + N 2 O -HN 2 O 22 reactions have been detected in rare-gas matrices, in which the hydrogen atoms are produced by photolysis of suitable precursors (e.g., C 2 H 2 , HCN, HNCO, and HBr) and the formed radicals are characterized by IR spectroscopy.An interesting case is the H + SO 2 reaction. 23The higher-energy HSO 2 isomer is formed in this reaction instead of the lower-energy isomer cis-HOSO, which is explained by a low barrier for the HSO 2 formation and a high barrier for the cis-HOSO formation.HSO 2 can be isomerised into cis-HOSO by visible light, one proposed mechanism of which is the direct isomerisation via intramolecular hydrogen atom transfer.
In the present work, we study the reactions of hydrogen atoms with FA in a Kr matrix and characterize the radical products using IR spectroscopy.Hydrogen (deuterium) atoms are produced by UV photolysis of HBr and HCl (D 2 C 2 ).Thermal mobilization of the hydrogen atoms promotes the H + FA reaction, leading first to the formation of an intermediate radical trans-H 2 COOH, which converts via hydrogen atom tunneling to the more stable radical trans-cis-HC(OH) 2 .The structural assignment is aided by extensive deuteration experiments and ab initio calculations at the UMP2 and UCCSD(T) levels of theory.

Computational details and results
The optimization of the geometrical structures and calculation of the relative energies and vibrational spectra were performed using the UCCSD(T) and UMP2 methods with the aug-cc-pVTZ basis set.The Molpro program was used for the UCCSD(T) calculations and Gaussian09 (version C.01) for the UMP2 calculations.The minima on the potential energy surface (PES) were verified on the basis of harmonic vibration analysis which yielded no imaginary frequencies.In the case of the transition state (TS), one imaginary frequency was obtained.The minima on the PES associated with the TS were verified by means of the intrinsic reaction coordinate (IRC) analysis.The relative energy (DE) between the optimized geometrical structures was corrected for the zero-point vibrational energy difference (DZPVE).
The reaction of a hydrogen atom with FA may proceed by three different channels: (i) reaction with the oxygen atom of the hydroxyl group (OH), (ii) reaction with the oxygen atom of the carbonyl group (CQO), and (iii) reaction with the carbon atom.These channels have been studied by the UMP2/aug-cc-pVTZ and UCCSD(T)/aug-cc-pVTZ calculations.For channel (i), no minimum was found on the PES.This process results in the formation of a new H-O bond, the breaking of the C-O bond and the formation of a water molecule, which finally leads to the HCOÁ Á ÁH 2 O complex.Two other channels produce radicals shown in Fig. 1.For channel (ii), the HC(OH) 2 radical is formed that accepts two conformers.The two OÀH bonds can be oriented in the same direction as the CÀH bond (cis-cis form) or one OÀH bond can be oriented in the opposite direction to the CÀH bond (trans-cis form).These structures present the simplest geminal diol radicals.The trans-trans form of HC(OH) 2 is a transition state with one imaginary frequency, which is above the trans-cis form by about 1.55 kcal mol À1 (UMP2).Channel (iii) leads to the H 2 COOH radical with two CÀH bonds and one OÀH bond which may be oriented in two directions as shown in Fig. 1.The trans conformer has an out-of-plane OH bond (+OÀCÀOÀH = 561) and two non-equivalent CH bonds.The analysis shows that the cis form has practically no stabilization barrier; thus, it can be excluded from the consideration.The optimized geometries of the species under consideration are presented in Table S1, ESI.† The energetic effect associated with the formation of the new bonds, HÀO in HC(OH) 2 and HÀC in H 2 COOH, is characterized by the reaction energy (DE r ) calculated for the reactions: H + HCOOH -HC(OH) 2 (1) The value of DE r is the difference between the total energy of the HC(OH) 2 or H 2 COOH radical and the sum of the total energies of the hydrogen atom and FA (the lowest energy form) corrected by DZPVE (see Table 1).The formation of the HC(OH) 2 radical is energetically favourable for both UMP2     1.The isomerization from trans-cis-HC(OH) 2 to cis-cis-HC(OH) 2 requires an activation energy (E a ) of 2.57 kcal mol À1 (UMP2), which is several times smaller than the activation energy of the trans-to-cis reaction of FA (11.7 kcal mol À1 ).The optimized geometrical structure of the transition state (TS1) is shown in Fig. 1.The conformational change from trans-H 2 COOH to transcis-HC(OH) 2 requires the transfer of the hydrogen atom from the carbon atom to the carbonyl oxygen.The energy barrier associated with this process (TS2, Fig. 1) is 21.31 kcal mol À1 (UMP2).The energy barrier of the opposite process is substantially higher and equals 35.73 kcal mol À1 (UMP2).The change in the total energy for hydrogen atom transfer from trans-H 2 COOH to trans-cis-HC(OH) 2 along the IRC path is shown in Fig. S1 (ESI †).
The characteristic frequencies of the proposed radicals and their deuterated analogues obtained at the UCCSD(T)/aug-cc-pVTZ level of theory are presented in Table 2.The full spectra at the UMP2/aug-cc-pVTZ levels are presented in Table S2, ESI.†
The vibrational spectra of the ground-state (trans) conformer of HCOOH, DCOOH, and HCOOD in a Kr matrix recorded in our experiments agree well with the literature. 24,25193 nm photolysis (1500 pulses) of a HCOOH/HBr/Kr (1/2/1000) matrix decomposes about 85% HBr, producing H and Br atoms.About 30% FA is also decomposed producing mainly the COÁ Á ÁH 2 O complex. 24nnealing of a photolyzed HCOOH/HBr/Kr matrix at 31 K mobilizes hydrogen atoms in a Kr matrix, 26 which leads to the formation of new absorption bands (trace 1 in Fig. 2).Annealing of the photolyzed matrices also leads to a loss of FA.Fig. 3 shows that the amount of the products after annealing increases with the corresponding loss of FA.
The new absorption bands can be separated into two products marked by I U and I S (labels U and S indicate unstable and stable, respectively).Immediately after annealing, product I U dominates over product I S .The bands of product I U slowly decrease in intensity in the dark whereas the bands of product I S rise (trace 1 in Fig. 4).This process is accelerated by broad band IR light of the spectrometer (trace 2 in Fig. 2).In both cases, the amounts of I U and I S change with the same rates.The amount of I U goes to zero and the amount of I S tends to saturate.These facts feature the I U -to-I S conversion.This conversion process is also accelerated by vibrational excitation of product I U at 3614 cm À1 by the OPO (trace 3 in Fig. 2) whereas the vibrational excitation of product I S (at 3636.0 and 3598.4 cm À1 ) does not lead to the opposite process.The same two products are observed after photolysis and annealing of HCOOH/HCl/Kr matrices (trace 4 in Fig. 2).
It is noted that several bands of the products are split, which is probably due to trapping of the molecules in different matrix sites.The sub-bands of product I U are bleached with different efficiencies when exposed to narrow-band IR light.For excitation at about 3614 cm À1 , the band at 3613.7 cm À1 is efficiently bleached together with the bands at 1116.0 and 962.0 cm À1 .In contrast, the decrease of the bands at 3615.0, 1115.1, and 964.5 cm À1 is much less efficient for excitation at about 3615 cm À1 .
The experiments with deuterated species were also performed under the same experimental conditions.Annealing of a photolyzed HCOOD/HBr/Kr matrix initially leads to product II U that then converts in the dark to two products II S and II S 0 (trace 2 in Fig. 4, see discussion below).It should be noted that FA in this matrix is, due to incomplete deuteration, a mixture of HCOOD and HCOOH (70/30), which explains the formation of some amount of products I U and I S .In trace 2, the bands of these products are subtracted.Most indicative, the new characteristic bands are observed in the OD stretching region (2690-2650 cm À1 ), which are absent in the case of HCOOH/HBr/Kr matrices (trace 1 in Fig. 4).Additional new strong bands are observed at 1180-1140 cm À1 and 901 cm À1 .Annealing of a photolyzed DCOOH/HBr/Kr matrix leads to product III U which converts in the dark to III S (trace 3 in Fig. 4).No bands appear in this case in the OD stretching region whereas the CD stretching bands are observed for both stable and unstable products (2189.2and 1952.0 cm À1 , respectively).The most characteristic bands are observed at 1288/1285 cm À1 for III S and at 1093/1091 cm À1 for III U .The same characteristic features are detected in a HCOOH/C 2 D 2 /Kr matrix, however, in quite small amounts presumably due to the relatively low efficiency of the C 2 D 2 photolysis.The formation of DKrCCD (920 cm À1 ) 27 confirms the presence and mobility of D atoms in the matrix.In the HCOOH/C 2 D 2 /Kr case, the products formed in HCOOH/ HBr/Kr matrices are also observed and subtracted in trace 4. The formation of these background products is explained by the presence of H atoms probably produced by photolysis of the COÁ Á ÁH 2 O complex, which is a photodecomposition product of FA. 24 The presence of hydrogen atoms in the matrix is evidenced by the formation of HDCCD radicals (ca.1220 and 795 cm À1 ) originated from the H + DCCD reaction. 19In the deuteration experiments, the band intensities of the products correlate well with the annealing-induced loss of FA.The decay of the initial products in the dark and the formation of the final products have similar rates independent of deuteration described above (characteristic time B7.5 hours at 4.3 K, Fig. 5).The conversion is accelerated at elevated temperatures (characteristic time B5.0 hours at 27 K).The sum of the relative intensities of the decreasing and rising bands does not change in time.The situation is absolutely different for a DCOOH/C 2 D 2 /Kr matrix.In this case, an annealing-induced band is observed at 3616.6 cm À1 , which is the OH stretching mode of the ''unstable'' product; however, this band does not decrease in the dark even in 45 hours.The decay of this band can be promoted by broadband light of the spectrometer.
The weakness of the bands in this experiment prevents a detailed analysis in other spectral regions.

Spectral assignments
The new products appear as a result of annealing at temperatures that mobilize hydrogen atoms in a Kr matrix. 26It follows that the initial (unstable) product most probably originates from reactions of mobile hydrogen atoms with another fragment.For the assignment, several possibilities have been considered.First of all, the COÁ Á ÁH 2 O complex is present in the matrix as a result of the decomposition of FA. 24 Mobile hydrogen atoms can react with this complex producing the HCOÁ Á ÁH 2 O complex, and this possibility should be discussed.In fact, the HCOÁ Á ÁH 2 O complex obtained from the H + COÁ Á ÁH 2 O reaction has been recently identified by us, 28 and it is also observed in the present work.The formation of the HCOÁ Á ÁH 2 O complex is the most efficient when FA is mostly decomposed, i.e. when the amount of the COÁ Á ÁH 2 O complex is large.As expected, the bands of the HCOÁ Á ÁH 2 O complex do not change in the dark, in contrast to the new products studied in the present work.The spectra of these products differ significantly from that of the HCOÁ Á ÁH 2 O complex.For example, the HCOÁ Á ÁH 2 O complex in a Kr matrix has a strong CO stretching band at 1853.3 cm À1 . 28he correlation between the amount of the products and the loss of FA (Fig. 3) strongly suggests that the new species originate from the H + FA reaction.The small deviation from the exact proportionality can be caused, for example, by the annealinginduced formation of FA dimers and by the consumption of hydrogen atoms in reactions with HBr.In principle, the H + FA reaction can lead to the formation of a H 2 molecule and a radical species via H abstraction.However, the possible radicals cannot simultaneously have the OH and CH stretching modes, which rules out this scenario.
The addition of a hydrogen atom to different positions in FA should also be considered.For example, a hydrogen atom may react with the oxygen atom of the OH group forming HCOOH 2 .Our calculations suggest that this species is not an energy minimum, and it spontaneously decomposes to a pair of H 2 O and HCO.Thus, this reaction (if occurred) would lead to the formation of the HCOÁ Á ÁH 2 O complex in a matrix.
The addition of a hydrogen atom to the oxygen atom of the CQO group is also possible, which would result in the HC(OH) 2 radical.The calculations show that the HC(OH) 2 radical has two local minima differing by the orientation of the OH groups, the higher-energy cis-cis form and the lower-energy trans-cis form (Fig. 1).Similar conformers exist for carbonic acid, although the energetics is opposite, with the cis-cis form being more stable than the trans-cis form (by 1.0-1.6 kcal mol À1 ). 29,30Both of these conformers of carbonic acid have been observed in an Ar matrix. 29For this reason, we have to consider the initial formation of cis-cis-HC(OH) 2 which converts to trans-cis-HC(OH) 2 .However, the calculated barrier between these two species is only 0.44 kcal mol À1 , which is much smaller than  that of carbonic acid (9.2 kcal mol À1 ) 30 and definitely too small to stabilize the higher-energy conformer at the timescale of hours.For comparison, the cis form of FA converts to trans-FA much faster than in the present study despite the calculated stabilization barrier of 7.7 kcal mol À1 .The exclusion of the assignment of I U to cis-cis-HC(OH) 2 is confirmed by the deuteration experiments.In particular, no OD stretching band is seen as a result of the D + HCOOH reaction (see trace 4 in Fig. 4).
After these exclusions, the most probable mechanism of the H + FA reaction is the addition of the hydrogen atom to the carbon atom, which results in the H 2 COOH radical.For this radical, only one true minimum, the trans conformer, is found, and its total energy is higher than that of trans-cis-HC(OH) 2 by 8.10 kcal mol À1 (UCCSD(T)).Thus, we assign the decreasing bands to the trans-H 2 COOH radical and the rising bands to the trans-cis-HC(OH) 2 radical.Our spectroscopic data show that the trans-H 2 COOH radical dominates immediately after annealing (trace 1 in Fig. 2 and Fig. 3 and 5) and then it slowly converts to trans-cis-HC(OH) 2 .
The amount of hydrogen atoms in the 193 nm irradiated HCOOH/HBr/Kr (1/2/1000) matrix is roughly twice as much as the amount of remaining FA molecules.It follows that the initial product of the H + FA reaction (trans-H 2 COOH) may in principle react with a second hydrogen atom to form the closedshell species CH 2 (OH) 2 (methanediol).This molecule has been observed in an Ar matrix and its vibrational spectrum is different from the spectra of our species. 31For example, the strong CO stretching band of CH 2 (OH) 2 is not detected in the present experiments.Now, we analyse the spectra of the reaction products (Table 2).The OH stretching region is very informative.In theory, the trans-cis-HC(OH) 2 radical has two OH stretching bands of similar intensity separated by about 33 cm À1 whereas the trans-H 2 COOH radical has one intense band located between the trans-cis-HC(OH) 2 bands.These calculations are in perfect agreement with the experiment.In the H + HCOOH experiment, product I S has two bands separated by 39 cm À1 and product I U has one OH stretching band located in between.In contrast, the calculated cis-cis-HC(OH) 2 radical has the OH stretching band higher in frequency than the two bands of trans-cis-HC(OH) 2 , which also supports our assignment of the initial product to trans-H 2 COOH.
For the H + HCOOD reaction, the OH stretching region exhibits similar rising bands as in the H + HCOOH case whereas no decreasing band appears in between.This fact strongly confirms that the unstable product is H 2 COOD that has no OH stretching mode but only OD.The two rising bands in the OD stretching region are separated by 26 cm À1 and the decreasing band appears in between (trace 2 in Fig. 4b).The appearance of four (not two) rising bands is explained by the fact that the trans-cis-HC(OH)(OD) radical has two isotopic analogues.In this case, the bands at 3635 and 2658/2654 cm À1 are assigned to the species with the cis-OH and trans-OD groups.The second isotopic analogue has the trans-OH and cis-OD groups, absorbing at 3598 and 2686/2684 cm À1 , respectively.
For the H + DCOOH and D + HCOOH reactions, the OH stretching bands (traces 3 and 4 in Fig. 4) are similar to the H + HCOOH case for both initial and final products (trace 1).An important observation is the absence of the OD stretching bands in these two cases, which indicates that the H (not D) atom is involved in the conversion process (see also later).For the D + HCOOH reaction, the lack of the OD stretching absorption provides an additional support that the initial product is formed by the addition of a D atom to the carbon atom of FA rather than to the oxygen atom.
In the nCH and nCD regions, weak but characteristic bands are observed.For the H + HCOOH and H + HCOOD reactions, both products absorb in the CH stretching region (3100-2800 cm À1 ) whereas the CD stretching region (2300-2000 cm À1 ) shows no bands.For the H + DCOOH reaction, the initial product exhibits the CH and CD stretching bands at 2872.0 and 1952.0 cm À1 , which further confirm the addition of a hydrogen atom to the carbon atom of FA.
The spectrum in the deformation region (1400-900 cm À1 ) also supports the proposed assignment.For the H + HCOOH reaction, the positions of the two decaying (B1116.0/1115.0 and 964.5/962.0cm À1 ) and three rising bands (B1326, 1186.0, and B1129 cm À1 ) are in good agreement with the calculations (1113.7 and 1002.2 cm À1 for trans-H 2 COOH; 1368.3, 1205.3, and 1155.1 cm À1 for trans-cis-HC(OH) 2 ).For the H + HCOOD reaction, the spectrum significantly changes in qualitative agreement with the calculations.An intense band at 1288 cm À1 is observed for the final product for the H + DCOOH and D + HCOOH reactions but it is absent in other cases.According to the calculations, this band is characteristic of trans-cis-DC(OH) 2 (calculated frequency 1329.7 cm À1 ).It should be admitted that some of the experimental shifts in this region are not accurately predicted.An example is the order of two rising bands and one decaying band observed in the H + HCOOD reaction at around 1150 cm À1 .In the experiment, the decaying band is between the two rising bands whereas the calculations predict it at a higher frequency.We connect this small discrepancy with the difficulties to describe open-shell species even at the UCCSD(T) level of theory.Very accurate calculations of vibrational spectra of open-shell species require special approaches, 32 but this exceeds the scope of the present work.Furthermore, this is not a surprise that some band intensities are not accurately predicted by the MP2 calculations.

Formation of the radicals and intramolecular hydrogen atom transfer
As shown above, new radicals trans-H 2 COOH and trans-cis-HC(OH) 2 are formed in the H + FA reaction.Our spectroscopic data in a Kr matrix suggest that this reaction primarily produces the trans-H 2 COOH radical by the addition of a hydrogen atom to the carbon atom of FA and this species dominates immediately after annealing.The formation of the trans-H 2 COOH radical in the H + FA reaction is energetically favourable (by 1.91 kcal mol À1 ) at the UCCSD(T) level of theory.This radical is higher in energy than the lowest-energy trans-cis-HC(OH) 2 species by 8.10 kcal mol À1 (UCCSD(T)).The situation when a reaction produces a metastable isomer is not unique to matrix-isolation studies.One example is the efficient formation of the higher-energy dimers of FA after deposition of matrices under some conditions. 33As another example, the H + SO 2 reaction mainly leads to a higher-energy isomer HSO 2 but not to the ground-state isomer HOSO. 23he experimental results show that the initial product trans-H 2 COOH converts to the trans-cis-HC(OH) 2 radical in the dark at cryogenic temperature (4.3 K).We explain this isomerization process by intramolecular hydrogen-atom tunneling from carbon to oxygen atom.This process agrees with the calculated energetics of the system showing that the trans-H 2 COOH radical is higher in energy than the trans-cis-HC(OH) 2 .The trans-H 2 COOH stabilization barrier is high enough (21.31 kcal mol À1 at the UMP2 level of theory) to prevent the over-barrier process at low temperatures; however, quantum tunneling of a hydrogen atom is possible through this barrier.The tunneling mechanism is supported by the strong H/D isotope effect: (i) no OD stretching bands appear as a result of the conversion of HDCOOH and (ii) D 2 COOH is practically stable.The strong H/D isotope effect has been repeatedly shown for quantum tunneling. 34,35he cis-cis-HC(OH) 2 radical can be formed in this tunneling process but quickly decays to the trans-cis conformer due to the very low stabilization barrier (0.44 kcal mol À1 ).It has been recently reported that methylhydroxycarbene (H 3 C-C-OH) isomerizes to acetaldehyde (H 3 C-CHQO) with a lifetime ca. 1 hour in an Ar matrix at 11 K, and the process was explained by a facile hydrogen tunneling through a barrier of 28.0 kcal mol À1 . 36he isomerization from trans-H 2 COOH to trans-cis-HC(OH) 2 is accelerated by vibrational excitation of nOH fundamental by the narrow-band IR light but the opposite process does not occur.This fact fully agrees with the computational energies showing a significantly higher isomerization barrier (35.73 kcal mol À1 ) for the opposite process.It should be noted that the photon energy (3536 cm À1 or 10.1 kcal mol À1 ) is smaller than the calculated stabilization barrier of trans-H 2 COOH; however, light-induced under-barrier tunneling was demonstrated. 37The conversion process is accelerated by broad-band IR light even for D 2 COOH, which is probably contributed by an over-barrier isomerization.Furthermore, this hydrogen-atom tunneling is consistent with the acceleration of the conversion process at elevated temperatures (ca.7.5 hours at 4 K and 5.0 hours at 27 K), which is typical for quantum tunneling. 34,35adical isomerization via hydrogen atom tunneling can compete with radical decomposition.In the present case, the possible decomposition channel would lead to H 2 + COOH (or HOCO).However, no evidence of these species is provided by the experiments as mentioned above.Another remarkable fact is that the present situation is quite different from the Cl + FA reaction leading to HCl. 9 It can also be noted that the weak trans-cis-HC(OH) 2 bands are seen immediately after annealing (trace 1 in Fig. 2).This observation seemingly suggests that some amount of this species is formed directly in the H + FA reaction.However, more probably, these weak bands appear due to the isomerization during the annealing and the first measurement.

Conclusions
This work presents the first report on the reaction of atomic hydrogen with formic acid.The experiments in a Kr matrix show that this reaction proceeds via addition of the hydrogen atom to the carbon atom of FA to form the trans-H 2 COOH radical as the initial product.The trans-H 2 COOH radical isomerizes via intramolecular hydrogen atom tunneling to the lower-energy trans-cis-HC(OH) 2 radical, which is the simplest geminal diol radical.These two radicals are directly observed for the first time and characterized by using IR spectroscopy and ab initio calculations at the UMP2 and UCCST(T) levels of theory.The assignment of the products and the hydrogen atom tunneling mechanism of the isomerization are confirmed by experiments with deuterated species.
The simplest geminal diol radical trans-cis-HC(OH) 2 obtained and identified in the present work is presumably very reactive, and further reaction channels are of particular interest.For instance, these species may constitute new channels for the initiation and propagation of more complex organic species in the interstellar clouds. 38To our knowledge, the new reactions and species observed in the present work have not been considered in these interstellar processes.Another natural extension of this work is the study of reactions between atomic hydrogen and more complex molecules such as amino and fatty acids.

Fig. 1
Fig.1Schematic potential energy diagram of the radicals that may be formed in the H + HCOOH reaction.The energy values (in kcal mol À1 ) are calculated at the UMP2 level (for the UCCSD(T) energies of the optimized structures see Table1).The structures shown in rectangles are transition states (TS1 and TS2); cis-H 2 COOH (shown in a dashed circle) is energetically unfavourable.In trans-H 2 COOH, H1 and H4 are on the same side of the heavy-atom skeleton.

Fig. 3
Fig. 3 Amount of products I U and I S versus the corresponding loss of FA (3536.5 cm À1 ) in the H + FA reaction.The spectra were recorded at 4.3 K immediately after annealing at 31 K of a photolyzed HCOOH/HBr/Kr (1/2/ 1000) matrix.The data points correspond to the integrated intensities of the absorption bands in different experiments.The lines are linear fits.

Fig. 4 (
Fig. 4 (a and c) Calculated IR spectra of trans-H 2 COOH (negative bands) and trans-cis-HC(OH) 2 (positive bands) and their deuterated analogues.The dashed line in the calculated spectrum (2) corresponds to trans-cis-HC(OH)(OD) with the cis-OD group.The splitting of the negative band in the calculated spectrum (3) of trans-HDCOOH is caused by different locations of the CD group.(b and d) Difference IR spectra showing the conversion processes after different reactions: (1) H + HCOOH, (2) H + HCOOD, (3) H + DCOOH, and (4) D + HCOOH.Thick and thin lines in the experimental spectra correspond to 4 and 8 hours in the dark at 4.3 K after annealing, respectively.Bands of FA are marked by asterisks.The FA/HBr (C 2 D 2 )/Kr matrices were first photolyzed at 193 nm and then annealed for 5 min at 31 K.The spectra were recorded at 4.3 K.Note the different absorbance scales in the left and right panels.

Fig. 5
Fig. 5 Relative concentrations of the products as a function of time at different temperatures in the dark.(a) The decay of the initial products trans-H 2 COOH (3614 cm À1 ), trans-H 2 COOD (2666 cm À1 ), and trans-HDCOOH (3614 cm À1 ); (b) The formation of the final products transcis-HC(OH) 2 (3636 cm À1 ), trans-cis-HC(OH)(OD) (2684 cm À1 ), and transcis-DC(OH) 2 (3636 cm À1 ).The sum of the relative concentrations of the unstable and stable products is shown by stars in panel a.

Table
).The structures shown in rectangles are transition states (TS1 and TS2); cis-H 2 COOH (shown in a dashed circle) is energetically unfavourable.In trans-H 2 COOH, H1 and H4 are on the same side of the heavy-atom skeleton.

Table 1
Energetics (in kcal mol À1 ) for the radicals that may be formed in the H + HCOOH reactions calculated at the UCCSD(T) and UMP2 (in parentheses) levels of theory This journal is © the Owner Societies 2014

Table 2
Experimental (I U -III U and I S -III S ) and calculated characteristic bands of trans-H 2 COOH and trans-cis-HC(OH) 2 and their deuterated analogues a