Abstract
Combining strong magnetic anisotropy with strong exchange interaction is a long standing goal in the design of quantum magnets. The lanthanide complexes, while exhibiting a very strong ionic anisotropy, usually display a weak exchange coupling, amounting to only a few wavenumbers. Recently, an isostructural series of mixed (Lnâ=âGd, Tb, Dy, Ho, Er) have been reported, in which the exchange splitting is estimated to reach hundreds wavenumbers. The microscopic mechanism governing the unusual exchange interaction in these compounds is revealed here by combining detailed modeling with density-functional theory and ab initio calculations. We find it to be basically kinetic and highly complex, involving non-negligible contributions up to seventh power of total angular momentum of each lanthanide site. The performed analysis also elucidates the origin of magnetization blocking in these compounds. Contrary to general expectations the latter is not always favored by strong exchange interaction.
Similar content being viewed by others
Introduction
The effects of strong magnetic anisotropy, traditionally investigated in magnetic insulators, especially, in f-electron systems1,2,3, recently attracted renewed interest in connection with molecular magnetic materials4. The investigation of molecular nanomagnets gave birth to new objects such as single-molecule magnets (SMMs)5,6 and single-chain magnets7 and initiated studies in the domain of molecular spintronics8,9 and quantum computation10,11,12. Among them, in the last years the accent moved towards lanthanide complexes which have already demonstrated several exciting properties13,14,15,16,17,18.
The key feature of lanthanide ions in materials is their strong magnetic anisotropy caused by strong spin-orbit coupling effects19, which often leads to highly axial ground and low-lying excited doublet states even in the lack of axial symmetry20. Due to small radius of electronic f-shells, the exchange interaction in lanthanide complexes is much weaker than the crystal-field splitting on lanthanide ions21. As a result, only individual doublet states on lanthanide sites, described by pseudospins , participate in the magnetic interaction. The latter is described by a Hamiltonian bilinear in pseudospins ( and ) in the case of two interacting lanthanide ions, or a pseudospin and a true spin (S2) in the case of a lanthanide ion interacting with a transition metal or a radical when the spin-orbit coupling in the second site is negligible. For strongly axial doublet states on the lanthanide sites (Ln) these Hamiltonians basically become of Ising type22:
either collinear (z1âââz2) or non-collinear depending on geometry21 and details of interaction23. The exchange parameter is contributed by magnetic dipolar and exchange interaction between the sites, , the former being usually stronger in net lanthanide complexes21.
This paradigm was recently challenged by a series of -radical bridged dilanthanide complexes [K(18âââcrownâââ6)]{[(Me3Si)2N]â(THF)Ln}2 (μâââη2:η2âââN2) (Lnâ=âGd (1), Tb (2), Dy (3), Ho (4), Er (5), THFâ=âtetrahydrofuran), shown in Fig. 1a14,24. In some of these compounds the exchange interaction was found to be two orders of magnitude stronger than in any known lanthanide system. This is of the same order of magnitude as the crystal-field splitting of J-multiplets on the lanthanide sites, implying that the picture of exchange interaction involving individual crystal-field doublets, Eq. (1), is no longer valid for these compounds. Moreover, the terbium complex from this series exhibits a magnetic hysteresis at 14âK and a 100âs blocking time at 13.9âK (one of the highest blocking temperatures among existing SMMs24), suggesting a possible implication of the giant exchange interaction in this SMM behavior.
The purpose of the present work is to reveal the mechanism of giant exchange interaction and the origin of the magnetization blocking of the series of the complexes based on adequate theoretical treatment. We apply an approach combining ab initio and density-functional theory (DFT) calculations with microscopic model description to unravel the nature of this exchange interaction. We also elucidate the origin of blocking barriers in these compounds and discuss the effect of strength of exchange interaction on magnetization blocking in strongly anisotropic complexes.
Results
Origin of giant exchange interaction
To understand the origin of such strong exchange interaction, we consider the simplest complex of the series, the gadolinium one. In this system the isotropic spins of Gd3+ ions (SGdâ=â7/2) interact with the radical spin of the bridge via Heisenberg exchange interaction, , described by a single parameter due to the inversion symmetry of the complex (Fig. 1a). Broken-symmetry DFT calculations25 give the value in close agreement with the experimental one, 14 and the previous DFT calculations26,27.
To get insight into the mechanism responsible for the obtained huge value of , we projected a series of DFT calculations into the effective tight-binding and Hubbard models acting in the space of interacting magnetic orbitals of two Gd ions and the radical (see the Supplemental Material for details). Because of the D2h symmetry of the exchange core (Fig. 2a), the antibonding Ï* orbital accommodating the unpaired electron of radical overlaps with only one of the 4f orbitals on each Ln site, the xyz one (Fig. 2b). The corresponding transfer parameter t was derived for the Gd complex as tâ=â1407âcmâ1. The value of t is obtained large because the radicalâs magnetic orbital Ï* resides on nearest-neighbor atoms (nitrogens) to both lanthanides. Most important, this orbital is found to lie higher than the 4fxyz orbitals by as much as Îâ=â5.2âÃâ104âcmâ1 (Fig. 2b). Because of this huge energy gap, small electron promotion energy is expected for the electron transfer from the Ï* to the 4fxyz orbitals: the Coulomb repulsion energy between the transferred electron and the f electrons is cancelled at large extent by Î. On the other hand, because of the same large gap Î, the promotion energy of electron transfer from 4f to Ï* orbital is at least one order of magnitude larger. Therefore, the contribution of this process to the exchange coupling can be neglected. Indeed, our analysis using the Hubbard model gives the experimental for the Gd complex with (averaged) promotion energy of , a value many times smaller than typical âHubbard Uâ in metal complexes28. Taking into account only the dominant virtual electron transfer, (4f)7â(Ï*)1âââ(4f)8â(Ï*)0âââ(4f)7â(Ï*)1, the kinetic contribution to the Gd3+-exchange parameter is written in a good approximation as 29,30.
Compared to this mechanism, the other contributions such as the direct exchange, the delocalization of unpaired electron of into the empty 5d orbitals of Gd3+ (Goodenoughâs mechanism31), the spin polarization and the magnetic dipolar interaction between Gd3+ ions are expected to be 1â2 orders of magnitude smaller. The reason is that all these contributions are expected to be of the same order of magnitude as in other lanthanide-radical compounds. Indeed, the direct exchange integral depends only on the shape of the 4f an radicalâs orbitals, which is not expected to be much different from other complexes. The Goodenoughâs contribution arises from higher (third) order of the perturbation theory compared to the usual kinetic exchange and involves the excitation energy into a higher 5d orbital on the Ln site. Both these contributions are usually neglected unless the conventional kinetic exchange appears to be small29,30. The spin polarization mechanism starts to play a role when the ligand bridging the magnetic centers contains a spectrum of low-lying orbital excitations, which is certainly not the case of . As for magnetic dipolar interaction, it is estimated for Gd3+- to be ~0.25âcmâ1.
The same physical situation is realized in the other complexes of the series. As Table 1 shows, the transfer parameters only slightly decrease with the increase of Ln atomic number. On the other hand, the gap Î between the 4f and the Ï* orbital levels is obtained as huge as in the Gd complex (Table 1), leading again to small promotion energy and, consequently, to the dominant role of the kinetic mechanism in the Ln3+- exchange coupling of complexes 2â5. Given the small change of t through 1â5, the strong variation of the strength of exchange interaction in this series of complexes, testified by the experimental magnetic susceptibilities (Fig. 1b), is expected to be due to the variation of the promotion energy.
Anisotropic exchange interaction
Contrary to the Gd complex, the other members of the series are characterized by strong magnetic anisotropy on the Ln sites induced by the crystal-field (CF) splitting of their atomic J multiplets. These CF-split multiplets are described by multi-configurational wave-functions, therefore, they should be treated by explicitly correlated ab initio approaches32,33 rather than DFT. The ab initio fragment calculations show that the CF split J multiplet on Tb3+ ions (â700âcmâ1) is of the same order of magnitude as the estimated isotropic exchange splitting in 1 (â400âcmâ1). Therefore, in sharp contrast with the common situation in lanthanides, the exchange coupling in the anisotropic 2â5 does not reduce to the interaction between individual (lowest) CF doublets on Ln sites with the Sâ=â1/2 spin of the radical, Eq. (1), but will intermix the entire CF spectrum arising from the ground atomic J multiplet at lanthanide ions. Then, such exchange interaction should be formulated in terms of the total angular momenta (iâ=â1,â2) on the lanthanide sites.
Extending the Andersonâs superexchange theory29,30 to strong spin-orbit coupled systems, the tensorial form of the kinetic (covalent) interaction has been recently derived from the microscopic electronic Hamiltonian34. The kinetic interaction between the lanthanide and radical centers contains besides the exchange part also the Ln3+- covalent contribution (arising from Ln3+- electron delocalization) to the CF splitting at the Ln3+ sites :
Here, and are the unit and the spin operators, respectively, of the radicalâs spin Sâ=â1/2, are the Stevens operators35 of rank k and component q and and are the exchange parameters34. The Stevens operator is a polynomial of (αâ=âx,ây,âz) of kth degree, in which |q| (=±q) corresponds to the order of . The maximal rank of k is 7 for the considered Ln3+ ions, whereas the maximal |q| is 5 in the present case because only the 4fxyz magnetic orbitals at the lanthanide sites contribute to the kinetic exchange. The summation over k in Eqs (2) and (3) is confined to even and odd ranks, respectively, which is required by the invariance of these Hamiltonians with respect to time-inversion. As it is seen from the form of these Hamiltonians, only contributes to the CF splitting of J multiplets on individual metal sites, whereas describes the interaction between powers of total angular momenta at the metal sites with components of spin Sâ=â1/2 of the radical. For comparison, the weak anisotropic exchange interaction between two spins (pseudospins) is described by the exchange Hamiltonian , where D is the 3âÃâ3 exchange matrix, containing one isotropic, five symmetric anisotropic and three antisymmetric (Dzyaloshinsky-Moriya) exchange parameters36. This Hamiltonian corresponds to the first rank contribution (kâ=â1) in Eq. (3), where are just the nine components of the above exchange matrix D. The expression for the exchange parameter 34 includes all virtual electron transfer processes, (4f)nâ(Ï*)1âââ(4f)n+1â(Ï*)0âââ(4f)nâ(Ï*)1, where n is the number of 4f electrons in Ln3+. The multiplet electronic structure of Ln2+ is fully included in the electron promotion energy U0â+âÎEα and the wave functions of the intermediate states, where by U0 we further denote the smallest promotion energy, α numbers the intermediate J-multiplets and ÎEα is the excitation energy of the multiplet α with respect to the ground one in Ln2+.
The highly complex tensorial form of the exchange Hamiltonian is inevitable for orbitally degenerate systems with strong spin-orbit coupling, as was pointed out long time ago37,38. Although all exchange parameters are in principle required for adequate description of the exchange interaction, it is hardly possible to extract a sufficient large number of them from experiment in a unique way. However, once are expressed via microscopic electronic parameters34, the latter can be determined from up-to-date quantum chemistry calculations. Thus the transfer parameter t is obtained here from DFT calculations, expected to be accurate enough39,40, whereas the excitation energies ÎEα and the CF states are obtained by fragment state-of-the-art ab initio calculations including spin-orbit coupling32,33. The only parameter that might be inaccurate when extracted from DFT or ab initio calculations is U0. Indeed, the former gives at most an averaged value over multiplets and the latter systematically overestimates it due to insufficient account of dynamical correlation.
In this way we construct the full microscopic Hamiltonian, , containing only one unknown parameter U0, where is the ab initio CF Hamiltonian for mononuclear Ln fragments (see Supplemental Materials for details). Diagonalizing this Hamiltonian, the magnetic susceptibility Ï for the entire series of compounds has been simulated as described elsewhere33. Figure 1b shows that the experiment is well reproduced for the values of minimal promotion energy U0 listed in Table 1. The calculated exchange parameters for the series of the complexes are shown in Table 2. We can see from the table that the exchange interaction involves non-negligible contributions up to the rank kâ=â7.
The low-lying exchange spectrum for the Tb complex is shown in Fig. 3a. The ground (1±) and the first two excited (2±, 3±) exchange Kramers doublets (KDs) mainly originate from the ground CF doublets on the Tb ions (94%, 87% and 88%, respectively). However, the third and fourth excited exchange KDs (4±, 5±) represent almost equal mixtures of the ground and the first excited CF doublets on the Tb3+ sites. This is remarkable because the mixed CF states are separated by 166âcmâ1 (Fig. 3a). Similar scenario is realized in 3 and 4, whereas in 5 the exchange interaction and the resulting mixing of CF states is relatively weak. The magnetic structure of the ground exchange KD is shown in Fig. 1a. The magnetic moments on Tb3+ sites are parallel due to inversion symmetry and almost coincide with the directions of the main magnetic axes in the ground local KDs (Fig. 1a). The magnetic moment of the radical, corresponding to isotropic Sâ=â1/2, is rotated with respect to the magnetic moments on Tb sites by small angle θ (Table 1) due to the non-Heisenberg contributions to the exchange interaction23.
One may notice that the dominant first rank term of the exchange interaction is of isotropic Heisenberg type despite the strong spin-orbit coupling in Ln3+ ions (Table 2). This looks surprising because even weak spin-orbit coupling makes the first-rank exchange interaction anisotropic36. The analysis of the expression for the first-rank exchange parameters 34 shows that they are in general of non-Heisenberg type, whereas the present case is the only possible exception (see Supplemental Material). Indeed, the isotropy of the first-rank exchange contribution requires involvement of only f orbitals with the projections mâ=â±2. This can only arise for high symmetry of the exchange bridge (Fig. 2a) and for situations with one single electron transfer path, as in the present case. If any other orbital (or more of them) contribute to the electron transfer, the first-rank exchange interaction becomes strongly anisotropic.
Magnetization blocking barriers
Table 1 shows that the transverse g-factors (gx and gy) in the ground exchange KD, the squares of which characterize the rate of quantum tunneling of magnetization (QTM)21, are the largest for 4 and the smallest for 2 and 3 complexes. This explains why large magnetization hysteresis is seen at low temperatures in the latter two compounds, while not seen at all in the former and only weakly observed in the complex 5â14,24. The path characterizing the activated magnetic relaxation in high-temperature domain is shown for the Tb complex in Fig. 3a by red arrows. The height of the activation barrier Ebarrier corresponds to the first excited exchange KD, because its two components (2± in Fig. 3a) are connected by a large magnetic moment matrix element which causes a large temperature-assisted QTM. Blocking barriers of similar structure (Fig. 3a) arise in 3 and 4, their calculated activation energies comparing well with the experimental ones (Table 1).
The unusually large matrix elements between the ground and the first excited exchange KDs are entirely due to the exchange mixing of the ground and the first excited CF doublets on the Ln sites. Indeed, if one quenches the exchange admixture of excited CF doublets to the ground ones, this matrix element becomes three orders of magnitude smaller (Fig. 3b). Then the activated relaxation will proceed via a higher exchange doublet, thereby doubling the height of the blocking barriers (Fig. 3b). Thus in the case of very strong exchange interaction, which is able to intermix the CF states on Ln sites, the axiality of the ground and excited exchange doublets is diminished dramatically and the blocking barriers do not exceed the energy of the first excited exchange KD. In other words, the strength of exchange interaction after reaching a certain value starts playing a destructive role for the magnetization blocking. Therefore, to exploit the effect of strong exchange interaction for achieving high magnetization blocking, an even stronger axial CF field on the Ln sites, precluding the exchange admixture of excited CF states, seems to be indispensable.
Discussion
The mixed lanthanide complexes 1â5 investigated in this work are unique because they show an exchange interaction up to two orders of magnitude stronger than in conventional lanthanide complexes. Due to such strong exchange interaction, a qualitatively new situation arises when the exchange coupling starts to intermix the CF multiplets on the Ln sites. In all previous lanthanide complexes only the ground CF doublets on Ln sites were involved, which led to conventional Ising-type exchange interactions. In the present case, due to the involvement of all CF doublets belonging to the atomic J-multiplet, the exchange interaction becomes highly complex, requiring a tensorial description and involving many parameters.
By combining DFT and ab initio calculations with the microscopic modeling of the exchange interaction, we were able to unravel the mechanism of giant exchange interaction in these complexes. This exchange interaction is found to be kinetic and highly complex, involving non-negligible contributions up to seventh power of total angular momentum of each Ln site. Based on the calculated exchange states, the mechanism of the magnetization blocking is revealed. Contrary to general expectations the latter is not always favored by strong exchange interaction. The accuracy of our approach is proved by the close reproduction of experimental magnetic susceptibility and magnetization blocking barrier for all investigated compounds.
The theoretical analysis proposed in this work opens the way for the investigation of highly complex exchange interaction in materials with strongly anisotropic magnetic sites. Given the large number of involved exchange parameters and the obvious difficulties of their experimental determination, such an approach can become a powerful tool for the study of magnetic materials of primary interest.
Methods
DFT calculations
All DFT calculations were carried out with ORCA 3.0.0. program41 using the B3LYP functional and SVP basis set. Scalar relativistic effects were taken into account within Douglas-Kroll-Hess Hamiltonian. The isotropic exchange parameter for the complex 1, , was derived by applying the broken-symmetry approach25. The obtained was divided by 2 to account for its overestimation due to the self-interaction error42,43. The 4f and the Ï* orbital levels and the transfer parameters t for all complexes 1â5 were derived by projecting the Kohn-Sham orbitals onto a tight-binding model. The averaged promotion energy for the complex 1 was derived by reproducing the energy difference between the high-spin and the broken-symmetry DFT states with a Hubbard model.
Ab initio calculations
Energies and wave functions of CF multiplets on Ln3+ sites in 1â5 have been obtained from fragment ab initio calculations including the spin-orbit coupling, using the quantum chemistry package Molcas 7.844. The calculations have been done for the experimental geometry of the complexes, in which one of the two Ln3+ ions was replaced by an isovalent closed-shell La3+ ions. The total number of electrons was reduced by unity in order to have a closed-shell electronic configuration on the dinitrogen bridge. To simulate the electrostatic crystal field from the removed radicalâs electron, two point charges of â0.5e were added on the nitrogen atoms. For this structural model of a Ln fragment, the complete active space self-consistent field (CASSCF) approach was used including all seven 4f orbitals of the Ln atom in the active space. The spin-orbit interaction was treated with the module SO-RASSI and the local magnetic properties were calculated with the SINGLE_ANISO module of Molcas45. Exchange energy spectrum and magnetic properties of the investigated polynuclear compounds were calculated using the POLY_ANISO program33,45, modified to treat the general form of exchange interaction, Eqs. (2), (3), within the kinetic exchange mechanism.
For further details, see Supplemental Material.
Additional Information
How to cite this article: Vieru, V. et al. Giant exchange interaction in mixed lanthanides. Sci. Rep. 6, 24046; doi: 10.1038/srep24046 (2016).
References
Santini, P. et al. Multipolar interactions in f-electron systems: The paradigm of actinide dioxides. Rev. Mod. Phys. 81, 807â863 (2009).
Zvezdin, A. K., Matveev, V. M., Mukhin, A. A. & Popov, A. I. Rare Earth Ions in Magnetically Ordered Crystals (Nauka, Moskow, 1985) In Russian.
Gingras, M. J. P. & McClarty, P. A. Quantum spin ice: a search for gapless quantum spin liquids in pyrochlore magnets. Rep. Prog. Phys. 77, 056501 1â26 (2014).
Gatteschi, D., Sessoli, R. & Villain, J. Molecular Nanomagnets (Oxford University Press, Oxford, 2006).
Sessoli, R., Gatteschi, D., Caneschi, A. & Novak, M. A. Magnetic bistability in a metal-ion cluster. Nature 365, 141â143 (1993).
Christou, G., Gatteschi, D., Hendrickson, D. N. & Sessoli, R. Single-molecule magnets. MRS Bulletin 25, 66â71 (2000).
Coulon, C., Miyasaka, H. & Clerac, R. Single-chain magnets: Theoretical approach and experimental systems. In Winpenny, R. (ed.) Single-Molecule Magnets and Related Phenomena vol. 122 of Struct. Bond. 163â206 (2006).
Heersche, H. B. et al. Electron transport through single Mn12 molecular magnets. Phys. Rev. Lett. 96, 206801 (2006).
Bogani, L. & Wernsdorfer, W. Molecular spintronics using single-molecule magnets. Nat. Mater. 7, 179â186 (2008).
Leuenberger, M. N. & Loss, D. Quantum computing in molecular magnets. Nature 410, 789â793 (2001).
Timco, G. A. et al. Engineering the coupling between molecular spin qubits by coordination chemistry. Nat. Nanotech. 4, 173â178 (2009).
Aromi, G., Aguila, D., Gamez, P., Luis, F. & Roubeau, O. Design of magnetic coordination complexes for quantum computing. Chem. Soc. Rev. 41, 537â546 (2012).
Ishikawa, N., Sugita, M., Ishikawa, T., Koshihara, S. & Kaizu, Y. Mononuclear lanthanide complexes with a long magnetization relaxation time at high temperatures: A new category of magnets at the single-molecular level. J. Phys. Chem. B 108, 11265â11271 (2004).
Rinehart, J. D., Fang, M., Evans, W. J. & Long, J. R. Strong exchange and magnetic blocking in -radical-bridged lanthanide complexes. Nat. Chem. 3, 538â542 (2011).
Woodruff, D. N., Winpenny, R. E. P. & Layfield, R. A. Lanthanide single-molecule magnets. Chem. Rev. 113, 5110â5148 (2013).
Ungur, L., Lin, S.-Y., Tang, J. & Chibotaru, L. F. Single-molecule toroics in Ising-type lanthanide molecular clusters. Chem. Soc. Rev. 43, 6894â6905 (2014).
Layfield, R. & Murugesu, M. (eds) Lanthanides and Actinides in Molecular Magnetism (Wiley, New Jersey, 2015).
Demir, S., Jeon, I.-R., Long, J. R. & Harris, T. D. Radical ligand-containing single-molecule magnets. Coord. Chem. Rev. 289, 290, 149â176 (2015).
Dieke, G. H. Spectra and Energy Levels of Rare-Earth Ions in Crystals (Academic Press Inc., New York, 1967).
Ungur, L. & Chibotaru, L. F. Magnetic anisotropy in the excited states of low symmetry lanthanide complexes. Phys. Chem. Chem. Phys. 13, 20086â20090 (2011).
Chibotaru, L. F. Theoretical understanding of Anisotropy in Molecular Nanomagnets. In Gao, S. (ed.) Molecular Nanomagnets and Related Phenomena vol. 164 of Struct. Bond. 185â229 (Springer: Berlin Heidelberg, 2015).
Chibotaru, L. F., Ungur, L. & Soncini, A. The origin of nonmagnetic Kramers doublets in the ground state of dysprosium triangles: evidence for a toroidal magnetic moment. Angew. Chem. Int. Ed. 120, 4194â4197 (2008).
Chibotaru, L. F. & Iwahara, N. Ising exchange interaction in lanthanides and actinides. New J. Phys. 17, 103028 1â15 (2015).
Rinehart, J. D., Fang, M., Evans, W. J. & Long, J. R. A -Radical-Bridged terbium complex exhibiting magnetic hysteresis at 14 K. J. Am. Chem. Soc. 133, 14236â14239 (2011).
Soda, T. et al. Ab initio computations of effective exchange integrals for HâH, HâHeâH and Mn2O2 complex: comparison of broken-symmetry approaches. Chem. Phys. Lett. 319, 223â230 (2000).
Rajeshkumar, T. & Rajaraman, G. Is a radical bridge a route to strong exchange interactions in lanthanide complexes? A computational examination. Chem. Commun. 48, 7856â7858 (2012).
Zhang, Y.-Q., Luo, C.-L., Wang, B.-W. & Gao, S. Understanding the Magnetic Anisotropy in a Family of Radical-Bridged Lanthanide Complexes: Density Functional Theory and ab Initio Calculations. J. Phys. Chem. A 117, 10873â10880 (2013).
van der Marel, D. & Sawatzky, G. A. Electron-electron interaction and localization in d and f transition metals. Phys. Rev. B 37, 10674â10684 (1988).
Anderson, P. W. New approach to the theory of superexchange interactions. Phys. Rev. 115, 2â13 (1959).
Anderson, P. W. Theory of magnetic exchange interactions: exchange in insulators and semiconductors. In Seitz, F. & Turnbull, D. (eds.) Solid State Physics vol. 14, 99â214 (Academic Press, New York, 1963).
Goodenough, J. B. Magnetism and the Chemical Bond (John Wiley & Sons, New York, 1963).
Chibotaru, L. F. & Ungur, L. Ab initio calculation of anisotropic magnetic properties of complexes. I. Unique definition of pseudospin Hamiltonians and their derivation. J. Chem. Phys. 137, 064112 1â22 (2012).
Ungur, L. & Chibotaru, L. F. Computational Modelling of the Magnetic Properties of Lanthanide Compounds. In Layfield, R. & Murugesu, M. (eds.) Lanthanides and Actinides in Molecular Magnetism 153â184 (Wiley, New Jersey, 2015).
Iwahara, N. & Chibotaru, L. F. Exchange interaction between J multiplets. Phys. Rev. B 91, 174438 1â18 (2015).
Stevens, K. W. H. Matrix Elements and Operator Equivalents Connected with the Magnetic Properties of Rare Earth Ions. Proc. Phys. Soc. London, Sec. A 65, 209â215 (1952).
Moriya, T. Anisotropic Superexchange Interaction and weak Ferromagnetism. Phys. Rev. 120, 91â98 (1960).
Elliott, R. J. & Thorpe, M. F. Orbital effects on exchange interactions. J. Appl. Phys. 39, 802â807 (1968).
Hartmann-Boutron, F. Interactions de superéchange en présence de dégénérescence orbitale et de couplage spin-orbite. J. Phys. France 29, 212â214 (1968).
Imada, M. & Miyake, T. Electronic structure calculation by first principles for strongly correlated electron systems. J. Phys. Soc. Jpn. 79, 112001 1â42 (2010).
Solovyev, I. V. Combining dft and many-body methods to understand correlated materials. J. Phys.: Condens. Matter 20, 293201 1â33 (2008).
Neese, F. The ORCA program system. WIREs Comput. Mol. Sci. 2, 73â78 (2012).
Polo, V., Gräfenstein, J., Kraka, E. & Cremer, D. Long-range and short-range Coulomb correlation effects as simulated by HartreeâFock, local density approximation and generalized gradient approximation exchange functionals. Theor. Chem. Acc. 109, 22â35 (2003).
Ruiz, E., Alvarez, S., Cano, J. & Polo, V. About the calculation of exchange coupling constants using density-functional theory: The role of the self-interaction error. J. Chem. Phys. 123, 164110 1â7 (2005).
Aquilante, F. et al. Molcas 7: The next generation. J. Comput. Chem. 31, 224â247 (2010).
Chibotaru, L. F. & Ungur, L. The computer programs SINGLE_ANISO and POLY_ANISO. University of Leuven (2006).
Acknowledgements
N.I. would like to acknowledge the financial support from the the Fonds Wetenschappelijk Onderzoek - Vlaanderen (FWO) and the GOA grant from Katholieke Universiteit (KU) Leuven. L.U. is a postdoctoral fellow of FWO and also acknowledges INPAC and Methusalem grands of KU Leuven.
Author information
Authors and Affiliations
Contributions
V.V. did DFT and CASSCF calculations, N.I. developed the methodology to derive the tight-binding and Hubbard parameters, L.U. implemented the methodology of calculation of anisotropic multi-rank exchange interaction in the module Poly_Aniso and L.F.C. conceived and supervised the work. All authors participated to the discussion of the results and the writing of the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Vieru, V., Iwahara, N., Ungur, L. et al. Giant exchange interaction in mixed lanthanides. Sci Rep 6, 24046 (2016). https://doi.org/10.1038/srep24046
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep24046
This article is cited by
-
Electrically driven spin resonance of 4f electrons in a single atom on a surface
Nature Communications (2024)
-
Computational design of magnetic molecules and their environment using quantum chemistry, machine learning and multiscale simulations
Nature Reviews Chemistry (2022)
-
Air-stable redox-active nanomagnets with lanthanide spins radical-bridged by a metalâmetal bond
Nature Communications (2019)
-
Giant coercivity and high magnetic blocking temperatures for N2 3â radical-bridged dilanthanide complexes upon ligand dissociation
Nature Communications (2017)
-
Single molecule magnet with an unpaired electron trapped between two lanthanide ions inside a fullerene
Nature Communications (2017)