Global aromaticity at the nanoscale.

Aromaticity can be defined by the ability of a molecule to sustain a ring current when placed in a magnetic field. Hückel's rule states that molecular rings with [4n + 2] π-electrons are aromatic, with an induced magnetization that opposes the external field inside the ring, whereas those with 4n π-electrons are antiaromatic, with the opposite magnetization. This rule reliably predicts the behaviour of small molecules, typically with fewer than 22 π-electrons (n = 5). It is not clear whether aromaticity has a size limit, or whether Hückel's rule extends to much larger macrocycles. Here, we present evidence for global aromaticity in porphyrin nanorings with circuits of up to 162 π-electrons (n = 40); aromaticity is controlled by changing the constitution, oxidation state and conformation. Whenever a ring current is observed, its direction is correctly predicted by Hückel's rule. The largest ring currents occur when the porphyrin units have fractional oxidation states.

The extent of electronic delocalisation in linear molecules is limited by the onset of symmetry-breaking transitions, which can be viewed as Peierls-type electron-vibration interactions or as shifts in mixed-valence behaviour[1]. For example, cyanine dyes feature a linear chain of C–C bonds with bond order 1.5 and negligible bond length alternation (like the C–C bonds in benzene) resulting in charge delocalisation, but if the chain exceeds a critical length, the symmetry collapses, localising the charge[2,3]. It is not clear whether similar effects limit the size of an aromatic ring, or whether molecular ring currents can extend into the domain of mesoscopic phenomena such as Aharonov–Bohm oscillations[4]. Many new globally aromatic macrocycles have been reported during the last few years[5-18], but apart from the porphyrin nanorings discussed here, there are no reports of aromatic rings with more than 62 π-electrons[7].

Here we explore circuits of up to 162 π-electrons in a large family of nanorings in a wide range of oxidation states. In these nanorings, each porphyrin contributes 10 electrons to the Hückel π-electron count, and each linking alkyne contributes 2 electrons, so a nanoring cation has an electron count of 10N + 4x + 2y – Q (where N is the number of porphyrin units; x and y are the number of butadiyne and ethyne links respectively). Recently, we reported that the butadiyne-linked six-porphyrin nanoring displays global aromaticity when oxidised[8], reduced[17] or electronically excited[18]. According to Schleyer’s terminology[19], all the nanorings discussed here are ‘trannulenes’ rather than ‘annulenes’, since the p-orbitals making up the aromatic π-system are oriented in the plane of the ring, whereas in annulenes the p-orbitals are perpendicular to the ring. Trannulenes follow Hückel’s rule in the same way as annulenes, although Hückel did not anticipate this type of system.[20]

Results and discussion

The six-porphyrin nanoring complexes , , and -P6[b provide a homologous series of compounds in which we systematically vary the number of π-electrons by changing the number of -C≡C- units, while preserving the circular geometry (which is locked by the template, T6* or T6, Fig. 1)[21,22]. The 1H NMR spectra of all four complexes were recorded in the 2+, 4+ and 6+ oxidation states, revealing the presence of aromatic or antiaromatic ring currents. The directions of these ring currents agree perfectly with Hückel’s rule. Thus, (neutral: 72 πe) and (neutral: 84 πe) are both aromatic in the 2+ and 6+ oxidation states and antiaromatic in the 4+ state, whereas (neutral: 74 πe) and (neutral: 82 πe) are both antiaromatic in the 2+ and 6+ states and aromatic in the 4+ state.

Fig. 1

1H NMR spectra of the aromatic and antiaromatic six-porphyrin nanoring template complexes in oxidation states 2+, 4+ and 6+.

a, *; b, *; c, and d, . A grid plot of the NICS(0) value in the xy plane of the nanoring, calculated without template, is shown for each oxidation state of each complex (5 × 5 nm; LC-ωhPBE/6-31G*, ω = 0.1; colour axis truncated above 20 and below –20 ppm; contours drawn every 5 ppm). 1H NMR spectra recorded at 500 MHz in CD2Cl2; oxidised states are generated by titration with thianthrenium hexafluoroantimonate. # and ‡ denote CHDCl2 and thianthrene, respectively. Detailed spectra are shown in Supplementary Figs. 9–12 and 19–34. Dashed vertical lines indicate 10-fold magnifications.

The most obvious evidence for these ring currents comes from the chemical shifts of the template α and β pyridyl 1H resonances; for example, these protons are strongly shielded (δH = –7.3 and –3.6 ppm, vs. 8.7 and 7.5, respectively in the free template) in aromatic 2+ and strongly deshielded (δH = 36.0 and 30.1 ppm, respectively) in antiaromatic 4+. Further evidence is provided by the trihexylsilyl (THS) 1H and 13C signals. In each spectrum, we observe one group of THS signals that is essentially unshifted, at 0–2 ppm as in the neutral compounds (peaks coloured orange in Fig. 1) which arises from the external THS group located near the zero-shielding cone of the nanoring (THSout, Fig. 2), and one group of THS signals that is shielded or deshielded, depending on the direction of the ring current (THSin, coloured green in Fig. 1 and 2). Interconversion of THSin and THSout is slow on the NMR timescale, and the assignment of THSin signals is confirmed by the observation of NOEs to protons of the template (T6 or T6*). There is an excellent linear correlation between changes in the chemical shift of the α-pyridyl template 1H and 13C signals and the THS CH2-Si and CH3 1H and 13C signals, showing that all six signals report on the same global ring currents (Supplementary Figs. 70–73).

Fig. 2

The magnetic shielding plotted in the xz plane perpendicular to the plane of the nanoring.

NICSiso contours for ( from –10 to 10 ppm. Internal trihexylsilyl groups (THSin, green) are sensitive to the global ring current whereas external ones (THSout, orange-brown) are not.

The observed ring currents in this set of 12 species (four nanorings in three oxidation states) were compared with the results of nucleus-independent chemical shift (NICS) calculations[23,24]. We screened a range of density functional theory (DFT) functionals for these NICS calculations (see Supplementary Figs. 80 and 83). The B3LYP functional[25] gave good agreement with the NMR results for most species, but it failed to reproduce the anti-aromaticity of the dications and . The LC-ωhPBE (ω = 0.1) functional[26] gave the best agreement with the observed NMR shifts (see plots on the right of each spectrum in Fig. 1) and we chose this functional for all the NICS calculations shown in this article.

The evolution of the ring current with increasing ring size is illustrated by the 1H NMR spectra of the 8-porphyrin nanoring complex , which has a circuit of 96 πe when neutral (Fig. 3)[21]. The THS and template protons show clear evidence for aromaticity in the 2+ and 6+ oxidation states, whereas the 4+ and 8+ oxidation states are antiaromatic. These results match the predictions of Hückel’s rule and agree with calculated NICS values (Fig. 3). The magnitude of the ring current varies substantially between different oxidation states; thus the mean change in the chemical shift of the α-protons of the template (relative to unbound T4*, δα = 8.67 ppm) is Δδα = +11.8 ppm in the 4+ state, and –15.6 ppm in the 6+ state, but it dwindles to +1.9 ppm in the 8+ state.

Fig. 3

Hückel behaviour in a template-bound eight-porphyrin ring.

Molecular structure of , Ar = 3,5-bis(trihexylsilyl))phenyl, and 1H NMR spectra in oxidation states 2+, 4+, 6+ and 8+. Labels denote the key resonances THSin (green), THSout (orange), and template (α, β); # and ‡ denote CHDCl2 and thianthrene, respectively. Detailed spectra are shown in Supplementary Figs. 14 and 45–47. NICS(0) grids (LC-ωhPBE/6-31G*, ω = 0.1) in the xy plane for each state (without template).

Fluorinated templates allow the aromaticity of nanorings to be evaluated using 19F NMR, as exemplified by the extended six-legged template T6ef (Fig. 4a). Two molecules of this template stack to form a stable 2:1 complex with the nanoring, (circuit electron count: 168 πe when neutral[27]), in which the CF3 groups are positioned to probe the global ring current. The CF3 19F resonance is shielded in the aromatic 6+ and 10+ oxidation states, but deshielded in the antiaromatic 8+ oxidation state (Fig. 4b). The shifts of the CF3 19F resonance are fully consistent with shifts of the THSin 1H signals (shaded green, Fig. 4b) and they agree with the predictions of NICS calculations. A range of DFT functionals were tested for and, as for the 6-porphyrin rings, LC-ωhPBE (ω = 0.1) was found to give the best agreement with experimental ring current shifts (see Supplementary Figs. 88–90 and Supplementary Tables 7 and 8).

Fig. 4

Ring currents in topologically distinct twelve-porphyrin ring complexes.

a, Molecular structure of with two stacked bound template units. b, 1H and 19F NMR spectra of oxidised ; green shading indicates the interior THSin resonance; backgrounds indicate a global aromatic (blue) or antiaromatic (orange) state. c, A circular conformation is predicted to exhibit a ring current and thus an induced magnetic moment in an external magnetic field. d, Molecular structure of figure-of-eight complex. e, 1H and 19F NMR spectra of oxidised . f, In the lemniscate, the two loops induce opposite currents which cancel, and no ring current is expected. Detailed spectra are shown in Supplementary Figs. 17, 18 and 58–67. Lk, Wr and Tw are the linking number, writhe and twist that specify the topology of the loop[32–34]. Axy, Axz and Ayz are the net cross section areas; for , Axy = 0 because the lemniscate has D2 symmetry and the areas of the two loops cancel.

No aromatic or antiaromatic ring current was detected for the 12+ oxidation state, and this result is reproduced by the NICS calculations, although this 12+ state is expected to be antiaromatic (156 πe; 4n; n = 39). The 19F NMR titrations also show CF3 signals attributed to the open-shell 7+, 9+ and 11+ oxidation states, at similar chemical shifts to the neutral compound — thus the shielding or deshielding effects in the open-shell cations are small compared with those in closed-shell species. EXSY NMR experiments show that the odd-electron oxidation states are in chemical exchange with the even-electron states, on a timescale of seconds. It is not surprising that mixtures of oxidation states are formed during these titrations, because the oxidation potentials are closely spaced (for calculated speciation curves, see Supplementary Fig. 79), but it is remarkable that these open-shell species give spectra that are sharp enough to be observed.

We used the 12-porphyrin nanoring to explore the relationship between 3D conformation and aromaticity. The magnitude of the ring current induced in a macroscopic ring of metal wire depends on the total magnetic flux passing through the ring. If the ring has a figure-of-eight shape, with two equal lobes such that the magnetic flux passing through each loop induces equal and opposite currents, then there will be no net ring current. We sought to test whether this principle applies on the molecular scale, so we synthesised a small six-legged template, T6f, which forms a figure-of-eight shaped[27] 1:2 complex (which is doubly twisted, not Möbius[28,29], Fig. 4d) and investigated the ring currents in this system by 1H and 19F NMR as a function of oxidation state, under identical conditions to those used for the circular . The resulting 1H, 13C and 19F NMR spectra show the absence of any detectable ring currents (Δδ <0.1 ppm) in the figure-of-eight nanoring, confirming that aromaticity can be switched on/off by geometry (Fig. 4e). This result is reproduced by the NICS calculations (Supplementary Fig. 91).

The suppression of ring currents in figure-of-eight shaped annulenes has been predicted theoretically[28,30], but it has not been observed experimentally in other figure-of-eight shaped aromatic systems[6,7,14,28,29,31], probably because they were not the right shape to achieve cancellation of the ring current. The topology of a closed ribbon can be described by the linking number Lk, the writhe Wr and the twist Tw (ref. 32–34). Total cancellation of the ring current is expected for an ideal lemniscate geometry with Wr = 2, Tw = 0 and Lk = 2 with D2 symmetry (Fig. 4f)[28], which is close to the geometry of (Wr = 1.8, Tw = 0.2 and Lk = 2, from the crystal structure of , ref: 35). The ring has a large cross sectional area (Axy = 22 nm2; Axz = Ayz = 0; Fig. 4c) resulting in a substantial ring current, whereas the lemniscate has a small net cross section (Axy = 0, Axz = 2.3 nm2, Ayz = 0.4 nm2 and Supplementary Fig. 93) resulting in a weak response to magnetic field. The ring current in is blocked by the global topology (Tw ≈ 0 and Wr ≈ 2), not by any local break in π-conjugation.

The ring currents observed in this whole family of nanorings are summarised in Fig. 5, which plots the shift in the THSin resonances (ΔδTHSin) as a function of oxidation state. Whenever a ring current is observed, its direction (aromatic or antiaromatic) matches the prediction from Hückel’s rule. The magnitude of the ring current varies with the average oxidation state of the porphyrin units The largest ring currents are observed in mixed-valence systems, where Global ring currents are not observed in the neutral rings where the local porphyrin ring current dominates; in the larger rings, the ring current also vanishes when (see square points for and in Fig. 5). The formation of a mixed valence state appears to be essential for efficient nanoscale charge delocalisation, just as the presence of a partially filled band is essential for conductance in an extended lattice[36].

Fig. 5

Summary of the shielding and deshielding of the trihexylsilyl groups across eight different nanorings.

Schematic representation of related porphyrin rings; red indicates ethyne, black indicates butadiyne links; black dots indicate zinc-porphyrins; templates in blue. Plot of the observed 1H chemical shift difference ΔδTHS = δTHS(in) – δTHS(out) between inner and outer THS protons for each oxidation state of the nanoring. Positive or negative shifts indicate a global antiaromatic or aromatic current, respectively. Empty circles indicate neutral rings (ΔδTHS = 0); squares indicate oxidation states with . Vertical dashed lines denote oxidation states with 4n +2 π-electrons. Lines connecting the points are for visual guidance only. ΔδTHS is plotted for the CH3 (dark points) and SiCH2 signals (faded points) of the THS chains.

Hückel’s rule was originally formulated to explain the unusual properties of benzene, and other molecules with 6 π-electrons[37]. It is remarkable that this simple rule correctly predicts the magnetic response of large oxidised nanorings with circuits of up to 162 π-electrons. This work shows that electronic delocalisation can extend coherently around molecular rings with circumferences of 16 nm. These supramolecular rings allow the magnitude of the ring current to be controlled by topology in a way that has not yet been demonstrated for small molecules.

Methods

All 1H NMR oxidation titrations were carried out by adding a well-stirred suspension of thianthrenium hexafluoroantimonate to a solution of the porphyrin nanoring in CD2Cl2 at –60 °C to –20 °C under a counterflow of argon (see Supplementary Material for details). These oxidation experiments were conducted in NMR tubes fitted with J. Young greaseless PTFE stopcocks using CD2Cl2 stored over molecular sieves, using standard Schlenk line techniques to exclude moisture. Exposure to water is immediately deleterious to porphyrin polycations, but they are stable to oxygen (O2). The neutral nanorings are only sparingly soluble in CD2Cl2 at low temperatures (< –40 °C), but solubility improves dramatically upon oxidation. The NMR measurements were performed on a Bruker AVII 500 (5 mm BBFO probe). 1H NMR chemical shifts were calibrated to residual proton signals of the solvent (CHDCl2 5.32 ppm). 13C NMR chemical shifts were referenced to the solvent peak (CD2Cl2 54.00 ppm). 19F NMR spectra were referenced to hexafluorobenzene (–164.8 ppm). At the end of each titration, decamethylferrocene (FeCp2*) was added to reduce the porphyrin nanoring back to its neutral form; the whole oxidation process is highly reversible.

Supplementary Material

Supplementary information is available for this paper at https://doi.org/10.?????