Structures of Annulenes and Model Annulene Systems in the Ground and Lowest Excited States

The paper introduces general considerations on structural properties of aromatic, antiaromatic and non-aromatic conjugated systems in terms of potential energy along bond length alternation and distortion coordinates, taking as examples benzene, cyclobutadiene and cyclooctatetraene. Pentalene, formally derived from cyclooctatetraene by cross linking, is also considered as a typical antiaromatic system. The main interest is concerned with [n]annulenes and model [n]annulene molecular systems, n ranging from 10 to 18. The rich variety of conformational and configurational isomers and of dynamical processes among them is described. Specific attention is devoted to bridged [10]and [14]annulenes in the ground and lowest excited states as well as to s-indacene and biphenylene. Experimental data obtained from vibrational and electronic spectroscopies are discussed and compared with ab initio calculation results. Finally, porphyrin, tetraoxaporphyrin dication and diprotonated porphyrin are presented as annulene structures adopting planar/non-planar geometries depending on the steric hindrance in the inner macrocycle ring. Radiative and non-radiative relaxation processes from excited state levels have been observed by means of time-resolved fluorescence and femtosecond transient absorption spectroscopy. A short account is also given of porphycene, the structural isomer of porphyrin, and of porphycene properties.


Introduction
Annulenes, i.e., conjugated hydrocarbon monocycles of general formula C 2n H 2n , form a class of organic molecules for which the experimental ingenuity of synthetic chemistry [1][2][3] and the internal robustness of quantum chemistry [4,5] have fully displayed their best qualities. For instance, benzene has represented along the years one of the preferred reference systems in order to validate new ab initio theoretical procedures. Higher annulenes, on the other hand, have allowed to investigate about isomeric rearrangements by means of advanced experimental techniques [6] and theoretical methods [7][8][9]. Comprehensive treatments of annulene preparation, structure and properties are available in textbooks [4,10,11], thematic issues [12,13] and review articles [14,15].
According to the Hückel rule conjugated monocycles with (4n + 2) π electrons are more stable than the corresponding open-chain polyenes and should be associated with a set of structural and chemical properties (bond equalized geometries, substitution reactivity) collectively known as aromatic [4,10]. However, the experimental difficulties under which the next higher homologues of benzene, [10]-and [14]annulene, have been synthesized [1] and the ease of thermal rearrangements [6] are in stark contrast with the stability of the first component of the (4n + 2) series. Additional steric factors related to the cycle geometry and size not explicitly considered by the Hückel theory obscure the energetic and structural manifestations of aromaticity. Therefrom two main strategies have been apparently developed. On one hand, aromaticity has been recovered lowering as much as possible the molecular flexibility by the introduction of chemical bridges between opposite ends of the annulene ring [16] or by replacement of double with triple bonds along the cycle [1]. On the other, the theoretical foundations responsible for the delocalized arrangement of the π electrons in aromatic systems have been critically reexamined leading to the view of π propensity to distortion in benzene and cyclobutadiene [17]. In addition, it has been found that the transition from bond equalized to bond alternating structures in (4n + 2) π systems, theoretically predicted long ago [18], may occur for rings of 30 C atoms [19] or lower [20].
Monocycles with 4n π electrons are the second large group of conjugated hydrocarbons, the antiaromatic systems [21,22], characterized by the instability with respect to acyclic polyenes and the tendency toward bond alternating geometries [4]. Cyclobutadiene has been considered the reference antiaromatic system [23,24] though recently criticism has been advanced to this view [25,26]. Going up with the ring size, while cyclooctatetraene is characterized by a non planar arrangement of the four π bonds [4,10,11], pentalene, where the ring is forced to planarity by crosslink [27,28], is equally well suited to model antiaromatic properties [29]. On the contrary, the more complex and formally antiaromatic structure of s-indacene peculiarly shows a bond equalized ring geometry [30,31].
It is the purpose of this review to report on the experimental and computational activity of our group in the field of aromatic as well as antiaromatic systems in the last years with particular attention to molecular structures of the ground and lowest excited states. Further, it has been thought to be useful to associate to our results an updated outline of the aromaticity/antiaromaticity theory [17] and a short description of isomeric structures up to [18]annulene [10,15,32,33]. The paper is organized in two Sections, one dedicated to basic aromatic/antiaromatic systems and the second to higher annulenes and model annulene systems to which the research activity of our group has mostly contributed.

Archetypal Conjugated Monocycles: Cyclobutadiene, Benzene and Cyclooctatetraene
This Section deals with the three simplest monocyclic conjugated hydrocarbons, cyclobutadiene, benzene and cyclooctatetraene. These molecules may be considered as reference systems to introduce the structural properties of ring geometries characterizing also higher annulenes. Review articles on cyclobutadiene [23,24,26], benzene [34,35] and cycloctatetraene [36,37] are well known references of the past years.

π-electron Approximation
In the π-electron approximation both the Hückel molecular orbital (HMO) method and the simple valence bond (VB) theory predict that displacing the planar cyclic structures from the bond equalized to the bond alternating configurations the ground state energy is lowered [38,39].
For instance, in the VB treatment of benzene [38] the two equivalent Kekule' structures K A and K B at the bond equalized D 6h geometry are degenerate with energy Q + (3/2)J o , Q and J o being the classical VB coulomb and exchange integrals [40], respectively. The interaction between the K A and K B functions split their energies by (9/10)J o and −(3/2)J o , thus giving origin to the ground and lowest excited states. The evolution from the bond equalized to the bond alternating structure, i.e., the D 6h → D 3h process described by the displacement parameter δ, increases the energy gap between the two states. For small δ values the two levels further shift by the amount ±(9/4)J o δ 2 . Therefore, since J o < 0 the total π energy of benzene in D 3h geometry is lower in the ground and higher in the excited state than in the D 6h geometry. Similar considerations apply to cyclobutadiene when distorted along the D 4h → D 2h trajectory. As it is summarized in Table I, top, the negative resonance energies of the two molecules, J o and (9/10)J o , respectively, imply that the π-electron systems are stabilized with respect to the Kekule' structures. Nevertheless, the π energies are lower in the bond alternating than in the bond equalized configuration.
In the Hückel theory the MO energies for monocycles C 2n H 2n with conjugated double/single bonds [41] were adapted to the bond alternation case [38]  where β o < 0 is the resonance integral. In this expression the ± and j labels are combined to give rise to the 2n orbital energies. At δ = 0 the traditional expression for the orbital energies of fully symmetrical monocycles is recovered. All (4n + 2) π-electron rings have a closed shell ground state. On the other hand, 4n π-electron rings are open shell systems, due to the occurrence of two non bonding orbitals and to their occupancy with two electrons. The alternation shifts all bonding orbitals (except the lowest) downward and all antibonding orbitals (except the highest) upward so that the π energy at δ = 0 is lower than at δ = 0. To first order in δ, the shift dependence on bond alternation of the bonding and antibonding energies (except the lowest and the highest) is quadratic in δ for both 4n and (4n + 2) π systems. However, the non bonding orbitals of the 4n π systems shift linearly with δ. As a result the variation of the π energy, ∆E π , with δ is different for cyclobutadiene and benzene as reported in Table 1, bottom. Table 1. Valence bond (VB) and molecular orbital (MO) energies of cyclobutadiene and benzene [38,39]. Upper, VB theory: under K A , K B ; energies of the Kekule' structures at the most symmetrical geometries, D 4h for cyclobutadiene and D 6h for benzene; under G; E π (0) ground state π energy after K A , K B mixing at the undistorted geometries, D 4h and D 6h , respectively, ∆E π (δ) π energy change as a function of the displacement parameter δ; under Exc; the same quantities for the excited state. Q, J 0 , coulomb and exchange integrals, both negative [40]. Lower, MO theory: under E b.eq. , π energies of cyclobutadiene and benzene in the Hückel approximation; under ∆E π (δ), π energy change as a function of the displacement parameter δ; under ∆E σ (δ), σ energy change as a function of the displacement parameter δ. β 0 , β 0 resonance integral and slope at the symmetrical configuration; k , force constant of the CC σ bond.

The Role of σ Electrons
Comparing the VB and MO data on benzene it is seen that the π energies depend in both cases on δ 2 . In the harmonic approximation the quadratic dependence holds also for the σ energy [39]. The relative weight of the σ and π contributions to the total energy determines the condition under which the delocalized structure is the ground state minimum. Since in the Hückel approximation E tot (δ) = E σ (δ) + E π (δ), the following expression is obtained for benzene [39] E tot (δ) = E tot (0) where the subscript index o refers to the equilibrium geometry of the delocalized structure, k > 0 is the stretching force constant of the CC σ bond and β o = (dβ/dδ) | o is the slope of β at equilibrium. From this expression a criterion may be established: if k | β o | /(β o ) 2 −2 > 0 benzene takes on the delocalized D 6h geometry while in the opposite, k | β o | /(β o ) 2 − 2 < 0, benzene shifts to the localized D 3h geometry. Since for benzene this quantity is positive [39], the D 6h structure is the equilibrium geometry. Thus, for a given β o and β o the ultimate reason of the bond equalized geometry of benzene is the high rigidity of the σ skeleton. Applying the same treatment to cyclobutadiene, the distortion to geometries less symmetrical than D 4h follows directly within the HMO model. In fact, the total ground state energy at δ = 0 of cyclobutadiene and, more in general, of 4n π-electron systems contains a linear term in δ as seen in Table 1. This means that all these systems are prone to distortion in their fully symmetrical geometries. The distortion in (4n + 2) π-electron rings is less obvious. Increasing the annulene size n, the HOMO-LUMO energy gap becomes smaller and as a result also the energy difference between the ground and the lowest excited state. Due to second order Jahn-Teller effect [42], electronic states close to each other may interact and for n sufficiently large a structure with bond alternation may be energetically advantageous even for the higher annulenes satisfying the Hückel rule [43,44].
Following early recognitions of the propensity of π electrons to distortivity [18,45], the rather surprising statement was advanced in a series of important papers [17,[46][47][48][49][50][51][52][53]. The delocalized structure of benzene was seen as intermediate configuration between the localized K A and K B structures in the bond exchange process K A → K B along the distortion mode where CC stretchings alternate with CC shrinkings. As illustrated in Figure 1, there are two possibilities for this process. On the left, the high energy gap between ground K A , K B and excited K * B , K * A states leads to localization of double bonds. If the gap is small as in Figure 1b, the delocalized geometry is preferred. The vertical gap (V G) is related to the energy required to unpair the π electrons, which is roughly three times the singlet-triplet excitation energy of the C=C system [49]. Then, a strong bond favors localization while a weak bond favors delocalization. Since the C=C bond is relatively strong it was concluded [48] that the π electrons are distortive in benzene and that the crucial factor for symmetrization is the σ propensity to bond equalization. Figure 1. The bond exchange process between the two alternating structures of benzene: K A , K B ground state (G) structures; K * A , K * B excited state (E) structures conserving electron pairing. Ψ G and Ψ E : wave functions after K A and K B mixing. Light and heavy lines: unperturbed and perturbed energy profiles. For the sake of clarity the correlation diagrams between the Kekule' structures have been assumed to be linearly dependent on the bond alternating coordinate δ. Bond alternating (1a) and bond equalized (1b) structures are stabilized according to the height of the vertical gap (V G). QM RE: quantum mechanical resonance energy, equal in both cases.
Several ab initio strategies have been applied in order to confirm the qualitative model [17,50,51,53]. The total distortion energy ∆E πσ was partitioned into components, ∆E πσ = ∆E σ + ∆E π , displacing the C atoms in such a manner that the nuclear repulsion term between these atoms is independent of distortion. The pure σ contribution ∆E σ was determined either for the π-septuplet state of benzene (where the π delocalization is intentionally eliminated) or for the bare σ framework (C + ) 6 H 6 . Since ∆E πσ of benzene is easily calculated, the comparison shows that ∆E σ > ∆E πσ > 0 and therefore that ∆E π adds negatively to ∆E σ . The restoring tendency to symmetrization decreases when π electrons are present [17]. The same conclusion was reached also starting from different partition schemes where the nuclear repulsion term is not kept constant [54]. As to cyclobutadiene, ∆E πσ is negative and ∆E σ positive [49]. It follows that ∆E π is negative and | ∆E π |> ∆E σ . The rigidity of the σ framework is not sufficient to constrain cyclobutadiene to a square geometry and the molecule shifts toward a rectangular geometry in the ground state. The quantities ∆E πσ , ∆E σ and ∆E π for the two molecules are collected in Table 2, shown later on for sake of clarity. The response to distortion is qualitatively illustrated in Figure 2. For benzene, the negative ∆E π curvature at the symmetrical equilibrium configuration is overcompensated by the positive curvature of ∆E σ ; the opposite occurs for cyclobutadiene. Figure 2. Schematic representation of the σ and π response to distortion in benzene (left) and cyclobutadiene (right). ∆E σ and ∆E π are the energy changes as a result of the bond alternating distortion and ∆E πσ = ∆E σ + ∆E π .

Quantum Mechanical and Thermochemical Resonance Energy
The stability of the (4n + 2)π-electron conjugated monocycles and instability of the 4nπ counterparts is commonly expressed by the Hückel rule [4,10]. A strong support to this rule is given by the chemical stability of benzene and by the elusive nature of cyclobutadiene. The Hückel rule has therefore become one of the most powerful paradigms of chemistry: aromatic species are associated with (i) delocalized arrangements of the π electrons, (ii) bond equalized geometries and (iii) special stability relative to an open-chain polyenic reference with the same number of C atoms [55][56][57]. In the converse, antiaromatic species are associated with bond localization/alternation and high reactivity [21]. A correlation has been assumed between π delocalization and geometry and this has led to the concept that geometries with equal bond lengths, such as in benzene, are due to the inherent tendency of the π electrons to delocalization. As this view is opposed by theoretical evidences about the π propensity to distortion in benzene and in cyclobutadiene [17,[48][49][50][51]53], arguments were advanced to reconcile the new structural approach with the Hückel rule and with well-known aromatic properties including thermochemical stability and substitution reactivity. This essential point has been discussed in detail [17] and will be shortly reviewed in the following paragraphs. To this purpose let us first recall the types of resonance energies pertinent to the discussion.
The quantum mechanical resonance energy (QM RE) is the energy lowering due to the linear mixing of the K A and K B wave functions with respect to the energy of the Kekule' structures at the same nuclear geometry. As made evident in Figure 3, no π conjugation among C=C bonds is associated with the Kekule' structures. The QM RE's of cyclobutadiene and benzene are largely different, being for instance QM RE(C 4 H 4 ) ≈ (1/3) QM RE(C 6 H 6 ) in a VBSCF/6-31G calculation [17]. The difference reflects the essence of the Hückel rule and indicates that π delocalization in benzene is favored over π delocalization in cyclobutadiene. A second resonance energy is adiabatic (ARE) where the Kekule' structures at the bond-equalized geometries are substituted by the optimized bond-alternating geometries in absence of all π conjugation. Taking as example benzene in Figure 4, ARE and QM RE are related according to the equation [58] QM RE = ARE + E c where E c is the compression energy from the non-conjugated optimized to the non-conjugated bond-equalized geometry.
Finally, the thermochemical, Dewar-type, resonance energy (T RE) is defined as the energy difference between the conjugated monocycle and an appropriate reference system R [59][60][61][62][63]. In the present case R is the hypothetical cyclic hexatriene or cyclic butadiene with the same amount of π-conjugation as in the open-chain polyene, as sketched in Figure 5. T RE is referred to also as aromatic stabilization energy (ASE) for benzene [58] while it is a sum of two contributions for cyclobutadiene, the (negative) ASE term due only to the antiaromatic destabilization and the angular strain term [64]. Since the resonance energies of open polyenes depend linearly on the number n of conjugative interactions between formal double bonds [65], ASE may be evaluated by subtracting the contribution of the hypothetical reference molecule from ARE [4,58] ASE = AREn x ARE (butadiene) where n = 2, 3 for cyclic butadiene and hexatriene, respectively.  Empirical estimates of resonance energies are obtained by means of isodesmic and homodesmotic reactions [66,67] or heats of formation [68,69]. The isodesmic reaction C 6 H 6 + 6CH 4 → 3(H 2 C = CH 2 ) + 3(H 3 C − CH 3 ), having equal numbers of formal single and double CC bonds in reactants and products, gives an approximate determination of ARE (C 6 H 6 ) since it compares the energy of benzene with that of three ethylene molecules. In a homodesmotic reaction the previous condition is strengthened, further demanding for reactants and products equal numbers of (i) C atoms in the same hybridization state and (ii) sp 2 and sp 3 C atoms with one, two or three hydrogen atoms attached. Examples of homodesmotic reactions for benzene and cyclobutadienes are reported in Figure 6 where the cyclic standards are 3 and 2 molecules, respectively, of s-transor s-cis-butadiene.
Due to the proper number of conjugative interaction in the acyclic reference, the enthalpies of these reactions are easily related to ASE. In the case of benzene, ASE is 21.6 kcal/mol and 28.8 kcal/mol with respect to s-transand s-cis-butadiene, respectively [58,68]. ASE (cyclobutadiene) has been recently investigated in detail and determined to be -42 kcal/mol and -37 kcal/mol with respect to s-transand s-cis-butadiene [64]. It has been further noted that the intrinsic ASE value of cyclobutadiene, i.e., resulting from separating the π repulsion effect from the total ASE with respect to the s-cis-butadiene is -12 kcal/mol [64]. Going back to the point in discussion, T RE may be derived [17] coupling the π-distortive model of benzene and cyclobutadiene with a thermochemical cycle involving benzene/cyclobutadiene (B; CB), distorted benzene / distorted cyclobutadiene (DB; DCB), Kekule' structures of distorted benzene /cyclobutadiene (DKB; DKCB) and the reference molecule R, i.e., cyclic hexatriene/cyclic butadiene.
In the analogous cycle for cyclobutadiene, Figure 7 lower, DKCB and R ( see Figure 5) are substantially different structures for in one case, DKCB, (i) the π bonds face one another and (ii) the ring is strained while in the second, R, the first effect is greatly reduced and the second is absent. In the third step not only partial π restoration occurs with resonance energy 17.9 kcal/mol [48] but in addition a relaxation term must be considered due to the decrease of the two effects. All recently reported values of angular strain are between 32 and 35 kcal/mol [64,[70][71][72]. ∆E rel is estimated to be largely negative, −59 kcal/mol from the strain energy of 34 kcal/mol [64] and the repulsion energy between π bonds of 25 kcal/mol c. Then, ∆E 3 = −QM RE(R) + ∆E rel = −76.9 kcal/mol. The other two terms are ∆E 1 = ∆E πσ = −3.4 kcal/mol (from Table 2) and ∆E 2 = QM RE(DCB) = 8.3 kcal/mol [17]. It results T RE(CB) = ∆E 1 + ∆E 2 + ∆E 3 = −72.0 kcal/mol, in good agreement with the reported value of T RE [17]. Further, if the strain and π repulsion energies are left aside, the portion of T RE(CB) associated with delocalization is only -13 kcal/mol, closely matching the intrinsic ASE(CB) value [64]. The thermochemical instability (with respect to the reference molecule) of CB is related to the instability of the Kekule' structure and not to absence of QM RE at the square geometry. The distortive π propensity of benzene and cyclobutadiene is consistent with the stability/instability of the hexagonal/square geometry with respect to the reference molecules [49]. ; R reference molecule, cyclic hexatriene and cyclic butadiene. T RE and QM RE: thermal and quantum mechanical resonance energies. ∆E rel : relaxation energy due to angle strain and π overlap in cyclobutadiene; this term is not considered for benzene [17]. The energy values (kcal/mol) characterizing the thermochemical cycle of benzene and cyclobutadiene are indicated.
Because of the σ-imposed rigidity of the C framework, benzene has geometries restricted to those close to perfect hexagons even at non-equilibrium configurations. Since the resonance energy is only moderately affected by small displacements from D 6h symmetry, the aromatic stability is conserved in all experiments probing near equilibrium properties such as reactivity, UV spectroscopy and magnetism. Much less common are experiments on π-distortivity. The unusually low infrared frequency of the b 2u bond alternating mode was an early evidence of π-distortivity related to the fact that the b 2u vibrating molecule spends a non vanishing fraction of time distorted as cyclohexatriene [45]. More recently, two-photon spectroscopy on the lowest excited B 2u state has supported this conclusion [73][74][75]. As the frequency of the b 2u mode, which is strongly active as vibronically induced two-photon B 2u × b 2u transition, is considerably increased with respect to the ground state value, this result fits nicely the model of Figure 1 since the π curvature adds to the σ curvature in the B 2u state while subtracts in the ground state. This behaviour is documented also for other aromatic hydrocarbons [76][77][78][79][80][81][82][83].

Cyclobutadiene in Square Geometry
In the MO theory the four π orbitals of square D 4h cyclobutadiene are classified into the a 2u , e g and b 1u symmetry species in order of increasing nodal planes [84,85]. All the MO's are determined by symmetry with the two degenerate e g orbitals characterized by one nodal plane and distributed either over the whole ring, φ 2 and φ 3 , or, equivalently, localized on opposite ring atoms, φ 2 and φ 3 . The two alternatives are shown in the schematic representation of Figure 8. The four π electrons fill the lowest a 2u orbital and partly the e g pair. With two orbitals containing altogether two electrons, six electronic configurations are possible, Since the a 2u orbital is doubly occupied, the a 2 2u × e 2 g representation is reduced to 1 A 1g + 3 A 2g + 1 B 1g + 1 B 2g . The four states, expressed in terms of only the φ 2 and φ 3 orbitals, i.e., omitting for simplicity the inner σ electrons and the lowest a 2u orbital, and the corresponding HF energies are [85] is the M s = 0 component of the triplet state. Since coulomb and exchange two-electron repulsion integrals are always positive, is the lowest state and Φ( 1 A 1g ) the highest among the four states. It may be also verified [85] that E( 1 B 1g ) < E( 1 B 2g ) when φ 2 and φ 3 are expressed in terms of the localized orbitals φ 2 and φ 3 , i.e., φ 2, In the VB formalism the two lowest singlet states are described by mixing the two Kekule' configurations K A and K B . Given the a, b, c, and d atomic orbitals, K A and K B are expressed as where A is the antisymmetry operator. Their variational mixing leads to the 1 B 1g state as the minus and to the 1 A 1g state as the plus combination Four equivalent configurations with one π bond broken (see Figure 8) contribute to the two lowest triplet states. For instance, the T 1 configuration is defined as and similarly for the other three triplet configurations T 2 , T 3 , T 4 . Two states, 3 A 2g and 3 E u , result from the linear combination of the four configurations since only three linearly independent combinations are possible for a system of four electrons in four orbitals [40]. Overall, the singlet pair, i.e., one resonant 1 B 1g and one antiresonant 1 A 1g state, and the triplet pair, i.e., one resonant 3 A 2g and one antiresonant 3 E u state, are obtained. Considering that E(T 1 ) > E(φ ab,cd ) and reasonably assuming that the resonance energy of the 3 A 2g triplet is not much greater than that of the 1 B 1g singlet, the ground state is predicted to be singlet and the first excited state to be triplet [85]. Qualitatively, both theories account for low-lying 1 B 1g and 3 A 2g states, however in reverse order. The relative energies of Table 3 obtained from ab initio SCF calculations [85] are in agreement with these considerations and indicate that at the HF/VB calculation level the lowest state of square cyclobutadiene is 1 3 Let us now consider the effect of configuration interaction on the MO energies of the singlet 1 1 B 1g and triplet 1 3 A 2g states. Singly-excited configurations within the ππ * space may arise from three different promotions, a 2 2u → a 1 2u b 1 1u , a 2 2u e 2 g → a 1 2u e 3 g and e 2 g → e 1 g b 1 1u . The only excited states interacting with the Out of the total 24 a 1 2u × e 2 g × b 1 1u configurations, 16 have a single occupancy of each π orbital and 8 a double occupancy of one of the two e g orbitals (and zero of the other) [92]. The 16 configurations are combined to six states, 5 B 1g , 2 3 B 1g , 3 A 2g , 1 B 1g , 1 A 2g while the other 8 configurations to four states, 2 1 A 1g and 2 3 B 2g . The lowest triplet, 1 3 A 2g , interacts with one singly excited state 3 A 2g and the lowest singlet, 1 1 B 1g , with one singly excited 1 B 1g state. Since the corresponding matrix elements differ approximately by a √ 3 factor in favor of the lowest singlet [93], the stabilization energy of 1 1 B 1g is three times that of 1 3 A 2g and may be larger than the unperturbed | E(1 1 B 1g ) − E(1 3 A 2g ) | energy gap. In principle, the singlet 1 1 B 1g may become the ground state for square cyclobutadiene, thus settling the difference with the VB results. As to the other pair of states, 1 1 B 2g and 1 1 A 1g , no 1 B 2g within the ππ * space is available for interaction while 1 1 A 1g may interact with the pair of singly-excited 1 A 1g states. The 1 A 1g energy may decrease until inversion with respect to SCF results. In fact, there is a large number of calculations made with increasingly refined methods, from full π-CI to MC/CASSCF to MRCC, which corroborate the two predictions [85][86][87][88][89][90][91]93,94]. It may be seen from Table 3 that 1 1 B 1g lies 6 -10 kcal/mol below 1 3 A 2g and that 1 1 A 1g is energetically lower than 1 1 B 2g . Table 3. Energies (kcal/mol) of the lowest excited states of square cyclobutadiene relative to that of 1 1 B 1g . (not considered in [85]). (c) limited π-CIS; 3-21G basis set [86]. (d) full π-CI; 6-31G * basis set [87]. (e) full π-CI; Dunning double-ζ basis set [85]. (f ) Dunning double-ζ basis set [88]. (g) CAS(4,4); 6-31G basis set [89]. (h) Multireference coupled cluster calculation with single and double excitations augmented non-iteratively by triple excitations; DZP basis set [90]. (i) Multireference average-quadratic coupled cluster calculation; aug -cc-pVTZ basis set [91].
The stabilization of 1 1 B 1g state with respect to 1 3 A 2g in the configuration interaction approach is also known as dynamic spin polarization of cyclobutadiene [93]. Shortly, the spin polarization of the two paired a 2 2u electrons depends on the spin state of the two unpaired e g electrons. Resorting to the φ 2 and φ 3 description of the e g orbitals it may be seen [93] that when two antiparallel electrons occupy φ 2 and φ 3 the a 2u electrons are easily polarized by mixing a 2u with the antibonding b 1u MO. On the contrary, with parallel electrons in φ 2 and φ 3 virtual orbitals of much higher energy are required for mixing, giving a smaller polarization contribution to the lowering of the 1 3 A 2g energy.

Structural Instabilities
Molecular geometries of non-degenerate electronic states may be unstable due to coupling with higher excited states through non-totally symmetric distortions [42,43]. The effect has been called second-order (or vibronic) Jahn-Teller effect and is distinguished from the first-order effect, relative to degenerate states [95]. For instance, the 1 1 B 1g state of square cyclobutadiene couples with 1 1 A 1g when the molecule is distorted along the rectangular b 1g mode consisting of adjacent in-and out-of-phase CC bond stretchings around the ring. For a sufficiently strong coupling term the square (D 4h ) → rectangular (D 2h ) structural change occurs. The consequent state correlation is This process may be viewed also as a bond shifting reaction interconverting two equivalent D 2h structures (see Figure 9) through the transition structure of D 4h symmetry [87]. The barrier height has been measured to be in the range 1.6-10 kcal/mol [96,97]. The stabilization energy ∆E(1 1 B 1g − 1 1 A g ) going from square to rectangular geometry has been calcultated [85][86][87][88][89][90][91]93,94,98] between 5 and 13 kcal/mol according to the results of Table 4.
Recent estimates including ZP E correction [89][90][91] lower the height to 3-6 kcal/mol. In these calculations accurate ground state structures of cyclobutadiene have been also reported exhibiting geometrical parameters "very close to the true equilibrium bond distances" [91]. Ab initio normal mode calculations on the 1 1 B 1g (D 4h ) state confirm its saddle point character given that imaginary frequencies, 990i, 803i, 1579i cm −1 , are found for the bond alternating b 1g vibration [89][90][91]. On the other hand, the 1 1 A g (D 2h ) state is associated with a minimum structure since all the vibrational modes have real frequencies. Geometries and frequencies [89][90][91]99] are reported in Table 5 while the calculated [90] potential energy profile is shown in Figure 10, left. Table 4. ∆E(1 1 B 1g −1 1 A g ) energy barrier (kcal/mol) for the interconversion process square D 4h → rectangular D 2h (first two rows) and for the process square D 4h → rhombic D 2h process (last row). Values in parentheses include the Zero Point Energy (ZP E) correction.
The distortions of cyclobutadiene (other than rectangular) to geometries lower than square may be guessed considering that a square X 4 ring has, in addition to the b 1g mode, three non-totally symmetric normal displacements, 1b 2g + 1b 1u + 1e u . The b 1u out-of-plane mode cannot couple the lowest excited ππ * states of cyclobutadiene which are all symmetric with respect to the molecular plane. The b 2g CCC bending mode couples 1 1 A 1g with 1 1 B 2g and is responsible of the square (D 4h ) → rhomboidal (D 2h ) distortion. The consequent state correlation is 1 1 A 1g → 1 1 A g ; 1 1 B 2g → 2 1 A g . The 1 1 A g state of rhombus cyclobutadiene is stabilized with respect to 1 1 A 1g of square cyclobutadiene by ≈ 3 kcal/mol. The rhomboidal stabilization (see Figure 10, right) is triggered by the b 2g bending mode with imaginary frequency, 9423i cm −1 (overestimated [89]) and 1104i cm −1 [90]. As to the 1 1 A 1g state (square geometry) it should be added [89] that in correspondence of the ground state stabilization along the rectangular distortion, this state is destabilized along the same coordinate. In other words, the b 1g frequency is considerably strengthened in 1 1 A 1g (square geometry), up to 2125/1783 cm −1 [89,90]. Second, the rhomboidal structure of the 1 1 A g state is not a true minimum structure. Two modes of b 2g and b 1u symmetry retain imaginary frequencies (see Table 5) and suggest the occurrence of a minimum structure with symmetry lower than rhomboidal [89,91]. Closely related potential energy curves have been documented for the cyclopentadienyl cation along the D 5h → C 2v distortion mode [100].
symmetry in order of increasing energy) as a function of the rectangular and rhomboidal distortion. Adapted MRCCSD(T) results [90]. The state correlation for the Finally, the e u mode, where one CC bond lengthens and the opposite shortens, may cause a square However, with the E u states high in energy, the coupling is expected to be weak and a trapezoidal equilibrium geometry is not reached in any of these states [87]. By the same argument applied to 1 3 A 2g it follows that, due to the large energy gap with respect to the next nearest state, 1 3 E u , the lowest triplet is stable toward trapezoidal and a fortiori all types of planar deformations [94]. The minimum D 4h structure of the 1 3 A 2g state in Table 5 is confirmed by the absence of imaginary frequencies in all ab initio vibrational calculations [89][90][91]. The square triplet structure of cyclobutadiene agrees with the reversal of the (4n + 2)/4n rule [22,101] according to which 4n π-electron rings with bond alternating geometry in the lowest singlet have bond equalized structures in the lowest triplet. Triplet aromaticity of cyclobutadiene and more in general of 4n π-electron annulenes up to the cyclononatetraenyl cation has been quantitatively estimated by high-level ab initio calculations in terms of aromatic stabilization energies and of geometric and magnetic properties [102]. Table 6. Bondlengths (Å) and bond angles (degrees) of benzene in the D 6h , D 2h (quinoidal and antiquinoidal) and C s structures. Structural data indicated but not reported are determined by symmetry. Imaginary frequencies ν (cm −1 ) of the 1 3 B 1u (D 6h ) and 1 3 B 1u (D 2h ) triplet states. ∆E (kcal/mol) is the energy barrier between minima.

Benzene: Ground and Lowest Excited States
calc (e−l)   [103]. (b) neutron diffraction data [104]. (c) electron diffraction results [105]. (d) rotational Raman results [106]. (e) MP2/6-311++G * * results [113]. (f ) MCSCF/CI/6-31G results [115]. (g) calculated values, Feynman path integral method [116]. (h) HF/MO ab initio results [111]. (i) HF/6-31G results [112]. (l) MP2/TZ2P+f calculation [114]. (m) MCSCF/6-31G results [119]. (n) MRDCI results [120]. (o) experimental determination [121]. (p) experimental determination [122]. (q) MCSCF/CAS(6,6)/STO-3G results [123]. (r) MCSCF/CIS/6-31-3G results [115]. (s) SCF/MO results [124]. (t) SINDO/CI calculation [125]. (u) MRDCI results [120]. (v) SCF/CI/4-31G calculation. Two angles are in addition defined for the C s structure: α = 44 • , β = 67 • , α and β being the angle between planes C 1 C 2 C 6 and C 3 C 4 C 5 with the plane All the structural determinations from X-ray [103], neutron [104] and electron [105] diffraction studies to rotational Raman spectroscopy [106] converge on a regular hexagonal D 6h geometry for the structure of benzene in the ground state. However, it has been acutely noted [107] that the X-ray structural data are compatible not only with an ordered crystal of D 6h molecules but also with one disordered with D 3h molecules overimposed and rotated 60 • . With a 0.10Å difference between C-C and C=C bondlengths the crystal disorder would be far below the resolving power of the X-ray diffraction technique. Since all other diffraction techniques are equally unable to resolve the D 6h /D 3h dilemma [107], this means that there are strong experimental indications in favor of the D 6h geometry which nevertheless do not allow to discard unequivocally the D 3h alternative. In spite of this, support to the D 6h geometry has been more recently provided by the analysis of the temperature dependence of the neutron diffraction data relative to the C 6 D 6 crystal [108,109]. Based on the atomic displacement parameters of the C atoms at 15 and 123 K [109], the model shows that the small in-plane temperature independent contribution to these parameters is due to the zero-point intramolecular motion and not to the disorder between Kekule' structures.
The theoretical calculations take on great significance [110]. Ab initio calculations [111][112][113][114][115][116][117] point to the D 6h structure as that of minimum energy for benzene on the ground state energy surface. The data [111][112][113][114][115][116] are collected in Table 6 and compared with the most important structural results. However, anomalous results such as the occurrence of imaginary frequencies in D 6h symmetry have been reported when incomplete basis sets are used in conjunction with correlated calculation methods [118].
In order of increasing number of nodal planes the six π orbitals of D 6h benzene are classified into the a 2u , e 1g , e 2u and b 2g symmetry species. The ππ * HOMO → LUMO excitation, i.e., e 4 1g → e 3 1g e 1 2u , gives origin to six ππ * states, three singlets, 1 B 2u , 1 B 1u and 1 E 1u , and three triplets, 3 B 2u , 3 B 1u and 3 E 1u . In addition, both the a 2 2u → a 1 2u e 1 2u and the e 4 1g → e 3 1g b 1 2g promotions contribute to the 1 E 2g and 3 E 2g ππ * states. The energy level diagram established on the basis of the assignment of the observed excitation energies [35,127] is summarized in Table 7. All the excited state calculations [116,[128][129][130][131][132] are consistent with the energy ordering 1 From the point of view of the VB theory these states may be classified into two groups, covalent and ionic [128,129,133]. States with covalent character are 1 B 2u , 3 B 1u and 3 E 1u described by Kekule', Dewar and quinoidal structures ac of Figure 11. Structures with charge separations [for instance, (d) in Figure 11] predominantly contribute to 1 B 1u , 1 E 1u and 3 B 2u , thus imparting ionic character to these states [133]. The excitation energies to covalent states are fairly well reproduced even in less extended ab initio calculations while more sophisticated calculation methods are needed for the calculation of the excitation energies to ionic states [129].
The equilibrium geometries of the lowest ππ * states have been the object of several studies [119,120,126,134]. For the lowest triplet, 3 B 1u , a structure with symmetry lower than D 6h is expected on the basis of the reversal of the (4n + 2)/4n rule [101] . Due to the coupling of 3 B 1u with the close lying 3 E 1u state, benzene shifts from hexagonal to structures of quinoidal Q and antiquinoidal AQ character sketched in Figure 11 (four long, two short bonds for Q; four short, two long bonds, for AQ). Both have D 2h symmetry.
The calculation indicates that 3 B 1u (D 6h ) is the transition state between AQ and Q minima whose energies are ≈ 2.3 kcal/mol below the transition state [120,134]. The saddle point nature of the 3 B 1u (D 6h ) state is confirmed by the occurrence of one imaginary ν 8 (e 2g ) frequency 924i cm −1 in the normal mode calculation [119]. The structures of 3 B 1u (D 6h ) and of the AQ and Q states are reported in Table 6. Evidence of 3 B 1u / 3 E 1u coupling comes from low-temperature EPR experiments according to which the lowest triplet of benzene is non hexagonal [135,136]. Further, the phosphorescence excitation spectrum [137] shows an intense doublet 239/252 cm −1 above (0-0) assigned [138] as the lowest vibrational transition of the ν 8 (e 2g ) vibration. This vibration has frequency 1584 cm −1 in the ground state. The strongly anharmonic potential well around the Q and AQ minima sets a value of 220 cm −1 for the vibrational transition 0 → 1 of the ν 8 (e 2g ) mode [134]. Table 7. Lowest singlet and triplet excited states of benzene: experimental S 0 → S n and S 0 → T n band maxima (eV) and calculated transition energies (eV).  According to experiment [35,121,122] and calculation data [115,123,124] the 1 1 B 2u state retains hexagonal symmetry with uniform expansion of the carbon ring and CC bondlength increase of ≈ 0.04Å with respect to the ground state. As to the second excited singlet, 1 1 B 1u , a SINDO calculation predicts a mere increase of CC bond lengths [125]. In contrast, ab initio calculations indicate [120] that 1 1 B 1u , similarly to 1 3 B 1u , distorts toward a quinoidal/antiquinoidal structure. The behaviour was justified noting the similarity of energy gaps between the singlet 1 1 B 1u /1 1 E 1u and triplet 1 3 B 1u /1 3 E 1u pairs and the independence of the vibronic coupling on the spin state. A biradicaloid structure of C s symmetry has been also reported for 1 1 B 1u on the basis of ab initio, limited CI, calculations, as intermediate structure in the photochemical conversion benzene → Dewar benzene [126].

Cyclooctatetraene: A Non Planar Cyclic Conjugated System
Due to angular strain cyclooctatetraene, the next higher homologue in the annulene series, is a non aromatic molecule with a non planar structural arrangement of four weakly interacting π bonds [4,10,11]. It is known from electron [139][140][141] and X-ray diffraction [142,143] data and recent results of femtosecond spectroscopy [144] that the molecule adopts a tub-shaped structure of D 2d symmetry, as shown in Figure 12.
Two basic dynamical processes characterize the molecule [145] (see Figure 13), (i) ring inversion  The four equienergetic D 2d conformations a, b, c, d are interconverted passing through transition states of D 4h symmetry for RI and D 8h symmetry for BS. In cyclooctatetraene only BS is probed by N M R spectroscopy [6]. The transition state is ≈ 14 kcal/mol above the tub structure [6,146]. In monoalkyl derivatives both processes are observed and the energy barriers are ≈ 11 kcal/mol for RI and 15 kcal/mol for BS [6]. The larger energy involved in the BS process is due to the fact that ring flattening is accompanied by bond length changes in order to have a D 8h transition state. The RI/BS barriers have been calculated for cyclooctatetraene with ab initio methods, ≈ 10/16.4 kcal/mol and 17.1/19.3 kcal/mol, respectively [147,148]. The mode responsible for ring inversion is the b 2u out-of-plane mode with frequency 139i cm −1 in D 4h symmetry while that for bond shift is the bond alternating b 1g mode with frequency 1962i/2363i cm −1 in D 8h symmetry [147,148].
The structures and energies (relative to the tub minimum) of the D 4h and D 8h stationary points and of one selected isomer and two conformers (bicyclo[4.2.0]octa-2,4,7-triene, crown and chair; C s , D 4 and C 2h symmetry) are reported in Figure 14. The chair and crown conformers are stable species 56.5 and 171.1 kcal/mol above the tub conformation [148]. They are obtained from the tub by rotation of two and four CC double bonds so that the π orbitals are directed toward the internal part of the cyclic structure, thus increasing the repulsion among π electrons [148]. The tub and bicyclo isomers are only a small fraction of the total number of (CH) 8 valence isomers, 21 including stereoisomers [149,150]. Disregarding H atoms, they correspond to regular graphs of degree three [150], i.e., their formulas are drawn under the condition that (i) three lines (bonds) always meet at each vertex (C atom) and (ii) no two of these lines cross. If these conditions are relaxed additional (CH) 8 isomers can be derived. On the left, D 2d cyclooctatetraene and BS, RI transition states of D 8h and D 4h symmetry; on the right cyclooctatetraene isomers (crown, chair and bicyclo[4.2.0]octa-2,4,7-triene, D 4 ,C 2h and C s symmetry, respectively). Structural parameters (Å) and energies (kcal/mol) relative to D 2d minimum, reported results [148] (MP2-CASSCF//CASSCF/6-31G* calculations). Experimental energy values (kcal/mol) of the D 8h and D 4h transition states are given in parentheses [151].
The properties of these isomers have been reviewed [36] and the isomerization reactions, such as those shown in Figure 15, have been discussed by means of extensive computational studies [152][153][154][155]. In particular, the relaxation pathways relevant to the photochemical reactivity of cycloctatetraene after excitation into S 2 [1E(D 2d )] have been described in terms of radiationless decay from S 2 into the S 1 -excited planar minimum of D 8h symmetry and internal conversion to S 0 through S 1 /S 0 conical intersections [156]. As to the transition structures of D 4h and D 8h symmetry, their occurrence has been provided by photodetachment experiments on the radical cyclooctatetraene anion [151]. Two states have been observed, 1 1 A 1g (D 4h ), i.e., the ground state of the D 4h transition structure and 1 3 A 2u (D 8h ), the lowest triplet state of the D 8h transition structure. Their energies with respect to the D 2d minimum have been determined, 23.5 kcal/mol. Since the difference between the RI and BS energy barriers is ≈ 4 kcal/mol [6], this is experimental evidence [151] that the lowest D 8h singlet, i.e., 1 1 B 1g (D 8h ), lies ≈ 8 kcal/mol below the lowest D 8h triplet, Pentalene Pentalene ( Figure 16) is formally derived from the regular octagonal π perimeter of cyclooctatetraene by cross-linking opposite atomic centers [157]. Due to the reduction of molecular symmetry from D 8h to D 2h , all the π orbitals of pentalene are non degenerate and the lowest four are occupied by eight π electrons. Simple arguments based on the second-order bond fixation in conjugated molecules [43,44] indicate that the symmetrical bond equalized D 2h structure corresponds to a maximum on the ground state energy surface and that the molecule distorts to C 2h symmetry with alternating single and double bonds. Bond alternation of the order of ≈ 0.15Å has been calculated by all ab initio methods for ground state C 2h pentalene [28,29,158] (see Table 8).
The highly reactive molecule is generated by photochemical cleavage of pentalene dimers with 254-nm light at 20 K [28]. Infrared and electronic spectra of pentalene have been measured and interpreted on the basis of ab initio calculations on C 2h molecular symmetry [28]. The three absorption bands at 3.48, 3.94 and 4.80 eV are assigned to the lowest allowed S 0 → nB u (n = 1, 2, 3) transitions.
At lower energy two states of A g symmetry are predicted with spectral activity too weak to be observed in the unsubstituted system. Their occurrence is due to the fact that in pentalene the HOM O → LU M O singly and doubly excited configurations are of a g symmetry and comparable energy [28]. Table 8. Upper: bond lengths (Å) of pentalene (D 2h and C 2h equilibrium structures, numbering as in Figure16) and of 1,3,5-tri-tert-butylpentalene (C s equilibrium structure) in the ground state. Structural data indicated but not reported are determined by symmetry. Imaginary frequency ν (cm −1 ) of the 1 1 B 1g (D 2h ) state. ∆E (kcal/mol) is the energy barrier between minima of C 2h symmetry. Lower: bond lengths (Å) of pentalene in the lowest excited states, from S 1 to S 4 . Imaginary frequency ν (cm −1 ) of the 2 1 A g (D 2h ) state. ∆E (kcal/mol) is the energy barrier between minima of C 2v symmetry.   [29]. (b) B3LYP/6-31G * results [28]. (c) B3LYP/6-31G * results [158]. (d) X-ray structural data of dimethyl-4,6-di-tert-butyl-1,2-dicarboxylate [159] (pair of values relative to symmetrically equivalent distances assuming a C 2h structure). (e) X-ray data of 1,3,5-tri-tert-butylpentalene [159] (pair of values relative to symmetrically equivalent distances assuming a C 2h structure). The first stable derivative at room temperature, 1,3,5-tri-tert-butylpentalene, was synthesized hindering sterically the dimerization [27]. The ring geometry was determined by X-ray diffraction [159]. Low energy absorption bands have been observed at 1.50 and 2.14 eV in the low temperature spectrum of the tert-butyl derivative and assigned to the two A g states previously mentioned [161]. Further, the derivative fluoresces from the lowest B u state at 3.17 eV (≈ 25550 cm −1 ) with quantum yield 0.002 and decay time 13 ps [29,160]. This band has been assigned as the emission counterpart of the 3.67 eV (≈ 29600 cm −1 ) absorption band (see Figure 17). Other two absorption bands are observed with maxima at 4.52 (≈ 36500 cm −1 ) and 5.76 eV (not shown in Figure 17). Calculated transition energies of the tri-tert-butyl derivative are in good agreement with experiment [160].  Similarly to cyclobutadiene and cyclooctatetraene, pentalene interconverts between C 2h isomers. In 1,3,5-tri-tert-butylpentalene the energy barrier was found to be ≈ 4 kcal/mol [162]. The calculated barrier height amounts to 9.7 kcal/mol [29] and 6.6 kcal/mol (including ZP E) [158]. Another structural rearrangement from D 2h to C 2v symmetry occurs in the second excited state due to the 2A g /1B 2u , i.e., S 2 /S 3 , vibronic coupling promoted by a b 2u vibration [29]. Overall, the energy diagram of the ground and of the lowest three excited states of pentalene is shown in Figure 18 as a function of the deviation from D 2h geometry. The excited state equilibrium geometries are reported in Table 8. Figure 18. Energy profiles of the lowest four states of pentalene ( 1 B 1g , 1 1 A g , 2 1 A g , 1 B 2u in D 2h symmetry) as a function of the displacement from D 2h geometry (vertical dotted line). Adapted MCSCF/CAS(8,8)/3-21G results [29]. For the sake of clarity the S 0 and S 2 structures of pentalene displaced from D 2h symmetry have strongly exaggerated differences. Actually, the S 0 minimum has C 2h symmetry (C 2 normal to the molecular plane) while the S 2 minimum has C 2v symmetry (C 2 along the short in-plane axis).

Higher Annulenes
In this Section the structural properties of fully conjugated monocyclic hydrocarbons, from [10]-to [18]annulenes, will be reviewed together with those of their most significant derivatives. Historical landmarks in the chemistry of these molecules are the synthesis of monocycles up to [30]annulene [163,164] and of the first bridged [10]annulene [165,166]. Due to angle strain and/or steric interactions between hydrogen atoms higher annulenes are characterized by a rich variety of non planar conformational and configurational isomers and of interconversion processes between them [6]. As a consequence, cis/trans isomerization takes place among configurations and several isomers usually exist in solution. A second important factor on which the annulene structures depend is the decreasing energy difference between localized and delocalized geometries as the ring size increases [18]. For annulenes containing more than 30 C atoms the localized arrangement of double bonds has been predicted to become energetically favored [19]. Limiting the attention to medium size annulenes, the insertion of atoms (or groups of atoms) as bridges between opposite C atoms or the introduction of triple bonds along the cycle in place of the alkene unit stiffens the ring framework. Bridged and dehydroannulenes are then formed with planar (or approximately planar) ring geometries and considerably reduced bond alternation.

Conformational Isomers
[10]-annulene was first trapped in a reaction mixture [167] and later isolated as two distinct crystalline isomers [168,169]. The first, (A), showed a single line in the 13 C and 1 H N M R spectra even down to −160 • C; on cooling the second, (B), below −100 • C a five line 13 C N M R spectrum was observed.
Based on these experiments and chemical evidence, A was assigned the all-cis and B the mono-trans configuration, in both cases in non planar conformation [169]. From the theoretical point of view the determination and the energy ordering of the C 10 H 10 isomers have represented a difficult task for the most advanced ab initio treatments [112,[170][171][172][173][174][175][176][177]. Five conformational structures have been identified as the lowest energy forms of [10]annulene, three of them non planar and two quasi planar, schematically shown in Figure 19. The three non planar isomers are all-cis or boat, mono-trans or twist, di-trans or naphthalene-like and the two quasi planar are azulene-like and heart, with the indication of the symmetry group. Figure 19. The low energy isomers of [10]annulene. The symmetry group, usual names and bondlengths (Å, CCSD/DZd and MP2/DZd results [174]) are indicated. Boat, twist and naphthalene-like isomers are non planar, azulene-like and heart are quasi planar.
At the HF calculation level these structures are local minima [170,171] except the heart isomer which originally [169] was proposed as the transition state between equivalent twist structures. The HF results of Table 9 are in good agreement with the assignment, being the twist and boat isomers the lowest energy conformers. On the other hand, correlated calculation methods such as density functional and MP2 assign the lowest minimum to the heart isomer [172]. This approach, however, has been shown to overestimate [173] the stability of the heart with respect to the twist isomer: In fact, the calculation at the correlated geometries using the high-level CCSD(T) method shifts the twist structure below heart by 3-7 kcal/mol. The result was subsequently confirmed starting with twist, naphthalene-like and heart structures optimized at the CCSD(T)/DZd level [174]. The final reference values, i.e., CCSD(T) single-point energies with the triple-ζ polarized TZ2P basis set at these geometries are 0, 1.4 and 4.24 kcal/mol, respectively (see Table 9). The 13 C N M R spectrum calculated with this method supports the assignment of the N M R spectrum of B to the twist form while the calculation does not help in the assignment of the A spectrum [175]. Replacing one CH pair with one nitrogen atom in aza [10]annulenes and thus lowering the repulsive interaction of the inner hydrogen atom with the ring in the heart isomer (refer to Figure 19 for clarity), heart becomes energetically competitive with twist. Nevertheless, this is not enough: the CCSD(T) calculation indicates that twist aza [10]annulene remains more stable than heart aza [10]annulene by ≈ 2.1 kcal/mol [178]. Table 9. Energy [174] (kcal/mol) of the [10]annulene isomers of Figure 19.  Among [10]annulene isomers the all-trans isomer, i.e., all-trans-cyclodecapentaene not shown in Figure 19 but similar to D 4 (crown) cyclooctatetraene (or [8]trannulene according to this classification) of Figure 14, has the remarkable property of p orbital overlap parallel to the average ring plane ("in-plane" conjugation). For a schematic depiction of the trannulene principle, see "Structure B" in Figure 42. While in planar (or nearly planar) [10]annulene geometries the p orbitals are oriented perpendicular to the ring plane, the all-trans configuration forces the p orbitals in the direction toward the ring center [179]. At the B3LYP/6-31G* calculation level, the all-trans [10]annulene, the first member of the group of all-trans[n]annulenes or [n] trannulenes, has D 5d symmetry with equal C-C bond lengths, 1.412Å [179]. It has been verified that [4n + 2] trannulenes have aromatic character on the basis of geometric and magnetic criteria.
The fully delocalized D 10h and the bond alternating D 5h planar structures are local maxima lying 31.9 and 31.2 kcal/mol above the twist isomer (SCF/DZP values) on a quite flat energy surface in the vicinity of the two maxima [170,171]. The instability of these structures is due to the fact that the energy gain of the planar arrangement is overwhelmed by the energy required to deform the CCC angle to 144 • . A strategy has been proposed [185] to generate planar [10]annulenes analogues of benzene where cyclopropene or cyclobutene units are fused with the D 10h ring. Density functional results indicate that derivatives with attached five cyclopropene or cyclobutene units arrange as planar D 5h structures that are local minima on the energy surface [185].
Introducing one alkyne moiety, 1,2-didehydro [10]annulene is formed, characterized by several conformational minima and transition states [186]. The isomer resembling [10]annulene heart is the lowest energy conformer with a nearly planar and approximately bond equalized structure. Evidently, the replacement of one double with one triple bond relaxes the angle strain in this form. The CCSD(T) calculation at the MP2/cc-pVDZ optimized geometries identifies other three low energy isomers. The next-to-heart isomer, 6.8 kcal/mol higher, has a non planar C 2 structure similar to the [10]annulene twist. A second planar C 2v isomer with CC distances between 1.372 and 1.434Å lies 16.6 kcal/mol above heart. The third is a boat structure of C s symmetry with energy ≈ 17.7 kcal/mol above heart. All other conformational minima are at much higher energies.

Bridged [10] Annulenes
Experimentally, three types of tautomeric structures have been reported for bridged [10]annulenes, aromatic, polyolefinic and norcaradienic (see Figure 20). According to X-ray diffraction studies [187][188][189][190] 1,6-methano-and 1,6-epoxy [10]annulene have bond equalized C 2v geometries. As to 1,6-methano [10]annulene in particular the maximum value of the torsional angles along the ring is 35 • indicating that the structure is appreciably displaced from planarity [187]. Substituting the hydrogen atoms of the > CH 2 group with π acceptor units like −CN the 11,11-dicyano derivative takes on the norcaradienic structure [191][192][193]. Finally, the bulky trimethylsilyl substituent, −M e 3 Si, in the 2,5,7,10 positions favors the polyolefinic tautomer [194]. The structural data on 1,6-methano-, 1,6-epoxy [10]annulene and on the 11,11-dicyano derivative are in agreement with density functional calculations [195]. The geometries of Table 10 are minima on the potential energy surface, all the vibrational modes having real frequencies. In particular, the norcaradienic minimum of the dicyano derivative has a C 1 · · · C 6 distance 1.542Å (1.558Å, exp [193]) and enhanced bond alternancy ≈ 0.13Å (0.14Å, exp [193]). The delicate balance between aromatic and norcaradienic geometries depends on the nature of the bridging groups [196]. As the strain energy decreases from the > CH 2 , > O, > N H triplet to the > P H and > S pair, due to the reduction of the bridge angle, the transannular interaction increases and stabilizes the norcaradienic geometry. The valence aromatic/norcaradienic tautomerism has been modelled as a function of the C 1 · · · C 6 distance considering different substituents of the bridge H atoms [197]. The results have been related to minima of the cyclopropane and trimethylene fragments. Table 10.
Spectroscopic measurements on 1,6-methano [10]annulene, including absorption, fluorescence and magnetic circular dichroism [200][201][202], were reported since the synthesis [165] and discussed in relation to the transannular interaction. The molecule was viewed as an example of "arrested transition state" along the electrocyclic path of the ring closure [201]. The UV spectrum was assigned to the aromatic L b , L a , B b transitions, the absorption bands being at 3.44, 4.16 and 4.87 eV [200]. In the last years the vibrational and electronic properties of 1,6-methano-and 1,6-epoxy [10]annulene have been subjected to additional scrutiny and the two molecules have been characterized in more detail [83,195,198,199,203]. The infrared and Raman spectra of both molecules have been measured and a complete vibrational assignment was reported on the basis of ab initio density functional calculations of frequencies and infrared and Raman intensities [195]. The good agreement between the observed and calculated Raman spectrum is seen in Figure 21 for 1,6-epoxy [10]annulene. Medium resolution fluorescence spectra of 1,6-methano [10]annulene have been observed by means of the site selection technique exciting on the onset of the L b absorption band (origin and the lowest three vibronic bands, see Figure 22) [198]. With electronic origin at 25081 cm −1 the most active vibronic bands of the fluorescence spectrum lie at 323, 614 and 778 cm −1 from the origin. From these bandheads overtones and combination bands extend with clear Franck-Condon profiles.
The analysis has allowed to make inferences on the L b equilibrium structure on the basis of the common feature of these modes, i.e., the large atomic displacements on C 1 and C 6 . Therefore, it was proposed that upon electronic excitation the C 1 · · · C 6 distance increases so that the transannular interaction as well as the bond alternancy along the ring is reduced. Similar conclusions hold for 1,6-epoxy [10]annulene considering the fluorescence excitation spectrum under supersonic jet expansion [199]. The allowed contribution to the spectrum is due to totally symmetric modes (290, 652, 677 cm −1 ) mostly involving the epoxide bridge. The calculated equilibrium structures of 1,6-methano-and 1,6-epoxy [10]annulene in the L b state are reported in Table 10. The fluorescence excitation spectrum of 1,6-methano [10]annulene has been measured under supersonic expansion and consists of a rich vibronic band structure with rotational substructure starting from the electronic origin at 25154 cm −1 up to 4000 cm −1 excess energy [203]. On the basis of the rotational envelope it has been experimentally verified that the S 0 → L b transition dipole oscillates along the molecular long x axis. By a careful analysis of the highly resolved spectrum the assignment has been made of 13 fundamental vibrations of a 1 symmetry and 11 of a 2 symmetry, these latter being active through L b − L a coupling. According to this analysis the observed/calculated frequencies of the ground [195] and lowest excited state form two linearly correlated sets of values, thus enforcing the assignment.
The two-photon fluorescence excitation spectrum of 1,6-methano [10]annulene gives additional information on excited states below 5 eV [83]. The strong one-photon B b band is almost degenerate with an equally strong two-photon band, as shown in Figure 23.
The latter however cannot be similarly assigned. The vertical excitation energies of Figure 24 have been obtained with MCSCF/CAS(10,10)/6-31G calculations including the MC/QDPT perturbation treatment [204]. The complete set of one-and two-photon data reported in Table 10 show the remarkably good agreement with experiment if the strong two-photon absorption band is attributed to a fourth ππ * excited state of A 1 symmetry, S 4 (A 1 ).  Calculated and observed vertical transition energies (eV) of 1,6-methano [10]annulene. Left: MCSCF/CAS(10,10)/6-31G results; center: results from MC/QDPT perturbation treatment [204]; right: observed one-and two-photon transitions.
A second interesting point is raised by the observation of the Kekule' mode of b 1 symmetry, favoring the bond alternancy in the excited state. This mode increases from 1355 cm −1 (ground state value [195]) to 1540 cm −1 , as expected from the π distortive model of cyclic conjugated systems [52]. A steeper potential occurs along the distortion mode in the excited state justifying the occurrence of the equilibrium C 2v geometry [198].
Note that the tri-trans isomer has two inner hydrogen atoms above and one below the average ring plane. Among the conformational processes occurring in the tri-trans isomer those of Figure 26, i.e., (a) the exchange 1a 1a of the three inner with the three outer trans hydrogen atoms and (b) the enantiomeric process of the three inner hydrogens 1a 1a , have been investigated in detail [208].
Two transition states have been identified, T S2 and T S3 of Figure 27, 4.5 and 2.7 kcal/mol above CT CT CT . T S2 is a C 2 -symmetric transition state with an energy in close agreement with the observed [206] activation energy of the (a) process, 4.1 kcal/mol.
The T S3 transition state has C 1 symmetry and is responsible of the (b) process.
The conformation 1b of Figure 27, obtained by simultaneous rotation around all single bonds (D 3 symmetry, three trans C=C bonds perpendicular to the ring plane), suggested earlier as the transition state for this process [206], is a rather shallow minimum which may be reached crossing the T S1 transition state.
Two possible alternatives have been considered for sand as-indacene: either a bond equalized structure of D 2h /C 2v symmetry, i.e., 3 and 5 of Figure 28, or two interconverting valence isomers of C 2h /C 2v symmetry with localized bonds, i.e., 2a,2b and 4a,4b. In the second case 3 and 5 are transition states connecting 2a with 2b and 4a with 4b, respectively. Note, however, that while 2a,2b are degenerate configurations, 4a,4b are structures with different energies. The two most recent ab initio calculations on s-indacene [210,211] give structural data in good agreement with each other (see Table 11) and small energy differences between the D 2h and C 2h structures, 0.7 and 0.1 kcal/mol in favor of the localized C 2h structure. In the B3LYP/6-31G * study [211], the D 2h structure was found to be the transition state for double bond isomerization 2a → 2b, the bond alternating b 1g mode having imaginary frequency, 349i cm −1 . If ZP E is included, the υ = 0 energy of the D 2h transition state is lower than the υ = 0 energy of the localized C 2h state by 0.6 kcal/mol. Accordingly, s-indacene was referred as a "quasi-delocalized" molecule. In the other study [210] single-point CASPT2 energy calculations at the MP2/6-31G * geometries lower D 2h s-indacene below C 2h by 3.1 kcal/mol. The s-indacene ground state has been described as due to bonding of two allyl radicals with a central benzene moiety (regarded in turn to be formed by two allyl radicals).
Although s-indacene is a highly reactive molecule, the 1,3,5,7 tert-butyl derivative is a stable compound at room temperature [30]. X-ray diffraction data have been taken at room as well as at 100 K [30,31] consistent with a symmetric D 2h ring geometry. Coupling these data with previous results it appears that tert-butyl substituents play a minor role in determining the ring structure, limited to the increase of the 1-8a bond [212]. Due to the hindered rotation of these groups around the C ring − C but bond several conformers occur some of which are responsible of the near-infrared S 0 → S 1 absorption spectrum when this molecule is dissolved in solution [212,213]. Table 11. Upper: ground state geometries (bondlengths,Å) of sand as-indacene and energy difference ∆E (kcal/mol) between localized and delocalized structures of Figure 28. Center: ground state geometries (bondlengths,Å) of s-indacene (I), 1,3,5,7-tetramethyl-s-indacene (TMI) and 1,3,5,7-tetra-tert-butyl-s-indacene (TTBI). Lower: experimental ground state geometry of 1,3,5,7-tetra-tert-butyl-s-indacene and calculated geometries of the S 2 and S 6 states of s-indacene. The atomic numbering is indicated in Figure 28. 1,3,5,7-tetra-tert-butyl-s-indacene fluoresces from the second excited state, a property shared by few others organic molecules among which the first and most famous is azulene [215], with fluorescence quantum yield 9×10 −5 and radiative decay constant τ r (S 2 ) ≈ 10 ns [216]. There is a large energy gap between the onset of the S 0 → S 1 absorption, 9270 cm −1 , and that of S 0 → S 2 , 18380 cm −1 and no other energy levels are calculated in the intermediate region [214,216,217]. The relaxation dynamics of the S 1 and S 2 levels has been studied by means of picosecond transient absorption experiments [216]. Time constants of non radiative deactivation from S 1 and S 2 , 18 and 2.5 ps respectively, were measured. These times are much smaller than the shortest time constants of intersystem crossing in conjugated cyclic hydrocarbons, τ ISC ≥ 10 ns [218]. Internal conversion is the dominant relaxation process in these cases. The S 1 → S 0 internal conversion has been further characterized by the occurrence of a S 1 /S 0 conical intersection close in energy and geometry to the S 1 minimum [219]. Due to this an efficient non radiative decay channel is promoted thus leading to a short S 1 lifetime and lack of S 1 → S 0 fluorescence.
Excited states of the 1,3,5,7-tetra-tert-butyl-s-indacene have been probed by means of one-photon absorption and resonance Raman spectroscopy at low temperature [214]. As the excitation energy approaches the S 0 → S 2 origin the totally symmetric Raman peaks 390, 572 and 742 cm −1 are enhanced together with the multiplet below 300 cm −1 (see Figure 29). When the excitation energy is tuned across the S 0 → S 2 absorption profile both resonance Raman and fluorescence signals are observed (see Figure 30).
The Raman enhancement profile of the 390 cm −1 peak (corresponding to the calculated 431 cm −1 mode) has been measured as a function of the exciting wave number and determined from the calculation of the ground and S 2 equilibrium structures and of the consequent Franck-Condon factor. The results are compared in Figure 31. Recalling the "allyl" representation of s-indacene, the vibrational mode (390 cm −1 , exp; 431 cm −1 , calc) describes the elongation of the two side allyl units with respect to the benzene moiety (composed of two allyl units whose CC bonds elongate), inducing the equilibrium geometry of the S 2 state. The structure of the S 2 state, having approximate D 2h symmetry according to semiempirical methods [214], is in good agreement with expectations. Figure 30. Resonance Raman and fluorescence spectra of tetra-tert-butyl-s-indacene as a function of the exciting wavenumber across the S 0 → S 2 absorption region (bottom to top, c = 10 −6 M in isopentane/ether, 77K). Solvents bands are indicated by asterisks. Figure 31. Experimental Raman enhancement factor (crosses and squares) of the 390 cm −1 peak scanning across the S 0 → S 2 origin band and the first vibronic band of tetra-tert-butyl-s-indacene and calculated factor of the corresponding 431 cm −1 mode for the model system mono-tert-butyl-s-indacene.
Biphenylene has been extensively investigated in the ground and lowest ππ * excited states [220][221][222][223][224][225][226]. The infrared spectrum of biphenylene embedded in a Ar matrix was measured at 12 K, allowing the assignment of 35 fundamental modes (out of the total 54) by spectral analysis and comparison with ab initio calculation data [220]. As to the ππ * electronic transitions, it should be recalled that oneand two-photon spectroscopy are mutually exclusive techniques since biphenylene is a centrosymmetric molecule with D 2h symmetry. Starting from the ground state, excited states of ungerade symmetry are active in one-photon spectroscopy. The second transition, S 0 → S 2 , is long-axis polarized [223] (see Figure 32 for the reference system) with origin at 27690 cm −1 . At higher energy a strong absorption, equally long-axis polarized, is found at 40240 cm −1 for biphenylene in solution [224]. Although the S 1 (B 1g ) state is one-photon inactive a careful analysis of the vibronically induced spectrum in the lowest absorption region [222] locates the electronic origin at 23884 cm −1 . Two-photon active states of ππ * character belong to the A g and B 1g symmetry species. The two-photon spectrum of biphenylene in Figure 32 has been measured in the spectral range 225-370 nm combining fluorescence excitation and thermal lensing detection [225]. While the transition to the S 1 state, in principle two-photon active, is too weak to be observed, at higher energy the first excited A g state, 2A g , is seen around 265 nm and the second excited B 1g , 2B 1g , has its onset at ≈ 225 nm. The vibronic band structure between 340 and 310 nm is assigned to vibronically induced transitions to the second excited state. The S 0 → S 1 two-photon transition has enhanced intensity in 2,3,6,7-tetra-ethylbiphenylene [226] due to the pseudoparity selection rules and to the inductive perturbation of the ring ππ * density [227]. The electronic origin in the derivative is observed at ≈ 23100 cm −1 .

[14]Annulene
The lowest energy isomer of [14]annulene has the structure of a distorted pyrene perimeter. The molecule is non planar and centrosymmetric but with approximate C 2h symmetry according to X-ray data [228], as indicated in Figure 33. The CC bond lengths are in the range 1.350-1.407Å without apparent alternation. Non planarity is due to the steric overcrowding of the four inner hydrogen atoms which are directed toward the center of the molecule. Figure 33. Experimental [228], left, and calculated [20] (B3LYP/6-311+G * * , center; BHLYP/6-311+G * * , right) structures of [14]annulene. Bondlengths inÅ.
Recent density functional calculations on [14]-, [18]-and [22]annulenes have shown that the C 2h structure of [14]annulene, though a local minimum at the B3LYP level, is a transition state when functionals like BHLYP and KMLYP, having a larger Hartree-Fock component, are chosen [20]. Using these functionals the C 2h transition state lies 2.2 -2.5 kcal/mol above the C s minimum. The calculated 1 H N M R chemical shifts of the inner and outer hydrogen atoms match well the experimental values [6] assuming the C s structure for [14]annulene while this is not verified when the X-ray geometry is used. Other experimental techniques such as vibrational and electronic absorption spectroscopies should provide in principle evidences in favour of the C s structure for [14]annulenes. It has been however observed that these methods are not sensitive to the C 2h → C s structural variation [20]. In the following subsection the UV absorption spectrum of [14]annulene will be compared with those of bridged [14]annulenes whose bond equalized structures are well established both on experimental and theoretical grounds.
The molecule undergoes in solution the bond shifting reaction of Figure 34 with the equilibrium state strongly displaced toward the configurational CT CT CT T isomer. The activation energy of this process was found to be 21 kcal/mol [6]. Since ∆ trans, i.e., the difference of trans double bonds between reactant and product, is 1, this means that one trans double bond isomerizes to cis. The reaction mechanism consists of three steps (see Figure 35) and requires in the final step a Möbius antiaromatic transition state [9].
Transition states of annulenes may be classified according to the Hückel or Möbius topology [229,230]. In the simple example of bond shifting in cyclobutadiene and cyclooctatetraene the number of cis double bonds is conserved and the interconversion process occurs through a Hückel antiaromatic transition state. In general, the number of trans double bonds are the reference index for bond shifting in annulenes: the reaction will occur with Hückel topology when ∆ trans is zero or an even number and with Möbius topology when ∆ trans is an odd number. For the reaction shown in Figure 35 ∆ trans = 1. The reaction requires a Möbius bond shifting step and since [14]annulene is a (4n + 2) π-electron system the bond equalized transition state is antiaromatic, T S1 in Figure 35. The three step reaction CT CT CT T → CCT CT CT has an energy barrier in the bond shifting step of 19.3 kcal/mol [9]. Figure 34. Cis/trans isomerization reaction of [14]annulene in solution at room temperature.

Bridged [14]Annulenes
Three types of bridged [14]annulenes have been synthesized, based on anthracene, pyrene and dicyclopentaheptalene perimeters [10]. Here we consider only the anthracene type molecules of Figure 36. The two methano bridges may be arranged in a synor in a antigeometry: in the synisomer 6 the ring is expected to be approximately planar while in the anti 7 the ring takes on a polyolefinic structure [16].
The Raman spectra depend on the excitation wavelength. In Figure 37 the 1064 nm excitation corresponds to the normal Raman spectrum of 9. With the 647.1 nm excitation all the Raman peaks above 800 cm −1 show decreased intensities with respect to those below 800 cm −1 . The peaks at 251, 420 and 615 cm −1 are assigned to modes of anthracene parentage, the lowest to an out-of-plane and the other two to in-plane modes, where the displacements of the bridgehead atoms enlarge the ring size [239]. Therefore, it has been assumed that the equilibrium geometry of the lowest excited state (L b ) has increased transannular distances. The analysis of high- [240] and medium-resolution [241] fluorescence spectra of the two systems confirms the hypothesis.  In the absorption spectra of 8 and 9, shown in Figure 38, the band maxima of the lowest electronic transition shift from 480 nm in 9 to 550 nm in 8. Reasonably, this effect is correlated with the increased ring planarity of 8 with respect to 9. Considering also the absorption spectra of [14]annulene [242] and of dimethyldihydropyrene [243], a bridged [14]annulene with pyrene perimeter, the band maxima further shift to ∼ 600 nm in the first case and to ∼ 640 nm in the second. Then, both these systems in solution are more planar than 8.
The absorption spectra [241,244] of 8 and 9 (see Figure 38) offer the good opportunity to investigate about the L a and B a geometries by means of resonance Raman spectroscopy. CC stretching modes of 8 are preferentially enhanced exciting on the onset of the L a transition (see Figure 39) while only the 234 cm −1 mode shows prominently with excitation at 350 nm on the foot of the B a transition. Qualitatively, the latter result points to the similarity between the B a and L b equilibrium structures. As indicated by the MCSCF/CAS(10,10) data of Table 12, the common feature of the two states is the increase of the transannular distances with respect to the ground state. This structural change is responsible for the strong activity of the 234 cm −1 mode in 8 and of the 251 cm −1 mode in 9. On the contrary, the slight decrease of these distances in the L a together with the enhancement of bond alternation favors an anthracene-like equilibrium structure for the L a state of both molecules. Table 12. Structural data (bondlengths,Å) of syn-1,6:8,13-bismethano [14]annulene (6), 1,6:8,13-ethanediylidene [14]annulene (8) and 1,6:8,13-propanediylidene [14]annulene (9). Upper row: ground state B3LYP/cc-pVDZ results [239]. Central and lower rows: geometries of the ground and of the lowest excited states (8, central; 9, lower); MCSCF/CAS(10,10)/6-31G results [241].  (a) X-ray data [233]. (b) X-ray data [235], from the 15,16-dimethyl derivative of 8. (c) X-ray data [236].

[16]Annulene
The synthesis of [16]annulene was first achieved starting from the dehydro [16]annulene derivative [245] and later from the cyclooctatetraene dimer [246]. At temperature below −67 • [16]annulene is a mixture of two configurational isomers, 10a and 10b of Figure 40, 10a being the major (80 per cent) component [6,247]. In the solid state 10a is the stable arrangement of [16]annulene in a nearly S 4 bond alternating structure [248,249]. The optimized structures of the two isomers, CT CT CT CT (10a) and CT CT T CT T (10b), at the BHHLYP/6-311+G * * calculation level are reported in Figure 40. [16]annulenes 10a and 10b are the first 4n π annulenes whose bond alternating structures are not overwhelmingly dominated by angle strain or steric hindrance of the inner hydrogens. In a sense, the UV-vis absorption spectrum of [16]annulene should be typical of 4n π systems. Two absorption bands are observed at 446 and 293 nm with strongly different extinction coefficients, both assigned to the 10a isomer [250]. The weakest band at 446 nm has been interpreted in S 4 symmetry as due to the vibronically induced 2b(HOM O) → 3b(LU M O) transition. The two orbitals 2b and 3b arise from splitting of the non bonding e 4u degenerate pair due to symmetry lowering from the fully planar D 16h to S 4 symmetry. The strongest band at 293 nm is due to the allowed 2b(HOM O) → 3e(LU M O + 1) transition. Comparing with more simple 4n π systems, the first two excited states of cyclobutadiene and pentalene are of g symmetry [23,26,28,29,99] (and, as a consequence, the transition to these states from the ground state not allowed) while in [16]annulene only the lowest excited state is totally symmetric (and then the S 0 → S 1 transition not allowed). The difference follows from a property of cyclobutadiene and pentalene not shared by [16]annulene, i.e., the occurrence of low energy doubly-excited configurations determining the second excited state in these molecules.  [16]annulene isomers calculated at the BHHLYP/6-311+G * * level [8]. Selected bondlengths (Å) and dihedral angles (•) are shown.
The 10a isomer lies below 10b by ≈ 1 kcal/mol, in agreement with experimental data at low temperature [254]. The conformational change 10a → 10a occurs in a stepwise mechanism following two different pathways with calculated activation energies 8 -10 kcal/mol. If the indicated trans double bonds of 10a are called T 1, T 2, T 3 and T 4 consecutively, in one case the rotation around T 1 is followed by that around the adjacent T 2 trans double bond, in the other by that around the opposite T 3 trans double bond. The configurational change 10a → 10b occurs with isomerization of one trans double bond to cis and then with Möbius topology. The key step of the process is a twist-coupled bond shifting during which a Möbius aromatic transition state is crossed [8]. The overall process then consists of (1) the conformational change around a trans double bond, (2) the Möbius step and (3) the final conformational change to 10b. The highest activation barrier is found for the second step and ranges from 16.7 kcal/mol (BHHLYP/6-311+G * * ) to 6.9 kcal/mol (B3LYP/6-311+G * * ). Figure 42. Schematic view of the synthetic strategy used to stabilize the C 2 Möbius structure of [16]annulene [255].
Recently, the synthesis of a stabilized [16]annulene isomer with Möbius topology has been reported by cycloaddition of a normal aromatic structure with a "belt-like" in-plane structure, as schematically depicted in Figure 42 [229,255,256]. While the Möbius rings of parent [16]annulene are at higher energies with respect to 10a and 10b [257], the stabilization imposed by the belt reverses the energy ordering placing two Möbius rings, of C 1 and C 2 symmetry, at the lowest and nearly equal energy [255]. The synthesized Möbius ring lies ∼ 2 kcal/mol above and has the structure anti − tZcZcZt shown in Figure 43 with C 2 symmetry [255]. It has been argued [258] about aromaticity of this system, on the basis of the bond alternation resulting from crystal data [255] and of large dihedral angles inhibiting the overlap between π orbitals. In contrast, the principal component analysis performed on the 153 most stable structures of [16]annulenes and the 25 most stable Möbius stabilized [16]annulenes indicates that the compound of Figure 43 belongs to the cluster of aromatic Möbius [16]annulene systems [256]. A detailed account of the discussion about aromaticity of the Möbius [16]annulene of Figure 43 can be found elsewhere [230]. Figure 43. Structure of anti − tZcZcZt [16]annulene in the Möbius topology. Bondlengths inÅ are from X-ray data [255]. Dihedral angles (•) are from DFT calculations (B3LYP/6-31G*) [255]. Nomenclature: the anti configuration is relative to the two H atoms at the quinoid double bonds; c and t denote the s-cis or s-trans geometry between two adjacent double bonds; E and Z the stereochemistry of the three central double bonds. For the sake of clarity the benzo groups fused with bonds marked with a b are omitted [230].

[18]Annulene
The molecule occupies a special position in the annulene series being the next higher homologue [259] of benzene which allows, though does not require, the highly symmetrical D 6h structure. According to X-ray diffraction data at room and 111 K [260,261] the planar strainless D 6h geometry of [18]annulene has twelve transoid and six cisoid CC bonds, 1.385 and 1.405Å long, respectively [261], without bond alternation (see Figure 44). Figure 44. The structure of [18]annulene: left, experimental [260,261], in parentheses, and calculated B3LYP/6-311+G * * [20] CC bond lengths; right, calculated [20] KMLYP/6-311+G * * CC bond lengths in C 2 symmetry. Bond lengths inÅ. The calculated structure on the right is slightly non planar with C 2 axis as indicated. No dihedral angle has been reported [20].
Along the years it has been debated about the stability of the bond alternating D 3h structure with respect to D 6h . In fact, semiempirical and ab initio HF calculations indicate a greater stability of the former structure [112,262,263]. On the contrary, using density functional and MP2 methods to take into account the correlation energy the D 6h structure is preferred with calculated bond lengths in substantial agreement with experimental data [19,264,265]. More recently, a bold statement has been advanced challenging the currently accepted structural information on [18]annulene [20]. On the basis of the gross disagreement between observed and calculated proton chemical shifts assuming the D 6h geometry and, in addition, considering that correlated methods overestimate delocalization in annulenes [173,174], functionals having a larger HF component, such as KMLYP and BHHLYP, have been used. With both functionals the D 6h geometry results to be a transition state between equivalent minima structures of C 2 symmetry (see Figure 44). The proton chemical shifts of the inner and outer hydrogen atoms calculated starting from the C 2 structure have been found in excellent agreement with experiment [6,266]. Static and dynamic crystal disorder may explain, according to this study [20], the difference with respect to the X-ray indications. For example, it has been noted that a 50/50 superposition of the two C 2 structures is close to the reported X-ray data [260,261] and therefore in principle this possibility cannot be excluded. Second, since the calculated barrier between the two interconverting minima is around 3 kcal/mol, the symmetrization of [18]annulene should take place rapidly even at 111 K, resulting in a time-average D 6h structure [20].
The proposed failure of crystallographic studies to determine the structure of [18]annulene has been subjected to critical reconsiderations [267]. First, the superposition of two D 3h (or C 2 ) bond alternating structures with ∆r = 0.1Å between C −C and C = C bondlengths and rotated 60 • gives an orientational disorder in the crystal of the order of 0.0025Å 2 for the twelve outer C atoms. In contrast with benzene, a disorder effect of this size can be detected, for instance performing variable temperature neutron diffraction studies [108]. Second, the experimental proton chemical shifts are relative to [18]annulene in solution. Given the small energy difference between the C 2 and D 6h structures, packing interactions in the crystal may favor the planar non alternating D 6h structure. It has been concluded [267] that further experimental and computational effort should be made in order to confirm the C 2 structure for [18]annulene.
Finally, although spectroscopic data do not allow the unequivocal assignment of the [18]annulene structure among geometrical alternatives [20], it is worth noting the better agreement with experimental results of infrared [164] and UV-vis [268] spectroscopy when calculations are made starting from the delocalized D 6h rather than from the localized D 3h structure [268,269]. Taking advantage of the occurrence of the inversion center for D 6h symmetry and of its absence for C 2 symmetry, a careful comparison of the infrared and near-infrared Raman spectra of the [18]annulene crystal at very low temperature appears to be highly recommended in order to solve the structural issue.

Porphyrins
Porphyrin and substituted porphyrins are molecules of enormous interest for their role in biological processes such as light absorption in green plants and oxygen transport in aerobic organisms [270]. They are also chemical structures of relevance by their own and accordingly in this Section results on structures and decay mechanisms of porphyrin and closely related species will be presented and discussed in relation to the annulene character of the macrocycle. Diprotonated porphyrin has a non planar structure mimicking that of non planar porphyrinic prosthetic groups [271]. Tetraoxaporphyrin dication, the oxygen analogue of porphyrin, has the highly symmetrical D 4h structure [272,273] which may be modulated through five redox stages [274]. The structural parameters of the three molecular species have been the object of a large number of computational studies all of which indicate that the stable forms of porphyrin and tetraoxaporphyrin have D 2h and D 4h geometries, respectively, when electron correlation is taken into account [275][276][277][278][279][280][281][282]. The experimental [272,283,284] and calculated [275] bondlengths of the three macrocycles are reported in Figure 45.
As to the diprotonated species, the structure of Figure 45 is the lowest energy conformer with D 2d symmetry. Other two stable conformers with approximate C 2h and C s symmetry, the so-called quasi-wave and quasi-ruffled isomers [285], have been found by calculations [275,286] ≈ 5 kcal/mol above the D 2d conformer. In qualitative discussions on porphyrins the concept of main conjugation path embedded in the macrocycle and responsible of the electronic absorption has been advanced [287,288]. The proposal may be useful also to estimate how the perturbation introduced by the >O, >NH and >N − groups affect the macrocycle structure. The oxo bridge, >O, has the weakest effect on the outer ring as evidenced by the fact that the CO distance is close to the single C-O bond. Tetraoxaporphyrin dication may thus be considered a slightly perturbed [20]annulene system with 18 π electrons. Also the D 2d structure of diprotonated porphyrin, though not planar, has ring bondlengths appreciably similar to those of tetraoxaporphyrin with a growing interaction of the macrocycle with the >NH bridges (the CN distance is roughly intermediate between single and double CN bonds). On the other hand, a quantitative determination of the main conjugation path in the porphyrin macrocycle has been proposed on the basis of structural data relative to a quite large number of porphyrin derivatives [289]. The porphyrin π system consists of the two pyrrole rings and of the "internal cross", i.e., the [16]annulene dianion with 18 π electrons [289]. Figure 45. Calculated (B3LYP/cc-pVDZ) [275] and (in parentheses) experimental bond lengths (Å) [272,283,284] of the stable structures of porphyrin (PH 2 , D 2h ), diprotonated porphyrin (PH 2+ 4 , D 2d ) and tetraoxaporphyrin dication (TOXP 2+ , D 4h ). The experimental bond lengths are in the three cases average values of distances symmetrically equivalent in D 2h , D 2d and D 4h symmetry, respectively. The reference axis system is also shown on the top of the Figure  A large amount of vibrational data have been acquired on porphyrin from normal and resonance Raman [290,291], infrared [292] and fluorescence/phosphorescence [293] spectroscopies and the assignment of the fundamental modes has been established in great detail [290,294]. The infrared and Raman spectra of the tetraoxa-and diprotonated porphyrin have been equally obtained and completely assigned making use of ab initio calculated frequencies and intensities and correlation with porphyrin [275]. The near resonance Raman spectra of porphyrin and diprotonated porphyrin excited at 406.7 nm and the normal Raman spectrum of tetraoxaporphyrin exciting at 785 nm are reported in Figure 46 and  The absorption spectra of the three macrocycles may be qualitatively interpreted on the basis of the main conjugation path. The conjugation path of the [20]annulene dication represents the suitable choice for the tetraoxa-and diprotonated porphyrin [288]. By means of the PMO theory [57] it is found that the >O and >NH perturbations split the doubly-degenerate HOMO orbitals of the regular [20]annulene polygon into two orbitals of different energy, HOMO and HOMO-1, while shift uniformly the doubly-degenerate LUMO orbitals to higher energies. Two transitions are expected, HOMO → LUMO and HOMO-1 → LUMO, which correspond to the weak Q and intense B (or Soret) bands of Figure 48.
The strong perturbation due to the >N − bridges further splits the LUMO orbitals into a doublet, thus giving the well known four-orbital model of porphyrins [295,296]. Four electronic transitions are expected, all of them active in one-photon spectroscopy, the two at lower energy denoted Q x , Q y while the two higher B x , B y coalesce in a single band (see Figure 48). A quantitative treatment of the vertical excitation energies for the three molecules can only be performed with ab initio calculation methods [297][298][299][300][301][302][303]. The results from three recent calculations relative to porphyrin and presented in Table 13 differ from experiment [299,304] by no more than 0.2 -0.3 eV not only for the Q x , Q y and B x , B y pairs but also for the weak transitions, usually denoted N x , N y and L y , and observed at higher energies with respect to the Soret band. Figure 48. Absorption spectra of PH 2 , PH 2+ 4 and TOXP 2+ solutions, c ≈ 10 −5 M, room temperature, from bottom to top. The solvents are benzene/cyclohexane 1:10 for PH 2 , benzene/cyclohexane 1:10 added with CF 3 COOH (5 % in volume) for PH 2+ 4 , concentrated HClO 4 for TOXP 2+ . Vertical arrows indicate pump wavelengths in our pump-probe experiments in the three cases.
Exciting into the Soret absorption band of tetraoxaporphyrin as indicated in Figure 48 radiative and non radiative decay processes are promoted [305,306]. The fluorescence emissions from the first and second excited state, Q → S 0 and B → S 0 , are easily observed with quantum yield 0.11 and 0.035, respectively [305]. The upper state fluorescence is favored by the sharpness of the Soret band having a considerably high extinction coefficient and by the absence of intermediate states in the large energy gap between B and Q. The B lifetime could only be estimated to be ≤ 100 ps, due to the instrumental resolution of the time-resolved fluorescence apparatus, and indirectly evaluated from quantum yield and radiative time constant to be ≈ 40 ps. Transient absorption pump-probe experiments have been set up in the femtosecond time domain [306] following experimental procedures described in detail [160,307,308] in order to characterize the non-radiative relaxation channels from the B state. In general, transient absorbances ∆A arise from four different contributions [307], excited state absorption (ESA), stimulated emission (SE), stimulated Raman gain (SRG) and bleaching (B). All these contributions are clearly seen in the transient spectrum of Figure 49.
The decay kinetics of the 405 and 425 nm ESA bands in Figure 50 identifies a single rise time of 25 ps, interpreted as the time constant of the population increase in the Q state. Successively, the transient absorbance decays with time constant 7.2 ns, i.e., the lifetime of the Q state. By comparison with the estimated B lifetime, ≈ 40 ps, it has been concluded that, within the experimental error, all the tetraoxaporphyrin molecules excited into the B-state decay by internal conversion into the Q state. Table 13. Transition energies and oscillator strengths (∆E, eV; f ) of porphyrin to the lowest excited states, ππ * unless otherwise specified. Exp: gas-phase [304] and solution band maxima ∆E [299] and f values [300]. Calculated transition energies under all other headings: DF/SCI results [299] (reference system as in Figure 33); TD/DFT, results [298]; CASPT2 results [297]. The transient absorption from the Q x ; Q y bands and the relaxation dynamics from the Soret band have been investigated in substituted porphyrins [309][310][311][312]. In tetraphenylporphyrin the internal conversions (B x , B y ) → Q y and Q y → Q x take place in less than 50 and 100 fs, respectively [310,311]. Later, three decay channels are activated; (a) intramolecular vibrational redistribution by dephasing and elastic collision with the solvent with time constants 100-200 fs and 1.4 ps, (b) thermal equilibration by energy transfer to the solvent in the temporal range 10-20 ps and (c) intersystem crossing Q x → T 1 with time constant 12 ns. In another experiment [312], exciting into the Q x ; Q y bands of the tetra-p-tolyl derivative of porphyrin and probing with infrared pulses, states of gerade symmetry in the vicinity of the Soret band have been observed by transient absorption spectroscopy.
Going to the unsubstituted molecule, the non radiative processes follow a decay pattern similar to that of tetraphenylporphyrin [306]. The femtosecond transient absorption results on porphyrin are presented in Figure 51 while Figure 52 shows the decay kinetics of the 412 nad 465 nm ESA bands.
After the unresolvably short (B x , B y ) → Q y internal conversion the Q y → Q x process is observed with time constant 140-150 fs, a rise/decay time for the 412/465 nm ESA band. As expected, this time constant is slightly longer than in substituted porphyrins. Also the subsequent fast processes with time constants 1.8 ps and 25 ps slow down in porphyrin with respect to the tetraphenyl derivative [311].   Diprotonated porphyrin is a case more complex than porphyrin and tetraoxaporphyrin. Although the energy gap between the B and Q states of the diprotonated species is comparable to that of tetraoxaporphyrin, efficient decay processes are promoted due either to the torsional flexibility of the non planar structure [313][314][315] or to the formation of the complex with the trifluoroacetic acid used for the preparation of the diprotonated species, i.e., (PH 2+ 4 )(CF 3 COO − ) 2 , as suggested in similar cases [316]. The complex has charge transfer states interspersed between B and Q and provides a deactivation mechanism more efficient than the direct B → Q internal conversion channel. According to the B3LYP/cc-pVDZ calculation the molecular complex isomerizes between the two equivalent A 1 and A 2 structures of Figure 53, expressed as (PH 2 )(CF 3 COOH) 2 , with an activation barrier of 0.05 eV (≈ 1.1 kcal/mol).
The UV-vis absorption spectra of porphycene and of a large number of porphycene derivatives have been comprehensively investigated applying the charged 18-electron [20]annulene planar model for a qualitative discussion of the Q and Soret bands [288]. Unlike the porphyrin case, by perturbing the [20]annulene ring with > N H and > N − bridges such that to obtain porphycene, one of the two degenerate LU M O orbitals of the model perimeter shifts to higher energy while the second has unaltered energy. Further, the effect of perturbation on the two degenerate HOM O orbitals is approximately equal so that they remain accidentally degenerate. Accordingly, the two Q transitions of porphycene,  The inner hydrogen migration between trans structures (see Figure 54, bottom) is common to porphyrin and porphycene [288]. Ground and excited state tautomerism in porphycene has been experimentally investigated by several techniques [323][324][325][326]. While this process has a stepwise mechanism in porphyrin [327], in porphycene the tautomerism has been described as a synchronous double hydrogen tunneling activated by a low frequency mode whose displacement affects the NH· · ·N distances [320]. By supersonic jet spectroscopy the tunneling splitting has been observed in the fluorescence excitation spectrum of the isolated molecule in the region of the (0-0) transition [325]. Also, the inner hydrogen transfer between the trans and the cis-B geometry is related to the efficient radiationless deactivation in excited porphycene [328]. It has been in fact proposed that crossing of the potential energy curves of the ground and lowest excited states for these two isomers occurs along the hydrogen transfer coordinate, thus giving origin to a conical intersection at a geometry close to that of the lowest cis-B excited state [328]. Additional information on excited state dynamics and relaxation in electronically excited porphycene is provided by transient absorption spectroscopy [329]. Similarly to the porphyrin case, the relaxation scheme for the molecule in solution excited into the Soret band involves ultrafast intramolecular decay processes and a slower thermal equilibration with the surrounding solvent.

Summary
In this paper the most significant structural aspects of annulene and model annulene systems relative to the ground and lowest excited states have been reviewed. The research activity of our group has been centered on one side on antiaromatic systems such as pentalene and s-indacene, and on the other on aromatic bridged annulenes and porphyrins. More recently, energy levels and decay mechanisms have been studied in planar and non planar porphyrins. The relaxation processes have been probed by means of time-resolved picosecond fluorescence and transient absorption femtosecond spectroscopy. Radiative and non radiative processes have been observed in unsubstituted porphyrin, diprotonated porphyrin and tetraoxaporphyrin dication describing the decay channels activated after excitation in the Soret band.
In an effort to present our results in the context of the aromaticity/antiaromaticity theory and to provide basic notions on the structural properties of annulenes the review has supplemented additional information also on these systems. A common feature of higher annulenes is the occurrence of non planar isomers and of thermal or photochemical valence isomerization among them. Two other fields of intense activity are worth mentioning. The design, almost entirely theoretical until now, of molecules exhibiting Möbius topology has stimulated growing interest in the chemistry of large annulenes. Möbius transition states have been proposed and computationally justified for the isomerization processes in these systems. The first Möbius ring has been synthesized, thus making feasible the direct comparison between experiment and theory. On the other hand, the structural transition from the delocalized to a localized valence tautomeric form in the series of [4n + 2] annulenes is a second issue giving evidence of the distortive nature of π electrons in conjugated systems. In this respect the discussion on the delocalized/localized ground state structure of [18]annulene takes on particular importance.
The review has covered in depth as much as possible the structural properties of aromatic/antiaromatic systems. Given the amplitude of the aromaticity subject, these properties do not exhaust the theme. As an example, the magnetic criteria of aromaticity have not been tackled, for which more complete reports are recommended [4,330,331].