Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 180.
Published online 2017 August 1. doi:  10.1038/s41467-017-00264-6
PMCID: PMC5539315

Robust spin correlations at high magnetic fields in the harmonic honeycomb iridates


The complex antiferromagnetic orders observed in the honeycomb iridates are a double-edged sword in the search for a quantum spin-liquid: both attesting that the magnetic interactions provide many of the necessary ingredients, while simultaneously impeding access. Focus has naturally been drawn to the unusual magnetic orders that hint at the underlying spin correlations. However, the study of any particular broken symmetry state generally provides little clue about the possibility of other nearby ground states. Here we use magnetic fields approaching 100 Tesla to reveal the extent of the spin correlations in γ-lithium iridate. We find that a small component of field along the magnetic easy-axis melts long-range order, revealing a bistable, strongly correlated spin state. Far from the usual destruction of antiferromagnetism via spin polarization, the high-field state possesses only a small fraction of the total iridium moment, without evidence for long-range order up to the highest attainable magnetic fields.


Spin systems with highly anisotropic exchange interactions have recently been proposed to host quantum spin-liquid states1. It has been suggested that the extreme exchange anisotropy required to achieve such states can be mediated by the strong spin–orbit interactions of transition-metal ions situated in undistorted and edge-sharing oxygen octahedra2. Both the two-dimensional (α) and three-dimensional (3D) (β and γ) polymorphs of the honeycomb iridates (A2IrO3, A = Na, Li) closely fulfill these geometric requirements and display intriguing low-field magnetic properties that indirectly indicate the presence of anisotropic exchange interactions between the iridium spins35. Kitaev’s spin-liquid ground state, however, results from the presence of exclusively bond-specific Ising exchange interactions1, which have not yet been verified experimentally. To date, all honeycomb iridates deviate sufficiently from this ideal, such that they transition to long-range magnetic order at finite temperature6, 7 drawing attention to the complex magnetic structures of their ground states813.

γ-lithium iridate features an incommensurate magnetic structure with non-coplanar and counter-rotating moments below T N = 38 K6. The magnetic anisotropy within the ordered phase was extensively characterized in our previous study4. With increasing magnetic field below T N, the magnetic torque τ, divided by the applied field H increases linearly up to an angle-dependent field H *, defining the phase boundary of the low-field ordered state4. Importantly, the sharp feature at H * is not accompanied by full saturation of the iridium moment4. Instead, the magnetic moment at H * is only ~0.1 μ B 4. The lack of a fully saturated moment at H * implies either the onset of another magnetically ordered phase above H *, or alternatively a transition into a paramagnetic state lacking long-range order. The latter implies that the spin correlations are controlled by exchange interactions much stronger than the applied magnetic field. Just as the strange metallic state near the quantum critical point in the hole-doped cuprates is revealed once superconductivity is suppressed with magnetic field14, we use high magnetic fields to destroy the antiferromagnetic order and expose the spin correlations in γ-lithium iridate.

In the following, we show that above H * γ-Li2IrO3 is characterized by a highly anisotropic nonlinear magnetic response with a small net magnetic moment up to the highest attainable magnetic fields. The absence of a sharp boundary with the high-temperature paramagnetic phase suggests that the high-field state does not break additional symmetries and does not exhibit long-range order. Furthermore, a hysteresis-free magnetic anisotropy that can be tied to the crystallographic directions evolves continuously with field out of the ordered state, uncharacteristic of glassy behavior. With these features in mind, and in light of recent studies that find dynamic evidence for liquid behavior in related transition-metal honeycomb structures15, 16, we will refer to the observed high-field state as a spin-fluid.


Experimental technique

We use torque magnetometry to study the high-field state of single crystal γ-lithium iridate. Specifically, we examine the nonlinear response at high magnetic fields. The anisotropy of the magnetic susceptibility α ij = χ i  χ j in the linear response regime leads to a smooth sin 2θ angle dependence with the torque vanishing when field is applied along the high symmetry directions (Supplementary Information of ref. 4). Deviations from this smooth angle dependence provide a direct probe of magnetic correlations. The extreme magnetic anisotropy of γ-lithium iridate necessitated that both smaller volume samples and stiffer cantilevers were employed for the pulsed high-field measurements compared to our prior low-field (DC magnet) measurements4. To this end, we utilized focused ion beam lithography to cut and precisely align samples on Seiko PRC120 piezoresistive levers. This had the added benefit of reducing the lever deflection with field and systematic angle offsets, particularly when field is aligned close to the magnetically hard-axis.

Magnetic anisotropy

Figure 1a, b show τ/H for field rotation in two planes that include the c-axis: the bc-plane (φ=0) and a plane that is perpendicular to one of the honeycomb planes, referred to as the hc-plane (φ=55) (Fig. 1c). We find that H * closely follows a 1∕∣cos θ cos φ angle dependence (Fig. 2a, b), where θ denotes the angle between the ab-plane and the applied field in both rotation sets. The collapse of 1/H * vs. the normalized b-component of magnetic field Hb∣∕∣H∣ = ∣cos θ, cos φ onto a straight line in both rotation planes (Fig. 2c) indicates that the torque at H * is dominated by the b-component of magnetization M b and that the absolute value of magnetization M(H *) is nearly independent of field orientation. Moreover, the absolute value in the denominator of the angle dependence of H * indicates bistable behavior, Mb=Mb*signHb. To determine whether the bistable behavior of M b persists to fields above H *, we turn to the angle dependence of τ/H across the entire field range.

Fig. 1
Crystal structure and anisotropic magnetic behavior of γ-lithium iridate. a, b Magnetic anisotropy, represented by τ/H, as a function of magnetic field H in the bc- and hc-planes at 4 K. A sharp and highly anisotropic feature defined ...
Fig. 2
Angle dependence of H *. a H * vs. rotation angle θ for both the bc- and b hc- rotational planes, with θ defined as the angle between the ab-plane and the magnetic field. c 1/H * plotted vs. |H b|/|H|. The magnitude of 1/H * collapses ...

Below H *, τ/H follows a sin 2θ dependence (Fig. 3a) that is characteristic of the linear response regime, where MiχijHj (Supplementary Information of ref. 4). The sin 2θ dependence below H * is in stark contrast with the angle dependence at high fields, where torque exhibits a sharp discontinuity as the magnetic field crosses the c-axis (Fig. 3b). We observe the nonlinear susceptibility at fields above H * as a sin θ component, accompanied by a sign(cos θ) factor that captures the observed discontinuity at 90°. The sin θ angle dependence indicates that the torque τMbHcMcHb in this high-field regime is dominated by the first term, where H c = H sin θ. The negligible contribution of the second term M c H b is confirmed by the observed torque in the hc-plane: if the discontinuity in τ/H vs. θ is driven exclusively by the saturated b-component of magnetization, then one would expect the amplitude of sin θ in the hc-plane to be reduced by a factor cos 55°  0.577 compared to the bc-plane. As expected, the reduced amplitude factor is 0.573. Thus, M b dominates the total magnetization at very high fields. Futhermore, the high-field magnetic response indicates that the Ising-like discontinuity Mb=Mb*signHb across the c-axis extends well beyond H *.

Fig. 3
Evolution of the angle dependence of τ/H with increasing magnetic field at low temperatures. a At 4 K, τ/H evolves from a characteristic sin angle dependence at small fields to b a correlation-driven sin θ sign(cos ...

In the ultra high-field limit, when all spins are nearly aligned with the applied magnetic field, the effective local anisotropy energy (per formula unit) is EμIr2β2cos2θ, where μ Ir is the magnetic moment of the iridium ion and β is the anisotropic spin stiffness. Therefore when all correlations are overcome by a very large applied magnetic field, the torque τ=dEdθμir2βsin2θ is expected to recover the sin 2θ angle dependence observed at low fields. By contrast, in γ-lithium iridate, the sin θ sign(cos θ) angle dependence of the torque persists up to the highest applied magnetic fields (Fig. 3c), indicating the presence of robust spin correlations throughout the entire field range and providing a lower bound for the magnitude of the exchange interactions.

To assess the extent of this correlated spin-fluid, we examine the thermal evolution of the nonlinear angle dependence of the torque at fields above H * (vertical line in Fig. 4b). Figure 4a shows the angle dependence of τ/H at 15 T for a range of temperatures, revealing a gradual decrease of the nonlinear response as temperature is raised. The nonlinearity becomes undetectable at this magnetic field above 70 K, the same temperature where paramagnetic behavior onsets in magnetic susceptibility measurements4, indicating a continuous crossover to a strongly correlated spin-fluid at low temperatures. This crossover is similarly captured by a smoothly evolving torque signal upon cooling through the transition temperature at fields above H * (Fig. 1d). We note that specific heat measurements also observe a broad peak upon cooling at high fields17, consistent with a decrease in entropy due to the onset of spin correlations18

Fig. 4
The temperature/magnetic field phase diagram and its relation to the spin-anisotropic exchange interactions of γ-lithium iridate. a The temperature evolution of the magnetic anisotropy represented by A/A 2. A and A 2 are the coefficients determined ...


The orientation of the spin-anisotropic exchange interactions with respect to the crystal directions gives rise to a strong magnetic anisotropy in the ordered state4. The extreme softening of χ b occurs because the b-axis is the only direction coaligned either parallel or perpendicular to all Ir-O2-Ir planes. None of the Ir-O2-Ir planes are parallel or perpendicular to either the a- or c-axis (Fig. 4c). In this study, we find that the b-component of magnetization continues to dominate the magnetic response beyond the ordered state to the highest measured fields (Fig. 4d). The high-field anisotropy bears resemblence to the broken symmetry state observed at low temperatures and low fields. This may suggest that the correlated object is a sub-unit of the complex magnetic structure observed at zero field: H * signifies the destruction of long-range order while leaving most of the local magnetic correlations intact. Perhaps the counter-propagating spin spirals seen by X-ray scattering in β-lithium iridate6, are decoupled by a small component of magnetic field along the b-direction, leading to a finite correlation length with only minimal polarization of the individual spins.

The unusual behavior of the spin-fluid directly indicates that the observed magnetic anisotropy is driven by exchange interactions, rather than g-factor anisotropy that is tied to the honeycomb planes (Supplementary Information of ref. 4). This observation, coupled with the anomalously small magnetic moment induced for all field orientations, identifies spin correlations that persist over a broad field and temperature range. In this context, we note that many conventional, as well as correlated, metals are unstable at low temperatures and undergo symmetry breaking that gaps their low-energy excitations19, 20. In frustrated magnets, long-range order can be stabilized by lattice distortions, crystal-field effects, and alternative exchange pathways. In the specific case of embedding the Kitaev model onto a 3D honeycomb lattice, it appears that long-range order is stabilized by intrinsic symmetry breaking—whereby only one component of the spin-anisotropic exchange can be coaligned with a crystal direction. Recent dynamic studies in related compounds15, 16 have found exotic gapless excitations persisting after magnetic order is suppressed with relatively small magnetic fields. Although the ordered state at low temperatures in γ-lithium iridate has quelled the possibility of a Kitaev spin-liquid ground state, we have shown that this system hosts an exotic spin state that is otherwise masked by the zero-field antiferromagnetic order. Other studies, such as nuclear magnetic resonance16, inelastic neutron scattering15 or linear thermal transport in the zero-temperature limit21 are necessary to determine whether the spin-fluid in γ-lithium iridate inherits any properties of a gapless spin-liquid.

Data availability

All relevant data is available from the authors.

Electronic supplementary material

Peer Review File(937K, pdf)


This work was performed at the National High Magnetic Field Laboratory, and was supported by the US Department of Energy BES ‘Science at 100 T’, the National Science Foundation DMR-1157490, and the State of Florida. R.M. acknowledges support from LANL LDRD DR20160085 ‘Topology and Strong Correlations’. J.A. and N.B. acknowledge support by the Department of Energy Early Career program, Office of Basic Energy Sciences under Contract No. DE-AC02-05CH11231. J.A. and N.B. also acknowledge support from the Gordon and Betty Moore Foundation’s EPiQS Initiative through Grant GBMF4374.

Author contributions

Author contributions

K.A.M., B.J.R., J.B.B., R.D.M. and A.S. performed the experiments at the National High Magnetic Field Laboratory-Pulsed Field Facility. K.A.M., B.J.R., R.D.M. and A.S. analyzed the data and wrote the manuscript with input from all authors. N.P.B. and J.G.A. synthesized and characterized high-quality single crystals of γ-lithium iridate at the University of California, Berkeley.


Competing interests

The authors declare no competing financial interests.


Electronic supplementary material

Supplementary Information accompanies this paper at doi:10.1038/s41467-017-00264-6.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Kitaev A. Anyons in an exactly solved model and beyond. Ann. Phys. 2006;321:2–111. doi: 10.1016/j.aop.2005.10.005. [Cross Ref]
2. 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. 2009;102:017205. doi: 10.1103/PhysRevLett.102.017205. [PubMed] [Cross Ref]
3. Chun SH, et al. Direct evidence for dominant bond-directional interactions in a honeycomb lattice iridate Na2IrO3. Nat. Phys. 2015;11:462–466. doi: 10.1038/nphys3322. [Cross Ref]
4. Modic KA, et al. Realization of a three-dimensional spin-anisotropic harmonic honeycomb iridate. Nat. Commun. 2014;5:4203. doi: 10.1038/ncomms5203. [PubMed] [Cross Ref]
5. Takayama T, et al. Hyperhoneycomb iridate β-Li2IrO3 as a platform for Kitaev magnetism. Phys. Rev. Lett. 2015;114:077202. doi: 10.1103/PhysRevLett.114.077202. [PubMed] [Cross Ref]
6. Biffin A, et al. Noncoplanar and counterrotating incommensurate magneticc order stabilized by Kitaev interactions in γ-Li2IrO3. Phys. Rev. Lett. 2014;113:197201. doi: 10.1103/PhysRevLett.113.197201. [PubMed] [Cross Ref]
7. Biffin A, et al. Unconventional magnetic order on the hyperhoneycomb Kitaev lattice in β-Li2IrO3: Full solution via magnetic resonant X-ray diffraction. Phys. Rev. B. 2014;90:205116. doi: 10.1103/PhysRevB.90.205116. [Cross Ref]
8. Rau JG, Kee H-Y. First-principles study of the magnetic structure of Na2IrO3. Phys. Rev. Lett. 2015;90:144417.
9. Sizyuk Y, Price C, Wolfle P, Perkins NB. Importance of anisotropic exchange interactions in honeycomb iridates. minimal model for zigzag antiferromagnetic order in Na2IrO3. Phys. Rev. B. 2014;90:155126. doi: 10.1103/PhysRevB.90.155126. [Cross Ref]
10. Chaloupka J, Jackeli G, Khaliullin G. On the origin of zigzag magnetic order in iridium oxide Na2IrO3. Phys. Rev. Lett. 2013;110:097204. doi: 10.1103/PhysRevLett.110.097204. [PubMed] [Cross Ref]
11. Lee EK-H, Kim YB. Theory of magnetic phase diagrams in hyperhoneycomb and harmonic-honeycomb iridates. Phys. Rev. B. 2015;91:064407. doi: 10.1103/PhysRevB.91.064407. [Cross Ref]
12. Kimchi I, Coldea R, Vishwanath A. Unified theory of spiral magnetism in the harmonic-honeycomb iridates α, β, γ Li2IrO3. Phys. Rev. B. 2015;91:245134. doi: 10.1103/PhysRevB.91.245134. [Cross Ref]
13. Reuther J, Thomale R, Rachel S. Spiral order in the honeycomb iridate Li2IrO3. Phys. Rev. B. 2014;90:100405. doi: 10.1103/PhysRevB.90.100405. [Cross Ref]
14. Ramshaw B, et al. Quasiparticle mass enhancement approaching optimal doping in a high-Tc superconductor. Science. 2015;348:317–320. doi: 10.1126/science.aaa4990. [PubMed] [Cross Ref]
15. Banerjee A, et al. Neutron scattering in a proximate quantum spin liquid. Science. 2017;356:1055–1059. doi: 10.1126/science.aah6015. [PubMed] [Cross Ref]
16. Baek, S.-H. et al. Observation of a field-induced quantum spin liquid in α-RuCl3, Preprint at (2017).
17. Ruiz, A. et al. Field-induced intertwined order in 3D Mott-Kitaev honeycomb β-Li2IrO3, Preprint at (2017).
18. Kohama Y, et al. Field-induced XY and Ising ground states in a quasi-two-dimensional S = 1/2 Heisenberg antiferromagnet. Phys. Rev. B. 2011;84:184402. doi: 10.1103/PhysRevB.84.184402. [Cross Ref]
19. Matthias BT, Geballe TH, Compton VB. Superconductivity. Rev. Mod. Phys. 1963;35:1. doi: 10.1103/RevModPhys.35.1. [Cross Ref]
20. Keimer B, Kivelson SA, Norman MR, Uchida S, Zaanen J. From quantum matter to high-temperature superconductivity in copper oxides. Nature. 2015;518:179–186. doi: 10.1038/nature14165. [PubMed] [Cross Ref]
21. Yamashita M, et al. Thermal-transport measurements in a quantum spin-liquid state of the frustrated triangular magnet κ-(bedt-ttf)2Cu2(CN)3. Nat. Phys. 2009;5:44–47. doi: 10.1038/nphys1134. [Cross Ref]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group