Matthias
Reiners
,
Dirk
Baabe
,
Kristoffer
Harms
,
Miyuki
Maekawa
,
Constantin G.
Daniliuc
,
Matthias
Freytag
,
Peter G.
Jones
and
Marc D.
Walter
*
Institut für Anorganische und Analytische Chemie, Technische Universität Braunschweig, Hagenring 30, 38106 Braunschweig, Germany. E-mail: mwalter@tu-bs.de
First published on 14th December 2015
Addition of N-heterocyclic carbenes (L = 1,3-di-tert-butylimidazol-2-ylidene (ItBu), 1,3-di-iso-propyl-4,5-dimethylimidazol-2-yildene (IiPr2Me2), 1,3-mesitylimidazol-2-yildene (IMes) and 1,3-di-(2,6-di-isopropylphenyl)imidazol-2-yildene (IPr)) to the iron half-sandwich complex [Cp′FeI]2 (Cp′ = η5-1,2,4-(Me3C)3C5H2, 1) forms the neutral, 16VE adducts [Cp′FeI(L)] (2–5) in moderate to excellent yields. These complexes were structurally characterised. The NHC ligand binds strongly to the Fe(II) atom, so that no exchange is observed on the NMR and chemical time scale. Fe(II) atoms in the starting material 1 adopt a high-spin configuration (S = 2) and are weakly antiferromagnetically coupled at low temperatures. Furthermore, in contrast to previous reports on related [(η5-C5Me5)FeCl(NHC)] systems, in which the Fe(II) atoms assume an intermediate spin (S = 1), no spin state change occurs upon coordination of the NHC ligand; the Fe(II) atoms in complexes 2–5 retain their high-spin state (S = 2) as shown by solid state magnetic susceptibility and zero-field 57Fe Mössbauer spectroscopy investigations. Density functional theory (DFT) studies at the B3LYP level of theory also agree with a well separated S = 2 ground state for compounds 2–5. Surprisingly for Fe(II) high-spin systems, compounds 1–5 exhibit slow paramagnetic relaxation in their Mössbauer spectra; this can be traced to spin–spin and spin–lattice relaxation processes with unusually large spin–lattice relaxation barriers. A structural model is proposed that associates these processes with crystal packing effects.
The paramagnetic adducts 2–5 were characterized by various spectroscopic techniques and elemental analyses. We were also keen to explore the stability of the NHC-adducts with respect to ligand exchange on the NMR and chemical time scales. To this end, free ItBu was added to C6D6 solutions of 2 and 5 and the 1H NMR spectra were recorded. However, the NMR resonances of 2 and 5 remained unperturbed and sharp resonances attributed to free ItBu were detected in the diamagnetic region. These solutions were then allowed to stand at ambient temperatures for 8 days and the 1H NMR spectra were recorded again, but unchanged spectra indicated that there is also no exchange on the chemical time scale. These observations imply that the NHC ligands bind strongly to the Fe(II) atom, which also raises questions regarding the spin state in these complexes. To address this aspect, the solid state molecular structures might provide some initial insights.
Compound reference | 2 | 3 | 4 | 5 |
---|---|---|---|---|
Chemical formula | C28H49FeIN2 | C28H49FeIN2 | C38H53FeIN2 | C44H65FeIN2 |
Formula mass | 596.44 | 596.44 | 720.57 | 804.73 |
Crystal system | Monoclinic | Monoclinic | Monoclinic | Orthorhombic |
a/Å | 26.809(5) | 9.35317(12) | 9.7691(4) | 18.4790(4) |
b/Å | 15.1557(16) | 18.27596(18) | 38.5982(14) | 16.7584(2) |
c/Å | 16.959(3) | 17.01417(16) | 9.6194(4) | 26.7961(4) |
α/° | 90.00 | 90.00 | 90.00 | 90.00 |
β/° | 120.93(2) | 95.7372(8) | 99.676(4) | 90.00 |
γ/° | 90.00 | 90.00 | 90.00 | 90.00 |
Unit cell volume/Å3 | 5910.8(15) | 2893.80(5) | 3575.6(2) | 8298.2(2) |
Temperature/K | 100(2) | 100(2) | 100(2) | 100(2) |
Space group | C2/c | P21/n | P21/c | Pbca |
No. of formula units per unit cell, Z | 8 | 4 | 4 | 8 |
Radiation type | Cu Kα | Mo Kα | Mo Kα | Cu Kα |
Absorption coefficient, μ/mm−1 | 12.397 | 1.605 | 1.312 | 8.970 |
No. of reflections measured | 39142 | 149287 | 139395 | 101789 |
No. of independent reflections | 5623 | 8415 | 8201 | 8484 |
R int | 0.1274 | 0.0370 | 0.0748 | 0.0607 |
Final R1 values (I > 2σ(I)) | 0.0544 | 0.0199 | 0.0582 | 0.0267 |
Final wR(F2) values (I > 2σ(I)) | 0.1212 | 0.0414 | 0.0957 | 0.0622 |
Final R1 values (all data) | 0.0736 | 0.0251 | 0.0684 | 0.0327 |
Final wR(F2) values (all data) | 0.1304 | 0.0431 | 0.0982 | 0.0650 |
Goodness of fit on F2 | 1.033 | 1.055 | 1.328 | 1.032 |
Δρ/e Å−3 | 2.986/−1.643 | 0.444/−0.328 | 1.299/−1.358 | 0.436/−0.440 |
2 | 3 | 4 | 5 | |
---|---|---|---|---|
Fe–C1 | 2.307(5) | 2.4095(11) | 2.334(4) | 2.3422(18) |
Fe–C2 | 2.327(5) | 2.3424(12) | 2.368(4) | 2.3769(19) |
Fe–C3 | 2.374(5) | 2.2792(12) | 2.357(4) | 2.3710(18) |
Fe–C4 | 2.436(5) | 2.3114(12) | 2.379(4) | 2.3410(18) |
Fe–C5 | 2.345(5) | 2.3660(12) | 2.329(4) | 2.2828(18) |
Fe–C(Cp′) (ave) | 2.3578 ± 0.0501 | 2.3417 ± 0.0500 | 2.3534 ± 0.0215 | 2.343 ± 0.0372 |
Fe–Cp′cent | 2.02 | 2.00 | 2.02 | 2.00 |
Fe–Cp′plane | 2.02 | 2.00 | 2.02 | 2.00 |
Fe–I | 2.7658(9) | 2.71039(19) | 2.7128(6) | 2.6883(3) |
Fe–C18 | 2.151(5) | 2.1244(12) | 2.162(4) | 2.1718(19) |
Cp′cent–Fe–I | 137.76 | 124.26 | 119.59 | 118.81 |
Cp′cent–Fe–C18 | 132.95 | 132.48 | 141.13 | 149.48 |
C18–Fe–I | 91.26(14) | 103.23(3) | 99.07(9) | 91.62(5) |
The molecular structures of several Fe half-sandwich complexes have been reported,10,13 and one important feature of these investigations is the correlation between the Cpcent–Fe distance and the spin state of these molecules. This distance is ca. 1.7 Å in low-spin Fe(II) complexes, such as Cp′2Fe (1.71 Å),10b and increases to ca. 1.9–2.1 Å for Fe(II) compounds with a high-spin configuration, e.g. [Cp′FeI]2 (1.93 Å).10b With increasing spin state, metal–ligand antibonding orbitals are occupied, thus weakening the Fe–C bonds and therefore also increasing the variations in the Cpcent–Fe distances. The coordination spheres around the Fe atom in the NHC-adducts 2–5 may be described as distorted trigonal planar, and the Cpcent–Fe distances range between 2.00–2.02 Å, consistent with a d6 high-spin configuration of the Fe(II) atom. It is noteworthy that for the related Cp* complexes such as [Cp*FeCl(IiPr2Me2)] and [Cp*FeCl(IMes)], the Cpcent–Fe distances are 1.78 and 1.93 Å, respectively.9a Despite this large spread of values, the authors assumed that the Fe(II) atoms in both systems adopt an intermediate spin configuration.9a Considering our observations on adducts 2–5 and assuming that the intermediate spin assignment is indeed correct, the Cp* ligand induces a stronger ligand field than the Cp′ ligand and therefore stabilizes the S = 1 spin configuration, which is not accessible for the Cp′-derived complexes 2–5. The increased electron-donating ability of Cp* compared to Cp′ was also indicated by electrochemical studies on various iron complexes with Cp* and Cp′ ligand systems.14 Another factor contributing to this apparent spin state change may be the significantly increased steric demand of the Cp′ relative to the Cp* ligand,15 which also destabilizes the intermediate spin state. The steric influence of alkyl substitution is also nicely demonstrated in the electronic properties of the corresponding manganocene systems; [Cp*2Mn] (S = 1/2) is a low-spin molecule,16 while [Cp′2Mn] (S = 5/2) adopts the high-spin state.17 The Fe–C(NHC) bond distances in complexes 2–5 are significantly longer than those found in the related [Cp*FeX(NHC)] (X = monoanionic ligand) derivatives,9a but are in the range of those observed for three-coordinate, high-spin Fe–NHC complexes (2.096(2) to 2.184(2) Å), that are not Cp-stabilized such as [(NHC)Fe{N(SiMe3)2}2] (NHC = IMes, IPr),18 [(IPr)Fe{N(SiMe3)2}{SePh}],19 [(IiPr2Me2)FeMes2],20 [(IPr)Fe{NHAr}2] (Ar = C6H3-2,6-Cl2, C6H3-2,6-iPr2),21 [(NHC)Fe(CH2SiMe3)2] and [(NHC)Fe(CH2SiMe3)(Cl)] (NHC = IPr and SIPr),22 [(NHC)Fe(CH2SiMe3)(Cl)] (NHC = IPr, IMes, SIPr, SiMes)23 and [(IiPr2Me2)Fe(σ-CPhCPh2)2].24
Fig. 2 Effective magnetic moment μeffvs. T plots for 1 and NHC-adducts 2–5 (recorded in an applied magnetic field of 0.1 T). |
Oxidation state | Compound | Spin state | δ (mm s−1) | ΔEQ (mm s−1) |
---|---|---|---|---|
a Recorded in solid state on powdered samples at T = 100 K unless otherwise stated. Isomer shifts are specified relative to metallic iron at room temperature. b See ref. 14b (recorded at 77 K). c See ref. 25 (recorded at 4.2 K). d See ref. 25 (recorded at 77 K). e See ref. 26. f See ref. 27. g See ref. 28 (recorded at 80 K). h This work. i See ref. 29 (recorded at 4.2 K). j See ref. 30. | ||||
+I | [Cp′Fe(μ-C10H8)FeCp*]b | S = 0 | 0.63 | 1.74 |
[(η5-C5H5)Fe][(η6-C6Me6)]c | S = 1/2 | 0.90 | 1.54 | |
+II | [(η5-C5H5)Fe][(η6-C6Me6)]+d | S = 0 | 0.45 | 2.00 |
[(η5-C5H5)2Fe]e | S = 0 | 0.53 | 2.45 | |
[(η5-C5H5)Fe(CO)2I]f | S = 0 | 0.23 | 1.83 | |
[Cp*Fe(dppp)][OTf]g | S = 1 | 0.55 | 1.75 | |
[Cp′FeI]2 (1)h | S = 2 | 1.03 | 1.91 | |
[Cp′FeI(ItBu)] (2)h | S = 2 | 0.95 | 3.11 | |
[Cp′FeI(IiPr2Me2)] (3)h | S = 2 | 0.86 | 1.53 | |
[Cp′FeI(IMes)] (4)h | S = 2 | 0.98 | 1.75 | |
[Cp′FeI(IPr)] (5)h | S = 2 | 0.99 | 1.59 | |
+III | [(η5-C5Me5)Fe(dppe)H][PF6]i | S = 1/2 | 0.26 | 0.84 |
[(η5-C5H5)2Fe][BF4]j | S = 1/2 | 0.61 | — | |
[Cp*Fe(dppe)(CO)H][PF6]i | S = 1/2 | 0.30 | 0.72 |
T (K) | δ (mm s−1) | ΔEQ (mm s−1) | H hf (T) | Γ HWHM (mm s−1) | τ c −1 (109 s−1) |
---|---|---|---|---|---|
a Recorded in solid state on powdered samples at various temperatures. Isomer shifts are specified relative to metallic iron at room temperature and were not corrected in terms of the second order Doppler shift. All spectra were analysed by a least-squares fitting routine based on the longitudinal relaxation model developed by Blume and Tjon.35 b The parameter was held constant for the fit at the value given. c Measured on an independently prepared sample. d The T = 300 K spectrum of complex 1 was analysed with a least-squares fit to a Lorentzian doublet, since the fit with the Blume–Tjon relaxation model was only possible with a constant numerical value for the relaxation rate of τc−1 > 50 × 109 s−1 (i.e., the relaxation rate approaches the fast relaxation limit). e The two sub-spectra consist of volume fractions of 35:65 and 63:37 at T = 20 and 30 K, respectively. | |||||
Complex 1 | |||||
4 | 1.116(22) | 2.034(44) | 53b | 0.190(7) | 0.32(13) |
10c | 1.087(22) | 1.962(44) | 53b | 0.195(8) | 0.37(19) |
30 | 1.116(18) | 2.036(34) | 53b | 0.190(6) | 0.35(13) |
40 | 1.097(13) | 2.006(26) | 53b | 0.191(6) | 0.44(15) |
50c | 1.062(10) | 1.933(20) | 53b | 0.175(5) | 0.60(21) |
75 | 1.066(4) | 1.973(7) | 53b | 0.158(3) | 1.38(41) |
100c | 1.027(4) | 1.907(8) | 53b | 0.148(6) | 3.44(30) |
125 | 1.042(2) | 1.970(3) | 53b | 0.151(3) | 7.49(50) |
150c | 1.006(3) | 1.906(6) | 53b | 0.137(5) | 14(3) |
175 | 1.011(2) | 1.954(3) | 53b | 0.139(3) | 24(6) |
200 | 0.994(2) | 1.945(4) | 53b | 0.141(4) | 47(29) |
300d | 0.923(3) | 1.904(6) | (53b) | 0.131(6) | (>50) |
Complex 2 | |||||
20 | 0.966(4) | 3.114(7) | 53b | 0.189(6) | 5.36(62) |
100 | 0.946(4) | 3.108(7) | 53b | 0.176(7) | 16(6) |
200 | 0.896(3) | 3.041(7) | 53b | 0.145(7) | 16(6) |
Complex 3 | |||||
20 | 0.878(4) | 1.332(7) | 50.2(1) | 0.097(26) | 0.0023(3) |
30 | 0.813(45) | 1.436(88) | 50.2b | 0.185(23) | 0.028(1) |
40 | 0.921(11) | 1.640(22) | 50.2b | 0.150(14) | 0.14(1) |
40c | 0.916(65) | 1.610(130) | 50.2b | 0.156(10) | 0.17(1) |
60c | 0.916(7) | 1.615(12) | 50.2b | 0.159(6) | 1.13(6) |
80 | 0.886(3) | 1.559(5) | 50.2b | 0.142(4) | 3.16(18) |
100 | 0.856(3) | 1.532(6) | 50.2b | 0.154(5) | 5.20(51) |
120 | 0.867(2) | 1.528(4) | 50.2b | 0.141(3) | 5.96(46) |
150 | 0.865(2) | 1.557(4) | 50.2b | 0.141(3) | 6.78(59) |
200 | 0.830(4) | 1.543(7) | 50.2b | 0.132(7) | 10(3) |
250 | 0.806(3) | 1.554(5) | 50.2b | 0.140(6) | 18(6) |
Complex 4 | |||||
20e | 0.960(48) | 1.692(88) | 56.1(4) | 0.250(22) | 0.006(3) |
0.861(20) | 1.540(40) | 56.1b | 0.328(50) | 0.051(7) | |
30e | 0.982(94) | 1.634(174) | 56.1b | 0.280(13) | 0.020(4) |
1.031(44) | 1.900(44) | 56.1b | 0.159(49) | 0.080(3) | |
40 | 1.341(200) | 2.480(410) | 56.1b | 0.197(22) | 0.085(6) |
60 | 1.016(50) | 1.890(98) | 56.1b | 0.212(17) | 0.35(3) |
100 | 0.970(8) | 1.790(15) | 56.1b | 0.195(8) | 1.68(11) |
100c | 0.977(1) | 1.705(12) | 56.1b | 0.170(6) | 1.57(8) |
200 | 0.919(7) | 1.673(14) | 56.1b | 0.159(12) | 10(3) |
200c | 0.912(4) | 1.660(8) | 56.1b | 0.149(7) | 7(1) |
300c | 0.856(8) | 1.616(15) | 56.1b | 0.192(15) | 14(8) |
Complex 5 | |||||
15c | 1.001(120) | 1.600(240) | 53b | 0.178(19) | 0.17(2) |
20c | 1.021(88) | 1.628(176) | 53b | 0.183(16) | 0.18(2) |
40 | 1.024(71) | 1.648(144) | 53b | 0.138(9) | 0.16(1) |
60c | 0.984(77) | 1.568(152) | 53b | 0.194(19) | 0.25(3) |
80 | 0.971(17) | 1.582(34) | 53b | 0.146(7) | 0.53(3) |
100 | 0.962(26) | 1.568(50) | 53b | 0.138(17) | 0.97(14) |
100c | 1.008(17) | 1.604(32) | 53b | 0.186(14) | 1.03(10) |
150 | 0.911(5) | 1.484(9) | 53b | 0.138(6) | 3.29(30) |
200c | 0.916(7) | 1.498(14) | 53b | 0.155(12) | 6(2) |
250 | 0.900(3) | 1.471(6) | 53b | 0.139(7) | 13(4) |
Nevertheless, all Mössbauer spectra obtained for compounds 1–5 clearly exhibit the characteristic line shape and temperature-dependent asymmetric line broadening (Fig. 3), indicating the presence of paramagnetic relaxation, which is slow or of the same order of magnitude as the 57Fe nuclear Larmor precession time in the magnetic hyperfine field. For example for complexes 3 and 4 the relaxation time (τc) approaches the experimental time scale τM of the Mössbauer spectroscopy (τM approx. 10−7 to 10−9 s) at a temperature above T = 200 K and reaches the slow relaxation limit at ca. T = 20 K, resulting in a single, magnetically split hyperfine pattern (Fig. 3 and 5). It should be noted that these fully developed six-line Mössbauer spectra are not attributed to a long-range magnetic ordering, which is further substantiated by the observation of Curie–Weiss behaviour for all NHC-adducts 2–5 with small Weiss temperatures in the magnetic susceptibility measurements between T = 2 and 300 K (vide supra).
Fig. 3 Zero-field Mössbauer spectra of compound 3 recorded at T between 40 and 250 K (a) and T between 20 and 40 K (b). |
This allows us to evaluate directly the magnitude of the magnetic hyperfine field at the 57Fe nucleus site (Hhf), which is of the order of Hhf = 50.2(1) and 56.1(4) T at T = 20 K for complex 3 and 4, respectively. The main contributions to the internal magnetic field in zero-field Mössbauer spectroscopy experiments are the Fermi contact (HF), the orbital (HL) and dipolar (HD) term.33,36 In highly ionic compounds, the Fermi contact term is predominantly determined by the spin state S, and HF can be estimated to be ca. 51 T,37 which is in good agreement with the experimentally determined magnetic hyperfine field for compounds 3 and 4. However, since covalence plays a significant role in complexes 1–5, tending to (significantly) reduce the HF value,37a,38 the measured hyperfine field is likely not determined solely by the Fermi contact term; orbital and dipolar terms may also contribute, and applied-field Mössbauer spectroscopy would be useful to quantify these contributions.
The presence of slow paramagnetic relaxation is unusual for Fe(II) high-spin complexes such as the series of compounds presented in this study. For symmetric Fe(II) high-spin complexes, a non-zero orbital momentum combined with spin–orbit coupling (SOC) commonly results in strong orbital-phonon coupling and, therefore, in short correlation times τc relative to the 57Fe nuclear Larmor precession time. Hence, we attribute the observation of slow paramagnetic relaxation for the complexes 1–5 to a more complete quenching of the orbital momentum and the presence of SOC in combination with a (large) negative axial zero-field splitting parameter D. This can be rationalized with a simple crystal field model. Because of the low molecular symmetry (C1) the five-fold degeneracy of the Fe d-orbitals (3d6, 5D) is completely removed and the orbital momentum is thoroughly quenched (Fig. 4).
Fig. 4 Schematic illustration of the crystal field splitting and spin–orbit coupling for the 5D (3d6) term. |
Nevertheless, spin–orbit coupling mixes the wave functions of the orbital ground state and of close-lying excited orbital states and therefore partially restores the orbital momentum and splits the ground state d-orbital singlet into Ms = 0, ±1 and ±2 sublevels, separated by D and 3D, respectively. Large negative D values (combined with a small or vanishing rhombic zero-field splitting parameter E/D) then result in a virtually isolated Ms = ±2 ground state, which is well separated from the Ms = ±1 and 0 state by 3D and 4D, respectively, and gives rise to a highly anisotropic magnetic ground state. This supposition is the starting point for the discussion of the Mössbauer results presented below; and it is worth mentioning that – in general – every Ms sublevel, when thermally populated, will contribute to the experimentally observed Mössbauer spectrum with a different magnetic hyperfine field and quadrupole splitting. The observation of a fully developed magnetically split hyperfine pattern at low temperatures for compounds 3 and 4 (vide supra) is consistent with a large negative D and a virtually isolated Ms = ±2 ground state as proposed above.
Paramagnetic relaxation processes observed in zero-field Mössbauer spectroscopy can be analysed by the stochastic longitudinal relaxation model developed by Blume and Tjon, which assumes a constant orientation of the electric field gradient relative to the direction of the magnetic hyperfine field.35 The magnetic hyperfine field then adopts two values, +Hhf and −Hhf, which corresponds to the virtually isolated Ms = ±2 ground state doublet derived by the proposed simple crystal field approach (vide supra). Although the Blume–Tjon model35 does not consider isotropic paramagnetic relaxation, we found a perfect agreement between theory and experimental data. As an example, we show the spectra of complex 3 for selected temperatures in Fig. 3, and for further comparison, we also plot the low-temperature spectra of complexes 1–5 in Fig. 5.
Fig. 5 Zero-field Mössbauer spectra of complex 1 measured at T = 4 K and NHC-adducts 2–5 recorded at T = 20 K. |
The numerical results of this analysis are summarised in Table 4. The strong correlation of the relaxation rate (τc−1) and the magnitude of the local magnetic hyperfine field also prevents both parameters from being simultaneously fitted when the relaxation rate is too fast. For complexes 3 and 4 the magnetic hyperfine field determined by the measurement at T = 20 K was held constant for the fit at elevated temperatures. Furthermore for complexes 1, 2 and 5 we applied, as an approximation, the mean hyperfine field derived from compounds 3 and 4 of approx. 53 T; this was also fixed for the fit. With the exception of complex 4 at T = 20 and 30 K (vide infra), all measured spectra of compounds 1–5 could consistently be fitted with a single 57Fe site (regardless of the temperature). Furthermore we found no indications of any iron-containing impurity phase. The line widths (ΓHWHM; half-width at half-maximum) at temperatures above the slow relaxation limit (T > 50 K) range between approx. 0.14 and 0.19 mm s−1 and are fairly close to the experimental line width of our spectrometer (ca. 0.13 mm s−1), which further substantiates the overall good sample homogeneity. Nevertheless, for the 20 and 30 K spectra of compound 4 we have to consider an additional sub-spectrum to account for the experimental data, whereas above 30 K the spectra of 4 can be analysed well by a model assuming a single 57Fe site. In the framework of the simple crystal field approach and the energy levels derived in combination with SOC (vide supra), the presence of a second sub-spectrum in the slow relaxation limit may be attributed to thermal occupation of the Ms = ±1 sublevel. Alternatively, we may associate the two sub-spectra with different volume fractions or two inequivalent 57Fe sites that are characterised by crystal packing. When both fractions exhibit different degrees of spin–lattice coupling within the solid state solution, different relaxation barriers and relaxation rates are expected. Unfortunately we are unable to distinguish between these alternative explanations at the current stage. However, we can exclude the possibility that the second sub-spectrum observed at low temperatures results from sample decomposition, since deliberate exposure of 4 to air for ca. five minutes reveals a prominent Fe(III) high-spin Mössbauer doublet at T = 20 K (see ESI‡ for details).
To further quantify the temperature-dependent behaviour of the relaxation rate for compounds 1–5, we considered different relaxation processes, i.e. temperature-independent spin–spin (or quantum tunnelling) processes and temperature-dependent direct one-phonon Orbach and indirect two-phonon Raman type spin–lattice relaxation mechanisms using eqn (1)39
1/τc = 1/τ0 + ATn + B/(exp(Δ/T) − 1) | (1) |
The 1/τ0 term denotes the spin–spin relaxation, while the ATn and B/(exp(Δ/T) − 1) terms describe the Raman and Orbach relaxation processes, respectively. The exponent of the Raman contribution was fixed to n = 7, which is appropriate for a non-Kramers doublet.39 The Orbach relaxation barriers Δ obtained from these analyses are summarised in Table 5, along with the parameters A and B of the simulation of eqn (1); Fig. 6 shows the simulation of the experimental data to the eqn (1).
Fig. 6 Temperature-dependent relaxation rates νc = τc−1 for complexes 1–5. The broken lines consider only the spin–spin and Orbach relaxation term of eqn (1) to illustrate the different contributions from Orbach and Raman type spin–lattice relaxation mechanisms. The full line represents the result of a least-squares fitting routine based on eqn (1). |
Complex | τ 0 −1 (109 s−1) | Δ (K) | A (109 s−1) | B (109 s−1) |
---|---|---|---|---|
a The values were obtained by a least-squares fitting routine based on eqn (1). For compounds 1 and 5, a significant contribution of the Raman relaxation term was found, while for the other adducts this contribution can be neglected. b Since the relaxation rate of adduct 2 is in the vicinity of the fast relaxation limit above T = 20 K, we used an Arrhenius ansatz, i.e. τc−1 = Bexp(−Δ/T), to find a rough estimate for the relaxation barrier Δ based on the data points available. c The parameter was held constant for the fit at the numerical value. | ||||
1 | 0.33(18) | 280(40) | 2.4(1) × 10−15 | 49(16) |
2 | — | <25b,c | — | 19b,c |
3 | <2 × 10−3c | 187(6) | — | 19c |
4 | <4 × 10−3c | 217(7) | — | 16c |
5 | 0.13(9) | 275(44) | 9.4(9) × 10−17 | 14(5) |
The Raman type relaxation is clearly relevant for compounds 1 and 5, whereas the temperature-dependent relaxation of compounds 2, 3 and 4 is sufficiently described solely by an Orbach type relaxation mechanism. Moreover, compound 2 shows relatively fast relaxation in the temperature range investigated. Therefore, the spin–spin and spin–lattice relaxation for this compound cannot be reliably quantified. In contrast, for 1 and 5 the relaxation processes associated with spin–spin interactions occur with relaxation rates in the range between approx. 0.13 × 109 and 0.33 × 109 s−1, while for 3 and 4 the spin–spin contribution with an upper limit of 1/τ0 < 0.002 × 109 and 0.004 × 109 s−1, respectively, is significantly smaller than for the other complexes. The Orbach relaxation barriers for 3 and 4 are of the order of Δ = 200 K (139 cm−1); in contrast adduct 2 shows a significantly smaller relaxation barrier with an upper limit of Δ < 25 K (17 cm−1) and complexes 1 and 5 feature a larger Δ of ca. 280 K (194 cm−1). To the best of our knowledge, the relaxation barriers evaluated in these systems are the largest spin–lattice relaxation barriers so far observed for any Fe(II) high-spin complexes.
Nevertheless, the differences in the relaxation processes observed for compounds 1–5 pose the questions whether these phenomena can be attributed to a structural property of these materials. In all cases the Fe(II) atoms are placed in a low-symmetry environment giving rise to a highly anisotropic magnetic Ms = ±2 ground state (vide supra). With the exception of the some variations for complexes 2–5 in the Fe–C18 and Fe–I bond distances and the Cp′–Fe–I, Cp′–Fe–C18 and I–Fe–C18 angles these molecules are structurally similar, so that the differences in the relaxation behaviour are presumably not only of molecular origin, but may also be associated with a different packing arrangement of these compounds. Within the series 1–5 two classes can clearly be distinguished, that is, those crystallising with four (compounds 3 and 4) or eight Fe(II) atoms (complexes 1,10b2 and 5) in the unit cell (Table 1). More importantly, complexes with four Fe(II) atoms per unit cell reach (or at least approach) the slow relaxation limit below T = 40 K, whereas those with eight Fe(II) atoms/unit cell do not. The packing diagrams of complexes 1–5 are shown in the ESI‡ along with the arrangement and intermolecular Fe⋯Fe distances in the respective complexes. Qualitatively the following trends emerge: (1) the Fe atoms of complexes 3 and 4 are arranged in planes to form rhombuses or zig-zag chains, respectively. Furthermore the intermolecular Fe⋯Fe distances are similar, lying in the range 10.1 to 11.4 Å. This corresponds to the observation of low spin–spin relaxation rates of 1/τ0 < 0.002 × 109 and 0.004 × 109 s−1 for these compounds. In contrast, for complexes with significant higher spin–lattice (compounds 1, 2 and 5) and spin–spin (compounds 1 and 5) relaxation rates, the intermolecular Fe⋯Fe distances vary significantly ranging from 7.7 to 13.4 Å, and for dimer 1 the intramolecular Fe⋯Fe distance is even shorter at 3.53 Å. (2) Remarkably, the more symmetric the arrangement of the Fe atoms within the unit cell of the NHC-adducts 2–5, the higher the Orbach relaxation barrier, reaching ca. 280 K (194 cm−1) for complex 5 in which the Fe atoms form a square prism. Although this is currently a rather crude phenomenological description, further work is in progress to unravel the correlation between slow paramagnetic relaxation and molecular structure and lattice contributions. This may also include applied-field Mössbauer spectroscopy, high-field EPR spectroscopy and AC magnetic susceptibility experiments.
[Cp*FeCl (IiPr2Me2)] | [Cp*FeCl(IMes)] | 3 | 4 | |
---|---|---|---|---|
a Computed at the B3LYP level of theory with the basis sets 6-311G(d,p) for Fe, C, H, N, Cl and SDD for I. Values given in parenthesis refer to the values obtained for B97D using the same basis sets. | ||||
S = 0 | 18.1 [4.3] | 20.3 [4.8] | 21.1 [5.9] | 25.5 [9.9] |
S = 1 | 2.5 [0.0] | 3.6 [0.0] | 3.0 [0.0] | 6.0 [0.0] |
S = 2 | 0.0 [5.2] | 0.0 [3.6] | 0.0 [3.1] | 0.0 [2.5] |
Computational methods in general are associated with uncertainties regarding the absolute values of the relative stabilities of the respective spin states; these are of the order of a few kcal mol−1.41 This aspect becomes especially significant when the energy difference between spin states is small, as in this case between S = 2 and S = 1 (Table 6). Table 6 also shows that B97D overestimates dispersion and non-covalent interactions and therefore it predicts an intermediate spin configuration for [Cp*FeCl(IiPr2Me2)], [Cp*FeCl(IMes)], 3 and 4, which is clearly inconsistent with the experimental data provided above. In contrast, B3LYP, which neglects dispersion effects, is known to (artificially) overstabilize the high-spin state.41,42
So while the electronic ground state might not be accurately determined computationally, DFT methods are known to predict reliably molecular structures of different spin states.43 Therefore a comparison of the computed and experimental structures might be more conclusive (Table 7). A closer inspection reveals that only for [Cp*FeCl(IiPr2Me2)] does the computed Fe–Cpcent distance agree with the experimental data in the intermediate spin configuration (S = 1), whereas for all other compounds this experimental Fe–Cpcent distance is much closer to those computed for the high-spin state.
S = 0 | S = 1 | S = 2 | X-ray data | |
---|---|---|---|---|
a Computed at the B3LYP level of theory with a 6-311G(d,p) basis set for Fe, C, H, N, Cl and a SDD basis set for I. Values given in parenthesis refer to the values obtained for B97D using the same basis sets. b Crystal structure data taken from ref. 9a. | ||||
[Cp*FeCl (IiPr2Me2)]b | ||||
Fe–Cpcent | 1.66 [1.60] | 1.81 [1.75] | 2.00 [1.98] | 1.78 |
Fe–Cl | 2.260 [2.256] | 2.283 [2.275] | 2.300 [2.277] | 2.2434(8) |
Fe–C18 | 1.978 [1.917] | 1.966 [1.890] | 2.140 [2.065] | 1.950(2) |
[Cp*FeCl(IMes)]b | ||||
Fe–Cpcent | 1.67 [1.62] | 1.82 [1.78] | 1.99 [2.00] | 1.93 |
Fe–Cl | 2.262 [2.258] | 2.285 [2.276] | 2.324 [2.310] | 2.2715(7) |
Fe–C18 | 2.006 [1.925] | 1.988 [1.911] | 2.187 [2.045] | 2.085(3) |
Complex 3 | ||||
Fe–Cpcent | 1.70 [1.63] | 1.85 [1.78] | 2.01 [1.96] | 2.00 |
Fe–I | 2.662 [2.628] | 2.693 [2.665] | 2.765 [2.729] | 2.71039(19) |
Fe–C18 | 2.025 [1.954] | 2.000 [1.927] | 2.158 [2.074] | 2.1244(12) |
Complex 4 | ||||
Fe–Cpcent | 1.71 [1.64] | 1.87 [1.80] | 2.02 [1.99] | 2.00 |
Fe–I | 2.693 [2.680] | 2.723 [2.712] | 2.837 [2.822] | 2.6883(3) |
Fe–C18 | 2.060 [1.968] | 2.027 [1.951] | 2.187 [2.086] | 2.1718(19) |
Overall, our DFT computations suggest that the energy difference between the intermediate (S = 1) and high-spin (S = 2) state is small and varies between 2.5–6.0 kcal mol−1 (with B3LYP), so that these 16VE adducts are also likely to adopt the maximum spin state. In addition, the qualitative crystal field splitting scheme of the d-orbitals (Fig. 4) is also found in our DFT analysis (see ESI‡ for details). While there is no ambiguity about the spin states of NHC-adducts 3 and 5 based on our experimental data, the situation differs for [Cp*FeCl(IiPr2Me2)] and [Cp*FeCl(IMes)].9a For the least sterically encumbered derivative [Cp*FeCl(IiPr2Me2)] the computed energy difference between intermediate (S = 1) and high-spin (S = 2) state is the smallest within the series and the computed molecular structure for S = 1 is also very close to the experimental one. All these observations lend some support to the original spin state assignment for [Cp*FeCl(IiPr2Me2)],9a but for [Cp*FeCl(IMes)] the computed high-spin geometry is closer to that of the experimental solid state structure. It appears probable that both spin states need to be considered for these adducts [Cp*FeX(NHC)] (X = monoanionic ligands), which might translate into substantial reactivity differences within this class of molecules.
Footnotes |
† Dedicated to Prof. Manfred Scheer on the occasion of his 60th birthday. |
‡ Electronic supplementary information (ESI) available: Crystallographic data as CIF, crystal packing diagrams, Curie–Weiss plots, UV/Vis and Mössbauer spectra and computational details. CCDC 1432917–1432920. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c5qi00235d |
This journal is © the Partner Organisations 2016 |