Schreiner and Grimme have examined a few compounds (see these previous posts) with long C-C bonds that are found in congested systems where dispersion greatly aids in stabilizing the stretched bond. Their new paper¹ continues this theme by examining 1 (again) and 2, using computations, and x-ray crystallography and gas-phase rotational spectroscopy and electron diffraction to establish the long C-C bond.

The distance of the long central bond in 1 is 1.647 Å (x-ray) and 1.630 Å (electron diffraction). Similarly, this distance in 2 is 1.642 Å (x-ray) and 1.632 Å (ED). These experiments discount any role for crystal packing forces in leading to the long bond.

A very nice result from the computations is that most functionals that include some dispersion correction predict the C-C distance in the optimized structures with an error of no more than 0.01 Å. (PW6B95-D3/DEF2-QZVP structures are shown in Figure 1.) Not surprisingly, HF and B3LYP without a dispersion correction predict a bond that is too long.) MP2 predicts a distance that is too short, but SCS-MP2 does a very good job.

Figure 1. PW6B95-D3/DEF2-QZVP optimized structures of 1 and 2.

References

1) Fokin, A. A.; Zhuk, T. S.; Blomeyer, S.; Pérez, C.; Chernish, L. V.; Pashenko, A. E.; Antony, J.; Vishnevskiy, Y. V.; Berger, R. J. F.; Grimme, S.; Logemann, C.; Schnell, M.; Mitzel, N. W.; Schreiner, P. R., "Intramolecular London Dispersion Interaction Effects on Gas-Phase and Solid-State Structures of Diamondoid Dimers." J. Am. Chem. Soc. 2017, 139, 16696-16707, DOI: 10.1021/jacs.7b07884.

InChIs

1: InChI=1S/C28H38/c1-13-7-23-19-3-15-4-20(17(1)19)24(8-13)27(23,11-15)28-12-16-5-21-18-2-14(9-25(21)28)10-26(28)22(18)6-16/h13-26H,1-12H2
InChIKey=MMYAZLNWLGPULP-UHFFFAOYSA-N

2: InChI=1S/C26H34O2/c1-11-3-19-15-7-13-9-25(19,21(5-11)23(27-13)17(1)15)26-10-14-8-16-18-2-12(4-20(16)26)6-22(26)24(18)28-14/h11-24H,1-10H2
InChIKey=VPBJYHMTINJMAE-UHFFFAOYSA-N

The Grimme group has examined all 1812 C₆₀ isomers, in part to benchmark some computational methods.¹ They computed all of these structures at PW6B95-D3/def2-QZVP//PBE-D3/def2-TZVP. The lowest energy structure is the expected fullerene 1 and the highest energy structure is the nanorod 2 (see Figure 1).

Figure 1. Optimized structures of the lowest (1) and highest (2) energy C₆₀ isomers.

About 70% of the isomers like in the range of 150-250 kcal mol^-1 above the fullerene 1, and the highest energy isomer 2 lies 549.1 kcal mol^-1 above 1. To benchmark some computational methods, they selected the five lowest energy isomers and five other isomers with higher energy to serve as a new database (C60ISO), with energies computed at DLPNO-CCSD(T)/CBS*. The mean absolute deviation of the PBE-D3/def2-TZVP relative energies with the DLPNO-CCSD(T)/CBS* energies is relative large 10.7 kcal mol^-1. However, the PW6B95-D3/def2-QZVP//PBE-D3/def2-TZVP method is considerably better, with a MAD of only 1.7 kcal mol^-1. This is clearly a reasonable compromise method for fullerene-like systems, balancing accuracy with computational time.

They also compared the relative energies of all 1812 isomers computed at PW6B95-D3/def2-QZVP//PBE-D3/def2-TZVP with a number of semi-empirical methods. The best results are with the DFTB-D3 method, with an MAD of 5.3 kcal mol^-1.

References

1) Sure, R.; Hansen, A.; Schwerdtfeger, P.; Grimme, S., "Comprehensive theoretical study of all 1812 C₆₀ isomers." Phys. Chem. Chem. Phys. 2017, 19, 14296-14305, DOI: 10.1039/C7CP00735C.

InChIs

1: InChI=1S/C60/c1-2-5-6-3(1)8-12-10-4(1)9-11-7(2)17-21-13(5)23-24-14(6)22-18(8)28-20(12)30-26-16(10)15(9)25-29-19(11)27(17)37-41-31(21)33(23)43-44-34(24)32(22)42-38(28)48-40(30)46-36(26)35(25)45-39(29)47(37)55-49(41)51(43)57-52(44)50(42)56(48)59-54(46)53(45)58(55)60(57)59
InChIKey=XMWRBQBLMFGWIX-UHFFFAOYSA-N

2: InChI=1S/C60/c1-11-12-2-21(1)31-41-32-22(1)3-13(11)15-5-24(3)34-43(32)53-55-47-36-26-6-16-17-7(26)28-9-19(17)20-10-29-8(18(16)20)27(6)37-46(36)54(51(41)55)52-42(31)33-23(2)4(14(12)15)25(5)35-44(33)58-56(52)48(37)39(29)50-40(30(9)10)49(38(28)47)57(53)59(45(34)35)60(50)58
InChIKey=AGZHNPDQKMDYHI-UHFFFAOYSA-N

Cyclobutadiene has long fascinated organic chemists. It is the 4e analogue of the 6e benzene molecule, yet it could hardly be more different. Despite nearly a century of effort, cyclobutadiene analogues were only first prepared in the 1970s, reflecting its strong antiaromatic character.

Per-trimethylsilylcyclobutadiene 1 offers opportunities to probe the properties of the cyclobutadiene ring as the bulky substituents diminish dimerization and polymerization of the reactive π-bonds. Kostenko and coworkers have now reported on the triplet state of 1.¹ They observe three EPR signals of 1 at temperatures above 350 K, and these signals increase in area with increasing temperature. This is strong evidence for the existence of triplet 1 in equilibrium with the lower energy singlet. Using the variable temperature EPR spectra, the singlet triplet gap is 13.9 ± 0.8 kcal mol^-1.

The structures of singlet and triplet 1 were optimized at B3LYP-D3/6-311+G(d,p) and shown in Figure 1. The singlet is the expected rectangle, with distinctly different C-C distance around the ring. The triplet is a square, with equivalent C-C distances. Since both the singlet and triplet states are likely to have multireference character, the energies of both states were obtained at RI-MRDDCI2-CASSCF(4,4)/def2-SVP//B3LYPD3/6-311+G(d,p) and give a singlet-triplet gap of 11.8 kcal mol^-1, in quite reasonable agreement with experiment.

singlet

triplet

Figure 1. Optimized geometries of singlet and triplet 1.

References

1. Kostenko, A.; Tumanskii, B.; Kobayashi, Y.; Nakamoto, M.; Sekiguchi, A.; Apeloig, Y., "Spectroscopic Observation of the Triplet Diradical State of a Cyclobutadiene." Angew. Chem. Int. Ed. 2017, 56, 10183-10187, DOI: 10.1002/anie.201705228.

InChIs

1: InChI=1S/C16H36Si4/c1-17(2,3)13-14(18(4,5)6)16(20(10,11)12)15(13)19(7,8)9/h1-12H3
InChIkey=AYOHYRSQVCLGKR-UHFFFAOYSA-N

Continuing their application of laser ablation molecular beam Fourier transform microwave (LA-MB-FTMW) spectroscopy and computational chemistry to biochemical molecules (see these previous posts), the Alonso group reports on the structure of the glycine-glycine dipeptide 1.¹ The microwave spectrum shows three different conformers. MP2/6-311++G(d,p) computations, the same method they have previously utilized for predicting geometries, revealed a number of different conformations. By matching the spectroscopic parameters obtained from the spectrum with those of the computed structures, they proposed the three conformations 1a, 1b, and 1c, shown in Figure 1.

1a

1b

1c

Figure 1. ωb97xd/6-31G(d) optimized structures of the three conformers of 1.
Note that the authors did not report their structures in their supporting materials(!) so I have optimized them.

The structures of conformers 1a and 1b are nearly planar. MP2 predicts a non-planar rotomer of 1a, which brings the carboxyl group out of plane, to be the lowest conformation in terms of electronic energy. With the M06-2x functional, this non-planar rotomer is about isoenergetic with 1a. With all computational levels 1a is the lowest in free energy. The barrier for rotation between the non-planar rotomer and 1a is very small, and this explains why it is not observed in the supersonic expansion.

References

1) Cabezas, C.; Varela, M.; Alonso, J. L., "The Structure of the Elusive Simplest Dipeptide Gly-Gly." Angew. Chem. Int. Ed. 2017, 56, 6420-6425, DOI: 10.1002/anie.201702425.

InChIs

1: InChI=1S/C4H8N2O3/c5-1-3(7)6-2-4(8)9/h1-2,5H2,(H,6,7)(H,8,9)
InChIKey=YMAWOPBAYDPSLA-UHFFFAOYSA-N

The synthesis of components of nanostructures (like fullerenes and nanotubes) has dramatically matured over the past few years. I have blogged about nanohoops before, and this post presents the recent work of the Itami group in preparing the nanobelt 1.¹

The synthesis is accomplished through a series of Wittig reactions with an aryl-aryl coupling to stitch together the final rings. The molecule is characterized by NMR and x-ray crystallography. The authors have also computed the structure of 1 at B3LYP/6-31G(d), shown in Figure 1. The computed C-C distances match up very well with the experimental distances. The strain energy of 1, presumably estimated by Reaction 1,² is computed to be about 119 kcal mol^-1.

1

Figure 1. B3LYP/6-31G(d) optimized structure of 1.

NICS(0) values were obtained at B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d); the rings along the middle of the belt have values of -7.44ppm and are indicative of normal aromatic 6-member rings, while the other rings have values of -2.00ppm. This suggests the dominant resonance structure shown below:

References

1) Povie, G.; Segawa, Y.; Nishihara, T.; Miyauchi, Y.; Itami, K., "Synthesis of a carbon nanobelt." Science 2017, 356, 172-175, DOI: 10.1126/science.aam8158.

2) Segawa, Y.; Yagi, A.; Ito, H.; Itami, K., "A Theoretical Study on the Strain Energy of Carbon Nanobelts." Org. Letters 2016, 18, 1430-1433, DOI: 10.1021/acs.orglett.6b00365.

InChIs:

1: InChI=1S/C48H24/c1-2-26-14-40-28-5-6-31-20-44-32(19-42(31)40)9-10-34-24-48-36(23-46(34)44)12-11-35-21-45-33(22-47(35)48)8-7-30-17-41-29(18-43(30)45)4-3-27-15-37(39(26)16-28)25(1)13-38(27)41/h1-24H
InChIKey=KJWRWEMHJRCQKK-UHFFFAOYSA-N

The Schreiner group has again reported an amazing experimental and computational study demonstrating a fascinating quantum mechanical tunneling effect, this time for the trifluoromethylhydroxycarbene (CF₃COH) 2.¹ (I have made on a number of posts discussing a series of important studies in this field by Schreiner.) Carbene 2 is formed, in analogy to many other hydroxycarbenes, by flash vapor pyrolysis of the appropriate oxoacid 1 and capturing the products on a noble gas matrix.

Carbene 2t is observed by IR spectroscopy, and its structure is identified by comparison with the computed CCSD(T)/cc-pVTZ frequencies. When 2t is subjected to 465 nm light, the signals for 2t disappear within 30s, and two new species are observed. The first species is the cis conformer 2c, confirmed by comparison with its computed CCSD(T)/cc-pVTZ frequencies. This cis conformer remains even with continued photolysis. The other product is determined to be trifluoroacetaldehyde 3. Perhaps most interesting is that 2t will convert to 3 in the absence of light at temperatures between 3 and 30 K, with a half-life of about 144 h. There is little rate difference at these temperatures. These results are quite indicative of quantum mechanical tunneling.

To aid in confirming tunneling, they computed the potential energy surface at CCSD(T)/cc-pVTZ. The trans isomer is 0.8 kcal mol^-1 lower in energy that the cis isomer, and this is much smaller than for other hydroxycarbenes they have examined. The rotational barrier TS1 between the two isomer is quite large, 26.4 kcal mol^-1, precluding their interchange by classical means at matrix temperatures. The barrier for conversion of 2t to 3 (TS2) is also quite large, 30.7 kcal mol^-1, and insurmountable at 10K by classical means. No transition state connecting 2c to 3 could be located. These geometries and energies are shown in Figure 1.

2c 0.8

TS1 26.4

2t 0.0

TS2 30.7

3 -49.7

Figure 1. Optimized geometries at CCSD(T)/cc-pVTZ. Relative energies (kcal mol^-1) of each species are listed as well.

WKB computations at M06-2X/6-311++G(d,p) predict a half-life of 172 h, in nice agreement with experiment. The computed half-life for deuterated 2t is 10⁶ years, and the experiment on the deuterated analogue revealed no formation of deuterated 3.

The novel component of this study is that tunneling is conformationally selective. The CF₃ group stabilizes the cis form probably through some weak H^…F interaction, so that the cis isomer can be observed, but no tunneling is observed from this isomer. Only the trans isomer has the migrating hydrogen atom properly arranged for a short hop over to the carbon, allowing the tunneling process to take place.

References

1) Mardyukov, A.; Quanz, H.; Schreiner, P. R., "Conformer-specific hydrogen atom tunnelling in trifluoromethylhydroxycarbene." Nat. Chem. 2017, 9, 71–76, DOI: 10.1038/nchem.2609.

InChIs

1: =1S/C3HF3O3/c4-3(5,6)1(7)2(8)9/h(H,8,9)
InChIKey=GVDJEHMDNREMFA-UHFFFAOYSA-N

2: InChI=1S/C2HF3O/c3-2(4,5)1-6/h6H
InChIKey=FVJVNIREIXAWKU-UHFFFAOYSA-N

3: InChI=1S/C2HF3O/c3-2(4,5)1-6/h1H
InChIKey=JVTSHOJDBRTPHD-UHFFFAOYSA-N

Use of computed NMR chemical shifts in structure determination is really growing fast. Presented here are a couple of recent examples.

Nguyen and Tantillo used computed chemical shifts with the DP4 analysis to identify the structure of three terpenes 1-3.¹ They optimized the geometries of all of the diastereomers of each compound, along with multiple conformations of each diastereomer, at B3LYP/6-31+G(d,p) and then computed the chemical shifts at SMD(CHCl₃)–mPW1PW91/6-311+G(2d,p). The chemical shifts were Boltzmann weighted including all conformations within 3 kcal mol^-1 of the lowest energy structure.

For 1, the DP4 analysis using just the proton shifts predicted a different isomer than using the carbon shifts, but when combined, DP4 predicted the structure, with 98.8% confidence, shown in the scheme above, and in Figure 1. For 2, the combined proton and carbon shift analysis with DP4 indicated a 100% confidence of the structure shown in the scheme and Figure 1. Lastly, for 3, which is more complicated due to the conformations of the 9-member ring, DP4 predicts with 100% confidence the structure shown in the scheme and Figure 1.

1

2

3

Figure 1. Optimized geometries of 1-3.

Feng, Davis and coworkers have examined a series of anthroquionones from Australian marine sponges.² The structure of one compound was a choice of two options: 4 or 5. Initial geometries were obtain by molecular mechanics and the low energy isomers were then reoptimized at B3LYP/6-31+G(d,p). The chemical shifts were computed using PCM/MPW1PW91/6-311+G(2d,p). Application of the DP4 method indicate the structure to be 4 with a 100% confidence level. The lowest energy conformer of 4 is shown in Figure 2.

Figure 2. Optimized geometry of 4.

References

1) Nguyen, Q. N. N.; Tantillo, D. J. “Using quantum chemical computations of NMR chemical shifts to assign relative configurations of terpenes from an engineered Streptomyces host,” J. Antibiotics 2016, 69, 534–540, DOI: 10.1038/ja.2016.51.

2) Khokhar, S.; Pierens, G. K.; Hooper, J. N. A.; Ekins, M. G.; Feng, Y.; Rohan A. Davis, R. A. “Rhodocomatulin-Type Anthraquinones from the Australian Marine Invertebrates Clathria hirsuta and Comatula rotalaria,” J. Nat. Prod., 2016, 79, 946–953, DOI: 10.1021/acs.jnatprod.5b01029.

InChIs

1: InChI=1S/C15H24/c1-10-5-6-15(4)8-11-7-14(2,3)9-12(11)13(10)15/h9-11,13H,5-8H2,1-4H3/t10-,11+,13-,15+/m1/s1
InChIKey=KVSCZIPUFBVHBM-OICBVUGWSA-N

2: InChI=1S/C15H24/c1-10-5-6-15(4)8-11-7-14(2,3)9-12(11)13(10)15/h5,11-13H,6-9H2,1-4H3/t11-,12-,13+,15-/m0/s1
InChIKey=ZLYGJLHCPYVGDA-XPCVCDNBSA-N

3: InChI=1S/C20H32/c1-14-6-9-18-19(3,4)10-11-20(18,5)13-17-15(2)7-8-16(17)12-14/h6,13,15-16,18H,7-12H2,1-5H3/b14-6-,17-13-/t15-,16-,18-,20+/m0/s1
InChIKey=JZGOFJIAHJJJDK-ICZJPRMTSA-N

4: InChI=1S/C18H14O7/c1-7(19)13-10(20)6-11(21)15-16(13)17(22)9-4-8(24-2)5-12(25-3)14(9)18(15)23/h4-6,20-21H,1-3H3
InChIKey=MPQMZEXRJVMYBT-UHFFFAOYSA-N

5: InChI=1S/C18H14O7/c1-7(19)13-10(20)6-11(21)15-16(13)14-9(17(22)18(15)23)4-8(24-2)5-12(14)25-3/h4-6,20-21H,1-3H3
InChIKey=WIKIUXNPFURKNF-UHFFFAOYSA-N

I have been delinquent in writing about the dehydro-Diels-Alder reactions, but really can’t put it off any further. These sets of reactions really deserve a fuller analysis than I am going to summarize here, but this post will provide a good jumping off point for anyone interested in further investigation.

So the Diels-Alder reaction is among the most famous and most important reactions in organic chemistry. The reaction creates a 6-member ring and sets up to four stereocenters. In the past couple of years many chemists have expressed interest in the variant where the four-carbon component is more highly unsaturated, i.e. enyne or diyne. I will summarize the results of three recent computational papers dealing with the reaction of a diyne with an yne.

The first paper is by Skraba-Joiner, Johnson, and Agarwal.¹ They discuss, among a number of interesting pericyclic reactions, the intramolecular Diels-Alder reaction of triyne 1 to give 2. They examined a concerted and stepwise pathway at (U)M05-2X/6-311+G(d,p) and find the concerted to be favored by 6.0 kcal mol^-1. CCSD(T) using these geometries increases the difference to 8.2 kcal mol^-1. The T₁ diagnostic is fairly large for both the concerted and stepwise transition states, so they also performed CCSD(T)/CBS computations, which had much lower T₁ values. The concerted TS remained favorable, but by only 2.7 kcal mol^-1.

In the same special issue of the Journal of Organic Chemistry, Cramer, Hoye, and Kuwata examined a reaction closely related to what Johnson examined above.² They looked at the reaction taking 3 into 4 via both experiments and computations. The M06-2x/6-311+G(d,p) geometries for the concerted and first TS along the stepwise path (with R1=R2=H) are shown in Figure 1. Evaluating the energies at SMD(o-dichlorobenzene)/B3LYP-D3BJ/6-311+G-(d,p)//M06-2X/6-311+G(d,p) find in this case (along with all of the other R1/R2 variants they examined) that the stepwise path has a lower barrier than the concerted path. In the case where R1=R2=H, the stepwise path is favored by 6.0 kcal mol^-1. Additionally, these stepwise barriers are in reasonable agreement with the experimentally-derived barriers.

Concerted TS

Stepwise TS

Figure 1. M06-2x/6-311+G(d,p) optimized geometries of the concerted and stepwise TSs for the reaction of 3H going to 4H.

It should be pointed out that the wavefunctions for the concerted TSs were all found to be unstable with regard to a restricted to unrestricted relaxation. Given this problem, they also performed a CASPT2 energy evaluation of the concerted and stepwise transition states for the case R1=R2=H. CASPT2 finds the stepwise barrier to be 3.7 kcal mol^-1 lower than the concerted barrier.

The last paper comes from the Houk lab, and examines the simplest set of intermolecular dehdro-Diels-Alder reactions.³ I will focus here on the most unsaturated analogue, the reaction of 1,3-butadiyne 5 with ethyne to give benzyne 6.

The concreted and stepwise transition states for this reaction (at (U)M06-2X/6-311+G(d,p)) are shown in Figure 2. The concerted barrier is 36.0 kcal moml^-1 while the stepwise barrier is slightly lower: 35.2 kcal mol^-1. The distortion energy for the concerted reaction is large (43.2 kcal mol^-1) due mostly to angle changes in the diyne. Its interaction energy is -7.2 kcal mol^-1, similar to the interaction energy in other similar Diels-Alder reactions. In contrast, the distortion energy for the stepwise pathway is 27.5 kcal mol^-1, but the interaction energy is +7.7 kcal mol^-1. These values are very similar to the distortion and interaction energy of the related (but less saturated DA reactions).

Concerted TS

Stepwise TS

Figure 2. (U)M06-2X/6-311+G(d,p) optimized concerted and stepwise TS for the reaction of 1,3-diyne with ethyne.

Molecular dynamics trajectories for both the concerted and stepwise paths reveal interesting differences. The concerted trajectories show an oscillatory behaviour of bending the angles at the C2 and C3 carbons prior to the TS, and then near synchronous formation of the new C-C bonds. The trajectories initiated at the stepwise TS show no systematic motion. Once the bond is formed, the biradical exhibits a long lifetime, on the order of picoseconds, much longer than the trajectory runs.

These three studies indicate the nature of the dehydro Diels-Alder reaction is very sensitive to reaction conditions, substituents, solvation, and all other manner of effects and will likely prove an area of interest for some time. It should keep a number of computational chemists busy for some time!

References

(1) Skraba-Joiner, S. L.; Johnson, R. P.; Agarwal, J. "Dehydropericyclic Reactions: Symmetry-Controlled Routes to Strained Reactive Intermediates," J. Org. Chem. 2015, 80, 11779-11787, DOI: 10.1021/acs.joc.5b01488.

(2) Marell, D. J.; Furan, L. R.; Woods, B. P.; Lei, X.; Bendelsmith, A. J.; Cramer, C. J.; Hoye, T. R.; Kuwata, K. T. "Mechanism of the Intramolecular Hexadehydro-Diels–Alder Reaction," J. Org. Chem. 2015, 80, 11744-11754, DOI: 10.1021/acs.joc.5b01356.

(3) Yu, P.; Yang, Z.; Liang, Y.; Hong, X.; Li, Y.; Houk, K. N. "Distortion-Controlled Reactivity and Molecular Dynamics of Dehydro-Diels–Alder Reactions," J. Am. Chem. Soc. 2016, 138, 8247-8252, DOI: 10.1021/jacs.6b04113.

InChIs

1: InChI=1S/C9H8/c1-3-5-7-9-8-6-4-2/h1-2H,5,7,9H2
InChIKey=IYZAZSVBWMMSLQ-UHFFFAOYSA-N

2: InChI=1S/C9H8/c1-2-5-9-7-3-6-8(9)4-1/h1,4H,3,6-7H2
InChIKey=PZJMTUKDGZUDBH-UHFFFAOYSA-N

3H: InChI=1S/C8H4O2/c1-3-5-6-7-10-8(9)4-2/h1-2H,7H2
InChIKey=MGXDIFXPYGGQLF-UHFFFAOYSA-N

4H: InChI=1S/C10H8O4/c1-6(11)14-8-2-3-9-7(4-8)5-13-10(9)12/h2-4H,5H2,1H3
InChIKey=GEFLHLNIKGXWCA-UHFFFAOYSA-N

5: InChI=1S/C4H2/c1-3-4-2/h1-2H
InChIKey=LLCSWKVOHICRDD-UHFFFAOYSA-N

6: InChI=1S/C6H4/c1-2-4-6-5-3-1/h1-4H
InChIKey=KLYCPFXDDDMZNQ-UHFFFAOYSA-N

Fullerenes can undergo the Diels-Alder reaction with some specificity: the diene preferentially adds across the bond shared by two fused 6-member rings over the bond shared by the fused 6- and 5-member rings. Garcia-Rodeja and colleagues have examined the analogous Diels-Alder reaction of cyclopentadiene with five curved aromatic compounds, 1-5.¹

The computations were performed at BP86-D3/def2-TZVPP//RI-BP86-D3/def2-SVP. Representative transition states for the addition of cyclopentadiene with 3 over the 6,6-bond and 5,6-bond are shown in Figure 1.

5,6-bond

6,6-bond

Figure 1. RI-BP86-D3/def2-SVP optimized transition states for the reaction of cyclopentadiene with 3.

For the reactions of cyclopentadiene with 2-5 the reactions with the 6,6-bond is both kinetically and thermodynamically favored, while with 1 the 6,6-bond is kinetically preffered and the 5,6-adduct is the thermodynamic product. As the molecules increase in size (from 1 to 5), the activation barrier decreases, and the barrier for the reaction with 5 is only 1.4 kcal mol^-1larger than the barrier with C₆₀. The reaction energy also becomes more exothermic with increasing size. There is a very good linear relationship between activation barrier and reaction energy.

Use of the distortion/interaction model indicates that the preference for the 6,6-regioselectivity come from better interaction energy than for the 5,6-reaction, and this seems to come about by better orbital overlap between the cyclopentadiene HOMO and the 6,6-LUMO of the buckybowl.

References

(1) García-Rodeja, Y.; Solà, M.; Bickelhaupt , F. M.; Fernández, I. "Reactivity and Selectivity of Bowl-Shaped Polycyclic Aromatic Hydrocarbons: Relationship to C60," Chem. Eur. J. 2016, 22, 1368-1378, DOI: .

InChIs

1: InChI=1S/C20H10/c1-2-12-5-6-14-9-10-15-8-7-13-4-3-11(1)16-17(12)19(14)20(15)18(13)16/h1-10H
InChIKey=VXRUJZQPKRBJKH-UHFFFAOYSA-N

2: InChIKey=ASIFYFRJNYNQLA-UHFFFAOYSA-N

3: InChI=1S/C26H12/c1-5-13-14-6-2-11-19-20-12-4-8-16-15-7-3-10-18-17(9-1)21(13)25(22(14)19)26(23(15)18)24(16)20/h1-12H
InChIKey=OUWFOTSXASFGQD-UHFFFAOYSA-N

4: InChI=1S/C30H12/c1-2-14-6-10-18-20-12-8-16-4-3-15-7-11-19-17-9-5-13(1)21-22(14)26(18)29(25(17)21)30-27(19)23(15)24(16)28(20)30/h1-12H
InChIKey=JEUCRZPADDQRKU-UHFFFAOYSA-N

5: InChI=1S/C36H12/c1-7-16-17-9-3-14-5-11-20-21-12-6-15-4-10-19-18-8-2-13(1)22-25(16)31-32(26(18)22)34-28(19)24(15)30(21)36(34)35-29(20)23(14)27(17)33(31)35/h1-12H
InChIKey=QMGQDOOJOCPYIA-UHFFFAOYSA-N

What is the size of a molecule that will stretch computational resources today? Chan and co-workers have examined some very large fullerenes¹ to both answer that question, and also to explore how large a fullerene must be to approach graphene-like properties.

They are interested in predicting the heat of formation of large fullerenes. So, they benchmark the heats of formation of C₆₀ using four different isodesmic reactions (Reaction 1-4), comparing the energies obtained using a variety of different methods and basis sets to those obtained at W1h. The methods include traditional functionals like B3LYP, B3PW91, CAM-B3LYP, PBE1PBE, TPSSh, B98, ωB97X, M06-2X3, and MN12-SX, and supplement them with the D3 dispersion correction. Additionally a number of doubly hybrid methods are tested (again with and without dispersion corrections), such as B2-PLYP, B2GPPLYP, B2K-PLYP, PWP-B95, DSD-PBEPBE, and DSD-B-P86. The cc-pVTZ and cc-pVQZ basis sets were used. Geometries were optimized at B3LYP/6-31G(2df,p).

Excellent results were obtained with DSD-PBEPBE-D3/cc-pVQZ (an error of only 1.8 kJ/mol), though even a method like BMK-D3/cc-pVTZ had an error of only 9.2 kJ/mol. They next set out to examine large fullerenes, including such behemoths as C₁₈₀, C₂₄₀, and C₃₂₀, whose geometries are shown in Figure 1. Heats of formation were obtained using isodesmic reactions that compare back to smaller fullerenes, such as in Reaction 5-8.

C180

C240

C320

Figure 1. B3LYP/6-31G(2df,p) optimized geometries of C₁₈₀, C₂₄₀, and C₃₂₀. (Don’t forget that clicking on these images will launch Jmol and allow you to manipulate the molecules in real-time.)

Next, taking the heat of formation per C for these fullerenes, using a power law relationship, they were able to extrapolate out the heat of formation per C for truly huge fullerenes, and find the truly massive fullerenes, like C₉₆₈₀, still have heats of formation per carbon 1 kJ/mol greater than for graphene itself.

References

(1) Chan, B.; Kawashima, Y.; Katouda, M.; Nakajima, T.; Hirao, K. "From C₆₀ to Infinity: Large-Scale Quantum Chemistry Calculations of the Heats of Formation of Higher Fullerenes," J. Am. Chem. Soc. 2016, 138, 1420-1429, DOI: 10.1021/jacs.5b12518.