Möbius aromaticity

In organic chemistry, Möbius aromaticity is a special type of aromaticity believed to exist in a number of organic molecules. In terms of molecular orbital theory these compounds have in common a monocyclic array of molecular orbitals in which there is an odd number of out-of-phase overlaps, the opposite pattern compared to the aromatic character to Hückel systems. The nodal plane of the orbitals, viewed as a ribbon, is a Möbius strip, rather than a cylinder, hence the name. The pattern of orbital energies is given by a rotated Frost circle (with the edge of the polygon on the bottom instead of a vertex), so systems with 4n electrons are aromatic, while those with 4n + 2 electrons are anti-aromatic/non-aromatic. Due to incrementally twisted nature of the orbitals of a Möbius aromatic system, stable Möbius aromatic molecules need to contain at least 8 electrons, although 4 electron Möbius aromatic transition states are well known in the context of the Dewar-Zimmerman framework for pericyclic reactions. Möbius molecular systems were considered in 1964 by Edgar Heilbronner by application of the Hückel method, but the first such isolable compound was not synthesized until 2003 by the group of Rainer Herges. However, the fleeting trans-C9H9+ cation, one conformation of which is shown on the right, was proposed to be a Möbius aromatic reactive intermediate in 1998 based on computational and experimental evidence.

Hückel-Möbius aromaticity
The Herges compound (6 in the image below) was synthesized in several photochemical cycloaddition reactions from tetradehydrodianthracene 1 and the ladderane syn-tricyclooctadiene 2 as a substitute for cyclooctatetraene.
 * [[Image:MobiusAromaticity.png|400px|Synthesis of a Mobius ring]]

Intermediate 5 was a mixture of 2 isomers and the final product 6 a mixture of 5 isomers with different cis and trans configurations. One of them was found to have a C2 molecular symmetry corresponding to a Möbius aromatic and another Hückel isomer was found with Cs symmetry. Despite having 16 electrons in its pi system (making it a 4n antiaromatic compound) the Heilbronner prediction was borne out because according to Herges the Möbius compound was found to have aromatic properties. With bond lengths deduced from X-ray crystallography a HOMA value was obtained of 0.50 (for the polyene part alone) and 0.35 for the whole compound which qualifies it as a moderate aromat.

It was pointed out by Henry Rzepa that the conversion of intermediate 5 to 6 can proceed by either a Hückel or a Möbius transition state.


 * [[Image:HuckelvsMobiusTransitionState.png|400px|Huckel vs Mobius Transition State]]

The difference was demonstrated in a hypothetical pericyclic ring opening reaction to cyclododecahexaene. The Hückel TS (left) involves 6 electrons (arrow pushing in red) with Cs molecular symmetry conserved throughout the reaction. The ring opening is disrotatory and suprafacial and both bond length alternation and NICS values indicate that the 6 membered ring is aromatic. The Möbius TS with 8 electrons on the other hand has lower computed activation energy and is characterized by C2 symmetry, a conrotatory and antarafacial ring opening and 8-membered ring aromaticity.

Another interesting system is the cyclononatetraenyl cation explored for over 30 years by Paul v. R. Schleyer et al. This reactive intermediate is implied in the solvolysis of the bicyclic chloride 9-deutero-9'-chlorobicyclo[6.1.0]-nonatriene 1 to the indene dihydroindenol 4. The starting chloride is deuterated in only one position but in the final product deuterium is distributed at every available position. This observation is explained by invoking a twisted 8-electron cyclononatetraenyl cation 2 for which a NICS value of -13.4 (outsmarting benzene) is calculated. A more recent study, however, suggests that the stability of trans-C9H9+ is not much different in energy compared to a Hückel topology isomer. The same study suggested that for [13]annulenyl cation, the Möbius topology penta-trans-C13H13+ is a global energy minimum and predicts that it may be directly observable.


 * Cyclotetranonaene cation.png[[Image:CyclononatetraenylCation.png|400px|the cyclononatetraenyl cation]]

In 2005 the same P. v. R. Schleyer questioned the 2003 Herges claim: he analyzed the same crystallographic data and concluded that there was indeed a large degree of bond length alternation resulting in a HOMA value of -0.02, a computed NICS value of -3.4 ppm also did not point towards aromaticity and (also inferred from a computer model) steric strain would prevent effective pi-orbital overlap.

A Hückel-Möbius aromaticity switch (2007) has been described based on a 28 pi-electron porphyrin system:
 * [[Image:HuckelMobiusSwitch.svg|Hückel-Möbius aromaticity switch]]

The phenylene rings in this molecule are free to rotate forming a set of conformers: one with a Möbius half-twist and another with a Hückel double-twist (a figure-eight configuration) of roughly equal energy.

In 2014, Zhu and Xia (with the help of Schleyer) synthesized a planar Möbius system that consisted of two pentene rings connected with an osmium atom. They formed derivatives where osmium had 16 and 18 electrons and determined that Craig–Möbius aromaticity is more important for the stabilization of the molecule than the metal's electron count.

Transition states
In contrast to the rarity of Möbius aromatic ground state molecular systems, there are many examples of pericyclic transition states that exhibit Möbius aromaticity. The classification of a pericyclic transition state as either Möbius or Hückel topology determines whether 4N or 4N + 2 electrons are required to make the transition state aromatic or antiaromatic, and therefore, allowed or forbidden, respectively. Based on the energy level diagrams derived from Hückel MO theory, (4N + 2)-electron Hückel and (4N)-electron Möbius transition states are aromatic and allowed, while (4N + 2)-electron Möbius and (4N)-electron Hückel transition states are antiaromatic and forbidden. This is the basic premise of the Möbius-Hückel concept.

Derivation of Hückel MO theory energy levels for Möbius topology
From the figure above, it can also be seen that the interaction between two consecutive $$p_z$$ AOs is attenuated by the incremental twisting between orbitals by $$\cos\omega$$, where $$\omega=\pi/N$$ is the angle of twisting between consecutive orbitals, compared to the usual Hückel system. For this reason resonance integral $$\beta^\prime$$ is given by
 * $$\beta^\prime=\beta\cos(\pi/N)$$,

where $$\beta$$ is the standard Hückel resonance integral value (with completely parallel orbitals). Nevertheless, after going all the way around, the Nth and 1st orbitals are almost completely out of phase. (If the twisting were to continue after the $$N$$th orbital, the $$(N+1)$$st orbital would be exactly phase-inverted compared to the 1st orbital). For this reason, in the Hückel matrix the resonance integral between carbon $$1$$ and $$N$$ is $$-\beta^\prime$$.

For the generic $$N$$ carbon Möbius system, the Hamiltonian matrix $$\mathbf{H}$$ is:

\mathbf{H}= \begin{pmatrix} \alpha& \beta' & 0 &\cdots& -\beta' \\ \beta' & \alpha& \beta' & \cdots & 0 \\ 0 & \beta' & \alpha & \cdots & 0 \\ \vdots &\vdots &\vdots &\ddots &\vdots \\ -\beta' &0& 0 & \cdots & \alpha \end{pmatrix} $$. Eigenvalues for this matrix can now be found, which correspond to the energy levels of the Möbius system. Since $$\mathbf{H}$$ is a $$N\times N$$ matrix, we will have $$N$$ eigenvalues $$E_k$$ and $$N$$ MOs. Defining the variable
 * $$x_k=\frac{\alpha-E_k}{\beta'}$$,

we have:

\begin{pmatrix} x_k& 1 & 0 &\cdots& -1 \\ 1 & x_k& 1 & \cdots & 0 \\ 0 & 1 & x_k & \cdots & 0 \\ \vdots &\vdots &\vdots &\ddots &\vdots \\ -1 &0& 0 & \cdots & x_k \end{pmatrix} \cdot \begin{pmatrix} c_1^{(k)} \\ c_2^{(k)} \\ c_3^{(k)} \\ \vdots\\ c_N^{(k)} \\ \end{pmatrix}=0 $$. To find nontrivial solutions to this equation, we set the determinant of this matrix to zero to obtain

x_k=-2\cos{\frac{(2k+1)\pi}{N}} $$. Hence, we find the energy levels for a cyclic system with Möbius topology,

E_k=\alpha+2\beta^\prime\cos{\frac{(2k+1)\pi}{N}}\quad (k=0,1,\ldots, \lceil N/2 \rceil-1) $$. In contrast, recall the energy levels for a cyclic system with Hückel topology,

E_k=\alpha+2\beta\cos{\frac{2k\pi}{N}}\quad (k=0,1,\ldots, \lfloor N/2 \rfloor) $$.