Does forming a Wheland intermediate disrupt all aromaticity?

Text books will announce that during aromatic electrophilic substitution, aromaticity is lost by the formation of a Wheland intermediate (and regained by eliminating a proton). Is that entirely true?

I will start by considering the simplest of all such intermediates[1],[2] the NMR of which was first reported by Olah and which is shown below (red). The values in black are from a ωB97XD/6-311G(d,p) calculation (chloroform)[3]. The discrepancy might be due to the difference in solvent (HF-SbF5/SO2ClF-SO2F2) and my failure to include a counter-ion in the calculation. Actually, this result is only obtained at -134°C; warming to -70°C increases the rate of [1,2] ring shifts and only one “aromatic” peak at 8.09 ppm is observed.

Click for 3D to see vibrational modes

 The computed MOs are shown below. MOs 21 and 15 are pretty much identical with those for benzene itself (the former because of symmetry),  so relatively little disruption of the conjugation there. An NBO analysis shows that each C-H bond from the sp3 centre acts as a donor to the π-system (E2 = 8.5 kcal/mol), which we know in another guise as hyperconjugation or σ-conjugation. In this sense, the two C-H bonds are acting as surrogates for a real p-AO on that carbon, and the E2 value is certainly not insignificant. Certainly they are nowhere near as good a donor as one real p-AO, but nonetheless, the two of them together are good enough to result in retention of a significant measure of cyclic π-conjugation, and hence probably aromaticity. The C-H bonds are weakened as a result of their electron donation, which reduces their normal mode wavenumber down by about 200 cm-1. The so-called Kekule-mode, which is depressed in benzene itself to about 1300 cm-1 because of the effect termed π-distortivity[4], is actually increased to 1442 cm -1 in the Wheland intermediate.[5] This mode is normally elevated by any weakening of the distortive π-system (in my previous post, I noted such an elevation induced by the quintet excited state of benzene) and so we many presume the same effect operates in the Wheland intermediate.

 

MO 20 Click for 3D

 

MO 21 Click for 3D

MO 15 Click for 3D

To get a measure of any aromaticity deriving from a ring current, I computed the NICS(1) value from the ring centroid. This latter was itself obtained from the coordinates of the ring-critical-point determined by a QTAIM analysis, and computed at that point and 1, 2 and 3Å above it. The NICS values are respectively -0.4, -6.1, -3.6 and -1.6 ppm. The maximum is indeed the NICS(1) point, and at that point, the value indicates modest, but real π-aromaticity.

So does forming a Wheland intermediate disrupt all aromaticity? The answer is a clear no! What might be interesting to compute (I have not tried doing so here) is an actual energy for the stabilisation resulting from the weak aromaticity present as a % of that present in benzene, so that the text-books can be amended.


It might also be that the exchange in the measured NMR spectrum was not entirely suppressed, and so some residual averaging of the peak positions remains.

References

  1. G.A. Olah, R.H. Schlosberg, D.P. Kelly, and G.D. Mateescu, "Stable carbonium ions. IC. Benzenonium ion (C6H7+) and its degenerate rearrangement", Journal of the American Chemical Society, vol. 92, pp. 2546-2548, 1970. http://dx.doi.org/10.1021/ja00711a057
  2. G.A. Olah, R.H. Schlosberg, R.D. Porter, Y.K. Mo, D.P. Kelly, and G.D. Mateescu, "Stable carbocations. CXXIV. Benzenium ion and monoalkylbenzenium ions", Journal of the American Chemical Society, vol. 94, pp. 2034-2043, 1972. http://dx.doi.org/10.1021/ja00761a041
  3. S. Shaik, A. Shurki, D. Danovich, and P.C. Hiberty, "A Different Story of π-DelocalizationThe Distortivity of π-Electrons and Its Chemical Manifestations", Chemical Reviews, vol. 101, pp. 1501-1540, 2001. http://dx.doi.org/10.1021/cr990363l
  4. Henry S. Rzepa., "Gaussian Job Archive for C6H7(1+)", 2013. http://dx.doi.org/10.6084/m9.figshare.870473
Henry Rzepa

Henry Rzepa is Emeritus Professor of Computational Chemistry at Imperial College London.

View Comments

  • The usage of the "ring critical point" comes imho a bit out of the blue, as there is no direct connection between the electron density topology and the the current topolgy _in general_ known as yet. I prefer in such cases the use of the central current stagnation point [r_0 for which j(r_0) = 0].

    Finding a general relation between the two would be high interest!

  • Re: "I prefer in such cases the use of the central current stagnation point [r_0 for which j(r_0) = 0]."

    Yes, that sounds a more pertinent point. But one does need specialised programs to compute this; QTAIM is rather more generally available.

    Perhaps you can suggest an easily available program for computing current stagnation points?

  • GIMIC interfaces with Turbomole and ACESII.
    Get it like that:

    git clone https://github.com/bakerjonas/gimic.git
    cd gimic
    ./setup
    cd build
    make
    sudo make install

    You can calculate vector fields like that one shown in the picture above. Searching the stagnation points "by hand" using "nested intervals" works pretty well for me1. Advanced computer users can also write some Newton-Raphson optimizer script for the inverse modulus of j (only a few lines in python f. i.).

  • Paul notes his 1999 investigation of cyclopentadiene, in which the two methylene hydrogens are replaced by SnR3 groups. To link into the above, I have done the same, but replacing the hydrogens with two Lithium atoms (doi: 10.6084/m9.figshare.872568 ). The resulting structure has all the C-C bond lengths almost equal, but a more interesting partitioning of the four electrons nominally associated with the two "C-Li" bonds. Two enter the ring to achieve 6-π-electron aromaticity for the ring, the other two form an in-plane σ-carbene. This species might be a very strong base as a result!

    The three highest occupied π-orbitals again strongly resemble those of benzene. The ELF analysis (below) shows the ring to contain five ~3-electron basins and the carbene lone pair (label 15, 2.6 electrons).

  • There is also a handfull of other programms capable of calculating currents. Keith Todd's AIMALL being one of them.

    This and othe issues concerning molecular currents are discussed in our facebook group called "Magnetically Induced Electronic Currents in Molecules". Anyone interested is cordially invited to join the group. Among our members are Jeppe Olsen, Trygve Helgaker, Kenneth Ruud, Jochen Autschbach and many more.

  • Here is another analogue, with BH2 replacing CH2 in cyclopentadiene, forming an anion. It too has some measure of apparent aromaticity (but less than some of the other systems discussed above).

  • What trends in terms of aromaticity are to be expected in the series cyclpropene, cp, cycloheptatriene, ... ?

  • Re: "What trends". Well, nominally, hyperconjugation should be suppressed in cyclopropene (4n anti-aromaticity), encouraged in cyclopentadiene (see above, and Paul's comment) and suppressed in cycloheptatriene. But the latter is bound to be non-planar, so it is not a fair progression.

    Certainly the little tricks played above (BH2 for CH2, or CLi2 for CH2) should not succeed for cyclopropene.

Recent Posts

Detecting anomeric effects in tetrahedral boron bearing four oxygen substituents.

In an earlier post, I discussed a phenomenon known as the "anomeric effect" exhibited by…

1 week ago

Internet Archeology: reviving a 2001 article published in the Internet Journal of Chemistry.

In the mid to late 1990s as the Web developed, it was becoming more obvious…

2 months ago

Detecting anomeric effects in tetrahedral carbon bearing four oxygen substituents.

I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…

2 months ago

Data Citation – a snapshot of the chemical landscape.

The recent release of the DataCite Data Citation corpus, which has the stated aim of…

2 months ago

Mechanistic templates computed for the Grubbs alkene-metathesis reaction.

Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…

3 months ago

3D Molecular model visualisation: 3 Million atoms +

In the late 1980s, as I recollected here the equipment needed for real time molecular…

3 months ago