Strong magnetic frustration and anti-site disordercausing spin-glass behavior in honeycomb Li{}_{2}RhO{}_{3}

Strong magnetic frustration and anti-site disorder
causing spin-glass behavior in honeycomb LiRhO

Vamshi M. Katukuri Satoshi Nishimoto Ioannis Rousochatzakis Institute for Theoretical Solid State Physics, IFW Dresden, Helmholtzstr. 20, 01069 Dresden, Germany Hermann Stoll Institute for Theoretical Chemistry, Universität Stuttgart, Pfaffenwaldring 55, 70550 Stuttgart, Germany Jeroen van den Brink Institute for Theoretical Solid State Physics, IFW Dresden, Helmholtzstr. 20, 01069 Dresden, Germany Liviu Hozoi Institute for Theoretical Solid State Physics, IFW Dresden, Helmholtzstr. 20, 01069 Dresden, Germany
July 12, 2019

With large spin-orbit coupling, the electron configuration in -metal oxides is prone to highly anisotropic exchange interactions and exotic magnetic properties. In iridates, given the existing variety of crystal structures, the magnetic anisotropy can be tuned from antisymmetric to symmetric Kitaev-type, with interaction strengths that outsize the isotropic terms. By many-body electronic-structure calculations we here address the nature of the magnetic exchange and the intriguing spin-glass behavior of LiRhO, a honeycomb oxide. For pristine crystals without Rh-Li site inversion, we predict a dimerized ground state as in the isostructural iridate LiIrO, with triplet spin dimers effectively placed on a frustrated triangular lattice. With Rh-Li anti-site disorder, we explain the observed spin-glass phase as a superposition of different, nearly degenerate symmetry-broken configurations.


For -metal compounds with localized magnetic moments, basic guidelines to soothsay the sign of the nearest-neighbor (NN) magnetic exchange interactions, i. e., the Anderson-Goodenough-Kanamori rules [1, 2, 3], were laid down back in the 1950’s. With one single bridging anion and half-filled states these rules safely predict antiferromagnetic (AF) exchange interactions, as is indeed encountered in numerous magnetic Mott insulators. It is however much harder to anticipate the sign of the couplings for geometries with two bridging ligands and bond angles close to . Illustrative recent examples are the honeycomb systems NaIrO and LiIrO. The signs of the NN Heisenberg and of the additional symmetric Kitaev anisotropy are intensely debated in these iridates, with both , [4, 5, 6] and , [7, 8, 9, 10, 11, 12] sets of parameters being used to explain the available experimental data.

The sizable Kitaev interactions, that is, uniaxial symmetric terms () that cyclically permute on the bonds of a particular hexagonal ring [13, 14, 15], are associated to strong frustration effects and unconventional magnetic ground states displaying, for example, noncollinear order, incommensurability, or spin-liquid behavior [4, 5, 7, 8, 15, 16, 17, 12, 6, 11]. Obviously, in the context of electronic-structure computational methods, such features cannot be thoroughly addressed by periodic total-energy calculations for a given set of predetermined spin configurations. A much more effective strategy is to first focus on individual pairs of NN -metal sites and obtain reliable values for the associated effective magnetic couplings by using ab initio many-body quantum chemistry (QC) machinery (for a recent review, see Ref. 18). The computed exchange parameters can be subsequently fed to effective spin Hamiltonians to be solved for larger sets of magnetically active lattice sites. Such an approach, earlier, allowed us to establish the signs plus the relative strengths of the Heisenberg and Kitaev interactions in both NaIrO and LiIrO and to additionally rationalize the qualitatively different types of AF orders in these two honeycomb iridates [16, 17].

The related honeycomb compound LiRhO is even more puzzling because it features no sign of long-range magnetic order. Instead, an experimental study suggests the presence of a spin-glass ground state [19]. While the spin-orbit couplings (SOC’s) are still sizable for the shell and may in principle give rise on the honeycomb lattice to compelling Kitaev physics, to date no conclusive evidence is in this respect available for LiRhO. To shed light on the nature of the essential exchange interactions in LiRhO we here carry out detailed ab initio QC calculations. We show that large trigonal splittings within the Rh shell, comparable with the strength of the SOC, dismiss a simple picture based on and effective states [20, 21, 22, 14]. The magnetic properties of the system can still be described, however, in terms of pseudospins. The calculations earmark LiRhO as a 4-electron system with remarkably large anisotropic magnetic couplings, in particular, FM Kitaev interactions of up to 10–15 meV. The isotropic Heisenberg exchange, on the other hand, features opposite signs on the two sets of structurally distinct links of Rh NN’s. This sign modulation of the NN Heisenberg interactions, with strong ferromagnetic (FM) ’s for one type of Rh-Rh bonds and weaker AF couplings for the other pairs of adjacent Rh sites, enables the initial hexagonal network to be mapped onto an effective model of spin-1 dimers on a frustrated triangular lattice. We further address the issue of Rh-Li anti-site disorder in samples of LiRhO. By exact-diagonalization (ED) calculations for an extended spin model that also includes second and third neighbor couplings, we show that the experimentally observed spin-glass behavior can be rationalized as a superposition of different nearly-degenerate symmetry-broken states arising at finite concentration of in-plane spin vacancies.


Rh electronic structure

The tetravalent rhodium ions in LiRhO display a 4 valence electron configuration, octahedral ligand coordination and bonding through two bridging ligands. In the simplest picture, i.e., for sufficiently large Rh splittings and degenerate levels, the ground-state electron configuration at each site is a effective spin-orbit doublet [20, 21, 22, 14]. For ions in a variety of three-dimensional, layered, and chain-like oxides, ab initio QC electronic-structure calculations yield excitation energies of 0.6–0.9 eV for the transitions between the and split levels [23, 24, 25, 26] and indicate values of 0.45–0.5 eV for the strength of the SOC , in agreement with earlier estimates [27]. Sharp features in the range of 0.6–0.9 eV are indeed found in the resonant x-ray scattering spectra [28, 25, 23, 24].

The validity of the approximation for the ground state of LiRhO is however questionable since the SOC is substantially weaker for elements. Indeed our QC calculations (see Table 1) indicate Rh splittings eV, close to values of 0.14–0.16 eV estimated for in various Rh oxides [29, 26]. For the ab initio QC investigation we employed multiconfiguration complete-active-space self-consistent-field (CASSCF) and multireference configuration-interaction (MRCI) calculations [30], see Supplemental Material (SM) and Refs. 26, 23. With and parameters of similar size, the and states are strongly admixed, as discussed in earlier work [20, 31] and illustrated in Table 1. In particular, for the relativistic ground-state wave function the hole is not equally distributed among the three Rh levels as for the “true” ground state [20, 21] but displays predominant character, close to 60%.

Magnetic couplings between two adjacent Rh ions

Interestingly, while the results for the on-site excitations are quite different as compared to the excitation energies [23], the computed intersite effective interactions are qualitatively similar to those obtained for the honeycomb iridate LiIrO [17]. The intersite exchange couplings were estimated by MRCI+SOC calculations for embedded fragments having two edge-sharing IrO octahedra in the active region. As described in earlier work [16, 17, 32] and in SM, the ab initio QC data for the lowest four spin-orbit states describing the magnetic spectrum of two NN octahedra is mapped in our scheme onto an effective spin Hamiltonian including both isotropic Heisenberg exchange and symmetric anisotropy. Yet the spin-orbit calculations, CASSCF or MRCI, incorporate all nine triplet and nine singlet states that arise from the two-site configuration.

For on-site Kramers-doublet configurations, the most general symmetry-allowed form of the effective spin Hamiltonian, for a pair of NN ions, is


where , are 1/2 pseudospin operators, is the isotropic Heisenberg interaction, the Kitaev coupling, and the coeeficients are offdiagonal elements of the symmetric anisotropic exchange matrix with . The antisymmetric anisotropic term vanishes since the crystallographic data reported in Ref. 33 indicate overall point-group symmetry for one block of NN RhO ocathedra, green (B1) bonds in Fig. 1, and only slight deviations from for the other type of NN’s, blue (B2 and B3) bonds in Fig. 1. For symmetry of the Rh-Rh link, . We note that in (1) and stand for components in the local, Kitaev bond reference frame. The axis is here perpendicular to the RhO plaquette (see SI and Refs. 14, 16, 17).

Relative energies for the four low-lying states describing the magnetic spectrum of two NN octahedra and the resulting effective coupling constants are listed in Table 2. For the effective picture of pseudospins assumed in Eq. (1), the set of four eigenfunctions contains the singlet and the triplet components , , . In symmetry, the “full” spin-orbit wave functions associated to , , and transform according to the , , and irreducible representations, respectively. Since two of the triplet terms may interact, the most compact way to express the eigenstates of (1) is then , , and . The angle parametrizes the amount of mixing, related to finite off-diagonal couplings. This degree of admixture is determined by analysis of the full QC spin-orbit wave functions. The effective parameters provided in Table 2 are obtained for each type of Rh-Rh link by using the , , , MRCI relative energies and the mixing coefficients (see SM).

For B1 links, we find that both and are FM. While by MRCI calculations always comes FM in spin-orbit coupled honeycomb systems [16, 17], the FM for the B1 bonds has much to do with the peculiar kind of dependence on the amount of trigonal squashing of the oxygen octahedra and consequently on the variation of the Rh-O-Rh angles of the RhO plaquette. The latter increase to values larger than 90 for finite trigonal compression. This dependence of the NN on the Rh-O-Rh bond angles is illustrated in Fig. 2 for a simplified structural model of LiRhO where the Rh-O bond lengths are all the same, set to the average bond length in the experimental crystal structure [33]. It is seen that displays a parabolic behavior, with a minimum of about –5 meV in the interval 92–93 and a change of sign to AF couplings around 96. For the B1 Rh-Rh links, the Rh-O-Rh bond angle is 93.4, close to the value that defines the minimum in Fig. 2. The difference between the –5 meV minimum of Fig. 2 and the –10 meV result listed in Table 2 comes from additional distortions of the O octahedra in the actual structure (see the footnotes in Table 2 and Ref. 33), not included in the idealized model considered for the plot in Fig. 2. An even stronger FM was computed for the B1 type bonds in the related compound LiIrO [17]. In NaIrO, on the other hand, the Ir-O-Ir bond angles are 97 and the NN turns AF on all short Ir-Ir links [16].

For the B2 and B3 links, we derive a FM Kitaev term and an AF Heisenberg interaction, again qualitatively similar to the QC data for LiIrO [17]. We assign the AF value of the NN on the B2/B3 units to the slightly larger Rh-O-Rh bond angle, which as shown in Fig. 2 pulls the towards a positive value, and most importantly to additional distortions that shift the bridging ligands on the Rh-O-Rh B2/B3 plaquettes in opposite senses parallel to the Rh-Rh axis [33]. The role of these additional distortions on the B2/B3 units was analyzed in detail in Ref. 16 and shown to enhance as well the AF component to the intersite exchange.

Effect of longer-range exchange interactions and occurence of spin-glass ground state

For further insights into the magnetic properties of LiRhO, we carried out ED calculations for an extended spin Hamiltonian that in addition to the NN terms of Eq.(1) incorporates longer-range second- and third-neighbor Heisenberg interactions and [9, 8, 7, 5]. We used clusters of 24 sites with periodic boundary conditions [4, 16, 17] and the quantum chemically derived NN coupling constants listed in Table 2. The static spin-structure factor was calculated as function of variable and parameters. For a given set of and values, the dominant order is determined according to the wave number Q = providing a maximum value of . The resulting phase diagram is shown in Fig. 3(a). Given the similar structure of the NN magnetic interactions, it is somewhat similar to that obtained in our previous study on LiIrO [17]. Six different regions can be identified for meV: FM, Néel, Kitaev spin liquid, stripy, diagonal zigzag and incommensurate (ICx) phases. The Kitaev spin liquid, stripy and incommensurate phases in strongly spin-orbit coupled honeycomb systems were analyzed in a number of earlier studies [4, 8, 7, 6, 16]. The detailed nature of the diagonal zigzag and incommensurate ICx ground states for large FM on one set of NN links was described in Ref. 17. Remarkably, for much larger than and , the initial hexagonal lattice can be mapped onto an effective triangular model of triplet dimers on the B1 bonds [17].

Since and are expected to be AF in honeycomb oxides [9, 8], the most likely candidate for the magnetic ground state of “clean” crystals of LiRhO, according to our results, is the diagonal zigzag state (see Fig. 3) and is found to be stable in a wide region of and . Experimentally, however, a spin-glass ground state was determined, with a spin freezing temperature of  K [19]. As possible cause of the observed spin-glass behavior in LiRhO we here investigate the role of Li-Rh site intercalation. Significant disorder on the cation sublattice is a well known feature in LiMO compounds. A typical value for the degree of LiM site inversion in these materials is 10–15 [34, 35]. Partial substitution of the “in-plane” Rh ions by nonmagnetic Li species introduces spin deffects in the honeycomb layer. On the 24-site cluster employed for our ED calculations, 10–15 site inversion translates in replacing two centers by vacancies. Hereafter, we denote the two spin defects as and .

The effect of spin vacancies on the static spin-structure factor in the diagonal zigzag phase () is shown in Fig. 3(d),(e). For comparison, the static spin-structure factor is also plotted in Fig. 3(c) for the ideal case without spin defects. In the absence of “defects”, the ground state is characterized in the bulk limit by symmetry-broken long-range order with either or wave vectors. Since the two symmetry-breaking states are degenerate, the structure factor displays four peaks, at and [see Fig. 3(c)]. However, by introducing spin vacancies, the degeneracy may be lifted via impurity pinning effects. For example, when the two defects occupy positions 17 and 20 [, see Fig. 3(b)] the spin structure defines one of the symmetry-breaking states with [Fig. 3(d)]; likewise, defects at yield a state with [Fig. 3(e)]. In other words, two different kinds of dominant short-range order can be obtained with anti-site disorder. The “locally” favored symmetry-breaking direction depends on the relative positions of the spin vacancies. In a macroscopic system, such “local” domains displaying different symmetry-breaking ordering directions are randomly distributed. Additional frustration is expected to arise because it is not possible to match two differently ordered domains without an emerging “string”. It is therefore likely that by creating some amount of spin defects the long-range zigzag order disappears and the resulting state is perceived as a spin glass at low temperature. A similar mechanism was proposed for the isotropic Heisenberg-Kitaev and -- models [5].

An early well known example of frustration induced through the competition between two different, degenerate spin configurations is the two-dimensional Ising model on a square lattice with randomly distributed, competing FM and AF bonds [36]. To investigate how the diagonal zigzag state is destroyed by increasing the concentration of spin defects, we also studied a simplified Ising model with for the B1 bonds, and all other interactions set to zero. This is a reasonable approximation for the honeycomb layers of LiRhO since the diagonal zigzag phase essentially consists of alternating spin-up and spin-down chains [see sketch in Fig. 3(a)]. Spin structures obtained this way for various spin-defect concentrations are shown in Fig. 3(f)-(i). For the symmetry-breaking diagonal zigzag state is realized, with degenerate and spin structures. At finite, low concentration those two configurations are no longer degenerate since one of them features slightly lower ground-state energy. We still have in this case a “macroscopically stable” ground state. At intermediate defect concentration the long-range order is in a strict sense destroyed. However, the large domain walls with either or seem to survive. At higher concentration the long-range order disappears completely. Moreover, we can now identify a mixture of local structures with different symmetry-breaking directions [see Fig. 3(i)].


In sum, we have calculated the microscopic neareast-neighbor magnetic interactions between effective 1/2 spins in LiRhO and uncovered a substantial difference between the two types of bonds that are present: one is dominated by Heisenberg and the other by Kitaev types of couplings. The latter give rise to strong frustration, even if the interactions are predominantly ferromagnetic. In this setting we have additionally considered the effect of the presence of anti-site disorder. Experimentally the in-plane spin-defect concentration in LiRhO has been estimated as 10–15 [34, 35]. Based on our theoretical findings it is likely that the observed spin-glass behavior arises from the combination of such anti-site disorder and strongly frustrating magnetic interactions, in particular, the different Kitaev/Heisenberg dominated magnetic bonds and the Ising-like physics associated with the triplet dimer formation that results from there.

Our combined ab initio and effective-model calculations on both LiRhO and related honeycomb iridates [23, 16, 17] indicate that a description in terms of on-site 1/2 pseudospins can well account for the diverse magnetic properties of these systems. While alternative models rely on the formation of delocalized, quasimolecular orbitals [37, 38] and for LiRhO downplay the role of spin-orbit interactions [38], here we show that the latter give rise in LiRhO to anisotropic Kitaev interactions the same magnitude as in iridates [16, 17, 39]. That happens in spite of having a Rh splitting and a spin-orbit coupling of similar magnitude, the same way similar sets of Ir and parameters in CaIrO [40] still generate symmetric anisotropic exchange terms in the range of 10 meV (work is in progress).


The molpro QC package was employed for all ab initio calculations [41]. To analyze the electronic ground state and the nature of the - excitations, a cluster consisting of one reference RhO octahedron plus three NN RhO octahedra and 15 nearby Li ions was used. The magnetic spectrum for two Rh ions was obtained from calculations on a cluster containing two reference and four NN RhO octahedra plus the surrounding 22 Li ions, see SI for details. The farther solid-state environment was in both cases modeled as a finite array of point charges fitted to reproduce the crystal Madelung field in the cluster region. The spin-orbit treatment was carried out according to the procedure described in Ref. 42, using spin-orbit pseudopotentials for Ir.


We thank V. Yushankhai, Y. Singh, N. A. Bogdanov, and U. K. Rößler for useful discussions. L. H. acknowledges financial support from the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG).

Author contributions

V.M.K. carried out the ab initio calculations and subsequent mapping of the ab initio results onto the effective spin Hamiltonian, with assistance from L.H., H.S. and I.R. S.N. performed the exact-diagonalization calculations. V.M.K., S.N., J.v.d.B. and L.H. analyzed the data and wrote the paper, with contributions from all other coauthors.

Additional information

Competing financial interests: Authors have no competing financial interests.


  • 1. Goodenough, J. B. An interpretation of the magnetic properties of the perovskite-type mixed crystals LaSrCoO. J. Phys. Chem. Sol. 6, 287 – 297 (1958).
  • 2. Kanamori, J. Superexchange interaction and symmetry properties of electron orbitals . J. Phys. Chem. Sol. 10, 87 – 98 (1959).
  • 3. Anderson, P. W. New approach to the theory of superexchange interactions. Phys. Rev. 115, 2 (1959).
  • 4. Chaloupka, J., Jackeli, G. & Khaliullin, G. Zigzag Magnetic Order in the Iridium Oxide NaIrO. Phys. Rev. Lett. 110, 097204 (2013).
  • 5. Andrade, E. C. & Vojta, M. Magnetism in spin models for depleted honeycomb-lattice iridates: Spin-glass order towards percolation. Phys. Rev. B 90, 205112 (2014).
  • 6. Rau, J. G., Lee, E. K.-H. & Kee, H.-Y. Generic Spin Model for the Honeycomb Iridates beyond the Kitaev Limit. Phys. Rev. Lett. 112, 077204 (2014).
  • 7. Kimchi, I. & You, Y.-Z. Kitaev-Heisenberg-- model for the iridates IrO. Phys. Rev. B 84, 180407 (2011).
  • 8. Singh, Y. et al. Relevance of the Heisenberg-Kitaev Model for the Honeycomb Lattice Iridates . Phys. Rev. Lett. 108, 127203 (2012).
  • 9. Choi, S. K. et al. Spin Waves and Revised Crystal Structure of Honeycomb Iridate . Phys. Rev. Lett. 108, 127204 (2012).
  • 10. Foyevtsova, K., Jeschke, H. O., Mazin, I. I., Khomskii, D. I. & Valentí, R. Ab initio analysis of the tight-binding parameters and magnetic interactions in NaIrO. Phys. Rev. B 88, 035107 (2013).
  • 11. Reuther, J., Thomale, R. & Rachel, S. Spiral order in the honeycomb iridate . Phys. Rev. B 90, 100405 (2014).
  • 12. Sela, E., Jiang, H.-C., Gerlach, M. H. & Trebst, S. Order-by-disorder and spin-orbital liquids in a distorted Heisenberg-Kitaev model. Phys. Rev. B 90, 035113 (2014).
  • 13. Kitaev, A. Anyons in an exactly solved model and beyond. Ann. Phys. 321, 2 – 111 (2006).
  • 14. Jackeli, G. & Khaliullin, G. Mott Insulators in the Strong Spin-Orbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev Models. Phys. Rev. Lett. 102, 017205 (2009).
  • 15. Chaloupka, J., Jackeli, G. & Khaliullin, G. Kitaev-Heisenberg Model on a Honeycomb Lattice: Possible Exotic Phases in Iridium Oxides AIrO. Phys. Rev. Lett. 105, 027204 (2010).
  • 16. Katukuri, V. M. et al. Kitaev interactions between moments in honeycomb NaIrO are large and ferromagnetic: insights from ab initio quantum chemistry calculations. New J. Phys. 16, 013056 (2014).
  • 17. Nishimoto, S. et al. Strongly frustrated triangular spin lattice emerging from triplet dimer formation in honeycomb LiIrO. arXiv:1403.6698 (unpublished) (2014).
  • 18. Malrieu, J. P., Caballol, R., Calzado, C. J., de Graaf, C. & Guihery, N. Magnetic Interactions in Molecules and Highly Correlated Materials: Physical Content, Analytical Derivation, and Rigorous Extraction of Magnetic Hamiltonians. Chem. Rev. 114, 429 (2014).
  • 19. Luo, Y. et al. LiRhO: A spin-glassy relativistic Mott insulator. Phys. Rev. B 87, 161121 (2013).
  • 20. Thornley, J. H. M. The magnetic properties of (IrX) complexes. J. Phys. C (Proc. Phys. Soc.) 1, 1024 (1968).
  • 21. Abragam, A. & Bleaney, B. Electron Paramagnetic Resonance of Transition Ions (Clarendon Press, Oxford, 1970).
  • 22. Kim, B. J. et al. Novel Mott State Induced by Relativistic Spin-Orbit Coupling in SrIrO. Phys. Rev. Lett. 101, 076402 (2008).
  • 23. Gretarsson, H. et al. Crystal-Field Splitting and Correlation Effect on the Electronic Structure of . Phys. Rev. Lett. 110, 076402 (2013).
  • 24. Liu, X. et al. Testing the Validity of the Strong Spin-Orbit-Coupling Limit for Octahedrally Coordinated Iridate Compounds in a Model System . Phys. Rev. Lett. 109, 157401 (2012).
  • 25. Hozoi, L. et al. Longer-range lattice anisotropy strongly competing with spin-orbit interactions in pyrochlore iridates. Phys. Rev. B 89, 115111 (2014).
  • 26. Katukuri, V. M. et al. Electronic Structure of Low-Dimensional 4 Oxides: Interplay of Ligand Distortions, Overall Lattice Anisotropy, and Spin–Orbit Interactions. Inorg. Chem. 53, 4833–4839 (2014).
  • 27. Andlauer, B., Schneider, J. & Tolksdorf, W. Optical Absorption, Fluorescence, and Electron Spin Resonance of Ir on Octahedral Sites in YGaO. Phys. Stat. Sol. B 73, 533 (1976).
  • 28. Kim, J. et al. Magnetic Excitation Spectra of Probed by Resonant Inelastic X-Ray Scattering: Establishing Links to Cuprate Superconductors. Phys. Rev. Lett. 108, 177003 (2012).
  • 29. Haverkort, M. W., Elfimov, I. S., Tjeng, L. H., Sawatzky, G. A. & Damascelli, A. Strong Spin-Orbit Coupling Effects on the Fermi Surface of SrRuO and SrRhO. Phys. Rev. Lett. 101, 026406 (2008).
  • 30. Helgaker, T., Jørgensen, P. & Olsen, J. Molecular Electronic-Structure Theory (Wiley, Chichester, 2000).
  • 31. Hill, N. J. Electron paramagnetic resonance of osmium-doped trichlorotris(diethylphenylphosphine)rhodium(III). J. Chem. Soc., Faraday Trans. 2 68, 427–434 (1972).
  • 32. Katukuri, V. M. et al. Mechanism of basal-plane antiferromagnetism in the spin-orbit driven iridate BaIrO. Phys. Rev. X 4, 021051 (2014).
  • 33. Todorova, V. & Jansen, M. Synthesis, Structural Characterization and Physical Properties of a New Member of Ternary Lithium Layered Compounds - LiRhO. Z. Anorg. Allg. Chem. 637, 37 (2011).
  • 34. Kobayashi, H. et al. Structure and lithium deintercalation of LiRuO. Solid State Ionics 82, 25–31 (1995).
  • 35. Kobayashi, H., Tabuchi, M., Shikano, M., Kageyama, H. & Kanno, R. Structure, and magnetic and electrochemical properties of layered oxides, LiIrO. J. Mater. Chem. 13, 957–962 (2003).
  • 36. Vannimenusi, J. & Toulouse, G. Theory of the frustration effect: II. Ising spins on a square lattice. J. Phys.: Solid State Phys. 10, L537–L542 (1977).
  • 37. Mazin, I. I., Jeschke, H. O., Foyevtsova, K., Valentí, R. & Khomskii, D. I. as a Molecular Orbital Crystal. Phys. Rev. Lett. 109, 197201 (2012).
  • 38. Mazin, I. I. et al. Origin of the insulating state in honeycomb iridates and rhodates. Phys. Rev. B 88, 035115 (2013).
  • 39. Gretarsson, H. et al. Magnetic excitation spectrum of NaIrO probed with resonant inelastic x-ray scattering. Phys. Rev. B 87, 220407 (2013).
  • 40. Bogdanov, N. A., Katukuri, V. M., Stoll, H., van den Brink, J. & Hozoi, L. Post-perovskite CaIrO: A quasi-one-dimensional antiferromagnet. Phys. Rev. B 85, 235147 (2012).
  • 41. Werner, H.-J., Knowles, P. J., Knizia, G., Manby, F. R. & Schütz, M. Molpro: a general-purpose quantum chemistry program package. WIREs Comput Mol Sci 2, 242–253 (2012).
  • 42. Berning, A., Schweizer, M., Werner, H.-J., Knowles, P. J. & Palmieri, P. Spin-orbit matrix elements for internally contracted multireference configuration interaction wavefunctions. Mol. Phys. 98, 1823–1833 (2000).
states Relative Wave-function composition
energies (normalized weights)
Table 1: Rh states in LiRhO, with composition of the wave functions (hole picture) and relative energies (meV). CASSCF results without and with SOC are shown. Only the three Rh orbitals were active [30] in CASSCF. By subsequent MRCI calculations, the relative energies of these states change to 0, 85, 95 without SOC and 0, 235, 285 meV with SOC included. Only one component of the Kramers’ doublet is shown for each spin-orbit wave function.
Energies & effective couplings B1 B2/B3

(Rh-O-Rh)=93.4, (Rh-Rh)=2.95 (2), (Rh-O)=2.03 Å.
(Rh-O-Rh)=94.1, (Rh-Rh)=2.95 (4), (Rh-O)=2.03, (Rh-O)=2.00 Å. O and O are the two bridging O’s.

Table 2: Relative energies of the four low-lying magnetic states and the associated effective exchange couplings (meV) for two NN RhO octahedra in LiRhO. Two distinct types of such [RhO] units, B1 and B2/B3 (see text), are found experimentally [33]. Results of spin-orbit MRCI calculations are shown, with a local coordinate frame for each Rh-Rh link ( along the Rh-Rh bond, perpendicular to the RhO plaquette). The form of the actual lattice spin model is detailed in the SM.
Figure 1: Layered network of edge-sharing RhO octahedra in LiRhO. The two distinct types, B1 and B2/B3, of NN two-octahedra units and the honeycomb lattice of Rh sites are evidenced.
Figure 2: Dependence of the NN and on the Rh-O-Rh bond angle for an idealized structural model where all Rh-O bond lengths are set to the average value in the experimental crystal structure [33]. MRCI+SOC results are shown. The variation of the Rh-O-Rh angles is the result of gradual trigonal compression of the O octahedra.
Figure 3: (color online) (a) Phase diagram for the effective model of Eq. (1) supplemented by 2nd- and 3rd-neighbor isotropic couplings , . The NN effective interaction constants are set to the QC values provided in Table II. The spin structure for the diagonal zigzag state is also shown. (b) Sketch of the cluster used in the ED calculations; the site index runs from 1 to 24. Spin structure factors for with either (c) no spin defects or two spin defects at (d) and (e) . Spin configurations for the simplified Ising model are shown for spin-defect concentrations of (f) , (g) , (h) , (i) . Filled dark and light circles indicate opposite spin directions. Open circles show the position of spin defects.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description