XML Generation Notice This is an automatically generated XML version of the publication: "NaRuO2: Kitaev-Heisenberg exchange in triangular-lattice setting - https://doi.org/10.34657/13640", originally published under the license CC BY 4.0 International. The XML version was generated on 28-08-2026 by Technische Informationsbibliothek (TIB). In spite of high standards in the development of the script, we cannot outrule the possibility of mistakes in the process of file conversion. The authors are not responsible for possibly existing mistakes or differences in this XML-version. The link above directs you to the authorised version of the publication. If you spot mistakes in the XML-version, please contact us at oarepository@tib.eu. This XML-Version of the publication "NaRuO2: Kitaev-Heisenberg exchange in triangular-lattice setting - https://doi.org/10.34657/13640" is provided "as is", without any guarantee (“ausdrücklich oder stillschweigend”). TIB does not take responsibility for potential mistakes, missing passages and consequences resulting from this use. Hinweis zur XML-Generierung Dies ist eine durch die Technische Informationsbibliothek – TIB am 28-08-2026 automatisiert erstellte XML-Version der Publikation "NaRuO2: Kitaev-Heisenberg exchange in triangular-lattice setting - https://doi.org/10.34657/13640", die unter der Lizenz CC BY 4.0 International veröffentlicht wurde. Trotz sorgfältiger Entwicklung des Skripts können Fehler im Konvertierungsprozess nicht ausgeschlossen werden. Die Autor*innen der Originalpublikation sind für etwaige Fehler/Differenzen in der XML-Version nicht verantwortlich. Die oben verlinkte Datei stellt die autoritative Version der Publikation dar. Hinweise auf Fehler on der XML-Version nehmen wir unter oarepository@tib.eu gern entgegen. Die vorliegende XML-Version der Publikation "https://doi.org/10.1038/s41535-023-00582-7" wird "wie es ist" bereitgestellt, ohne jegliche Gewährleistung oder Garantie, sei es ausdrücklich oder stillschweigend. Die TIB übernimmt keine Verantwortung für etwaige Fehler, Auslassungen oder Folgen, die aus der Nutzung dieser Datei resultieren. xml_created 28-08-2026 dc.identifier.uri https://doi.org/10.1038/s41535-023-00582-7 dc.rights.license CC BY 4.0 International dc.relation.doi https://doi.org/10.34657/13640 NaRuO2: Kitaev-Heisenberg exchange in triangular-lattice setting NaRuO 2 : Kitaev-Heisenberg exchange in triangular-lattice setting, Bhattacharyya, PritamBogdanov, NikolayNishimoto, SatoshiWilson, StephenHozoi, Liviu, , (). Kitaev exchange, a new paradigm in quantum magnetism research, occurs for 90° metal-ligand-metal links, t2g5 transition ions, and sizable spin-orbit coupling. It is being studied in honeycomb compounds but also on triangular lattices. While for the former it is known by now that the Kitaev intersite couplings are ferromagnetic, for the latter the situation is unclear. Here we pin down the exchange mechanisms and determine the effective coupling constants in the t2g5 triangular-lattice material NaRuO2, recently found to host a quantum spin liquid ground state. We show that, compared to honeycomb compounds, the characteristic triangular-lattice cation surroundings dramatically affect exchange paths and effective coupling parameters, changing the Kitaev interactions to antiferromagnetic. Quantum chemical analysis combined with subsequent effective spin model simulations provide perspective onto the nature of the experimentally observed quantum spin liquid—it seemingly implies fairly large antiferromagnetic second-neighbor isotropic exchange, and the atypical proximity to ferromagnetic order is related to ferromagnetic nearest-neighbor Heisenberg coupling. INTRODUCTIONExchange through 90°metal-ligand-metal bonds represents one of the limiting cases in (super)exchange theory 1,2 . In the simplest situation of half-filled d-metal orbitals, this geometry is associated with Heisenberg ferromagnetism. Away from half-filling and for degenerate orbitals, however, very intricate physics may arise, as pointed out by Jackeli and Khaliullin for t 5 2g magnetic centers with sizable spin-orbit coupling: highly anisotropic effective interactions involving only (pseudo)spin components normal to the M 2 L 2 square plaquette formed by two transition-metal (TM) ions and the two bridging ligands 3,4 . In crystals of NaCl type and derivative rhombohedral structures which imply three different possible orientations of the M 2 L 2 plaquettes (see Fig. 1), this translates into directional dependence of the nearest-neighbor exchange: on differently oriented adjacent M 2 L 2 units-i.e., normal to either x, y, or z-different components of the magnetic moments interact, either Sx i and Sx j (red arrows in Fig. 1b for x-type, normal to the xaxis plaquettes), Sy i and Sy k (blue arrows, y-type plaquettes), or Sz i and Sz l (green arrows, z type). The (super)exchange model of Jackeli and Khaliullin 3 thus formalizes Kitaev's effective Hamiltonian of bond-dependent anisotropic magnetic couplings 5 proposed initially more like a heuristic device. It launched a whole new subfield in the research area of quantum magnetic materials, that of 5d 5 and 4d 5 honeycomb magnets 4 , with subsequent extension to 3d 7 t 5 2g e 2 g hexagonal networks of edge-sharing ML 6 octahedra.Here we explore the nature of nearest-neighbor couplings in a 4d 5 triangular-lattice magnet, NaRuO 2 , and evidence the presence of sizable bond-dependent Kitaev interactions. Interestingly, those are antiferromagnetic, different from the case of honeycomb 4d 5 and 5d 5 magnetic lattices. Preliminary numerical tests highlight the important role of cation species around the ligands mediating (super)exchange: the sign change is related to electrostatics having to do with adjacent TM ions that in the honeycomb setting are missing. Sizable antiferromagnetic off-diagonal intersite couplings are also predicted, along with a somewhat larger ferromagnetic Heisenberg exchange. The latter seems to be consistent qualitatively with features seen in experiment: incipient ferromagnetism within a quantum-disordered ground state 6,7 . RESULTSOn-site multiplet structure and intersite couplings in NaRuO 2 NaRuO 2 is rather unique, a realization of j eff ≈ 1/2 t 5 2g moments on a triangular network of edge-sharing RuO 6 octahedra (see Fig. 2). The basic Ru 3+ 4d 5 multiplet structure in the specific delafossite crystalline setting of NaRuO 2 is illustrated in Table 1, as obtained by quantum chemical 8 complete-active-space self-consistent-field (CASSCF) 8,9 and post-CASSCF multireference configurationinteraction (MRCI) 8,10 embedded-cluster calculations (see Methods for computational details and Supplementary Note 1 for basis sets). To separately visualize the effects of crystal-field (CF) splittings and spin-orbit coupling (SOC), both scalar relativistic (CASSCF, MRCI) and MRCI+SOC 11 results are listed. The t 2g -e g splitting is larger than in e.g., RuCl 3 12 , due to larger ligand effective charges and shorter TMligand bonds in NaRuO 2 . It is seen that the MRCI correlation treatment brings significant corrections to some of the CASSCF relative energies, in particular for the 6 A 1g term. Also important is the SOC (see the last two columns in Table 1). However, the effects of trigonal fields are visible too: those remove the degeneracy of the 2 T 2g CF states (p-orbital-like l eff = 1/2 13 and s = 1/2 quantum numbers in ideal cubic environment) and speak for significant deviations from textbook 3,13 j eff = 1/2 spin-orbit moments.Even with sizable trigonal fields, existing estimates of several meV for the Kitaev interactions in honeycomb-lattice 4d t 5 2g compounds such as RuCl 3 (see e.g. 12,[14][15][16][17][18][19] and discussion in ref. 20 ) and Li 2 RhO 3 21 provide strong motivation for detailed ab initio investigation of the effective couplings in NaRuO 2 . To this end, further quantum chemical computations were performed for a block of two adjacent edge-sharing RuO 6 octahedra, following the procedure described in refs. 12,22 (see also Methods). Mapping the quantum chemical data onto the relevant effective spin model for C 2h symmetry of the [Ru 2 O 10 ] magnetic unit,H ðγÞ ij ¼ J Si Á Sj þ K Sγ i Sγ j þ X α≠β Γ αβ ð Sα i Sβ j þ Sβ i Sα j Þ;(1)by spin-orbit MRCI, we arrive to Kitaev, Heisenberg, and off-diagonal Γ coupling constants K = 2.0, J = −5.2, Γ ≡ Γ xy = 3.6, Γ 0 Γ yz ¼ Γ zx ¼ 1:1 (meV). Interestingly, both K and Γ are antiferromagnetic, with Γ > K, but the largest interaction parameter is the Heisenberg J. The latter being ferromagnetic, it brings us away from the antiferromagnetic ground states (zigzag antiferromagnetic order, in particular) experimentally found in j ≈ 1/2 honeycomb systems, as depicted for instance in the phase diagrams computed for J-K-Γ triangular-lattice models in ref. 23 . According to ref. 23 , quantum spin liquid (QSL) ground states can only be realized for relatively large, antiferromagnetic Γ. Our result for the strength of Γ represents the largest ab initio quantum chemical estimate so far in a real material setting. Yet, a nearest-neighbor J that is even larger in magnitude and ferromagnetic (see the values listed above and at the top of Table 2) implies that only longer-range Heisenberg couplings (analyzed in isotropic context in e.g. refs. [24][25][26][27] ) and/or higherorder interactions such as ring exchange [28][29][30] can pull the system to the QSL regime evidenced experimentally 6 .Signatures of a nearby ferromagnetic state in NaRuO 2 appear via the enhanced Van Vleck term in the magnetic susceptibility 6 , which approaches values observed in nearly magnetic metals like Pd. Furthermore, despite possessing a charge gap and insulating ground state, clear quasiparticle excitations generate a Sommerfeld-like term in the low-temperature heat capacity 6 . This suggests the presence of strong spin fluctuations associated with a nearby magnetic state. These seemingly gapless spin fluctuations are robust, and are directly resolved in neutron scattering as diffuse continuum-like modes about the nuclear zone center and near Q = 0 6 . Initial attempts at scaling S(Q, E) suggest the proximity of a nearby ferromagnetic quantum critical point, and small, nonmagnetic chemical perturbations to the network of Ru 3+ ions induce glass-like freezing of the magnetic moments 7 . What makes K antiferromagneticFinding an antiferromagnetic Kitaev coupling in NaRuO 2 , opposite to the ferromagnetic K generally found in honeycomb Kitaev- Heisenberg magnets (see 12,21,31 for quantum chemical results, refs. [14][15][16][17] for estimates based on effective-model (super)exchange theory, and refs. 18,19 for values obtained by fitting experimental spectra) is intriguing. To shed light on this aspect, we performed the following numerical experiment: in a new set of quantum chemical computations, the two magnetic-plane cations in the immediate neighborhood of the bridging ligands, forming a line perpendicular to the link of nearest-neighbor Ru sites (see Fig. 2c), were removed and their charge redistributed within the pointcharge array modeling the extended solid-state surroundings.Remarkably, without those nearby positive ions, three of the effective magnetic couplings change sign (see 31 , produced in that case by cations residing between the magnetic layers. Ru 2 O 2 -unit correlations. Extended magnetic latticeinto the interplay between spin-orbit couplings and electron correlation effects on a plaquette of two nearestneighbor Ru ions and two bridging ligands can be obtained by analyzing computational data obtained at lower levels of approximation, i.e., spin-orbit calculations including only the t 5 2gt 5 2g electron configuration and CASSCF+SOC calculations that also account for excited-state t 4 2g -t 6 2g configurations since all possible occupations are considered within the six-orbital (three t 2g orbitals per Ru site) active space in the latter case.As illustrated in Table 3, anisotropic effective intersite interactions comparable in size with the isotropic Heisenberg J are already found at the single-configuration (SC) level, i.e., when accounting for just t 5 2g -t 5 2g Coulomb exchange (sometimes referred to as direct exchange 2 ). Large anisotropic Coulomb exchange as found in the SC calculation represents physics not addressed so far in the literature-Kitaev magnetism is presently exclusively explained through TM-TM kinetic exchange and TM-L-TM superexchange.CASSCF, through the inclusion of intersite t 2g → t 2g excitations (i.e., TM-TM kinetic exchange), brings sizable corrections to K, J, and Γ in particular. By considering next, in the spin-orbit MRCI treatment, all possible single and double excitations out of the Ru t 2g and bridging-ligand orbitals (Ru(t 2g )→Ru(e g ) and O(2p)→Ru(4d) transitions, among others), the most significant post-CASSCF corrections are found to occur for K and Γ. This suggests that more sophisticated calculations, e.g., MRCI+SOC based on larger active spaces in the prior CASSCF, would bring significant additional corrections to K and Γ mostly, i.e., would only enhance antiferromagnetic fluctuations on the extended lattice.An interesting observation is that, adding the Coulombexchange contributions computed at SC level for K and J (−1.0 and −0.7 meV, respectively, first line in Table 3) to the K and J estimates derived from effective-model calculations where Coulomb exchange is neglected (K = 2.9 and J = −4.2 meV) 32 brings the latter in rather good agreement with the MRCI values (2.0 and −5.2 meV, respectively), although this is not the case for Γ and Γ 0 .To obtain first impressions on the role of longer-range Heisenberg interactions, second-neighbor (J 2 ) and third-neighbor (J 3 ), in particular, we performed density-matrix renormalization group (DMRG) calculations 33,34 for a J-K-Γ-Γ 0 -J 2 -J 3 model on fragments of 19, 27, and 37 sites of the triangular lattice. In order to prevent artifacts, e.g., artificial stabilization (or destabilization) of particular magnetic states, we applied open boundary conditions. The validity of this material model is discussed in Supplementary Note 2 and Supplementary Note 3. Setting J, K, Γ, and Γ 0 to -5.5, 4, 4, and 1.5 meV, respectively, we obtain a very rich phase diagram, see Fig. 3. Notably, a QSL phase is found for antiferromagnetic J 2 values and ferromagnetic J 3 . The corresponding structure factor has no Bragg peaks but exhibits a characteristic pattern. It indicates that the QSL emerges from the competition of adjacent ordered states. The most remarkable spot is the region where the QSL neighbors ferromagnetic order and a spiral phase, namely, around J 2 = 5, J 3 = −2-as shown in Fig. 3, at J 2 = 5, J 3 = −2 the structure factor implies a broad spectral feature for momenta near |Q| = 0, seemingly consistent with experimental observations 6 . DISCUSSIONWhile the way nearby cations are structurally arranged can strongly affect on-site electronic-structure features such as subshell level splittings 22,35,36 , single-ion anisotropies 35 , and g factors 22,36 , detailed ab initio investigations of the effect on intersite magnetic couplings are less numerous. Here we show that, compared to honeycomb compounds, the characteristic triangular-lattice cation surroundings dramatically affect superexchange paths and effective coupling parameters. In particular, with respect to honeycomb j ≈ 1/2 systems, the Kitaev interaction constant changes its sign in triangular-lattice NaRuO 2 . By providing insight into the signs and strengths of the nearest-neighbor magnetic interactions in this material, our work sets the frame within which the role of longer-range effective spin couplings should be addressed. Interestingly, while giving rise to antiferromagnetic order on honeycomb Kitaev-Heisenberg lattices, the latter appear to be decisive in realizing quantum disorder 6 in NaRuO 2 . Last but not least, we establish the role of anisotropic Coulomb exchange (in trianguar-lattice setting), a mechanism not addressed so far in Kitaev-Heisenberg magnetism. METHODS Ru-site multiplet structureAll quantum chemical computations were carried out using the MOLPRO suite of programs 37 . Crystallographic data as reported by Ortiz et al. 6 were utilized. For each type of embedded cluster, the crystalline environment was modeled as a large array of point charges which reproduces the crystalline Madelung field within the cluster volume; we employed the EWALD program 38,39 to generate the point-charge embeddings.To clarify the Ru-site multiplet structure, a cluster consisting of one 'central' RuO 6 octahedron, the six in-plane adjacent octahedra, and 12 nearby Na cations was designed. The quantum chemical study was initiated as CASSCF calculations 8,9 with an active orbital space containing the five 4d orbitals of the central Ru ion. Post-CASSCF correlation computations were performed at the MRCI level, with single and double excitations 8,10 out of the Ru 4d and O 2p orbitals of the central RuO 6 octahedron. SOCs were accounted for following the procedure described in ref. 11 . Intersite exchange in NaRuO 2Clusters with two edge-sharing RuO 6 octahedra in the central region were considered in order to derive the intersite effective magnetic couplings. The eight in-plane octahedra directly linked to the two-octahedra central unit and 22 nearby Na ions were also explicitly included in the quantum chemical computations but using much more compact basis sets. CASSCF computations were carried out with six (Ru t 2g ) valence orbitals and ten electrons as active (abbreviated as (10e,6o) active space); the t 2g orbitals of the adjacent TM ions were part of the inactive orbital space. In the subsequent MRCI correlation treatment, single and double excitations out of the central-unit Ru t 2g and bridging-O 2p levels were considered. We used the Pipek-Mezey methodology 40 to obtain localized central-unit orbitals.The CASSCF optimization was performed for the lowest nine singlet and lowest nine triplet states associated with the (10e,6o) setting. Those were the states for which SOCs were further accounted for 11 , either at SC, CASSCF, or MRCI level, which finally yields a number of 36 spin-orbit states in each case. The SC label in Table 3 in the main text stands for a CASCI in which intersite excitations are not considered. This is also referred to as occupation-restricted multiple active space (ORMAS) scheme 41 .The lowest four spin-orbit eigenstates from the MOLPRO output (eigenvalues lower by ~200 meV or more than the eigenvalues of higher-lying excited states, as illustrated for example in Table 1) were mapped onto the eigenvectors of the effective spin Hamiltonian (1), following the procedure described in refs. 12,22 : those four expectation values and the matrix elements of the Zeeman Hamiltonian in the basis of the four lowest-energy spinorbit eigenvectors are put in direct correspondence with the respective eigenvalues and matrix elements of (1). Having two of the states in the same irreducible representation of the C 2h point group 21 , such one-to-one mapping translates into two possible sets of effective magnetic couplings. The relevant array is chosen as the one whose g factors fit the g factors obtained for a single RuO 6 t 5 2g octahedron. We used the standard coordinate frame usually employed in the literature, different from the rotated frame employed in earlier quantum chemical studies 12,21,31 that affects the sign of Γ. This is the reason the sign of Γ for RuCl 3 in Table 2 in the main text is different from the sign in ref. 12 (see also footnote [48] in ref. 20 ). Fig. 3 Single-layer magnetism. Ground-state phase diagram obtained by DMRG computations for the J-K-Γ-Γ 0 -J 2 -J 3 model (right) and typical static spin structure factors for each phase (left). Using the quantum chemical analysis as a guide (see data in Table 3 and related discussion), J, K, Γ, and Γ 0 were set to to −5.5, 4, 4, and 1.5 meV, respectively. Fig. 1 Fig. 1 Orthogonal M 2 L 2 plaquettes in solids. a Rocksalt lattice, the prototype for 90°M-L-M exchange 1,2 . b If two different cation varieties-alkaline (A, yellow) and magnetic TM (B) ions-form successive layers normal to the [111] axis, a rhombohedral ABL 2 structure results, with a triangular network of edge-sharing octahedra (either AL 6 or BL 6 ) in each layer. On each B 2 L 2 plaquette (atoms not explicitly shown), the Kitaev interaction couples only (pseudo)spin components perpendicular to the respective plaquette. Red, blue, and green indicate x, y, and z projections, respectively; only two spin components are shown at each magnetic site. Fig. 2 Fig. 2 NaRuO 2 crystal structure. a Sketch of successive layers of edge-sharing RuO 6 and NaO 6 octahedra. Gray, red, and yellow spheres represent Ru, O, and Na ions, respectively. b The triangular magnetic network of RuO 6 units. c To understand how the nearby surroundings affect magnetic couplings, Ru neighbors in the median plane (sites 3, 4) of a given Ru-Ru magnetic link (sites 1, 2) were removed in one set of quantum chemical computations. Table 1 . Table 1 . Ru 3+ 4d 5 multiplet structure, with all five 4d orbitals active in CASSCF.Each MRCI+SOC value refers to a Kramers doublet (KD); just the lowest and highest KDs are shown for each group of t 4 2g e 1 g spin-orbit states. Only the CF terms listed in the table entered the spin-orbit calculations. For simplicity, notations corresponding to O h symmetry are used. Ru 3+ 4d 5 CASSCF MRCI MRCI splittings (eV) +SOC 2 T 2g (t 5 2g ) 0 0 0 0.11 0.12 0.21 0.11 0.12 0.26 4 T 1g (t 4 2g e 1 g ) 1.42 1.58 1.61 1.45 1.62 | 1.45 1.62 1.74 6 A 1g (t 3 2g e 2 g ) 1.42 1.91 2.08 (×3) 4 T 2g (t 4 2g e 1 g ) 2.16 2.28 2.37 2.16 2.28 | 2.20 2.30 2.42 Table 2 . Table 2 . Nearest-neighbor effective magnetic couplings (MRCI+SOC) for experimentally determined crystal structures of NaRuO 2 /RuCl 3 and for adjusted configurations with modified median plane (MMP) cation distribution (see main text and Fig. 2c for details). MRCI+SOC (meV) K J Γ Γ 0 NaRuO 2 2.0 -5.2 3.6 1.1 NaRuO 2 , MMP -19 12 1.2 -1.1 RuCl 3 12 -5.6 1.2 1.2 -0.7 RuCl 3 , MMP 1.9 -12 1.1 -0.1 Table 2 ) Table 2 ) . This suggests that strong orbital polarization effects at the ligand sites, induced by the extra nearby positive charge on the triangular lattice (+3 effective nearby charges in NaRuO 2 versus +1 in honeycomb '213' compounds such as Na 2 IrO 3 and Li 2 RhO 3 or 0 in RuCl 3 ) have dramatic consequences on hopping matrix elements and superexchange processes. Similar additional tests for RuCl 3 , where +3 ions were placed in the centers of two edge-sharing Ru 6 hexagonal rings, confirm the trend: the extra positive charge in the neighborhood of the Ru 2 O 2 'magnetic' plaquette (i.e., [Ru 2 O 10 ] cluster of edge-sharing octahedra) invert the sign of the Kitaev coupling constant, from ferromagnetic in actual RuCl 3 12 to antiferromagnetic in the presence of positive charge in the centers of the two hexagons sharing the magnetic Ru-Ru link (see lowest lines in Table 2). Strong polarization effects of this type were earlier pointed out in the case of the Kitaev-Heisenberg honeycomb iridate H 3 LiIr 2 O 6 Table 3 . Table 3 . Nearest-neighbor magnetic couplings (meV) at different levels of theory: SC (only the t 5 2g -t 5 2g configuration considered), CASSCF (also t 4 2g -t 6 2g states included), and MRCI (single and double excitations out of the Ru 4d t 2g and bridging-O 2p levels on top of CASSCF). Magnetic couplings K J Γ Γ 0 SC+SOC -1.0 -0.7 0.4 0.8 CASSCF+SOC 1.3 -5.2 2.2 1.0 MRCI+SOC 2.0 -5.2 3.6 1.1 ACKNOWLEDGEMENTS P. B. and L. H. acknowledge financial support from the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG), Project ID 441216021, and technical assistance from U. Nitzsche. S. N. acknowledges financial support via project A05 of the Collaborative Research Center 1143 of the DFG (Project ID 247310070). S. D. W. acknowledges support via DOE, Office of Science, Basic Energy Sciences under Award DE-SC0017752. P. B. thanks M. S. Eldeeb for discussions. We also acknowledge instructive discussions with I. Rousochatzakis, U. K. Rößler, and G. Khaliullin. DATA AVAILABILITYQuantum chemical data (atomic coordinates, point-charge data, MOLPRO output) have been deposited in the RADAR database under the https://doi.org/10.22000/1736. CODE AVAILABILITYRelevant codes are available from the corresponding author upon reasonable request. AUTHOR CONTRIBUTIONS FUNDINGOpen Access funding enabled and organized by Projekt DEAL. COMPETING INTERESTSThe authors declare no competing interests. ADDITIONAL INFORMATION Supplementary informationThe online version contains supplementary material available at https://doi.org/10.1038/s41535-023-00582-7.Correspondence and requests for materials should be addressed to Pritam Bhattacharyya or Liviu Hozoi.Reprints and permission information is available at http://www.nature.com/ reprints Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Superexchange interaction and symmetry properties of electron orbitals, Kanamori, JJ. Phys. Chem. Solids, 10, 87 (1959). New approach to the theory of superexchange interactions, Anderson, PPhys. Rev, 115, (1959). Mott insulators in the strong spin-orbit coupling limit: from heisenberg to a quantum compass and kitaev models, Jackeli, GKhaliullin, GPhys. Rev. Lett, 102, 17205 (2009). Concept and realization of kitaev quantum spin liquids, Takagi, HTakayama, TJackeli, GKhaliullin, GNagler, SNat. Rev. Phys, 1, 264 (2019). Anyons in an exactly solved model and beyond, Kitaev, AAnn. Phys, 321, (2006). Quantum disordered ground state in the triangular-lattice magnet NaRuO 2, Ortiz, BNat. Phys, 19, 943 (2023). Defect control in the heisenberg-kitaev candidate material NaRuO 2, Ortiz, BSarte, PAvidor, AWilson, SPhys. Rev. Mater, 6, 104413 (2022). , Molecular electronic structure theory, , (2000). MCSCF optimization revisited. II. combined first-and second-order orbital optimization for large molecules, Kreplin, DKnowles, PWerner, H.-JJ. Chem. Phys, 152, 74102 (2020). Internally contracted multiconfiguration-reference configuration interaction calculations for excited states, Knowles, PWerner, H.-JTheor. Chim. Acta, 84, 95 (1992). Spin-orbit matrix elements for internally contracted multireference configuration interaction wavefunctions, Berning, ASchweizer, MWerner, H.-JKnowles, PPalmieri, PMol. Phys, 98, 1823 (2000). Kitaev exchange and field-induced quantum spin-liquid states in honeycomb α-RuCl 3, Yadav, RSci. Rep, 6, 37925 (2016). , Electron Paramagnetic Resonance of Transition Ions, , (1970). Crystal structure and magnetism in α-RuCl 3 : An ab initio study, Kim, H.-SKee, H.-YPhys. Rev. B, 93, 155143 (2016). Challenges in design of kitaev materials: Magnetic interactions from competing energy scales, Winter, SLi, YJeschke, HValenti, RPhys. Rev. B, 93, 214431 (2016). Theoretical investigation of magnetic dynamics in α-RuCl 3, Wang, WDong, Z.-YYu, S.-LLi, J.-XPhys. Rev. B, 96, 115103 (2017). Deriving models for the kitaev spin-liquid candidate material α-RuCl 3 from first principles, Eichstaedt, CPhys. Rev. B, 100, 75110 (2019). Identification of magnetic interactions and high-field quantum spin liquid in α-RuCl 3, Li, HNat. Commun, 12, 4007 (2021). Extraction of interaction parameters for α-RuCl 3 from neutron data using machine learning, Samarakoon, APhys. Rev. Res, 4, L022061 (2022). Magnetization processes of zigzag states on the honeycomb lattice: Identifying spin models for α-RuCl and Na 2 IrO 3, Janssen, LAndrade, EVojta, MPhys. Rev. B, 96, 64430 (2017). Strong magnetic frustration and anti-site disorder causing spin-glass behavior in honeycomb Li 2 RhO 3, Katukuri, VSci. Rep, 5, 14718 (2015). Orbital reconstruction in nonpolar tetravalent transitionmetal oxide layers, Bogdanov, NNat. Commun, 6, 7306 (2015). Comprehensive study of the global phase diagram of the J − K − Γ model on a triangular lattice, Wang, SPhys. Rev. B, 103, 54410 (2021). Theory of a competitive spin liquid state for weak mott insulators on the triangular lattice, Mishmash, RGarrison, JBieri, SXu, CPhys. Rev. Lett, 111, 157203 (2013). Spin liquid phase of the S ¼ 1 2 J 1 À J 2 Heisenberg model on the triangular lattice, Zhu, ZWhite, SPhys. Rev. B, 92, 41105 (2015). Spin liquid nature in the heisenberg J 1 −J 2 triangular antiferromagnet, Iqbal, YHu, W.-JThomale, RPoilblanc, DBecca, FPhys. Rev. B, 93, 144411 (2016). Dirac spin liquid on the spin-1/2 triangular heisenberg antiferromagnet, Hu, SZhu, WEggert, SHe, Y.-CPhys. Rev. Lett, 123, 207203 (2019). Variational study of triangular lattice spin-1/2 model with ring exchanges and spin liquid state in k-(ET) 2 Cu 2 (CN) 3, Motrunich, OPhys. Rev. B, 72, 45105 (2005). Four-spin terms and the origin of the chiral spin liquid in mott insulators on the triangular lattice, Cookmeyer, TMotruk, JMoore, JPhys. Rev. Lett, 127, 87201 (2021). Thermal properties and instability of a U(1) spin liquid on the triangular lattice, Zhao, Q.-RLiu, Z.-XPhys. Rev. Lett, 127, 127205 (2021). Engineering kitaev exchange in stacked iridate layers: impact of inter-layer species on in-plane magnetism, Yadav, RChem. Sci, 10, 1866 (2019). , A j eff = 1/2 kitaev material on the triangular lattice: the case of NaRuO 2 . npj Quantum Mater, 8, 36 (2023). Density matrix formulation for quantum renormalization groups, White, SPhys. Rev. Lett, 69, 2863 (1992). , The ITensor Software Library for Tensor Network Calculations, , 4 (2022). , Bogdanov, NMaurice, RRousochatzakis, IVan Den Brink, JHozoi, LPhys. Rev. Lett, 110, 127206 (2013). Single-site magnetic anisotropy governed by interlayer cation charge imbalance in triangular-lattice AYbX 2, Zangeneh, ZAvdoshenko, SVan Den Brink, JHozoi, LPhys. Rev. B, 100, 174436 (2019). Molpro: a generalpurpose quantum chemistry program package, Werner, H.-JKnowles, PKnizia, GManby, FSchütz, MWIREs Comput. Mol. Sci, 2, 242 (2012). Accurate crystal fields for embedded cluster calculations, Klintenberg, MDerenzo, SWeber, MComp. Phys. Commun, 131, 120 (2000). Determining point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations, Derenzo, SKlintenberg, MWeber, MJ. Chem. Phys, 112, 2074 (2000). A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions, Pipek, JMezey, PJ. Chem. Phys, 90, 4916 (1989). Direct configuration interaction and multiconfigurational self-consistentfield method for multiple active spaces with variable occupations. I. Method, Ivanic, JJ. Chem. Phys, 119, 9364 (2003).