Marcin A.
Majewski
a,
Piotr J.
Chmielewski
a,
Alan
Chien
b,
Yongseok
Hong
c,
Tadeusz
Lis
a,
Maciej
Witwicki
a,
Dongho
Kim
*c,
Paul M.
Zimmerman
*b and
Marcin
Stępień
*a
aWydział Chemii, Uniwersytet Wrocławski, ul. F. Joliot-Curie 14, 50-383 Wrocław, Poland. E-mail: marcin.stepien@chem.uni.wroc.pl; Web: http://www.mstepien.edu.pl
bDepartment of Chemistry, University of Michigan, 930 N. University Ave, Ann Arbor, MI 48109, USA. E-mail: paulzim@umich.edu
cDepartment of Chemistry, Yonsei University, 50 Yonsei-ro, Seoul 120-749, Korea. E-mail: dongho@yonsei.ac.kr
First published on 7th February 2019
A diindenophenanthrene biradicaloid, formally derived from Chichibabin's hydrocarbon, is obtained in a short, scalable synthesis. The present system is electron-rich and devoid of conjugated substituents, and still exhibits very good stability under ambient conditions. The introduction of the diindeno[1,2-a:2′,1′-i] phenanthrene ring framework results in a singlet biradicaloid system with an easily accessible triplet state (ΔES–T = −1.30 kcal mol−1) and a small electronic bandgap (1.39 V). The stability limits of the title hydrocarbon were explored systematically in the solid state, to reveal an unusual thermally initiated hydrogen-scrambling oligomerization process.
Here we show that by a judicious yet simple modification of the classical structure of Chichibabin's hydrocarbon 1, it is possible to create a much more stable biradicaloid system. This effect is achieved by (a) fusion of an extra benzene ring in the biphenyl section of 1, (b) pentannulation of two phenyl substituents, and (c) steric blocking of the remaining two phenyl groups in the form of bulky mesityl substituents. The resulting molecule, 5,10-dimesityldiindeno[1,2-a:2′,1′-i]phenanthrene 2 features a fused ring system (DIPh) isomeric to the DIAn framework of 3. Advantageously, 2 achieves ambient stability in the absence of bulky TIPS-acetylene substituents. This structural simplification results in a pure hydrocarbon system with a significantly lower molecular weight. The type of pentannulation introduced in 2 is also more effective than earlier modifications of 1 in which adjacent phenyl groups were fused to form terminal fluorene subunits. Those systems would either spontaneously polymerize25–27 or required multiple benzo fusion to achieve ambient stability.28
In the solid state, the diindenophenanthrene core of 2 is found to deviate slightly from planarity, apparently as a result of steric repulsions in the two cove-like regions of the fused ring system (Fig. 1). DFT-optimized geometries of 2 are however planar, indicating that the experimentally observed distortion might be promoted by packing forces in the crystal. No π-stacking of DIPh cores was observed in the solid-state structure, in line with the bulky Mes substitution of the fused ring system. The “quinomethane” bonds in 2 (1.395(4) and 1.408(4) Å) are longer than expected of formally double CC bonds, a feature that correlates with the increasing biradicaloid character of diindenoarenes.7 The above bond lengths are similar to the corresponding distances determined for 3 (ref. 5) and 1 (ref. 4) (1.406 and 1.415(3) Å, respectively), in line with the presumed biradicaloid character of 2.
Fig. 1 Molecular structure of 2 determined using X-ray crystallographic analysis. An n-hexane solvent molecule and hydrogen atoms are omitted for clarity. |
In the electronic spectrum of 2, the most intense band (λmax = 600 nm, ε = 4.1 × 104 M−1 cm−1, Fig. 2) is responsible for the deep blue color of the compound, and appears at a longer wavelength than that reported for 1 (λmax = 574 nm).4 While the main absorption in 2 is hypsochromically shifted relative to the corresponding band of 3 (690 nm), the absorption edge is observed at longer wavelengths (1050 nm vs. ca. 900 nm for 3), indicating a reduced optical bandgap of 2. According to a broken-symmetry TD-DFT calculation (ESI†), the first transition, dominated by the HOMO–LUMO excitation, is dipole-forbidden, in accord with the weak intensity of the corresponding experimental absorption.
In cyclic voltammetry measurements (Fig. 3), 2 undergoes three one-electron oxidations (0.26, 0.66, and 1.11 V vs. Fc+/Fc), of which only the first one is fully reversible. Good reversibility was however observed for two reductions (−1.13 and −1.42 V). All of the above redox events are shifted by up to 0.15 V to higher potentials relative to those determined for 3. The electrochemical bandgap of 1.39 V determined for 2 is smaller than the already uniquely narrow gap of 3 (1.45 V). Quinoidal molecules containing isomeric ring fusion-patterns often yield large variations of optical signatures and biradicaloid characters.12,14,32 In particular, the number of Clar sextets in the open-shell valence structure of DIPh is increased relative to DIAn (four in the open-shell formula 2′vs. three in 3′, cf.Scheme 1), and such a change has been demonstrated to enhance open-shell characters of hydrocarbon biradicaloids.33 In the present case, the bandgap change is likely induced not only by the different topologies of the ring systems in 2 and 3, but also by the difference in substitution patterns. TIPS-ethynyl groups were previously shown to reduce the optical bandgap of pentacene by 0.27 eV, an effect attributed to disproportionate lowering of the LUMO energy.34 The influence of substituents on the bandgap of 3 has not been investigated experimentally; however, if the response of 3 to ethynyl substitution is qualitatively similar to that observed in pentacene, the bandgap-reducing effect of replacing anthracene with phenanthrene in 2 would be even greater than indicated by electrochemistry data.
Fig. 3 Cyclic voltammetry (top) and differential pulse voltammetry scans (bottom) obtained for 2 (dichloromethane, tetrabutylammonium hexafluorophosphate as a supporting electrolyte). |
At 220 K, 2 produces a sharp 1H NMR spectrum in the CDCl3 solution, enabling unambiguous assignment of all signals (Fig. 4). When the sample was heated to 360 K, the spectrum showed progressive yet reversible line broadening, accompanied by shortening of T1 relaxation times, observed already at 300 K (Fig. S10 and S11†). These spectral changes were indicative of chemical exchange with a paramagnetic species, which was presumed to be the thermally populated triplet state of 2, in analogy to the reported behavior of 3 and related systems.
Fig. 4 Temperature-dependent 1H NMR spectra of 2 (CDCl3, 600 MHz). The signals were assigned at 220 K using 2D spectroscopic methods (see the ESI†). |
Solid-state samples of 2 yielded a well-defined ESR line (g = 2.00248), however, without resolved triplet components. The lack of fine structure is apparently consistent with the relatively large distance between the two quinomethane carbons in 2 (9.84 Å vs. 10.32 Å in 3 (ref. 5) and 10.05 Å in 1 (ref. 4)). In the point-dipole approximation,35 this distance gives a value of the zero-field splitting (ZFS) parameter |D| = 29 G (0.003 cm−1). This estimate ignores electron delocalization36,37 and can be taken as the lower bound of the actual ZFS interaction. DFT calculations performed for 2 at various levels of theory produced somewhat variable, but noticeably higher estimates of D (50 to 100 G, 0.005 to 0.009 cm−1, see ESI† for computational details). These higher values reflect the effect of spin delocalization but nevertheless suggest that the ZFS in 2 may indeed be at the limit of measurability.38 Changes of the doubly integrated ESR signal intensity recorded in the 143–333 K range could be fitted with the Bleaney–Bowers equation to provide a singlet–triplet gap of ΔES–T = −1.30 kcal mol−1. The measurement of bulk magnetic susceptibility using the Evans method (for a toluene-d8 solution), showed a similar increase of the triplet contribution, yielding a ΔES–T value of −1.19 kcal mol−1, respectively. It was however observed that prolonged heating of the samples above 340 K resulted in further non-reversible increase of the magnetism, which was linked to the thermally induced reactivity of 2, discussed below.
2 can be chemically oxidized to a stable radical cation 2˙+ using [N(p-C6H4Br)3][SbCl6] (magic blue, Scheme 3, Fig. 2 and S20†). Further oxidation was observable when a stronger oxidant, (NO)[SbF6], was employed, but the resulting product was not stable. When treated with an excess of sodium metal in THF-d8, 2 cleanly produced the brownish dianion 22−,6,14 which was extremely unstable to traces of moisture and oxygen. The dianion yielded clearly defined, diamagnetic 1H and 13C NMR spectra, which were analyzed in detail with the help of 2D correlation methods (Fig. 5, S13 and S14†). The 1H resonance pattern of 22− is characterized by significant downfield relocations of the a, b, and c resonances relative to their positions in the spectrum of 2, suggesting that the aromaticity of the phenanthrene core is restored in the dianion. 13C NMR signals of quinomethane carbons in 22− were unambiguously identified at 99.0 ppm on the basis of 1H–13C HMBC correlations. This upfield shift, reproduced in a GIAO calculation (95.4 ppm), is characteristic of π-conjugated carbanions,39 and provides strong evidence for the charged nature of 22−. Interestingly, the bond length pattern obtained for 22− in a DFT calculation is very similar to that observed experimentally in the fluorenide anion,40 indicating a possible similarity of the electronic structures of these two systems. When traces of water were present during the Na-induced reduction of 2, the protonated monoanion (2-H)− was obtained instead of 22−. This new species, identified in situ on the basis of its 1D and 2D 1H NMR spectra (Fig. S12†), could be partly converted back to the neutral 2 by exposure to air.
Scheme 3 Redox reactivity of 2. Reagents and conditions: (a) [N(p-C6H4Br)3][SbCl6] (magic blue). (b) Na metal, THF-d8. (c) Na metal, THF-d8, traces of moisture. (d) Air (partial recovery). |
Fig. 5 1H NMR spectrum (THF-d8, 600 MHz, 300 K, bottom) and DFT-optimized geometry of 22− (PCM(THF)/B3LYP/6-31G(d,p), top). |
Femtosecond transient absorption (fs-TA) measurements performed for 2 (Fig. 6A) revealed ultrafast relaxation dynamics, similar to those observed in other quinodimethane biradicaloids.7,8 Evolution-associated spectra and time traces of population ratios, obtained by global analysis (Fig. 6B and S18†), showed an initial decay with a time constant τ1 = 200 ± 20 fs corresponding to an internal conversion to the S1 state, followed by rapid deactivation to the ground state (τ2 = 7.0 ± 0.7 ps), possibly admixed with vibrational cooling. The short excited-state lifetime of 2 is comparable with those of other diindenoarenes,41 and may be similarly linked to non-radiative decay via an S1/S0 conical intersection.
Fig. 6 fs-TA spectra (A) and evolution-associated spectra (B) of 2 in THF with the excitation of 600 nm. |
The biradicaloid nature of 2 was explored computationally using the spin-flip approach.42 A RAS(2,2)-SF/cc-pVTZ calculation43,44 was performed on the substituent-free DIPh core of 2 optimized at the CAM-B3LYP/6-311++G(d,p) level of theory. A value of the ΔES–T gap of −1.43 kcal mol−1 was obtained, in remarkably good agreement with the results of ESR and Evans measurements. The strong biradicaloid character of the singlet state is indicated by the RAS-SF natural orbital occupancies, which gave a diradical index γ0 of 0.69. In both the S and T states, the odd-electron density, obtained from the RAS-SF calculation, is localized principally at the five-membered rings (Fig. 7B), in line with the open-shell configuration 2′. The validity of 2′ was further verified using 1D and 2D NICS scans performed on the DIPh core (Fig. 7A and S17†), which showed that the aromaticity of the central phenanthrene unit increases in the triplet state. In the singlet state, the terminal six-membered rings of DIPh remain aromatic, whereas the five-membered rings show a weak antiaromatic ring current. This current is apparently responsible for the upfield relocation of the o-Me resonance of 2 (signal i, Fig. 4).
The stability of 2 in solution depends on the choice of solvent and on storage environment (Fig. S1 and S2†). Remarkably, a 20 μM solution of 2 in dichloromethane showed no changes of the absorption spectrum after 25 days under ambient conditions. In THF, DMF, DMSO, and toluene, complete bleaching typically occurred within one day, when the sample was exposed to light. In the dark, however, THF solutions could be stored for over a week without significant decomposition, even under aerobic conditions. Solid 2 can be heated at ca. 70–110 °C for several hours without noticeable degradation. However, at these temperatures, the overall concentration of radicals in the sample, estimated using ESR, is seen to increase irreversibly, indicating a very slow, thermally activated transformation of 2. The onset of decomposition is additionally marked by the appearance of hyperfine structure in the ESR signal, which is particularly well resolved for solution samples. Line shape simulations suggest that this signal may correspond to a monoradical (or a family of structurally related species) formed by addition to one of the quinomethane carbons of 2 (ESI†). Formation of doublet radicals via addition to sp2 centers in biradicals was postulated in earlier work.38,45,46
Above 150 °C, the decomposition of 2 becomes faster, yielding an exotherm in a TGA–DTA scan (Tonset = 178 °C, Fig. S3†), with no significant mass loss. To drive the process to completion in a controlled way, solid 2 was sealed under argon and heated at 240 °C for 14 h. When so treated, 2 produced a greenish, fluorescent material (7, ca. 60% isolated yield), which was accompanied by ca. 5% of recovered 2 (displaced in the reaction tube by sublimation) and, unexpectedly, by ca. 25% of the dihydro precursor 6 (Scheme 4). While indene-based biradicaloids can be prone to hydrogenation,47 the present instance is of interest because of the absence of an external hydrogen source, and the bulky substitution at the five-membered rings. The 1H NMR spectrum of 7 is poorly defined but could be used to confirm the presence of saturated C(H)Mes moieties analogous to those in 6 (Fig. S4†). This structural feature, and the concurrent formation of 6 indicate that considerable scrambling of hydrogens occurs in the above thermally induced transformation.
Scheme 4 Thermally induced oligomerization of 2. Reagents and conditions: (a) argon, pressure tube, 240 °C, 14 h. |
The oligomeric nature of 7 was inferred from 1H DOSY maps and mass spectrometry (Fig. S6, S7 and S32†). Oligomer sizes of up to n = 15 were detectable in linear-mode MALDI analyses. High-resolution mass spectra of the lower oligomers (n ≤ 6) revealed a complex m/z pattern, which corresponded to a mixture of singly charged molecular ions with the general formula [(C46H36)n – xH + yO]+, with the major components characterized by 0 ≤ x ≤ n and y = 0 to 5. Since the oligomerization takes place in the absence of dioxygen, the oxidized species are presumed to form either during workup or in the ion source. The absorption spectrum of 7 (Fig. 2) resembles that of 2, with the absorption onset at ca. 1000 nm, except that the vis-NIR part is broadened and has a diminished intensity relative to the UV part of the spectrum. All these spectral features of 7 support an oligomer structure combining quinoidal and non-quinoidal sections (Scheme 4).
We initially hypothesized that the oligomerization may start with the formation of σ-dimers such as 8 (Scheme 5), which would undergo subsequent hydrogen scrambling to produce the oligomeric products. Analogous additions have been documented for open-shell systems for both inter-48 and intramolecular cases,49–51 with the resulting sp3-linked products undergoing facile dehydrogenation (if possible). Thermal induction of the oligomerization process was at first presumed to indicate that the triplet state of 2 (T-2)is the initial reactive species.52 Using DFT, all possible symmetrical σ-dimers of 2 were optimized (a–a, b–b, etc.), taking into account the different spin states (singlet vs. triplet) and stereoisomers (racemo vs. meso). The e–e dimer 8 was found to have much lower energy than all other structures (26.84 kcal mol−1 for the meso stereoisomer, triplet state), indicating that position e might be the preferred addition site. However, it was also found that the barriers to dimerization are very high on both the broken-symmetry singlet (bsS) and triplet hypersurfaces (ΔG‡ = 39.97 and 40.20 kcal mol−1, respectively), indicating that direct dimerization should be thermally inaccessible.
Looking for a more feasible pathway, we considered addition of the neutral 2 to the radical cation 2˙+, which is expected to form in small amounts in bulk samples of 2. Such an addition would produce the radical-cation intermediate 8˙+ (Scheme 5), which yields broken-symmetry geometries characteristic of a Robin–Day class-I or -II mixed-valence (MV) system.53 The calculated barrier for this process (ΔG‡ = 31.82 kcal mol−1, doublet hypersurface) is significantly lower than those obtained for dimerization of the neutral 2. While the value is still relatively high, it does not bear a quantitative significance because the accuracy of the calculation is potentially limited by the MV character of the system,53 the neglect of solvation and dispersion phenomena, and the choice of the standard state.
The calculated ΔG value for the electron transfer process 8˙+ + 2 → 2˙+ + 8 are low (0 to 2 kcal mol−1) indicating that in principle, the electron can be transferred between chains of various lengths. The oligomerization can thus propagate via a sequence of hydrogen, proton, and electron transfer steps. Importantly, while the initial adduct 8 is high in energy, the formation of its rearranged isomers, such as 9a–b, which contain hydrogenated 5-membered rings, is predicted to be significantly exergonic. A similar stabilization is predicted for 9a˙+ (ΔG of ca. −25 kcal mol−1vs. S-2 + 2˙+). The overall reactivity described above would be distinct from reported cases of biradicaloid oligo- and cyclooligomerizations, in which the addition typically involves sites with highest spin densities, sometimes relatively hindered,26,27,52,54–56 or conjugated acetylene substituents.57 DFT shows that, while the initial barrier to dimerization may be relatively high, the proposed oligomerization is a strongly exergonic process, and should be feasible for neat 2 at high temperatures. Further growth of the oligomer chain can produce a variety of species containing random sequences of closed-shell and open-shell sections. It is worth noting that in such sequences, terminal fluorenyl radicals may form (as seen in the idealized oligomer 11), which can be expected to produce a doublet ESR spectrum, qualitatively similar to that observed for 2 at the beginning of its decomposition (vide supra).
Footnote |
† Electronic supplementary information (ESI) available: Experimental procedures, spectra of new compounds, and details of computational studies. CCDC 1861241. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c9sc00170k |
This journal is © The Royal Society of Chemistry 2019 |