Triplet State Baird Aromaticity in Macrocycles: Scope, Limitations, and Complications

The aromaticity of cyclic 4nπ-electron molecules in their first ππ* triplet state (T1), labeled Baird aromaticity, has gained growing attention in the past decade. Here we explore computationally the limitations of T1 state Baird aromaticity in macrocyclic compounds, [n]CM’s, which are cyclic oligomers of four different monocycles (M = p-phenylene (PP), 2,5-linked furan (FU), 1,4-linked cyclohexa-1,3-diene (CHD), and 1,4-linked cyclopentadiene (CPD)). We strive for conclusions that are general for various DFT functionals, although for macrocycles with up to 20 π-electrons in their main conjugation paths we find that for their T1 states single-point energies at both canonical UCCSD(T) and approximative DLPNO-UCCSD(T) levels are lowest when based on UB3LYP over UM06-2X and UCAM-B3LYP geometries. This finding is in contrast to what has earlier been observed for the electronic ground state of expanded porphyrins. Yet, irrespective of functional, macrocycles with 2,5-linked furans ([n]CFU’s) retain Baird aromaticity until larger n than those composed of the other three monocycles. Also, when based on geometric, electronic and energetic aspects of aromaticity, a 3[n]CFU with a specific n is more strongly Baird-aromatic than the analogous 3[n]CPP while the magnetic indices tell the opposite. To construct large T1 state Baird-aromatic [n]CM’s, the design should be such that the T1 state Baird aromaticity of the macrocyclic perimeter dominates over a situation with local closed-shell Hückel aromaticity of one or a few monocycles and semilocalized triplet diradical character. Monomers with lower Hückel aromaticity in S0 than benzene (e.g., furan) that do not impose steric congestion are preferred. Structural confinement imposed by, e.g., methylene bridges is also an approach to larger Baird-aromatic macrocycles. Finally, by using the Zilberg–Haas description of T1 state aromaticity, we reveal the analogy to the Hückel aromaticity of the corresponding closed-shell dications yet observe stronger Hückel aromaticity in the macrocyclic dications than Baird aromaticity in the T1 states of the neutral macrocycles.


■ INTRODUCTION
In recent years, it has become accepted that many compounds with 4nπ-electron cycles are aromatic in their lowest electronically excited states, a characteristic now labeled as Baird aromaticity. 1−19 Different types of compounds have been examined; monocyclic, polycyclic, heterocyclic and macrocyclic. Yet, the limitations and complications of the excited state Baird aromaticity concept need to be assessed in more detail so that we can avoid using it on molecules for which it is not applicable. A critical evaluation is also crucial for the successful and efficient implementation of Baird aromaticity in design of molecules with targeted properties, e.g., for applications in organic electronics or spintronics. Herein we explore macrocyclic compounds which are oligomers of four different monocyclic compounds; two of the monocycles used are Huckel-aromatic in their electronic ground states (S 0 ) while the other two are nonaromatic with 1,3-diene segments. Our focus is on the potential T 1 state Baird aromaticity of the macrocycles. What type of small monocycles as monomer units provide for strong Baird aromaticity in large macrocycles ? Peeks, Anderson and co-workers explored the potential excited-state Baird-aromatic and -antiaromatic characters of large π-conjugated porphyrin nanorings with 70−112 πelectrons in the main macrocyclic conjugation pathway (MCP), including both Huckel-antiaromatic ([4n]) and -aromatic ([4n + 2]) cases. These nanorings are exceptional as they are composed of a series of Huckel-aromatic porphyrins linked via butadiynyl segments. Yet, it was observed that the extent of T 1 state Baird aromaticity assessed through computations varies strongly with choice of DFT functional. 20 On the other hand, for [n]cycloparaphenylene nanohoops ([n]CPP's) 21−24 in their T 1 states it has been considered that long-range corrected functionals are preferred. 25 When going to gradually larger [n]CPP's in their T 1 states calculated with UCAM-B3LYP, structures changed from circular with fully delocalized triplet excitations to egg-shaped with the excitations localized to a subset of the p-phenylene units at n > 5.
Herein, we argue that in the T 1 state there is a tug-of-war between dominance of macrocyclic Baird-aromatic character of cyclic 4nπ-electron MCP's and dominance of closed-shell Huckel-aromatic character of one or a few monocycles which forces the triplet diradical character to localize to a part of the macrocycle ( Figure 1A). The outcome of this tug-of-war should vary with the type and size of a macrocycle and its monocyclic constituents; e.g., macrocyclic Baird aromaticity should be sustained more easily in [n]cyclooligo(2,5-furan)s ([n]CFU) than in [n]CPP's as the closed-shell Huckel aromaticity of furan is lower than that of benzene. 26 Thus, one may argue that [n]CFU's are better suited for macrocyclic Baird aromaticity than [n]CPP's. Here it is notable that a [n]CFU derivative has been synthesized and studied experimentally, 27 and the parent [n]CFU's have been explored through quantum chemical computations revealing that those with n < 8 display planar or nearly planar S 0 state geometries and have electronic properties which make them interesting for organic electronics. 28 The effect of heteroatoms on the T 1 state Baird aromaticity of [6]CPP has been reported in a computational study where the phenylene units were changed successively to 1,4-azaborines, 29 a monocycle which is less Huckel-aromatic in S 0 than benzene. 30 Yet, such replacements in the cases investigated led to attenuated Baird aromaticity according to nucleus independent chemical shifts (NICS's) calculated with the PBE0 functional. Can other heterocycles lead to increased Baird aromaticity?
An interesting description of Baird aromaticity is the description of Zilberg and Haas in which the T 1 state of a neutral 4nπ-electron cycle is viewed as that of the corresponding Huckel-aromatic dication with 4n − 2 spinpaired π-electrons plus two nonbonding same-spin πelectrons. 31 This description was utilized to rationalize the aromaticity/antiaromaticity switch that occurs when going from S 0 to T 1 in a [26]-and [28]porphyrinoid pair. 13,32 Indeed, the potential Baird aromaticity of a neutral [n]CPP in its T 1 state could therefore closely resemble that of the corresponding [n]CPP dication ( Figure 1B), species for which in-plane Huckel aromaticity with 4n − 2 π-electrons has been demonstrated. 33 Yet, is the aromatic character stronger or weaker in the T 1 state of [n]CPP than in the closed-shell singlet state [n]CPP 2+ , or is it similar? Furthermore, is there a similarity to the dianions which are Huckel-aromatic with 4n + 2 π-electrons? The [n]CPP dianions have been explored earlier, both experimentally and computationally, 34−36 and are found to exhibit some aromaticity according to the computations. 34 Another issue is how to relate [n]CPP and [n]CFU with suitable large monocycles? As found by Beans, Fowler and Soncini for cyclic C 16 , cyclo[4n]carbons are interesting because their two orthogonally oriented 4nπ-electron cycles allow them to be doubly Baird-aromatic in their lowest quintet states (Qu 1 ), 37 in contrast to the S 0 state where they are doubly Huckel-antiaromatic until a certain size. 38−40 For the T 1 state, it was concluded that the out-of-plane π-conjugated path of C 16 displayed Baird aromaticity while the in-plane path displayed Huckel antiaromaticity. A situation in the T 1 state with two similarly strong counteracting ring-currents, one diatropic (Baird-aromatic) and one paratropic (Huckelantiaromatic), could lead to a NICS value that corresponds to a nonaromatic situation. Indeed, such seemingly nonaromatic characters were recently observed from computations of cyclic C 16 , C 20 . and C 24 . 41 Furthermore, after the recent experimental generation of cyclo [18]carbon by Anderson and co-workers, 42 [9]CPP from [9]CPP 2+ by addition of two nonbonding same-spin π-electrons to [9]CPP 2+ , which has a Huckel-aromatic MCP with 4 × 9−2 (i.e., 4 × 8 + 2) spin-paired πelectrons.
The Journal of Physical Chemistry A pubs.acs.org/JPCA Article macrocyclic MCP ( Figure 2) one arrives at a [n]CPP with n para-phenylene (PP) units and at a [n]CFU with n 2,5-linked furan units (FU). How will the Baird aromaticity change upon this gradual transformation? It is noteworthy that the MCP involves in-plane orbitals in [n]CPP while it involves out-ofplane orbitals in [n]CFU's with n < 8 (cf., trannulenes and annulenes, respectively). As [5]CPP is the smallest [n]CPP that has been synthesized and explored experimentally, 45 we start our analysis at the C 20 cycle and replace one, two, or five 1,3-butadiyne segments with M units, reaching, e.g., [5]CPP. This [n]CPP has the narrowest HOMO−LUMO gap among the [n]CPP's synthesized and it has a benzenoid structure in its S 0 state. 33 A NICS(0) iso value of −1.3 ppm in this state indicates a nonaromatic character with respect to the macrocyclic conjugation path. In our analysis we consider four different M units ( Figure 2); the 1,4-linked benzene (p-phenylene, PP), 1,4-linked cyclohexa-1,3-diene (CHD), 2,5-linked furan (FU), and 1,4-linked cyclopentadiene (CPD). The compounds are labeled as C 20−4n [n]CM with M as either PP, CPD, CHD, or FU, and n as their number. The [n]CCPD and [n]CCHD are macrocycles with just single 4nπ-electron paths and, accordingly, monocycles from an electronic perspective. In macrocycles with two M units these were placed at maximal distance from each other.
A number of questions can now be posed: It is known that the extent of Baird aromaticity varies with computational method, 20 but how does it vary between the different macrocycles? Also, to what extent is there an agreement between various types of aromaticity indices on the magnitude of Baird aromaticity of macrocycles? Discrepancies between magnetic and electronic indices have earlier been highlighted, and it was concluded that (anti)aromaticity assessments exclusively based on the magnetic aspects are not justified. 46 Furthermore, based on the Zilberg−Haas model 31 one may ask how similar in aromatic character are the closed-shell S 0 state dications (e.g., [n]CPP 2+ and [n]CFu 2+ ) and the corresponding neutral triplet states? One may also consider the similarity to the corresponding S 0 state dianions. The overarching goal of this study is to find routes to large macrocycles which to high probability are triplet state Baird-aromatic, not only in silico but also in reality.

■ RESULTS AND DISCUSSION
The investigation can be divided into two parts. The first and larger one provides an in-depth analysis of similarities, differences and trends among macrocycles with 20π-electron MCP's when based on different aromaticity indices (geometric, electronic, energetic, and magnetic) and functionals; how does the extent of T 1 state Baird aromaticity change when more M units are incorporated into the 20π-electron path? The second part focuses on how the knowledge gained from the small macrocycles with 20π-electron paths can be utilized to  20 and a monocyclic unit that participates with the same number of π-electrons in macrocyclic conjugation. Rings with X = O and CH 2 are named CFU and CCPD, respectively, and rings with X−X = CH=CH and CH 2 CH 2 are labeled as CPP and CCHD, respectively.
The Journal of Physical Chemistry A pubs.acs.org/JPCA Article understand larger macrocycles with potential T 1 state Baird aromaticity. Yet, to assess which functional is preferred for the T 1 states of the macrocycles investigated we first carried out single-point energy calculations with canonical UCCSD(T)/6-311G(d,p), DLPNO-UCCSD(T)/cc-pVTZ 5 5 , 5 6 and CASPT2/ANO-RCC-VDZP using the T 1 geometries at UB3LYP, 47,48 UM06-2X, 49 and UCAM-B3LYP 50 levels. Still, as we strive for conclusions that are functional independent, the discussion is often based on results from more than one functional as this provides a broader base for conclusions on general trends and differences among the macrocycles. Two idealized Baird-aromatic conformational types can be drawn for the C 20−4n [n]CM molecules in their T 1 states; one with the M units lying in the plane of the macrocycle leading to an out-of-plane π-conjugated circuit and one with the monocycles arranged perpendicularly to the macrocycle plane so that the π-conjugated circuit runs in-plane. In the T 1 state, the macrocycles with M = PP and CHD have the M units arranged perpendicularly to the plane of the macrocycle for each n, while those with M = FU and CPD adopt essentially planar conformations, even though 3 [5]CFU has a slightly bowl-shaped structure. Noteworthy, the compounds adopt the same conformational type in the S 0 and T 1 states (see Supporting Information).
Assessment of Computational Methods. Studies on [4n + 2]annulenes and expanded porphyrins in their S 0 states have shown that functionals with a medium or low amount of HF exchange, such as B3LYP (20%), tend to overstabilize bond-length equalized structures compared to coupled cluster calculations. 76−80 Peeks, Anderson, and co-workers analyzed the effect of different density functionals on the triplet delocalization of porphyrin nanorings and found that the delocalization is sensitive to the computational method; functionals that normally overdelocalize gave larger NICS(0) values.
For expanded porphyrins in the S 0 state, functionals with higher amount of HF exchange, in particular M06-2X (54%) and CAM-B3LYP (19−65%), are in better agreement with geometries from canonical CCSD(T) and DLPNO−CCSD-(T), and the extent of aromaticity as measured by electronic indices is also lower with these functionals. 76 As these results do not necessarily generalize to the T 1 state, we performed calculations of 3 C 8 H 8 , 3 C 12 H 12 , and 3 C 16 H 16 and also on 3 C 1 6 [1]CPP, 3 C 12 [2]CPP, 3 [5]CPP, 3 C 16 [1]CFU, 3 C 12 [2]CFU, and 3 [5]CFU. Each of these species adopted just one stable conformer with each of the three functionals, also for the two [5]CM's. We applied both canonical UCCSD(T) and the DLPNO-UCCSD(T). Due to large spin contamination of some compounds in T 1 with UHF, we opted to use both a restricted open-shell HF (ROHF) reference and a UBLYP DFT reference wave function for the UCCSD(T) calculation. These two references give no (ROHF) or only minute (UBLYP) spin contamination. The T 1 diagnostics values were in the range 0.011−0.031 (Table S4), and varying with the geometry used (UB3LYP, UM06-2X, or UCAM-B3LYP), they increased in the order UB3LYP < UM06-2X < UCAM-B3LYP. The values were also consistently lower in CCSD calculations with the ROHF reference wave function. The fractional occupation number weighted density (FOD) 62 values reveal that the multiconfigurational character in the T 1 state is highest in cyclic C 20 (Tables S5 -S6), and that it decreased with number of M units in 3 C 20−4n [n]CM. For this reason, we disregard cyclic C 20 hereafter. Finally, for 3 C 16 [1]CPP, 3 C 16 [1]CFU and the three [4n]annulenes in T 1 , we applied CASPT2//UDFT.
For the [4n]annulenes in their T 1 states, the lowest coupled cluster energies are generally found when the calculations are carried out on UB3LYP optimized geometries, followed by UM06-2X optimized geometries. This trend is found at both canonical UCCSD(T) and DLPNO-UCCSD(T) levels (Tables 1 and 2), and it is supported by results at CASPT2/ ANO-RCC-VDZP//UDFT/6-311G(d,p) level (Table 1). However, the UB3LYP geometry of C 16 H 16 gives the lowest UCCSD(T) energy with an ROHF initial guess while the UM06-2X geometry gives the lowest energy with a UBLYP guess. For 3 C 16 [1]CFU and 3 C 12 [2]CFU, UB3LYP geometries lead to the lowest coupled cluster energies regardless of initial guess, while for 3 C 16 [1]CPP and 3 C 12 [2]CPP it depends on the choice of guess (ROHF or UBLYP). For 3 [5]CPP and 3 [5]CFU, DLPNO-UCCSD(T) calculations still provisionally indicate that the UB3LYP geometry is the best, but for these species, the calculations had to be done with a small basis set. A caveat here is the observation that only one conformer each exist for the two [5]CM's, also with UCAM-B3LYP. For the larger [n]CPP's, several conformers can exist and our conclusions on method preference are therefore not fully unambiguous as we were unable to carry out DLPNO-UCCSD(T) calculations on these.
Noteworthy, DLPNO−CCSD(T) calculations of systems with multireference character should be carried out with the TightPNO keyword. 81 Such calculations for 3 C 16 [1]CFU and 3 C 16 [1]CPP, the two species with the highest T 1 diagnostics and FOD values (Table S6), give relative energies which differ by up to 2.2 kcal/mol from the default DLPNO−CCSD(T) results, yet, always in the direction favoring coupled cluster energies based on the B3LYP geometries. This finding is corroborated by results on the two compounds at CASPT2// UDFT level as the UB3LYP T 1 geometries again lead to the lowest energies (Table S7). Thus, UB3LYP on average performs best, followed by UM06-2X, while UCAM-B3LYP performs worst. Yet, as indicated above we base our conclusions in the following on at least two functionals despite the provisional preference for UB3LYP.
The differences between the T 1 geometries of the compounds at UB3LYP and UM06-2X levels are significant as the CC bond lengths differ by up to 0.081 Å and on average by 0.015 Å (see Figure S2 and S5). They differ on average even more between UB3LYP and UCAM-B3LYP (0.021 Å), although the maximal difference is smaller (0.051 Å). For Table 1. Relative Energies at Canonical UCCSD(T)/6-311G(d,p)//UDFT/6-311G(d,p) (Normal Print) and CASPT2/ANO-RCC-VDZP//UDFT/6-311G(d,p) (Italics) Level with UB3LYP, UM06-2X, and UCAM-B3LYP  (Tables S3 and S6), in line with the observation by Matito and co-workers. 82 Here it can be noted that spin delocalization has earlier been observed experimentally in porphyrin nanorings by ENDOR spectroscopy, 83 in line with the UB3LYP result. Yet, the time scale of the EPR measurements (up to 2 μs) may not capture the situation with a rapid interconversion between conformers with semilocalized spin densities. Another aspect of relevance for the computational method choice is the presence of repulsive H−H interactions between adjacent M units and the potential need for dispersioncorrected DFT. However, through a comparison of geometries calculated with and without Grimme's dispersion correction, 84 we find that the effects are only modest and that they do not alter the results significantly (see Supporting Information).
Geometrical Impacts of Macrocyclic Baird Aromaticity. A high degree of aromaticity normally goes with highly symmetric structures. Thus, one can expect that 3 [5]CM macrocycles with marked Baird aromaticity in their T 1 states possess D nh or C nv symmetry (or nearly so), whereby each monocycle in the macrocycle is structurally equivalent. A Baird-aromatic character should also result in bond length equalization along the CC bonds in the perimeter. To probe this we applied the geometry-based HOMA index, 70 although it should be noted that it can only be applied for the perimeters in the full [n]CM's as it will severely underestimate the T 1 aromatic character of species with polyynic segments. Additionally, one monocycle (furan) is known to be a difficult case for HOMA. 85 For these reasons, we also estimated the impact of T 1 state Baird aromaticity from geometrical parameters that function as approximate indicators.
The HOMA values for the all-carbon perimeters in the [5]CM's in their S 0 states are low, corresponding to nonaromatic situations (Table 3). Yet, they increase in the T 1 states with each of the three functionals with the largest increase (ΔHOMA(T 1 −S 0 )) in [5]CFU leading to HOMA values (0.75−0.90) that unambiguously correspond to macrocyclic Baird aromaticity. Noteworthy, despite that 3 [5]CFU with UM06-2X and UCAM-B3LYP lacks symmetry it is still highly aromatic according to HOMA. The increases in the perimetric HOMA values upon excitation to T 1 are also significant in [5]CPP, but for this compound the actual values in T 1 are significantly lower (0.44−0.53) than for 3 [5]CFU, and correspond to weaker Baird aromaticity or nonaromaticity. In contrast, the ΔHOMA(T 1 −S 0 ) for [5]CCHD and [5]CCPD are more modest, although with UB3LYP they both have higher HOMA values in their T 1 states than The geometrical impact of T 1 state Baird aromaticity is also visible in several other geometric parameters, particularly if one considers the geometric changes when going from S 0 to T 1 . One such parameter, the C−C−C−C dihedral angle between adjacent monomer units (Λ(C−C−C−C)) in 3 [5]CM's, should decrease upon excitation for systems that adopt Baird aromaticity. Indeed, when going from S 0 to T 1 the average Λ(C−C−C−C) angles of 3 [5]CPP, 3 [5]CCPD, and 3 [5]CCHD decrease at all three DFT levels, indicative of attainment of some macrocyclic conjugation in the T 1 state. However, for [5]CFU, this angle is already small in S 0 , and there is therefore either no change or only a slight increase upon excitation to T 1 . Another useful geometrical parameter is Relative energies obtained when using the TightPNO keyword in the DLPNO-UCCSD(T) calculations. b Calculations with the cc-pVDZ basis set. Table 3. Perimetric HOMA Values, Λ(C−C−C−C) Dihedral Angle Ranges, and Average Differences in Λ(C−C−C−C) between the S 0 and T 1 States at (U)B3LYP (Normal Print), (U)M06-2X (italics), and (U)CAM-B3LYP (Parentheses) Levels The Journal of Physical Chemistry A pubs.acs.org/JPCA Article the length of the C−C bonds between the monocycles as these bonds should shorten in Baird-aromatic macrocycles. The shortening in these bonds is apparent in both 3 [5]CPP and 3 [5]CFU with all three functionals (  [1]CM to 3 C 12 [2]CM at UB3LYP level, one can note a significant weakening of the T 1 Baird-aromatic character of Electronic Impact of Macrocyclic Baird Aromaticity. In an idealized (strongly) Baird-aromatic cycle with a D nh or C nv symmetric structure, the spin density distribution should be such that there is a continuous excess of α-spin (or β-spin) density along the complete cycle. Excess of β-spin density at any atom(s) along the path may suggest a discontinuous (acyclic) distribution of the two unpaired α-electrons of the triplet diradical, and consequently, a lowered or negligible Baird aromaticity. At UB3LYP level, this spin density criterion is satisfied by eight of the 12 20-π-electron macrocyclic compounds examined herein; the four exceptions are 3 [5]CCPD and each of the three 3 C 20−4n [n]CCHD's (n = 1, 2, or 5). Two of these compounds ( 3 [5]CCPD and 3 [5]CCHD) have geometries where the cyclic π-conjugation in the MCP is reduced or broken due to steric congestion between adjacent M units while the other two have planar structures that nominally could provide strong Baird aromaticity, yet they are nonetheless nonaromatic. However, the functional has a decisive impact on the spin density distribution as seen for 3 [5]CFU since the excess in α-spin density gets interrupted by β-spin excess at points along the macrocycle when going from UB3LYP to UM06-2X and UCAM-B3LYP (Figure 3), and this also occurs for 3 [5]CPP ( Figure S9). Thus, a continuous and evenly distributed α-spin excess around a complete cycle reveals idealized Baird aromaticity, although species with modest disruptions are still Baird-aromatic as seen for 3 [5]CFU where the HOMA value is well above 0.5 with UCAM-B3LYP.
The extent of Baird aromaticity was also assessed through the electronic FLU index 71     Energy Impacts of Macrocyclic Baird Aromaticity. As an energy-based aromaticity index we applied isomerization stabilization energies (ISEs) in versions where the 1,3hydrogen shift(s) takes place at one single (di)methyl substituted M unit. The methyl/methylene groups are placed so that steric congestion is avoided on either side of the ISE reaction ( Figure 4). For the compounds with PP monocycles we used dimethyl substitution at one PP unit as the macrocyclic conjugation otherwise cannot be broken.
For each M, and with both UB3LYP and UM06-2X, one can note gradually decreasing ISE values when n increases in 3 C 20−4n [n]CM (n = 1, 2 or 5), indicating reduced Baird aromaticity. For 3 [5]CPP and 3 [5]CFU the ISE values are still negative, revealing stabilization upon formation of Bairdaromatic macrocycles in T 1 , but they are positive for 3 [5]CCHD and 3 [5]CCPD indicating a lack of aromatic stabilization. Also, the ISE values for 3 C 20−4n [n]CFU's are always higher than those of the analogous 3 C 20−4n [n]CPP. As a comparison, the ISE values of 3 C 8 H 8 and 3 C 20 H 20 at (U)B3LYP ((U)M06-2X) levels are −16.9 (−10.8) and −19.1 (−14.4) kcal/mol, and for benzene in the S 0 state (the aromatic archetype), it is −33.9 (−34.2) kcal/mol, resembling earlier computed ISE values for small aromatic monocycles in T 1 and S 0 . 72,86,87 The T 1 state aromatic stabilization in a polycyclic compound with a cyclooctatetraene core has also been determined experimentally as 21−22 kcal/ mol. 14 Thus, the Baird-aromatic stabilization of 3 [5]CPP is already low with UB3LYP, and with UM06-2X it is merely 23% of the ISE value of S 0 state benzene at the same level. For 3 [5]CFU, on the other hand, the ISE values are essentially equal to those of triplet state C 20 H 20 and approximately half the aromatic stabilization of benzene in S 0 . Hence, it is apparent at both UB3LYP and UM06-2X levels that the aromatic stabilization drops as the number of PP and FU units increase in a 3 C 20−4n [n]CM. Yet, the results also suggest that a macrocycle 3 [n]CM composed of moderately Huckel-aromatic M units such as furans has a higher likeliness to sustain macrocyclic Baird aromaticity to larger sizes in its T 1 state. With UCAM-B3LYP the ISE values of 3 [5]CPP and 3 [5]CFU are −7.8 and −12.5 kcal/mol, respectively, resembling those at UM06-2X level ( Table 6).

Magnetic Aspects of Macrocyclic Baird Aromaticity.
Magnetic aromaticity indices (NICS and induced current density maps) are the most frequently applied computational tools for aromaticity assessments. 88 However, it has been pointed out that these indicators should not be used alone, 46 and it has been shown that NICS can be an ambiguous measure of magnetically induced paratropic or diatropic ring currents in molecules. 89 NICS values should therefore be combined with current density maps. Now, to what extent do the results of magnetic indices on potential Baird aromaticity in the present macrocycles differ from the results of other types of indices?
With regard to NICS, all 3 C 20−4n [n]CM (n = 1, 2 or 5) with M = PP or FU calculated with UB3LYP or UM06-2X show    [5]CPP than in 3 [5]CFU. This ordering in relative degree of aromaticity is contrary to the findings made with the other index types. It is also noteworthy that the strongest aromatic character according to NICS(0) zz is found for the 3 C 12 [2]CM species. There is a small attenuation when going to 3 [5]CPP while 3 [5]CFU is as aromatic as  20 there is a distinct variation between different functionals to what extent they describe large (macro)cycles as Baird-aromatic. CAM-B3LYP strongly disfavors Baird aromaticity, yet based on this feature, one may also single out which macrocycles with high certainty are Baird-aromatic in reality. On the other hand, our coupled cluster computations indicate that the UB3LYP geometries are the best for the species studied herein. With regard to the aromaticity indices, they report essentially similar trends between the various 3 C 20−4n [n]CM compounds, except for NICS indicating that the Baird aromaticity in 3 C 20−4n [n]CPP's to be more markedly Baird-aromatic than the Baird-aromatic 3 C 20−4n [n]CFU's.
Finally, with regard to the first aspect listed above, 3 [5]CFU has a highly Baird-aromatic character according to all aromaticity indices applied herein and according to two of the three functionals applied. Furthermore, the attenuation in Baird-aromatic character when going from 3 C 16 -[1]CM to 3 [5]CM is smallest when M = FU indicating that this unit is suitable for designing larger [n]CM macrocycles. In contrast, the [n]CPP's seem not to be optimal if one seeks macrocycles with strong T 1 state Baird aromaticity, unless further tailored.
T 1 State Baird Aromaticity of Larger and Altered [n]CMs. We now focus on larger 3 [n]CPP and 3 [n]CFU (n > 5), and their derivatives, in order to assess at what sizes n these macrocycles cease to be Baird-aromatic. The 3 [n]CFU's are particularly promising as 3 [5]CFU was Baird-aromatic even with UCAM-B3LYP, in contrast to 3 [5]CPP. Conversely, since 3 [5]CCPD and 3 [5]CCHD are nonaromatic even with UB3LYP, larger 3 [n]CCPD's and 3 [n]CCHD's with n > 5 are not further considered. For macrocycles for which we found Baird-aromatic minima we have checked for nonaromatic minima by starting from several different geometries that were markedly distorted.
According to UCAM-B3LYP, the 3 [n]CPP's with n > 5 have egg-shaped structures which are nonaromatic, 25 and similar attenuations of the Baird aromaticity occur also with UB3LYP   Table S14.
The Journal of Physical Chemistry A pubs.acs.org/JPCA Article and UM06-2X but at higher n. The distortion to nonaromatic 3 [n]CPP's is complete at n = 8 with UB3LYP and at n = 6 or 7 with UM06-2X, seen both with NICS and HOMA (Tables 8, S17, and S18). Yet the FLU values, in contrast to NICS and HOMA, reveal no significant weakening of the electron delocalization in the 3 [n]CPP's when going from n = 6 to n = 8; neither with B3LYP nor M06-2X. That is, also the unsymmetric 3 [8]CPP has FLU values that correspond to Baird aromaticity with both functionals despite that the interring ΛC−C−C-C angles are found in the ranges 2.2−29.2°( UB3LYP) and 6.9−35.7°(UM06-2X), and with several dihedral angles larger than 20°(four angles with UB3LYP and six angles with UM06-2X). Additionally, the delocalization index (DI), on which FLU depends, also does not help to explain the incongruencies between this electronic index, on the one side, and NICS and HOMA on the other. By taking the mean of the DI of all bonds that participate in the circuit, no significant variation is observed when going to larger macrocycles (the mean DI = 1.296, 1.297, and 1.299 e for 3 [6]CPP, 3 [7]CPP, and 3 [8]CPP, respectively).
For the larger macrocycles, the energy differences between the highly symmetric (Baird-aromatic) T 1 state structures and the distorted (nonaromatic) ones vary extensively with functional. With UB3LYP, the C 4v and D 8h symmetric structures of 3 [8]CPP, which are second-and higher-order saddle-points, are merely 2.1 and 3.7 kcal/mol higher in energy than the unsymmetric one. In contrast, with UM06-2X and UCAM-B3LYP these structures are 22−30 kcal/mol higher in energy. In the singlet excited states Kim, Nakano, Yamago, and co-workers recently revealed rapid interconversion between self-trapped exciton and delocalized states, 90 yet it has earlier been found through time-dependent CAM-B3LYP calculations that the S 1 state of [n]CPP's is more delocalized than the T 1 state. 25 Similar observations in the T 1 state of various [n]CPP's could give experimental support as to which functional gives the right description: UB3LYP or UCAM-B3LYP?
Now, how to design large [n]CM's (or derivatives) that are Baird-aromatic in their T 1 states even with UCAM-B3LYP? To enhance the Baird-aromatic character in 3 [n]CPP's it should be crucial to reduce the steric congestion between adjacent PP units when aligned. Itami and co-workers recently reported synthesis and further experimental studies of methylenebridged [6]CPP (MB[6]CPP), 91 a compound that displays magnetically induced paratropic ring currents in the S 0 state. 92 We now explored 3 MB[n]CPP's with n = 6−8 using UB3LYP, UM06-2X, and UCAM-B3LYP, and also 3 MB [12]CPP with UCAM-B3LYP. Indeed, with each of the three functionals the 3 MB[n]CPP's exhibit fully delocalized yet unevenly distributed α-spin density (see Table S9). The HOMA and NICS(0) zz values of 3 MB[8]CPP with UCAM-B3LYP (0.78 and −40.3 ppm, respectively) clearly reflect Baird-aromatic character. The larger 3 MB [12]CPP calculated with UCAM-B3LYP, however, goes to a structure which similar to 3 [12]CPP has the spin density distribution localized to a part of the macrocycle and Huckel-aromatic benzene rings in the other part.
The Baird-aromatic character also impacts on the adiabatic T 1 energies (E(T 1 )) as seen in Table 9. For the Baird-aromatic MB[n]CPP's one can note a gradual increase in E(T 1 ) when going from n = 6 to n = 8 with each of the three functionals. In contrast, the E(T 1 ) for the unconstrained [n]CPP's, calculated with UM06-2X and UCAM-B3LYP, drop below the ones of the corresponding MB[n]CPP's, particularly for n = 8. This indicates that the geometrical relaxation available in the latter, leading to semilocalized triplet diradical character and a few PP units with closed-shell Huckel-aromatic character, providing a stabilization of the T 1 state with the loss of Baird-aromatic character.
To tailor large unconstrained macrocycles which are Bairdaromatic in T 1 one may argue that an approach is to utilize M units that can retain Huckel aromaticity in one part of the unit simultaneously as another part is involved in the Bairdaromatic macrocyclic circuit. The 1,4-linked naphthalene could be such a unit, but at UCAM-B3LYP level the [n]cycloparanaphthalenes ([n]CN's) loose Baird-aromatic character already at n = 6. Steric congestion is likely a contributing factor as the repulsive nonbonded H−H interactions should increase in [n]CN's when compared to that in [n]CPP's. Additionally, a Baird-aromatic 3 [6]CN can be described as a triplet state [24]annulene, a species where the lowest energy conformer has a distorted structure with attenuated aromaticity according to UCAM-B3LYP.
Other M units than those based on strongly Huckelaromatic benzene rings are required in order to achieve unconstrained large [n]CM macrocycles with strong Baird aromaticity. For this aim, it should be favorable to utilize M units that are moderately Huckel-aromatic in S 0 and/or impose Indeed, the existence of a local Baird-aromatic minimum on the T 1 potential energy surface of 3 [8]CFU with UB3LYP is interesting as it supports a tug-of-war between a macrocyclic T 1 state Baird-aromatic conformer and one conformer having a semilocalized triplet diradical together with one or a few closed-shell Huckel-aromatic FU units. This feature impacts on the E(T 1 ) energies because, similar as for the [n]CPP's, one can note that T 1 state Baird-aromatic [n]CFU's have higher E(T 1 ) than the ones with nonaromatic global minima in T 1 ( Table 9). This feature also points to a pitfall in the explorations of a potentially Baird-aromatic macrocycle as the Baird-aromatic minimum located may not be the global minimum on the T 1 PES.
Based on the observation that weakened S 0 state Huckel aromaticity together with reduced steric congestion brings increased T 1 state Baird aromaticity one may ask if changing the PP units in [n]CPP's to either of the three diazines (pyridazine, pyrimidine, and pyrazine) can provide for hetero-[n]CPP's which are Baird-aromatic at larger n's than the allcarbon [n]CPP's? These heterocycles have a lower Huckel aromaticity in S 0 , 26 and they are less sterically demanding in nanohoops than the PP unit. 93 Hence, we explored the hetero- [n]CPP's with UCAM-B3LYP as this functional should underestimate the Baird-aromatic character. Yet, we find that each of the three different hetero-[n]CPP's, n = 6 and 8, have unsymmetric global minima in T 1 and semilocalized spin densities indicating nonaromatic character. Several of the macrocycles also have Baird-aromatic local minima 12−26 kcal/mol higher in energy. Interestingly, some of these higherenergy local minima vanish with UB3LYP and UM06-2X so that only the nonaromatic structures exist with those functionals.
Neutral Macrocycles in T 1 vs Dicationic and Dianionic Ones in S 0 . Interesting perspectives on macrocycles with potentially Baird-aromatic T 1 states can also be gained through comparisons with the corresponding dications and dianions. Several aromatic [n]CPP 2+ have earlier been reported, 33,94−97 and based on the Zilberg−Haas valence bond description of triplet state aromaticity ( Figure 1B), 33 one may ask to what extent they resemble 3 [n]CPP's? On the other hand, from a molecular orbital perspective it could be more obvious to compare with the dianions as the same orbitals are occupied. Also the dianions [n]CPP 2− have been explored earlier, and they are aromatic according to computations and have significant quinoidal character in the X-ray crystal structures. 34−36 Indeed, the dications, dianions, and neutral triplets when viewed with Huckel MO theory all have the same resonance energy as an equal occupancy of the orbitals which are HOMO and LUMO in the neutral S 0 state will cancel out the contributions from these orbitals according to the Pairing theorem for even alternant hydrocarbons. Thus, the dianion and the neutral triplet should both resemble the dication in terms of resonance stabilization and aromaticity.
The dication [5]CPP 2+ in its closed-shell S 0 state has a similar NICS value as 3 [5]CPP (with (U)B3LYP: −41.8 and −40.5 ppm, respectively), and the similarity is underlined by nearly identical HOMA values (0.571 and 0.572, respectively). The dianion [5]CPP 2− has NICS and HOMA values of −38.7 ppm and 0.523, respectively. Yet, the analogy between the three species disappears when going to gradually larger   [7]CPP. For the latter, the ISE values are much smaller than those of [7]CPP 2+ , also with UB3LYP, revealing that the macrocyclic Baird aromaticity in the neutral triplet is negligible.
According to both NICS and HOMA at B3LYP and M06-2X levels, the aromatic character of [n]CFU 2+ and [n]CPP 2+ diminishes at, respectively, n > 10 and 12 (Tables S17 and S18), and similar observations are made for the dianions. Indeed, the observation that dicationic macrocycles provide for stronger aromaticity than the neutral T 1 state is in line with findings made by the Anderson group on the much larger porphyrin nanorings. 98 For these, magnetically induced diatropic ring-currents were most pronounced in oxidized nanorings having a charge per porphyrin unit in the range 0.5− 0.7. In another finding, which similar to that for 3 [n]CFU and 3 [n]CPP, a [n]CFU 2+ macrocycle is slightly more aromatic than a [n]CPP 2+ with the same n according to HOMA, ISE, and FLU while it is the opposite between the two species with NICS ( Figure 6 and Tables S17 and S18).
Taken together, the macrocyclic dications [n]CM 2+ and dianions [n]CM 2− provide upper limits in n for Baird-aromatic character in the T 1 state of a [n]CM. If the macrocyclic dication or dianion, which potentially is readily investigated experimentally, is not Huckel-aromatic, then the analogous neutral triplet will not be Baird-aromatic.

■ SUMMARY AND OUTLOOK
Herein we used computational tools to explore a series of macrocycles composed of either benzene, furan, cyclohexa-1,3diene or cyclopentadiene units linked together so that macrocyclic conjugated 4nπ-electron paths are enabled. These macrocycles may all exhibit T 1 state Baird-aromatic character, yet, a number of aspects limit the attainment of such character. First, steric congestion distorting the structures away from the ones with maximal macrocyclic π-conjugation is highly detrimental for T 1 state Baird aromaticity. A second factor is the tug-of-war between macrocyclic Baird aromaticity, on the one side, and closed-shell Huckel aromaticity of one or a few M units together with semilocalized tripled diradical character on the remaining M units, on the other. The 2,5linked furan unit (FU) is ideal as it does not induce extensive steric congestion, and it has a lower Huckel-aromatic character in S 0 than benzene. Hence, regardless of functional used we observe that the T 1 state Baird aromaticity is less rapidly attenuated in 3 [n]CFU's than in 3 [n]CPP's. On the other side, the steric congestion between adjacent PP units in [n]CPP's can be overcome by methylene bridges between the PP units.
The functional used has extensive impact on the breakpoint n at which a [n]CM macrocycle class goes from T 1 state Bairdaromatic to nonaromatic character, leading to different breakpoints n at which the T 1 structures switch character. Our comparative computations with coupled cluster lead to a provisional preference for UB3LYP geometries, but the trends are similar to the three functionals used. Other functionals may provide different threshold values n (e.g., GGA functionals at somewhat higher n than those here), yet the trend will remain. Here, it should also be stressed that B3LYP is known to overestimate aromaticity of large π-conjugated cycles in the S 0 state, 82 and one should be careful when assessing the degree of triplet state aromaticity of a macrocycle. There is also a variation between different aromaticity indices, although the overall trend is that described above, except for NICS which describes the few 3 [n]CPP's that are Baird-aromatic to have stronger such character than the analogous 3 [n]CFU's. Combined, as the choice of functional and aromaticity index, both impact on whether a certain macrocycle is described as triplet state Baird-aromatic or not; one should be cautious in the assessment. As the evaluation of Baird-aromatic character relies heavily on computational results, one should preferably speak of likeliness that a certain macrocycle is Baird-aromatic (or not). There are certainly pitfalls in the area of excited state aromaticity and antiaromaticity.
We also explored the analogy between the neutral triplets and the macrocycle dications, following the valence bond theoretical model by Zilberg and Haas, 31 and find that the neutral triplets turn nonaromatic at lower n than the dications. Comparisons were also made with the dianions, which in terms of aromaticity resemble the corresponding dications. Hence, if a macrocyclic dication (or dianion) is nonaromatic then the neutral triplet will be as well.
To extrapolate from our findings, we postulate that the tugof-war can tip to the favor of macrocyclic Baird aromaticity if π-conjugated linkers are inserted between the M units. This also eliminates steric congestion between adjacent M units. A few such molecules have been prepared and explored experimentally. 99,100 The structurally more complex porphyrinoid macrocycles could also be interesting as the T 1 state Baird aromaticity in such compounds has been explored extensively. 7,13,101−107 For the S 0 state, Aihara and co-workers concluded, based on a graph theoretical analysis, that Huckel's rule is applicable to the main conjugation paths of porphyrinoids having 4n + 2 π-electrons. 108−114 Yet, the macrocyclic aromaticity in S 0 was significantly weaker than the aromaticity of the individual monocycles. How does this change in the T 1 state? These aspects of macrocyclic T 1 Baird aromaticity will be reported in forthcoming studies.