Archive for the 'Molecules' Category

Triplet cyclobutadiene

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.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.



Figure 1. Optimized geometries of singlet and triplet 1.


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.


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

Aromaticity &cyclobutadiene Steven Bachrach 11 Sep 2017 1 Comment

Structure of GlyGly

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.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.




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.


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.


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

amino acids Steven Bachrach 17 Jul 2017 1 Comment


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.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,2 is computed to be about 119 kcal mol-1.


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

Rxn 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:


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.


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

Aromaticity &nanohoops Steven Bachrach 22 May 2017 No Comments

Conformationally selective tunneling

The Schreiner group has again reported an amazing experimental and computational study demonstrating a fascinating quantum mechanical tunneling effect, this time for the trifluoromethylhydroxycarbene (CF3COH) 2.1 (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.






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 106 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 CF3 group stabilizes the cis form probably through some weak HF 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.


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.


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

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

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

carbenes &Schreiner &Tunneling Steven Bachrach 07 Feb 2017 2 Comments

More examples of structure determination with computed NMR chemical shifts

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.1 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(CHCl3)–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.




Figure 1. Optimized geometries of 1-3.

Feng, Davis and coworkers have examined a series of anthroquionones from Australian marine sponges.2 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.


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.


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

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

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

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

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

NMR &terpenes Steven Bachrach 25 Oct 2016 No Comments

Dehydro-Diels-Alder Reactions

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.1 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 T1 diagnostic is fairly large for both the concerted and stepwise transition states, so they also performed CCSD(T)/CBS computations, which had much lower T1 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.2 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.3 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!


(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.


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

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

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

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

5: InChI=1S/C4H2/c1-3-4-2/h1-2H

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

benzynes &Cramer &Diels-Alder &Houk Steven Bachrach 25 Jul 2016 No Comments

Diels-Alder reaction of buckybowls

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.1

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.



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 C60. 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.


(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: <10.1002/chem.201502248.


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


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

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

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

Diels-Alder &fullerene Steven Bachrach 23 May 2016 No Comments

Calculating large fullerenes

What is the size of a molecule that will stretch computational resources today? Chan and co-workers have examined some very large fullerenes1 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 C60 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).

C60 + 10 benzene → 6 corannulene

Reaction 1

C60 + 10 naphthalene → 8 corannulene

Reaction 2

C60 + 10 phenanthrene → 10 corannulene

Reaction 3

C60 + 10 triphenylene → 12 corannulene

Reaction 4

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 C180, C240, and C320, 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.

C70 + 5 styrene → C60 + 5 naphthalene

Reaction 5

C180 → 3 C60

Reaction 6

C320 + 2/3 C60 → 2 C180

Reaction 7




Figure 1. B3LYP/6-31G(2df,p) optimized geometries of C180, C240, and C320. (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 C9680, still have heats of formation per carbon 1 kJ/mol greater than for graphene itself.


(1) Chan, B.; Kawashima, Y.; Katouda, M.; Nakajima, T.; Hirao, K. "From C60 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.

fullerene Steven Bachrach 22 Feb 2016 4 Comments

Highly efficient Buckycatchers

Capturing buckyballs involves molecular design based on non-covalent interactions. This poses interesting challenges for both the designer and the computational chemist. The curved surface of the buckyball demands a sequestering agent with a complementary curved surface, likely an aromatic curved surface to facilitate π-π stacking interactions. For the computational chemist, weak interactions, like dispersion and π-π stacking demand special attention, particularly density functionals designed to account for these interactions.

Two very intriguing new buckycatchers were recently prepared in the Sygula lab, and also examined by DFT.1 Compounds 1 and 2 make use of the scaffold developed by Klärner.2 In these two buckycatchers, the tongs are corranulenes, providing a curved aromatic surface to match the C60 and C70 surface. They differ in the length of the connector unit.

B97-D/TZVP computations of the complex of 1 and 2 with C60 were carried out. The optimized structures are shown in Figure 1. The binding energies (computed at B97-D/QZVP*//B97-D/TZVP) of these two complexes are really quite large. The binding energy for 1:C60 is 33.6 kcal mol-1, comparable to some previous Buckycatchers, but the binding energy of 2:C60 is 50.0 kcal mol-1, larger than any predicted before.



Figure 1. B97-D/TZVP optimized geometries of 1:C60and 2:C60.

Measurement of the binding energy using NMR was complicated by a competition for one or two molecules of 2 binding to buckyballs. Nonetheless, the experimental data show 2 binds to C60 and C70 more effectively than any previous host. They were also able to obtain a crystal structure of 2:C60.


(1) Abeyratne Kuragama, P. L.; Fronczek, F. R.; Sygula, A. "Bis-corannulene Receptors for Fullerenes Based on Klärner’s Tethers: Reaching the Affinity Limits," Org. Lett. 2015, ASAP, DOI: 10.1021/acs.orglett.5b02666.

(2) Klärner, F.-G.; Schrader, T. "Aromatic Interactions by Molecular Tweezers and Clips in Chemical and Biological Systems," Acc. Chem. Res. 2013, 46, 967-978, DOI: 10.1021/ar300061c.


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

2: InChI=1S/C66H36O2/c1-67-65-51-24-45-43-23-44(42-20-38-34-16-12-30-8-4-27-3-7-29-11-15-33(37(38)19-41(42)43)59-55(29)53(27)56(30)60(34)59)46(45)25-52(51)66(68-2)64-50-26-49(63(64)65)47-21-39-35-17-13-31-9-5-28-6-10-32-14-18-36(40(39)22-48(47)50)62-58(32)54(28)57(31)61(35)62/h3-22,24-25,43-44,49-50H,23,26H2,1-2H3/t43-,44+,49+,50-

Aromaticity &fullerene &host-guest Steven Bachrach 30 Nov 2015 No Comments

Bistetracene is a biradical singlet

Feng, Müller and co-workers have prepared a bistetracene analogue 1.1 This molecule displays some interesting features. While a closed shell Kekule structure can be written, a biradical structure results in more closed Clar rings, suggesting that perhaps the molecule is a ground state singlet biradical. The loss of NMR signals with increasing temperature along with an EPR signal that increases with temperature both support the notion of a ground state singlet biradical with a triplet excited state. The EPR measurement suggest as singlet-triplet gap of 3.4 kcal mol-1.

The optimized B3LYP/6-31G(d,p) geometries of the biradical singlet and triplet states are shown in Figure 1. The singlet is lower in energy by 6.7 kcal mol-1. The largest spin densities are on the carbons that carry the lone electron within the diradical-type Kekule structures.

singlet 1

triplet 1

Figure 1. B3LYP/6-31G(d,p) optimized geometries of the biradical singlet and triplet states of 1.


(1) Liu, J.; Ravat, P.; Wagner, M.; Baumgarten, M.; Feng, X.; Müllen, K. "Tetrabenzo[a,f,j,o]perylene: A Polycyclic Aromatic Hydrocarbon With An Open-Shell Singlet Biradical Ground State," Angew. Chem. Int. Ed. 2015, 54, 12442-12446, DOI: 10.1002/anie.201502657.


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

Aromaticity &diradicals Steven Bachrach 16 Nov 2015 No Comments

Next Page »