Summary: Conformations of the mono, di, tri and tetra hetero [12] annulenes which have C2 or pseudo C2 symmetry are explored using ab initio theory for their potential Möbius aromatic properties, as quantified using the nucleus independent chemical shift (NICS) values at the ring centroid. This is found to be greatest for the sulphur heterocycles, and a number of known examples of tri and tetrahetero systems show C2 symmetric conformations for the crystal structures and modestly aromatic NICS values. The structures of several other annulenes show distorsions resulting from Möbius aromaticity, most promimently the tetrathiocine and tetraselenocine ring systems. We also propose forms of the elements nitrogen and sulfur which have some Möbius-aromatic character, the latter in particular being only 11 kcal mol-1 higher than the most stable conformation.
Heilbronner1 proposed in 1964 that a 180o phase shift in the pp-pp overlap of a cyclic array of such orbitals occupied by 4n p electrons and delocalised over the entire ring could result in a so-called Möbius conformation having aromatic properties, rather than the Hückel antiaromaticity associated with planarity. A characteristic feature noted subsequently is the chirality of these systems, often associated with a formal or approximate C2 axis of symmetry passing through the aromatic ring. In 1998, Schleyer and co-workers2 reported that the lowest energy conformation of the known conjugated 4n pp C9H9+ cation was of this type, and we followed this with the suggestion of a Möbius form with C2 symmetry for the neutral [16] annulene3. However, this conformation did not exhibit any significant aromaticity, as quantified using the theoretical NICS (nucleus independent chemical shift) probe.4 More aromatic Möbius conformations for [12], [16] and [20] annulenes were subsequently published by Schleyer and co-workers5, and although their predicted energies were around 6 kcal mol-1 higher than the non-aromatic forms corresponding to the X-Ray structure of e.g. [16] annulene, it was suggested that these conformations might nevertheless play a role in the chemistry of [4n]annulenes. We had earlier shown6 that Möbius minima can also be located for heteroatom derivatives of [8] annulenes, but these were predicted to be > 20 kcal mol-1 higher than non-Möbius forms. To date therefore, no Möbius aromatic neutral annulene rings have been identified which are more stable than their non or anti-aromatic isomers, and hence capable of structural characterisation and study. We report here our findings for 12p rings, including those derived from [12] annulene by replacing between one to four double bonds with heteroatoms (1-6), along with several other heteroannulene ring systems.
Molecular structures were normally optimised without symmetry constraints using Gaussian987 at the closed shell B3LYP/6-31G(d) level. The aromaticity of these systems was quantified using a NICS4 probe placed at the ring centroid. Coordinates are available via the electronic supplemental data (ESI) at http://www.rsc.org/ESI/...
The hetero[12]annulene series reveals several entries for compounds of the general type 1-5. In contrast to a number of crystalline examples of 7-ring heteropine [8]annulenes, only a single example of the [12] analogue 1, X=S, has been reported,14 but with a methano-bridge constricting it to a plane of symmetry. Likewise one example of 2 is known; an oxastibecene (X=O, Y=SbPh), although no 3D structure is reported. More examples occur for the ring system3 including the benzo-substituted examples X=Y=O13 and X=Y=S13 which exhibit a formal C2 symmetry axis rather than a plane. One example containing X=S, Y=N-acetyl15 can be identified as having a pseudo C2 symmetry axis which passes through one S and the opposing C=C bond. These latter species can all be described as having conformations formed by replacing three C=C units in the recently identified Möbius form of [12] annulene5 with the corresponding heteroatom.
There are also several interesting heteroatom rich systems. Six examples of the 12p ring system 4 (X=S, a tetrathiocin) are known17 to exhibit chiral D2 ring symmetry, although occurances of this ring with Cs symmetry are more common. The only known selenium analog 4 (X=Se, tetraselenocine) exhibits D2 ring symmetry,18 as do all the 4,5,6,7 isomers of tetrathiocin (5, X=S) which are also not rare.19 The unusual hexa-aza[12] annulene 616 has distinct C2 like distorsions from planarity.
We conclude that a search of the Cambridge database does seem to reveal tantalising evidence of stable conformations of [12] heteroannulenes which could exibit to some degree Möbius characteristics. To understand how the geometrical preferences of these systems arise, we turned to ab initio modelling at the B3LYP/6-31G(d) theoretical level.
For all three heteroatoms an unsymmetric (C1) conformation was the most stable, although the C2 isomers were significantly more accessible in energy (4.0 kcal mol-1 higher for X=S, 7.4 for X=O, 12.0 for X=N) than we had previously observed for the smaller ring [8] heteropines.6 The NICS(0) value for the C2 Möbius forms continued the trend of increased aromaticity along the series 1-3, as do the C=C ring bond lengths which reveal a small degree of increasing delocalisation (1-3 X=S; 1.339, 1.343 and 1.346/1.355/1.363Å respectively).
For the homodihetero systems, two C2 symmetric conformations were located, the slightly lower having an axis perpendicular to the approximate plane of the ring and the other passing directly through the two heteroatoms. The Cs symmetric conformation shows all cis alkenes and the plane again passing through the two heteroatoms. With X=O or NH, the latter conformation was the more stable, whereas with X=S, the C2 form was. The NICS(0) of the latter in particular corresponds to a modest diatropic ring current corresponding to mild aromaticity.
The conformational surface of these 9-membered rings is relatively complex, and no attempt here has been made to identify all the possible conformations. We instead concentrate on comparing those with a proper or psuedo C2 symmetry axis with conformations with a Cs plane of symmetry. Where searches have revealed conformations with no (C1) symmetry these are also included (Table 2). Thus the triaza system 1, X=Y=NH reveals a number of significant conformations (Table 2), the lowest in energy of which has an approximate C2 symmetry axis passing through one NH bond and the mid point of the opposite C-C bond, broken by the orientation of the two remaining NH bonds, respectively axial and equatorial with respect to the ring itself. Isomers with two axial or two equatorial NH bonds exhibit proper C2 symmetry but are 2.4-3.0 kcal mol-1 higher in energy. A Cs symmetric isomer with cis C=C double bonds is 11 kcal mol-1 higher in energy; a form with two trans C=C double bonds being 41.6 higher. Despite the greater stability of the C2 forms, the NICS(0) value at the ring centroid of the most stable conformation (+1.3 ppm) indicates neither aromatic or anti-aromatic character. No examples of this simple ring system have been hitherto reported, but experimental characterisation does appear desirable. The phosphorus analogue shows a subtle difference to nitrogen. The more stable form reveals a pseudo C2 axis passing from a C=C mid point through an axis defined by the P atom and its lone pair of electrons, indicating perhaps a lower degree of conjugation of this lone pair to the ring. Again, this ring is not known experimentally.
The analogous oxygen heterocycle (X=Y=O) is conformationally simpler with the C2 form 6.4 kcal mol-1 lower than Cs, and this mode of course corresponds to the observed crystal structure of the tribenzo derivative.13 As with the triaza heterocycle, the NICS(0) value indicates no aromaticity. We have previously suggested3 that ring fluorination has the effect of enhancing Möbius characteristics, and in this instance the perfluoro derivative does have a more aromatic NICS(0) value, but at the expense of stability, since this conformer is now 3.0 kcal mol-1 less stable than the corresponding Cs form. The trithia heterocycle (X=Y=S) is the most interesting of the series. Attempts to locate the all cis alkene Cs form all converged to a C2 isomer; the only Cs symmetric form that could be located contained two trans alkenes and was 18.8 kcal mol-1 higher in energy. As with oxygen, perfluorination tends to a more aromatic NICS(0) but less stable C2 isomer. Here again a known tribenzo derivative is known to adopt the C2 form.14
We also explored several other variations of X or Y. Replacing Y=O by Y=CH- forms a charged dihetero ring for which the C2 form is no longer the more stable (by 3.1 kcal mol-1), but for this system perfluorination makes the C2 form more stable by 1.7 kcal mol-1. The system X=O, Y=S has two C2-like isomers differing by 2.0 kcal mol-1 and both were much lower than the only Cs form locatable. With X=S, Y=CH-, the stable conformation was again C2, with ring perfluorination reducing the NICS(0) value significantly. The system X=S, Y=NR is of interest, since an example of this ring is known.15 With Y=NH, a pseudo C2 structure corresponding to the crystal structure can be located; a Cs in this instance is computed to be 3.2 kcal mol-1 more stable, but may well be the result of NH...S hydrogen bonding in the latter orientation. With X=S, Y=O, again C2-like isomers differ by 1.6 kcal mol-1, with the two trans C=C Cs form some 26.4 kcal mol-1 higher.
In the trihetero ring system therefore, we see evidence of the preference for conformations with a C2 rather than a Cs symmetry, and this preference does seem to be enhanced by sulfur in particular. The magnetic properties of these rings as quantified by the NICS(0) probe do not reveal any prominent aromaticity however, but as 4n (n=3) p systems, they are clearly not Hückel antiaromatics.
Species 6 is a heteroatom rich [12] annulene with the unusual out of plane C2-like distorsion as noted above. The calculations confirm this, and indicate also modest aromaticity (Table 2). This interpretation however must be moderated by the potential "trannulene" effect as described by Schleyer and co-workers22 which may also contribute to the diatropic ring current. The tetrathiocine ring 9, a=1, b=0 has a slightly aromatic NICS value, and this is more true of the alkene equivalent 7a, X=S, which also reveals a 5.7 kcal mol-1 preference for D2 symmetry. The molecular orbitals (Figure 1) for 7a, X=S are characteristic of systems with Möbius character3, revealing a catenane-like phase structure for the electronic Möbius strip. Ring fluorination in this case decreases the aromaticity and the stability of the Möbius form. The alternate isomers 7b, X=S are in absolute energies lower than 7a, X=S, and show similar mild aromaticity. The selenium heterocycle 7a, X=Se, R=H is predicted to be more twisted and with greater aromaticity than the sulfur analogue, an effect remarked upon previously in the context of spiroaromatic systems for the series N, P, As.23. The known18 example (7a, X=Se, R=CF3) has a measured C=C bond length of 1.385Å, which is approaching the fully aromatic delocalised value in e.g. benzene. The NICS value of -7.6 indicates this system is the most aromatic of the various discussed. The equivalent bond length in 7a, X=S is 1.33Å. The homologue tetrathia[20]annulene 9, a=2,b=1 reveals a strong D2 twist in common with the series 4, but only a very slightly negative NICS value. This may in part be because the distance from the ring centroid to the nearest carbon atom is 2.64Å, which may be too far to induce any significant magnetic effect. The molecular orbitals of 9, a=2,b=1 have characteristically Möbius photogenic form!
One can remove all the carbon atoms entirely from the ring, resulting in S8, for which a C2 symmetric form 10 equivalent to a [16] annulene can be readily located, being only 11 kcal mol-1 higher than the known stable form. This Möbius sulfur has modest NICS aromaticity, and again characteristic molecular orbitals (Figure 1). Continuing the removal of carbon atoms, the three benzene rings in 6 can be replaced by three N=N units, giving the N12 allotrope 11. Unlike 6, it exhibits a true C2 axis and modest aromaticity. Its energy relative to gaseous nitrogen (353 kcal mol-1) suggests this is a purely theoretical species (we do not urge its synthesis!) but it does serve to illustrate the generality of the Möbius distorsion in the [12] and higher annulenes.
The prototypic molecule from which the concept of aromaticity and anti-aromaticity is largely derived is the monocyclic annulene, and although the suggestion that a Möbius aromatic form might be possible was made almost forty years ago, experimental characterisation of any such annulenes has hitherto been largely unaccomplished or unrecognised. Properties such as the molecular geometry, NMR characteristics and reactivity have thus far only been inferred by theoretical calculation. Following the recent suggestion that low energy Möbius conformations of medium sized 4np annulenes may be accessible, we here reinforce this by drawing attention to various hetero-[12]annulenes, particularly those containing sulfur, whose Möbius characteristics in stable (C2 symmetric) conformations may not have hitherto been recognised. We feel it appropriate to conclude that the prospects for the synthesis of new stable forms of Möbius aromatic molecules must be good.
| Table 2. Calculated energies (B3LYP/6-31G(d), Hartree) or relative energies(kcal mol-1) for heteroannulenes 1-7 and NICS(0) values (ppm). | ||||
|---|---|---|---|---|
| X, Y, Ra | Relative Energy of Cs (C1) [Planar] | Total Energy of C2 (pseudo C2) | NICS(0) [Planar] | |
| 1, X=NH, R=H | -8.8 (-12.0) | -442.31797 | -5.4 | |
| 1, X=O, R=H | -1.5 (-7.4) | -462.17742 | -5.6 | |
| 1, X=S, R=H | 1.2 (-4.0) | -785.15962 | -6.3 | |
| 1, X=S, R=F | -1.1 (-3.1) | -1777.42173 | -4.8 | |
| 2, X=Y=N | -11.3 | 5.3, -420.27282 | -2.9, -1.2 | |
| 2, X=Y=O | -2.8 | 8.4, -459.99477 | -0.1, -0.6 | |
| 2, X=Y=S | 3.2 | 1.8, -1105.94349 | -4.8, -1.4 | |
| 2, X=Y=S, R=F | -10.2 | -1899.74067 | -4.6 | |
| 2, X=O, Y=SbH, R=H | -390.75207 | -390.72924 | -0.6 | |
| 2, X=O, Y=SbF, R=F |
-1283.91234 | -1283.91401 |
-2.7 | |
| 3, X=Y=NH, R=H | 13.5, [-398.17027, -397.90839] | 2.5, 3.1 (-398.20690) | 1.3 [not calculated, +17.3] | |
| 3, X=Y=O, R=H | 6.4, [-457.73639] | -457.78250 | 1.5 [+7.7] |
|
| 3, X=Y=O, R=F | -3.1 | -1053.15767 | -2.3 | |
| 3, X=Y=PH, R=H | 4.3 | (-1258.00411) | 0.1 |
|
| 3, X=Y=S, R=H | 6.7, [-1426.65943] | -1426.72236 | 0.0 [+8.1] |
|
| 3, X=Y=S, R=F | -5.1 | -2022.06902 | -3.1 | |
| 3, X=O, Y=CH-, R=H | -3.2 | -421.26886 | 2.8 |
|
| 3, X=O, Y=CF-, R=F | 1.7 | -1115.91871 | -1.0 | |
| 3, X=O, Y=S, R=H | 5.6 | -780.76745 (2.0) | 0.7 |
|
| 3, X=O, Y=S, R=F | 28.2 | (-1376.13160) | -2.8 | |
| 3, X=S, Y=CH-, R=H | 9.1 | -1067.25508 | 5.3 |
|
| 3, X=S, Y=CF-, R=F | 4.4 | -1761.88127 | -0.2 | |
| 3, X=S, Y=NH, R=H | -3.2 | ( -1083.88695) | 0.2 | |
| 3, X=S, Y=O, R=H | 26.4 | -1103.74840 (1.6) | 0.0 |
|
| 3, X=S, Y=O, R=F | -3.8 | -1699.10275 | -3.3 | |
| 3, X=NH, Y=S, R=H | -741.04548 | -741.03715, -741.04513,-741.04512 | 0.9, -1.7, 0.2 | |
| 3, X=NF, Y=S, R=F |
-1534.73962 |
-1534.72736 | -3.1 |
|
| 4, X=S, R=H | -1747.52621 | -1747.53529 | -4.7 | |
| 5, X=S, R=H | - | -1747.55077 | -3.3 | |
| 4, X=S, R=F | -2144.43257 | -2144.42289 | -3.9 | |
| 4, X=Se, R=H | -9752.33901 | -9752.34808 | -6.2 | |
| 5, X=Se, R=H | (-9752.37322) | -9752.37739 | -1.6 | |
| 4, X=Se, R=F | -10149.251472 | -10149.24279 | -4.9 | |
| 4, X=Se, R=CF3 | - | -11100.4724 | -7.6 | |
| 6 | - | (-1021.50351) | -5.9 | |
| 7 | - | (-656.58174) | -4.4 | |