The mesityl-substituted TFF derivative 4a was obtained in a four-step procedure described in the Supplementary Information and completely characterized. In the solid state (Fig. 2A), the two diindenofluorenylidene (DIF) subunits of 4a form a twist angle θeq of 50.5°, which is larger than the corresponding torsion in 1a (42°), yet smaller than in the sterically congested 1b and 2. However, with a length of 1.431(7) Å, the central alkene bond a is significantly elongated in comparison with unhindered 9,9'-bifluorenylidene (BF) derivatives (1.36–1.38 Å),[35, 37, 48] and is even longer than the corresponding bond in 2 (1.40–1.41 Å).[40] Since the b bonds in 4a (1.422(3) Å) are also shorter than those in the reported unhindered BF systems (1.46–1.48 Å), the central double bond in 4a is weakened by conjugation within the DIF subunits. Interestingly, distance c, which corresponds to a formal single bond in the closed-shell configuration of 4a, is shorter than the formally double bond d (1.398(4) Å vs. 1.456(4) Å). Thus, the solid-state bond length pattern suggests considerable open-shell contributions to the electronic structure of TFF.
Compound 4a, which is NMR-silent in solution, revealed a persistent ESR signal in the solid samples of 4a, whose intensity decreased with the decreasing temperature. The temperature dependence of the ESR signal was modeled using the Bleaney–Bowers equation, yielding an estimate of the singlet–triplet gap ΔES–T of − 1.5(1) kcal/mol. A SQUID analysis confirmed that 4a is a ground-state singlet, with ΔES–T = − 3.1(1) kcal/mol calculated using the Bleaney-Bowers model. Absorption spectra of the deep-blue 4a contained several strong maxima in the 400 to 800 nm range and a weak band tailing beyond 1000 nm, corresponding to a relatively small energy gap of ca. 1.2 eV. In line with the latter observation, 4a displayed pronounced redox amphoterism in electrochemical experiments (Supplementary Fig. 1). In differential pulse voltammetry, three oxidation events were found at 0.01 V, 0.16 V, and 1.03 V (vs. Fc+/Fc in dichloromethane). Reduction of 4a occurred at − 1.07 V (2e), − 1.90 V (1e), and − 2.01 V (1e), respectively.
Titration of 4a with a one-electron oxidant, tris(4-bromophenyl)ammoniumyl hexachloroantimonate (BAHA, Eox = 0.7 V in DCM) revealed consecutive formation of two NIR-absorbing species, with each step yielding near-perfect isosbestic points (Fig. 3A). The neutral 4a could be quantitatively recovered by reduction with KO2. The two oxidation products were also generated electrochemically (Supplementary Figs. 15–16), and identified, respectively, as the radical cation [4a]•+ (λmax = 860 nm) and dication [4a]2+ (λmax = 985 nm). The absorption of the radical cation features a shoulder at ca. 1200 nm, implying a smaller energy gap in [4a]•+ than in [4a]2+. In contrast to the neutral 4a, the dication [4a]2+ is diamagnetic. In its 13C NMR spectrum, the C-11 resonance was identified at a highly downfield position of 184 ppm, indicating partial localization of the positive charge at this site. The structure of the dication was further elucidated in an X-ray diffraction analysis of a single crystal of the [4a][SbCl6]2 salt (Fig. 2B). In comparison with the neutral 4a, the a bond is somewhat lengthened (to 1.462(6) Å), while the θeq torsion increases to 54.5–55.8°. These changes suggest that the two-electron oxidation results in a moderate decrease of the a bond order.
Interestingly, even though a two-electron reduction event had been revealed in electrochemical measurements, titration of 4a with cobaltocene (CoCp2, Ered = − 1.3 V) provided evidence for stepwise electron transfer (Fig. 3B'). Specifically, upon addition of up to 4 equiv of CoCp2, we observed initial formation of the radical anion [4a]•– (λmax = 685, 920, 1155, ~ 1355 nm). Further addition produced the dianion [4a]2– (λmax = 760, 1120 nm), which could be oxidized back to 4a using diiodine. A broader range of anionic states of 4a was achievable by reduction with sodium naphthalenide (NaN), in the presence of 15-crown-5 (Figs. 3B and 3C). Initial spectra, observed in the range of 0–9 formal equivalents of added NaN corresponded to the sequential formation of [4a]•– and [4a]2–. On further addition of NaN (ca. 17 equiv), we observed the formation of a new species, with NIR absorptions extending beyond 2000 nm, which was presumed to be the radical trianion [4a]•3–. When an even larger excess of NaN was added (up to 33 equiv), these characteristic bands disappeared, and the final spectrum had an absorption onset at ca. 1200 nm. Thus, the ultimate reduced product had a larger energy gap than all the preceding forms and was tentatively assumed to be the tetraanion [4a]4–. Both [4a]4– and the dianion [4a]2– could be selectively generated on a larger scale and characterized using 1H NMR (Supplementary Figures S2 and S7–S11).
Single crystals containing [4a]2– and [4a]4– anions were grown from THF solutions of 4a reduced with sodium metal in the absence of the crown ether additive. In particular, the tetraanion structure [Na(THF)3]4[4a] revealed a highly regular pattern of Na cations coordinated directly to the π system (Fig. 2D). The cations are bound near the fused edges of the BF core, with the shortest Na···C distances of 2.636(6) Å. Apparently, binding to the five-membered rings, which are presumed to carry a significant portion of the negative charge, is not feasible because of the steric protection by the Mes substituents. The a bond in the tetraanion is stronger than in other oxidation levels of 4a, as evidenced by its short length of 1.35(1) Å and the smaller θeq torsion of 40.9°. Partial reduction of 4a yielded crystals with a stoichiometry of [Na(THF)6][Na(THF)5]0.74[4a]·8.3THF, indicating a mixed-valence character of 4a. Sodium occupancies indicate that the crystal contains mostly the dianion [4a]2– with a ca. 26% admixture of the radical anion [4a]•–. The structure is notable for the lack of Na···π coordination, reflecting the smaller negative charge residing in the π system. Specifically, two non-equivalent sodium sites were found: octahedral [Na(THF)6]+, and disordered square-pyramidal [Na(THF)5]+. The apparent geometry of 4a is averaged over the two contributing redox states (–1 and − 2) and features a relatively long bond a (1.436(6) Å) and a large θeq torsion (58.1°). These parameters may indicate weaker conjugation between the DIF subunits in [4a]2– and [4a]•– than in [4a]4–.
Computational analysis. DFT calculations performed for the substituent-free TFF molecule 4b (R = H, Fig. 1) revealed significant variations of the key geometrical parameters as a function of the charge and multiplicity of the system (Supplementary Table 1). At the UCAM-B3LYP/6-31G(d,p) level of theory (hereafter denoted CAM), the equilibrium torsion, θeq, and the central bond distance a ranged from 33.7° and 1.372 Å in 5[4b] to 90.0° 1.477 Å in 3[4b]4–, respectively. The a distance shows good correlation with θeq, indicating that both parameters can be used to quantify the strength of inter-subunit interaction in TFF. Relaxed potential energy surface (PES) scans along the θ coordinate revealed a complex dependence of the energy profile on the charge and multiplicity of 4b (Fig. 4, Supplementary Table 1). Specifically, all singlets and doublets have a twisted equilibrium geometry, characterized by θeq < 60°, and a transition state at θ = 90°. The singlet dication 1[4b]2+ features the lowest twist barrier (ΔEtwist = ΔErel(90°) – ΔErel(θeq) = 0.8 kcal/mol) and the least acute θeq angle (59.2°). The ΔEtwist barriers increase in the order 1[4b]2+ < 1[4b]2– < 2[4b]+ < 1[4b] < 2[4b]– < 2[4b]3– < 1[4b]4– and correlate with a decrease of the respective θeq angles. The above sequence can thus be assumed to reflect an increase of the bond order a. PES scans of the triplets lie above the corresponding singlet scans at all θ angles: a single-well potential with θeq = 90° is predicted for 3[4b]4–, while double-well potentials are found for 3[4b]2+ and 3[4b]2–. The latter two species are significantly destabilized relative to the respective singlets, however, their twist barriers ΔEtwist are actually higher, implying that the a bond becomes stronger in these two triplet states. Conversely, the neutral triplet, 3[4b], has a shallow minimum at θeq = 58.2° with a low ΔEtwist barrier of 1.1 kcal/mol, implying a particularly weak a bond. This behavior could be considered typical of an alkene, however, in the quintet state 5[4b], the bond is predicted to be very strong, with a much higher twist barrier of 15.7 kcal/mol. This unusual response of the inner alkene in the neutral 4b to changes of spin multiplicity is further confirmed by the torsional dependence of a and b bond lengths (Supplementary Fig. 16). The triplet 3[4b] and quintet 5[4b] approach the limits of a pure single and double a bond, respectively, whereas the singlet 1[4b] contains a strongly conjugated alkene with an intermediate bond order.
Adiabatic singlet–triplet gaps predicted for even-electron ions of 4b at the CAM level are relatively large (ca. − 9 to − 15 kcal/mol, Fig. 4, Supplementary Table 1), in line with the observed diamagnetism of [4a]2+, [4a]2–, and [4a]4–. For the neutral 4b, a significantly smaller gap was obtained (ΔEST = − 3.91 kcal/mol), relatively close to the experimental SQUID value, whereas the quintet state 5[4b] was predicted to have a much higher energy (ΔESQ = − 14 kcal/mol). Since CAM, as a single-reference method, is not fully suitable for quantitative assessment of spin-state energetics, we evaluated energies of the neutral m[4b] (m = 1, 3, 5) using two active-space methods, i.e. CAS-SCF(6,6)/6-31G(d,p) (denoted CAS), and the spin-flip approach[49, 50] at the RAS(4,4)-SF-srB3LYP/cc-pVDZ level of theory (denoted RAS). RAS calculations, performed for the CAM-optimized minima and PES scans, indicate that the ST and SQ gaps may be smaller than predicted by the CAM level (ca. − 3 and − 7 kcal/mol, respectively, Supplementary Figure S18).
An analysis of natural orbital occupation numbers (NOONs) showed the neutral singlet 1[4b] had a pronounced tetraradicaloid character (y0 ≥ 0.98 and y1 ≥ 0.29), with a possible smaller hexaradicaloid contribution (y2CAM = 0.12). The number of unpaired electrons obtained from the CAM-derived NOONs (nUCAM, Supplementary Table 1) is non-zero for all states except for 1[4b]4–, which is the only one with a purely closed-shell configuration. The nUCAM value of 3.15 obtained for the neutral singlet 1[4b] is in fact higher than in the corresponding triplet state. Extensive mixing of open-shell configurations is indicated by the high values of nUCAM (i.e., exceeding m – 1) obtained for 2[4b]+ and 2[4b]–. Since twisting of an alkene normally leads to π bond breaking, one could intuitively expect the nU values to increase with increasing θ. Remarkably, however, no such general relationship is found for 4b (Supplementary Fig. 17). Paradoxically, nU decreases with θ for the neutral singlet 1[4b], as well as for 2[4b]–, 2[4b]+, 1[4b]2–, and 1[4b]2+, suggesting that in these species, electron pairing is actually enhanced by decoupling of the DIF subunits.
Nucleus-independent chemical shifts (NICS) revealed striking variations of magnetism in 4b caused by changes of its oxidation and spin state (Fig. 5 and Supplementary Fig. 23). Remarkably, while the triplet state 3[4b] is less aromatic than the singlet 1[4b], the quintet 5[4b] state shows enhanced aromaticity. The NICS map obtained for the triplet is essentially identical with the map obtained for the diindenofluorenyl radical 2[DIF-H], consistent with the weak interaction between DIF subunits in 3[4b] found in the PES scan. Thus, the enhancement of aromaticity in the singlet and quintet originates from the stronger inter-subunit coupling in each of these two spin states. This conclusion is supported by the harmonic oscillator model of aromaticity (HOMA, Supplementary Table 2), which produced significantly higher indexes of rings B in 5[4b] (0.89) and 1[4b] (0.71) than in 3[4b] (0.63). The doubly charged 1[4b]2+ and 1[4b]2–, show opposite changes of their magnetism, being respectively para- and diatropic. The tetraanion 1[4b]4– experiences a dramatic increase in the aromaticity of rings C–D, confirming that it can be treated as a union of four fluorenyl anions. In line with this interpretation, the HOMA indices for rings B, C, and D, are 0.76, 0.42 and 0.76.
Valence structure of TFF. A unified description of the valence structure of TFF, which is valid for all oxidation levels, can be developed using the five partial canonical structures A through E*** shown in Fig. 6. These structures differ in (a) the number of Clar sextets, (b) the number of formally non-bonding sites (denoted with an asterisk), and (c) the order of the linking bond (single or double). For instance, there are no Clar sextets and no non-bonding sites in structure A, whereas structure B* features three sextets and one non-bonding site, i.e., either a cation, a radical, or an anion. Valence structures of TFF can be constructed from pairs of partial canonical contributors with matching linking bond orders, e.g. A + A or B* + B*. Each such structure can be characterized by the total number of sextets (NCS) and the number of unpaired electrons (NUE). Given that the stability of the structures is expected to increase for high NCS values and low NUE values, one can consider a simple stability metric ΔN = NCS – NUE.
Data obtained for the neutral TFF indicate that the triplet state 3[4b] is well approximated by the diradical structure B•–B•, containing a single bond between the DIF units. Similarly, the tetraradicaloid form D••=D•• provides an accurate representation of the quintet 5[4b]. Mixing of these two contributions in 1[4b] can be proposed to explain the intermediate aromaticity, inter-subunit bond order, and high tetraradicaloid character of the singlet state. Interestingly, B•–B• and D••=D•• have the highest ΔN = 4 among all canonical forms, which explains their relative importance. Doubly charged TFF ions 1[4b]2+ and 1[4b]2– can be similarly characterized with singly bonded structures B+–B+ and B––B–, respectively, however, small contributions of the doubly bonded forms D•–=D•– and D•+=D•+ need to be invoked to explain the non-zero nU values and non-vanishing inversion barriers of these two species. The latter two forms should become dominant in the structures of respective triplets, 3[4b]2+ and 3[4b]2–, explaining the high ΔEtwist values predicted for these species. Analogous contributions become even more relevant in the singly charged 2[4b]+ and 2[4b]– (D••=D•– and D••=D•+, respectively), both of which have significant triradicaloid character (nU > 2.5). Finally, the doubly bonded contributors fully dominate in the higher anions, 2[4b]3– and 1[4b]4–, which feature the highest ΔEtwist barriers, and very low nU values.