Tunnelling in carbonic acid †

The cis,trans-conformer of carbonic acid (H2CO3), generated by nearinfrared radiation, undergoes an unreported quantum mechanical tunnelling rotamerization with half-lives in cryogenic matrices of 4–20 h, depending on temperature and host material. First-principles quantum chemistry at high levels of theory gives a tunnelling half-life of about 1 h, quite near those measured for the fastest rotamerizations.

Although carbonic acid (H 2 CO 3 , 1, Fig. 1) is constantly on the lips of millions in the form of carbonated beverages, evidence for the molecule in the gas phase was obtained only in 1987 by thermolysis of NH 4 HCO 3 and mass spectrometric detection. 1or many years this profoundly important molecule, which is considered essential for many biological processes, 2 escaped unequivocal identification in the gas phase.Indeed, it was vigorously debated whether H 2 CO 3 could persist long enough to be spectroscopically identifiable, as the molecule decomposes readily into CO 2 and H 2 O. 3 At the same time, evidence was found for the existence of 1 in atmospheric clouds both on earth and several terrestrial planets. 4Finally, in 2009 5 and 2011, 6 Fouriertransform microwave spectroscopy provided the first spectroscopic characterization of gaseous H 2 CO 3 , produced by passing a pulsed electric discharge of CO 2 in a supersonic jet through a water reservoir.The cis-trans (1ct) 5 and cis-cis (1cc) 6 rotamers (Fig. 1) were identified unequivocally, while the trans-trans (1tt) isomer was deemed unattainable because of its high relative energy and small rotational barrier (1.8 kcal mol À1 at CCSD(T)/CBS) 7 to 1ct.
[10] Shortly thereafter, the protonation of potassium bicarbonate with excess HCl in glassy methanolic solutions led to a solid, later termed the a-H 2 CO 3 polymorph. 11An alternative b-H 2 CO 3 polymorph was generated in an analogous fashion in aqueous glassy solutions in the presence of excess HBr, and conversion of b-H 2 CO 3 to a-H 2 CO 3 was found upon treatment with HCl in methanol. 12Overturning previous understandings, we proved recently 7 through comparison of matrix-isolation 13,14 infrared (IR) spectra of independently prepared compounds that the substance named a-H 2 CO 3 is in fact carbonic acid monomethyl ester, 15,16 while only the b-polymorph 17,18 is indeed solid H 2 CO 3 . 7,19n doing so, we developed a novel synthetic route to gas-phase H 2 CO 3 utilizing the ester pyrolysis of either Boc 2 O or di-t-butyl carbonate, a conceptually simpler approach 20 generating much higher concentrations in the gas phase than previously achieved.
Carbonic acid is generally considered a fleeting molecule at ambient temperatures due to rapid decomposition into CO 2 and H 2 O in the gas phase 3 (Fig. 1) or deprotonation in solution.However, the gas-phase decomposition of 1 has a large barrier in excess of 40 kcal mol À1 , 21,22 because it is a thermally forbidden [2 s +2 p ] cycloreversion.In the presence of water, an (auto)catalytic [23][24][25] process appears to be operative instead: computations reveal that complexation of 1 with explicit water molecules allows passage through more favourable cyclic transition states of six or more members that are lower in energy by 16-19 kcal mol À1 . 21Here we outline hitherto unreported conformational tunnelling of 1, thereby capturing one of its key reactivities.The detailed delineation of the conformational, reactivity, and spectroscopic landscape of 1 will aid in its identification in the atmosphere and extraterrestrial environments.
Our highly efficient gas-phase preparation of 1 via ester pyrolysis (Scheme 1) and its trapping in noble gas matrices under cryogenic conditions allows the unimolecular kinetic behaviour to be probed. 7Di-t-butyl carbonate (2) readily undergoes double ester fragmentation via high vacuum flash pyrolysis (HVFP) with concomitant loss of two isobutene equivalents.Immediately after passing 2 through the pyrolysis tube, the products are trapped in noble gas matrices, kept at low temperatures (3-20 K), and monitored by means of IR spectroscopy.Bimolecular reactions do not occur because the trapped molecules are present in low concentrations and are fully surrounded by noble gas atoms (ratio ca. 1 : 10 3 ); the solid noble gas matrices do not allow molecular diffusion.
Infrared bands for most vibrational modes of the two observable isomers of 1 were identified before. 7However, we were able to assign several new absorptions (Tables S16 and  S17 After enrichment of the 1ct conformer by narrow band NIR excitation (see ESI † for details) of matrices containing 1, we observe its decay in solid Ne (3 K), Ar (3, 12.5, 17.5, 20 K), Kr (12.5 K), and Xe (12.5 K) with effective half-lives (t eff ) of 4-20 h.During the kinetic measurements a long-pass IR filter cutting off all wavelengths below 4.5 mm (ca.2200 cm À1 ) was used to prevent unwanted photochemistry induced by the globar radiation of the IR spectrometer.Decomposition of 1 into CO 2 and water (Fig. 1A) can be excluded, as we monitor the concomitant disappearance of 1ct and appearance of 1cc.Since the activation barriers for the other two possible processes (Fig. 1B and C) are too large to be overcome at cryogenic temperatures, quantum mechanical tunnelling (QMT) 26 must be responsible for the disappearance of 1ct.The [1,3]H-shift of 1ct to 1cc (Fig. 1B) via TS2 exhibits a 34.6 kcal mol À1 barrier at the CCSD(T)/CBS level; this compares well with the QCISD(T)/6-311++G(d,p)//MP2(full)/ 6-31++G(d,p) barrier of 34.3 kcal mol À1 reported earlier. 27arriers of this magnitude are not too large in principle to prohibit H-tunnelling on the time scales observed here, but in this case the potential energy profile is not ''thin'' [28][29][30][31][32][33][34] enough to make hydrogen transfer viable (Fig. 1).Unsurprisingly, we computed a prohibitive half-life of t = 4.4 Â 10 6 y for this reaction along a CCSD(T)/cc-pVQZ//MP2/aug-cc-pVTZ minimum energy path within a Wenzel-Kramers-Brillouin (WKB) tunnelling model 35,36 that has proved very accurate for H-transfer reactions. 28,37he isomerization of 1ct through C-O single bond rotation has a barrier of only 7.9 kcal mol À1 (Fig. 1), which is in the 2-13 kcal mol À1 range typical of carboxylic acids [38][39][40][41][42][43][44] but still much too large to be overcome thermally at the temperatures of our experiments.9][40][41][42][43][44] We executed a detailed kinetic analysis of this tunnelling phenomenon in different noble gases and at various temperatures (Table 1).
The intricate kinetic behaviour of matrix-isolated 1 was treated within a multi-exponential model we recently developed to describe a similar conformational tunnelling process in oxalic acid, termed domino tunnelling. 45To achieve highquality fits of the data within the experimental error bars, the assumption of three classes of matrix-isolated carbonic acid molecules was required: a frozen class, in which tunnelling is entirely precluded by the environment; a slow class with seemingly hindered conformational tunnelling and half-lives (t 1 ) of 11-26 h; and a fast class, characterized by half-lives (t 2 ) of 2-5 h, for which the matrix atoms appear to interfere to only a limited extent with the conformational tunnelling.We emphasize here that our three-class model is phenomenological in nature, although it is appealing to interpret the experimental results in terms of a hindering matrix environment that softens at higher temperatures.The final kinetic analysis employed two half-lives (t 1 , t 2 ) and three initial population ratios (1ct fast : slow, 1ct frozen : slow, and total 1ct : 1cc) as fitting parameters in a simultaneous, nonlinear least-squares fit of all suitable infrared spectral bands of both 1ct and 1cc.Complete details of the kinetic analysis are given in the ESI † along with plots of the data and fits.Deuteration completely suppresses the tunnelling Scheme 1 High-vacuum flash pyrolysis of di-t-butyl carbonate (2) and trapping of 1 in noble gas matrices.
Table 1 Tunnelling half-lives t (in h) in solid noble gas matrices at several temperatures.X frozen , X fast , and X slow are the mole fractions for the frozen (t = N), slow (t 1 ), and fast (t 2 ) classes of trapped molecules of 1ct, respectively a Matrix (T/K) t isomerization, and we compute t = 2400 y for monodeuterated 1 (see ESI †); no tunnelling of deuterated 1 was observed over one day at 20 K in Ar.The Ar (20 K) results shown in Fig. 2 are representative of our kinetic analysis of conformational tunnelling.A single-exponential fit to the 1ct decay profile yields an effective half-life of 4.4 h, but the poor quality of the fit demonstrates the multi-exponential character of the kinetics.The initial ratio of fast : slow sites is near two, and in this case the data do not support the presence of frozen sites (X frozen = 0).The slow and fast tunnelling of 1ct exhibits the half-lives t 1 = 14 AE 3 h and t 2 = 3.0 AE 0.2 h, respectively.
To theoretically support our experiments, a distinguished reaction path (DRP) was generated for 1ct -1cc by a series of MP2/aug-cc-pVTZ constrained geometry optimizations, whereby the torsion angle t OQC-O-H was varied in 51 increments between 01 and 1801 while all other internal coordinates were optimized.Along this path, CCSD(T)/CBS energies were computed to obtain an accurate potential energy curve for the isomerization.Finally, a DRP projected frequency analysis 46 was employed to evaluate harmonic vibrational frequencies (o i ) and zero-point vibrational (ZPVE) corrections for each value of t OQC-O-H .Tunnelling probabilities (k) for 1ct -1cc were computed using exact numerical integration 45 and the WKB approximation. 35Tunnelling rates were evaluated as the product of the transmission probability k and the classical rate (o 0 ) at which the reactant hits the barrier, assuming that the collision energy under cryogenic conditions is the reaction-mode ZPVE of o 0 /2.The computed intrinsic tunnelling half-lives (t 0 ) for the 1ct -1cc isomerization in the gasphase amount to 56 min employing exact numerical integration and to 72 min within the WKB approximation.In support of the methodology underlying these results, a systematic computational study on acid tunnelling has found that ''the minimum energy path can be considered as a good approximation for the cis-to-trans tunneling path with maximum transparency''. 47lthough our initial tunnelling computations do not incorporate reaction-path curvature effects that would further enhance the rate, 48,49 the results compare well with the fast components (t 2 ) obtained from the multi-exponential kinetic analysis of our experiments.Importantly, the X fast mole fractions increase toward 1 and the t eff values diminish toward t 0 upon elevating the temperature of the Ar matrix (Table 1 and Fig. 3).While a rigorous assessment of various theoretical methods for computing tunnelling rates is beyond the scope of this communication, we report additional predictions of half-lives obtained with the ZCT and SCT approaches implemented in POLYRATE in the ESI.† [50][51][52] The complicated observed tunnelling behaviour is presumably attributable to varying degrees in which the matrix environment hinders the rotamerization of 1.The thermal energy of our highesttemperature experiments (20 K) is not sufficient to significantly excite the C-O torsional vibration (491 cm À1 , CCSD(T)/cc-pVTZ) to accelerate the tunnelling rotamerization, which would only become important at about 60 K; however, there is enough thermal energy to excite the solid noble gas host material, whose Debye frequencies are in the range of 40-70 cm À1 . 53As the matrix softens at higher temperatures, it can better accommodate the structural changes required for tunnelling of 1ct.Such temperature dependences of conformational tunnelling processes have already been observed experimentally 54 and rationalized through a theoretical study of the noble gas environment, 55 which showed that the temperature dependence of the tunnelling event is most likely attributable to a reorganization of the matrix environment.The circumstances are very different for the highly exothermic tunnelling of hydroxycarbenes to the corresponding aldehydes. 32,37he intricate dependence of QMT on the noble gas host is currently not fully understood (Fig. 3).Argon seems in this respect well behaved as tunnelling accelerates with increasing temperature in the 3-20 K temperature regime; this corresponds well with a temperature-dependent activation of the noble gas atoms during the QMT rotamerization.Krypton is comparable, while the most polarizable Xe-matrix impedes tunnelling.Comparison with the Ne-matrix is difficult, but its low polarizability is likely to reduce the interaction with the interconverting conformers relative to Ar at the same temperature.Such behaviour has been reported before, e.g., in the conformational tunnelling of formic acid that exhibits increasing tunnelling rates in going from Xe to Ne, while its O-deuterated isotopologue shows the opposite trend. 56n terms of competing tunnelling mechanisms, it is clear that while the two possible pathways from 1ct to 1cc via TS2 and TS3, respectively, have approximately the same width (Fig. 1), the tunnelling probabilities are primarily determined by barrier height, as found in classic kinetically controlled reactions.This situation is in marked contrast to the phenomenon of tunnelling control, 32 where barrier width trumps barrier height in favour of the thermodynamic product.The 26.7 kcal mol À1 activation barrier difference between the two transition structures translates into a tunnelling rate ratio of approximately 1 : 10 10 , i.e., essential quantitative selectivity.
We have unveiled a novel conformational tunnelling isomerization of 1 (in the range of 4-20 h) that is supported by high-level quantum chemical computations.The tunnelling process strongly depends on the temperature and nature of the matrix environment.The three phenomenological classes of tunnelling sites demonstrate that tunnelling is affected by the environment and therefore can in principle be controlled externally.On the basis of these findings and the theoretical predictions, we expect the tunnelling half-life for the isomerization of 1ct to 1cc would be even faster in the gas phase.
This work was supported by a collaboration grant from the U.S. National Science Foundation (CHE-1124885) and the Deutsche Forschungsgemeinschaft (Schr 597/18-1) to W. D. A. and P. R. S.; V. H. thanks the Liebig College at Justus-Liebig University, sponsored by the State of Hesse.

Fig. 1
Fig. 1 Potential energy curves at CCSD(T)/CBS for disappearance of 1ct by decomposition into H 2 O and CO 2 (A), [1,3]hydrogen transfer (B), or internal C-O bond rotation (C).The profiles are drawn in terms of computed arc length (in u 1/2 bohr) along the intrinsic reaction path.
, ESI †), which assist the detection of 1 in extraterrestrial and other environments.The additional assignments include the n s (O-C-O) symmetric stretching and d(O-C-O) deformation modes of both conformers, as well as the d(O-CQO) deformation of 1ct.With decreasing polarizability of the noble gas host from Xe to Ne, we usually observe blue-shifted spectral features that slowly approach the gas-phase values for the vibrational transitions (Fig. S1, ESI †).On this basis it was also possible to assign the very close lying d(H-O-C) + n as (O-C-O) vibrations of 1cc and 1ct (computed at 1183.8 and 1182.0 cm À1 , respectively) in neon in the theoretically predicted ordering.

Fig. 2
Fig. 2 Time evolution of H 2 CO 3 conformer concentrations in an Ar matrix at 20 K after NIR irradiation.Blue and green curves result from the global, nonlinear least-squares fit of the multi-exponential kinetic model to three and five IR bands of 1ct and 1cc, respectively.The dashed curve is the best single-exponential fit for 1ct.

Fig. 3
Fig. 3 Mole fractions of frozen, slow, and fast tunnelling sites in different noble gas matrices at several temperatures.