Correlation between Slow Magnetic Relaxations and Molecular Structures of Dy(III) Complexes with N5O4 Nona-Coordination

A series of Dy(III) mononuclear complexes [DyA2L] (L denotes Schiff base N5 ligand that occupies equatorial positions and A− denotes bidentate anionic O-donor ligands such as NO3− (1), AcO− (2), and acac− (3)) were synthesized to investigate the correlation between the slow magnetic relaxation phenomena and the coordination structures around Dy(III). The Dy(III) ion in each complex is in a nona-coordination with the anionic O-donor ligand occupying upand down-side positions of the N5 equatorial plane. 2 and 3 show slow magnetic relaxation phenomena under a zero bias-field condition, and all complexes showed slow magnetic relaxation under the applied 1000-Oe bias-field conditions. Arrhenius analyses revealed that the ∆E/kB, the barrier height for magnetization flipping, increases in this order, with the values of 24.1(6), 85(3), and 140(15) K. The effects of the exchanging axial ligands on the magnetic anisotropy were discussed together with the DFT calculations.


Introduction
Single-molecule magnets (SMMs) are fascinating molecule-based nanomaterials, which are characterized by slow relaxation of magnetization at low temperatures [1][2][3][4][5][6][7][8][9][10].Magnetic anisotropy plays an essential role in preventing the magnetization flipping; as the orbital angular momentum of 4f electrons is unquenched in the complex formations, each lanthanide(III) (Ln(III)) ion possesses a large magnetic moment correlated with the total angular momentum J, which is defined by the length of the vector summation of the spin angular momentum S and the orbital angular momentum L. The Dy(III) ion is the most fascinating lanthanide ion due to a large total angular momentum of J = 15/2, accompanied by the Kramers characteristic and an oblate type electronic distribution [11][12][13][14][15][16][17][18][19].The magnetic anisotropy of lanthanide ions is strongly correlated with the electronic repulsion with the crystal field of an appropriate anisotropy; this means an axially stressed crystal field is advantageous for realizing an easy axis anisotropy of the oblate type lanthanide ion.The J ground state splits into 2J + 1 numbered substates with different components along the z axis, i.e., J z .In an appropriate anisotropic crystal field, the pair of substates with the highest Ising character was relatively stabilized compared with those of less Ising type pairs to give an easy axial magnetic anisotropy.To achieve such an anisotropic crystal field, the combination of neutral and anionic ligands, with the former located at equatorial positions and the latter located along the z axis, is simple but powerful [18,[20][21][22][23].This strategy is effective both for light and heavy lanthanide ions with oblate type electronic distributions, such as Ce(III), Nd(III) [24,25], Tb(III), and Dy(III).To this effect, we have previously reported several mononuclear Ln(III) complexes (Ln = Ce, Pr, and Nd) incorporated with neutral and anionic ligand pairs, such as 18-crown-6 and NO 3 − , 1,10-diaza-18-crown-6, and NO 3 − [26], and a single helical N 6 ligand (Figure 1, L ) and NO 3 − [21], respectively.In all these complexes, the central Ln(III) ion is surrounded by an (aza)crown ether or a helical ligand L in an equatorial manner, and axial positions are occupied by two or three nitrate anions; this leads to an easy-axial magnetic anisotropy and slow magnetic relaxation phenomena of Ce(III) and Nd(III) as Kramers ions.The radii of the macrocyclic and helical ligands are slightly larger than the ionic radii of the lanthanide ions, and hence the ligands show distortion [20,26] or helication [27] in the complex formation, hence the coordination distances are rather long and the coordination itself is weak.The mismatching between ligand radius and lanthanide radius becomes larger for heavy lanthanide ions and, for the case of macrocyclic ligands, we have not succeeded in synthesizing heavy lanthanide complexes with macrocyclic ligands to reveal the correlation between slow magnetic relaxation phenomena and the molecular structure accompanied with a series of axial ligands.In this study, we decided to introduce a smaller N 5 ligand L as an equatorial ligand, shown in  state splits into 2J + 1 numbered substates with different components along the z axis, i.e., Jz.In an appropriate anisotropic crystal field, the pair of substates with the highest Ising character was relatively stabilized compared with those of less Ising type pairs to give an easy axial magnetic anisotropy.To achieve such an anisotropic crystal field, the combination of neutral and anionic ligands, with the former located at equatorial positions and the latter located along the z axis, is simple but powerful [18,[20][21][22][23].This strategy is effective both for light and heavy lanthanide ions with oblate type electronic distributions, such as Ce(III), Nd(III) [24,25], Tb(III), and Dy(III).To this effect, we have previously reported several mononuclear Ln(III) complexes (Ln = Ce, Pr, and Nd) incorporated with neutral and anionic ligand pairs, such as 18-crown-6 and NO3 − , 1,10-diaza-18-crown-6, and NO3 − [26], and a single helical N6 ligand (Figure 1, L′) and NO3 − [21], respectively.In all these complexes, the central Ln(III) ion is surrounded by an (aza)crown ether or a helical ligand L′ in an equatorial manner, and axial positions are occupied by two or three nitrate anions; this leads to an easy-axial magnetic anisotropy and slow magnetic relaxation phenomena of Ce(III) and Nd(III) as Kramers ions.The radii of the macrocyclic and helical ligands are slightly larger than the ionic radii of the lanthanide ions, and hence the ligands show distortion [20,26] or helication [27] in the complex formation, hence the coordination distances are rather long and the coordination itself is weak.The mismatching between ligand radius and lanthanide radius becomes larger for heavy lanthanide ions and, for the case of macrocyclic ligands, we have not succeeded in synthesizing heavy lanthanide complexes with macrocyclic ligands to reveal the correlation between slow magnetic relaxation phenomena and the molecular structure accompanied with a series of axial ligands.In this study, we decided to introduce a smaller N5 ligand L as an equatorial ligand, shown in Figure 1.The molecular radius of L is smaller than those of the macrocyclic ligands and the helical ligand L′.Complexes of L with heavy lanthanide ions, such as Tb(III) or Dy(III), will be sufficiently stable to exchange axial ligands with different donating abilities.Moreover, it is expected that L can occupy equatorial positions in a flatter manner than L′, being free from helication.In this study, we synthesized a family of Dy(III) complexes with the general formula of

Synthesis and Characterization
We have newly synthesized three Dy(III) complexes by the reaction of an appropriate dysprosium salt with L in an organic solvent, such as Dy(NO3)3•5H2O in MeCN for 1, Dy(AcO)3•4H2O and Dy(CF3SO3)3 in a 2:1 molar ratio in EtOH for 2, and Dy(acac)3 and Dy(CF3SO3)3 in a 2:1 molar ratio in 2-propanol for 3, to obtain crystals suitable for single X-ray crystallography.Their structures were revealed from the crystallographic analyses both for single crystals and

Synthesis and Characterization
We have newly synthesized three Dy(III) complexes by the reaction of an appropriate dysprosium salt with L in an organic solvent, such as Dy(NO 3 ) 3 •5H 2 O in MeCN for 1, Dy(AcO) 3 •4H 2 O and Dy(CF 3 SO 3 ) 3 in a 2:1 molar ratio in EtOH for 2, and Dy(acac) 3 and Dy(CF 3 SO 3 ) 3 in a 2:1 molar ratio in 2-propanol for 3, to obtain crystals suitable for single X-ray crystallography.Their structures were revealed from the crystallographic analyses both for single crystals and microcrystalline samples (Powder X-ray diffraction (PXRD) data are given as Figure S1, in the Supplementary Information).
Crystal structures of the cationic part of three complexes are shown in Figure 2 and Figures S2-S5.Crystallographic data, accompanied by the selected distances and angles, are summarized in Tables S1  and S2. 1 and 2 crystallized in a triclinic crystal system with a P-1 space group, whereas 3 crystallized in monoclinic P2 1 /n.In the complexes, a neutral L ligates as a pentadentate ligand which occupies equatorial positions forming four pentagonal chelate rings.Two anionic ligands, such as NO 3 − in 1, AcO − in 2, or acac − in 3, occupy the positions above and below the lanthanide ion in a bidentate fashion, to complete the N 5 O 4 nona-coordination of the Dy(III) ion.In all cases, one whole molecule is crystallographically independent; however, each has pseudo two-fold symmetry along the axis passing through Dy and N4 of the central pyridine ring.L shows twist distortion around the pseudo two-fold axis, which can be evaluated from the values of the dihedral angles between the two terminal pyridine rings involving N1 and N7.The estimated angles were 25.45(11)  and it also enforced a slight elongation of the Dy-N distances.The anisotropic ligand field is defined by both the difference in coordination distances and the difference in negative charge distributions of donor atoms.The overall ligand field geometries of 1-3 were similar, and could lead to an axially stressed ligand field by sandwiching the Dy(III) ion between negatively charged O-donor ligands.
The observed structural trend suggests that the magnetic anisotropy is enhanced in the order of 1, 2, and 3 because the difference in the coordination distances become larger in this order.The larger negative charges and the shorter coordination distances of O-donors lead to a larger stress along z-axis, and this is likely to enhance the easy-axial magnetic anisotropy for Dy(III) as an oblate ion.
In 2, crystalline solvent molecules showed heavy disordering.They were assigned as two EtOH and one H 2 O, with 50%, 25%, and 50% occupancies according to the elemental analyses.These solvent molecules were closely located to O3 (~2.85 Å) suggesting the hydrogen bond formation, however, π electrons of the acetate ligand are symmetrically delocalized (O3-C22 = 1.263(3)Å and O4-C22 = 1.265(3)Å), and hence the negative charge distributions of O3 and O4 would be similar.Therefore we assumed that the influence of the hydrogen bond formations on the ligand field anisotropy around Dy(III) ion would be negligible.
The shortest Dy•••Dy distances were estimated as 7.2871(13) Å (symmetry code: 1 -x, 1 -y, 1-z), 7.9634(10) Å (symmetry code: 1 -x, -y, 1 -z), and 9.1817(8) Å (symmetry code: 1 -x, 1 -y, -z).The Dy•••Dy separations were long enough to avoid through-space magnetic interactions.The packing diagrams in the unit cell are given as Figures S2-S4   DC susceptibility data were recorded at various temperatures for all complexes under an application of 1000 Oe DC field.Figure 3 shows the temperature dependence of χ M T products.The Curie constant for Dy(III) ion in an isotropic ligand field was estimated as 14.2 emu K mol −1 (J = 15/2, g = 4/3), and the observed values at 300 K, 13.8, 16.0, and 14.6 emu K mol −1 , respectively, were close to the expected values.This may indicate that the magnetic anisotropy is small in these complexes.The χ M T values were almost constant down to 150 K (for 1) or 100 K (for 2 and 3), and then χ M T gradually decreased as the temperature was cooled, due to the thermal depopulation among the substates arising from the 6 H 15/2 ground state which split under the anisotropic crystal field.
ligand field anisotropy around Dy(III) ion would be negligible.

DC Susceptibility of the Complexes
DC susceptibility data were recorded at various temperatures for all complexes under an application of 1000 Oe DC field.Figure 3 shows the temperature dependence of χMT products.The Curie constant for Dy(III) ion in an isotropic ligand field was estimated as 14.2 emu K mol −1 (J = 15/2, g = 4/3), and the observed values at 300 K, 13.8, 16.0, and 14.6 emu K mol −1 , respectively, were close to the expected values.This may indicate that the magnetic anisotropy is small in these complexes.The χMT values were almost constant down to 150 K (for 1) or 100 K (for 2 and 3), and then χMT gradually decreased as the temperature was cooled, due to the thermal depopulation among the substates arising from the 6 H15/2 ground state which split under the anisotropic crystal field.

AC Susceptibility of the Complexes
The slow magnetic relaxation of the complexes was revealed by measuring the alternating current (AC) magnetic susceptibility under the applied conditions of zero field or direct current (DC) field.For all complexes, slow magnetic relaxation phenomena were initially analyzed under a zero DC bias field (Figure S6), which exhibited no out-of-phase signals χM″ for 1, and weak signals for 2 and 3 owing to the quantum-tunneling magnetization (QTM) relaxation process, which was faster for 1 and 3 than the flipping of the magnetic field.For 2, χM″ shows a peak at ~6000 Hz, which is almost constant in the temperature range up to 6.0 K.This indicates that the magnetization relaxation occurs via QTM in this temperature range.The frequency dependences of χM″ were

AC Susceptibility of the Complexes
The slow magnetic relaxation of the complexes was revealed by measuring the alternating current (AC) magnetic susceptibility under the applied conditions of zero field or direct current (DC) field.For all complexes, slow magnetic relaxation phenomena were initially analyzed under a zero DC bias field (Figure S6), which exhibited no out-of-phase signals χ M " for 1, and weak signals for 2 and 3 owing to the quantum-tunneling magnetization (QTM) relaxation process, which was faster for 1 and 3 than the flipping of the magnetic field.For 2, χ M " shows a peak at ~6000 Hz, which is almost constant in the temperature range up to 6.0 K.This indicates that the magnetization relaxation occurs via QTM in this temperature range.The frequency dependences of χ M " were analyzed on the basis of the Cole-Cole equation (see below), and the QTM relaxation rate τ QTM was estimated as 2.70(9) × 10 −5 s.Upon applying the DC bias field, fast relaxation via the QTM process was suppressed and slow magnetic relaxations were observed for all complexes.To reveal the effects of the bias field on slow magnetic relaxations, the bias field dependences of the out-of-phase susceptibilities were measured under the DC field in the ranges of 0 to 3000 Oe and 0 to 5000 Oe for several temperatures (Figure 4).
analyzed on the basis of the Cole-Cole equation (see below), and the QTM relaxation rate τQTM was estimated as 2.70(9) × 10 −5 s.Upon applying the DC bias field, fast relaxation via the QTM process was suppressed and slow magnetic relaxations were observed for all complexes.To reveal the effects of the bias field on slow magnetic relaxations, the bias field dependences of the out-of-phase susceptibilities were measured under the DC field in the ranges of 0 to 3000 Oe and 0 to 5000 Oe for several temperatures (Figure 4).In all complexes, slow magnetic relaxation occurred when a weak DC field of 200 Oe was applied, and it became slower as the applied field increased.The relaxation rate reached a minimum under a DC field of 1000 to 2000 Oe, and was enhanced again when the field was further strengthened, mainly due to the direct process of relaxation becoming predominant.Hence, the dynamic magnetic property of each complex was revealed under the application of 1000 Oe DC, where the QTM and direct processes were effectively suppressed.Under this condition, all complexes exhibited frequency-dependent in-phase (χM′) and out-of-phase (χM″) susceptibilities at temperatures up to 4.2 K, 7.8 K, and 14 K, for 1, 2, and 3 respectively, with an AC frequency up to 10,000 Hz. Figure 5 shows the frequency dependence of the χM′T products and χM″ values.In all complexes, slow magnetic relaxation occurred when a weak DC field of 200 Oe was applied, and it became slower as the applied field increased.The relaxation rate reached a minimum under a DC field of 1000 to 2000 Oe, and was enhanced again when the field was further strengthened, mainly due to the direct process of relaxation becoming predominant.Hence, the dynamic magnetic property of each complex was revealed under the application of 1000 Oe DC, where the QTM and direct processes were effectively suppressed.Under this condition, all complexes exhibited frequency-dependent in-phase (χ M ) and out-of-phase (χ M ") susceptibilities at temperatures up to 4.2 K, 7.8 K, and 14 K, for 1, 2, and 3 respectively, with an AC frequency up to 10,000 Hz. Figure 5 shows the frequency dependence of the χ M T products and χ M " values.
1 and 2 showed bilaterally symmetric shapes of out-of-phase signals, and hence the AC susceptibility data were analyzed with the Cole-Cole Equation (1) given below [28].
Here χ T and χ S denote the isothermal and the adiabatic susceptibilities, respectively, τ denotes the relaxation time at each temperature, and α denotes the distribution of τ (Tables S3 and S4).At first, the frequency dependence of χ M T products was fitted using these four parameters, and the frequency dependences of χ M " were reproduced using the estimated parameters, which showed good agreement with the observations.The estimated α values were small enough, and hence the slow magnetic relaxations occurred via a single process.3 showed slightly asymmetric signals both for χ M T and χ M ", which obeyed the Cole-Davidson Equation (2) given below [29].
Here β denotes the distribution of τ.The analyses were carried out similarly to those for 1 and 2; the frequency dependences of χ M T products were fitted with Equation (2) first, of which the estimated parameters were listed in Table S5, and then the frequency dependence of χ M " were reproduced.The β values were large at low temperature, with the largest value of 0.53 at 3.0-4.0K; these became smaller at higher temperatures, reaching down to 0.41 at 13 K.Both α and β parameters describe the distributions of τ from the ideal value; however, their scales were slightly different.Figure 6 shows the theoretical calculations of χ" for the different α and β values.In both cases, α = 0 and β = 0 mean that the distribution of relaxation time is zero and all molecules flip with a unique relaxation time.This obeys the Debye equation.The peaks of both plots were normalized as 1.0 at α = 0 and β = 0.When the distribution of τ becomes wider, the number of molecules which flip in the different relaxation time increases, and hence the peak height becomes lower depending on α and β.For the same height of the χ" plots, β takes a value 2.0-2.8 times larger than that of α.Considering this relationship, the β values for 3 correspond to the α values at around 0.2, and hence the distribution of τ is regarded to be sufficiently small.Hence, we assumed 3 to be a field-induced SMM.

+ 𝑖ωτ
Here χT and χS denote the isothermal and the adiabatic susceptibilities, respectively, τ denotes the relaxation time at each temperature, and α denotes the distribution of τ (Tables S3 and S4).At first, the frequency dependence of χM′T products was fitted using these four parameters, and the frequency dependences of χM″ were reproduced using the estimated parameters, which showed good agreement with the observations.The estimated α values were small enough, and hence the slow magnetic relaxations occurred via a single process.
Here β denotes the distribution of τ.The analyses were carried out similarly to those for 1 and 2; the frequency dependences of χM′T products were fitted with Equation ( 2) first, of which the estimated parameters were listed in Table S5, and then the frequency dependence of χM″ were reproduced.The β values were large at low temperature, with the largest value of 0.53 at 3.0-4.0K; these became smaller at higher temperatures, reaching down to 0.41 at 13 K.Both α and β parameters describe the distributions of τ from the ideal value; however, their scales were slightly different.Figure 6 shows the theoretical calculations of χ″ for the different α and β values.In both cases, α = 0 and β = 0 mean that the distribution of relaxation time is zero and all molecules flip with a unique relaxation time.This obeys the Debye equation.The peaks of both plots were normalized as 1.0 at α = 0 and β = 0.When the distribution of τ becomes wider, the number of molecules which flip in the different relaxation time increases, and hence the peak height becomes lower depending on α and β.For the same height of the χ″ plots, β takes a value 2.0-2.8 times larger than that of α.Considering this relationship, the β values for 3 correspond to the α values at around 0.2, and hence the distribution of τ is regarded to be sufficiently small.Hence, we assumed 3 to be a field-induced SMM.Next we examined the temperature dependence of relaxation time τ using two methods.The left column in Figure 7 shows the Arrhenius plots for three complexes measured under a 1000 Oe DC applied field.Bent plots were obtained for all complexes over the entire temperature range.The QTM and direct processes were effectively suppressed under these conditions (Figure 4), and hence the bent shaped plots were due to the presence of Raman and Orbach processes; these two processes dominate in different temperature ranges.The data were first analyzed using the linear Arrhenius equation for the data in high temperature regions, where the Orbach process is predominant.This gave the best fit parameters of ∆E/k B = 19.The right column in Figure 7 shows the temperature dependence of τ in double logarithm plots.It is known that the Raman process dominates in low temperature regions and that it obeys Equation (3) below.In plots of ln(τ) vs. ln(T), relaxation via the Raman process is easily identified as a linear region found at low temperatures.The data here were analyzed using Equation (3), which gave the best fit parameters of C = 96(15) s −1 K −4.0 and n = 4.0(2) for 1, C = 0.17(3) s −1 K −5.42 and n = 5.42(13) for 2, and C = 0.048(4) s −1 K −5.12 and n = 5.12(4) for 3.
Finally, the whole dataset was analyzed using Equation ( 4) which considers both Raman and Orbach processes [30,31].
Magnetochemistry 2019, 5, x FOR PEER REVIEW 9 of 15 Next we examined the temperature dependence of relaxation time τ using two methods.The left column in Figure 7 shows the Arrhenius plots for three complexes measured under a 1000 Oe DC applied field.Bent plots were obtained for all complexes over the entire temperature range.The QTM and direct processes were effectively suppressed under these conditions (Figure 4), and hence the bent shaped plots were due to the presence of Raman and Orbach processes; these two processes dominate in different temperature ranges.The data were first analyzed using the linear Arrhenius equation for the data in high temperature regions, where the Orbach process is predominant.This gave the best fit parameters of ΔE/kB The right column in Figure 7 shows the temperature dependence of τ in double logarithm plots.It is known that the Raman process dominates in low temperature regions and that it obeys Equation (3) below.In plots of ln(τ) vs. ln(T), relaxation via the Raman process is easily identified as a linear The estimated values are summarized in Table 1.Some parameters differ from the values based on the simple Raman or Orbach process.This may be due to the narrow range of temperatures where slow magnetic relaxation was observed.From the structural analyses, the crystal field anisotropy as well as the magnetic anisotropy is enhanced in the order of 1 to 3, and the kinetic parameters in Table 1 clearly support the idea.This result indicates that magnetic relaxation phenomena are strongly dependent on the nature of axial ligands, and that the increase of negative charge distributions on O-donor atoms plays a primitive role in enhancing the magnetic anisotropy.However, the structural changes from complex 1 with NO 3 − axial ligands to 2 with AcO − axial ligands were small, while the increase of ∆E/k B values was rather drastic.This may suggest the presence of other factors enhancing the barrier height of the magnetization flipping; hence we examined the electronic structures of the complexes with the DFT technique shown below.

Molecular Orbital Analyses
To elucidate the electronic structures of these complexes, we performed density functional theory (DFT) calculations for three model structures m1-m3, pertaining to complexes 1-3, respectively, as illustrated in Figure 8. Calculated atomic spin densities of the Dy(III) ions in each model summarized in Table 2 are confirmed to be ~5.0,showing that the appropriate electronic configurations are obtained.Kohn-Sham orbital shapes and the energies of f orbitals are summarized in Figure S7 in the Supplementary Information, together with orbital energies of the highest occupied molecular orbital (HOMO).The energy level of occupied f electrons is approximately −10 to −13 eV, much lower than those of HOMO; this suggests that those are independent of the effects of frontier orbitals.On the other hand, those energy levels are comparable to the levels of coordinating σ orbitals and the stable π orbitals of ligands.As a consequence, some f orbitals can be hybridized with those ligand orbitals, although it is usually assumed that the f orbitals do not interact with ligand orbitals.As depicted in Figure S7, the f orbitals of m1 are almost localized, except for HOMO-39(a) and HOMO-46(a) that are also found in equatorial ligand, while those of m2 and m3 are more delocalized by hybridization especially with σ and π orbitals of axial AcO − and acac − ligands, respectively.In addition, the atomic spin densities of the Dy(III) ions in m2 and m3 are slightly smaller than those in m1.The result also indicates that the f spins are slightly delocalized to the ligands in m2 and m3.Hence, the f electrons are considered to interact with axial ligand orbitals in m2 and m3 but not in m1, suggesting that the magnetic anisotropy of m2 and m3 can be affected by such hybridization between f and ligand orbitals.

General Procedures and Methods
All chemicals and reagents were of reagent grade and used without further purification.All chemical reactions and sample preparations for physical measurements were performed in an ambient atmosphere.Variable temperature magnetic susceptibility measurements were performed on PPMS-9 and MPMS-XL magnetometers (Quantum Design).Diamagnetic corrections for each sample were applied using Pascal's constants.Elemental analyses were carried out with the help of the Research and Analytical Centre for Giant Molecules, Graduate School of Science, Tohoku University.

Synthesis of Complexes
Synthesis of L: The ligand L was prepared according to a previously reported method [32].10 mmol of 2,6-diacetylpyridine and 20 mmol of 2-pyridylhydrazine were dissolved into 25 mL of EtOH and one drop of conc.HCl solution was added.The resulting solution was stirred at 60 • C for 1 h to give a pale-yellow precipitate.The solution was left to stand for one night at room temperature to complete the reaction, and then the precipitate was filtrated under reduced pressure, washed with EtOH, and dried in vacuo.Yield 3.01 g (8.7 mmol, 87%).

Crystallography
A single crystal of each complex was mounted on a Rigaku Varimax Saturn area detector for data collection using confocal monochromated MoKα radiation at low temperature (153 K).Intensity data were corrected for absorption using an empirical method included in the Crystal Clear software [33].

Figure 1 .
The molecular radius of L is smaller than those of the macrocyclic ligands and the helical ligand L .Complexes of L with heavy lanthanide ions, such as Tb(III) or Dy(III), will be sufficiently stable to exchange axial ligands with different donating abilities.Moreover, it is expected that L can occupy equatorial positions in a flatter manner than L , being free from helication.In this study, we synthesized a family of Dy(III) complexes with the general formula of [DyA 2 L] + (A − = bidentate anionic O-donor ligand), such as [Dy(NO 3 ) 2 L]NO 3 (1), [Dy(AcO) 2 L]CF 3 SO 3 (2), and [Dy(acac) 2 L]CF 3 SO 3 (3), to investigate the correlation between the axial ligand nature and slow magnetic relaxation phenomena.Syntheses, crystal structures, magnetic properties, and DFT calculation results of these complexes are discussed.

Figure 1 .
Figure 1.Structure of N 6 ligand L and N 5 ligand L.

( 2 ) Å for 3 .
The longest Dy-N distances in 3 resulted from the twisting of L. The coordination distances of O-donors and N-donors in 1 were almost similar, whereas in 2, Dy-O distances were slightly shorter than those of Dy-N.The molecular structures of 1 and 2 were very similar, with twisting of L, and the difference in coordination distances reflects the difference in negative charge distributions between nitrate and acetate O-donor atoms.Because a negative charge in the nitrate delocalized among three oxygen atoms, the charge density of O-donor atoms in the nitrate is slightly smaller than those of the acetate anion.Larger negative charges in the O-donor enforced a slight shortening of the coordination distance in 2, .

Figure 2 .
Figure 2. ORTEP drawing of the cationic parts of 1 (top), 2 (central), and 3 (bottom) at the 50% probability level.Hydrogen atoms are omitted for clarity: (a,c,e) Top view and (b,d,f) side view of the molecules.

Figure 3 .
Figure 3. Temperature dependence of the χMT product for each complex measured under a 1000 Oe DC field applied condition.

Figure 3 .
Figure 3. Temperature dependence of the χ M T product for each complex measured under a 1000 Oe DC field applied condition.

Figure 4 .
Figure 4. DC field dependence of the relaxation rate τ −1 of (a) 1, (b) 2, and (c) 3 measured in the temperature ranges of 2 to 4 K, 4 to 7 K, and 6 to 12 K, respectively.

Figure 4 .
Figure 4. DC field dependence of the relaxation rate τ −1 of (a) 1, (b) 2, and (c) 3 measured in the temperature ranges of 2 to 4 K, 4 to 7 K, and 6 to 12 K, respectively.

Figure 5 .
Figure 5. AC susceptibility data of (a) 1, (b) 2, and (c) 3 measured under an application of 1000 Oe bias field: frequency dependence of the products of temperature and in-phase susceptibility (closed circles) and out-of-phase susceptibility (open circles) measured at several temperatures.Solid curves represent theoretical calculations on the basis of the Cole-Cole equation for 1 and 2, or the Cole-Davidson equation for 3, of which the estimated parameters are listed in Tables S3-S5.

Figure 6 .
Figure 6.Theoretical curves of the out-of-phase susceptibility χ″ as a function of ωτ for several α values (left) or β values (right), calculated on the basis of the Cole-Cole and Cole-Davidson equations.At α = 0 and β = 0, both plots obey Debye equation.The value of χ″ was normalized as 1.0 at ωτ = 0 and α = 0 for the left plot and at ωτ = 0 and = 0 for the right plot.

Figure 6 .
Figure 6.Theoretical curves of the out-of-phase susceptibility χ" as a function of ωτ for several α values (left) or β values (right), calculated on the basis of the Cole-Cole and Cole-Davidson equations.At α = 0 and β = 0, both plots obey Debye equation.The value of χ" was normalized as 1.0 at ωτ = 0 and α = 0 for the left plot and at ωτ = 0 and = 0 for the right plot.

Figure 7 .
Figure 7. Temperature dependence of the flipping time of magnetization τ measured for (a) 1, (b) 2, and (c) 3 under an application of 1000 Oe DC field.(left) Arrhenius plots and (right) ln(τ) vs. ln(T) plots.The blue lines represent theoretical calculations based on the linear Arrhenius equation for the Orbach process (left) or based on Equation (3) for the Raman process (left).The blue curves represent theoretical calculations considering both Orbach and Raman processes, expressed as Equation (4).

Figure 7 .
Figure 7. Temperature dependence of the flipping time of magnetization τ measured for (a) 1, (b) 2, and (c) 3 under an application of 1000 Oe DC field.(left) Arrhenius plots and (right) ln(τ) vs. ln(T) plots.The blue lines represent theoretical calculations based on the linear Arrhenius equation for the Orbach process (left) or based on Equation (3) for the Raman process (right).The blue curves represent theoretical calculations considering both Orbach and Raman processes, expressed as Equation (4).

Table 1 .
Best fit kinetic parameters for magnetization flipping estimated for complexes 1-3 on the basis of the Equation (4).

Table 2 .
Calculated atomic spin densities assigned to each f orbitals.

Table 2 .
Calculated atomic spin densities assigned to each f orbitals.