Chemists are constantly checking the limits of theories, and the limits of bonding is one that has been subject to many tests of late. I have posted on two recent papers (here, here) that probe just how long a C-C bond can be, and now Li, Miller, and co-workers report a structure that pushes that limit even further out.1

They prepared and obtained the x-ray structure of five
derivatives of o-carborane,
namely compounds 1, 2a, 3a, 3b and 4. In all of these, the C-C bond in the carborane is stretched well beyond that of a typical C-C
bond (see Table 1). The longest case is in 3b
where the C-C bond length is a whopping 1.931 Å
(see Figure 1), which obliterates the previous record holder at 1.798
Å.2 B3PW91-D3/cc-pVTZ computations corroborate these structures and the long
C-C bond.

Scheme 1: Carboranes with long C-C bonds (highlighted in blue)

Table 1. C-C bond distance (Å)

cmpd

r(C-C) expt

r(C-C) DFT

1

1.829

1.839

2a

1.720

1.710

3a

1.893

1.917

3b

1.931

1.936

4

1.627

1.607

Figure 1. B3PW91-D3/cc-pVTZ optimized structure of 3b.

Topological electron density analysis locates a bond path
between the two carbons in all five structures. The Wiberg
bond index is small, with a value of only 0.34 in 3b. Natural bond orbital (NBO) analysis identifies a
negative hyperconjugation interaction between the nitrogen lone pair
and the σ*C-C orbital. This
rationalizes both the very long C-C bond and the very short C-N bonds, and the
trends associated with the variation between 1°
amine, 2° amine and imine.

Pericyclic reactions remain a fruitful area of research despite the seminal publication of the Woodward-Hoffmann rules decades ago. Here are two related papers of pericyclic reactions that violate the Woodward-Hoffmann rules.

First, Solomek, Ravat, Mou, Kertesz, and Jurícek reported on the thermal and photochemical electrocyclization reaction of diphenylcetherene 1a.1 Though they were not able to directly detect the intermediate 2, through careful examination of the photochemical reaction, they were able to infer that the thermal cyclization goes via the formally forbidden conrotatory pathway (see Scheme 1).

Scheme 2.

Kinetic studies estimate the activation barrier is 14.1 kcal mol-1. They performed DFT computations of the parent 1b using a variety of functionals with both restricted and unrestricted wavefunctions. The allowed pathway to 2syn is predicted to be greater than 27 kcal mol-1, while the formally forbidden pathway to 2anti is estimated to have a lower barrier of about 23 kcal mol-1. The two transition states for these different pathways are shown in Figure 1, and the sterics that force a helical structure to 1 help make the forbidden pathway more favorable.

TS(1b→2b-syn)

TS(1b→2b-anti)

Figure 1. (U)B3LYP/6-31G optimized geometries of the transition states taking 1 into 2.

Nonetheless, all of the DFT computations significantly overestimate the activation barrier. The authors make the case that a low-lying singlet excited state results in an early conical intersection that reduces the symmetry from C2 to C1. In this lower symmetry pathway, all of the states can mix, leading to a lower barrier. However, since DFT is intrinsically a single Slater configuration, the mixing of the other states is not accounted for, leading to the overestimated barrier height.

In a follow up study, this group examined the thermal and photo cyclization of 13,14-dimethylcethrene 4.2 The added methyl groups make the centhrene backbone more helical, and this precludes the formal allowed disrotatory process. The methyl groups also prohibit the oxidation that occurs with 1, driven by aromatization, allowing for the isolation of the direct product of the cyclization 5. This anti stereochemistry is confirmed by NMR and x-ray crystallography. The interconversion between 4 and 5 can be controlled by heat and light, making the system an interesting photoswitch.

Also of interest is the singlet-triplet gap of 1 and 4. The DFT computed ΔEST is about 10 kcal mol-1 for 4, larger than the computed value of 6 kcal mol-1 for 1b. The EPR of 1b does show a signal while that of 4 has no signal. To assess the role of the methyl group, they computed the singlet triplet gaps for 1b and 4 at two different geometries: where the distance between the carbons bearing the methyl groups is that in 1b (3.03 Å) and in 4 (3.37 Å). The lengthening of this distance by the methyl substituents is due to increased helical twist in 4 than in 1b. For 1b, the gap increases with twisting, from 7.1 to 8.3 kcal mol-1, while for 4 the gap increases by 1.8 kcal mol-1 with the increased twisting. This change is less than the effect of methyl substitution, which increases the gap by 2.2 kcal mol-1 at the shorter distance and 2.8 kcal mol-1 at the longer distance. Thus, the electronic (orbital) effect of methyl substitution affects the singlet-triplet gap more than the geometric twisting.

Interesting 18 π-electron systems involving cyclooctadecanonenetriyne rings have been synthesized and examined by computations.1 The mono-, di- and tri-C18
ring compounds 1, 2, and 3 were prepared and the x-ray structure of 2 was obtained. The B3PW91/6-31G(d,p) optimized geometries of 1-3 and of the tetra ring 4 are shown in Figure 1.

1

2

3

4

Figure 1. B3PW91/6-31G(d,p) optimized geometries of 1-4.

Since the rings are composed of 18 π-electrons in the π-system perpendicular to the nearly planar ring, the natural question is to wonder if the ring is aromatic. The authors computed NICS(0) and NICS(1) values at the center of the C18 rings. For all four compounds, both the NICS(0) and NICS(1) values are negative, ranging from -12.4 to -14.9 ppm, indicating that the rings are aromatic.

They find that vinyl triflates and cyclic vinyl triflates will react with [Ph3C]+[HCB11Cl11]– and triethylsilane to generate vinyl cations that can then be trapped through a C-H insertion reaction. For example, cyclohexenyl triflate 1 reacts in a cyclohexane solvent to give the insertion product 2.

The reactions of isomers 3 and 4 give different ratios of the two products 5 and 6. In both cases, the cyclohexyl is trapped predominantly at the site of the triflate substituent. This means that the mechanism cannot involve a cyclohexene intermediate, since then the two ratios should be identical.

They performed molecular dynamic trajectory analysis at the M062X/6-311+G(d,p) level, starting with the two transition states leading from 3 (TS3) and 4 (TS4), the only transition states located for the insertion reaction. The structures of these TSs are shown in Figure 1.

TS3

TS4

Figure 1. M062X/6-311+G(d,p) optimized geometries of TS3 and TS4.

The trajectories end up in two product basins associated with 5 and 6 starting with either TS3 or TS4. Thus, these transition states are ambimodal, and typical of reactions where dynamic effects dominate. For the reaction of 3, the majority of the trajectories starting at TS3 end up as 5, consistent with the experiments. Similarly, for the trajectories that start at TS4, the majority end up as 6, consistent with experiments.

Once again, we see that relatively simple organic reactions do not follow simple reaction mechanisms, that a single transition state leads to two different products and the product distributions are dependent on reaction dynamics. This may not be too surprising for the vinyl cation insertions given the many examples provide by the Tantillo group of cation rearrangements that are controlled by reaction dynamics (see for examples, this post and this post).

I have recently been interested in curved aromatic systems – see my own paper on double helicenes.1 In this post, I cover four recent papers that discuss non-planar aromatic molecules.

The first paper2 discusses the warped aromatic 1 built off of the scaffold of depleiadene 3. The crystal structure of 1 shows the molecule to be a saddle with near C2v symmetry. B3LYP/6-31G computations indicate that the saddle isomer is 10.5 kcal mol-1 more stable than the twisted isomer, and the barrier between them is 16.0 kcal mol-1, with a twisted saddle intermediate as well.

The PES is significantly simpler for the structure lacking the t-butyl groups, 2. The B3LYP/6-31G PES of 2 has the saddle as the transition state interconverting mirror images of the twisted saddle isomer, and this barrier is only 1.8 kcal mol-1. Figure 1 displays the twisted saddle and the saddle transition state. Clearly, the t-butyl groups significantly alter the flexibility of this C86 aromatic surface. One should be somewhat concerned about the small basis set employed here, especially lacking polarization functions, and a functional that lacks dispersion correction. However, the computed geometry of 1 is quite similar to that of the x-ray structure.

2 twisted saddle (ground state)

2 saddle (transition state)

Figure 1. B3LYP/6-31G optimized geometries of the isomer of 2.

The second paper presents 4, a non-planar aromatic based on [8]circulene 6.3 (See this post for a general study of circulenes.) [8]circulene has a tub-shape, but is flexible and can undergo tub-to-tub inversion. The expanded aromatic 4 is found to have a twisted shape in the x-ray crystal structure. A simplified model 5 was computed at B3LYP/6-31G(d) and the twisted isomer is 4.1 kcal mol-1 lower in energy than the saddle (tub) isomer (see Figure 2). The barrier for interconversion of the two isomers is only 6.2 kcal mol-1, indicating a quite labile structure.

The third paper presents a geodesic molecule based on 1,3,5-trisubstitued phenyl repeat units.4 The authors prepared 7, and its x-ray structure shows a saddle-shape. The NMR indicate a molecule that undergoes considerable conformational dynamics. To address this, they did some computations on the methyl analogue 8. The D7h structure is 309 kcal mol-1 above the local energy minimum structure, which is way too high to be accessed at room temperature. PM6 computations identified a TS only 0.6 kcal mol-1 above the saddle ground state. (I performed a PM6 optimization starting from the x-ray structure, which is highly disordered, and the structure obtained is shown in Figure 3. Unfortunately, the authors did not report the optimized coordinates of any structure!)

Figure 3. PM6 optimized structure of 8.

The fourth and last paper describes the aza-buckybowl 9.5 The x-ray crystal structure shows a curved bowl shape with Cs symmetry. NICS(0) values were computed for the parent molecule 10 B3LYP/6-31G(d). These values are shown in Scheme 1 and the geometry is shown in Figure 4. The 6-member rings that surround the azacylopentadienyl ring all have NICS(0) near zero, which suggests significant bond localization.

Scheme 1. NICS(0) values of 10

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

Our understanding of what aromaticity really means is constantly being challenged!

Sarpong and Tantillo have examined the acid-catalyzed Prins/semipinacol rearrangement of hydroxylated pinenes, such as Reaction 1.1

Rxn 1

Interestingly, only the fenchone scaffold products, like 1, are observed and the camphor scaffold products, like 2, are not observed. Cation intermediates are likely, and this means that a primary alkyl shift is taking place in preference to a tertiary alkyl shift, see Scheme 1.

Scheme 1.

Primary alkyl shift

Tertiary alkyl shift

They proposed the following key steps in the reaction mechanism:

ωB97X-D/6-31+G(d,p) computations find a flat surface around cation intermediate 4: the TS leading to 5 and 6 are only 1.3 and 3.3 kcal mol-1, respectively. Since these small barriers are quite susceptible to changes in basis set and functional, and since Tantillo has found many examples of post-transition state bifurcations in cation systems, the authors reasonably decided to conduct molecular dynamics trajectories originating at the TS connecting 3 and 4. The geometries of the critical points are shown in Figure 1.

The trajectory study shows all the usual characteristics of reactions that are under dynamic control. A third of the trajectories show recrossing of the barrier, typical of very flat surfaces. Nearly all of the remaining trajectories led to 5, with only 2 trajectories (~1%) leading to 6. The dynamics are understandable in terms of favoring the primary alkyl shift over the tertiary since a significantly smaller mass needs to move in the former case.

It never hurts to promote one’s science through clever names – think cubane, buckminsterfullerene, bullvalene, etc. Host-guest chemistry is not immune to this cliché too, and this post discusses the latest synthesis and computations of a nano-Saturn; nano-Saturns are a spherical guest molecule captured inside a ring host molecule. I discussed an example of this a number of years ago – the nano-Saturn comprised of C60 fullerene surrounded by [10]cycloparaphenylene.

Yamamoto, Tsurumaki, Wakamatsu, and Toyota have prepared a nano-Saturn complex with the goal of making a flatter ring component.1 The inner planet is modeled again by C60 and the ring is the [24]circulene analogue 1. The x-ray crystal structure of this nano-Saturn complex is shown in Figure 1.

1: R = 2,4,6-tri-iso-propylphenyl2: R = H

Figure 1. X-ray crystal structure of the nano-Saturn complex of 1 with C60.

Variable temperature NMR experiments gave the binding values of ΔH = -18.1 ± 2.3 kJ mol-1 and TΔS = 0.8 ± 2.2 kJ mol-1 at 298 K. To gauge this binding energy, they computed the complex of C60 with the parent compound 2 at B3LYP/6-1G(d)//M05-2X/6-31G(d), unfortunately without publishing the coordinates in the supporting materials. The computed binding enthalpy is ΔH = -50.6 kJ mol-1, but this computation is for the gas phase. The computed structure shows close contacts of 0.29 nm between the fullerene and the C9-proton of the anthracenyl groups, in excellent agreement with the x-ray structure. These weak C-H…π interactions undoubtedly help stabilize the complex, especially given that the fullerene carries a very tiny Mulliken charge of +0.08 e.

I recently posted on a paper proposing 1,2-benzoquinone and related compounds as the diene component for bioorthogonal labeling. Levandowski, Gamache, Murphy, and Houk report on tetrachlorocyclopentadiene ketal 1 as an active ambiphilic diene component.1

1 is sterically congested to diminish self-dimerization and will react with both electron-rich and electron-poor dienes. To test it as an active diene in bioorthogonal labeling applications, they optimized the structures of the transition states at CPCM(water)/M06-2X/6-311+G(d,p)//CPCM(water)/M06-2X/6-31G(d) for the Diels-Alder reaction of 1 with a variety of dienophiles including trans-cyclooctene 2 and endo-bicyclononyne 3. These transition states are shown in Figure 1. The activation free energy is quite low for each: 18.1 kcal mol-1 with 2 and 18.9 kcal mol-1 with 3.

TS(1+2)

TS(1+3)

Figure 1. CPCM(water)/M06-2X/6-31G(d) optimized geometries for the TSs of the reaction of 1 with 2 and 3.

Experiments were successfully run using 1 as a label on a neuropeptide.

Eckhardt and Schreiner have spectroscopically characterized the aminomethylene carbene 1.1 Their characterization rests on IR spectra, with comparison to the computed AE-CCSD(T)/cc-pCVQZ anharmonic vibrational frequencies, and the UV-Vis spectra, with comparison to the computed B3LYP/6–311++G(2d,2p) transitions.

1 can be converted to 2 by photolysis. Interestingly, 1 does not convert to 2 after 5 days on the matrix in the dark. This is in distinct contrast to hydroxycarbene and related other carbene which undergo quantum mechanical tunneling (see this post and this post). Examination of the potential energy surface for the reaction of 1 to 2 at AE-CCSD(T)/cc-pCVQZ (see Figure 1) identifies that the lowest barrier is 45.8 kcal mol-1, about 15 kcal mol-1 larger than the barrier for the hydroxycarbene rearrangement. Additionally, the barrier width for 1 → 2 is 25% larger than for the hydroxycarbenes. Both of these suggest substantially reduced tunneling, and WKB analysis predicts a tunneling half-life of more than a billion years. The stability of 1 is attributed to the strong π-donor ability of nitrogen to the electron-poor carbene. This is reflected in a very short C-N bond (1.27 Å).

Figure 1. Structures and energies of 1 and 2 and the transition states that connect them. The relative energies (kcal mol-1) are computed at AE-CCSD(T)/cc-pCVQZ.

This is a third post in a series dealing with very short or very long distances between atoms. Ishigaki, Shimajiri, Takeda, Katoono, and Suzuki have prepared three related analogues of hexaphenylethane that all have long C-C bonds.1 The idea is to create a core by fusing two adjacent phenyls into a naphylene, and then protect the long C-C bond through a shell made up of large aryl groups, 1. Fusing another 5-member ring opposite to the stretched C-C bond (2) creates a scissor effect that should stretch that bond further, even more so in the unsaturated version 3.

Their M062-x/6-31G* computations predict an increasing longer C-C bond (highlighted in blue in the above drawing): 1.730 Å in 1, 1.767 Å in 2, and 1.771 Å in 3. The structure of 3 is shown in Figure 3.

Figure 1. M06-2x/6-31G(d) optimized structure of 3.

These three compounds were synthesized, and characterized by IR and Raman spectroscopy. Their x-ray crystal structures at 200 K and 400K were also determined. The C-C distances are 1.742 Å (1), 1.773 Å (2) and 1.798 Å (3) with distances slightly longer at 400 K. These rank as the longest C-C bonds recorded.