2. Results and Discussion
The syntheses of 1-(3,4-difluorophenyl)- (
1a) and 1-(2,4-difluorophenyl)-3-phenyl-1,4-dihydrobenzo[e][
1,
2,
4]triazin-4-yl (
1b) were carried out via the classic route involving oxidation of amidrazones
3a,b in situ, giving 1,2,4-triazabutadienes
6a,b, which then underwent electrocyclic ring closure thus affording benzotriazines, which were next oxidized to desired 1,2,4-benzotriazinyls
1a,b (
Scheme 1). Initial amidrazones
3a,b were prepared by a reaction of corresponding arylhydrazines
2a,b with
N-phenylbenzimidoyl chloride [
48]. The yields of radicals
1a,b were moderate because the syntheses of amidrazones
3a,b were accompanied by the formation of corresponding backside products since the hydrazino group can attack via either the α or β nitrogen atom.
Analytically pure samples of 1,2,4-benzotriazinyls 1a,b were obtained by column chromatography on silica gel, followed by procedures of recrystallization from a mixture of CH2Cl2 with n-heptane. Both 1,2,4-benzotriazinyls were comprehensively studied in solution and in the solid state.
Determination of electrochemical parameters of the processes of oxidation and reduction of 1,2,4-benzotriazinyls
1a,b was carried out by cyclic voltammetry (CV) on a glassy carbon disc electrode in acetonitrile. It was found that radicals
1a,b are oxidized and reduced chemically and electrochemically reversibly (Figures 1 and S1), which is typical of 1,2,4-benzotriazinyls [
49,
50,
51,
52]. This is evidenced by our analysis of the CV curves of oxidation and reduction obtained at potential sweep rates of 0.05–1.00 V/s (Figures 2, 3, S2 and S3). It was shown that for each process, the ratio of currents in the reverse and forward peaks does not depend on the potential sweep rate and is close to 1.0 (chemical reversibility), and that the value of the interval between the forward and reverse peaks lies in the range of 0.054–0.058 V and differs slightly from theoretical value 0.059 V for electrochemically reversible processes [
53]. It is worth noting that 1,2,4-benzotriazinyl
1b is oxidized and reduced at higher potentials (E
1/2ox = 0.336 V and E
1/2red = –0.766 V) as compared with derivative
1a (0.295 and –0.709 V, respectively).
Figure 1.
CV curves of oxidation and reduction of 1a (3 × 10−3 M) in 0.1 M Bu4NBF4/MeCN on a glassy carbon disk electrode at a potential sweep rate of 0.1 V/s.
Figure 1.
CV curves of oxidation and reduction of 1a (3 × 10−3 M) in 0.1 M Bu4NBF4/MeCN on a glassy carbon disk electrode at a potential sweep rate of 0.1 V/s.
Figure 2.
(left) CV curves of the oxidation of 1a (3 × 10−3 M) in 0.1 M Bu4NBF4/MeCN on a glassy carbon disk electrode at potential sweep rates of 0.05, 0.10, 0.20, 0.50, and 1.00 V/s. (right) Dependences of the potentials of the peaks of oxidation and reciprocal reduction on the current at the peak for the corresponding process.
Figure 2.
(left) CV curves of the oxidation of 1a (3 × 10−3 M) in 0.1 M Bu4NBF4/MeCN on a glassy carbon disk electrode at potential sweep rates of 0.05, 0.10, 0.20, 0.50, and 1.00 V/s. (right) Dependences of the potentials of the peaks of oxidation and reciprocal reduction on the current at the peak for the corresponding process.
Figure 3.
(left) CV curves of the reduction of 1a (3 × 10−3 M) in 0.1 M Bu4NBF4/MeCN on a glassy carbon disc electrode at potential sweep rates of 0.05, 0.10, 0.20, 0.50, and 1.00 V/s. (right) Plots of reduction and reciprocal oxidation peak potentials versus the peak current for the respective process.
Figure 3.
(left) CV curves of the reduction of 1a (3 × 10−3 M) in 0.1 M Bu4NBF4/MeCN on a glassy carbon disc electrode at potential sweep rates of 0.05, 0.10, 0.20, 0.50, and 1.00 V/s. (right) Plots of reduction and reciprocal oxidation peak potentials versus the peak current for the respective process.
Continuous-wave ESR spectra of 1,2,4-benzotriazinyls 1a,b were recorded at room temperature in degassed toluene solutions of various concentrations. Dilute solutions of 1a,b produced EPR spectra typical of Blatter radicals (Figures 4 and S5), which consist of seven-line multiplets resulting from hyperfine interactions between an unpaired electron and three similar but slightly nonequivalent 14N nuclei. The hyperfine coupling (HFC) constants and giso determined by the modeling of the EPR spectra and calculated at the UPBE0/def2-TZVP level are aN(1) = 0.73 / 0.57, aN(2) = 0.49/ 0.41, and aN(4) = 0.51 / 0.45 mT with g = 2.0040 / 2.0044 in 1a and aN(1) = 0.72 / 0.57, aN(2) = 0.49 / 0.41, and aN(4) = 0.51 / 0.45 mT with g = 2.0040 / 2.0044 in 1b. A comparison of the calculated HFC constants and Mulliken spin populations of N atoms indicated that they do not correlate, which can be explained by a slight deviation from the planarity of the triazinyl ring. The deviation leads to a small difference in the hybridization of the N atoms and therefore in populations of their s-type AOs.
Figure 4.
A) The EPR spectrum of radical 1a in a toluene solution (~10–5 M) at 298 K. B) Optimized structure of 1a with the numbering of nitrogen atoms and their Mulliken spin populations.
Figure 4.
A) The EPR spectrum of radical 1a in a toluene solution (~10–5 M) at 298 K. B) Optimized structure of 1a with the numbering of nitrogen atoms and their Mulliken spin populations.
Figures 5 and S6 depict electronic absorption spectra of radicals
1a,b in acetonitrile at room temperature. One can see that both paramagnets
1a,b show rather strong absorption bands in the ultraviolet region with maxima at 268 and 266 nm, respectively. Other absorption bands in the near-UV and visible regions have substantially lower intensity; in
1a, the maxima of the bands are at 314, 369, 425, and 491, and there is a shoulder (sh.) at ~540 nm; in
1b, they are at 314, ~326 (sh.), 367, ~411 (sh.), ~454 (sh.), 486, and ~535 (sh.) nm. On the whole, the observed spectra are characteristic of related 1,2,4-benzotriazinyl radicals.[
54]
Figure 5A shows that the results of the time-dependent DFT calculations are in good agreement with the experiment. According to the calculations, a series of 40 electronic transitions contributes to the UV-Vis spectra of radicals
1a and
1b in the range 17000−43000 cm
−1 (~230–600 nm). At the same time, the near-UV and visible region (17000–34000 cm
−1 or ~295–600 nm) corresponds to excitations in approximately 20 excited states. Thus, the structure of this moderate-intensity absorption region is mainly related to the superposition of absorption bands onto different excited states. Nonetheless, manifestation of vibrational structures of some bands could not be ruled out either.
Figure 5B displays a number of occupied and unoccupied molecular orbitals (MOs) in the ground state of
1a. All MOs related to these diagrams are presented in SI (
Table S1). The letters a–d denote transitions of electrons between higher occupied and lower unoccupied MOs. According to the calculations, the main contribution (92%) to the lowest-energy excited state (
Figure 5A) comes from the configuration with the promotion of the -electron from SOMO to LUMO. The main contribution (82%) to the second excited state comes from the promotion of the -electron from HOMO to the SOMO counterpart (
Figure 5B, excitation b). A similar description was also performed for transitions to the 4th and 6th excited states with the main contributions of promotions c and d (
Figure 5A,B), respectively. Of note, the most intense transition at 275 nm (f = 0.75) matches excitations to the 21st excited state, whose wavefunction is the sum of a large number of configurations, the predominant of which correspond to the promotions of both (26%) and (20%) electrons from the HOMO to the LUMO.
The structures of radicals
1a,b were successfully confirmed by single-crystal XRD.
Figure 6 shows the structure of molecules
1a and
1b as a projection onto the plane of the bicyclic moiety. In structure
1a, two crystallographic molecules
1aA and
1aB are similar in bond lengths but differ in the rotation of the plane of the nonfluorinated phenyl ring relative to the bicycle plane (
Table 1); in
1aA, this angle (–22.8(3)°) is much larger than that in
1aB (5.0(3)°), whereas in
1b, this angle has an intermediate value (10.5(3)°). In
1a, the fluorinated ring forms angles 51–52° with the heterocyclic plane. By contrast, in
1b, it is turned at almost the same angle in the opposite direction (–58.3(2)°).
In radical
1a, the shortest intermolecular contacts in the structures (less than 3.5 Å) are shown in
Figure 7. Alternating centrosymmetric dimers
1aA⋅⋅⋅
1aA and
1aB⋅⋅⋅
1aB form stacks in the structure. Dimers
1aA…
1aA arise due to H-bonds C–H⋅⋅⋅N [C⋅⋅⋅N 3.427(3) Å]. In
1aB⋅⋅⋅
1aB dimers, the shortest distances between atoms of parallel bicyclic parts are 3.351(3) and 3.444(3) Å. Short N⋅⋅⋅N contacts of 3.646(2) Å can be seen between dimers (
Figure 7). In structure
1b, the parallel arrangement of bicyclic parts and phenyl rings favors the appearance of π-stacking interactions between them [the C⋅⋅⋅C and C⋅⋅⋅N distances are 3.413(3) and 3.388(2) Å, respectively], which are complemented by C–H⋅⋅⋅N hydrogen bonds [C⋅⋅⋅N 3.377(3) Å]. Considering these contacts, stacks of molecules can be distinguished in the structure of
1b (
Figure 7).
To estimate the strength of exchange interactions between radicals and to select model magnetic motifs, parameters of exchange interactions within and between stacks were calculated using the spin-unrestricted broken-symmetry approach. This approach has previously been shown to work well for evaluating intermolecular exchange interactions.[
55,
56] In the case of
1a, the stacks consist of two types of radicals,
1aA and
1aB, with three types of contacts (…
1aA…
1aA…
1aB…
1aB…). Therefore, in the calculations for
1a, three considerably different parameters
J were obtained (
Table 2,
J1aA…1aA,
J1aA…1aB, and
J1aB…1aB) at a ratio of 1.00:0.12:0.61. In the case of
1b, the structure consists of uniformly bound stacks with a single parameter,
J1b…1b.
Table 2 indicates that
J1aA…1aA,
J1aB…1aB, and
J1b…1b are of the same order of magnitude, while
J1aA…1aB is much smaller. For both
1a and
1b, our calculations predicted negligible exchange interactions between the radicals of neighboring stacks (
0.2 cm
−1).
The temperature dependences of χT for radicals
1a and
1b are shown in
Figure 8. For
1a, χT at 300 K is 0.23 cm
3mol
−1K and decreases with lowering temperature, reaching a plateau of ~0.013 cm
3mol
−1K below 60 K. The high-temperature value of χT is much lower than the theoretical spin-only value of 0.375 cm
3mol
−1K for the noninteracting radicals at g = 2, and this finding indicates strong antiferromagnetic exchange interactions; the latter result is consistent with our calculations (
Table 2). We simulated the dependence of χT on T by means of a magnetic motif consisting of 12-membered radical stacks rolled into rings. To avoid overparametrization, we fixed the ratio of
J values (1:0.12:0.61) obtained from the calculations. The following parameters correspond to the best fit: g = 2.01 ± 0.05 and
J1aA…1aA/k
B = –292 ± 10 K. Most likely, the nonzero value of χT below 60 K is due to a small amount of monomeric radicals ( ≈ 2.0%).
For
1b, χT at 300 K is also 0.23 cm
3mol
−1K and diminishes with decreasing temperature to 0.009 cm
3mol
−1K at 2 K (
Figure 8). Simulation of the dependence of χT on T with 12-membered looped uniform stacks yielded the following best-fit parameters:
J1b…1b/k
B = –222 ± 17 K, ρ ≈ 2.0%, and g = 2.04 ± 0.01.
Table 2 suggests that the results of BS-DFT calculations are in semiquantitative agreement with the parameters extracted from the simulation of the experimental temperature dependences of the molar magnetic susceptibility. This finding confirms the proposed magnetic motifs of crystals of radicals
1a and
1b.
Given that it was possible to solve crystal structures of difluoro derivatives
1a and
1b, it was reasonable to compare the structures of their solid phases (as well as magnetic properties) with those of nonfluorinated analog
1c (CSD refcode: TICMOH[
57,
58]). According to single-crystal X-ray analysis, radicals
1c are stacked along the
b axis. In the columns, the phenyl ring at C3 and the 1,2,4-triazin-4(1
H)-yl ring lie alternately one on top of the other (
Figure 9). Within the array, the mean interplanar distance is 3.449 Å, which can be considered the mean distance between adjacent π systems. Nevertheless, the centers of the 1,2,4-triazin-4(1
H)-yl rings, bearing most of the spin population, are 5.50 Å apart. Consequently, there are only very weak intermolecular interactions along the columns, in agreement with the results of static magnetic susceptibility measurements (
J1c…1c/k
B = –2.2 ± 0.2 K).[
59]
4. Experimental Part and Computational Details
4.1. Reagents and General Methods
All commercial reagents were used without further purification. Solvents were freshly distilled. All reaction mixtures and column eluents were monitored by thin-layer chromatography (TLC) using commercial aluminum-backed TLC plates (Merck Kieselgel 60 F254). The plates were visualized under UV radiation at 254 nm. Chromatography was performed on silica gel (0.063–0.200 mm) for column chromatography. Melting points were measured by means of Stuart melting point apparatus SMP 30; solvents used for recrystallization are indicated after the melting points.
Measurement of absorption spectra in the UV and visible range was carried out on an Agilent 8453 spectrometer for 10-4 M solutions of compounds in acetonitrile, using a 10 mm quartz cuvette with a Teflon cap. IR spectra of samples in KBr were recorded on a BRUKER Vertex-70 FTIR spectrometer, and strong, medium, and weak peaks are labeled as s, m, and w, respectively. 1H and 13C NMR spectra were acquired on a Bruker Avance 300 spectrometer at 300 and 75 MHz, respectively. Deuterated solvents were used to achieve a homonuclear lock, and the signals are presented in reference to the deuterated solvent peaks. Masses of molecular ions were determined by high-resolution mass spectrometry (HRMS) by means of a DFS Thermo Scientific instrument (EI, 70 eV). Elemental analyses were performed using a Euro EA 3000 elemental analyzer.
4.2. Syntheses
N-phenylbenzimidoyl chloride.[
60] In a round-bottom Schlenk flask,
N-phenylbenzamide (3.0 g, 15.23 mmol) and SOCl
2 (8.85 ml, 8 eq.) were heated at 70 °C. After 4 h, the excess of SOCl
2 was removed to obtain the corresponding imidoyl chloride as a greenish solid (3.15 g, 96% yield), which was subjected to the next step without additional purification. The NMR results are in agreement with previously reported data.
1H NMR (300 MHz, DMSO-d
6) δ 7.06–7.15 (m, 1 H), 7.30–7.41 (m, 2 H), 7.49–7.64 (m, 3 H), 7.78–7.86 (m, 2 H), 7.95–8.03 (m, 2 H).
13C NMR (76 MHz, DMSO-d
6) δ 120.77, 124.07, 128.16, 128.82, 129.03, 131.99, 135.41, 139.59, 165.93.
(3,4-Difluorophenyl)hydrazine. To a solution of 3,4-difluoroaniline (2 g, 15.59 mmol) in a mixture of water (6 ml) and concentrated HCl (12 M, 10 ml), we added dropwise a solution of NaNO2 (1.18 g, 17.04 mmol) in water (2 ml) at 0–5 °C, and the reaction mixture was stirred at this temperature for 1 h. Then, a solution of SnCl2 (8.74 g, 39.73 mmol) in concentrated HCl (12 M, 7 ml) was introduced dropwise at 0–5 °C with vigorous stirring. The mixture was then allowed to warm up to room temperature and incubated for 30 min. The precipitate was filtered off, washed with Et2O (5 ml), and dried in vacuo at 40 °C for 5 h, thus giving (3,4-difluorophenyl)hydrazine hydrochloride as a pale brown solid. Next, the salt was treated with a 40% aqueous NaOH solution (10 ml), and the product was extracted with CH2Cl2 (2 × 5 ml). The combined organic layers were dried over Na2SO4 and concentrated in vacuo to obtain the desired products. Yield 1.69 g (75%), yellow crystals. The product was used at the next reaction step without purification.
(2,4-Difluorophenyl)hydrazine.[
61] To a solution of aniline (2 g, 15.59 mmol) in a mixture of water (6 ml) and concentrated HCl (12 M, 10 ml), a solution of NaNO
2 (1.18 g, 17.04 mmol) in water (2 ml) was added dropwise at 0–5 °C, and the reaction mixture was stirred at this temperature for 1 h. Then, a solution of SnCl
2 (8.74 g, 39.73 mmol) in concentrated HCl (12 M, 7 ml) was introduced dropwise at 0–5 °C with vigorous stirring. The mixture was then warmed to room temperature and incubated for 30 min. The precipitate was filtered off, washed with Et
2O (5 ml), and dried
in vacuo at 40 °C for 5 h, thereby affording (2,4-difluorophenyl)hydrazine hydrochloride as a violet solid. The salt was treated with a 40% aqueous NaOH solution (10 ml), and the mixture was extracted with CH
2Cl
2 (2 × 5 ml) to obtain the title product. Yield 1.69 g (75%), dark yellow crystals, mp 59–61 °C. The product was subjected to the next reaction step without purification.
N′-(3,4-Difluorophenyl)-N′′-phenylbenzohydrazonamide (3a) or N′-(2,4-difluorophenyl)-N′′-phenylbenzohydrazonamide (3b). N-phenylbenzimidoyl chloride (0.3 g, 0.98 mmol) was added to a stirred solution of (3,4-difluorophenyl)hydrazine or (2,4-difluorophenyl)hydrazine (0.478 g, 2.22 mmol) and Et3N (0.336 g, 3.32 mmol) in dry THF (3 ml) at –10 °C. After that, the reaction mixture was stirred at 0 °C for 8 h. The solvent was evaporated under reduced pressure, the residue was treated with a 2% aqueous acetic acid solution (10 ml), and the organic products were extracted with CH2Cl2. The combined extracts were dried with Na2SO4, and the solvent was evaporated. Crude products 3a or 3b were partially purified by flash chromatography (SiO2, hexane/CH2Cl2, 3:1), and due to limited stability, were quickly subjected to subsequent reactions.
1-(3,4-Difluorophenyl)-3-phenyl-1,4-dihydrobenzo[e][
1,
2,
4]
triazin-4-yl (
1a). MnO
2 (0.81 g, 9.32 mmol) was added to a stirred solution of amidrazone
3a (0.3 g, 0.93 mmol) in dry CH
2Cl
2 (4 mL), then the reaction mixture was stirred for 6 h at room temperature. After evaporation of the solvent, the crude product was purified by column chromatography (petroleum ether/ethyl acetate [9:1],
Rf = 0.72). Yield 0.086 g (29%), crystals of dark-brown color, mp 110.2–110.7 °C (from a mixture of CH
2Cl
2 with
n-heptane). IR (KBr): 3447 (m), 3058 (w), 2958 (w), 2927 (w), 2857 (w), 1728 (m), 1613 (m), 1513 (s), 1484 (s), 1450 (m), 1393 (s), 1328 (m), 1268 (m), 1203 (m), 1175 (w) cm
−1. EPR (toluene):
g = 2.0040;
aN(1) = 0.73,
aN(4) = 0.51,
aN(2) = 0.49 mT. HRMS (ESI): m/z [M]
+ calcd for C
19H
12F
2N
3 320.0999, found 320.1002. Found, %: С, 71.37; H, 3.64; F, 11.39; N, 13.07. Calcd. for C
19H
12F
2N
3, %: C, 71.24; H, 3.78; F, 11.86; N, 13.12.
1-(2,4-Difluorophenyl)-3-phenyl-1,4-dihydrobenzo[e][
1,
2,
4]
triazin-4-yl (
1b). MnO
2 (0.81 g, 9.32 mmol) was added to a stirred solution of amidrazone
3b (0.3 g, 0.93 mmol) in dry CH
2Cl
2 (4 mL), then the reaction mixture was stirred for 6 h at room temperature. After evaporation of the solvent, the crude product was purified by column chromatography (petroleum ether/ethyl acetate [9:1],
Rf = 0.68). Yield 0.098 g (33%), swamp-colored crystals, mp 128.6–129.2 °C (from a mixture of CH
2Cl
2 with
n-heptane). IR (KBr): 3471 (m), 3062 (w), 3030 (w), 1947 (w), 1877 (w), 1610 (s), 1508 (s), 1483 (s), 1452 (m), 1431 (m), 1395 (s), 1347 (w), 1269 (m), 1229 (m), 1203 (w), 1141 (m) cm
−1. EPR (toluene):
g = 2.0040;
aN(1) = 0.72,
aN(4) = 0.51,
aN(2) = 0.49 mT. HRMS (ESI): m/z [M]
+ calcd for C
19H
12F
2N
3 320.0999, found 320.0998. Found, %: С, 71.64; H, 3.52; F, 11.41; N, 13.10. Calcd. for C
19H
12F
2N
3, %: C, 71.24; H, 3.78; F, 11.86; N, 13.12.
4.3. EPR Spectroscopy
EPR measurements were carried out on a Jeol JES-FA200 X-band spectrometer using a Jeol X-Band Microwave Unit (9.8 GHz) at 290 K in dilute (down to approximately 10
–5 M) toluene solutions degassed by argon bubbling. The spectra were recorded during one slow (~1 h) scan with a modulation of 0.2 mT at 100 kHz and a power of 4 MW. Isotropic g-factor values were measured experimentally using MgO doped with Mn(II) ions as a standard placed in the resonator simultaneously with the test solution. The spectra were simulated in Winsim v.0.96 software [
62]. Estimated accuracy of determination of HFS constants and of the g value was 0.005 mT and 0.0001, respectively.
4.4. CV
CV measurements were performed in a dry glove box with an argon atmosphere and humidity and oxygen levels not exceeding 1 ppm. The analyzed compounds dissolved in acetonitrile (10−4 М) were subjected to electrochemical transformations in a standard three-electrode glass cell at a potential sweep rate of 0.05–1.00 V⋅s−1. The working electrode was a glassy carbon disc electrode with a diameter of 2.7 mm. Before use, its surface was polished with sandpaper and then by means of chromium (III) oxide paste to achieve a mirror finish. The auxiliary electrode was a platinum wire calcined in a gas burner flame to remove oxides and other possible contaminants from the surface. Potentials of the studied processes were measured relative to a reference electrode, which was a silver wire coated with a layer of silver chloride, separated from the bulk of the electrolyte by an electrolytic bridge filled with a solution of the supporting electrolyte. The reference electrode was calibrated against a ferrocene/ferrocenium couple (E0 = 0.640 V with respect to SHE). The supporting electrolyte was a 0.1 M solution of Bu4NBF4 (99%, Sigma-Aldrich) in acetonitrile (≥99.9%, HPLC gradient grade, Fisher Chemical) with a water content not exceeding 20 ppm, according to Karl Fischer titration; Mettler-Toledo titrator C10SD was employed.
4.5. Single-Crystal X-ray Diffractometry
Intensity data on single crystals of
1 and
2 were collected on a Bruker AXS diffractometer, an Apex Duo (Cu Kα, λ = 1.54178 Å, room temperature). Absorption correction was applied by means of Bruker SADABS (version 2.10) [
63]. The structure was solved by direct methods and refined by the full-matrix least-squares method in an anisotropic approximation for all nonhydrogen atoms. Positions of the H atoms were calculated geometrically and included in the refinement in a riding model. All calculations for structure solution and refinement were performed in SHELXL-2018/3.[
64,
65] Crystallographic data were deposited with the Cambridge Crystallographic Data Centre and can be obtained free of charge via
http://www.ccdc.cam.ac.uk/conts/retrieving.html (or from the CCDC, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: +44 1223 336033; E-mail: deposit@ccdc.cam.ac.uk).
Crystallographic data for 1: C19H12F2N3, M = 320.32, T = 296 K, triclinic space group P-1, a = 10.090(3), b = 13.144(3), c = 13.207(4) Å, α = 75.170(15), β = 89.178(15), γ = 67.551(14)°, V = 1558.0(7) Å3; Z = 4, calc = 1.366 g·cm–3, μ(Cu Kα) = 0.829 mm–1, θ range 3.477–68.110°, Ihkl collected/unique 13490/5451, Rint = 0.0582, 3797 Ihkl with I > 2σ(I), 542 refined parameters, GooF = 1.059, R1 = 0.0458, wR2 = 0.1357, CCDC 2253805.
Crystallographic data for 2: C19H12F2N3, M = 320.32, T = 296 K, monoclinic space group C2/c, a = 24.2137(14), b = 4.3811(3), c = 29.4449(19) Å, β = 106.403(4)°, V = 2996.5(3) Å3; Z = 8, calc = 1.420 g·cm–3, μ(Cu Kα) = 0.863 mm–1, θ range 3.806–66.727°, Ihkl collected/unique 9767/2575, Rint = 0.0464, 1800 Ihkl with I > 2σ(I), 217 refined parameters, GooF = 1.030, R1 = 0.0451, wR2 = 0.1235, CCDC 2253804.
4.6. Magnetic Measurements
Magnetic susceptibility of the polycrystalline samples was measured with a Quantum Design MPMSXL SQUID magnetometer in the temperature range 2–300 K in a magnetic field of up to 5 kOe. Diamagnetic corrections were made using the Pascal constants.
4.7. Computational Details
Parameters of exchange interactions (
) between radicals were calculated via the spin-unrestricted broken-symmetry approach[
66] and the Yamaguchi formula [
67]
where the energies of the high-spin (HS) and low-spin (LS) states were calculated at the UB3LYP/def2-TZVP (with ECP for I)[
68,
69] level for XRD geometry of radical pairs.
Geometries of
1a and
1b optimized at the B97-D3/def2-SVP [
70,
71,
72] level in acetonitrile and toluene were used to calculate the EPR parameters (g
iso and HFC constants with
14N) as well as the energies and oscillator strengths of electronic transitions in the UV-Vis spectra. The solvent was taken into account using the CPCM model [
73]. The EPR parameters were computed at the PBE0/def2-TZVP level.[
74] The UV/Vis spectra of radicals
1a and
1b were calculated via the time-dependent DFT approach[
75] with the B3LYP functional and the def2-TZVP basis set.
All calculations were performed using the ORCA 4.0.1 software package [
76].
To simulate the temperature dependence of the molar magnetic susceptibility, χ(
T), in the form of product χ(
T) × T, we utilized custom-designed software.[
77] For both radicals, the magnetic motifs consist of stacks of exchange-coupled radicals. In the case of
1b, these are uniformly bound stacks. For
1a, the situation is more complicated: the stacks consist of two types of molecules,
1aA and
1aB, with three types of contacts (…
1aA…
1aA…
1aB…
1aB…), and consequently there are three different
J parameters (
J1aA…1aA,
J1aA…1aB, and
J1aB…1aB). To avoid overparametrization, the ratio of these parameters was fixed based on the calculation results (
Table 2).