Magnetic Interactions in a Series of Homodinuclear Lanthanide

Nov 20, 2015 - CIF files for [Y2L(OAc)4]PF6, [Gd2L(OAc)4]PF6, [Tb2L(OAc)4]PF6, [Dy2L(OAc)4]PF6, [Ho2L(OAc)4]PF6, [Er2L(OAc)4]PF6 and [Lu2L(OAc)4]PF6 a...
0 downloads 5 Views 3MB Size
Article pubs.acs.org/IC

Magnetic Interactions in a Series of Homodinuclear Lanthanide Complexes Peter Comba,*,† Michael Großhauser,† Rüdiger Klingeler,‡ Changhyun Koo,‡ Yanhua Lan,§,∥ Dennis Müller,† Jaena Park,‡ Annie Powell,§,∥ Mark J. Riley,⊥ and Hubert Wadepohl† †

Anorganisch-Chemisches Institut und Interdisziplinäres Zentrum für Wissenschaftliches Rechnen (IWR), Universität Heidelberg, INF 270, D-69120 Heidelberg, Germany ‡ Kirchhoff-Institut für Physik, Universität Heidelberg, INF 227, D-69120 Heidelberg, Germany § Institut für Anorganische Chemie, Karlsruhe Institute of Technology (KIT), Engesserstrasse 15, Geb. 30.45, D-76131 Karlsruhe, Germany ∥ Institute of Nanotechnology (INT), P.O. Box 3640, D-76021 Karlsruhe, Germany ⊥ School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, Queensland 4072, Australia S Supporting Information *

ABSTRACT: A series of seven isostructural homodinuclear lanthanide complexes are reported. The magnetic properties (ac and dc SQUID measurements) are discussed on the basis of the Xray structural properties which show that the two lanthanide sites are structurally different. MCD spectroscopy of the dysprosium(III) and neodymium(III) complexes ([DyIII2(L)(OAc)4]+ and [NdIII2(L)(OAc)4]+) allowed us to thoroughly analyze the ligand field, and high-frequency EPR spectroscopy of the gadolinium(III) species ([GdIII2(L)(OAc)4]+) showed the importance of dipolar coupling in these systems. An extensive quantum-chemical analysis of the dysprosium(III) complex ([DyIII2(L)(OAc)4]+), involving an ab initio (CASSCF) wave function, explicit spin−orbit coupling (RASSI-SO), and a ligand field analysis (Lines model and Stevens operators), is in full agreement with all experimental data (SQUID, HF-EPR, MCD) and specifically allowed us to accurately simulate the experimental χT versus T data, which therefore allowed us to establish a qualitative model for all relaxation pathways.



INTRODUCTION

and eventually lead to a basis which may enable us to optimize lanthanide-based molecular magnets.18−22 We report a series of eight isostructural homodinuclear lanthanide(III) complexes with the general formula [Ln2L(OAc)4]PF6, where LnIII = Y, Nd, Gd, Tb, Dy, Ho, Er, and Lu and L = 2,6-bis((bis(pyridin-2-ylmethyl)amino)methyl)-4methylphenolate (see Chart 1). The aim was to thoroughly investigate the magnetic exchange between the lanthanide centers, in particular to fully understand the electronic structure

Single molecular magnetism (SMM) was first described some 25 years ago and has since attracted much attention,1−4 not only due to a number of possible applications5,6 but also because the fundamental reasons for bistability and relaxation, specifically the anisotropy barrier and quantum tunneling, are known to be related to the electronic structure of the metal sites, magnetic exchange, and magnetic anisotropy.7−9 These are related to molecular structure, and therefore in theory are predictable.10 While the investigation of SMMs in the 1990s focused on oligonuclear 3d metal complexes,11 the larger spin− orbit coupling of second and third row transition metal ions and the lanthanides, specifically the large intrinsic anisotropy of the lanthanides, more recently started to attract attention and resulted in an increasing number of oligonuclear lanthanide and mixed-metal complexes with SMM behavior.12−16 Major problems in this area still are the predictability of structures and their correlation with electronic properties,12,17 but interesting developments in theory may help to more deeply understand experimental data on coupled lanthanide systems © XXXX American Chemical Society

Chart 1. Structure of the Ligand L

Received: July 30, 2015

A

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry and to quantify the anisotropy. We have used experimental structural, magnetic, optical spectroscopic, and high-field electron paramagnetic resonance data and have also analyzed the system computationally with ab initio CASSCF calculations and introduced spin−orbit coupling via the restricted active space state interaction method (RASSI).23 The magnetic coupling was then analyzed on the basis of the “Lines model”,24 combined with a thorough ligand field analysis, based on extended Stevens operators,25 using the quantumchemical data. This allowed us to fully interpret the experimental data and, importantly, to carefully check for consistency between the various experimental and computational data sets. Therefore, the combined experimental and computational analysis serves as a benchmark for the theoretical model, which may now be used for the design and optimization of lanthanide-based SMMs.



RESULTS AND DISCUSSION Syntheses and Structural Properties. The preparation of the homodinuclear lanthanide complexes involved lanthanide(III) acetate salts and the dinucleating ligand L (Chart 1). The complexes of YIII, NdIII, GdIII, TbIII, DyIII, HoIII, ErIII, and LuIII were isolated and characterized;, from all but NdIII X-ray crystal structures were obtained. All structures are isomorphous and contain approximately three molecules of methanol per complex in the crystal lattice. Figure 1 shows the molecular structure of the bisdisprosium(III) complex, and similar structural plots of all other complexes are given as Supporting Information; selected structural parameters of all structures are listed in Table 1. The cations are asymmetric with different coordination numbers at the two lanthanide centers (Ln1 is 8-coordinate with an N3O5 chromophore, Ln2 is 9-coordinate with an N3O6 donor set). L provides two pyridine and a tertiary amine donor to each site; the two metal cations are bridged by the μ2-OPh phenolate of the dinucleating ligand and by two acetate anions. The latter have different bonding modes: one is terminal at the 9coordinate center 2 and bridges to the other center via a μ1-O5Ac oxygen atom. The other acetate at site 2 is a μ2-CH3COOκ2O:O′ moiety and coordinates to the two centers with a μ1-O3Ac and an O4Ac atom. Both metal centers are also coordinated to terminal acetates. The average LnIII−O distances vary significantly, and the distance between the two metal ions is between 3.767 Å (ErIII) and 3.812 Å (TbIII, see Table 1). The inner coordination sphere of the 9-coordinate site 2 is best described as a monocapped distorted antiprism, and the 8coordinate site 1 is a triangular dodecahedron; both centers have highly distorted icosahedral symmetry (see Figure 1). Intermolecular magnetic exchange often plays an important role in magnetic and spectroscopic properties of solids.26 Therefore, as an example the structure of the dinuclear DyIII2 complex was evaluated for intermolecular interactions. It appears that there is only one significant intermolecular contact, i.e., a π−π interaction (see Figure 2): the bridging phenolate ring (Ph) aligns with a pyridine donor of the 8coordinated site 1 (py1) as well as a pyridine moiety of the 9coordinated site 2 (py2). The distances and angles between the aromatic rings are also summarized in Table 1. Electronic Properties. i. Magnetism. Direct current (dc) magnetic susceptibilities were measured on powdered crystals of the samples used for crystallography at 500 G over a temperature range 2−300 K. χMT versus T plots of all [Ln2(L) (OAc)4]PF6 complexes are shown in Figure 3.

Figure 1. Chemical structure, ORTEP plot of the full structure (hydrogen atoms and anions omitted) and of the inner coordination sphere of [DyIII2(L)(OAc)4]+, and coordination polyhedra of the two lanthanide sites (left, blue, 8-coordinate Dy1 site; right, red, 9coordinated Dy2 site).

The χMT value of the dinuclear NdIII2 species at 300 K is 2.94 cm3 K mol−1; for two independent NdIII ions with a 4I9/2 ground state a somewhat higher χMT value of 3.28 cm3 K mol−1 is expected at ambient temperature. In a known alkoxy-bridged dinuclear NdIII complex, a similar value for χMT = 2.84 cm3 K mol−1 was observed at 300 K.27 The χMT value decreases steadily with decreasing temperature to a value of 0.94 cm3 K mol−1 at 2.3 K, in agreement with antiferromagnetic intramolecular coupling. The χMT value of the GdIII2 complex is 15.10 cm3 K mol−1 at 300 K, somewhat smaller than the expected value 15.82 cm3 K mol−1 of two uncoupled 8S7/2 centers. χMT slowly decreases to 14.16 cm3 K mol−1 at 23.3 K and drops to 9.71 cm3 K mol−1 at 2.0 K. The almost temperature independent shape for a dinuclear GdIII species28 indicates a weak antiferromagnetic B

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry Table 1. Selected Inter- and Intramolecular Distances (Å) and Angles (deg) of [LnIII2(L)(OAc)4]PF6·nMeOH [LnIII2(L)(OAc)4]+ Ln1−Ln2 Ln1−O1AC1 Ln1−O2AC Ln1−O3AC Ln1−O5AC Ln2−O4AC Ln2−O5AC Ln2−O6AC Ln2−O7AC Ln2−O8AC Ln1−O9Ph Ln2−O9Ph Ln1−N1 Ln2−N2 Ln1−N1Py Ln1−N2Py Ln2−N3Py Ln2−N4Py ∠Ln1−O9Ph−Ln2 Py2′−Py1 Py1−Ph Ph−Py2 Py2−Py1′ ∠Py2′−Py1 ∠Py1−Ph ∠Ph−Py2 ∠Py2−Py1′

Y

Gd

Tb

Dy

Ho

Er

Lu

3.775(3) 2.363(4) 2.373(4) 2.284(4) 2.276(4) 2.347(4) 2.450(4) 2.456(4) 2.402(4) 2.408(4) 2.312(4) 2.342(4) 2.549(4) 2.585(5) 2.536(5) 2.535(5) 2.552(5) 2.554(4) 108.44(15) 4.073 3.843 3.471 4.073 21.86 37.87 20.64 21.86

3.791(1) 2.373(2) 2.407(2) 2.313(2) 2.317(1) 2.375(2) 2.464(2) 2.508(2) 2.438(2) 2.432(2) 2.334(1) 2.360(1) 2.570(2) 2.633(2) 2.575(2) 2.570(2) 2.557(2) 2.558(2) 107.73(6) 4.040 3.857 3.453 4.040 20.45 38.40 20.93 20.45

3.812(2) 2.381(3) 2.400(3) 2.305(3) 2.307(3) 2.373(3) 2.473(3) 2.493(3) 2.425(3) 2.420(3) 2.340(3) 2.358(3) 2.572(3) 2.603(4) 2.557(4) 2.548(4) 2.562(3) 2.578(3) 108.50(11) 4.085 3.867 3.477 4.085 20.80 38.36 20.46 20.80

3.794(1) 2.372(2) 2.385(2) 2.299(2) 2.304(2) 2.354(2) 2.459(2) 2.477(2) 2.415(2) 2.410(2) 2.313(2) 2.356(2) 2.563(2) 2.594(2) 2.548(2) 2.536(2) 2.552(2) 2.562(2) 108.73(6) 4.084 3.858 3.473 4.084 21.68 38.24 20.77 21.68

3.773(2) 2.360(3) 2.372(3) 2.283(3) 2.285(3) 2.345(3) 2.445(3) 2.459(3) 2.399(3) 2.401(3) 2.314(3) 2.339(3) 2.549(3) 2.579(4) 2.535(4) 2.517(4) 2.538(3) 2.557(3) 108.35(11) 4.081 3.837 3.464 4.081 22.03 38.28 20.53 22.03

3.767(1) 2.353(2) 2.372(2) 2.271(2) 2.277(1) 2.340(2) 2.439(2) 2.446(2) 2.392(2) 2.402(2) 2.308(2) 2.330(1) 2.540(2) 2.580(2) 2.523(2) 2.518(2) 2.532(2) 2.550(2) 108.61(6) 4.135 3.847 3.462 4.135 23.04 38.62 20.24 23.04

3.729(1) 2.324(2) 2.353(2) 2.241(2) 2.251(2) 2.310(2) 2.419(2) 2.415(3) 2.372(2) 2.385(2) 2.273(2) 2.300(2) 2.513(2) 2.553(3) 2.499(3) 2.498(3) 2.508(3) 2.525(2) 109.25(8) 4.090 3.813 3.458 4.090 23.37 37.86 20.29 23.37

Figure 2. Mercury plot for the [DyIII2(L)(OAc)4]+ and its neighboring molecules. Black spheres show ring centroid for the interacting pyridines and phenols. The dashed lines show the direct contact to the next aromatic ring.

where g is the Landé factor, β the Bohr magneton, N the Avogadro number, and k the Boltzmann constant; a good fit of the experimental data is found with J = −0.065 cm−1 and g = 1.97. The small coupling strength J obtained from the fit is indeed in agreement with the hypothesis from the thermal evolution of the χMT curve. Other comparable dinuclear GdIII complexes have only slightly larger antiferromagnetic coupling constants.30 The χMT value of the TbIII2 compound at 300 K with 24.10 cm3 K mol−1 is higher than the expected value of 22.78 cm3 K mol−1 for a 7F6 ground state. With a decrease in the temperature to 60 K, there is a small decrease to 22.86 cm3 K mol−1, with 16.42 cm3 K mol−1 at 2 K.

interaction between the GdIII centers. With the known isotropic ̂ ŜGd2 results in the interaction, the Hamiltonian Ĥ = −JSGd1 formula for the temperature dependence of the molar magnetic susceptibility in eq 129 χM =

2g 2Nβ 2 kT × [140e56J / kT + 91e 42J / kT + 55e30J / kT + 30e 20J / kT + 14e12J / kT + 5e6J / kT + e 2J / kT ]/[15e56J / kT + 13e 42J / kT + 11e30J / kT + 9e 20J / kT + 7e12J / kT + 5e6J / kT + 3e 2J / kT + 1]

(1) C

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 3. χMT vs T plot of all [LnIII2(L)(OAc)4]PF6 complexes together with a fit of the GdIII compound (eq 1, g = 1.97, J = −0.065 cm−1).

Figure 4. Resonance frequencies vs magnetic field of the dinuclear GdIII complex, at 4 K (●). Representative EPR spectra are shown as gray lines. The black solid line displays the simulated frequency dependence of the ground state resonance for magnetic fields parallel to the easy axis. Parameters are given in the text. Note that the linear fit to the data mentioned in the text nearly perfectly coincides with the simulated data and hence is not shown.

The χMT value of 27.34 cm3 K mol−1 for the DyIII2 complex at 300 K is lower than the expected value of 28.36 cm3 K mol−1 for a 6H15/2 ground state. This decreases slightly down to a minimum of 23.52 cm3 K mol−1 at 2.3 K and then slightly rises to 23.73 cm3 K mol−1 at 2 K. The χMT value for the HoIII2 species at 300 K is 27.54 cm3 K mol−1, i.e., lower than the theoretical value of 28.10 cm3 K mol−1 for a 5I8 ground state. With a continuous decrease it reaches 7.66 cm3 K mol−1 at 2 K. The χMT value of the dinuclear ErIII2 complex is 21.36 cm3 K mol−1 at 300 K, significantly lower than the expected 22.95 cm3 K mol−1 for two uncoupled centers with a 4I15/2 ground state. After a continuous decrease, χMT reaches 11.11 cm3 K mol−1 at 2 K. As described above, the general feature of the χMT curves for the four heavy lanthanide complexes is very similar. The continuous decrease of χMT with the decrease of temperature could partially result from thermal depopulation of the Stark levels, i.e., from the splitting of the free-ion ground state of LnIII by the crystal field. On the basis of the analysis of the data of the GdIII compound, the exchange coupling between the lanthanides in these four isostructural dinuclear compounds is probably also antiferromagnetic in nature, but should be very weak to the point of being virtually nonexistent. For the DyIII2 complex the temperature dependent ac susceptibility data were also collected (see Supporting Information). The isolines are split in the out-of-phase susceptibility. No maximum is observed, implying that this may occur below 1.8 K, the lowest possible temperature of the setup used. This indicates only very weak SMM behavior. Very rough estimates of Ueff and τ0 may be obtained from eq 231,17 ln(χ ″ /χ ′) = ln(ωτ0) + Ueff /kT

axis.32,33 We note that this behavior does not appear in randomly oriented powder spectra, but it is typical for the particular case of the aligned powder, i.e., where all local anisotropy axes are aligned along the applied magnetic field. As will be described below, our analysis suggests a dipolar origin of the anisotropy. The temperature dependence of the spectra at 130 GHz supports the easy axis scenario as spectral weight is shifted from lower to higher magnetic field upon heating (see Supporting Information). The anisotropy is, however, small so that the excited state transitions and the low-lying state transitions are not clearly separated. We attribute the broad shoulder seen in the resonance features to only partially separated resonances. The resonance fields of the ground state transition at 4 K are shown in the resonance frequency versus magnetic field diagram in Figure 4. Linear fitting of the resonance branches provides the g-factor and the zero-field splitting (ZFS), i.e., g = 2.01 ± 0.01 and Δ = 5.9 ± 1.0 GHz. We mention the nonvanishing size of the ZFS which in the case of GdIII ions with a half-filled 4f-shell and negligible spin−orbit coupling does not originate from single-ion anisotropy. We conclude that magnetic dipolar coupling accounts for the observed ZFS. Quantitatively, the magnetic anisotropy tensor induced by magnetic dipolar coupling for weakly coupled pairs of ions is written as in eq 3.34−36 In eq 3, βr = μB2/r3; μB is the Bohr

(2) −6

leading to Ueff ∼ 1.97 K and τ0 = 1.62 × 10 s. This is comparable to barriers of similar compounds.31 ii. High-Frequency EPR Spectroscopy. Representative HFEPR spectra of the GdIII2 complex obtained at 4 K display broad resonance features (see Figure 4). The shape of a particular resonance feature allows qualitative conclusions about the presence and the sign of anisotropy as maximum absorption appears in the lower field region of the resonances while the resonance intensity smoothly decreases as the magnetic field increases. Such a behavior in general suggests an easy-axis-type of anisotropy. To be specific, maximum absorption at low field and at low temperature indicates the ground state transition for the magnetic field parallel to the easy

magneton, and r is the distance between two ions (r = 3.791 Å for the GdIII dinuclear complex at hand). ξ(η) is the angle between the ionic z (y)-axis and the distance vector r̂. All the terms here are referred from ref 35. The spin Hamiltonian including the anisotropy term induced by magnetic dipolar coupling can be written as in eq 4 2

⃡ Ĥ = −JS1̂ ·S2̂ + gμB ∑ Sî ·B̂ + S1̂ ·Ddipolar ·S2̂ i=1

(4)

with the intradimer exchange integral J extracted from the magnetic susceptibility. The Hamiltonian of eq 4 is applied to D

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 5. (a) Overview of the [NdIII2(L)(OAc)4]+ MCD spectrum at 8.0 K and 5.0 T. (b) As above, scaled to show the weaker transitions. The stick spectra are the calculated values for the free ion (black), the Nd1 center (red), and the Nd2 center (blue). The transitions to the 2P1/2 and 2D1/2 states are seen for both centers as indicated. The observed and calculated peak positions are given as Supporting Information.

the discussion of the spectra is guided by ligand field calculations. Ligand Field Calculations. The calculations of the ligandfield-split states of the dinuclear complexes are made assuming that the interaction between the centers is weak, i.e., that the optical spectra can be interpreted in terms of a superposition of single ion spectra. The Hamiltonian to be solved is the usual series of atomic terms together with the ligand field of the atomic environment. The atomic terms include the electron− electron repulsion (F2, F4, F6), spin−orbit coupling (ζ), and the higher order terms that include the two body configuration interaction parameters (α, β, γ), the three body parameters (Tk, k = 2, 3, 4, 6, 7, 8), the magnetic parameters (Mk, k = 0, 2, 4), describing the spin−spin and spin−other orbit interactions, and the electrostatically correlated spin−orbit interaction (Pk, k = 2, 4, 6). These 19 atomic parameters for NdIII37 and DyIII38were taken from the literature, where they were obtained from the fitting to LnCl3 spectra. These data sets are considered reliable and not plagued by errors due to the use of incorrect reduced matrix elements as in other calculations. Details of the ligand field calculations are tabulated in the Supporting Information. For a lanthanide ion in a ligand field with no formal symmetry, there are 27 nonzero Bkq parameters,39 far more than can reasonably be determined with the available experimental data. Here, we make use of the relationship between the angular overlap model (AOM) and the Bkq parameters.39 In this manner, the ligand field is determined by the geometry of the first coordination sphere and adjustable eσ, eπ parameters. Details of the ligand field parametrization are given in the Supporting Information. The calculation of the f state energy levels is well-known.39−41 Good agreement was found with

simulate the ZFS for the magnetic field parallel to the easy axis. η does not affect the results significantly for the magnetic field parallel to the easy axis. The best simulation parameters with ξ = 30°, η = 0°, and isotropic g = 2.01 for a Gd3+ ion yields Δ = 6.00 GHz, which matches the experimental results of Δ = 5.9 ± 1.0 GHz. Our analysis hence confirms the dipolar origin of ZFS observed in the HF-EPR data of the GdIII dinuclear complex. We note that when the magnetic field determines the easy axis direction, the observed easy-axis-like behavior for the aligned case is straightforwardly explained. iii. MCD Spectroscopy. Figures 5a and 6a show an overview of the MCD spectra of the [Nd2(L)(OAc)4]+ and [Dy2(L)(OAc)4]+ complexes, respectively, together with the assignment of the bands to transitions between appropriate multiplets; detailed spectra of particular transitions are shown in Figures 7 and 8, respectively. The spectra were obtained from the complexes dissolved in a 3:1 ethanol/methanol glass at 5 T and 8 K. The absorption electronic spectra were measured simultaneously, but were generally of low quality due to the low extinction coefficients and a background that was difficult to subtract. The phase-sensitive detection employed in the MCD measurement is more sensitive to small signals, has minimal background, and allows the f−f transitions to be clearly observed. The spectra can be interpreted in terms of a superposition of the single ion spectra of f−f transitions. For Nd (Dy) there are a total of 182 (1001) Kramers doublets, of which there are 96 (116) expected to occur in the 9000−30 000 cm−1 spectral range. As each lanthanide site in the [Ln2(L)(OAc)4]+ complexes experiences a different ligand field, their superimposed spectra give rise to a large number of transitions, and E

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 6. (a) Overview of the [DyIII2(L)(OAc)4]+ MCD spectrum at 8.0 K and 5.0 T. (b) As above, scaled to show the weaker transitions. The stick spectra are the calculated values for the free ion (black), the Nd1 center (red), and the Nd2 center (blue). The observed and calculated peak positions are given as Supporting Information.

calculations, where the spin−spin interaction part of Mk parameters was included.37 As we are interested in the interaction of the energy levels with a magnetic field we have also calculated the g-values, assuming isolated Kramers doublets. We have used an established methodology42 to calculate and diagonalize the g2-tensor, which gives the principle values and directions (see Supporting Information). Discussion of the Spectra. Figures 5b and 6b show a more detailed view of the spectra, and these also include stick spectra of the calculated energy positions of the free ion (black) as well as for each lanthanide center Ln1 (red) and Ln2 (blue). The height of the stick spectra is the fraction of the calculated wave function that has the same spin as the ground state, which gives some indication of the relative intensity of the transitions. The experimental and calculated energies for each lanthanide center are given as Supporting Information. It is important to note that no fitting of the parameters to the experimental spectra has been attempted, as the spectra of the superimposed centers need 73 (19 + 2 × 27) parameters. However, some insights can be gained from using the same atomic parameters for each center, with different ligand fields specific to Ln1 and Ln2 being used. Therefore, for [NdIII2(L)(OAc)4]+ and [DyIII2(L)(OAc)4]+ the same ligand fields are used for the Ln1 and Ln2 coordination geometry. Comparing the calculated stick spectra for Nd1 and Nd2 in Figure 5b, one can see the shift of the free-ion levels (black) to higher energy (red/blue) when the ligand field is added to give good agreement with the experimental levels. The spread of the calculated levels and that of the experiment also agree quite well. In general, the experimental intensities correlate with the calculated fractional spin composition.

Figure 7 displays the temperature dependence of the [NdIII2(L)(OAc)4]+ spectra in the region of the 4I9/2 → 4G5/2 and 4I9/2 → 4D5/2, 4D3/2,4D1/2 transitions, clearly showing the C-term behavior expected for this odd electron system. With a fit of the integrated intensity to a simple S = 1/2 function in a VTVH plot, one finds effective g-values of 4.850 and 3.197 from the temperature dependence of the 4I9/2 → 4G5/2 and 4I9/2 → 4 D5/2,4D3/2, 4D1/2 groups of transitions, respectively. In a comparison of the calculated ground state g-values of Nd1 (0.405, 0.977, 4.958) and Nd2 (0.006, 1.618, 3.561, see Supporting Information), this suggests that the former transition is xy polarized while the latter may have xy and z polarized components. The negative MCD signals of the 4I9/2 → 4D1/2 transitions in Figure 7 are interesting, as there is only a single predicted transition. The ∼190 cm−1 separation of these peaks must therefore be due to the 4I9/2 → 4D1/2 transitions of each metal site. The 4D1/2 state in Nd1 is calculated to be 60 cm−1 higher than that in Nd2. A calculation of the energy levels as a function of the ligand field strength is shown in the Supporting Information, and the 4D1/2 state remains separate from the nearby bands. Another calculated isolated single state is 2P1/2 at ∼23 300 cm−1 (Supporting Information), and indeed, two peaks are observed in this region (see Figure 5b and Supporting Information). Again, the higher energy peak is predicted to be due to Nd1, but the calculated separation of 15 cm−1 is significantly smaller than that observed (123 cm−1), indicating that the difference in the ligand fields between the two sites may be larger than assumed here. The [Dy2(L)(OAc)4]+ spectrum shown in Figure 6a is dominated by the 6H15/2 → 6Fn (n = 9/2, 7/2, 5/2, 3/2) transitions in the near-infrared and the 6H15/2 → 6P5/2, 6P7/2 F

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 7. (a) MCD spectrum of [NdIII2(L)(OAc)4]+, showing the transitions of the 4I9/2 → 4G5/2 multiplet as a function of temperature. (b) MCD spectrum of [NdIII2(L)(OAc)4]+, showing the transitions of the 4I9/2 → 4D5/2 + 4D3/2 multiplets as a function of temperature.

transitions in the blue region. At energies above 20 000 cm−1, the quartet and sextet states are very mixed (see Supporting Information) although it is easy to identify the 3 components of the 6P5/2 and the 4 components of the 6P7/2 states embedded in the quartet states. The transitions to the components of these two states appear as a single positive band and a single negative band, respectively. The 6H15/2 → 6F1/2 transition is calculated to be an isolated single peak, but as shown in Figure 8a, it does not appear experimentally. This is due to the ΔJ > 6, involved for this transition which results in it being forbidden.43 However, the small splitting observed in the 6H15/2 → 6F3/2 transition is likely due to the two lanthanide sites as the calculated separation of the two 6F3/2 components within each Dy site is 8 (Dy1) and 2 (Dy2) cm−1, while the separation of the mean energies between the two sites is 40 cm−1, which is reasonably close to the observed 35 cm−1. The observed 23 cm−1 splitting in the 6H15/2

→ 6F5/2 transition is more ambiguous, as the splitting between the 6F5/2 components is comparable with energy differences between the two sites. Figure 8b shows the temperature dependence of the 6H15/2 → 6F7/2 transition. This is less than that observed in the [NdIII2(L)(OAc)4]+ complex, implying a larger g-value. At the lower edge of the band, derivative shaped A-terms develop, implying a difference in the ground and excited state g-values, so that transitions of the opposite sign form an A-term rather than cancel as seen in the [NdIII2(L)(OAc)4]+ spectra of Figure 7. The additional A-term that appears at higher temperatures 81 cm−1 to lower energy is then a hot band due to the population of higher levels of the split ground state. The calculated first energy levels of 39 (Dy1) and 44 (Dy2) cm−1 again may indicate that the ligand field is larger than estimated. Computation of the Electronic Properties. The calculation of the wave function was done with ab initio G

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 8. (a) MCD spectrum of [DyIII2(L)(OAc)4]+, showing the transitions of the 6H15/2 → 6F5/2, 6F3/2 multiplets. Note the complete absence of the 6H15/2 → 6F1/2 transition. The calculated positions of the Dy1 (red) and Dy2 (blue) centers are shown as stick spectra. (b) MCD spectrum of [DyIII2(L)(OAc)4]+, showing the transitions of the 6H15/2 → 6F7/2 multiplet as a function of temperature.

only have a small contribution. After obtaining the spin-free wave functions through CASSCF,47 spin−orbit coupling was introduced via the restricted active space state interaction method (RASSI).23 Due to the high number of states arising from the given active space, only a limited number can be mixed by spin−orbit coupling. All sextet, 50 doublet, and 50 quartet states were considered. The resulting RASSI-SOcoupled multiplet states were introduced to the Molcas SINGLE-ANISO module44 to calculate the local electronic and magnetic properties of the single magnetic center fragments. To further account for the magnetic coupling between the two lanthanide centers the “Lines model” was applied.24 Even though this is a solid state isotropic exchange model, one can apply the effective Hamiltonian on the anisotropic basis (spin− orbit coupled multiplets from RASSI), calculated for the single magnetic center fragments beforehand. This procedure is welldescribed and implemented in the POLY_ANISO mod-

calculations on the crystal structural coordinates with Molcas 8.0.44 An established fragmentalized approach was used, where the neighboring magnetic centers are substituted with their closest diamagnetic analogue; i.e., DyIII was substituted by LuIII with an ionic radius closer than that of LaIII.19 This results in two fragments with only one magnetic center to calculate. All atoms were described by contracted ANO-RCC basis sets, included in Molcas.45 The active space of the employed CASSCF method included nine electrons in seven 4f orbitals. Additional dynamic correlation through second-order perturbation theory was not considered, since the effect is known to be small for lanthanides.46 From the ligand field manifold for a 4f9 electronic configuration arise 21 sextet, 224 quartet, and 490 doublet states. Sufficient results can be obtained with all sextet, 50 doublet, and 50 quartet states. Considering more roots increases calculation time and demand drastically and has very little influence on the energy of the ground and low-lying excited states, since they correspond to high-energy states that H

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry ule.46,48−50 One core feature of the Lines model to bear in mind is that it becomes exact in the case of two Ising-type ions (strong axial magnetization). With the resulting wave function for the dinuclear complex the magnetic susceptibility from dcSQUID data can be simulated (see Figure 10). The only parameters are the exchange coupling constant J and a parameter zJ, which accounts for intermolecular interactions in the powder by a mean-field approach.

wavenumber. The fast relaxation exhibited by the dinuclear DyIII complex, as seen in the experimental magnetic data indicates a very low blocking barrier, which is well-reproduced by the computational results, and this also emerges from the computed exchange spectrum given in Figure 11. Apart from the noncolinearity of the local magnetic moments on the DyIII centers, the nonaxiality of the ligand field proves to be an equally big problem. This becomes evident when looking at the ligand field in terms of the extended Stevens operators (see Supporting Information).25 Since the system has C1symmetry no main quantization axis by standard convention can be given. The standard interpretation of the ESO parameters remains intact, even in nonaxial systems, but a total of 27 parameters are needed to describe the exact energy splitting. The z-axis was chosen to align with the orientation of the magnetic moment of the ground state Kramers doublet (Zm, see Figure 9). The axiality of the ground state KD (see Table 2 KD1(gx, gy, gz)) is due to the perturbation by a close ligand atom. The Bn0-values for the Dy2 center are all negative. This is the result of the axial component of the ligand field being stronger than the equatorial one; a strong axiality is needed to stabilize the MJ = ±15/2 state.12 Especially, the Dy1 center shows an odd term B21 which is non-negligible and points to a rather weak axiality of the system. The two DyIII ions experience a vastly different ligand field, which is evident in the disjoined local anisotropy axes. This strengthens the belief that lanthanides are not that inert to ligand field effects as usually assumed and emphasizes the importance of a strongly axial ligand field.51 Without a strong axiality of the ground and several excited Kramers doublets, one will always observe a very fast relaxation of one magnetic center, because of the possibility of reversal of magnetization by quantum tunneling (QTM) and Orbach or Raman relaxation pathways. Tunneling comes more into effect as the gx and gy components increase.52 There are indications that the deviation of the orientation of the magnetic moments of the ground state Kramers doublet to the excited ones plays a major part in supporting spin−phonon relaxation pathways.53 As of now there is no method to quantify the three possible relaxation pathways (thermal, QTM, Orbach), but a qualitative insight can be gained by the transition magnetic moment connecting the respective states (Figure 11).54 This leads to the requirements that have to be met by a Dy-based SMM. A strong axial ligand field that leads to an energy separation of the ground state and excited state Kramers doublet by at least 100 cm−1, a gx, gy < 0.2 and gz > 18 and a colinearity of the orientation of magnetic axes between the ground and excited state Kramer doublets.

Figure 9. Orientation of the ab initio calculated main magnetic axes of the ground state Kramer doublets of [DyIII2(L)(OAc)4]+ on the magnetic centers Dy1 (left) and Dy2 (right).



CONCLUSION A large range of experimental and computational methods have been applied to a complete series of isostructural dinuclear lanthanide complexes. The ligand field analysis and ab initio calculations of the DyIII2 complex are in excellent agreement with the experimental results, and the quantum-chemical calculations provide a good understanding of the fundamental reasons for the fast relaxation. The weak antiferromagnetic coupling of the DyIII ions and the ill-matched local anisotropy, due to variations in local ligand fields, are the main reasons for the relaxation pathways. The results show a strong anisotropy but a very small barrier, leading to fast relaxation of magnetization. This highlights the importance of the coordination geometry, when designing small lanthanidebased SMMs. Although spin−orbit effects are of a much larger

Figure 10. Black: Measured χT vs T data of [DyIII2(L)(OAc)4]+. Red: Simulation of the χT vs T data on the basis of the ab initio calculated wave function with the exchange coupling constant J = −1.2 cm−1 and the parameter zJ = −0.02 cm−1, to account for intermolecular coupling in the measured powder.

The directions of the local anisotropy axes for the dinuclear DyIII complex are shown in Figure 9; the computed parameters are given in Table 2. This indicates that the local axes do not align, resulting in an easy plane and not an easy axis for the overall anisotropy of the complex. The simulation of the experimental χT versus T data yields an energy spectrum of the four lowest exchange states (Figure 11). The overall height of the potential is very low. Especially, the energy gap between the first excited exchange state and the ground state is less than a I

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Table 2. Energies and g-Tensor Components of the RASSI-SO Coupled Eight Kramers Doublets of the DyIII Ground State Multipleta E [cm−1] Kramers doublets

a

Dy1

∠KD1−Dy−KDX [deg] Dy2

1

0.000

0.000

2

82.678

47.380

3

126.352

101.130

4

168.99

135.012

5

190.307

190.559

6

244.019

227.811

7

292.388

250.966

8

505.778

308.520

gx gy gz gx gy gz gx gy gz gx gy gz gx gy gz gx gy gz gx gy gz gx gy gz

Dy1

Dy2

0.0513 0.0807 19.6774 0.0105 0.2256 17.9927 1.0239 2.1238 14.5144 2.3042 3.6608 12.5689 2.4158 5.6653 9.1808 1.9934 2.8726 13.4074 0.4709 0.5829 18.5805 0.0087 0.0115 19.7885

0.1523 0.1785 19.5093 0.1085 0.2098 19.3994 2.2305 3.2754 13.0475 9.5218 8.1274 2.3023 1.6465 3.9979 10.9756 1.8470 5.7434 10.9760 0.9296 1.9883 15.3747 0.2281 0.6409 18.2012

Dy1

Dy2

81.3

47.8

133.7

17.5

85.9

23.5

84.1

68.3

96.3

80.5

101.3

98.0

58.8

83.7

The angle KD1−Dy−KDX gives the deviation of the main magnetic axis of each excited Kramers doublet with respect to the ground state KD1.

complexes, provided the structures of the complexes are known.



EXPERIMENTAL SECTION

Measurements. Magnetic Measurements. Magnetic measurements were carried out on a MPMS-XL 5T (Quantum Design) SQUID magnetometer. Samples were powdered and pressed in PTFE tape to avoid field-induced orientation. Data were corrected for diamagnetism of the sample holder, and Pascal’s constants were used for diamagnetic corrections of the sample. Dynamic magnetic ac SQUID measurements were performed with a Quantum Design MPMS-XL7 SQUID. Polycrystalline samples were embedded in vaseline to prevent torqueing. High-Frequency Electron Paramagnetic Resonance (HF-EPR). These measurements were performed on a powder sample of the dinuclear GdIII complex in the frequency range ν = 32−330 GHz, and in a temperature range from 4 to 150 K. A millimeter vector network analyzer (MVNA) was used as a source and detector of stable microwave, and a superconducting magnet was applied to provide a magnetic field up to 15 T.55 Magnetic Circular Dichroism (MCD). Spectra were obtained in the range 350−1400 nm using a Jobin Yvon 750s monochromator with a Hinds photoelastic modulator to generate the circularly polarized modulation with either a Si-avalanche photodiode or an InGaAs detector for the visible and near-IR regions, respectively. The signals were acquired with Stanford SR830 lock-in amplifiers at chopped and circularly polarized modulations allowing the simultaneous measurement of absorption and MCD spectra. The sample was held in a 5 Telsa magnetic field of an Oxford SM-4 Spectramag at temperatures 2−10 K using an Oxford ITC504 controller. The MCD spectrum was calibrated using Co(en)33+ as a standard for circular dichroism.56 X-ray Crystal Structure Determinations. Crystal data and details of structure determinations are given in the Supporting Information. CIF files for [Y2L(OAc)4]PF6, [Gd2L(OAc)4]PF6, [Tb2L(OAc)4]PF6, [Dy2L(OAc)4]PF6, [Ho2L(OAc)4]PF6, [Er2L(OAc)4]PF6 and [Lu2L-

Figure 11. Ab initio calculated exchange spectrum of [DyIII2(L)(OAc)4]+ with four basis functions per magnetic center; the numbers at each arrow indicate the mean average matrix elements [(|μx| + |μy| + |μz|)/3] connecting the stationary points to give a qualitative idea of the most probable relaxation pathway. Tunelling beween respective states is expected to be relevant above a factor of approximately 10−3.

dimension, the energy barriers for relaxation are very small, and the small ligand field contributions therefore play a crucial role as well. The important result of the present study is that the ab initio analysis used leads to excellent agreement with the experimental data and therefore allows for a thorough understanding of the relaxation pathways. Importantly, this will now allow us to predict the SMM behavior of this type of J

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry (OAc)4]PF6 are available as Supporting Information. CCDC files contain the crystallographic data of the structures reported here. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/data_request/ cif. Elemental analyses were obtained from the microanalytical laboratory of the Chemical Institutes of the University of Heidelberg. Syntheses. General. Chemicals were used as supplied. Technical grade solvents were distilled prior to use. The ligand L (2,6bis((bis(pyridin-2-ylmethyl)amino)methyl)-4-methylphenolate) was prepared as described.57,58 General Procedure for [Ln2(L)(OAc)4]PF6. To a solution of L (200 mg, 0.4 mmol) in methanol (15 mL) was added a solution of Ln(OAc)3 (0.8 mmol) in methanol/water (1:1, 10 mL). The solution was stirred at 55 °C for 1 h. NH4PF6 (160 mg, 1 mmol) was added, and the solution was stirred for an additional 30 min and was then left for evaporation in air. [Y2(L)(OAc)4]PF6. After 1 day, colorless crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Y2(L)(OAc)4]PF6·H2O (%): C, 44.50; H, 4.28; N, 7.59. Found: C, 44.30; H, 4.21; N, 7.70. HR ESI-MS m/z: 915.1439 ([Y2(L)(OAc)3(OMe)]+). [Nd2(L)(OAc)4]PF6. After 1 day, yellow crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Gd2(L)(OAc)4]PF6 (%): C, 40.06; H, 4.02; N, 6.84. Found: C, 39.38; H, 4.42; N, 10.55. HR ESIMS m/z: 1080.1789; calcd 1080.1724 ([Gd2L(OAc)4]+). [Gd2(L)(OAc)4]PF6. After 1 day, yellow crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Gd2(L)(OAc)4]PF6 (%): C, 40.06; H, 4.02; N, 6.84. Found: C, 40.31; H, 3.94; N, 6.65. HR ESIMS m/z: 1080.1789; calcd 1080.1724 ([Gd2(L)(OAc)4]+). [Tb2(L)(OAc)4]PF6. After 1 day, yellow crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Tb2(L)(OAc)4]PF6 (%): C, 40.08; H, 3.69; N, 6.84. Found: C, 39.84; H, 3.94; N, 6.96. HR ESIMS m/z: 1083.1772; calcd 1083.1750 ([Tb2(L)(OAc)4]+). [Dy2(L)(OAc)4]PF6. After 1 day, yellow crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Dy2(L)(OAc)4]PF6 (%): C, 39.85; H, 3.67; N, 6.46. Found: C, 39.62; H, 3.71; N, 6.46. HR ESIMS m/z: 1090.1859; calcd 1090.1804 ([Dy2(L)(OAc)4]+). [Ho2(L)(OAc)4]PF6. After 1 day, red crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Ho2(L)(OAc)4]PF6· NH4PF6 (%): C, 34.82; H, 4.18; N, 6.46. Found: C, 34.95; H, 4.05; N, 6.30. HR ESI-MS m/z: 1027.1885 ([Ho2(L)(OAc)2(OH)(MeO)]+). [Er2(L)(OAc)4]PF6. After 1 day, red crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Er2(L)(OAc)4]PF6 (%): C, 39.54; H, 3.64; N, 6.75. Found: C, 39.65; H, 3.70; N, 6.62. HR ESIMS m/z: 1055.1825 ([Er2(L)(OAc)3(OH)]+). [Lu2(L)(OAc)4]PF6. After 1 day colorless crystals suitable for X-ray diffraction were obtained. Anal. Calcd for [Lu2(L)(OAc)4]PF6· 2NH4PF6 (%): C, 30.96; H, 3.61; N, 7.04. Found: C, 30.45; H, 3.56; N, 7.02. HR ESI-MS m/z: 1087.2141 ([Lu 2 (L)(OAc)3(OMe)]+).



Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Financial support by the German Science Foundation (DFG), the Deutscher Akademischer Austauschdienst (DAAD), the Graduate School HGS MathComp of the IWR, and the University of Heidelberg are gratefully acknowledged.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.5b01673. Crystallographic data (CIF) Details of the crystallographic experiments, the complete list of magnetic data, details of the ligand field analysis, the MCD data, and quantum-chemical analysis (PDF)



REFERENCES

(1) Caneschi, A.; Gatteschi, D.; Sessoli, R.; Barra, A. L.; Brunel, L. C.; Guillot, M. J. Am. Chem. Soc. 1991, 113, 5873−5874. (2) Sessoli, R.; Tsai, H.-L.; Schake, A. R.; Wang, S.; Vincent, J. B.; Folting, K.; Gatteschi, D.; Christou, G.; Hendrickson, D. N. J. Am. Chem. Soc. 1993, 115, 1804−1816. (3) Sessoli, R.; Gatteschi, D.; Caneschi, A.; Novak, M. A. Nature 1993, 365, 141−143. (4) Arom, G.; Aubin, S. M. J.; Bolcar, M. A.; Christou, G.; Eppley, H. J.; Folting, K.; Hendrickson, D. N.; Huffman, J. C.; Squire, R. C.; Tsai, H.-L.; Wang, S.; Wemple, M. W. Polyhedron 1998, 17, 3005−3020. (5) Gatteschi, D.; Sessoli, R. Angew. Chem., Int. Ed. 2003, 42, 268− 297. (6) Long, J. R. Molecular Cluster Magnets. In Molecular Cluster Magnets in Chemistry of Nanostructured Materials; Yang, P., Ed.; World Scientific: Hong Kong, 2003; pp 291−315. (7) Glaser, T. Chem. Commun. 2011, 47, 116−130. (8) Atanasov, M.; Comba, P.; Hausberg, S.; Martin, B. Coord. Chem. Rev. 2009, 253, 2306−2314. (9) Neese, F.; Pantazis, D. A. Faraday Discuss. 2011, 148, 229−238. (10) Comba, P.; Kerscher, M. Coord. Chem. Rev. 2009, 253, 564− 574. (11) Aromi, G.; Brechin, E. K. In Struct. Bonding (Berlin); SpringerVerlag: Berlin, 2006; Vol. 122, pp 1−67. (12) Rinehart, J. D.; Long, J. R. Chem. Sci. 2011, 2, 2078−2085. (13) Sorace, L.; Benelli, C.; Gatteschi, D. Chem. Soc. Rev. 2011, 40, 3092−3104. (14) Guo, Y.-N.; Xu, G.-F.; Guo, Y.; Tang, J. Dalton Trans. 2011, 40, 9953−9963. (15) Hoeke, V.; Stammler, A.; Bögge, H.; Schnack, J.; Glaser, T. Inorg. Chem. 2014, 53, 257−268R. (16) Sessoli, R.; Powell, A. K. Coord. Chem. Rev. 2009, 253, 2328− 2341. (17) Guo, Y.-N.; Chen, X.-H.; Xue, S.; Tang, J. Inorg. Chem. 2011, 50, 9705−9713. (18) Ceulemans, A.; Chibotaru, L. F.; Heylen, G. A.; Pierloot, K.; Vanquickenborne, L. G. Chem. Rev. 2000, 100, 787−806. (19) Chibotaru, L. F.; Ungur, L. J. Chem. Phys. 2012, 137, 064112-1− 064112-22. (20) Singh, S. K.; Rajaraman, G. Chem. - Eur. J. 2014, 20, 113−123. (21) Atanasov, M.; Busche, C.; Comba, P.; El Hallak, F.; Martin, B.; Rajaraman, G.; van Slageren, J.; Wadepohl, H. Inorg. Chem. 2008, 47, 8112−8125. (22) Atanasov, M.; Comba, P.; Helmle, S.; Müller, D.; Neese, F. Inorg. Chem. 2012, 51, 12324−12335. (23) Malmqvist, P.-A.; Roos, B. O.; Schimmelpfennig, B. Chem. Phys. Lett. 2002, 357, 230−240. (24) Lines, M. E. J. Chem. Phys. 1971, 55, 2977−2984. (25) Rudowicz, C. J. Phys. C: Solid State Phys. 1985, 18, 1415−1430. (26) Atanasov, M.; Comba, P.; Förster, S.; Linti, G.; Malcherek, T.; Miletich, R.; Prikhod'ko, A.; Wadepohl, H. Inorg. Chem. 2006, 45, 7722−7735. (27) Miyashita, Y.; Sanada, M.; Islam, M. M.; Amir, N.; Koyano, T.; Ikeda, H.; Fujisawa, K.; Okamoto, K. Inorg. Chem. Commun. 2005, 8, 785−788. (28) Li, Y.; Zheng, F.-K.; Liu, X.; Zou, W.-Q.; Guo, G.-C.; Lu, C.-Z.; Huang, J.-S. Inorg. Chem. 2006, 45, 6308−6316. (29) Niu, S. Y.; Jin, J.; Jin, X. L.; Yang, Z. Z. Solid State Sci. 2002, 4, 1103−1106.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Fax: +49-6221546617. K

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry (30) Costes, J.-P.; Dahan, F.; Dupuis, A.; Lagrave, S.; Laurent, J.-P. Inorg. Chem. 1998, 37, 153−155. (31) Luis, F.; Bartolome, J.; Fernandez, J. F.; Tejada, J.; Hernandez, J. M.; Zhang, X. X.; Ziolo, R. Phys. Rev. B: Condens. Matter Mater. Phys. 1997, 55, 11448−11456. (32) Abragam, A.; Bleaney, B. Electron Paramagnetic Resonance of transition Ions; Clarendon Press: Oxford, 1970. (33) An easy plane-type system can also have a maximal absorption appearing in the lower field region of the resonance at low temperature. However, the simulation of ZFS induced by the magnetic dipolar coupling for an easy plane-type system with g ∼ 2.01 yields 2.5 GHz, which is not confirmed by the experimental results. (34) Benchini, A.; Gatteschi, D. EPR of Exchange Coupled Systems; Dover: Mineola, NY, 2012. (35) Smith, T. D.; Pilbrow, J. R. Coord. Chem. Rev. 1974, 13, 173− 278. (36) Baggio, R.; Calvo, R.; Garland, M. T.; Pena, O.; Perec, M.; Rizzi, A. Inorg. Chem. 2005, 44, 8979−8987. (37) Yeung, Y.-Y.; Tanner, P. A. J. Alloys Compd. 2013, 575, 54−60. (38) Rana, R. S.; Shertzer, J.; Kaseta, F. W.; Garvey, R.; Rana, D.; Feng, S. Y. J. Chem. Phys. 1988, 88, 2242−2248. (39) Görller-Walrand, C.; Binnemans, K. In Handbook on the Physics and Chemistry of Rare Earths; Gschneidner, K. A., Eyring, L., Eds.; Elsevier: New York, 1996; Vol. 23, pp 121−283. (40) Wybourne, B. G. Spectroscopic Properties of Rare Earths; Interscience: New York, 1965. (41) Judd, B. R. Operator Techniques in Atomic Spectroscopy; University Press: Princeton, NJ, 1998. (42) Gerloch, M.; McMeeking, R. F. J. Chem. Soc., Dalton Trans. 1975, 2443−2451. (43) Görller-Walrand, C.; Binnemans, K. Spectral Intensities of f-f Transitions. In Handbook on the Physics and Chemistry of Rare Earths; Gschneidner, K. A., Eyring, L., Ed.; Elsevier: New York, 1998; Vol. 25, pp 101−264. (44) Aquilante, F.; De Vico, L.; Ferré, N.; Ghigo, G.; Malmqvist, P.A.; Neogrády, P.; Pedersen, T. B.; Pitonak, M.; Reiher, M.; Roos, B. O.; Serrano-Andrés, L.; Urban, M.; Veryazov, V.; Lindh, R. J. Comput. Chem. 2010, 31, 224−247. (45) Roos, B. O.; Lindh, R.; Malmqvist, P.-A.; Veryazov, V.; Widmark, P.-O. J. Phys. Chem. A 2005, 109, 6575−6579. (46) Ungur, L.; Van den Heuvel, W.; Chibotaru, L. F. New J. Chem. 2009, 33, 1224−1230. (47) Roos, B. O.; Taylor, P. R.; Siegbahn, P. E. M. Chem. Phys. 1980, 48, 157−173. (48) Chibotaru, L. F.; Ungur, L.; Soncini, A. Angew. Chem., Int. Ed. 2008, 47, 4126−4129. (49) Chibotaru, L. F.; Ungur, L.; Aronica, C.; Elmoll, H.; Pilet, G.; Luneau, D. J. Am. Chem. Soc. 2008, 130, 12445−12455. (50) Ungur, L.; Chibotaru, L. F.; Thewissen, M.; Costes, J.-P.; Wernsdorfer, W. Inorg. Chem. 2013, 52, 6328−6337. (51) Ungur, L.; Chibotaru, L. F. Phys. Chem. Chem. Phys. 2011, 13, 20086−20090. (52) Friedman, J. R.; Sarachik, M.; Tejada, J.; Ziolo, R. Phys. Rev. Lett. 1996, 76, 3830−3833. (53) Leuenberger, M. N.; Loss, D. Phys. Rev. B: Condens. Matter Mater. Phys. 2000, 61, 1286−1302. (54) Ungur, L.; Le Roy, J. J.; Korobkov, I.; Murugesu, M.; Chibotaru, L. F. Angew. Chem., Int. Ed. 2014, 53, 4413−4417. (55) Golze, C.; Alfonsov, A.; Klingeler, R.; Büchner, B.; Kataev, V.; Mennerich, C.; Klauss, H.-H.; Goiran, M.; Broto, J.-M.; Rakoto, H.; Demeshko, S.; Leibeling, G.; Meyer, F. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 73, 224403-1−224403-8. (56) McCaffery, A. J.; Mason, S. F. Mol. Phys. 1963, 6, 359−371. (57) Suzuki, M.; Kanatomi, H.; Murase, I. Chem. Lett. 1981, 10, 1745−1748. (58) Borovik, A. S.; Papaefthymiou, V.; Taylor, L. F.; Anderson, O. P.; Que, L. J. Am. Chem. Soc. 1989, 111, 6183−6195.

L

DOI: 10.1021/acs.inorgchem.5b01673 Inorg. Chem. XXXX, XXX, XXX−XXX