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, poly<span class="Chemical">cyclic, 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 Hückel-aromatic in their electronic ground states (S0) while the other two are nonaromatic with 1,3-diene segments. Our focus is on the potential T1 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 macro<span class="Chemical">cyclic conjugation pathway (MCP), including both Hückel-antiaromatic ([4n]) and -aromatic ([4n + 2]) cases. These nanorings are exceptional as they are composed of a series of Hückel-aromatic porphyrins linked via butadiynyl segments. Yet, it was observed that the extent of T1 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 T1 states it has been considered that long-range corrected functionals are preferred.[25] When going to gradually larger [n]CPP’s in their T1 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 T1 state there is a tug-of-war between dominance of macrocyclic Baird-aromatic character of <span class="Chemical">cyclic 4nπ-electron MCP’s and dominance of closed-shell Hückel-aromatic character of one or a few monocycles which forces the triplet diradical character to localize to a part of the macrocycle (Figure A). 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 Hückel 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 S0 state geometries and have electronic properties which make them interesting for organic electronics.[28] The effect of heteroatoms on the T1 state Baird aromaticity of [6]CPP has been reported in a computational study where the phenylene units were changed successively to 1,4-azaborines,[6]Cycloparaphenylenes. J. Mater. Chem. C. 2017 ">29] a monocycle which is less Hückel-aromatic in S0 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?

Figure 1

(A) Two (resonance) structures of [9]CPP, one representing macrocyclic Baird aromaticity and the other with a locally Hückel-aromatic PP unit and with the triplet diradical distributed on the other eight PP units. (B) Conceptual formation of [9]CPP from [9]CPP by addition of two nonbonding same-spin π-electrons to [9]CPP, which has a Hückel-aromatic MCP with 4 × 9–2 (i.e., 4 × 8 + 2) spin-paired π-electrons.

An interesting description of Baird aromaticity is the description of Zilberg and Haas in which the T1 state of a neutral 4nπ-electron cycle is viewed as that of the corresponding Hückel-aromatic dication with 4n – 2 spin-paired π-electrons plus two nonbonding same-spin π-electrons.[31] This description was utilized to rationalize the aromaticity/antiaromaticity switch that occurs when going from S0 to T1 in a [26]- and [28]porphyrinoid pair.[13,32] Indeed, the potential Baird aromaticity of a neutral [n]CPP in its T1 state could therefore closely resemble that of the corresponding [n]CPP dication (Figure B), species for which in-plane Hückel aromaticity with 4n − 2 π-electrons has been demonstrated.[33] Yet, is the aromatic character stronger or weaker in the T1 state of [n]CPP than in the closed-shell singlet state [n]CPP, or is it similar? Furthermore, is there a similarity to the dianions which are Hückel-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 C16, 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 (Qu1),[37] in contrast to the S0 state where they are doubly Hückel-antiaromatic until a certain size.[38−40] For the T1 state, it was concluded that the out-of-plane π-conjugated path of C16 displayed Baird aromaticity while the in-plane path displayed Hückel antiaromaticity. A situation in the T1 state with two similarly strong counteracting ring-currents, one diatropic (Baird-aromatic) and one paratropic (Hückel-antiaromatic), could lead to a NICS value that corresponds to a nonaromatic situation. Indeed, such seemingly nonaromatic characters were recently observed from computations of cyclic C16, C20. and C24.[41] Furthermore, after the recent experimental generation of cyclo[18]carbon by Anderson and co-workers,[18]carbon. Science. 2019 ">42,18]carbon via Debromination of C. J. Am. Chem. Soc.. 2020 ">43] the aromaticity of cyclo[4n + 2]carbons has gained increased attention,[39,44] yet these are outside our explicit focus herein.

Now, if every 1,3-butadiyne unit in a <span class="Chemical">cyclo[4n]carbon is replaced by monocycles contributing four π-electrons to the macrocyclic MCP (Figure ) 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-of-plane orbitals in [n]CFU’s with n < 8 (cf., trannulenes and annulenes, respectively).

Figure 2

1,3-Butadiyne segment in cyclic C20 and a monocyclic unit that participates with the same number of π-electrons in macrocyclic conjugation. Rings with X = O and CH2 are named CFU and CCPD, respectively, and rings with X–X = CH=CH and CH2CH2 are labeled as CPP and CCHD, respectively.

1,3-Butadiyne segment in <span class="Chemical">cyclic C20 and a monocyclic unit that participates with the same number of π-electrons in macrocyclic conjugation. Rings with X = O and CH2 are named CFU and CCPD, respectively, and rings with X–X = CH=CH and CH2CH2 are labeled as CPP and CCHD, respectively.

As [5]CPP is the smallest [n]CPP that has been synthesized and explored experimentally,[45] we start our analysis at the C20 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 <span class="Chemical">benzenoid structure in its S0 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 ); 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[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 S0 state dications (e.g., [n]CPP and [n]CFu) and the corresponding neutral triplet states? One may also consider the similarity to the corresponding S0 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.

Computational Methods

The S0 and T1 geometries were optimized using the B3LYP,[47,48] M06-2X,[49] and CAM-B3<span class="Gene">LYP[50] functionals and the 6-311G(d,p) basis set.[51] Dispersion-corrected UB3LYP with Becke–Johnson damping (GD3BJ)[52] was also computed in some cases. Relative energies were calculated using UDFT (B3LYP, M06-2X, and CAM-B3LYP) and coupled cluster with single and double substitutions and triple excitations with single point calculations at the various triplet UDFT geometries (UCCSD(T)/6-311G(d,p)//UDFT/6-311G(d,p)).[53,54] Canonical CCSD(T) was only computed for C8H8, C12H12 and C16H16. Domain pair natural orbitals (DLPNO) were employed in most of the CCSD(T) calculations (DLPNO-UCCSD(T)/cc-pVTZ//UDFT/6-311G(d,p)).[55−57] In the DLPNO-UCCSDT computations, cc-pVTZ was chosen as an auxiliary basis set. In the canonical UCCSD(T) and DLPNO-UCCSD(T) calculations, three initial guesses were used, i.e., UHF, BLYP and ROHF. For C8H8 and C12H12, CASPT2/ANO-RCC-VDZP//UDFT/6-311G(d,p) calculations were carried out using the atomic natural orbital (ANO-RCC-VDZP) basis set, while for C16H16 CASPT2/6-31G(d)//UDFT/6-311G(d,p) was instead used due to the large active space used.[58−60] Regarding the latter computations, the active space used was 8 electrons in 8 orbitals (8in8), 12in12 and 16in16 for C8H8, C12H12, and C16H16, respectively. Most of the calculations (including the T1 diagnostic using the “T1diag” keyword) were performed with Gaussian 16 revision B.01[61] except for DLPNO-UCCSD(T), fractional occupation number weighted density (FOD)[62] and CASPT2 computations for which ORCA[63] and OpenMolcas packages were used,[60] respectively. Aromaticity was evaluated in terms of the nucleus independent chemical shift (NICS)[64,65] with the GIAO method,[66] ACID plots,[67,68] the harmonic oscillator model of aromaticity (HOMA),[69,70] the aromatic fluctuation index (FLU)[71] and isomerization stabilization energies (ISEs),[72] computed at the B3LYP and M06-2X optimized geometries. ACID plots were computed using the software by Herges,[67] HOMA was computed using Multiwfn[73] and FLU was obtained with the ESI-3D program[74] and AIMAll.[75]

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 T1 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 understand larger macrocycles with potential T1 state Baird aromaticity. Yet, to assess which functional is preferred for the T1 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[55,56] and CASPT2/ANO-RCC-VDZP using the T1 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[n]CM molecules in their T1 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 T1 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 [5]CFU has a slightly bowl-shaped structure. Noteworthy, the compounds adopt the same conformational type in the S0 and T1 states (see Supporting Information).

Assessment of Computational Methods

Studies on [4n + 2]annulenes and expanded <span class="Chemical">porphyrins in their S0 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 S0 state, functionals with higher amount of HF exchange, in particular M06-2X (54%) and CAM-B3<span class="Gene">LYP (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 T1 state, we performed calculations of 3C8H8, 3C12H12, and 3C16H16 and also on C[1]CPP, C[2]CPP, [5]CPP, C[1]CFU, C[2]CFU, and [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 T1 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 T1 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 T1 state is highest in cyclic C20 (Tables S5 - S6), and that it decreased with number of M units in C[n]CM. For this reason, we disregard cyclic C20 hereafter. Finally, for C[1]CPP, C[1]CFU and the three [4n]annulenes in T1, we applied CASPT2//UDFT.

For the [4n]annulenes in their T1 states, the lowest coupled cluster energies are generally found when the calculations are carried out on UB3<span class="Gene">LYP optimized geometries, followed by UM06-2X optimized geometries. This trend is found at both canonical UCCSD(T) and DLPNO-UCCSD(T) levels (Tables and 2), and it is supported by results at CASPT2/ANO-RCC-VDZP//UDFT/6-311G(d,p) level (Table ). However, the UB3LYP geometry of C16H16 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 C[1]CFU and C[2]CFU, UB3LYP geometries lead to the lowest coupled cluster energies regardless of initial guess, while for C[1]CPP and C[2]CPP it depends on the choice of guess (ROHF or UBLYP). For [5]CPP and [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.

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 Geometriesa

geometries	B3LYP	CAM-B3LYP	M06-2X
3C8H8	0.0, 0.0	0.8, 0.4	0.4, 0.2
3C12H12	0.0, 0.0	(4.4)b, 1.5	(4.4)b, 0.0
3C16H16	0.0, 0.0c	(7.6)b, 2.4c	1.4, 0.9c

UHF reference wavefunction in UCCSD(T) calculation. Relative energies in kcal/mol.

Values in parentheses are results for which the spin contamination in the reference exceeds 20%, as seen in Table S5.

CASPT2 calculations for 3C16H16 run with the 6-31G(d) basis set.

Table 2

Relative DLPNO-UCCSD(T)/cc-pVTZ Energies (in kcal/mol) in T1 State at Geometries Optimized with UB3LYP, UCAM-B3LYP, and UM06-2X and with UBLYP and ROHF Reference Wavefunctions (ref. wfn.)

ref. wfn.	BLYP	ROHF
geometry	B3LYP	CAM-B3LYP	M06-2X	B3LYP	CAM-B3LYP	M06-2X
3C8H8	0.0	0.3	0.1	0.0	0.3	0.1
3C12H12	0.0	1.1	0.2	0.0	1.6	1.0
3C16H16	0.0	0.3	–0.2	0.0	4.2	3.7
3C16[1]CFU	0.0, 0.0a	7.3, 6.7a	0.5, 1.3a	0.0, 0.0a	10.2, 9.4a	0.3, 1.4a
3C12[2]CFU	0.0	2.1	0.6	0.0	1.6	5.0
3[5]CFU	0.0b	1.4b	0.8b	0.0b	3.2b	2.8b
3C16[1]CPP	0.0, 0.0a	0.6, 2.8a	–0.6, 0.0a	0.0, 0.0a	0.6, 2.5a	–0.5, 0.0a
3C12[2]CPP	0.0	–1.5	–2.3	0.0	2.4	1.4
3[5]CPP	0.0b	2.6b	0.7b	0.0b	6.0b	2.7b

Relative energies obtained when using the TightPNO keyword in the DLPNO-UCCSD(T) calculations.

Calculations with the cc-pVDZ basis set.

Noteworthy, DLPNO–CCSD(T) calculations of systems with multireference character should be carried out with the TightPNO keyword.[81] Such calculations for C[1]CFU and C[1]CPP, the two species with the highest T1 diagnostics and FOD values (Table S6), give relative energies which differ by up to 2.2 kcal/mol from the default DLPNO–<span class="Chemical">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 T1 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 T1 geometries of the compounds at UB3LYP and <span class="Chemical">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 the [4n]annulenes in their T1 states, the bond length alternations is smallest with UB3LYP and largest with 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 dispersion-corrected 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 [5]CM macrocycles with marked Baird aromaticity in their T1 states possess D or C 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 T1 aromatic character of species with polyynic segments. Additionally, one monocycle (<span class="Chemical">furan) is known to be a difficult case for HOMA.[85] For these reasons, we also estimated the impact of T1 state Baird aromaticity from geometrical parameters that function as approximate indicators.

The HOMA values for the all-<span class="Chemical">carbon perimeters in the [5]CM’s in their S0 states are low, corresponding to nonaromatic situations (Table ). Yet, they increase in the T1 states with each of the three functionals with the largest increase (ΔHOMA(T1–S0)) in [5]CFU leading to HOMA values (0.75–0.90) that unambiguously correspond to macrocyclic Baird aromaticity. Noteworthy, despite that [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 T1 are also significant in [5]CPP, but for this compound the actual values in T1 are significantly lower (0.44–0.53) than for [5]CFU, and correspond to weaker Baird aromaticity or nonaromaticity. In contrast, the ΔHOMA(T1–S0) for [5]CCHD and [5]CCPD are more modest, although with UB3LYP they both have higher HOMA values in their T1 states than [5]CPP.

Table 3

Perimetric HOMA Values, Λ(C–C–C–C) Dihedral Angle Ranges, and Average Differences in Λ(C–C–C–C) between the S0 and T1 States at (U)B3LYP (Normal Print), (U)M06-2X (italics), and (U)CAM-B3LYP (Parentheses) Levels

	HOMA			Λ(C–C–C–C) range
compound	S0	T1	ΔHOMA(T1–S0)	S0	T1	ΔΛavg(C–C–C–C; T1–S0)
[5]CPP	0.183, 0.248 (0.255)	0.535, 0.529 (0.443)	0.352, 0.287 (0.188)	6.3–25.2, 10.3–25.3 (4.2–26.8)	0.0, 1.3–9.3 (0.9–12.5)	–16.7, –14.6 (−9.8)
[5]CCHD	0.492, 0.247 (0.246)	0.628, 0.347 (0.315)	0.136, 0.100 (0.069)	0.6–11.3, 0.4–17.9 (0.4–12.9)	0.6–12.0, 0.7–15.6 (0.5–17.1)	–3.2, –5.4 (−3.1)
[5]CFU	0.485, 0.393 (0.423)	0.900, 0.833 (0.749)	0.415, 0.440 (0.326)	0.0, 0.1–0.5 (0.2–8.9)	0.0, 0.8–8.9 (0.8–13.9)	0.0, 6.5 (4.9)
[5]CCPD	0.388, 0.125 (0.135)	0.552, 0.276 (0.279)	0.164, 0.151 (0.144)	0.0–27.9, 0.0–35.4 (0.0–34.3)	0.0–18.9, 0.0–29.1 (0.0–26.5)	–13.5, –14.3 (−14.6)

The geometrical impact of T1 state Baird aromaticity is also visible in several other geometric parameters, particularly if one considers the geometric changes when going from S0 to T1. One such parameter, the C–C–C–C dihedral angle between adjacent monomer units (Λ(C–C–C–C)) in [5]CM’s, should decrease upon excitation for systems that adopt Baird aromaticity. Indeed, when going from S0 to T1 the average Λ(C–C–C–C) angles of [5]CPP, [5]CCPD, and [5]CCHD decrease at all three DFT levels, indicative of attainment of some macrocyclic conjugation in the T1 state. However, for [5]CFU, this angle is already small in S0, and there is therefore either no change or only a slight increase upon excitation to T1. Another useful geometrical parameter is 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 [5]CPP and [5]CFU with all three functionals (Table S8a), but it is only with (U)B3LYP that the shortening is of equal magnitude in the five C–C bonds between the M units (the structures are not symmetric with the other functionals). Furthermore, with each of the three functionals these C–C bonds are slightly shorter in [5]CFU than in [5]CPP (1.415 vs. 1.437 Å with UB3LYP and 1.392–1.430 Å vs. 1.405–1.469 Å with UM06-2X), indicating more distinct Baird-aromatic character in [5]CFU.

For the C[1]CM and C[2]CM compounds, we examined approximate HOMA values (HOMAappr) determined based on only the diene segments of the M units that participate in the macrocyclic Baird aromaticity (Table ). At the UB3LYP level, each C[1]CM and C[2]CM is Baird-aromatic while with UCAM-B3LYP this is the case only for C[1]CPP and C[1]CFU. Interestingly, when going from C[1]CM to C[2]CM at UB3LYP level, one can note a significant weakening of the T1 Baird-aromatic character of C[2]CM in the case when M = PP, a minute weakening when M = FU, and a strengthening when M = CPD and CHD. This suggests that T1 state Baird-aromatic character is hampered by a higher proportion of PP units in a macrocycle. It can also be noted that when based on the HOMAappr values at UB3LYP level, C[n]CFU (n = 1 and 2) display slightly higher Baird aromaticity than C[n]CPP (n = 1 and 2). Finally, the C–C(M) bond lengths are shorter in C[1]CFU and C[2]CFU (1.376 and 1.371 Å with UB3LYP) than in C[1]CPP and C[2]CPP (1.401 and 1.386 Å), revealing a stronger macrocyclic conjugation in the first compound class.

Table 4

Approximate HOMA Values (HOMAappr) of the Diene Units of the Monocycles of C[n]CM’s in the S0 and T1 Statesa

	HOMAapprb
compound	S0	T1	ΔHOMA(T1–S0)
C16[1]CPP	0.924, 0.965, (0.975)	0.803, 0.917, (0.953)	–0.121, –0.048 (−0.022)
C12[2]CPP	0.924, 0.965, (0.976)	0.636, 0.917, 0.098 (0.953, 0.118)	–0.288, –c (−)c
C16[1]CCHD	0.624, 0.739, (0.406)	0.451, 0.054, (−0.007)	–0.173, –0.685 (−0.413)
C12[2]CCHD	0.618, 0.404, (0.395)	0.721, 0.556, 0.039, (−0.021, 0.503)	0.103, –c (−)c
C16[1]CFU	0.891, 0.833, (0.823)	0.960, 0.985, (0.933)	0.069, 0.152, (0.110)
C12[2]CFU	0.876, 0.814, (0.804)	0.931, 0.871, (0.934, 0.494)	0.055, 0.057, (−)c
C16[1]CCPD	0.515, 0.355, (0.344)	0.844, 0.213, (0.064)	0.329, –0.142 (−0.280)
C12[2]CCPD	0.504, 0.318, (0.332)	0.925, 0.469, 0.157, (0.016, 0.416)	0.421, –c (−)c

Results with (U)B3LYP (normal print), (U)M06-2X (italics), and (U)CAM-B3LYP (parentheses).

The HOMA values determined for the diene segments in the monocycles that contribute to the macrocyclic Baird aromaticity.

As the HOMA values differ for the different M units within the compound in its T1 state there is no distinct ΔHOMA(T1–S0).

Results with (U)B3LYP (normal print), (U)M06-2X (italics), and (U)CAM-B3<span class="Gene">LYP (parentheses).

The HOMA values determined for the <span class="Chemical">diene segments in the monocycles that contribute to the macrocyclic Baird aromaticity.

Electronic Impact of Macrocyclic Baird Aromaticity

In an idealized (strongly) Baird-aromatic cycle with a Dh or Cv 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 [5]CCPD and each of the three C[n]CCHD’s (n = 1, 2, or 5). Two of these compounds ([5]CCPD and [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 [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 ), and this also occurs for [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 [5]CFU where the HOMA value is well above 0.5 with UCAM-B3LYP.

Figure 3

Spin densities of [5]CFU at UB3LYP, UM06-2X, and UCAM-B3LYP levels at the top row and [5]CPP, [5]CCPD, and [5]CCHD at UB3LYP level at the bottom row. The isodensity surface values in the plots are 0.0040.

Spin densities of [5]CFU at UB3LYP, <span class="Chemical">UM06-2X, and UCAM-B3LYP levels at the top row and [5]CPP, [5]CCPD, and [5]CCHD at UB3LYP level at the bottom row. The isodensity surface values in the plots are 0.0040.

The extent of Baird aromaticity was also assessed through the electronic FLU index[71] for the [4n]annulenes and in the perimeters of the [5]CM macrocycles (Table ). Compounds with polyyne segments were not considered as their FLU values deviate from the others due to their CC triple bonds. According to FLU all [4n]annulenes are T1 Baird-aromatic at both the B3LYP and M06-2X levels. The differences between the FLU values in the S0 and T1 states decrease when going to larger cycles but this is primarily a result of reductions in the values for the S0 state when the annulenes become larger. In T1, the extent of aromaticity is more or less constant according to FLU at UB3LYP level while there is a variation with UM06-2X with 3C8H8 as the most and 3C20H20 as the least Baird-aromatic, in line with similar findings by Matito and co-workers on [4n]annulenes in their T1 states up until C16H16.[82]

Table 5

FLU Values of Compounds in Their S0 and T1 States at (U)B3LYP/6-311G(d,p) (Normal Print) and (U)M06-2X/6-311G(d,p) (Italics) Levels

Compound	S0	T1	ΔFLUT1-S0
C8H8	0.056, 0.062	0.001, 0.001	–0.055, –0.062
C12H12	0.042, 0.051	0.002, 0.008	–0.040, –0.043
C16H16	0.029, 0.041	0.002, 0.006	–0.027, –0.036
C20H20	0.023, 0.034	0.001, 0.011	–0.022, –0.023
[5]CCHD	0.037, 0.046	0.015, 0.026	–0.023, –0.019
[5]CCPD	0.036, 0.046	0.016, 0.028	–0.020, –0.018
[5]CFU	0.018, 0.022	0.009, 0.012	–0.009, –0.009
[5]CPP	0.017, 0.016	0.016, 0.017	–0.000, 0.001

With regard to the [5]CM’s, the one with the most Baird-aromatic character according to FLU at both UB3<span class="Gene">LYP and UM06-2X levels is [5]CFU. With UM06-2X, the FLU value of this macrocycle is similar to that of 3C20H20 while it is slightly less Baird-aromatic than the latter with UB3LYP. The other three [5]CM’s are somewhat less aromatic than [5]CFU at both levels.

Energy Impacts of Macrocyclic Baird Aromaticity

As an energy-based aromaticity index we applied isomerization stabilization energies (ISEs) in versions where the 1,3-hydrogen 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 ). For the compounds with PP monocycles we used dimethyl substitution at one PP unit as the macrocyclic conjugation otherwise cannot be broken.

Figure 4

1,3-Hydrogen shift reactions used to assess the isomerization stabilization energies (ISEs). Blue bonds represent π-conjugated segments: acyclic on the reactant side and cyclic (potentially Baird-aromatic) on the product side.

1,3-Hydrogen shift reactions used to assess the isomerization stabilization energies (ISEs). Blue bonds represent π-conjugated segments: a<span class="Chemical">cyclic on the reactant side and cyclic (potentially Baird-aromatic) on the product side.

For each M, and with both UB3LYP and <span class="Chemical">UM06-2X, one can note gradually decreasing ISE values when n increases in C[n]CM (n = 1, 2 or 5), indicating reduced Baird aromaticity. For [5]CPP and [5]CFU the ISE values are still negative, revealing stabilization upon formation of Baird-aromatic macrocycles in T1, but they are positive for [5]CCHD and [5]CCPD indicating a lack of aromatic stabilization. Also, the ISE values for C[n]CFU’s are always higher than those of the analogous C[n]CPP. As a comparison, the ISE values of 3C8H8 and 3C20H20 at (U)B3LYP ((U)M06-2X) levels are −16.9 (−10.8) and −19.1 (−14.4) kcal/mol, and for benzene in the S0 state (the aromatic archetype), it is −33.9 (−34.2) kcal/mol, resembling earlier computed ISE values for small aromatic monocycles in T1 and S0.[72,86,87] The T1 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 [5]CPP is already low with UB3LYP, and with UM06-2X it is merely 23% of the ISE value of S0 state benzene at the same level. For [5]CFU, on the other hand, the ISE values are essentially equal to those of triplet state C20H20 and approximately half the aromatic stabilization of benzene in S0. 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 C[n]CM. Yet, the results also suggest that a macrocycle [n]CM composed of moderately Hückel-aromatic M units such as furans has a higher likeliness to sustain macrocyclic Baird aromaticity to larger sizes in its T1 state. With UCAM-B3LYP the ISE values of [5]CPP and [5]CFU are −7.8 and −12.5 kcal/mol, respectively, resembling those at UM06-2X level (Table ).

Table 6

Isomerization Stabilization Energies (ISE’s) of C[n]CM in Their T1 States at UB3LYP/6-311G(d,p) (Normal Print) and UM06-2X/6-311G(d,p) (Italics) Levels

n	PP	FU	CPD	CHD
1	–22.9, –14.1	–25.4, –22.3	–14.2, –12.8	–7.4, –4.4
2	–17.4, –10.3	–24.3, –24.1	–9.5, –1.2	–7.3, –5.1
5	–13.1, –7.8	–19.1, –14.2	8.6, 12.1	8.1, 12.5

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?

The ACID plots for the T1 states, calculated with UB3LYP, reveal diatropic ring-currents along the perimeters of all except two macrocycles when going from C[1]CM to [5]CM (Figure and Figures S8–S92), the exceptions being [5]C<span class="Chemical">CPD and [5]CCHD. Slightly fewer macrocycles exhibit diatropic ring currents with UM06-2X as they are observed for C[n]CFU and C[n]CPP (n = 1, 2 and 5), and with UCAM-B3LYP only C[2]CFU and [5]CFU exhibit such ring currents.

Figure 5

ACID plots of (A) C[n]CFU (n = 2 and 5) and (B) C[n]CPP (n = 2 and 5) at UB3LYP/6-311+G(d,p) and UM06-2X/6-311+G(d,p) levels. For the ACID plots in the S0 state and all other compounds, see the Supporting Information.

With regard to NICS, all C[n]CM (n = 1, 2 or 5) with M = PP or FU calculated with UB3LYP or <span class="Chemical">UM06-2X show moderately to highly negative NICS(0) values (Table ). Thus, the NICS results are in line with the induced diatropic ring currents displayed in the ACID plots. Similar as observed with the other indices, the NICS values from UM06-2X indicate slightly lower Baird aromaticity than those at UB3LYP level, except for the NICS of all macrocycles with M = CHD and CPD where UM06-2X indicate lack of aromaticity. Such nonaromatic situations are even found with UB3LYP for [5]CCHD and [5]CCPD. In contrast, the T1 aromatic character according to NICS is preserved in both [5]CPP and [5]CFU according to UB3LYP and UM06-2X while they are just modestly Baird-aromatic according to NICS at UCAM-B3LYP level (NICS(0) = −11.3 and −7.4 ppm, respectively). Interestingly, the NICS values with all three functionals are more negative in [5]CPP than in [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) is found for the C[2]CM species. There is a small attenuation when going to [5]CPP while [5]CFU is as aromatic as C[1]CFU.

Table 7

NICS(0) Values of C[n]CM Compounds in Their T1 States at GIAO/UB3LYP/6-311+G(d,p) (Normal) and GIAO/UM06-2X/6-311+G(d,p) (Italics) Levelsa

n of CM units	PP	FU	CPD	CHD
1	–48.3, −33.8	–34.5, –22.1	–31.9, –10.2	–28.5, –8.8
2	–49.9, –29.5	–45.5, –41.7	–42.9, –5.9	–33.3, –8.4
5	–40.5, –28.3	–31.6, –22.0	–4.8, –3.2	–8.2, 0.5

NICS(0)iso values found in Table S14.

Taken together, the variations in Baird-aromatic character based on magnetic indicators follow largely the same trends as observed by the other indices, except for the opposite order in the extent of aromatic character of C[n]CFU vs. C[n]CPP when based on NICS values.

Synopsis of Results on 3C20–4[n]CM Macrocycles

Conclusions on three aspects can be made based on the computational results above: (i) on the extent of Baird aromaticity in the various C[n]CM species, (ii) on the (dis)similarities in the aromaticity assessments based on the various indices, and (iii) on the differences between the three functionals used. Starting with the last aspect, as already reported,[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 C[n]CM compounds, except for NICS indicating that the Baird aromaticity in C[n]CPP’s to be more markedly Baird-aromatic than the Baird-aromatic C[n]CFU’s.

Finally, with regard to the first aspect listed above, [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 C to [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 T1 state Baird aromaticity, unless further tailored.

T1 State Baird Aromaticity of Larger and Altered [n]CMs

We now focus on larger [n]CPP and [n]CFU (n > 5), and their derivatives, in order to assess at what sizes n these macrocycles cease to be Baird-aromatic. The [n]CFU’s are particularly promising as [5]CFU was Baird-aromatic even with UCAM-B3LYP, in contrast to [5]CPP. Conversely, since [5]C<span class="Chemical">CPD and [5]CCHD are nonaromatic even with UB3LYP, larger [n]CCPD’s and [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 [n]CPP’s with n > 5 have egg-shaped structures which are nonaromatic,[25] and similar attenuations of the Baird aromaticity occur also with UB3<span class="Gene">LYP and UM06-2X but at higher n. The distortion to nonaromatic [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 , S17, and S18). Yet the FLU values, in contrast to NICS and HOMA, reveal no significant weakening of the electron delocalization in the [n]CPP’s when going from n = 6 to n = 8; neither with B3LYP nor M06-2X. That is, also the unsymmetric [8]CPP has FLU values that correspond to Baird aromaticity with both functionals despite that the inter-ring Λ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 [6]CPP, [7]CPP, and [8]CPP, respectively).

Table 8

HOMA, NICS(0), and FLU Values for [n]CM (M = PP or FU, n = 6–8) and at UB3LYP/6-311G(d,p) (Normal Print) and UM06-2X/6-311G(d,p) (Italics) Levels

compound	sym	HOMA T1	NICS(0)zz T1	FLU T1
3[6]CPP	C2h, C1	0.588, 0.424	–41.3, –8.6	0.0150, 0.0184
3[7]CPP	C1, C2	0.504, 0.407	–15.1, –6.2	0.0163, 0.0183
3[8]CPP	C1, C1	0.457, 0.403	–4.7, –5.8	0.0169, 0.0181
3[6]CFU	C6v, C6v	0.949, 0.858	–33.1, –16.6	0.0088, 0.0122
3[7]CFU	D7h, D7h	0.969, 0.841	–36.3, –10.3	0.0085, 0.0131
3[8]CFU	D8h, C1	0.974, 0.808	–38.2, –3.2	0.0084, 0.0143

For the larger macrocycles, the energy differences between the highly symmetric (Baird-aromatic) T1 state structures and the distorted (nonaromatic) ones vary extensively with functional. With UB3LYP, the C4 and D8 symmetric structures of [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 <span class="Chemical">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 S1 state of [n]CPP’s is more delocalized than the T1 state.[25] Similar observations in the T1 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 T1 states even with UCAM-B3LYP? To enhance the Baird-aromatic character in [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 methylene-bridged [6]CPP (MB[6]CPP),[91] a compound that displays magnetically induced paratropic ring currents in the S0 state.[92] We now explored MB[n]CPP’s with n = 6–8 using UB3<span class="Gene">LYP, UM06-2X, and UCAM-B3LYP, and also MB[12]CPP with UCAM-B3LYP. Indeed, with each of the three functionals the MB[n]CPP’s exhibit fully delocalized yet unevenly distributed α-spin density (see Table S9). The HOMA and NICS(0) values of MB[8]CPP with UCAM-B3LYP (0.78 and −40.3 ppm, respectively) clearly reflect Baird-aromatic character. The larger MB[12]CPP calculated with UCAM-B3LYP, however, goes to a structure which similar to [12]CPP has the spin density distribution localized to a part of the macrocycle and Hückel-aromatic benzene rings in the other part.

The Baird-aromatic character also impacts on the adiabatic T1 energies (E(T1)) as seen in Table . For the Baird-aromatic MB[n]CPP’s one can note a gradual increase in E(T1) when going from n = 6 to n = 8 with each of the three functionals. In contrast, the E(T1) 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 Hückel-aromatic character, providing a stabilization of the T1 state with the loss of Baird-aromatic character.

Table 9

Adiabatic Triplet State Energies E(T1) in kcal/mola

macrocycle	(U)B3LYP	(U)M06-2X	(U)CAM-B3LYP
[6]CPP	31.5	43.5	39.1
[7]CPP	37.6	47.5	43.2
[8]CPP	42.2	52.0	46.3
MB[6]CPP	34.1	46.1	43.9
MB[7]CPP	39.3	50.7	47.7
MB[8]CPP	42.2	60.0	62.1
[6]CFU	30.5	43.0	39.0
[7]CFU	30.8	44.1	54.9
[8]CFU	28.8	41.2	37.3

Triplet energies calculated as the absolute free energy differences between the minima in the S0 and T1 states.

To tailor large unconstrained macrocycles which are Baird-aromatic in T1 one may argue that an approach is to utilize M units that can retain Hückel aromaticity in one part of the unit simultaneously as another part is involved in the Baird-aromatic macrocyclic circuit. The <span class="Chemical">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 [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 Hückel-aromatic 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 Hückel-aromatic in S0 and/or impose less steric congestion, and the FU unit is such a unit. Indeed, the [6]CFU and [7]CFU macrocycles have clear Baird-aromatic character with all three functionals, and the breakpoint to nonaromatic structures occurs at [8]CFU. For the latter macrocycle there is just one conformer according to UB3<span class="Gene">LYP (a Baird-aromatic with D8 symmetry) but two with UM06-2X and UCAM-B3LYP; the D8 symmetric conformer and an unsymmetric conformer which are lower in energy by 18.9 and 17.9 kcal/mol, respectively. The unsymmetric conformer of [8]CFU at UM06-2X and UCAM-B3LYP levels has NICS(0) values that reveal complete lack of aromaticity (−3.2 and 1.7 ppm, respectively), and the α-spin density is unevenly distributed. Thus, [8]CFU should likely be considered as nonaromatic despite that the HOMA values with UM06-2X and UCAM-B3LYP (0.81 and 0.78, respectively) are situated within the aromatic range.

Indeed, the existence of a local Baird-aromatic minimum on the T1 potential energy surface of [8]CFU with UB3LYP is interesting as it supports a tug-of-war between a macrocyclic T1 state Baird-aromatic conformer and one conformer having a semilocalized triplet diradical together with one or a few closed-shell Hückel-aromatic FU units. This feature impacts on the E(T1) energies because, similar as for the [n]CPP’s, one can note that T1 state Baird-aromatic [n]CFU’s have higher E(T1) than the ones with nonaromatic global minima in T1 (Table ). 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 T1 PES.

Based on the observation that weakened S0 state Hückel aromaticity together with reduced steric congestion brings increased T1 state Baird aromaticity one may ask if changing the PP units in [n]CPP’s to either of the three diazines (<span class="Chemical">pyridazine, pyrimidine, and pyrazine) can provide for hetero-[n]CPP’s which are Baird-aromatic at larger n’s than the all-carbon [n]CPP’s? These heterocycles have a lower Hückel aromaticity in S0,[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 T1 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 higher-energy local minima vanish with UB3LYP and UM06-2X so that only the nonaromatic structures exist with those functionals.

Neutral Macrocycles in T1 vs Dicationic and Dianionic Ones in S0

Interesting perspectives on macrocycles with potentially Baird-aromatic T1 states can also be gained through comparisons with the corresponding dications and dianions. Several aromatic [n]CPP have earlier been reported,[33,94−97] and based on the Zilberg–Haas valence bond description of triplet state aromaticity (Figure B),[33] one may ask to what extent they resemble [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 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 Hückel MO theory all have the same resonance energy as an equal occupancy of the orbitals which are HOMO and LUMO in the neutral S0 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 in its closed-shell S0 state has a similar NICS value as [5]CPP (with (U)B3LYP: −41.8 and −40.5 ppm, respectively), and the similarity is underlined by nearly identical <span class="Chemical">HOMA values (0.571 and 0.572, respectively). The dianion [5]CPP 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 [n]CPP’s because the dication [7]CPP and dianion [7]CPP with, respectively, 26 and 30 π-electrons in their MCP’s, are still Hückel-aromatic whereas the neutral triplet [7]CPP approaches a nonaromatic situation (Figure ). The same situation is found with M06-2X. In contrast, the trend is different for the [7]CFU, [7]CFU, and [7]CFU as the degree of aromaticity according to each of the three indices used decreases with number of electrons. One can also see that the inter-ring distance is slightly longer in [7]CFU than in [7]CFU and [7]CFU. When going to larger [n]CFU’s, the [n]CFU reduces its aromaticity before the corresponding dication and dianion, and the trend is the same with (U)B3LYP and (U)M06-2X when considering NICS and HOMA values.

Figure 6

HOMA, ISE and NICS(0) values of [7]CFU and [7]CPP in their dicationic and dianionic S0 states and neutral T1 states at (U)B3LYP/6-311G(d,p) and GIAO/(U)B3LYP/6-311+G(d,p) levels. Inter-ring C–C bonds in Å are given below each compound. Symmetries in the optimized geometries are given within each compound.

HOMA, ISE and NICS(0) values of [7]CFU and [7]CPP in their dicationic and dianionic S0 states and neutral T1 states at (U)B3<span class="Gene">LYP/6-311G(d,p) and GIAO/(U)B3LYP/6-311+G(d,p) levels. Inter-ring C–C bonds in Å are given below each compound. Symmetries in the optimized geometries are given within each compound.

It is further notable that the inter-ring CC bonds are shorter in [7]CFU, [7]CFU and [7]CFU than in the corresponding [7]CPP, [7]CPP, and [7]CPP suggesting stronger macro<span class="Chemical">cyclic aromaticity in the former macrocycle class regardless if dication, neutral triplet, or dianion. The ISE values are, however, high for all [7]CM’s except for [7]CPP. For the latter, the ISE values are much smaller than those of [7]CPP, also with UB3LYP, revealing that the macrocyclic Baird aromaticity in the neutral triplet is negligible.

According to both NICS and HOMA at B3<span class="Gene">LYP and M06-2X levels, the aromatic character of [n]CFU2+ and [n]CPP 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 T1 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 [n]CFU and [n]CPP, a [n]CFU macrocycle is slightly more aromatic than a [n]CPP with the same n according to HOMA, ISE, and FLU while it is the opposite between the two species with NICS (Figure and Tables S17 and S18).

Taken together, the macrocyclic dications [n]CM and dianions [n]CM provide upper limits in n for Baird-aromatic character in the T1 state of a [n]CM. If the macrocyclic dication or dianion, which potentially is readily investigated experimentally, is not Hückel-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, <span class="Chemical">furan, cyclohexa-1,3-diene or cyclopentadiene units linked together so that macrocyclic conjugated 4nπ-electron paths are enabled. These macrocycles may all exhibit T1 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 T1 state Baird aromaticity. A second factor is the tug-of-war between macrocyclic Baird aromaticity, on the one side, and closed-shell Hückel aromaticity of one or a few M units together with semilocalized tripled diradical character on the remaining M units, on the other. The 2,5-linked furan unit (FU) is ideal as it does not induce extensive steric congestion, and it has a lower Hückel-aromatic character in S0 than benzene. Hence, regardless of functional used we observe that the T1 state Baird aromaticity is less rapidly attenuated in [n]CFU’s than in [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 T1 state Baird-aromatic to nonaromatic character, leading to different breakpoints n at which the T1 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 S0 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 [n]CPP’s that are Baird-aromatic to have stronger such character than the analogous [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 tug-of-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 <span class="Chemical">porphyrinoid macrocycles could also be interesting as the T1 state Baird aromaticity in such compounds has been explored extensively.[7,13,101−107] For the S0 state, Aihara and co-workers concluded, based on a graph theoretical analysis, that Hückel’s rule is applicable to the main conjugation paths of porphyrinoids having 4n + 2 π-electrons.[108−114] Yet, the macrocyclic aromaticity in S0 was significantly weaker than the aromaticity of the individual monocycles. How does this change in the T1 state? These aspects of macrocyclic T1 Baird aromaticity will be reported in forthcoming studies.