Robust quantum spin liquid state in the presence of giant magnetic isotope effect in D3LiIr2O6

Robust quantum spin liquid state in the presence of giant magnetic isotope effect in D3LiIr2O6

Introduction

Quantum spin liquid (QSL) represents a novel quantum state of matter where interacting spins remain disordered down to zero temperature and harbors a topologically nontrivial state with fractionalized excitations1,2. In two-dimensional magnets, since the proposal of the resonating valence bond (RVB) state by P. W. Anderson for the S = 1/2 Heisenberg antiferromagnet on a triangular lattice3, a variety of theoretical models has been put forward to realize a QSL state. Among those, the S = 1/2 bond-dependent Ising model on a honeycomb lattice proposed by A. Kitaev4 is distinct from the other two-dimensional models in that it is exactly solvable. By introducing two kinds of Majorana operators4, the ground state of the Kitaev model was shown to be a QSL consisting of an ensemble of itinerant and localized Majorana fermions. The Kitaev model and its QSL ground state have been attracting considerable interest recently.

The essential ingredient of the Kitaev model is the competing Ising interactions with their easy-axis orthogonal to each other on the three 120° bonds of honeycomb lattice4. Such bond-dependent Ising couplings have been theoretically proposed to be realized in iridium (Ir4+) oxides (iridates) and a ruthenium (Ru3+) chloride with an edge-shared honeycomb network of IrO6 and RuCl6 octahedra5. In the family of Ir4+ oxides with five 5d electrons, the cubic crystal field overcomes Hund’s first rule, and all the five 5d-electrons are accommodated in the t2g manifold. Strong spin-orbit coupling λSO of Ir, as large as 0.5 eV, entangles spin and orbital degrees of freedom, yielding a fully-filled lower Jeff = 3/2 quartet and a half-filled upper Jeff = 1/2 doublet out of the triply-degenerate t2g orbitals6. A spin-orbit-entangled Mott state with localized Jeff = 1/2 pseudospins can be stabilized only with a modest on-site Coulomb interaction of a few eV between 5d electrons7. When the two neighboring IrO6 octahedra are connected by sharing their edges, bond-dependent ferromagnetic Ising coupling between the two Jeff = 1/2 pseudospins appears through a super-exchange process in the Ir-O2-Ir planar bonds5. If the other magnetic couplings than this super-exchange coupling can be ignored, the Kitaev model is realized in the edge-shared honeycomb network of IrO6 octahedra, leading to a quantum liquid state of spin-orbit-entangled pseudospins (denoted QSL hereafter for brevity).

The layered iridates Na2IrO3 and α-Li2IrO38,9 and α-RuCl310 with edge-shared honeycomb networks of IrO6 (RuCl6) octahedra, and later their three-dimensional analogs β– and γ-Li2IrO311,12, have been considered as candidates to realize the Kitaev model and the QSL state. However, all those candidate materials were found to undergo either a zigzag or a spiral magnetic ordering at low temperatures13,14,15,16. The appearance of magnetic ordering rather than a QSL state is most likely attributable to the superposition of other types of magnetic interactions on the Kitaev coupling, such as antiferromagnetic Heisenberg coupling due to direct dd overlap, symmetry-allowed off-diagonal exchange and further-neighbor interactions17,18,19. Nevertheless, the predominance of Kitaev coupling and the proximity to the Kitaev QSL state has been argued both experimentally and theoretically20,21,22.

In contrast to the aforementioned candidates, the protonated honeycomb iridate H3LiIr2O6 was discovered to host a QSL ground state23. H3LiIr2O6 is obtained by a soft-chemical substitution on α-Li2IrO324,25. H+ ions replace all of the Li+ ions in the interlayer sites but not those in the honeycomb layers of α-Li2IrO3, which leaves the LiIr2O6 layer units intact. The Curie-Weiss behavior of magnetic susceptibility at high temperatures evidences the presence of localized and interacting spin-orbit-entangled Jeff = 1/2 moments. The negative Curie-Weiss temperature, θCW ~ −100 K, points to the sizable antiferromagnetic couplings23, very likely the superposition of more than one kind of interaction, such as Kitaev, Heisenberg, and off-diagonal interactions. Despite those magnetic interactions, H3LiIr2O6 does not show any sign of magnetic ordering or spin-glass freezing in the specific heat and the nuclear magnetic resonance (NMR) measurements down to 50 mK, which is indicative of a QSL ground state. However, any clear evidence for the Kitaev-type QSL, including the signature of localized and itinerant Majorana fermions, has been lacking so far. The specific heat at low temperatures below 1 K was found to be dominated only by a T3-term, which very likely represents the lattice contribution, indicating a gap in the spin excitations23. This suggests that the QSL state of H3LiIr2O6 is at least away from the pure Kitaev limit, and some mechanism suppressing competing long-range magnetic order should be operative behind the QSL state.

The formation of the QSL state in H3LiIr2O6, in sharp contrast to the other honeycomb iridates and the ruthenium chloride, brought our attention to the role of interlayer protons, a structural ingredient distinct from the other honeycomb iridates. We examined the H/D isotope effect to experimentally capture the role of the OH bond by synthesizing the deuterium-substituted compound D3LiIr2O6. The structural refinements based on the neutron diffraction measurements revealed the presence of covalent OD bonds with oxygen atoms either above or below the deuterium plane, as well as the large structural isotope effect in the OD and OH bonds. A gigantic isotope effect on the magnetic interactions was discovered, indicating the unexpectedly large influence of OD/OH bonds on the magnetic couplings. The antiferromagnetic Curie-Weiss temperature |θCW| of D3LiIr2O6 was found to be enhanced substantially to ~170 K from ~100 K of H3LiIr2O6. Nevertheless, D3LiIr2O6 was found to host a QSL state of Jeff = 1/2 pseudospins as in H3LiIr2O6. The contrast of the OD/OH bonds should account for the large isotope effect on the magnetic interactions governed by the super-exchange process through the oxygen 2p orbitals. Random distribution of the O-D-O (or O-H-O) bond structures in the two materials very likely gives rise to bond-disorder in the in-plane magnetic interactions of the Jeff = 1/2 moments, which we argue may be one of the important ingredients to stabilize the QSL state in D3LiIr2O6 and H3LiIr2O6.

Results

Crystal structure

The powder sample of D3LiIr2O6 was synthesized from α-Li2IrO3 by an ion-exchange technique as described in “Methods”. The powder X-ray diffraction (PXRD) indicated that the product comprised a single phase of D3LiIr2O6 crystallizing in the same layered honeycomb structure as H3LiIr2O6 (see Supplementary Fig. S1 for the PXRD pattern). The presence of strong stacking fault disorder was identified in the PXRD pattern as in H3LiIr2O624. All the interlayer Li+ ions of α-Li2IrO3 are replaced by D+ ions within the given resolution of our experiments. Unlike α-Li2IrO3, O2- ions in the top and the bottom layers of the adjacent LiIr2O6 honeycomb layer face each other and form nearly straight O-D-O units24.

The comparison of the PXRD pattern with that of H3LiIr2O6 indicates a large structural isotope effect. The shift of 001 reflections is clearly seen in the PXRD pattern (Fig. 1b), which corresponds to the increased distance between the Ir-honeycomb planes (d001) by ~1% via deuteration. The lattice parameter c at room temperature estimated from the neutron diffraction measurements described below is c = 4.884(2) Å for H3LiIr2O6 and c = 4.916(6) Å for D3LiIr2O6, respectively. A similar gigantic structural isotope effect from the interlayer H+/D+ ions was observed in layered oxides such as H(D)CrO226,27, which has been ascribed to the distinct character of the interlayer O-H-O and O-D-O bonds. The H+ or D+ ion placed between the two O2- ions experiences double-well potential. While the D+ ion is very likely trapped at one of the potential wells, the H+ ion is argued to delocalize over the two possible positions due to the pronounced quantum tunneling effect26,27,28.

Fig. 1: Crystal and electronic structure of D3LiIr2O6.
Robust quantum spin liquid state in the presence of giant magnetic isotope effect in D3LiIr2O6

a Crystal structure of D3LiIr2O6 (space group C2/m) refined from the neutron diffraction data. The gray, blue, light green, and yellow spheres show iridium, oxygen, lithium, and deuterium atoms, respectively. Two different O-D-O bonds are indicated by dotted circles. There are two possible positions for each D atom, and the occupancy of each position is 0.5. The refined structural parameters are listed in Table 1. b Powder X-ray diffraction (PXRD) patterns for D(H)3LiIr2O6 near 001 reflection, highlighting the structural isotope effect. The PXRD patterns were recorded using Cu-Kα radiation. The distance of (001) planes, d001, estimated from PXRD, is ~4.559 Å for D3LiIr2O6 and ~4.526 Å for H3LiIr2O6, respectively. c Difference Fourier map of O-H-O interlayer bonds in H3LiIr2O6 (left) and the O-D-O bonds in D3LiIr2O6 (right) from the neutron diffraction measurements. While the D atoms have a large occupational probability close to the oxygen atoms, the H atoms reside near the center of O-O bonds. This may reflect the more dynamic character of H+ ions due to the quantum tunneling effect. d Ir L3-edge RIXS spectra for D(H)3LiIr2O6 recorded at T = 10 K. The dotted line shows the spectrum of α-Li2IrO3 obtained in the same experimental condition for comparison. The spectra are normalized by the intensity of elastic lines.

Full size image

Such distinct characteristics of the interlayer bonds may also be present in H3LiIr2O6 and D3LiIr2O6. The experimental determination of the position of D+ or H+ ions, however, was not feasible in the PXRD because of its weak X-ray scattering factor. To determine their positions, the neutron diffraction measurement was conducted on the powder sample enriched with 7Li- and 193Ir-isotopes29. The use of 7Li and 193Ir reduces the absorption of neutrons by a factor of ~4. The refined structural parameters for D3LiIr2O6 are listed in Table 1 (see Supplementary Fig. S2 for the powder neutron diffraction patterns), and the corresponding crystal structure is illustrated in Fig. 1a. There are two crystallographic sites for deuterium, D1 and D2, forming the O1-D1-O1 and O2-D2-O2 bonds between the honeycomb layers with the two inequivalent oxygen sites O1 and O2, respectively. The O1-D1-O1 bond is nearly perpendicular to the honeycomb plane. The configuration of the O2-D2-O2 bond is essentially the same as that of O1-D1-O1 but more tilted with respect to the honeycomb plane. We will not distinguish the two bonds explicitly in this paper. The refinement indicates that each D ion has two possible off-center positions, above and below the center, between the two oxygen atoms facing each other (see Supplementary Note 1 for details of the refinements). This contrasts with the symmetric positions of interlayer cations in other delafossite-type oxides such as Ag3LiIr2O630. The estimated occupancies of the two off-center positions are 50% each, which naturally means that only one of the two off-center positions between two O atoms is randomly occupied by D atoms. The simultaneous occupation of the two off-center positions is highly unlikely due to electrostatic repulsion. D atoms form a short OD bond of ~1 Å with one of the O atoms, as short as that in D2O31. The distance to the other oxygen atom is as large as ~1.6 Å. These strongly suggest that a D atom between the two oxygen atoms forms a rigid covalent bond with one and a hydrogen bond with the other, such as O-D···O or O···D-O, where – and ··· denote a covalent and a hydrogen bond, respectively. The O-D···O and O···D-O bonds are distributed randomly, and no signature of layer ordering of O-D···O and O···D-O bonds in the interlayer plane was found in the refinement.

Table 1 Structural parameters of D3LiIr2O6 refined by the neutron diffraction at room temperature
Full size table

The refinement of neutron diffraction data on the isotope-enriched H37Li193Ir2O6 (Supplementary Fig. S3) also indicates the two off-center positions for H+ ions as in the deuterated analog. The refined structural parameters are listed in Supplemental Table S1. Nonetheless, we found signatures of the different characters of O-H-O bonds. The difference Fourier maps for the O-H-O and O-D-O bonds depicted in Fig. 1c highlights the contrast of the two bonds. The Fourier map for the O-D-O bond shows a large occupational probability of the D ions close to either of two O atoms, whereas that of the O-H-O bonds showed that H+ ions reside close to the center of O-O bonds. This difference is manifested in the slightly shorter length of the O-D bonds (d(O1-D1) = 1.001(8) Å, d(O2-D2) = 1.024(12) Å) than the O-H bonds (d(O1-H1) = 1.017(17) Å, d(O2-H2) = 1.026(10) Å) and naturally explains the large structural isotope effect. In addition, the refined isotropic displacement parameter for H atom, Uiso ~ 0.0177 Å2, is larger than those of D atom (Uiso ~ 0.0124 Å2). These suggest that H+ ions have more dynamical character and are possibly tunneling between two off-center positions, as discussed in HCrO226,27. In the case of H3LiIr2O6, however, the strong stacking fault disorder24 should disturb the coherent alignment of oxygen atoms above and below, which may hinder the tunneling and trap a H+ ion at one of the off-center positions inhomogeneously.

Compared to the large isotope effect on the interlayer bonds, the impact of deuteration on the LiIr2O6 honeycomb layer is not as significant as the interlayer distance. Both H3LiIr2O6 and D3LiIr2O6 have strongly compressed IrO6 octahedra. The Ir-O-Ir angles of ~100° are nearly the same in both compounds, as the elongation of the c-axis in D3LiIr2O6 is mostly attributed to the enlarged O-O distance between the LiIr2O6 honeycomb layers by ~2%. The in-plane lattice parameter b(a) increases(decreases) only by ~0.1(0.05)% by deuteration, respectively. Therefore, the local electronic structure should not be modified substantially by the isotope exchange, which was confirmed by resonant inelastic X-ray scattering (RIXS) measurement described below.

Electronic and magnetic properties

The resistivity measurement showed that D3LiIr2O6 is an insulator with an activation gap of ~0.1 eV as shown in the inset of Fig. 2a, a typical value for an Ir4+ Mott insulator23. Magnetic susceptibility χ(T) displays a Curie-Weiss behavior at high temperatures (Fig. 2a). The Curie-Weiss fit above 250 K, shown in Fig. 2b, yields an effective moment of ~1.69 μB/Ir, which is very close to that of Jeff = 1/2 moment. These provide the support that D3LiIr2O6 is a spin-orbit-entangled Mott insulator with localized Jeff = 1/2 pseudospins as in H3LiIr2O6.

Fig. 2: Magnetic and thermodynamic properties of D3LiIr2O6.
figure 2

a Temperature dependence of magnetic susceptibility χ(T) measured at a magnetic field of 1 T. The inset shows the Arrhenius plot of resistivity ρ(T). b Inverse of χ(T) for D3LiIr2O6 and H3LiIr2O6 displaying the isotope effect on the magnetic interactions. The Curie-Weiss temperature changes from ~ −100 K for H to ~ −170 K for D. c 7Li Knight shift K(T) of D3LiIr2O6 and H3LiIr2O6 measured at a magnetic field of 2 T. The broad hump of K(T) appears at a higher temperature in D3LiIr2O6 than that in H3LiIr2O6. d Specific heat divided by temperature C(T)/T for D3LiIr2O6 at zero magnetic field. The nuclear Schottky contributions at low temperatures were removed by the method described in ref. 23. No magnetic order was found down to 50 mK.

Full size image

The Ir L3-edge RIXS spectrum shown in Fig. 1d is consistent with the Jeff = 1/2 picture. D3LiIr2O6 exhibits two excitation peaks at ~0.78 eV and ~1.04 eV, which can be ascribed to an intra-atomic excitation from Jeff = 1/2 to 3/2 states (in the hole description). A similar peak has been observed around ~0.7 eV in the RIXS spectra of a number of iridates32,33, which corresponds to the split of Jeff = 1/2 and 3/2 states, namely 3/2λSO. The double-peak structure observed in D3LiIr2O6 is very likely ascribed to the split of Jeff = 3/2 states by the trigonal crystal field associated with the compression of IrO6 octahedra along the c-axis. The comparison of the double-peak structure with those in α-Li2IrO3 and H3LiIr2O6 indeed indicates the intimate correlation between the peak split in the RIXS spectra and the trigonal structural distortion. While D3LiIr2O6 with a large trigonal distortion shows a clear double-peak structure with the splitting of ~0.25 eV, α-Li2IrO3 with a moderate trigonal distortion shows only broadening within the experimental resolution. Comparing the spectra of D3LiIr2O6 and H3LiIr2O6, the splitting appears to be not appreciably different, at most only marginally larger for D3LiIr2O6 (see Supplementary Note 2 and Fig. S4). This is consistent with the comparable or only slightly larger trigonal distortion for D3LiIr2O6.

While the deuterium isotope exchange does not have any appreciable impact on the electronic structure on eV scale, it gives rise to an unexpectedly large change of the low-energy magnetic interactions, which can be seen clearly in the comparison of magnetic susceptibility χ(Τ) of D3LiIr2O6 with that of H3LiIr2O6 (Fig. 2b). From the Curie-Weiss fit above 250 K, an antiferromagnetic Curie-Weiss temperature, θCW ~ −170 K, is estimated for D3LiIr2O6, which is almost a factor of two larger than that of H3LiIr2O6, θCW ~ −100 K23, pointing to much enhanced antiferromagnetic interactions. With decreasing temperature below 250 K, χ(T) displays a broad hump at around 180 K, which can be regarded as a temperature scale for the development of antiferromagnetic correlations. Consistent with the enhancement of |θCW | , the hump temperature of 180 K is much higher than ~130 K for H3LiIr2O6. The enhancement of hump temperature in the uniform susceptibility can be recognized more clearly in the Knight shift K(T) from 7Li-NMR (Fig. 2c), where Curie-like contributions from spin defects are absent.

The information about spin dynamics obtained by the spin-lattice relaxation time T1 of NMR corroborates the gigantic H/D isotope effect on χ(T) and K(T) (Fig. 3c, d). On cooling from room temperature, (T1T)-1 for 7Li shows a slight increase in accord with the Curie-Weiss behavior of χ(T) and then shows a broad hump around 230 K followed by a decrease. The (T1T)-1 hump temperature is appreciably higher than that of H3LiIr2O6, ~180 K23. As 7Li atoms reside at the center of the Ir-honeycomb lattice and the nuclei experience hyperfine fields from the six surrounding Ir4+ ions, the decrease of (T1T)-1 very likely derives from the cancelation of the hyperfine fields by the development of antiferromagnetic correlations over an Ir hexagon. Hence, the higher (T1T)-1 hump temperature of D3LiIr2O6 also indicates enhanced antiferromagnetic interactions upon deuteration.

Fig. 3: Spin-liquid state of D3LiIr2O6 evidenced by NMR measurement.
figure 3

a 7Li- and b 2D-NMR spectra. The solid vertical line is a guide for spectral positions without any internal fields. The split of spectra in 2D-NMR originates from the nuclear electric quadrupole splitting for a nuclear spin I = 1. The red solid line shows a powder pattern fitted with a single set of quadrupole parameters, νQ = 0.225 MHz and η = 0, where νQ equals to 3/2 × e2qQ/h for I = 1 (e2qQ/h is the deuterium quadrupole coupling constant) and η is the asymmetry parameter, respectively61. The shoulder-like structure is attributable to the quadrupolar interactions in the randomly oriented powder sample. c Temperature dependences of the inverse of spin-lattice relaxation time divided by temperature, (T1T)−1, for 7Li atom in D3LiIr2O6 and H3LiIr2O623. A broad hump of (T1T)-1 is seen at ~230 K for D3LiIr2O6, which is higher than ~180 K for H3LiIr2O6, suggesting the enhanced antiferromagnetic interactions. d Temperature dependence of (T1T)−1 for 2D. The (T1T)-1 data of 1H in H3LiIr2O623 are shown for comparison. Magnetic fields are applied parallel to the honeycomb plane of the field-oriented powder samples.

Full size image

A more significant isotope effect was found in (T1T)-1 of 2D/1H (Fig. 3d). Each D(H) atom is connected to an oxygen atom forming an Ir-O2-Ir plaquette. (T1T)-1 of D(H) should be, therefore, sensitive to the nearest neighbor Ir-Ir magnetic correlation. In H3LiIr2O6, (T1T)-1 of 1H shows a monotonic increase below 100 K, unlike that of 7Li, implying a much weaker antiferromagnetic or possibly ferromagnetic correlation between the nearest Ir atoms23. In contrast, (T1T)-1 of 2D in D3LiIr2O6 exhibits a decrease with lowering temperature below 150 K and shows a dip at ~50 K as in 7Li, suggestive of the appreciable nearest neighbor antiferromagnetic correlations, which may consistently reflect the enhanced antiferromagnetic interactions in D3LiIr2O6. We note here that the positions of H/D atoms are slightly different in H(D)3LiIr2O6 (Fig. 1c), which could influence the hyperfine field at the H/D site and make the contrast more pronounced.

Spin-liquid behavior of D3LiIr2O6

Despite the gigantic isotope effect on the magnetic interactions between the Jeff = 1/2 moments, the ground state of D3LiIr2O6 was found to be a QSL as in H3LiIr2O6. On cooling further below the hump temperature, χ(T) in Fig. 2a shows a Curie-like upturn at low temperatures below 20 K. There was no signature of magnetic order or spin-glass freezing down to 2 K. No magnetic order was observed in specific heat C(T) down to 50 mK either, as displayed in Fig. 2d. The absence of magnetic order/spin-glass freezing is corroborated by the NMR measurements. Figure 3a, b show the 7Li and 2D-NMR spectra down to 2 K, respectively. Both spectra remain sharp without showing any splits or appreciable broadening on cooling, excluding the appearance of magnetic order or a spin-glass state. These results indicate that the QSL state is robust against the large modification of magnetic couplings owing to the H/D isotope effect.

K(T) is nearly temperature-independent below 60 K, as seen in Fig. 2c, indicating that the low-temperature upturn of χ(T) originates predominantly from magnetic impurities or defects. We note that the finite K(0) does not necessarily indicate gapless excitations in the presence of strong spin-orbit coupling and Dzyaloshinskii-Moriya interaction34,35. The analysis of low-temperature specific heat C(T) for H3LiIr2O6 provides a strong indication of the existence of a magnetic excitation gap23.

The low-temperature specific heat C(T, B) with and without applying magnetic field for D3LiIr2O6 shows strikingly similar behavior to those of H3LiIr2O6 as shown in Fig. 4. It is hard to distinguish the data for the two materials, in contrast to the gigantic H/D isotope effect on K(T) and χ(T). At zero magnetic field, C(T)/T displays a weak power-law increase C/TT-1/2, which points to abundant low-energy excitations. By application of magnetic fields, C(T)/T at low temperatures is strongly reduced, suggesting an opening of a gap. The C(T, B) of H3LiIr2O6 at low temperatures was found to follow the B/T scaling behavior, B1/2 C(T)/T = F(T/B)23, which should hold for that of D3LiIr2O6 as well. After the subtraction of the scaled contributions, C(T) of H3LiIr2O6 shows almost field-independent T3-behavior, which likely originates from the lattice, indicative of the absence of low-energy magnetic excitations and thereby a sizeable excitation gap23. The analogous C(T) behavior of D3LiIr2O6, regardless of distinct magnetic exchange parameters, supports that the low-temperature C(T, B) is dominated by spin defects of ~1% isolated from the QSL background and the magnetic contribution to C(T) from the bulk spin-liquid state is negligibly small below ~5 K. This also suggests that D3LiIr2O6 and H3LiIr2O6 share the similar disorder despite the contrast of OD and OH bonds. The low-temperature C(T, B) of H3LiIr2O6 has been argued to be in line with those of Kitaev spin liquid with bond-disorder or spin vacancies36,37,38,39.

Fig. 4: Low-lying magnetic excitation of D3LiIr2O6.
figure 4

Low-temperature specific heat C(T, B)/T of D3LiIr2O6 under magnetic fields B of 0, 2, and 7 T. The data for H3LiIr2O6 are also shown23. The low-temperature C(T)/T is strongly suppressed by B in both D3LiIr2O6 and H3LiIr2O6, indicating a similar form of low-lying excitations.

Full size image

Discussion

The gigantic isotope effect on magnetic couplings is clearly linked to that on the O-D(H)-O interlayer structural unit. The formation of a covalent O-D(H) bond, unique to D(H)3LiIr2O6 among the honeycomb iridates, shall reduce the Ir 5 d-O 2p hybridization and hence suppress the super-exchange coupling between the Jeff = 1/2 moments40. If the magnetic exchange coupling is dominated by the dpd hopping processes, the distinct character of OD and OH bonds (Fig. 1c) naturally accounts for the giant isotope effect on the magnetic interaction. The theoretical calculations incorporating the effect of OD(OH) bonds indeed support the impact of OD(OH) bonds on magnetic couplings40,41,42 as well as the large difference between the OD and OH bonds40. We argue that this is firm evidence that the super-exchange process through oxygen is dominant in the interaction of Jeff = 1/2 moments of honeycomb iridates, which is a necessary condition but not a sufficient condition to realize the Kitaev model.

D3LiIr2O6 and H3LiIr2O6 are distinct from the other honeycomb Jeff = 1/2 magnets in that the ground state is a QSL instead of a magnetically ordered state. The early theoretical estimate of interaction parameters predicted that Kiteav interaction becomes largest when the Ir-O-Ir bond angle is ~100°43,44. However, the estimated magnetic exchange parameters for H3LiIr2O6 seem to be away from the very narrow region of the spin-liquid phase located near the pure Kitaev limit in the framework of the extended Kitaev model, including Heisenberg, off-diagonal, and further-neighbor interactions18,40,41,42,45. Considering the large difference in the exchange parameters between D3LiIr2O6 and H3LiIr2O6 due to the gigantic isotope effect, it is highly unlikely that the parameters for both D3LiIr2O6 and H3LiIr2O6 fall into the narrow parameter region for the QSL state. Additional ingredient(s) must be incorporated to stabilize the QSL state in these two compounds.

The disorder effect associated with OD bonds should be one of the important ingredients for the suppression of long-range magnetic order and the realization of the QSL state. Within the experimental resolution, the D+ ions fully substitute the interlayer Li+ ions of the parent α-Li2IrO3, and no signature of Li-D inter-site mixing was identified in the honeycomb layers. The dominant source of disorder in D3LiIr2O6 should be the random distribution of OD bonds. If we have a mixture of four different configurations of Ir-O2-Ir edge-shared bonds, namely Ir-O(OD)-Ir, Ir-(OD)O-Ir, Ir-O2-Ir and Ir-(OD)2-Ir, we expect to have essentially three different kinds of exchange paths, which distribute randomly and give rise to strong bond-disorder in the magnetic interactions. The random distribution of OD bonds above and below a honeycomb layer gives rise to disorder in the local crystal field acting on the Ir site and, hence, site-disorder as well.

How can such bond- and site-disorder lead to the QSL state in D3LiIr2O6? The theoretical studies pointed out that the Kitaev spin liquid is robust against the presence of disorders such as spin vacancies and bond-disorder40,41,42,46. However, this does not mean that disorder stabilizes the Kitaev spin-liquid state, and competing long-range orders that mask the spin-liquid state in other Kitaev materials need to be evaded. The appearance of a spin-liquid-like state by suppressing long-range order with the disorder has been discussed in geometrically frustrated systems47,48,49. It is conceivable that bond/site-disorder associated with OD bonds destabilizes competing long-range ordering in D3LiIr2O6. The disorder of OD bonds modulates not only the Kitaev exchange but also affects non-Kitaev interactions that would stabilize a long-range order, such as off-diagonal Γ coupling and particularly the Heisenberg coupling; when the two Ir-O-Ir super-exchange paths are asymmetric in the Ir-O(OD)-Ir bond, the destructive interference in the super-exchange paths is invalidated, and the sizeable Heisenberg coupling may be restored40. It is tempting to infer that such disordered non-Kitaev interactions destabilize competing magnetic orders while the emergent spin-liquid state is robust against the bond-disordered Kitaev interactions. The recent RIXS study on H3LiIr2O6 showed momentum-independent magnetic continuum excitations, supportive of a spin-liquid state associated with Kitaev-type coupling50. We, therefore, argue that the bond-disorder and frustrated Kitaev coupling concertedly play a role in realizing the spin-liquid state. Indeed, the theoretical Kitaev-Heisenberg model, including long-range interactions, indicates that the phase range of Kitaev spin liquid is extended in the presence of bond-disorder51.

For H3LiIr2O6, if H+ ions are fluctuating between the two off-center positions due to the quantum tunneling effect, all the Ir-Ir bonds would be regarded as identical, and no bond- or site-disorder would be present. However, the time scale of H+ motions was reported to be much slower compared to that of magnetic exchange52, and H+ ions should be considered static in the magnetic exchange process. In addition, strong stacking faults of the honeycomb layers may trap and localize part of H+ ions. Analogous bond/site-disorder should be present in H3LiIr2O6 as well.

To conclude, the deuterium-substituted honeycomb iridate D3LiIr2O6 realizes a quantum liquid state of spin-orbit-entangled Jeff = 1/2 pseudospin as in H3LiIr2O6. The isotope exchange from hydrogen to deuterium induces a pronounced structural modification of the O-D(H)-O interlayer unit and results in the drastic change of magnetic interactions through the super-exchange process involving O 2p state. Nevertheless, the QSL ground state is robust against the isotope effect. We argue that the disorder of OD (or OH) bonds, incorporated in a frustrated Kitaev material, plays a key role in realizing the quantum liquid state. We also expect that the variant super-exchange interactions associated with OH(D) bonds, as well as the bond/site-disorder effect, may be exploited to realize exotic magnetic ground states not only in Kitaev materials but in a wide variety of magnets53,54.

Methods

Sample synthesis and characterizations

Powder sample of D3LiIr2O6 was prepared by an ion-exchange reaction on α-Li2IrO3. The powder of α-Li2IrO3 was sealed in a Teflon-lined autoclave together with 4 M D2SO4 solution diluted by D2O, and heated at 125 °C for 120 h. The obtained powder was rinsed with D2O and dried in air. For the neutron diffraction study, α7Li2193IrO3 was synthesized from 7Li2CO3 and metallic 193Ir, and the ion-exchange reaction was conducted on it. Transport, magnetic, and thermodynamic properties above 0.4 K were measured on a cold-pressed pellet by commercial apparatus (Quantum Design, MPMS-XL, and PPMS). The heat capacity measurement was performed by a relaxation method.

Powder X-ray diffraction

The crystal structure of the product was analyzed by powder X-ray diffraction. The diffraction patterns were collected at room temperature with Bruker D2 Phaser with Cu-Kα radiation. The detailed measurement shown in Fig. S1 was performed by using a Stoe Stadi-P transmission diffractometer (primary beam Johann-type Ge (111) monochromator for Ag-Kα1 radiation (λ = 0.55941 Å), Mythen-Dectris PSD with 12° 2θ opening) with the sample sealed in a glass capillary of 0.5 mm diameter (Hilgenberg, glass No. 50). The sample was spun during measurement for better particle statistics in both cases.

Neutron diffraction

The powder neutron diffraction measurement was performed at room temperature on the HRPD beamline at the ISIS Neutron and Muon Source. The powder was contained in a 6 mm diameter cylindrical vanadium can with data being collected using the 10–110 ms time-of-flight window. Isotope-enriched powder samples of D37Li193Ir2O6 and H37Li193Ir2O6 were prepared for neutron diffraction to reduce the absorption of neutrons to allow the refinement of fine structural details. The Rietveld refinement of crystal structure was performed using the GSAS-II software55. The crystal structure is visualized by using VESTA software56.

Low-temperature specific heat measurement

Specific heat C between 0.05 and 1 K was measured by using a cell installed in a dilution refrigerator57. The sample was wrapped by an annealed silver foil (residual resistivity ratio, RRR > 1000) and fixed by silver paste to ensure enough thermal link between the sample and the sample stage. A commercial setup (Quantum Design PPMS) was used above 0.4 K. The data was analyzed by the relaxation method. The subtraction of nuclear Schottky contributions was performed with the procedure described in ref. 23.

NMR measurements

NMR measurements were carried out with a standard coherent pulse spectrometer, cryo-cooled preamplifier, double-rotation probe, and 3He-circulation lines. The NMR frequency-swept spectra were acquired using a Fourier-step-sum technique, where each Fourier-transformed spectrum shifted with its center frequency are accumulated during the sweep. The quadrupole splitting for 7Li nuclei (I = 3/2) was very small, less than 20 kHz, throughout the experiment. The spin-lattice relaxation rate 1/T1 for 7Li was measured by a comb-shaped pulse recovery method with an empirical stretched-exponential recovery function, 1 − exp{−(t/ T1)β}, while that for 2D (I = 1) was obtained with the formula, 1 − {0.75exp(−(3t/T1)β) + 0.25exp(−(t/T1)β)}. The field-oriented samples were used for the measurements23.

Resonant inelastic X-ray scattering

(RIXS): The RIXS experiment at the Ir L3-edge was performed at BL11XU of SPring-8. The energy of incident X-rays was tuned at 11.214 keV, which is about 3 eV below the peak energy of the white line in the X-ray absorption spectrum at the Ir L3-edge. This energy corresponds to the excitation of Ir 2p3/2 to 5d t2g. The incident X-ray beam was monochromatised by a double-crystal Si(111) monochromator and a secondary Si(844) four-crystal monochromator. The scattered X-rays were analyzed by a diced and spherically bent Si(844) crystal. The scattering angle was kept at 90 degrees, and π-polarized incident X-rays were used to minimize the elastic scattering. The measurement was performed at 10 K, and the total energy resolution was about 50 meV. The cold-pressed pellets of D(H)3LiIr2O6 powder were used for the measurement. The spectrum of α-Li2IrO3 sintered pellet was also measured for comparison. The obtained spectra can be regarded as a q-averaged one. Since the dd excitations do not show sizable q-dependence31, the spectra should reflect the essential features of dd excitations. Further discussions about the RIXS spectra are given in Supplementary information, which includes refs. 58,59,60.

Related Articles

Emergent quantum Majorana metal from a chiral spin liquid

We propose a mechanism to explain the emergence of an intermediate gapless spin liquid phase in the antiferromagnetic Kitaev model in an externally applied magnetic field, sandwiched between the well-known gapped chiral spin liquid and the gapped partially polarized phase. We propose that, in moderate fields, π-fluxes nucleate in the ground state and trap Majorana zero modes. As these fluxes proliferate with increasing field, the Majorana zero modes overlap creating an emergent ({{mathbb{Z}}}_{2}) quantum Majorana metallic state with a “Fermi surface” at zero energy. We further show that the Majorana spectral function captures the dynamical spin and dimer correlations obtained by the infinite Projected Entangled Pair States method, thereby validating our variational approach.

Topological excitations at time vortices in periodically driven systems

We consider two-dimensional periodically driven systems of fermions with particle-hole symmetry. Such systems support non-trivial topological phases, including ones that cannot be realized in equilibrium. We show that a space-time defect in the driving Hamiltonian, dubbed a “time vortex,” can bind π Majorana modes. A time vortex is a point in space around which the phase lag of the Hamiltonian changes by a multiple of 2π. We demonstrate this behavior on a periodically driven version of Kitaev’s honeycomb spin model, where ({{mathbb{Z}}}_{2}) fluxes and time vortices can realize any combination of 0 and π Majorana modes. We show that a time vortex can be created using Clifford gates, simplifying its realization in near-term quantum simulators.

Microscopic roadmap to a Kitaev-Yao-Lee spin-orbital liquid

The exactly solvable spin-1/2 Kitaev model on a honeycomb lattice has drawn significant interest, as it offers a pathway to realizing the long-sought-after quantum spin liquid. Building upon the Kitaev model, Yao and Lee introduced another exactly solvable model on an unusual star lattice featuring non-abelian spinons. The additional pseudospin degrees of freedom in this model could provide greater stability against perturbations, making this model appealing. However, a mechanism to realize such an interaction in a standard honeycomb lattice remains unknown. Here, we provide a microscopic theory to obtain the Yao-Lee model, on a honeycomb lattice by utilizing strong spin-orbit coupling of anions edge-shared between two eg ions in the exchange processes. This mechanism leads to the desired bond-dependent interaction among spins rather than orbitals, unique to our model, implying that the orbitals fractionalize into gapless Majorana fermions and fermionic octupolar excitations emerge. Since the conventional Kugel-Khomskii interaction also appears, we examine the phase diagram, including these interactions, using classical Monte Carlo simulations and exact diagonalization techniques. Our findings reveal a broad region of disordered states that break rotational symmetry in the bond energy, suggesting intriguing behavior reminiscent of a spin-orbital liquid.

Z2 flux binding to higher-spin impurities in the Kitaev spin liquid

Stabilizing Z2 fluxes in Kitaev spin liquids (KSLs) is crucial for both characterizing candidate materials and identifying Ising anyons. In this study, we investigate the effects of spin-S magnetic impurities embedded in the spin-1/2 KSL. Utilizing exact diagonalization and density matrix renormalization group methods, we examine the impurity magnetization and ground-state flux sector with varying impurity coupling and spin size. Our findings reveal that impurity magnetization exhibits an integer/half-integer spin dependence, which aligns with analytical predictions, and a flux-sector transition from bound-flux to zero-flux occurs at low coupling strengths, independent of the impurity spin. Notably, for spin-3/2 impurities, we observe a reentrant bound-flux sector, which remains stable under magnetic fields. By considering fermionic representations of our spin Hamiltonian, we provide phenomenological explanations for the transitions. Our results suggest a novel way of binding a flux in KSLs, beyond the proposals of vacancies or Kondo impurities.

Efficient computation using spatial-photonic Ising machines with low-rank and circulant matrix constraints

Spatial-photonic Ising machines (SPIMs) have shown promise as an energy-efficient Ising machine, but currently can only solve a limited set of Ising problems. There is currently limited understanding on what experimental constraints may impact the performance of SPIM, and what computationally intensive problems can be efficiently solved by SPIM. Our results indicate that the performance of SPIMs is critically affected by the rank and precision of the coupling matrices. By developing and assessing advanced decomposition techniques, we expand the range of problems SPIMs can solve, overcoming the limitations of traditional Mattis-type matrices. Our approach accommodates a diverse array of coupling matrices, including those with inherently low ranks, applicable to complex NP-complete problems. We explore the practical benefits of the low-rank approximation in optimisation tasks, particularly in financial optimisation, to demonstrate the real-world applications of SPIMs. Finally, we evaluate the computational limitations imposed by SPIM hardware precision and suggest strategies to optimise the performance of these systems within these constraints.

Responses

Your email address will not be published. Required fields are marked *