⇦Chinese Physics Letters, 2021, Vol. 38, No. 4, Article code 047502Express Letter Rare-Earth Chalcohalides: A Family of van der Waals Layered Kitaev Spin Liquid Candidates Jianting Ji (籍建葶)1†, Mengjie Sun (孙梦杰)1,2†, Yanzhen Cai (蔡焱桢)3, Yimeng Wang (王艺朦)1,2, Yingqi Sun (孙颖琦)1,2, Wei Ren (任玮)3, Zheng Zhang (张政)1,2, Feng Jin (金峰)1, and Qingming Zhang (张清明)3,1* Affiliations 1Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 2Department of Physics, Renmin University of China, Beijing 100872, China 3School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China Received 9 March 2021; accepted 27 March 2021; published online 2 April 2021 Supported by the National Key Research and Development Program of China (Grant Nos. 2017YFA0302904 and 2016YFA0300504), the National Natural Science Founation of China (Grant Nos. U1932215 and 11774419), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB33010100).
†These authors contributed equally to this work.
*Corresponding author. Email: qmzhang@ruc.edu.cn Citation Text: Ji J T, Sun M J, Cai Y Z, Wang Y M, and Sun Y Q et al. 2021 Chin. Phys. Lett. 38 047502    Abstract The Kitaev spin liquid (KSL) system has attracted tremendous attention in recent years because of its fundamental significance in condensed matter physics and promising applications in fault-tolerant topological quantum computation. Material realization of such a system remains a major challenge in the field due to the unusual configuration of anisotropic spin interactions, though great effort has been made before. Here we reveal that rare-earth chalcohalides REChX (RE = rare earth; Ch = O, S, Se, Te; X = F, Cl, Br, I) can serve as a family of KSL candidates. Most family members have the typical SmSI-type structure with a high symmetry of $R\bar{3}m$, and rare-earth magnetic ions form an undistorted honeycomb lattice. The strong spin-orbit coupling of $4f$ electrons intrinsically offers anisotropic spin interactions as required by the Kitaev model. We have grown the crystals of YbOCl and synthesized the polycrystals of SmSI, ErOF, HoOF and DyOF, and made careful structural characterizations. We carry out magnetic and heat capacity measurements down to 1.8 K and find no obvious magnetic transition in all the samples but DyOF. The van der Waals interlayer coupling highlights the true two-dimensionality of the family which is vital for the exact realization of Abelian/non-Abelian anyons, and the graphene-like feature will be a prominent advantage for developing miniaturized devices. The family is expected to act as an inspiring material platform for the exploration of KSL physics. DOI:10.1088/0256-307X/38/4/047502 © 2021 Chinese Physics Society Article Text Quantum spin liquid (QSL) refers to a novel quantum spin disordered state with strong spin entanglement, which hosts many exotic excitations and phenomena.[1] Generally it can be driven by geometrical spin frustrations plus quantum spin fluctuations. Kitaev proposed an exactly solvable spin model on honeycomb lattice with the frustration of exchange coupling rather than a geometrical one,[2] which offers an alternative way to realize QSL. The most intriguing feature of the model is that it has a gapless QSL phase which allows excitations of non-Abelian anyons when a spin gap is stabilized by external magnetic fields. This opens the door to the scheme of topological quantum computation based on Kitaev spin liquid. On the other hand, the bond-dependent anisotropic exchange coupling is really unusual and challenges general experience in material realization of quantum spin systems. Strong spin-orbit coupling (SOC) is believed to be a feasible mechanism yielding such a kind of anisotropic spin interactions.[3] Jackeli and Khaliulin proposed an attractive scenario based on SOC,[4] in which both spin-orbit-coupled effective moments and “spin” interactions between them are highly anisotropic. The isotropic term can be canceled out through a delicate arrangement of mediating anions and the bond-dependent term becomes dominant. This greatly stimulates and clearly guides the exploration of $4d/5d$ Kitaev materials. The $5d$ iridates A$_2$IrO$_3$ (A = Li, Na, Cu) are the first Kitaev materials predicted by the Jackeli–Khaliulin mechanism.[4–7] The Ir$^{4+}$ (5$d^5$) ions are centered in the edge-sharing IrO$_6$ octahedra and form a honeycomb network normal to the diagonals. A crystalline electric field (CEF) ground state of Kramers doublet with an effective spin $J_{\rm eff}$=1/2 is created in the ions by the strong SOC and the CEF effect. The iridates are magnetically ordered at low temperatures, indicating a non-KSL ground state. The related compound A$_3$LiIr$_2$O$_6$ (A = Ag, Cu, H), which is structurally similar to the original iridates, was also reported.[8–10] Particularly in H$_3$LiIr$_2$O$_6$, no magnetic transition was observed down to 50 mK. Whether or not this is an evidence for KSL ground state still remains to be clarified, because of the substantial existence of quenched disorders in the compound.[10] The $4d^5$ Li$_2$RhO$_3$ with the similar structure and $J_{\rm eff}=1/2$ was also investigated and a partial spin-freezing was observed at 6 K.[11] Recently, the $4d$ ruthenium halide $\alpha$-RuCl$_3$ was extensively studied as a KSL candidate.[12–19] Similar to the Ir$^{4+}$ ions in the iridates, Ru$^{3+}$($4d^5$) ions have an effective spin $J_{\rm eff}=1/2$ from Kramers doublet guaranteed by the strong SOC and the CEF effect, but sit on a slightly distorted honeycomb lattice at room temperature. It undergoes an antiferromagnetic (AFM) transition around 8 K, also indicating a non-KSL ground state. Neutron scattering experiments reported the unusual signals possibly from Majorana excitations, suggesting that $\alpha$-RuCl$_3$ is proximate to a KSL system.[18] With the suppression of the magnetic transition by in-plane fields, Kasahara et al. reported the observation of a half-integer quantum Hall effect, which was considered as an indication of Majorana excitations in thermodynamic channel.[19] Along the way of trichlorides, more KSL candidates have been explored very recently. YbCl$_3$, having a monoclinic structure and an effective spin $J_{\rm eff}$=1/2, was suggested to be proximate to the KSL physics.[20] Also, the short- and long-range magnetic transitions were observed at low temperatures in the compound. The $5d$ osmium trichloride Os$_x$Cl$_3$ was reported to be a KSL candidate very recently.[21] Just like the cases mentioned above, Os$^{3+}$ ($5d^5$) ions have an effective spin $J_{\rm eff}=1/2$, but originally sit on a triangular rather than honeycomb lattice. It was suggested that there exist nano-domains in which Os$^{3+}$ ions are placed on a honeycomb network. The deficiency of osmium remains a crucial issue in the material. As we can see, the current effort for searching for KSL candidates mainly focuses on $4d/5d$ systems, as inspired by the original Jackeli–Khaliulin scenario. Actually, most rare-earth $4f$ ions typically have strong SOC and an effective spin 1/2 guaranteed by Kramers doublet, and hence are natural building blocks for KSL candidates. Theoretical considerations in this aspect have been made,[22–24] while rare-earth-based KSL materials are little explored at present. In our previous studies, we have revealed the geometrically frustrated rare-earth-based quantum magnet YbMgGaO$_4$ and the family of rare-earth chalcogenides,[25–27] which serve as an ideal platform for the exploration of QSL physics. In the rare-earth materials, exchange coupling is relatively small (normally several kelvins) because of the screening of the outer $5s/5p$ shells. However, this also brings us some advantages, for example, one can safely consider the nearest-neighbor spin interactions without worrying about further neighbors and the low-energy spin Hamiltonian can be well investigated by the application of magnetic fields generated in the lab. In this study, we point out that rare-earth chalcohalides REChX (RE = rare earth; Ch = O, S, Se, Te; X = F, Cl, Br, I) may be a promising family of KSL materials. The general crystal structure of the family belongs to the SmSI-type or its polytypes which basically show a high symmetry of $R\bar{3}m$. This allows rare-earth magnetic ions to form an undistorted honeycomb lattice. The honeycomb magnetic layers are well separated and there is only a weak van der Waals coupling between the layers. We have grown the crystals of YbOCl and synthesized the polycrystals of SmSI, ErOF, HoOF and DyOF. They are structurally determined by x-ray diffraction measurements and careful Rietveld refinements. We carry out magnetic susceptibility and heat capacity measurements down to 1.8 K, showing no obvious phase transitions in the samples. We further discuss anisotropic spin interactions in the compounds. Synthesis of Crystals/Polycrystals and Experiments. The crystals of YbOCl have been grown by a standard flux method.[28] The operations related to raw materials are carried out in a glove box. Yb$_2$O$_3$ (99.9$\%$, Alfa) and anhydrous YbCl$_3$ (99.9$\%$, Alfa) are mixed with the ratio of $1\!:\!10$, and the mixed raw materials are sealed in a quartz tube filled with argon gas. The sealed tubes are placed in a box furnace and heat up to 1050 ℃ at which the furnace stay for four days. Then the furnace is cooled down to room temperature with the rate of 8℃ per hour. The products are filtered and cleaned in distilled water. Finally the nearly colorless YbOCl crystals are obtained with the maximum size of $\sim$$5.0 \times 5.0 \times 0.1$ mm$^3$ (see the inset in Fig. 1). The SmSI polycrystals are synthesized with a standard solid state reaction method.[29,30] Samarium and sulfur powder together with iodine balls are mixed up with the ratio of $1\!:\!1$:1 in a glove box and sealed in a quartz tube. The mixtures are placed in a box furnace and heated up to 450–500℃ at which the furnace stays for some time. Then temperature is slowly increased to 830–900℃ at which the furnace stays for one week. The polycrystals are obtained after the furnace is naturally cooled down to room temperature. It should be noted here that SmSI is quite unstable in air. The sample is sealed with sticky tape when making structural characterizations. The polycrystals of ErOF, DyOF and HoOF are also synthesized by a solid state reaction method.[31] The raw materials are Ln$_2$O$_3$ powder and LnF$_3$ powder (99.99$\%$, Aladdin, Ln = Er, Dy, Ho), which are mixed up in a glove box with the molar ratios of Er$_2$O$_3$:ErF$_3=0.998$:1, Dy$_2$O$_3$:DyF$_3=1$:1, and Ho$_2$O$_3$:HoF$_3=1$:1, respectively. The mixtures are sealed in a quartz tube and slowly heated up to 800℃ in a box furnace. After staying at 800℃ for 60 min, the furnace is slowly cooled down to 600℃ and then naturally cooled down to room temperature. The YbOCl single crystal is structurally analyzed by a single crystal x-ray diffractometer (XRD) (Bruker D8 Venture) and the structure was solved and refined using the Bruker SHELXTL software package. The polycrystals are measured using a powder XRD (Bruker D8 Da Vinci) with the mode of step scanning. The crystal structures are refined using the TOPAS program. Magnetic susceptibility, magnetization and heat capacity measurements from 1.8 K to 300 K are performed with a vibrating sample magnetometer (Quantum Design, physical property measurement system, PPMS). For magnetic measurements, the polycrystalline samples which are pressed into thin plates, or YbOCl crystals are fixed on the sample holder with GE varnish. For heat capacity measurements, the samples are mounted with N grease for better thermal contact. The background heat capacity from the sample holder and grease is measured before sample is mounted, and subtracted in the analysis of heat capacity from samples.
cpl-38-4-047502-fig1.png
Fig. 1. Crystal structure of YbOCl/SmSI and their structural refinements. The upper middle shows the polyhedron surrounding Yb/Sm and the lower middle is the top view along the $c$ axis showing the honeycomb lattice formed by Yb/Sm ions. The elevated XRD background at small angles in SmSI comes from the sticky tape which was used to protect the air-sensitive sample.
cpl-38-4-047502-fig2.png
Fig. 2. Crystal structure of ErOF/HoOF/DyOF and their structural refinements. The upper middle shows the polyhedron surrounding Er/Ho/Dy and the lower middle is the top view along the $c$ axis showing the honeycomb lattice formed by Er/Ho/Dy ions.
Results and Discussions. We have determined the crystal structure of YbOCl and SmSI (Fig. 1). Both belong to the SmSI-type structure with a high symmetry of $R\bar{3}m$. A central rare-earth ion is coordinated by three halide anions and four chalcogen anions (see upper middle of Fig. 1). It is quite rare since such a rare-earth polyhedron contains two kinds rather than one kind of anions. Though the YbO$_4$Cl$_3$/SmS$_4$I$_3$ polyhedron looks a little more complicated, interestingly it retains the symmetry of $C_{3v}$. The configuration of CEF splitting is, in principle, determined by the point-group symmetries of crystals. This suggests that the CEF splitting in YbOCl and SmSI is formally similar to the cases in YbMgGaO$_4$ and rare-earth chalcogenides.[25–27,32] In other words, one can still expect the configurations of four Kramers doublets for Yb$^{3+}$, eight Kramers doublets for Er$^{3+}$, three Kramers doublets for Sm$^{3+}$, etc. The basic structural feature is that two adjacent magnetic layers form a double-layer sheet which, as a whole, are well separated from the neighboring double-layer ones with a weak van der Waals inter-sheet coupling. The A/B sublattices of a bipartite honeycomb lattice are precisely assigned to two layers in a double-layer sheet. In other words, the spin sites in a double-layer sheet, though not exactly in a flat plane, form a regular honeycomb lattice without any bond distortion (see lower middle of Fig. 1). The spin configuration well matches the original Kitaev model. Furthermore, the nearest-neighbor Dzyaloshinskii–Moriya (DM) interaction is symmetrically prohibited and the next nearest-neighbor one is allowed. This is expected to bring more interesting spin physics. We will have more discussions on it later.
cpl-38-4-047502-fig3.png
Fig. 3. Different stacking order in the compounds studied here. Adjacent magnetic layers form charge-neutral double-layer sheets which are well separated with a weak van der Waals interaction. YbOCl/SmSI has an ABC stacking order (left panel) while an ACB order for ErOF/HoOF/DyOF (right panel).
cpl-38-4-047502-fig4.png
Fig. 4. Magnetic susceptibility for (a)$\&$(b) YbOCl/SmSI and (c)$\&$(d) ErOF/HoOF/DyOF. The temperature ranges for C-W fitting [(b)$\&$(d)] are estimated and listed in Table 2 and the factors for the estimations are discussed in the context. In DyOF, there exists a drop around 3 K, which suggests an AFM-like transition and is consistent with heat capacity results (Fig. 6). The other compounds show no obvious transition down to 1.8 K.
The structure of Er/Dy/HoOF is also determined using powder x-ray diffraction and Rietveld refinements (Fig. 2). Basically, it is isostructural to the SmSI-structure and the high symmetry of $R\bar{3}m$ is still kept. The slight difference is the stacking order of the double-layer sheets, which is an A–B–C stacking in YbOCl/SmSI and becomes an A–C–B one in Er/Ho/DyOF (Fig. 3). Compared to YbOCl/SmSI, Er/Ho/DyOF serial compounds exhibit larger ratios of intra-sheet to inter-sheet distances, which offers a potential way to manipulate the inter-layer coupling in both charge and spin channels. Some key distances have been listed in Table 1. We would like to highlight that the weak van der Waals coupling guarantees the true two-dimensionality, which is a dispensable prerequisite for the exact realization of Abelian/non-Abelian anyons. In principle, even a small deviation from the two-dimensionality family will invalidate the braiding statistics for anyons and hence is potentially disadvantageous to the scheme of topological computation.
Table 1. Some key distances in the refined crystal structures. Here $a$ and $c$ are lattice constants and the other ones are illustrated in Figs. 1 and 2.
Compounds Symmetry $a$ (nm) $c$ (nm) $d$ (nm) $D$ (nm) $d_1$ (nm) $d_2$ (nm)
YbOCl $R\bar{3}m$ 0.37199 2.7791 0.28203 0.6443 0.3545 0.372
SmSI $R\bar{3}m$ 0.45575 3.26631 0.33464 0.75413 0.42570 0.45575
ErOF $R\bar{3}m$ 0.38741 1.88293 0.28068 0.34696 0.35569 0.37841
HoOF $R\bar{3}m$ 0.38015 1.89104 0.28189 0.34846 0.35726 0.38016
Table 2. Magnetic parameters extracted from the Curie–Weiss analysis $\chi=C/(T-\theta_{\rm cw}$). Here $\mu_{\rm eff}$ refers to the effective moments in units of Bohr magnetons. The temperature ranges for C-W fitting are listed.
Compounds Direction $C$ $\theta_{\rm cw}$(K) $\mu_{\rm eff}$/$\mu_{_{\rm B}}$ Temperature range (K)
YbOCl $c$-axis 0.21785 $-1.22$ 1.32 25–55
YbOCl $ab$-plane 0.82372 $-6.15$ 2.567 25–55
SmSI 0.13716 $-18.94$ 1.0475 10–25
ErOF 11.72 $-13.0$ 9.683 25–55
HoOF 13.265 $-5.74$ 10.3 25–55
Before we go to magnetic properties of the compounds, it should be pointed out that the rare-earth ions with an odd number of $4f$ electrons generally exhibit an effective spin $J_{\rm eff}=1/2$ related to Kramers doublets guaranteed by time reversal symmetry, while the ones with an even number of electrons having no Kramers doublets need to be discussed circumstantially. In our case, YbOCl, SmSI and ErOF with $J_{\rm eff}=1/2$ should be more relevant to KSL, while Ho$^{3+}$ and Dy$^{3+}$ with larger moments may involve dipolar interactions and behave more classically. Figure 4 shows magnetic susceptibilities down to 1.8 K, which fall into two categories [Figs. 4(a) and 4(c)] according to the structural similarity. The Curie–Weiss (C-W) analysis has been made based on the susceptibilities [Figs. 4(b) and 4(d)]. It should be noted that a reasonable temperature range for the C-W analysis strongly depends on CEF excitations of rare-earth magnetic ions. More specifically, it requires that the upper limit is far below the first CEF excitation and the lower one is far above exchange coupling. Though the CEF excitations in the compounds are not investigated yet, generally a compound with better ionicity is expected to have higher CEF excitations. Therefore, one can presume that the compounds studied here have larger first CEF excitation levels comparable to those of other rare-earth oxides like YbMgGaO$_4$.[33] Considering the exchange coupling of several kelvins for rare-earth magnetic ions due to the screening, we are able to estimate the effective C-W temperature ranges for the compounds. The C-W fitting details have been summarized in Table 2. It should be pointed out that the reverse susceptibility of SmSI in Fig. 4(b) shows a less linear behavior which outputs a relatively high C-W temperature of $\sim $18 K. At present it is unclear if the deviation from a linear form is related to the intrinsic spin and CEF excitations in SmSI, or extrinsically caused by the air-sensitivity. More careful investigations are required to clarify the issue in the future.
cpl-38-4-047502-fig5.png
Fig. 5. Magnetization for (a) YbOCl/SmSI and (b) ErOF/HoOF/DyOF. SmSI reaches no magnetic saturation up to 7 T while the other four are saturated around 1–3 T.
Magnetic susceptibilities demonstrate that all the compounds but DyOF show no obvious phase transitions at least down to 1.8 K. This is further evidenced by magnetization (Fig. 5) and heat capacity (Fig. 6) measurements. DyOF undergoes an AFM-like transition around 3 K. It remains to be uncovered that the transition is driven by spin-spin interactions or by dipolar interactions. YbOCl shows a magnetic anisotropy of $\sim $2:1 (Table 2) which is not very high and may favor the $K$ term rather than anisotropic off-diagonal terms.[34] There is an unusual upturn in heat capacity of YbOCl with lowering temperatures down to 1.8 K [Fig. 6(a)]. Further measurements down to milikelvins are highly required to figure out whether or not the upturn is an indication of a phase transition. ErOF and HoOF show broad features in heat capacity at temperatures above 10 K [Fig. 6(b)], which are higher than the C-W temperatures. The features are possibly related to a large number of CEF excitation levels in the two rare-earth ions and detailed discussions can be found somewhere. SmSI has a smaller moment but shows an unusual magnetization which is not saturated up to 7 T [Fig. 5(a)]. As comparison, the other compounds are easily saturated around 1–3 T. Its magnetic susceptibility and the extracted C-W temperature are also remarkably distinct from those of the other ones. Whether or not this suggests a rare-earth spin system with much stronger spin interactions, needs to be explored in the future.
cpl-38-4-047502-fig6.png
Fig. 6. Heat capacity measurements for (a) YbOCl/SmSI and (b) ErOF/HoOF/DyOF. The AFM-like transition in DyOF is also evidenced by the peak in heat capacity around 3 K (b). The broad features above 10 K in ErOF and HoOF may be related to rich CEF excitation levels in the two compounds. The origin of the upturn below 10 K in YbOCl (a) is not clear yet and needs to be identified by further experiments down to milikelvins.
As mentioned above, the strong SOC in rare-earth magnetic ions intrinsically provides highly anisotropic spin interactions, which have been well investigated in the triangular lattice QSL candidates.[26] When the anisotropic spin interactions are transplanted on an undistorted honeycomb lattice, we can fully expect the bond-dependent Kitaev interactions coming out of them. In fact, the generic nearest neighbor spin Hamiltonian of rare-earth ions on a honeycomb lattice can be equivalently written as the $J$–$K$–$\varGamma$–$\varGamma'$ model,[24,35] in which the Kitaev term $K$ is clearly highlighted. The transformation relations of exchange parameters between the two models are derived.[24] As discussed above, the next nearest neighbor DM term is symmetrically allowed in our case. For accuracy, the DM term should be taken into account circumstantially and it may bring non-trivial effects in some cases. It has been predicted that the DM term in Kitaev materials can contribute a non-quantized thermal Hall conductance under magnetic fields,[36] which is understood in term of the motion of the emergent spinons twisted by an internal gauge flux. In summary, we have reported the discovery of a rare-earth based family of KSL candidates. Most family members share the SmSI-type structure with the symmetry of $R\bar{3}m$ and rare-earth magnetic ions form a regular honeycomb lattice. The combination of the honeycomb lattice and strongly spin-orbit coupled moments of $4f$ electrons generates bond-dependent anisotropic exchange couplings required by the Kitaev model. The family features a weak van der Waals interlayer coupling which is of particular advantage for exact anyonic excitations and potential applications. We have synthesized some representative family members including the crystals of YbOCl and the polycrystals of SmSI, ErOF, HoOF and DyOF. Their structures are carefully determined by Rietveld refinements. Thermodynamic measurements including magnetic susceptibility, magnetization and heat capacity have been carried out down to 1.8 K, and no obvious phase transitions are found in the samples but DyOF. Thanks to the diversity of the family members, the family is expected to be a promising and manipulable platform for the study of KSL and may host more intriguing spin physics beyond Kitaev physics. Q. M. Zhang acknowledges the support from Users with Excellence Program of Hefei Science Center and High Magnetic Field Facility, CAS. The refined structural data are available on request. References Resonating valence bonds: A new kind of insulator?Anyons in an exactly solved model and beyondSpin-orbit effects in Na 4 Ir 3 O 8 : A hyper-kagome lattice antiferromagnetMott Insulators in the Strong Spin-Orbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev ModelsStructure, and magnetic and electrochemical properties of layered oxides, Li2IrO3Antiferromagnetic Mott insulating state in single crystals of the honeycomb lattice material Na 2 IrO 3 Cu 2 IrO 3 : A New Magnetically Frustrated Honeycomb IridateOn AgRhO2, and the new quaternary delafossites AgLi1/3M2/3O2, syntheses and analyses of real structuresIridium containing honeycomb Delafossites by topotactic cation exchangeA spin–orbital-entangled quantum liquid on a honeycomb latticeLocal magnetism and spin dynamics of the frustrated honeycomb rhodate Li 2 RhO 3 Moessbauer spectroscopy and magnetization studies of .alpha.- and .beta.-ruthenium trichloride α − RuCl 3 : A spin-orbit assisted Mott insulator on a honeycomb latticeMagnetic order in α − RuCl 3 : A honeycomb-lattice quantum magnet with strong spin-orbit couplingAnisotropic Ru 3 +   4 d 5 magnetism in the α − RuCl 3 honeycomb system: Susceptibility, specific heat, and zero-field NMRLow-temperature crystal and magnetic structure of α − RuCl 3 Magnetization processes of zigzag states on the honeycomb lattice: Identifying spin models for α − RuCl 3 and Na 2 IrO 3 Neutron scattering in the proximate quantum spin liquid α-RuCl 3Majorana quantization and half-integer thermal quantum Hall effect in a Kitaev spin liquidNéel-type antiferromagnetic order and magnetic field–temperature phase diagram in the spin- 1 2 rare-earth honeycomb compound YbCl 3 Kitaev Spin Liquid Candidate Os x Cl 3 Comprised of Honeycomb Nano-DomainsKitaev materials beyond iridates: Order by quantum disorder and Weyl magnons in rare-earth double perovskitesGround-state phase diagram of the Kitaev-Heisenberg model on a kagome latticeHoneycomb rare-earth magnets with anisotropic exchange interactionsGapless quantum spin liquid ground state in the two-dimensional spin-1/2 triangular antiferromagnet YbMgGaO4Rare-Earth Triangular Lattice Spin Liquid: A Single-Crystal Study of YbMgGaO 4 Rare-Earth Chalcogenides: A Large Family of Triangular Lattice Spin Liquid CandidatesPreparation, powder data and crystal structure of YbOClStructure cristalline du sulfoiodure de samarium, SmSIZur Hochdruckpolymorphie der Seltenerdsulfidiodide LnSIProduction of the oxyfluorides of holmium, erbium, and ytterbium and study of their crystal-chemical propertiesCrystalline electric field excitations in the quantum spin liquid candidate NaYbSe 2 Crystalline Electric-Field Randomness in the Triangular Lattice Spin-Liquid YbMgGaO 4 Ferromagnetic Kitaev interaction and the origin of large magnetic anisotropy in α-RuCl3Generic Spin Model for the Honeycomb Iridates beyond the Kitaev LimitThermal Hall signatures of non-Kitaev spin liquids in honeycomb Kitaev materials
[1] Anderson P W 1973 Mater. Res. Bull. 8 153
[2] Kitaev A 2006 Ann. Phys. 321 2
[3] Chen G and Balents L 2008 Phys. Rev. B 78 094403
[4] Jackeli G and Khaliullin G 2009 Phys. Rev. Lett. 102 017205
[5] Kobayashi H, Tabuchi M, Shikano M, Kageyama H and Kanno R 2003 J. Mater. Chem. 13 957
[6] Singh Y and Gegenwart P 2010 Phys. Rev. B 82 064412
[7] Abramchuk M, Ozsoy-Keskinbora C, Krizan J W, Metz K R, Bell D C and Tafti F 2017 J. Am. Chem. Soc. 139 15371
[8] Todorova V, Leineweber A, Kienle L, Duppel V and Jansen M 2011 J. Solid State Chem. 184 1112
[9] Roudebush J H, Ross K A and Cava R J 2016 Dalton Trans. 45 8783
[10] Kitagawa K, Takayama T, Matsumoto Y, Kato A, Takano R, Kishimoto Y, Bette S, Dinnebier R, Jackeli G and Takagi H 2018 Nature 554 341
[11] Khuntia P, Manni S, Foronda F R, Lancaster T, Blundell S J, Gegenwart P and Baenitz M 2017 Phys. Rev. B 96 094432
[12] Kobayashi Y, Okada T, Asai K, Katada M, Sano H and Ambe F 1992 Inorg. Chem. 31 4570
[13] Plumb K W, Clancy J P, Sandilands L J, Shankar V V, Hu Y F, Burch K S, Kee H Y and Kim Y J 2014 Phys. Rev. B 90 041112
[14] Sears J A, Songvilay M, Plumb K W, Clancy J P, Qiu Y, Zhao Y, Parshall D and Kim Y J 2015 Phys. Rev. B 91 144420
[15] Majumder M, Schmidt M, Rosner H, Tsirlin A A, Yasuoka H and Baenitz M 2015 Phys. Rev. B 91 180401
[16] Cao H B, Banerjee A, Yan J Q, Bridges C A, Lumsden M D, Mandrus D G, Tennant D A, Chakoumakos B C and Nagler S E 2016 Phys. Rev. B 93 134423
[17] Janssen L, Andrade E C and Vojta M 2017 Phys. Rev. B 96 064430
[18] Banerjee A, Yan J, Knolle J, Bridges C A, Stone M B, Lumsden M D, Mandrus D G, Tennant D A, Moessner R and Nagler S E 2017 Science 356 1055
[19] Kasahara Y, Ohnishi T, Mizukami Y, Tanaka O, Ma S, Sugii K, Kurita N, Tanaka H, Nasu J, Motome Y, Shibauchi T and Matsuda Y 2018 Nature 559 227
[20] Xing J, Feng E, Liu Y, Emmanouilidou E, Hu C, Liu J, Graf D, Ramirez A P, Chen G, Cao H and Ni N 2020 Phys. Rev. B 102 014427
[21] Kataoka K, Hirai D, Yajima T, Nishio-Hamane D, Ishii R, Choi K Y, Wulferding D, Lemmens P, Kittaka S, Sakakibara T, Ishikawa H, Matsuo A, Kindo K and Hiroi Z 2020 Phys. Soc. Jpn. 89 114709
[22] Li F Y, Li Y D, Yu Y, Paramekanti A and Chen G 2017 Phys. Rev. B 95 085132
[23] Morita K, Kishimoto M and Tohyama T 2018 Phys. Rev. B 98 134437
[24] Luo Z X and Chen G 2020 SciPost Phys. Core 3 004
[25] Li Y, Liao H, Zhang Z, Li S, Jin F, Ling L, Zhang L, Zou Y, Pi L, Yang Z, Wang J, Wu Z and Zhang Q 2015 Sci. Rep. 5 16419
[26] Li Y, Chen G, Tong W, Pi L, Liu J, Yang Z, Wang X and Zhang Q 2015 Phys. Rev. Lett. 115 167203
[27] Liu W, Zhang Z, Ji J, Liu Y, Li J, Wang X, Lei H, Chen G and Zhang Q 2018 Chin. Phys. Lett. 35 117501
[28] Brandt G and Diehl R 1974 Mater. Res. Bull. 9 411
[29] Savigny N, Laruelle P and Flahaut J 1973 Acta Crystallogr. B 29 345
[30] Beck H P and Strobel C 1986 Z. Anorg. Allg. Chem. 535 229
[31] Podberezskaya N V, Batsanova L R and Egorova L S 1966 J. Struct. Chem. 6 815
[32] Zhang Z, Ma X, Li J, Wang G, Adroja D T, Perring T P, Liu W, Jin F, Ji J, Wang Y, Kamiya Y, Wang X, Ma J and Zhang Q 2021 Phys. Rev. B 103 035144
[33] Li Y, Adroja D, Bewley R I, Voneshen D, Tsirlin A A, Gegenwart P and Zhang Q 2017 Phys. Rev. Lett. 118 107202
[34] Sears J A, Chern L E, Kim S, Bereciartua P J, Francoual S, Kim Y B and Kim Y J 2020 Nat. Phys. 16 837
[35] Rau J G, Lee E K H and Kee H Y 2014 Phys. Rev. Lett. 112 077204
[36] Gao Y H, Hickey C, Xiang T, Trebst S and Chen G 2019 Phys. Rev. Res. 1 013014