Open Access Article
Shinichi
Ogawa
,
Takuya
Kanetomo
* and
Masaya
Enomoto
*
Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan. E-mail: kanetomo@rs.tus.ac.jp
First published on 2nd April 2024
Spirobiacridine diradicals exhibit intramolecular ferromagnetic coupling owing to their unique cruciform structure (D2d molecular symmetry). In this study, novel spirodiradicals 2,2′,7,7′-tetra(tert-butyl)-9,9′(10H,10′H)-spirobisilaacridine-10,10′-dioxyl (1-Si) and a Ge-centred derivative (1-Ge) were synthesized. These compounds exhibited a ground triplet spin state (S = 1), which was confirmed by spectroscopic and magnetic studies. Interestingly, the exchange coupling constant of 1-Ge (+19.3(2) K) was larger than that of 1-Si (+14.2(3) K), although the distance between intramolecular radical centres was increased owing to its larger bonding lengths. This result implies that σ*(Si–Cα)–π* hyperconjugation has a negative effect on spin–spin coupling.
The compound 9,9′(10H,10′H)spirobiacridine-10,10′-dioxyl (SBDO-C) has two orthogonal π systems, resulting in D2d symmetry.11 Each dihydroacridine moiety possesses one unpaired electron, and singly occupied molecular orbitals (SOMOs) are symmetric to the bisecting plane, resulting in the ground triplet state (S = 1). In fact, the spin state of this spirobiacridine diradical is dictated by its tetra(tert-butyl) derivative, 2,2′,7,7′-tetra(tert-butyl)-9,9′(10H,10′H)spirobiacridine-10,10′-dioxyl (1-C).12 While the structural analysis of 1-C recognises the distortion in these acridine moieties, deviating from a strict D2d symmetry, magnetic studies confirm the ground triplet state (2J/kB = +23(1) K in
). Spirobiacridine diradicals, which can manifest the triplet state despite structural distortion, represent an ideal class of substances for exploring molecular designs based on molecular symmetry.13
A spiro centre plays a crucial role in the development of spiro compounds. Group 14 elements such as C, Si and Ge have been intensively studied.14–16 In particular, silicon compounds exhibit a significant reduction in lowest unoccupied molecular orbital (LUMO) energy levels in π-conjugated systems, which is attributed to a robust σ*(Si–C) and π* hyperconjugation effect.17 The effect of spiro-centre elements on intramolecular spin–spin interaction through a spiro-skeleton has not been sufficiently investigated.18 In this study, we synthesized Si and Ge substitutes of 1-C, namely 2,2′,7,7′-tetra(tert-butyl)-9,9′(10H,10′H)spirobisilaacridine-10,10′-dioxyl (1-Si) and 2,2′,7,7′-tetra(tert-butyl)-9,9′(10H,10′H)spirobigermaacridine-10,10′-dioxyl (1-Ge). Based on the theoretical and experimental results, both compounds exhibited ground triplet states. Interestingly, the intramolecular exchange coupling constant of 1-Ge was larger than that of 1-Si, despite the longer distance between the radical centres.
The frozen-solution ESR spectra of 1-Si and 1-Ge in degassed toluene at 100 K are shown in Figs. 1a and b, respectively. The experimental spectra exhibit the sum of zero-field splitting (ZFS) structures based on D2d (D ≠ 0 and E = 0) and D2 (D ≠ 0 and E ≠ 0) symmetries. The small central signal can be assigned to monoradicals as an impurity. These findings provide a good reproduction of the results for 1-C. In addition, a forbidden signal with |Δms| = 2 appeared at half field (Figs. S3 and S4†), indicating the presence of the triplet spin state (S = 1). ZFS parameters for 1-Si and 1-Ge were estimated using EasySpin software.20 For 1-Si, the experimental spectrum was reproduced by simulating D2/D2d/monoradical components with an intensity ratio of 89.6%/6.6%/3.8%. The simulation gave gxx = gyy = 2.0048, gzz = 2.0087, |D| = 7.37 × 10−3 mT and |E| = 7.67 × 10−4 mT for D2 symmetry; gxx = 2.0061, gyy = 2.0036, gzz = 2.0035 and |D| = 7.38 × 10−3 mT for D2d symmetry; and gxx = 2.0046, gyy = 2.0039 and gzz = 2.0147 for monoradical impurity. For 1-Ge, the experimental spectrum was reproduced by simulating three components with the intensity ratio of 91.2%, 5.6% and 3.2%. The ZFS parameters were gxx = 2.0046, gyy = 2.0048, gzz = 2.0081, |D| = 6.64 × 10−3 mT and |E| = 7.83 × 10−4 mT for the D2 symmetry; gxx = gyy = 2.0048, gzz = 2.0019 and |D| = 6.67 × 10−3 mT for D2d symmetry; and gxx = 2.0047, gyy = 2.0045 and gzz = 2.0082 for monoradical impurity. The experimental |D| and g values of D2 symmetry for 1-Si and 1-Ge led to dipole–dipole distances of 7.1 and 7.3 Å, respectively, from the point dipole approximation 2D = 3g2μB/r3. These values reproduced the crystallographically determined distances of 7.4 Å for 1-Si and 7.5 Å for 1-Ge (more details are given below). Compound 1-C had a dipole–dipole distance of 6.9 Å and an experimentally determined distance of 7.2 Å.12
space group (Table S1†). The crystal structure and molecular packing diagram are shown in Fig. 2. Two crystallographically independent molecules were present in a unit cell. The two independent molecules were denoted as 1-Si-a (the Si1 atom) and 1-Si-b (the Si2 atom). The N1–O1, N2–O2, N3–O3 and N4–O4 bond lengths are 1.288(1), 1.289(1), 1.288(1) and 1.288(1) Å, respectively, which are close to the typical values for nitroxides.21 The bond lengths between Si and neighbouring carbon (Cα) atoms are 1.850(1)–1.866(1) Å around the Si1 atom and 1.850(1)–1.868(1) Å around the Si2 atom. These values reproduce those of typical Si–Csp2 bonds.14c,16b,22,23 The Si atom connects two dihydroacridines, forming a cruciform structure with dihedral angles of 80.33° and 82.02° between the two dihydroacridine moieties. In addition, the dihedral angles between the terminal phenyl rings are 28.78° and 6.42° for the N1 and N2 sides, respectively, in 1-Si-a and 5.17° and 38.67° for the N3 and N4 sides, respectively, in 1-Si-b. High-planarity skeletons (small dihedral angles of 6.42 and 5.17°) result from the intermolecular π–π contact between the dihydroacridines, including the N2 and N3 atoms, as shown in Fig. 2a. In this contact, the N2–O2 and N3–O3 sites are positioned above the C41–C46 and C21–C26 sites, respectively. The N2⋯C41 and N3⋯C21 distances are 3.558(1) and 3.546(1) Å, respectively, slightly exceeding the sum of the van der Waals (vdW) radii (C/N: 3.3 Å),24 indicating weak intermolecular contacts. Alternatively, low-planarity skeletons (large dihedral angles of 28.78 and 38.67°) appeared to be related to the intermolecular contacts around the tert-butyl groups (C33 and C53 atoms) and the dihydroacridine moieties (Fig. 2b). The intermolecular O1⋯O4 distance is 3.828(1) Å, exceeding the sum of the vdW radii (O/O: 3.04 Å).24 Nevertheless, this can still contribute to strong intermolecular antiferromagnetic coupling owing to the direct contacts between the radical centres.
The compound 1-Ge crystallized in the monoclinic P21/n space group (Table S1†). The crystal structure and molecular packing diagram are shown in Fig. 3. The N1–O1 and N2–O2 bond lengths are 1.2868(15) and 1.2846(16) Å, respectively, which indicates the presence of a nitroxide group.21 Ge–Cα bond lengths are 1.9345(14)–1.9445(14) Å, which reproduce typical Ge–Csp2 bonds.16b,23,25 The Ge atom connects two dihydroacridines, forming a cruciform structure with a dihedral angle of 84.99°. In addition, the dihedral angles between the terminal phenyl rings are 40.20° and 4.61° for the N1 and N2 sides, respectively. The former is attributed to the contact with the neighbouring tert-butyl groups and 1-Si. The O1 atom is located at 5.653(2) Å from the nearest neighbouring O atom (O1*), suggesting weak intermolecular interaction. Alternatively, a high-planarity skeleton results from the intermolecular π–π contact between the dihydroacridine moieties, as shown in Fig. 3b. The N2–O2 site is positioned above the C21–C22 site, and the O2⋯C21 and N2⋯C22 distances are 3.340(2) and 3.396(2) Å, respectively. These values are close to the sum of the vdW radii (C/O, 3.2 Å; C/N, 3.3 Å),24 indicating strong interactions compared with 1-Si.
The exchange coupling constant 2J/kB of +13.81 K for 1-Ge was larger than that of +12.93 and +11.25 K (a mean value of 12.09 K) for 1-Si. Despite the increase in the distance between the intramolecular radical centres (according to the ESR and structural studies), the observed enhancement in 2J/kB is an intriguing result. However, the MOs for 1-Si-a, 1-Si-b and 1-Ge are complicated owing to the distortion of the spiro structure and the dihydroacridine moieties. Therefore, we confirmed the magnitude of the interaction by the optimized structure for model compounds, SBDO-C and the Si and Ge derivatives SBDO-Si and SBDO-Ge (for details, see ESI, Figs. S8–S10 and Tables S8–S13†). These calculations gave 2J/kB values of +30.37, +18.42 and +21.06 K for SBDO-C, -Si and -Ge, respectively. The values for the Si and Ge derivatives were larger than the above-mentioned constants, implying the suppression of distortion for the dihydroacridine moieties.
To investigate the relationship between exchange couplings and intramolecular spin-centres distances, we calculated second model compounds (SBDO-C′ and -C′′), in which the spiro centres were replaced by a carbon atom in the optimised structural parameters of SBDO-Si and -Ge. The N–N distances for SBDO-C, -C′ and -C′′ were 5.93, 6.44 and 6.52 Å, respectively. The energy gaps between the bonding and antibonding orbitals resulting from spiro-conjugation (ΔEspiro), as shown in Fig. 4a, showed a narrowing trend (Table 1). The exchange coupling constants (2J/kB) for SBDO-C, -C′ and -C′′ were +30.37, +29.00 and +27.05 K, respectively. These values show a proportional change with respect to the distance between the radical centres. These findings imply that spirodiradicals not only rely on molecular symmetry, but also require space/bond interactions within each unit to enhance spin–spin interaction. In this case, the spiro-conjugation of the occupied MOs could play a role in the interaction of the two dihydroacridine radical units.
| α spins | β spins | |||
|---|---|---|---|---|
| ΔEspiro | UMOs/eV | OMOs/eV | UMOs/eV | OMOs/eV |
| SBDO-C | 1.001 | 1.169 | 1.088 | 1.130 |
| SBDO-C′ | 0.438 | 0.429 | 0.461 | 0.423 |
| SBDO-C′′ | 0.380 | 0.358 | 0.401 | 0.354 |
| SBDO-Si | 0.646 | 0.564 | 0.719 | 0.546 |
| SBDO-Ge | 0.491 | 0.430 | 0.524 | 0.420 |
The difference in the interactions of SBDO-Si and -C′ (+18.42 K vs. +29.00 K) showed a more pronounced change compared to that of SBDO-Ge and -C′′ (+21.06 K vs. +27.05 K). For SBDO-C, -Si and -Ge, ΔEspiro values are summarized in Table 1. The tendency of the change in ΔEspiro values with the SBDO-C, -Si and -Ge agrees with the SBDO-C, -C′ and -C′′ cases, and thus the pronounced difference cannot be only explained by it. Therefore, we should consider the contribution of hyperconjugation between the molecular orbitals of the dihydroacridine moieties and the σ*(Si–Cα) orbital. Fig. 5 presents the energy diagrams of the α and β spins in the triplet state for SBDO-C, -Si and -Ge (for more details, see Tables S11–S13†). For SBDO-C, the highest occupied MO (HOMO) and the lowest unoccupied MO (LUMO) were assigned to the bonding and antibonding orbitals based on the spiro-conjugation in the OMO and UMO, respectively. However, for SBDO-Si and -Ge, the HOMOs were the antibonding orbitals from the spiro-conjugation, whereas the LUMOs were the π* orbital of the dihydroacridine units. The energy levels of this orbital for SBDO-Si and -Ge were lower than that of SBDO-C. Here, the LUMOs for SBDO-Si and -Ge have the same phase with lobes at two neighbouring carbon atoms to the spiro centre (Fig. 4b). This symmetry allows the σ*(X–Cα) orbital to interact with the π* orbital, namely σ*(X–Cα)–π* hyperconjugation, as shown in Fig. 4c. This facilitates electron transfer from the spiro centre to the dihydroacridine moieties. In fact, the calculated spin densities of the Si spiro centre of −0.030 was lower than those of the spiro centres for SBDO-C and -Ge, as shown in Table S14.† On the other hand, the O atom for SBDO-Si was +0.474, which was the largest for the three compounds. These results indicate that for SBDO-Si, radical spin is localized at the NO site, resulting in small exchange coupling.
For 1-Si, structural studies revealed two intermolecular contacts, namely indirect O2⋯O3 and direct O1⋯O4 interactions. Notably, to reduce computational costs, the dihydroacridine moiety, which is not involved in the contacts, and the tert-butyl groups were replaced by methyl groups. First, the indirect O2⋯O3 contact through C atoms was evaluated. The calculated spin densities of the T and BS states were mapped onto the molecular skeleton, as shown in Fig. S11.† The triplet and singlet energy states were ET = −2080.30652569 au with 〈S2〉T = 2.0024 and EBS = −2080.30652658 au with 〈S2〉BS = 0.2455. The intermolecular interaction 2J/kB was estimated to be −0.32 K, indicating weak antiferromagnetic coupling. Next, the direct O1⋯O4 contact was evaluated. The calculated spin densities were mapped onto the molecular skeleton, as shown in Fig. S12.† The triplet and singlet energy states were ET = −2080.31124125 au with 〈S2〉 T = 2.0005 and EBS = −2080.31133968 au with 〈S2〉 BS = 0.2159. The intermolecular interaction 2J/kB was estimated to be −34.81 K, indicating strong antiferromagnetic coupling.
For 1-Ge, an indirect O2⋯O2 contact was present, as shown in Fig. 3b. To reduce computational costs, the dihydroacridine moiety, which was not involved in the contacts, and the tert-butyl groups were replaced by methyl groups. The calculated spin densities of the T and BS states were mapped onto the molecular skeleton, as shown in Fig. S13.† The triplet and singlet energy states were ET = −5655.17316304 au with 〈S2〉T = 2.0006 and EBS = −5655.17317585 au with 〈S2〉BS = 0.2520. The intermolecular magnetic interaction 2J/kB was estimated to be −4.62 K, indicating antiferromagnetic coupling.
. This model assumes intramolecular coupling and the direct O1⋯O4 contact. The g value was fixed at 2.0064 from the ESR study. The best-fit curve was achieved with 2J/kB = +14.2(3) K, 2J′/kB = −41.2(2) K and TIP = −0.249(4) × 10−3 cm3 mol−1. These exchange coupling constants almost reproduced the DFT calculation results. The field dependence of the magnetization curve at 2 K is shown in the inset of Fig. 6a. The M value was saturated and reached 0.90 μB at 7 T. This value was lower than the expected M value of 2.0 μB (S = 1 and g = 2.0) and close to the expected 1.0 μB (S = 1/2 and g = 2.0). This result indicates strong intermolecular antiferromagnetic coupling (J′). Namely, in the dimer, described in the 4-spin model based on the O1⋯O4 contact, the spin state reflects the cancellation of two spins. Note the monomer exhibited the ground triplet state, as confirmed by the experimental (2J/kB = +14.2 K) and theoretical results.
![]() | ||
| Fig. 6 Temperature dependence of the product χmT measured at 5 kOe for (a) 1-Si and (b) 1-Ge. Solid lines represent the simulated curves. The inset shows the magnetization curve measured at 2.0 K. | ||
For 1-Ge (Fig. 6b), the χmT value at 300 K was 0.714 cm3 K mol−1, which was close to the expected value of 0.750 cm3 K mol−1. Upon cooling, the χmT value increased abruptly at 252 K, which was a reversible process according to the results of the cooling and second heating measurements (Fig. S15†). This unique behaviour has been further investigated and will be reported in another paper. Upon further cooling, the χmT value increased and reached 0.950 cm3 K mol−1 at 8.5 K, and then decreased to 0.699 cm3 K mol−1 at 2.0 K. The experimental data of 1-Ge were analysed using PHI software.27 The spin Hamiltonian is defined as the linear four S = 1/2 spin model:
. This model assumes intramolecular coupling and the indirect O2⋯O2 contact. The best-fit curve was achieved with 2J/kB = +19.3(2) K, 2J′/kB = −3.69(1) K, g = 2.059(1) and TIP = −0.323(7) × 10−3 cm3 mol−1. These exchange coupling constants almost reproduced the DFT calculation results. The field dependence of the magnetization curve at 2 K is shown in the inset of Fig. 6b. The M value was saturated and reached 2.00 μB at 7 T, which reproduces to be the expected M value (S = 1 and g = 2.0).
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 23349212334922. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4qo00419a |
| This journal is © the Partner Organisations 2024 |