Schaefer, Csaszar, and Allen have applied the focal point method towards predicting the energies and structures of cysteine.¹ This very high level method refines the structures that can be used to compare against those observed by Alonso² in his laser ablation molecular beam Fourier transform microwave spectroscopy experiment (see this post). They performed a broad conformation search, initially examining some 66,664 structures. These reduced to 71 unique conformations at MP2/cc-pvTZ. The lowest 11 energy structures were further optimized at MP2(FC)/aug-cc-pV(T+d)Z. The four lowest energy conformations are shown in Figure 1 along with their relative energies.

I (0.0)	II (4.79)
III (5.81)	IV (5.95)

Figure 1. MP2(FC)/aug-cc-pV(T+d)Z optimized geometries and focal point relative energies (kJ mol^-1) of the four lowest energy conformers of cysteine.¹

The three lowest energy structures found here match up with the lowest two structures found by Alonso and the energy differences are also quite comparable: 4.79 kJ and 5.81 mol^-1 with the focal point method 3.89 and 5.38 kJ mol^-1 with MP4/6-311++G(d,p)// MP2/6-311++G(d,p). So the identification of the cysteine conformers made by Alonso remains on firm ground.

References

(1) Wilke, J. J.; Lind, M. C.; Schaefer, H. F.; Csaszar, A. G.; Allen, W. D., "Conformers of Gaseous Cysteine," J. Chem. Theory Comput. 2009, DOI: 10.1021/ct900005c.

(2) Sanz, M. E.; Blanco, S.; López, J. C.; Alonso, J. L., "Rotational Probes of Six Conformers of Neutral Cysteine," Angew. Chem. Int. Ed. 2008, 4, 6216-6220, DOI: 10.1002/anie.200801337

InChIs

Cysteine:
InChI=1/C3H7NO2S/c4-2(1-7)3(5)6/h2,7H,1,4H2,(H,5,6)/t2-/m0/s1
InChIKey: XUJNEKJLAYXESH-REOHCLBHBU

Malonaldehyde 1 possesses a very short intramolecular hydrogen bond. Its potential energy surface has two local minima (the two mirror image hydrogen-bonded structures) separated by a C_2v transition state. Schaefer reports a high-level computational study for the search for even shorter hydrogen bonds that might even lead to a single well on the PES.¹

The hydrogen bond distance is characterized by the non-bonding separation between the two oxygen atoms. Table 1 shows the O^…O distance for a number of substituted malonaldehydes computed at B3LYP/DZP++. Electron withdrawing groups on C₂ reduce the O^..O distance (see trend in 1 → 4). Electron donating groups on C₁ and C₃ also reduce this distance (see 5 and 6). Bulky substituents on the terminal carbons also reduce the O^…O distance (see 7). Combining all of these substituent effects in 8 leads to the very short O^…O distance of 2.380 Å.

Table 1. Distance (Å) between the two oxygen atoms and the barrier for hydrogen transfer of substituted malonaldehydes .¹

^aFocal point energy. ^bFocal point energy and corrected for zero-point vibrational energy.

A shorter O^…O distance might imply a smaller barrier for hydrogen transfer between the two oxygens. The structures of 8 and the transition state for its hydrogen transfer are shown in Figure 1. The energies of a number of substituted malonaldehydes were computed using the focal point method, and the barriers for hydrogen transfer are listed in Table 1. There is a nice correlation between the O^…O distance and the barrier height. The barrier for 8 is quite small, suggesting that with some bulkier substituents, the barrier might vanish altogether, leaving only a symmetric structure. In fact, the barrier appears to vanish when zero-point vibrational energies are included.

8

8TS

Figure 1. B3LYP/DZP++ optimized geometries of 8 and the transition state for hydrogen transfer 8TS.¹

References

(1) Hargis, J. C.; Evangelista, F. A.; Ingels, J. B.; Schaefer, H. F., "Short Intramolecular Hydrogen Bonds: Derivatives of Malonaldehyde with Symmetrical Substituents," J. Am. Chem. Soc., 2008, 130, 17471-17478, DOI: 10.1021/ja8060672.

InChIs

1: InChI=1/C3H4O2/c4-2-1-3-5/h1-4H/b2-1-
InChIKey=GMSHJLUJOABYOM-UPHRSURJBI

2: InChI=1/C4H3NO2/c5-1-4(2-6)3-7/h2-3,6H/b4-2-
InChIKey=BHYIQMFSOGUTRT-RQOWECAXBC

3: InChI=1/C3H3NO4/c5-1-3(2-6)4(7)8/h1-2,5H/b3-1+
InChIKey=JBBHDCMVSJADCE-HNQUOIGGBS

4: InChI=1/C3H5BO2/c4-3(1-5)2-6/h1-2,5H,4H2/b3-1+
InChIKey=IQNKNZSFMBIPBI-HNQUOIGGBX

5: InChI=1/C3H6N2O2/c4-2(6)1-3(5)7/h1,6H,4H2,(H2,5,7)/b2-1-/f/h5H2
InChIKey=AOZIOAJNRYLOAH-KRHGAQEYDI

6: InChI=1/C5H8O4/c1-8-4(6)3-5(7)9-2/h3,6H,1-2H3/b4-3+
InChIKey=BYYYYPBUMVENKB-ONEGZZNKBI

7: InChI=1/C11H20O2/c1-10(2,3)8(12)7-9(13)11(4,5)6/h7,12H,1-6H3/b8-7-
InChIKey=SOZFXLUMSLXZFW-FPLPWBNLBX

8: InChI=1/C3H5N3O4/c4-2(7)1(3(5)8)6(9)10/h7H,4H2,(H2,5,8)/b2-1+/f/h5H2
InChIKey=IHYUFGCOUITNJP-CHFMFTGODK

Schaefer and Allen have applied their focal point method to the question of the benzylic effect in the S_N2 reaction.¹ S_N2 reactions are accelerated when the attack occurs at the benzylic carbon, a well-known phenomenon yet the reason for this remains unclear. The standard textbook-like argument has been that the negative charge built up in the S_N2 transition state can be delocalized into the phenyl ring. However, solution phase Hammett studies are often U-shaped, indicating that both electron donating and withdrawing group accelerate the substitution reaction. (This is usually argued as indicative of a mechanism change from S_N2 to S_N1.)

The focal point method involves a series of very large computations where both basis set size and degree of electron correlation are systematically increased, allowing for an extrapolation to essentially infinite basis set and complete correlation energy. The energy of the transition state (relative to separated reactants) for four simple S_N2 reactions evaluated with the focal point method are listed in Table 1. The barrier for the benzylic substitutions is lower than for the methyl cases, indicative of the benzylic effect.

Table 1. Energy (kcal mol^-1) of the transition state relative to reactants.¹

To answer the question of why the benzylic substitution reactions are faster, they examined the charge distribution evaluated at B3LYP/DZP++. As seen in Table 1, this method does not accurately reproduce the activation barriers, but the errors are not terrible, and the trends are correct.

In Figure 1 are the geometries of the transition states for the reaction of fluoride with methylflouride or benzylfluoride. The NBO atomic charges show that the phenyl ring acquired very little negative charge at the transition state. Rather, the electric potential at the carbon under attack is much more revealing. The potential is significantly more positive for the benzylic carbon than the methyl carbon in both the reactant and transition states.

VC = -405.156 V	VC = -404.379 V

Figure 1. MP2/DZP++ transition states for the reaction of fluoride with methylfluoride and benzylflouride. NBO charges on F and C and the electrostatic potential in Volts.¹

They next examined the reaction of fluoride with a series of para-substituted benzylfluorides. The relation between the Hammet σ constants and the activation energy is fair (r = 0.971). But the relation between the electrostatic potential at the benzylic carbon (in either the reactant or transition state) with the activation energy is excellent (r = 0.994 or 0.998). Thus, they argue that it is the increased electrostatic potential at the benzylic carbon that accounts for the increased rate of the S_N2 reaction.

References

(1) Galabov, B.; Nikolova, V.; Wilke, J. J.; Schaefer III, H. F.; Allen, W. D., "Origin of the S_N2 Benzylic Effect," J. Am. Chem. Soc., 2008, 130, 9887-9896, DOI: 10.1021/ja802246y.

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