Not particularly strong programming at the year’s spring ACS meeting – but one great session in the organic division yesterday. This was the awards session in honor of John Baldwin getting the James Flack Norris Award for physical organic chemistry.

First to speak was James Duncan, who discussed his recent CASSCF computations looking for pseudopericylic [3,3]-sigmatropic migrations. I will be commenting on his latest work in a post that will appear soon.

I had to skip the next talk, but came back to hear John Brauman discuss recent work on the solvation effect in the S_N2 reaction. This is an interesting case of where the screening of larger substituents is counterbalanced by geometric changes that lead to greater charge distribution. The net effect is that they cancel each other out, and the methyl,ethyl, iso-propyl, butyl β-effect is negligible.

Next was Peter Schreiner who discussed his carbene work, specifically the enormous tunneling effect observed in hydroxymethylene (see this post). He discussed some new work, that is if anything even more fantastic on methylhydroxycarbene – look for this work perhaps later in 2011.

Last to speak was John Baldwin – and he described his truly tour de force efforts in examining the [1,3]-rearrangements of vinylcyclopropane and vinylcyclobutane. The former work is described in my book, while the later study is still ongoing.

John’s work is amazingly painstaking and careful. I am truly in awe of his dedication in taking on extremely difficult studies that require enormous care. John has really taught us a lot – not just about these rearrangements (they involve diradicals on a flat plateau demanding dynamic analysis – but how to think about a study and then carry it out to fruition so that all details are assessed. A truly deserving recipient!

Following up on his previous studies of isotope effects on the ring opening of cyclopropylcarbinyl radical 1 to give 2 (see my previous post), Borden now reports on its kinetic isotope effect (KIE).¹

Using the small-curvature tunneling approximation along with structures and frequencies computed at B3LYP/6-31G(d), he finds a negligible KIE at C₁, consistent with little motion of C₁ in the transition vector. The KIE for substitution at C₄ is large (k(¹²C/¹⁴C)=5.46), also consistent with its large motion in the transition vector. What is surprising is the KIE for deuterium substitution at C₁: 0.37. This is a large inverse isotope effect!

Analysis of the vibrational frequencies that involve the C₁ hydrogens provides an explanation. In going to the TS for the ring opening, both the torsional motion about the C₁-C₂ bond (making the double bond) and the pyramidal motion increase in frequency. This leads to a higher activation barrier for H than D, and the inverse isotope effect.

References

(1) Zhang, X.; Datta, A.; Hrovat, D. A.; Borden, W. T., “Calculations Predict a Large Inverse H/D Kinetic Isotope Effect on the Rate of Tunneling in the Ring Opening of Cyclopropylcarbinyl Radical,” J. Am. Chem. Soc., 2009, 131, 16002-16003, DOI: 10.1021/ja907406q.

Following on the great study of hydroxycarbene¹ (see my blog post), Schreiner now reports on the synthesis and characterization of dihydroxycarbene 1.² It is prepared by high-vacuum flash pyrolysis of oxalic acid (Scheme 1).

Scheme 1

Dihydroxycarbene can exist in three different conformations characterized by the relationship about the C-O bond, either s-cis or s-trans. The three conformations are shown in Figure 1, and the s-trans,s-trans structure is the local energy minimum (computed at CCSD(T)/cc-pVTZ).

1tt (0.0)

1ct (0.1)

1cc (6.7)

Figure 1. CCSD(T)/cc-pVTZ optimized geometries and relative energies (kcal mol^-1) of the conformers of 1.²

Identification of the 1 is made through comparison of the experimental and computed IR vibrational frequencies. As an example, the experimental and computed frequencies for the s-trans,s-trans conformer are listed in Table 1. The agreement is excellent.

Table 1. Computed and experimental vibrational frequencies (cm^-1) and intensities (in parentheses) of the s-trans,s-trans conformation of 1.²

Unlike hydroxycarbene, dihydroxycarbene is stable. The amazing instability of hydroxycarbene is due to tunneling through a large barrier: nearly 30 kcal mol^-1.¹ The tunneling route for the decomposition of 1 is more difficult for two reasons. First, its C-O bond is quite strong; the C-O distance is quite short, 1.325 Å. This makes a long distance that must be traversed in the tunneling mode. (The strong bond is due to π-donation from the oxygen lone pair into the empty carbon p orbital; this is noted by the large rotational barrier about the C-O bonds of 17 kcal mol^-1!) Second, the activation barrier for decomposition is very high, at least 34 kcal mol^-1.

References

(1) Schreiner, P. R.; Reisenauer, H. P.; Pickard Iv, F. C.; Simmonett, A. C.; Allen, W. D.; Matyus, E.; Csaszar, A. G., "Capture of hydroxymethylene and its fast disappearance through tunnelling," Nature, 2008, 453, 906-909, DOI: 10.1038/nature07010.

(2) Schreiner, P. R.; Reisenauer, H. P., "Spectroscopic Identification of Dihydroxycarbene13," Angew. Chem. Int. Ed., 2008, 47, 7071-7074, DOI: 10.1002/anie.200802105

InChIs

1: InChI=1/CH2O2/c2-1-3/h2-3H

A nice summary of the tunneling behavior of hydroxymethylene¹ was just published by Bucher in Angewandte Chemie.² Bucher strongly points out that the really novel part of this work is the very large barrier through which the proton tunnels. My blog post on this topic is here.

References

(1) Schreiner, P. R.; Reisenauer, H. P.; Pickard IV, F. C.; Simmonett, A. C.; Allen, W.
D.; Matyus, E.; Csaszar, A. G., "Capture of hydroxymethylene and its fast disappearance through tunnelling," Nature, 2008, 453, 906-909, DOI: 10.1038/nature07010.

(2) Bucher, G.; “Hydroxycarbene: Watching a Molecular Mole at Work,” Angew. Chem. Int. Ed., 2008, 47, 6957 – 6958, DOI: 10.1002/anie.200803195

The very simple carbene hydroxymethylene, HOCH, has finally been prepared and characterized.¹ Glyoxylic acid CHOCO₂H is subjected to high-vacuum laser photolysis. It fragments into HOCH, which is then trapped into an argon matrix. The experimental IR frequencies match up very well with the CCSD(T)/cc-pVQZ harmonic frequencies of the trans isomer 1t that are also adjusted for anharmonic effects. The computed vertical excitation energy of 415 nm matches well with the experimental value of the maximum absorption in the UV/vis spectra of 427 nm.

The other very interesting experimental result is that HOCH has a lifetime of about 2 hours in the matrix, while the deuterated species DOCH is stable. To explain these results, Schreiner, Allen and co-workers optimized a number of structures on the PES at CCSD(T)/cc-pVQZ and computed their energies using the focal point technique. The optimized structures and their relative energies are given in Figure 1.

1t (0.0)	TS2 (29.7)	2 (-52.1)
TS1(26.8)
1c (4.4)

Figure 1. Optimized CCSD(T)/cc-pVQZ structures of HOCH isomers and their Focal Point relative energies (kcal mol^-1).¹

The barriers for rearrangement from 1t are both very high. Rearrangement to formaldehyde 2 requires crossing a barrier of 29.7 kcal mol^-1, while the barrier to convert to the cis isomer 1c is 26.8 kcal mol^-1. (Note that from 1c a cleavage into CO and H2 can occur, but this barrier is another 47.0 kcal mol^-1.) These barriers are too large to be crossed at the very low temperatures of the matrices. However, using the intrinsic reaction potential at CCSD(T)/cc-pVQZ and WKB theory, the tunneling lifetime of HOCH is computed to be 122 minutes, in excellent accord with the experiment. The lifetime for DOCH is computed to be over 1200 years. Thus, the degradation of hydroxymethylene is entirely due to tunneling through a very large classical barrier! This rapid tunneling casts serious doubt on the ability to ever identify any hydroxymethylene in interstellar space.

References

(1) Schreiner, P. R.; Reisenauer, H. P.; Pickard IV, F. C.; Simmonett, A. C.; Allen, W.
D.; Matyus, E.; Csaszar, A. G., "Capture of hydroxymethylene and its fast disappearance through tunnelling," Nature, 2008, 453, 906-909, DOI: 10.1038/nature07010.

InChI

1: InChI=1/CH2O/c1-2/h1-2H
2: InChI=1/CH2O/c1-2/h1H2

Wes Borden has been exploring reactions where tunneling is operational. These studies have been inspired by Bill Doering’s¹ statement regarding tunneling in 1,5-sigmatropic shifts: “The tunneling effect is likely, in the opinion of some, to remain relegated to the virtual world of calculations”. Borden’s first two papers dealt with the kinetic isotope effects for the [1,5]-H shift in 1,3-cyclopentadiene and 5-methyl-1,3-cyclopentadiene.^2,3

His latest article examines carbon tunneling,⁴ which, due to the much heavier mass of the carbon nucleus relative to a proton, is likely to play a minimal role at best. Borden looked at the ring opening of cyclopropylcarbinyl radical 1 to 3-butene-1-yl radical 2, passing through transition state TS1-2. The B3LYP/6-31G(d) optimized structures are shown in Figure 1.

1		2
1	TS1-2	2

Figure 1. B3LYP/6-31G(d) optimized geometries of 1, 2, and TS1-2.⁴

The predicted rate of the reaction at 298 K using canonical variational transition state theory is increased by about 50% when small-curvature tunneling is included. This predicted rate is a bit smaller than the experimental value. Experiments also shows a linear Arrhenius plot, and Borden’s calculations agree until one reaches very low temperatures. Below 150 K the Arrhenius curve begins to deviate from linearity, and below 20 K the curve is flat – the rate is no longer temperature dependent! Thus, at cryogenic temperatures, the tunneling rate far exceeds traditional crossing of the variational barrier. Borden hopes that experimentalists will reinvestigate this problem (and hopefully confirm his predictions).

References

(1) Doering, W. v. E.; Zhao, X., "Effect on Kinetics by Deuterium in the 1,5-Hydrogen Shift of a Cisoid-Locked 1,3(Z)-Pentadiene, 2-Methyl-10-methylenebicycloJ. Am. Chem. Soc., 2006, 128, 9080-9085, DOI: 10.1021/ja057377v.

(2) Shelton, G. R.; Hrovat, D. A.; Borden, W. T., "Tunneling in the 1,5-Hydrogen Shift Reactions of 1,3-Cyclopentadiene and 5-Methyl-1,3-Cyclopentadiene," J. Am. Chem. Soc., 2007, 129, 164-168, DOI: 10.1021/ja0664279.

(3) Shelton, G. R.; Hrovat, D. A.; Borden, W. T., "Calculations of the Effect of Tunneling on the Swain-Schaad Exponents (SSEs) for the 1,5-Hydrogen Shift in 5-Methyl-1,3-cyclopentadiene. Can SSEs Be Used to Diagnose the Occurrence of Tunneling?," J. Am. Chem. Soc., 2007, 129, 16115-16118, DOI: 10.1021/ja076132a.

(4) Datta, A.; Hrovat, D. A.; Borden, W. T., "Calculations Predict Rapid Tunneling by Carbon from the Vibrational Ground State in the Ring Opening of Cyclopropylcarbinyl Radical at Cryogenic Temperatures," J. Am. Chem. Soc., 2008, 130, 6684-6685, DOI: 10.1021/ja801089p.

InChIs

1: InChI=1/C4H7/c1-4-2-3-4/h4H,1-3H2

2: InChI=1/C4H7/c1-3-4-2/h3H,1-2,4H2