To realize the synthesis of [6]PA with sufficient stability and solubility, our strategy is to introduce four bulky and electron-withdrawing 2,6-dichlorophenyl groups at the most reactive sites of the two zigzag edges and to incorporate C(sp3)-bridged cyclopenta (CP) rings at the bay regions of [6]PA where each sp3 carbon substituted by two 3,5-di-tert-butylphenyl groups. In addition to these substituents, four additional 3,5-di-tert-butylphenyl groups are also attached onto the π-conjugated framework to further improve solubility. As depicted in Scheme 1, the key intermediate for the synthesis of the target molecule 1 is the partially fused dibromo- intermediate 7, which is prepared via multiple cross-coupling, regioselective cyclization and oxidative dehydrogenation reactions. The CP ring-fused perylene 2 was synthesized from parent perylene according to our reported procedre.11 The monobrom CP-fused perylene 3 was firsty synthesized by regio-selective bromination with bromine under controlled conditions at the bay-region, and then Miyaura borylation reaction of 3 with bis(pinacolato)diboron (B2(pin)2)in the presence of palladium catalyst provided monoboronic ester 4 in 63% yield. After that, the Suzuki coupling of 5 with 4 yielded the intermediate 6 in 68% yield. Next, the key intermediate 7 was obtained by the treatment of 6 with CuBr2 in the presence of base where a one-pot reaction of
benzannulation and intramolecular oxidative cyclodehydrogenation was achieved. The structure of 7 was confirmed by single-crystal analysis (Figure S9). Afterwards, lithium-halogen exchange reaction by treatment of 7 with n-BuLi proceeds and subsequent quenching with 2,6-dichlorobenzaldehyde to generate the diol, which underwent an intramolecular Friedel-Crafts cyclization in the presence of trifluoromethanesulfonic acid (TfOH) to afford the dihydro- precursor 8 in 53% yield over three steps from 7. Finally, the deprotonation of 8 with tBuOK and 18-crown-6 in THF followed by the oxidation with p-chloranil furnished the target peri-hexacene derivative 1 in 73% yield. Beacause two cyclopenta (CP) ring was attached onto the arm-chair edges and the four bulky and electron withdrawing 2,6-dichlorophenyl groups are attached onto the most reactive zig-zag edges, compound 1 shows good solubility and stability, and can be purified by triethylamine deactivated silica gel column. Stability test of the solution of 1 in carbon tetrachloride under ambient air and light conditions followed by UV/Vis-NIR absorption measurements revealed a half-life time (t1/2) of about 24 h when monitored at 617 nm (Figure S1 in the Supporting Information (SI)).
Single crystals of 1 suitable for X-ray crystallographic analysis were grown by slow diffusion of acetonitrile into its solution in toluene.12 The molecule has a centrosymmetric geometry with a slightly distorted skeleton (Figs. 2a,b). The dihedral angles between the π-conjugated skeleton and its 2,6-dichlorophenyl substituents range from 83.1 to 84.5o (Fig. 2b). Bond length analysis of the 1 backbone revealed that the bond lengths of the CC bonds (the bonds shown in red, Fig. 2c) linking the two hexacene units were in the range of 1.409–1.429 Å, which are shorter than that in [4]PA derivative (1.442–1.459 Å),6 indicating enhanced electronic coupling between two hexacene units. The CC bonds shown in blue (Fig. 2c) are much longer (1.464 Å) like a typical C(sp2)-C(sp2) single bond, which could be attributed to the elongation induced by the steric hindrance between the neighboring 2,6-dichlorophenyl substituenets. The calculated harmonic oscillator model of aromaticity (HOMA) values13 based on the X-ray structures revealed that the four benzenoid rings at the termini (rings A/F/L/Q) have largest HOMA values (0.87/0.84, Fig. 2c), indicating a large aromatic character. The rings C/D/N/O have much smaller HOMA values (0.38 ~ 0.44, Fig. 2c) while the other benzenoid rings have a HOMA value of 0.65 ~ 0.80 (Fig. 2c), indicating significant contribution of the open-shell diradical form B with five localized aromatic sextet rings shown in Fig. 1d to the ground-state structure. On the other hand, nucleus independent chemical shift (NICS)14 calculations show that these five benzenoid rings have large negative NICS(1)zz values (-20.49 ppm ~ -22.20 ppm), the rings B/E/H/J/M/P have moderate negative NICS(1)zz values (-11.93 ppm ~ -16.00 ppm) while the rings C/D/N/O are almost non-aromatic (Fig. 3a). All these suggest that [6]PA possesses a large diradical character, which can be stabilized by five aromatic sextet rings. In addition, no intermolecular π-π interaction was observed in the 3D packing structure (Figure S10).
To gain a deeper understanding of the electronic structure and aromaticity of compound 1, we carried out anisotropy of induced current-density (ACID)15 and isochemical shielding surface (ICSS)16 calculations (B3LYP/6-31G(d,p)). ACID plot shows obvious clockwise diatropic ring current circuit along the periphery (Fig. 3b), indicating that the molecule is globally aromatic and the two hexacene units are coupled well with each other, in accordance with the bond length analysis. Notably, the central anthracene units show clockwise diatropic ring current (Fig. 3b), in agreement with the HOMA value and NICS analysis. Furthermore, the calculated 2D ICSS map also shows that the rings A/F/L/Q have the strongest shielding, the rings B/E/H/J/M/P have moderate shielding, and there is the least shielding in other benzenoid rings (Fig. 3c), consistent with the NICS and ACID calculations. For comparison, although the open-shell diradical form of [4]PA exhibits a localized aromatic character (Fig. 1a), the central benzenoid ring connecting the two tetracene units in[4]PA is almost non-aromatic. The central benzenoid rings connecting the two hexacene/heptacene units in [6]PA and [7]PA are aromatic. This difference could be ascribed to the much larger diradical character of [6]PA and [7]PA compared to [4]PA. Therefore, the lateral extension in [6]PA and [7]PA provides the opportunity to form local aromatic sextets along the central row of the backbone, resulting in different electronic structure.
To further understand the ground-state electronic structure of compond 1, we conducted magnetic measurements such as variable temperature (VT) NMR, electron spin resonance (ESR) and superconducting quantum interference device (SQUID) measurements. The 1H NMR spectrum of 1 in THF-d8 at room temperature exhibited one almost flat baseline, and cooling of sample to -80 oC also did not obtain full resolution (Figure S3). The NMR broadening can be explained by the thermal population of triplet diradical species at elevated temperatures, as expected for species with strong diradical character. Further experimental evidence was obtained by ESR measurement on 1 in solid, showing strong ESR signal at room temperature with a g-factor of 2.003 (Figure S4). Furthermore, we also carried out superconducting quantum interference device (SQUID) measurement for the microcrystalline powder sample of 1 at 2–300 K. The magnetic susceptibility increased with increasing temperature, and the singlet-triplet energy gap (∆ES-T) was estimated to be -1.33 kcal/mol by careful fitting of the data by using the Bleaney–Bowers17 equation (Fig. 4). All these indicate that compound 1 is an open-shell singlet diradicaloid and exhibits a large diradical character with a small singlet-triplet energy gap. For comparison, the [4]PA derivative (-2.5 kcal/mol)6 showed a larger ∆ES-T value, in consistence with its smaller diradical character. On the other hand, broken-symmetry DFT calculations (UB3LYP/6-31G(d,p)) also predict that the open-shell singlet biradical state of 1 has lower energy compared to the triplet biradical state, and the singlet–triplet energy gap (∆ES-T) is calculated to be -2.08 kcal/mol, which is consistent with the above experimental results.
The UV-vis absorption spectra of 1 was recorded in carbon tetrachloride and compared to that of precursor 8 in DCM, as illustrated in Fig. 5a. Compound 8 exhibits a well-resolved absorption maximum at 632 nm along with a small shoulder band peaking at 580 nm. In contrast, the absorption spectrum of 1 displays a broad band in the near-infrared (NIR) region with a maximum at 789 nm, which is significantly redshifted compared with that of precursor 8. Moreover, 1 shows a weak absorption band in the NIR region with maximum (λmax) at 1090 nm along with a shoulder (λsh) at 1197 nm (Fig. 5a). 1 also displays multiple intense absorption bands in the visible region. The NIR absorption band has a similar structure to that in [4]PA derivative6 and [7]PA derivative,7 however, the intensity is much weaker. Furthermore, the characteristic long-wavelength shoulder peak for 1 is strong indication that it has an open-shell singlet ground state, and the band is originated from the ground-state HOMO,HOMO→LUMO,LUMO double excitation.18 The optical band gap (Egopt) of 1 from the onset of its UV/Vis absorption is roughly estimated to be 0.99 eV, which is smaller than that of [4]PA derivative (1.12 eV)6 due to the π-extension.
The electrochemical property of 1 was investigated by differential pulse voltammetry (DPV) and cyclic voltammetry (CV) measurements in dichloromethane (Fig. 5b). Compound 1 displays two reversible oxidation waves with halfwave potentials E1/2ox at 0.04 V and 0.32 V, and two reversible reduction waves with half-wave potentials E1/2red at − 1.03 V and − 1.24 V (vs Fc/Fc+). Compared with [4]PA derivative, the neighboring oxidation or reduction waves of 1 shows a smaller segregation, which can be explained by the larger π-conjugated backbone that reduces the intramolecular Coulomb charge repulsion. In addition, according to the onset potentials of the first oxidation/reduction waves, the HOMO and LUMO energy levels of 1 are estimated to be -4.82 eV and − 3.76 eV, respectively. Thus, the corresponding electrochemical energy gap (EgEC) is estimated to be 1.06 eV, which is consistent with the optical energy gap. On the other hand, compound 1 can be chemically oxidized into its corresponding radical cation and dication with NO•SbF6 in anhydrous DCM. The radical cation exhibits multiple intense absorption bands in the NIR region with λmax at 858 nm, 1210 nm and 1590 nm, respectively (Figure S2), and a strong ESR signal was observed for the radical cation (Figure S4). The dication shows an intense absorption band with λmax at 845 nm in the NIR region (Figure S2). All absorption spectra are in agreement with the time-dependent (TD) DFT calculations (see SI).