Jin
Xiong
^{a},
Hai-Yan
Ding
^{b},
Yin-Shan
Meng
^{a},
Chen
Gao
^{a},
Xue-Jing
Zhang
^{c},
Zhao-Sha
Meng
^{a},
Yi-Quan
Zhang
^{b},
Wei
Shi
^{c},
Bing-Wu
Wang
*^{a} and
Song
Gao
*^{a}
^{a}Beijing National Laboratory of Molecular Science, College of Chemistry and Molecular Engineering, State Key Laboratory of Rare Earth Materials Chemistry and Applications, Peking University, Beijing, 100871, P. R. China. E-mail: wangbw@pku.edu.cn; gaosong@pku.edu.cn; Web: http://www.chem.pku.edu.cn/mmm/ Fax: +8610-6275-1708
^{b}Jiangsu Key Laboratory for NSLSCS, School of Physical Science and Technology, Nanjing Normal University, Nanjing 210023, P. R. China
^{c}Department of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (MOE), State Key Laboratory of Elemento-Organic Chemistry, Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Nankai University, Tianjin 300071, China
First published on 3rd October 2016
A hydroxide-bridged centrosymmetric Dy^{III} dimer with each Dy^{III} being five-coordinated has been synthesized using bulky hindered phenolate ligands. Magnetic studies revealed that this compound exhibits a slow magnetic relaxation of a single-ion origin together with a step-like magnetic hysteresis of the magnetic coupled cluster. The thermal relaxation barrier of magnetization is 721 K in the absence of a static magnetic field, while the intramolecular magnetic interaction is very large among reported 4f-only dimers. CASSCF calculations with a larger active space were performed to understand the electronic structure of the compound. The thermal relaxation regime and the quantum tunneling regime are well separated, representing a good model to study the relaxation mechanism of SMMs with intramolecular Dy–Dy magnetic interactions.
Rare examples were reported regarding the synthesis and magnetic studies of Dy^{III} based low coordination number (CN) compounds.^{31–35} Among the low-CN Dy^{III}-based SMMs, {Dy_{5}},^{36} {Dy_{4}K_{2}}^{37} and [Dy(BIPMTMS)(BIPMTMSH)]^{26} molecules (all CN = 6) exhibited high relaxation temperatures and relatively strong Dy^{III}–ligand bonds, suggesting the existence of interesting magneto-structural correlations. One classical method to synthesize low CN Ln compounds is using bulky ligands to increase the steric hindrance.^{30} In addition, with proper hindered ligands, such compounds could accept small hindered ligands. This makes it possible to use a small ligand to bridge such monomers to introduce magnetic interactions and build 4f-only SMMs.
Herein, we present a hydroxide-bridged symmetric five-coordinate Dy^{III} dimer, [Dy(μ-OH)(DBP)_{2}(THF)]_{2} (DBP^{−} = 2,6-di-tert-butylphenolate), which was synthesized by the hydration of low-CN Dy^{III} monomers. This compound exhibits slow magnetic relaxation from a single-ion origin as well as step-like hysteresis curves induced by intramolecular antiferromagnetic Dy–Dy coupling. To the best of our knowledge, this compound represents the highest U_{eff} in lanthanide-only dimers.
Single-crystal X-ray diffraction experiments indicate compounds 1, 2 and 4 are almost isostructural, as they are all crystallized in an orthorhombic space group, Pbca, with close cell parameters (Table S1†). Powder X-ray diffraction experiments confirm that the magnetic diluted sample 3 is a pure single-phase which is isostructural with 2 (Fig. S1†). In 1, each molecule is a centrosymmetric dimer bridged by two hydroxide groups. No H-bonding acceptors can be found for these two hydroxide groups. This is consistent with the FT-IR spectrum (O–H stretching, 3678 cm^{−1}, m, sh, Fig. S2†). Five O atoms coordinate to each Dy^{III} ion, two from the DBP^{−} ligands, two from hydroxide and one from THF (Fig. 1a). The intramolecular Dy–Dy distance is 3.71 Å while the shortest intermolecular Dy–Dy distance is 10.76 Å, suggesting that the intermolecular dipole–dipole interaction could be negligible in comparison to the intramolecular interaction. The crystal field for each Dy^{III} ion exhibits strong anisotropy with four short Dy–O bonds and one longer Dy–O bond (Fig. 1b). The two phenolate O atoms are the most electronegative, leading to very short Ln–O bond lengths of 2.09 Å and 2.12 Å respectively. The two Ln–OH bonds are 2.24 Å and 2.28 Å. The longest Ln–O bond in the dimer is 2.38 Å, in which the O atom is from the electroneutral THF ligand. The bond lengths and angles of Ln_{2}O_{2} parallelograms are listed in Table 1. To evaluate the geometric symmetry of the first coordination sphere, CShM's values (Table S2†) of the Dy^{III} ion in comparison with all of the reference standard 5-coordinate polyhedrons were calculated using the SHAPE 2.0 program.^{41–43} The smallest value is 1.923 for C_{4v} symmetry (square pyramid), indicating no apparent geometric symmetry. (CShM's value is a kind of index to evaluate the geometric similarity between two coordination polyhedrons. During the calculations, the atom type is not considered. Smaller CShM values mean that the two polyhedrons are closer. A value of zero means that the two polyhedrons are geometrically identical).
Fig. 1 Structure details of 1. (a) The molecular structure of 1 with all of the H atoms omitted. (b) The first coordination shell of Dy^{III} in 1. (c) Polyhedrons of the Dy^{III} ions in 1. |
Compound | Bond length (Å) | Bond angles (°) | ||
---|---|---|---|---|
Ln_{1}–O_{4} | Ln_{1}–O′_{4} | O_{4}–Ln_{1}–O′_{4} | Ln_{1}–O_{4}–Ln′_{1} | |
1(Dy^{III}) | 2.28 | 2.24 | 69.5 | 110.5 |
2(Y^{III}) | 2.24 | 2.28 | 68.7 | 111.3 |
4(Gd^{III}) | 2.31 | 2.25 | 69.5 | 110.5 |
To probe the static magnetic behaviour, direct current (dc) magnetic measurements were performed on both of the polycrystalline samples of 1, 3 and 4. The χ_{m}T value of 1 at 300 K was 27.76 cm^{3} mol^{−1} K, which is close to the expected value (28.34 cm^{3} mol^{−1} K) for the two isolated Dy^{III} ions. The sharp decline of the χ_{m}T values at low temperatures may be attributed to three possible reasons: antiferromagnetic coupling, thermal depopulation of low lying crystal field states and magnetic anisotropy. When the χ_{m}T value of 3 was divided by the molar ratio of Dy (ω in eqn (1), 5.7%), the result (χ_{m}T(3′), blue circles in Fig. 2) could be considered as the value of the two isolated Dy^{III} ions with the same coordination environment of 1. Hence the difference defined in eqn (1) can be mainly attributed to the contribution of the intramolecular magnetic interactions. As the insert of Fig. 2 shows, Δχ_{m}T values are negative between 10 K and 16 K, implying an antiferromagnetic interaction between the Dy^{III} ions. The apparent difference between the zero-field-cooled (ZFC) and field-cooled (FC) magnetic susceptibility curves (Fig. S5 and S7†) indicates strong anisotropy as well as open magnetic hysteresis curves for both 1 and 3.
Δχ_{m}T = χ_{m}T(1) − χ_{m}T(3′) = χ_{m}T(1) − χ_{m}T(3)/ω | (1) |
Fig. 2 Plots of the temperature dependence of χ_{m}T(1) (black square) and χ_{m}T(3′) (blue circle). Insert: an enlarged version of the Δχ_{m}T as calculated from eqn (1). All measurements were performed on polycrystalline samples. |
As expected, in the magnetic hysteresis measurements of 1, step-like hysteresis curves were observed up until 8 K under a scanning rate of 200 Oe s^{−1} (Fig. 3). Such steps have been observed on other coupled Dy^{III} dimers, with higher scanning rates or at lower temperatures, before.^{21–23} The loop was kept open near to zero field with a coercive field (H_{c}) of about 2500 Oe at 2 K. Steps appeared at the field near to ±2500 Oe, then the loop was enlarged at higher field until it was finally closed at around ±20000 Oe. These steps could be attributed to a level-crossing between the low lying states of the dimer.^{21–23} Without intramolecular coupling, the hysteresis curve of 3 (Fig. 3 insert) was a typical butterfly-like loop for SIMs with a QTM process near to zero field.
Fig. 3 Plots of the magnetic hysteresis of polycrystalline samples of 1 (main) and 3 (insert). The lines are guides for the eyes. |
Alternative current (ac) magnetic measurements show the dynamic characteristics of SMMs. When 1 was subjected to an ac field from 100 Hz to 10000 Hz under zero static field, the maxima of both the in-phase (χ′_{m}) and out-of-phase (χ′′_{m}) components showed a clear frequency dependence, revealing the typical slow magnetic relaxation behaviour of SMMs (Fig. 4 top and Fig. S14–S16†). Resonance peaks of χ′′_{m} could be observed between 35 K (100 Hz) and 49 K (10000 Hz). The ac susceptibility measurements as a function of the dc field were also carried out. At 2 K, both χ′_{m} and χ′′_{m} remained very small under a dc field smaller than 2000 Oe, then sharply increased to a maximum around the field of 2500 Oe (Fig. S13†). This behaviour coincides with the level-crossing signal detected in the hysteresis measurement. The raised signal in the low temperature regime should be the consequence of the level-crossing process. The slow thermal relaxation behaviour remained the same under a static field of 2500 Oe (Fig. S22–S24†).
Fig. 4 Plots of the temperature dependent out-of-phase ac susceptibility (χ′′_{m}) of 1 (top) and 3 (bottom) at zero static field. The lines are guides for the eyes. |
For comparison, ac magnetic measurements under zero (Fig. 4 bottom and Fig. S30–S32†) and a 2500 Oe dc field (Fig. S38–S40†) were also performed on the polycrystalline samples of 3. As expected, two well-separated relaxation regimes were observed under zero static field (Fig. 4 bottom). In the high temperature regime, frequency-dependent relaxation peaks of χ′′_{m} could be observed between 36 K (100 Hz) and 48 K (10000 Hz) under both zero field and a 2500 Oe dc field. The resonance temperature changes are small after magnetic dilution, implying that the thermal relaxation in the high temperature regime may be a single-ion behaviour. In the low temperature regime, the χ′′_{m} peaks are raised when under zero field and disappear at a 2500 Oe field. This is a typical behaviour of the zero-field QTM process. The QTM process appeared under zero field due to the absence of a bias field that is generated by the intramolecular interaction, and could be suppressed by an external dc field. No resonance peaks were observed in the low temperature regime, as the QTM process is slower than the time limit of the measurement.
The ac measurements at a lower frequency range (1–1000 Hz) were also performed (Fig. S18–S20, S26–S28, S34–S36 and S42–S44†). Although the first coordination sphere of each Dy^{III} has no obvious geometric symmetry, the results of the low-frequency ac susceptibility measurements indicate that, even for the diluted sample 3, the QTM process is slower than 1 Hz (Fig. S35†). This is not usual for the Dy^{III} species, since broken axial symmetry will induce the mixing of the m_{J} states, hence accelerating the QTM process.
The resonance temperatures and extracted relaxation times between the overlapped frequency regimes (100–1000 Hz) of the two measurements (100–10000 Hz on PPMS, 1–1000 Hz on MPMS) were well-matched (see Table S3–S6†). As Fig. 5 shows, the plots of lnτ vs. T^{−1} for the thermal relaxation process are almost the same regardless of magnetic dilution or an applied dc field, confirming that such relaxation is the nature of a single Dy^{III} ion, which is consistent with the conclusions in previously reported work.^{23} Arrhenius fitting of the data of 1 at zero static field (1@ZF) between 40 K and 49 K gave a U_{eff} of 721 K (501 cm^{−1}) and a τ_{0} of 6.6 × 10^{−12} s. The Arrhenius fitting parameters of the other data (1@HF, 3@ZF and 3@HF, HF = 2500 Oe dc field) are listed in Table S7.† The parameters are very close to that of 1@ZF. However, the non-linear trend of the plots implies that a contribution of the Raman process cannot be omitted. All of the plots in the whole measurement regime could be fitted well when one Raman process term was added into the fitting equation. Dual process fitting parameters are listed in Table S8.† For 1@ZF, the dual process fitting gave a U_{eff} of 754 K (524 cm^{−1}) and a τ_{0} of 3.5 × 10^{−12} s. These two parameters are close to those from the Arrhenius fitting, and C is very small (2.7 × 10^{−4} s^{−1} K^{−4.00}), implying that the Orbach process is predominate.
Complete-active-space self-consistent field (CASSCF) calculations on individual Dy^{III} fragments of compound 1 (in which the other Dy^{III} is replaced by a diamagnetic Lu^{III}) on the simplified structure of the X-ray determined geometry have been carried out with MOLCAS 7.8 program packs^{44} (see ESI Section IV for detailed method and results†). However, the results of the conventional method (only choosing the 4f orbitals as the active space, 9 electrons spanning 7 orbitals, CAS(9,7)) could not well explain the effective energy barrier on the basis of a recently reported methodology.^{45} According to the methodology, the rate of the thermal-assisted tunnelling transition (TAT) through the second-excited Kramers doublet (KD3) is 1.7 times faster than that through the first-excited Kramers doublet (KD2) at 50 K. (Approximately, P ∝ exp(−E/kBT)^{2} where P is the transition probability, E is the energy of the doublet relative to the ground state doublet and is the average value of the magnetic moment matrix elements connecting the opposite components of doublets). Hence the TAT process will occur through the second-excited Kramers doublet (KD3 in Table S9†). However, the corresponding energy gap is only 661.8 K (Table S9† and Fig. 6 top), which is smaller than the U_{eff} obtained by the Arrhenius fitting. This is not reasonable since the Arrhenius U_{eff} is usually smaller or equal to the corresponding TAT energy gap. One possible reason for the smaller calculated barrier is the improper choice of active space in the computation. Having only seven f orbitals as the active space in the CASSCF part may not be suitable in treating the present compound, considering the existence of short Ln–O bonds, in which the non-electrostatic interaction may be larger. Thus expansion of the active space was attempted in further calculations, considering the possible contribution of the strong coordinated atoms. Limited by our hardware resources, CAS(11,8) calculations were carried out. The corresponding calculated energy gap was increased to 763 K and the relaxation still occurred through the second-excited state (Table S9† and Fig. 6 bottom). These results were consistent with the Orbach energy barrier that was obtained by the dual relaxation process fitting. It seems that the CAS(11,8) calculation may be more appropriate than the conventional CAS(9,7) calculation to understand the electronic structure of the present system, where the ligands strongly coordinate to Dy^{III}.
Fig. 6 Magnetization blocking barrier in the Dy^{III} fragment in the CAS(9,7) (top) and CAS(11,8) (bottom) cases. The KDs are arranged according to the values of their magnetic moments and noted with the same number in Table S9.† The values and corresponding arrows were marked with the same colours (orange for |i^{−}〉 → |i^{+}〉, green for |i^{−}〉 → |i + 1^{+}〉 and blue for |i^{−}〉 → |i + 1^{−}〉). The red arrows show the relaxation pathway. The purple dashed line is a guide curve of the U_{eff} (501 cm^{−1}). |
Both calculations gave a ground spin state very close to the Ising limit state (Table S10†). The easy axis is close to the perpendicular direction relative to the Dy–O3 (O3 is the O atom of THF, Fig. 1a) bond (Fig. S46†). O3 is the atom which has the least negative Mulliken charge distribution among the five O atoms of the first coordination sphere (Fig. S48†). This demonstrates that when the other ligands are negatively charged, the only one neutral ligand has a decisive role in whether to ‘pin’ the easy axis in the perpendicular direction.^{46}
The intramolecular Dy^{III}–Dy^{III} interaction was also analysed using theoretical and experimental methods (Table 2). The 4f–4f interaction comprises two components, the dipole–dipole interaction and the exchange interaction (coupling constants are represented as J_{dip} and J_{exch}). Ab initio calculations (CAS(11,8)) gave g_{z} = 19.86, g_{x} = 0.002 and g_{y} = 0.002, which are very close to the values of the Ising limit state of Dy^{III} with a pseudospin of ±1/2. Thus an Ising approximation in the Lines model^{47} is reasonable. The total coupling constant can be described as shown in eqn (2) using the Hamiltonian description as shown in eqn (3). The magnetic dipole–dipole contribution to the coupling constant was calculated as J_{dip} = 1.8 cm^{−1}, ferromagnetically (see ESI Section V for the detailed process†).
J = J_{dip} + J_{exch} | (2) |
(3) |
Method | J _{dip} | J _{exch} | J |
---|---|---|---|
Level-crossing | 1.8 | −6.4 | −4.6 |
POLY_ANISO | 1.8 | −5.4 | −3.6 |
Two methods were applied to obtain the magnetic coupling interaction constant. One method is based on the Zeeman effect. Considering that the steps on hysteresis curves are the consequence of level-crossing between the ground state and the first-excited state of the dimer, eqn (4) could be applied to calculate the J value.^{21} The first derivative of the magnetization of 1 reveals that the crossing field could be estimated as being 2.5 kOe (Fig. S10 and S11†). Thus J = −4.6 cm^{−1} was obtained when g = g_{z} = 19.85, hence J_{exch} = −6.4 cm^{−1}. The other method involves fitting the dc magnetic susceptibility of compound 1 during 2 K to 300 K using the program POLY_ANISO,^{48–50} in which all the spin orbitals of each Dy^{III} were considered using the Lines model (see Section V in ESI for details†). The fitting gave a result of J = −3.6 cm^{−1} and J_{exch} = −5.4 cm^{−1}. Both results indicate that the exchange coupling transferred through the two hydroxide bridges is strong among the bridged lanthanide coordination compounds, leading to an apparent antiferromagnetic interaction although the dipolar component is ferromagnetic.
H_{cross} = −J/2gβ | (4) |
For comparison studies, the intramolecular Gd^{III}–Gd^{III} interaction in 4 was also examined using both an analytical method and numerical method (Fig. S50, see ESI Section V for details†). After applying the van Vleck equation with an isotropic spin Hamiltonian, we obtained J_{iso}(Gd^{III}–Gd^{III}) = −0.231 ± 0.000 cm^{−1}, g = 2.01 ± 0.00 (S = 7/2 and R^{2} = 1.000). Numerically, we used the PHI 2.1.6 program^{51} to fit the magnetic susceptibility. The fitting gave J_{iso}(Gd–Gd) = −0.23 cm^{−1} and g = 2.01, revealing the large coupling constants among the reported Gd^{III}-only dimers.^{21,23,35,52} Therefore, the hydroxide-bridged structure is indeed a good candidate for transferring a strong magnetic exchange interaction.
Footnote |
† Electronic supplementary information (ESI) available: Synthetic procedures, structure, magnetism, ab initio calculation, magnetic interaction analysis, Fig. S1–S50, Table S1–S10, the crystal information files of compounds 1, 2 and 4. CCDC 1418657 for 1, 1418658 for 2 and 1506174 for 4. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6sc03621j |
This journal is © The Royal Society of Chemistry 2017 |