## Unconventional Magnetic Order on the Hyperhoneycomb Kitaev Lattice in β -Li IrO : 2 3 Full Solution via Magnetic Resonant X-ray Diffraction

A. Biffin, 1 R. D. Johnson, 1 Sungkyun Choi, 1 F. Freund, 2 S. Manni, 2

A. Bombardi, 3 P. Manuel, 4 P. Gegenwart 2 and R. Coldea 1

1 Clarendon Laboratory, University of Oxford, Parks Road, Oxford OX1 3PU, U.K.

2 EP VI, Center for Electronic Correlations and Magnetism,

Augsburg University, D-86159 Augsburg, Germany

3 Diamond Light Source Ltd., Harwell Science and Innovation Campus, OX11 0DE, U.K. and

4 ISIS Facility, Rutherford Appleton Laboratory, Harwell Oxford, OX11 0QX, U.K.

The recently-synthesized iridate β -Li 2 IrO 3 has been proposed as a candidate to display novel magnetic behavior stabilized by frustration effects from bond-dependent, anisotropic interactions (Kitaev model) on a three-dimensional 'hyperhoneycomb' lattice. Here we report a combined study using neutron powder diffraction and magnetic resonant x-ray diffraction to solve the complete magnetic structure. We find a complex, incommensurate magnetic order with non-coplanar and counter-rotating Ir moments, which surprisingly shares many of its features with the related structural polytype 'stripyhoneycomb' γ -Li 2 IrO , where dominant Kitaev interactions have been 3 invoked to explain the stability of the observed magnetic structure. The similarities of behavior between those two structural polytypes, which have different global lattice topologies but the same local connectivity, is strongly suggestive that the same magnetic interactions and the same underlying mechanism governs the stability of the magnetic order in both materials, indicating that both β -and γ -Li 2 IrO 3 are strong candidates to realize dominant Kitaev interactions in a solid state material.

PACS numbers: 75.25.-j, 75.10.Jm

## I. INTRODUCTION

Materials containing ions with a strong spin-orbit interaction are attracting much attention as candidates to display novel electronic states ranging from topological insulators with protected gapless surface states to quantum spin liquids with exotic excitations. 1 Complex magnetic behavior can arise when the strong spin-orbit coupling (as found for heavy transition metal 4 d and 5 d ions) and crystal field effects stabilize spin-orbit entangled magnetic moments, which may interact via anisotropic exchange interactions where the anisotropy axis depends on the bond orientation. On certain lattices these interactions may be strongly frustrated, leading potentially to novel forms of cooperative magnetic order and/or excitations, not found for magnets of 3 d ions (where the spin-orbit coupling is much weaker and the orbital moment is in general quenched). The magnetic ground state of Ir 4+ ions in a cubic IrO 6 octahedron is a doublet with a mixed spin-orbital character, 2 J eff = 1 2, and for edge-/ sharing IrO 6 octahedra it was proposed 3,4 that the leading super-exchange is a ferromagnetic Ising coupling between the Ir magnetic moment components perpendicular to the Ir-O -Ir 2 plane. For three-fold coordinated, edge-sharing octahedra this leads to orthogonal Ir-O 2 -Ir planes for the three Ir-Ir bonds and correspondingly orthogonal x y , , z components coupled along the three Ir-Ir bonds. For a honyecomb lattice this realizes the Kitaev model, 5 where the strong frustration effects stabilize an (exactly solvable) quantum spin liquid state with novel excitations (Majorana fermions and fluxes). 6 In search for realizations of such physics, α -Na IrO 2 3 [Refs. 7-11] and α -Li 2 IrO 3 [Refs. 12,13] have been actively explored experimentally, however no clear evidence of novel Kitaev phenomena has been observed.

Generalization of the Kitaev model to threedimensional (3D) lattices have also been shown to have quantum spin liquid ground states. 14-16 Furthermore, non-trivial magnetic behavior has been predicted when (finite) additional interactions suppress the spin liquid and stabilize magnetic order. 15,17 The recently synthesized 'hyperhoneycomb' 18 β -Li 2 IrO 3 [see Fig. 1] and 'stripyhoneycomb' 19 γ -Li 2 IrO , 3 which have the local connectivity of three-fold coordinated IrO 6 octahedra with near-orthogonal Ir-O -Ir planes meeting at each Ir 2 site, are prime candidates to realize 3D Kitaev physics. The availability of single crystal samples of both of those two distinct structural polytypes offers unique opportunities to perform comparative studies between them and gain insight into the underlying physics. In recent experiments 20 on the γ -polytype we have observed an unexpectedly complex, yet highly-symmetric incommensurate magnetic structure with non-coplanar and counterrotating moments, and theoretical calculations showed that this structure could be stabilized by a spin Hamiltonian with dominant Kitaev couplings and some additional interactions. Here we extend our experimental studies to the β -polytype, where we find some striking similarities in the magnetic structure. This shows that the key features of the magnetic order are robust, independent of the changes in the global lattice topology between the two distinct structural polytypes that share the same building blocks. This is strongly suggestive that the same underlying magnetic interactions govern the stability of the magnetic order in both structural polytypes, which would be an important constraint on any theoreti-

cal model of the cooperative magnetism in this family of materials.

The structural polytype β -Li 2 IrO 3 was discovered by Takayama et. al. 18 and the crystal structure is illustrated in Fig. 1(left). It has an orthorhombic unit cell (space group Fddd , for full structural details see Appendix A) and the same mass density and same fundamental building blocks - edge-sharing LiO 6 and IrO 6 octahedra (the latter shaded in red) as in the layered polytype α -Li 2 IrO 3 (Ref. 12), where edge-sharing IrO 6 octahedra form planar honeycomb layers separated by hexagonal Li layers. In the β -structure [see Fig. 1] the Ir lattice is three-fold coordinated as in a planar honeycomb, but now the links of the lattice form a 3D network; this 3D lattice connectivity has been named a 'hyperhoneycomb'. 18 It can be thought of as being made of zig-zag chains [shaded cyan and purple in Fig. 1(right)] stacked along c and alternating in direction between the two basal plane diagonals a ± b (the planar honeycomb structure is obtained when the zig-zag iridium chains are not alternating in orientation, but are all running parallel to either the a + , or b the a -b directions).

FIG. 1: (color online) Crystal structure of β -Li 2 IrO . 3 Two unit cells are shown: (left) full structure with Li (white balls), O (black) and Ir (red) located inside IrO 6 octahedra (shaded polyhedra), (right) Ir ions are arranged in a 'hyperhoneycomb' structure. Light (cyan) and dark (purple) colors indicate zig-zag chains stacked along c and directed alternatingly along a ± b .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-4d29713b-fafb-42f4-b55c-3cbff2d2db04.jpg)

Earlier studies 18 have shown that β -Li 2 IrO 3 is an insulator. The temperature-dependence of the magnetic susceptibility parameterized in terms of a Curie-Weiss law gave an effective magnetic moment µ eff = 1.61 µ B , consistent with localized, J eff = 1 2 moments at the /

iridium sites. Both specific heat and susceptibility data showed an anomaly near T N = 38 K attributed to the onset of magnetic order. Here we confirm that the lowtemperature phase has magnetic long-range order and we provide a full experimental magnetic structure solution. Our results give uniquely detailed information about the correlations that govern the cooperative magnetism and the relevant magnetic interactions that stabilize them. We first used neutron diffraction on a powder sample to confirm the presence of magnetic order at low temperatures and to obtain candidate magnetic propagation vectors. We then performed magnetic resonant x-ray diffraction (MRXD) experiments on a 17 µ m diameter single crystal at the L 3 edge of iridium, where experiments on other iridates 2,8,21 have reported a strong enhancement of the magnetic scattering cross-section. We exploit the polarization dependence of the MRXD intensity (probed via azimuth scans) to deduce that the magnetic structure has rotating magnetic moments and we determine the plane of rotation and relative phases between all 16 iridium sites in the structural unit cell.

The paper is organized as follows: Sec. II presents the neutron powder diffraction measurements and the analysis by which we obtain candidate magnetic propagation vectors. Sec. III presents the single-crystal MRXD measurements, which observe magnetic Bragg peaks with an incommensurate propagation vector along the (100) direction. The observed diffraction pattern is analyzed in terms of magnetic basis vectors in Sec. III A and III B. The azimuth dependence of the diffraction intensity is used to determine the polarization of all magnetic basis vectors in the ground state and the relative phase between them (in Sec. III C). Finally the absolute value of the ordered magnetic moment is extracted from the neutron powder data in Sec. III D. The resulting magnetic structure is presented in Sec. III E and similarities with the magnetic structure in the γ -polytype are discussed in Sec. IV. Finally, conclusions are summarized in Sec. V. Further technical details of the analysis are presented in the Appendices: (A) crystal structure refinement at low temperature from single-crystal x-ray data, (B) decomposition of the magnetic structure in terms of its Fourier components, (C) derivation of the selection rules for magnetic scattering for the various basis vectors, (D) the azimuth dependence of the MRXD intensity enabling determination of the relative phase between basis vectors, and (E) an equivalence mapping between the magnetic structures in the β - and γ -polytypes of Li 2 IrO . 3

## II. MAGNETIC NEUTRON POWDER DIFFRACTION

Neutron powder diffraction measurements to obtain information about the magnetic propagation vector were performed using the time-of-flight diffractometer WISH at ISIS. 0.71g of powder β -Li 2 IrO 3 (synthesized as described in Appendix A) was placed in an aluminium can

of annular geometry (to minimize the strong neutron absorption by the Ir nuclei) and was mounted on the cold finger of a closed cycle refrigerator. Diffraction patterns were collected at a selection of temperatures from base (5.6 K) to paramagnetic (70 K), well above the magnetic ordering transition T N = 38 K inferred from thermodynamic measurements. 18 The obtained neutron diffraction pattern did not allow for a full structural refinement (most likely due to neutron absorption effects in the presence of iridium). We therefore performed additional x-ray measurements on a small single crystal piece extracted from the same powder batch to determine precise internal atomic positions in the unit cell at a temperature of 100 K, considered to be cold enough to be representative of the crystal structure in the low-temperature limit. The x-ray diffraction was performed on a Mo-source Oxford Diffraction Supernova x-ray diffractometer under a N 2 gas flow, which gave a temperature of 100 K at the sample position. The data confirmed the expected space group and full structural refinement (for details see Appendix A) gave atomic positions consistent with those previously reported at room temperature. 18

Fig. 2 shows the neutron powder diffraction pattern in the lowest angle bank of detectors, which covers the region of large d -spacings where magnetic diffraction is expected. In the main panel the visible diffraction peaks are of structural origin. β -Li 2 IrO 3 peak positions are denoted by the upper row of vertical marks (below the data) and the two large peaks in the data at d -spacing near 2 ˚ (labelled 'Al') are the lowest-order diffraction peaks A from the Al can and were fitted as such. The pattern also showed two peaks labelled (*), originating from a small impurity phase identified as pure Ir, and those regions in the data were excluded from the analysis. We performed the structural fit (using FullProf 22 ) to the neutron data (both paramagnetic and base temperature) with the internal atomic positions kept fixed, and only the lattice parameters left free to vary. The structural fit was consistent over all detector banks - corresponding to a range of scattering angles and resolutions - and we were also able to consistently fit the instrument parameters (which provide the conversion from time-of-flight to d -spacing) over all the detector banks tested. Fig. 2 illustrates the very good agreement in describing the structural diffraction pattern.

## A. Candidate magnetic propagation vectors

The inset in Fig. 2 highlights the region of large d -spacing, where upon cooling below T N two new Bragg peaks appear, attributed to the onset of magnetic ordering. The peak widths are comparable to those of structural Bragg peaks, implying long-range magnetic order. With only two magnetic peaks the magnetic structure cannot be solved uniquely, however the data provides strong constraints on the magnetic propagation vector, q . Typically q will lie on a point, line or plane of high sym-

FIG. 2: (color online) Neutron powder diffraction at base temperature (5.6 K, red circles) and in the paramagnetic regime (70 K, green circles) in the lowest-angle detector bank. Solid black line shows the fit to the structural and magnetic contributions as discussed in the text. Positions of structural β -Li IrO 2 3 peaks, aluminium peaks and magnetic Bragg peaks are marked below the pattern in the upper, middle and lower rows, respectively and the blue line underneath represents the difference between data and fit. Labels 'Al' indicate aluminium peaks and (*) denote two peaks from an impurity phase identified as pure Ir, removed from the fit. Inset: zoom into the large d -spacing region showing two magnetic Bragg peaks, labels indicate the magnetic peaks index (symmetrically equivalent indexes are omitted).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-363f56f8-b429-4832-89a4-c9c5428c5211.jpg)

metry in the Brillouin zone. Using the k-search tool in FullProf 22 we systematically searched all such positions for the Fddd space group in order of decreasing symmetry to find propagation vectors that could reproduce the magnetic peak positions observed in the powder neutron data. All commensurate orderings with wavevectors at the Γ, Z Y , , T and L points (where the Miller-Love convention for labelling high-symmetry Brillouin zone points has been used 23 ) were ruled out by this analysis. Therefore we considered propagation vectors along lines of high symmetry with general points denoted by Σ, ∆, Λ. Only the first two were compatible with the neutron data and the obtained solutions are listed in Table I, along with possible solutions for other high symmetry directions.

TABLE I: Candidate magnetic propagation vectors along high symmetry directions in the Brillouin zone compatible with the magnetic neutron powder diffraction pattern. The ellipses denotes the fact that multiple solutions exist for a general propagation vector in the ab -plane.

| Position   | q fit            |
|------------|------------------|
| Σ          | (0.57,0,0)       |
| ∆          | (0,0.81,0)       |
| J          | (0.27,0, 1 / 2 ) |
| E          | (0,0.39, 1 / 2 ) |
| M          | (0.43,0.53,0)    |
|            | . . .            |

## III. MAGNETIC RESONANT X-RAY DIFFRACTION

To determine which of the candidate magnetic propagation wavevectors identified by the neutron data actually occurs, we have performed a magnetic resonant x-ray diffraction experiment in reflection geometry at the L 3 edge of Ir using the I16 beamline at Diamond. We used a single crystal of β -Li 2 IrO 3 (characterized via x-ray diffraction as described in Appendix A) with an orthorhombic morphology with its maximum dimension less than 17 µ m. The sample was placed on a Si (111) disk with the reciprocal lattice vector (-6,5,7) being approximately surface normal and it was cooled using a closedcycle refrigerator with Be domes. The x-ray energy was tuned to the L 3 edge of Iridium (11.215 keV) and on cooling below T N new Bragg peaks appeared at satellite positions of reciprocal lattice points at ( h k l , , ) ± q , with h k l , , integers and an incommensurate ordering wavevector q =(0.57(1), 0, 0). Scans through such a peak are shown in Fig. 3a-c). The scans emphasize the incommensurate position along h and centering at integer k and l values. In the language of Table I this is then a propagation vector along the Σ line of symmetry, and the magnitude found from the x-ray data is consistent with the positions of peaks in the neutron data (first line in Table I). Lorentzian squared fits to the peak shape profiles (dashed lines in Figs. 3a-c),e) gave widths comparable to those measured for the nearby structural Bragg peak (-6,6,12), indicating that the magnetic order is longrange in all three directions. To estimate a lower bound on the correlation lengths we use the inverse of the peak half-widths, which gives values in excess of 100, 200 and 300 ˚ along the A a , b and c axes, respectively (more than ≈ 20 unit cells in each direction), and in reality the correlation lengths will be much larger as the above estimate did not include the peak broadening due to the finite instrumental angular resolution.

Several independent tests can be performed to confirm the magnetic origin of the observed diffraction peaks. Firstly, in Fig. 3(d) the integrated intensity of (-6,6,12)q as a function of temperature (filled circles) shows a typical order-parameter behavior. A power-law fit (solid line) gives an onset temperature of T fit = 36.7 K, which is very close to the transition temperature deduced from specific heat measurements. 18 We attribute the apparent small offset of ∆ T glyph[similarequal] 1.3 K to the fact that the sensor where the temperature was measured in the x-ray experiments was at the bottom of the sample support and the actual sample temperature may have been higher by ∆ T due to local heating effects from the very intense synchrotron beam, a common occurrence in resonant x-ray experiments.

Further evidence of the magnetic origin of the incommensurate diffraction peaks can be obtained via polarization analysis of the scattered beam. For a scattering experiment with incident beam polarization ˆ glyph[epsilon1] normal to the ( k k , ′ ) scattering plane ( σ -polarized), mag-

FIG. 3: (color online) Magnetic Bragg peak at (-6,6,12) -q ; (a) Scan shows that the peak occurs at an incommensurate position along h , and is well-centered at integer (b) k , and (c) l . (d) Temperature-dependence of the integrated peak intensity (solid line is a power-law fit).(e) Scans with a polarizer in the scattered beam show that the peak is only present in the rotated polarization channel ( σ π -′ ) (light blue, filled) as expected for magnetic scattering. (f) Intensity as a function of the x-ray energy (filled points) shows that the signal is resonant on the edge of the fluorescence signal onset (dashed line), as expected for magnetic scattering.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-5873df71-8745-4f3f-8ae6-dbe80d555292.jpg)

netic scattering is expected to rotate the polarization by 90 ◦ compared to the incident beam direction (scattering in the σ π -′ channel, see diagram in Fig. 4a) (inset)), whereas conventional (charge) scattering would leave the polarization direction unchanged (scattering in the σ σ -′ channel). 25 In the experiments polarization analysis was achieved via the 3 rd harmonic reflection from the (1,1,1) planes of an Au crystal, the d -spacing of 0.785 ˚ being A ideally suited for polarization analysis at the L 3 edge of Ir (giving perfect filtering at 11.16 keV, see Ref. 24). Fig. 3(e) shows an h -scan across the incommensurate peak performed with polarizer parallel (open symbols, σ σ -′ ) and perpendicular (filled symbols, σ π -′ ) to the original incident beam polarization. The signal appears only in the rotated polarization channel, as expected for magnetic scattering. Note that the higher background in the σ σ -′ channel is due to incoherent x-ray scattering, which is largely filtered out by the polarization analysis in the σ π -′ channel revealing a pure magnetic signal.

A third test is the dependence on the x-ray energy as the magnetic scattering is only expected at resonance. Indeed, the incommensurate peak intensity has a characteristic resonant response as shown in Fig. 3(f) (filled

blue symbols) with a maximum at the resonance energy of 11.215 keV in agreement with previous resonant studies on other iridates, 8,21 and coincides with the edge of the fluorescence signal from the sample (black, dotted line).

## A. Magnetic Basis Vectors

The observation of incommensurate magnetic Bragg peaks indicates a magnetic structure with rotating moments or a spin-density wave, both with a wavelength that is incommensurate with the underlying crystal lattice periodicity. In order to distinguish between those two scenarios, and to determine fully the orientations of the magnetic moments at all magnetic sites in the unit cell, we consider below symmetry-allowed basis vectors for magnetic structures compatible with the observed propagation vector q . This provides a natural framework to directly link diffraction data with a magnetic structure model and allows one to develop a systematic strategy in experiments to determine the magnetic structure completely.

The magnetic ions, Ir 4+ , occupy a single crystallographic site with ions at (1 / 8, 1 / 8, z ) [ z = 0.70845(7) glyph[similarequal] 5/8+1/12] and symmetry equivalent positions, giving four Ir ions in the primitive unit cell (labelled 1-4 in Fig. 6 and with positions listed explicitly in Table IV in Appendix B). Taking into account that at each site the magnetic moment could have components along the x , y and z -axes (along the orthorhombic a , b , c -axes), gives a 12-component representation of the magnetic structure. Its irreducible representations and associated magnetic basis vectors (obtained using the BasIrreps tool, part of the FullProf suite 22 ) are listed in Table II.

TABLE II: Irreducible representations and basis vectors for a magnetic structure with propagation vector q =( , 0, 0). q

| Irreducible Representation   | Basis Vectors   |
|------------------------------|-----------------|
| Γ 1                          | F x , G y , A z |
| Γ 2                          | C x , A y , G z |
| Γ 3                          | G x , F y , C z |
| Γ 4                          | A x , C y , F z |

The basis vectors contain symmetry-imposed relations between the Fourier components of the magnetic structure between the four sites in the primitive cell. The types of basis vectors that can occur are:

$$F = \begin{bmatrix} 1 \\ 1 \\ \delta \\ \delta \end{bmatrix}, C = \begin{bmatrix} 1 \\ 1 \\ - \delta \\ - \delta \end{bmatrix}, A = \begin{bmatrix} 1 \\ - 1 \\ - \delta \\ \delta \end{bmatrix}, G = \begin{bmatrix} 1 \\ - 1 \\ \delta \\ - \delta \end{bmatrix}, \quad & \text{where} \\ ( 1 ) \quad & \text{where} \\ \text{where for each vector the four values are relative phase} \quad & ( h, k,$$

where for each vector the four values are relative phase factors between the Fourier components at the four sites.

δ = e -i q · ( r 3 -r 1 ) = e -iπq/ 2 is a displacement phase factor that takes into account the fact that sites 3 and 4 are displaced relative to sites 1 and 2 by a/ 4 in the direction of the propagation vector q . In an F basis vector the Fourier components at the 4 sites are related by M q ,1 = M q ,2 = δ -1 M q ,3 = δ -1 M q ,4 , i.e. there are no additional phase factors between the four sites apart from the natural displacement phase factor, and in the limit q → 0 ( δ → 1) one recovers ferromagnetic order. For the basis vectors C , A and G two out of the four sites have a change in sign in addition to the normal displacement phase factor, and in the limit q → 0 one finds an antiferromagnetic alignment between the four sites. To determine the full magnetic structure, we first identify which basis vectors are present, determine their polarization (i.e. along x , y or z ) and then find the relative phase between them, as follows.

## B. Selection Rules for Magnetic Scattering

Each magnetic basis vector contains a strict, symmetry imposed phase relation between moments on all Ir sites in the primitive unit cell. As such, one can derive selection rules for non-zero magnetic diffraction intensity, i.e., each basis vector will only contribute to magnetic satellite reflections of certain reciprocal lattice points and will have zero structure factor for others. It follows that simply the presence or absence of magnetic Bragg peaks at certain positions can already identify which basis vectors are present.

Explicitly, the structure factor for a magnetic superlattice reflection at Q = ( h k l , , ) ± q is

$$\mathcal { F } ( Q ) = \mathcal { F } ( ( h, k, l ) \pm q ) = f _ { F } \, \sum _ { n } M _ { \pm q, n } e ^ { i Q \cdot r _ { n } }, \quad ( 2 )$$

where the pre-factor f F = 1+ e iπ h ( + ) k + e iπ k ( + ) l + e iπ l ( + h ) is due to the F -centering of the orthorhombic structural unit cell. The sum extends over all sites in the primitive cell ( n =1 -4) and M ± q , n are the Fourier components of the magnetic moments at site n with position in the unit cell r n . For an F -basis vector we obtain the structure factor

$$\text{ions} \\ \text{ruc-} \\ \text{The} \quad \mathcal { F } ^ { F } ( Q ) & = \begin{cases} \ 1 6 M _ { \pm q, 1 } \cos \frac { \pi l } { 6 } e ^ { i \xi _ { \pm } }, \ h, k, l \, \text{all even and} \\ h + k + l = 4 p, \ p \, \text{integer} \\ \ 1 6 M _ { \pm q, 1 } \cos \frac { \pi l } { 6 } \cos \frac { \pi ( h + l - k ) } { 4 } e ^ { i \chi _ { \pm } }, \\ h, k, l \, \text{all odd} \\ 0, \quad \text{otherwise,} \end{cases} \\ \Big \rfloor \quad & \text{where} \, \xi _ { \pm } = \pi ( h + k - 3 l \, \pm q ) / 4 \, \text{and} \, \chi _ { \pm } = \pi ( 2 h - 2 l \, \pm }$$

where ξ ± = π h ( + k -3 l ± q / ) 4 and χ ± = π (2 h -2 l ± q / ) 4 and to obtain a closed-form analytic expression we used the 'ideal' iridium position z =5/8+1/12. From (3) we conclude that an F basis vector will contribute to the intensity of ± q satellites of reciprocal lattice points ( h k l , , ) with the selection rule h k l , , all odd, or all even with h + k + l = 4 , p p integer. This is the same selection

rule as for structural Bragg scattering from the iridium ion sublattice at the ( h k l , , ) position, as expected since in the limit q → 0 the F -type magnetic basis vectors recover ferromagnetic order, and this has the same selection rule as structural scattering. Structure factor expressions for the other basis vectors are given in Appendix C. The resulting selection rule for A is same as for F , but with no contribution from A for l = 6 m m , integer. The selection rule for C and G is h k l , , all odd, or all even and h + k + l = 4 p + 2, p integer with no contribution from G if l = 6 m m , integer.

The experimental geometry constrained the search for magnetic satellite reflections to regions of reciprocal space in the vicinity of the surface normal reflection, ≈ (-6,5,7). Using (3) and similar expressions for C G , and A (see Appendix C) we ascertained 3 positions of interest for identifying which basis vectors contribute to the magnetic structure, namely (-6,6,12) ± q (pure F ), (-8,6,12) ± q (pure C ) and (-5,5,9) ± q (mixed A and G ) and found magnetic peaks at all these positions. Therefore we deduce that the structure contains F , C and one (or both) of A and G . To determine the direction of the magnetic moment components for each of those basis vectors we make use of the polarization dependence of the magnetic x-ray diffraction intensity, as follows.

## C. Azimuth Scans to Determine the Moment's Direction and Relative Phases Between Basis Vectors

In the dipolar approximation the magnetic x-ray diffraction intensity at resonance is proportional to

$$L ( \theta ) \mathcal { A } \left | ( \hat { \epsilon } ^ { \prime } \times \hat { \epsilon } ) \cdot \mathcal { F } ( Q ) \right | ^ { 2 },$$

where L θ ( ) is the Lorentz factor at the scattering angle 2 , θ A is an absorption correction dependent upon the experimental geometry, F ( Q ) is the magnetic structure factor vector given in (2), and ˆ glyph[epsilon1] ′ and ˆ glyph[epsilon1] are unit vectors along the polarization of the electric field component of the scattered and incident x-ray beams, respectively. 25 For a σ -polarized incident beam magnetic resonant scattering occurs only in the σ π -′ channel [see diagram in Fig. 4a) inset], meaning that the product of the electric field polarization vectors is along the scattered beam direction, i.e. ˆ glyph[epsilon1] ′ × ˆ glyph[epsilon1] = ˆ k ′ , so only the component of the structure factor vector along the scattered beam direction, F ‖ = F · ˆ k ′ , contributes to the magnetic intensity. By keeping the instrument in the scattering condition and rotating the sample around the scattering wavevector, Q = k ′ -k , the projection F ‖ of the structure factor vector onto the (fixed) direction ˆ k ′ varies depending on the azimuth angle Ψ, with maximum magnetic intensity when the magnetic moments that give rise to the scattering make the smallest angle with ˆ k ′ and zero intensity when they are perpendicular.

FIG. 4: (color online) Intensity as a function of azimuth for three magnetic Bragg peaks a) C y , b) F z and c) A x . Top diagram illustrates the scattering geometry. Data points (filled circles) are integrated peak intensities (in sample rocking curve scans) corrected for absorption and Lorentz factor. Thick (red) lines show fits that include all contributions to the magnetic scattering intensity 26 for the magnetic structure model Γ 4 depicted in Fig. 6. Orange, green and blue lines denote other models (Γ 1,2,3 , respectively), which can be easily ruled out. The azimuth origin Ψ = 0 corresponds to the case when the (010) axis is in the scattering plane making the smallest angle with the incident beam direction ˆ k .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-92795362-a307-4429-b6f8-c085a9be0d8f.jpg)

Azimuth scans to test the polarization dependence are shown in Figs. 4a-c) where the data points are obtained by integrating the magnetic peak intensity in sample rocking curve scans for each value of the azimuth, Ψ. As discussed earlier, the selection rules identify the signal in panel a) at (-8,6,12)+ q to be of C -character and the observed maximum intensity at Ψ = 0 uniquely iden-

tifies the polarization to be along the y -axis (red thick line) as at Ψ = 0 the y -axis is in the scattering plane and makes the smallest angle with the ˆ k ′ direction and makes a larger angle with ˆ k ′ when the azimuth is displaced away from 0 in both directions; green and blue curves show the expected intensity dependence for other polarizations that can be clearly ruled out. Similarly, in panel b) the peak at (-6,6,12) -q is identified as being of F character and polarized along z (thick red line). Panel c) shows data at (-5,5,9) -q , clearly of x -polarization (thick red line). The structure factor is compatible with both an A or G character and we will show later that only a pure A character can describe this and other azimuth scans quantitatively, so we identify this peak as coming from an A x basis vector.

FIG. 5: (color online) Intensity as a function of azimuth for 'mixed' magnetic Bragg peaks (black circles) corrected for absorption and Lorentz factor, compared to calculations (red thick line) for the magnetic structure model in Fig. 6, as well as calculations for other possible models (a) (-6,6,8)+ q is a mixed F : A peak with contributions in the ratio 1 : 3, (b) (-5,5,7) -q is a mixed F : C : A peak with contributions in the ratio 3 : 3 : 1. Note that the data clearly distinguishes between different relative phases of basis vectors (cyan and purple).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-b169743c-ee75-4092-a7e2-5d6ce0be8acd.jpg)

The next step in solving the structure is to determine the relative phase between the A x , C y and F z basis vectors. A phase of 0 or π between basis vectors with orthogonal components, say along A x and F z , means that the moments vary sinusoidally along a direction in the xz plane, whereas a phase offset of π/ 2 means that the moments rotate in the xz plane. The relative phase between two basis vectors can be determined from an azimuth scan at a magnetic Bragg peak position where both basis vectors contribute as the intensity is the sum of the intensities due to the two separate basis vectors plus an additional cross-term that is sensitive to the phase offset (for more details see Appendix C).

Fig. 5(a) shows the azimuth scan at (-6,6,8)+ q where structure-factor calculations show that both F and A basis vectors contribute in the ratio 1:3. The intensity dependence uniquely identifies it as originating from the particular basis vector combination F z + iA x , i.e. the x components are π/ 2 out-of-phase with respect to the z -components, meaning that the magnetic moments rotate in the xz plane (in a particular sense on each site, note that the opposite sense of rotations encoded in the basis vector F z -iA x is not consistent with the data cyan line). The data clearly rules out the case when the moments are not rotating in the xz plane, but are varying sinusoidally along a line in the xz plane ( F z ± A x , purple line), a model with no A x component is also ruled out (not shown). This analysis shows that an A x basis vector must exist in the ground state and that it is π/ 2 out-of-phase with respect to F z . Similarly, to determine the phase offset between the F z and C y basis vectors we measured the azimuth scan at the magnetic reflection (5,5,7) -q of mixed F , C and A contributions in the ratio 3:3:1, shown in Fig. 5(b), and only a model with the basis vector contribution F z + iC y (thick red line) matches the data, with other combinations (blue and green curves) clearly ruled out (note that those calculations also include the iA x contribution as already determined).

Having established all the basis vectors that are present, their polarization, and their relative phases, the relative magnitudes are obtained by quantitatively fitting the azimuth intensity dependencies. Three free parameters were fit; the relative ratios M /M x z and M /M y z , and a single overall intensity scale factor for all the data. The simultaneous fit of this magnetic structure model to the data in Figs. 4a-c) and Fig. 5a-b) (shown by the thick red solid lines) gave the moment amplitudes in the relative ratio M x : M y : M z = 0.45(1) : 0.65(1) : 1. Note that the basis vector combination ( iA x , iC y , F z ) with the above moment amplitudes quantitatively explains the mixed F z , A x peak in Fig. 5a), and the same magnitude A x basis vector also quantitatively accounts for the azimuth dependence in Fig. 4c), further confirming the identification of the signal in this magnetic Bragg peak as coming from a pure A x basis vector (if a G x basis vector is also present and contributing to the signal in Fig. 4c), its magnitude is very small, below the accuracy of the present experiments). Given this we conclude that the magnetic structure is described by the basis vector combination

( iA x , iC y , F z ), which corresponds to a single irreducible representation, Γ 4 in Table II.

## D. Absolute Value of the Ordered Moment

The only remaining parameter still to be determined is the absolute magnitude of the magnetic moments. This is difficult to extract reliably from the magnetic resonant x-ray diffraction data as it requires accurate determination of scale factors between the magnetic and structural peaks (the latter being of order 10 5 -10 6 more intense). However, in neutron scattering experiments the magnetic and structural peaks are of comparable magnitudes and one can obtain a reliable determination of the relative scale factor. To this aim we return to the powder neutron data in Fig. 2 and simultaneously fit three contributions to the data: structural peaks of the sample calculated using the crystal structure with fixed internal parameters (as deduced from single crystal x-ray diffraction), aluminium structural peaks (from the sample can), and magnetic peaks calculated using the full magnetic structure model deduced from the resonant x-ray experiments. The fit is plotted as a solid black line in Fig. 2 and shows excellent agreement with the data. The resulting magnetic moment is of magnitude 0.47(1) µ B when aligned along the c -axis, and slightly reduced by ∼ 20% when in the ab plane. The propagation vector was also fitted and we find q =(0.5768(3), 0, 0), consistent with the x-ray measurements.

## E. Magnetic Structure

The obtained magnetic structure projected onto the ac plane is shown in Fig. 6a). The magnetic moments are counter-rotating between all nearest-neighbor sites and the plane of rotation alternates between the two sites of every vertical ( c -axis) bond as illustrated by the pattern of light and dark shaded elliptical envelopes. The planes of rotation are obtained from the ac plane by rotation around the c -axis by an angle ± φ , with φ = tan -1 M y M x = 55(1) ◦ . Each zig-zag chain has the top and bottom sites counter-rotating in one and the same plane, which then alternates between consecutive zig-zag chains verticallylinked along the c -axis. The zig-zag chains are directed alternatingly along the a ± b directions and are made up of iridium sites of type (2,4) and (1,3), respectively, with vertical ( c -axis) bonds coupling sites of types 1-2 and 34. The magnetic order pattern is such that for every vertical bond the spin components along the y -axis are ferromagnetically-aligned (see Fig. 7), in accordance with the basis vector C y , which has equal Fourier components at sites 1-2 and 3-4, respectively [see eq. (1)].

We note that the magnetic moment rotation at each site defines an elliptical envelope, distorted from circular with the moment slightly smaller when in the ab plane compared to when along the c -axis, √ M 2 x + M /M 2 y z = 0.80(1). This effect may be due to a larger g -factor along the c -axis compared to the ˆ a cos φ ± ˆ b sin φ directions in the ab plane. We note that single-crystal susceptibility data 19 in the related polytype γ -Li 2 IrO 3 do provide evidence for the presence of g -tensor anisotropy for the iridium moments (in the high-temperature limit the susceptibility, χ , is expected to be proportional to the squared g -factor along the applied field direction, and experiments observe χ c larger than both χ a and χ b , implying an anisotropic g -tensor). The very similar local environment around the iridium sites in β -Li 2 IrO 3 suggests that an anisotropic g -tensor is also likely here. Another possibility is that zero-point quantum fluctuations may prefer a non-fixed length ordered moment between sites, with a larger ordered moment along the c -axis, which is the only direction that is common to all planes of rotation. Evidence for the presence of zero-point quantum fluctuations in the ground state is provided by the fact that the absolute magnitude of the ordered moment (0.47(1) µ B ) is significantly reduced from what is believed to be the available full-moment value (estimated at gµ B J eff glyph[similarequal] 1 µ B assuming g glyph[similarequal] 2), so a structure with small modulations on an already significantly reduced ordered moment could be compatible with the experimental results.

## IV. DISCUSSION

It is interesting to note that the obtained magnetic structure has striking similarities with the magnetic structure in the related polytype γ -Li 2 IrO 3 [Ref. 20], with which β -Li 2 IrO 3 shares the same size orthorhombic unit cell ( a × b × c ) and also many key structural features. In both polytypes the iridium lattice is locally three-fold coordinated and is made up of vertically-linked zig-zag chains that alternate in orientation between the a ± b direction. The only difference between them is that in the hyperhoneycomb β -polytype this alternation occurs between consecutive zig-zag chains, whereas in the stripyhoneycomb γ -polytype the alternation occurs between pairs of parallel zig-zag chains (which form a honeycomb row). The magnetic structures in both polytypes are incommensurate with the same propagation vectors within experimental error, q =(0.57(1), 0, 0), the moments are counter-rotating between the sites of every nearest-neighbor bond, and the plane of rotation alternates between two orientations tilted away (by an angle φ ) from the ac plane between consecutive zig-zag iridium chains vertically-linked along the c -axis. The φ angle was found to be somewhat smaller in γ -Li 2 IrO , 3 but apart from this difference the magnetic structure in γ -Li 2 IrO 3 can be regarded as 'equivalent' to that in β -Li 2 IrO 3 but in a different lattice setting (for a formal mapping see Appendix E). Those similarities are strongly suggestive that the defining features of the magnetic structure, namely

FIG. 6: (color online) a) Projection of the magnetic structure on the ac plane showing the sites 1-4 of the primitive cell. Left curly arrows indicate counter-rotation of moments between consecutive sites along c . In unit cell 2 light and dark shaded elliptical envelopes indicate an alternating tilt of the plane of moments' rotation away from the ac face. In unit cell 3 the color of bonds shows the anisotropy axis of the exchange in a Kitaev model 15 (ferromagnetic Ising exchange for each bond, but with a different Ising axis x y , , z for the blue/green/red bonds). Right-hand labels (b)-(e) indicate where slices through the magnetic structure are taken at different heights in the unit cell and projected onto the ab plane to illustrate the direction of the zig-zag chains and the alternating tilt of the plane of rotation away from the ac plane by ± φ between adjacent zig-zag chains stacked along c .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-64d1e475-e6c9-4634-9b08-f02eff05015b.jpg)

FIG. 7: (Color online) Projection of the iridium lattice on the ac plane showing the ordering of the magnetic moment components along the b -axis only. These are shown rotated (for ease of visualization) from the b to the a -axis and indicated by horizontal red arrows (length of arrow indicates magnitude of M y ( r ) at each site. Note the ferromagnetic alignment between the two sites of each vertical ( c -axis) bond.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-9aa550bf-ef86-4bc9-b30b-5374c0e6753e.jpg)

non-coplanarity, the counter-rotation on every nearestneighbor bond, and the direction and even the magnitude of the incommensurate ordering wavevector, are determined by the same type of short-range magnetic interactions in both polytypes, and those features appear to be robust against changes in the global lattice connectivity (i.e. how zig-zag chains are vertically-linked in a long-range pattern) which leave the local connectivity of the zig-zag chains unchanged. In the γ -polytype it was shown 20 that the key features of the magnetic order can be stabilized by a spin Hamiltonian with dominant Kitaev couplings with some additional interactions; the observation of similar magnetic structures in the β -and γ -polytypes suggests the same underlying magnetic interactions occur in both cases. The relevant parent Kitaev model for the β -polytype has been considered in Refs. 15,16 and is illustrated in Fig. 6a) (unit cell 3) where the color of the bonds indicates the anisotropy axis of the exchange, i.e. a ferromagnetic Ising coupling between the components normal to the Ir-O -Ir plane of 2 each Ir-Ir bond, those directions define the cubic axes x , y and z related to the crystallographic axes by ˆ x = ( ˆ a + ) ˆ c / √ 2, ˆ y = ( ˆ a -ˆ c ) / √ 2 and ˆ z = ˆ b (where we have assumed the parent 'idealized' lattice with cubic IrO 6 octahedra and a b c : : =1: √ 2:3). We note that the magnetic interactions may also be very similar in magnitude between the β and γ -polytypes, as both materials show very similar values for the magnetic ordering temperature.

## V. CONCLUSIONS

We have reported a combined study using magnetic neutron powder diffraction and single-crystal magnetic resonant x-ray diffraction experiments at the L 3 edge of Ir to explore the magnetic structure of β -Li 2 IrO , which 3 is a unique realization of a hyperhoneycomb iridium lattice, a candidate to display 3D Kitaev physics of strongly-frustrated, bond-anisotropic interactions. Using the azimuth dependence of the magnetic Bragg peak intensities combined with magnetic symmetry analysis we have obtained a complete solution for the magnetic structure for all 16 iridium sites in the structural unit cell. We find an incommensurate, non-coplanar magnetic structure with moments counter-rotating on

- 1 for a review see W. Witczak-Krempa, G. Chen, Y. B. Kim, L. Balents, Annual Review of Condensed Matter Physics 5 , 57 (2014).
- 2 B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, H. Takagi, and T. Arima, Science 323 , 1329 (2009).
- 3 G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 102 , 017205 (2009).
- 4 J. Chaloupka, G. Jackeli and G. Khaliullin, Phys. Rev. Lett. 105 , 027204 (2010).
- 5 A. Kitaev, Ann. Phys. (N.Y.) 321 , 2 (2006).
- 6 J. Knolle, D. L. Kovrizhin, J. T. Chalker, and R. Moessner, Phys. Rev. Lett. 112 , 207203 (2014).
- 7 Y. Singh and P. Gegenwart, Phys. Rev. B 82 , 064412 (2010).
- 8 X. Liu, T. Berlijn, W.-G. Yin, W. Ku, A. M. Tsvelik, Young-June Kim, H. Gretarsson, Yogesh Singh, P. Gegenwart, and J. P. Hill, Phys. Rev. B 83 , 220403(R) (2011).
- 9 S. K. Choi, R. Coldea, A. N. Kolmogorov, T. Lancaster, I. I. Mazin, S. J. Blundell, P. G. Radaelli, Yogesh Singh, P. Gegenwart, K. R. Choi, S.-W. Cheong, P. J. Baker, C. Stock, and J. Taylor, Phys. Rev. Lett. 108 , 127204 (2012).
- 10 F. Ye, S. Chi, H. Cao, B. C. Chakoumakos, J. A. Fernandez-Baca, R. Custelcean, T. F. Qi, O. B. Korneta, and G. Cao, Phys. Rev. B 85 , 180403(R) (2012).
- 11 H. Gretarsson, J. P. Clancy, Y. Singh, P. Gegenwart, J. P. Hill, J. Kim, M. H. Upton, A. H. Said, D. Casa, T. Gog, and Y.-J. Kim, Phys. Rev. B 87 , 220407(R) (2013).
- 12 M. J. O'Malley, H. Verweij and P.M. Woodward, J. Solid State Chem. 181 , 1803 (2008).
- 13 Y. Singh, S. Manni, J. Reuther, T. Berlijn, R. Thomale, W. Ku, S. Trebst, and P. Gegenwart, Phys. Rev. Lett. 108 , 127203 (2012).
- 14 S. Mandal and N. Surendran, Phys. Rev. B 79 , 024426 (2009).
- 15 E. K.-H. Lee, R. Schaffer, S. Bhattacharjee, and Y. B. Kim, Phys. Rev. B 89 , 045117 (2014).
- 16 I. Kimchi, J. G. Analytis, A. Vishwanath, arXiv:1309.1171 (2013).
- 17 E. K.-H. Lee and Y. B. Kim, arXiv:1407.4125 (2014).
- 18 T. Takayama, A. Kato, R. Dinnebier, J. Nuss, and H. Takagi, arXiv:1403.3296 (2014).
- 19 K. A. Modic, T. E. Smidt, I. Kimchi, N. P. Breznay, A. Biffin, S. Choi, R. D. Johnson, R. Coldea, P. Watkins-Curry,

every nearest-neighbor bond. The magnetic structure shows striking similarities to the magnetic structure of the related structural polytype γ -Li 2 IrO , 3 suggesting that the same underlying Hamiltonian (with dominant Kitaev interactions) stabilizes the defining features of the magnetic structure in both structural polytypes.

## VI. ACKNOWLEDGEMENTS

Work at Oxford was supported by EPSRC (UK) and at Diamond and ISIS by STFC (UK). Work at Augsburg was supported by the Helmholtz Virtual Institute 521 ('New states of matter and their excitations').

- G. T. McCandless, J. Y. Chan, F. Gandara, Z. Islam, A. Vishwanath, A. Shekhter, R. D. McDonald, and J. G. Analytis, Nature Comm. 5 , 4203 (2014).
- 20 A. Biffin, R. D. Johnson, I. Kimchi, R. Morris, A. Bombardi, J. G. Analytis, A. Vishwanath, and R. Coldea, arXiv:1407.3954 (2014).
- 21 S. Boseggia, R. Springell, H. C. Walker, H. M. Rønnow, Ch. R¨egg, H. Okabe, M. Isobe, R. S. Perry, S. P. Collins, u and D. F. McMorrow, Phys. Rev. Lett. 110 , 117207 (2013).
- 22 J. Rodriguez-Carvajal, Physica B. 192 , 55 (1993).
- 23 S.C. Miller and W.F. Love, Tables of Irreducible Representations of Space Groups and Co-Representations of Magnetic Space Groups (Pruett, Boulder, 1967)
- 24 Resonant Elastic X-ray Scattering , Aussois, France (2011) 25 J. P. Hill and D. F. McMorrow, Acta Cryst. A 52 , 236, (1996).
- 26 http://forge.ill.eu/svn/magnetix

## Appendix A: Sample synthesis, characterization and low-temperature structural refinement

The synthesis route to obtain the samples used in the neutron and x-ray experiments was as follows. First, a powder sample of the layered honeycomb phase α -Li IrO 2 3 was prepared as described in Ref. 13. Then repetitive annealing at 1100 C transformed the sample ◦ into hyperhoneycomb β -Li 2 IrO , 3 confirmed via powder x-ray diffraction. Bulk magnetic susceptibility and heat capacity measurements on the powder sample indicated evidence for magnetic ordering at 38 K, in full agreement with Ref. 18. The polycrystalline powder contained small shiny crystallites. Several such crystallites were extracted from the powder and their diffraction pattern measured using a Mo source single-crystal SuperNova xray diffractometer. Representative diffraction patterns in the ( hk 0) and ( h l 1 ) planes obtained from the sample used in the magnetic resonant x-ray diffraction experiments (diameter 17 µ m) are shown in Fig. 8. The patterns observe sharp Bragg peaks, indicating a highquality single crystal, with selection rules consistent with

the Fddd space group expected for β -Li 2 IrO . 3 The re-

FIG. 8: (color online) X-ray diffraction pattern at 300 K in the ( hk 0) and ( h l 1 ) planes from the crystal used in the resonant x-ray diffraction experiments.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-9b6de051-c252-424a-80df-deebb5b1d258.jpg)

sults of the full structural refinement of a complete xray diffraction data set collected at 100 K from another crystal of comparable dimensions from the same batch are given in Table III where we have used the same convention for the origin of the unit cell as in Ref. 18. In the table A.D.P. stands for Atomic Displacement Parameters, which were assumed to be isotropic. The reliability factors of the Rietveld refinement were R F ( 2 ) = 4.50%, wR F ( 2 ) = 5.24%, R F ( ) = 2.42% and χ 2 = 0.318.

TABLE III: β -Li 2 IrO 3 structural parameters at 100 K

## Cell parameters

Space group: Fddd (#70, origin choice 2) Z = 16

a b c , , ( ˚ ): A 5.8903(2) 8.4261(3) 17.7924(7)

$$\alpha, \beta, \gamma \colon \quad 9 0 ^ { \circ } \quad \ 9 0 ^ { \circ } \quad \ 9 0 ^ { \circ }$$

Volume ( ˚ A ): 883.08(6) 3

Atomic fractional coordinates and isotropic A.D.P.'s

| Atom   | x        | y          | z          | U iso (A ˚ 2 )   |
|--------|----------|------------|------------|------------------|
| Ir     | 0.125    | 0.125      | 0.70845(7) | 0.0025(3)        |
| Li1    | 0.125    | 0.125      | 0.04167    | 0.00633          |
| Li2    | 0.125    | 0.125      | 0.875      | 0.00633          |
| O1     | 0.855(7) | 0.125      | 0.125      | 0.002(5)         |
| O2     | 0.621(8) | 0.3669(19) | 0.0384(7)  | 0.002(3)         |

## Data collection

# measured reflections: 3770

Data reduction R int :

6.99%

(Criterion for observed reflections: I &gt; 2.0 σ I ( ))

# observed independent reflections: 298

# fitted parameters: 9

## Appendix B: Magnetic Structure and Fourier Decomposition

This section lists the coordinates of the four iridium ions in the primitive cell and also gives explicitly the expressions for the magnetic moments at all sites in the unit cell and their Fourier decomposition. The iridium ions occupy a single crystallographic site with generating position (1 / 8, 1 / 8, z ) with z = 0.70445(7) and the four sites per primitive cell (labelled 1-4 in Fig. 6a) (left panel) have the coordinates listed in Table. IV.

TABLE IV: Fractional atomic coordinates of the iridium sites in the primitive cell and corresponding magnetic basis vector components in the determined magnetic structure.

|   Site | Coordinates                | v x   | v y   | v z   |
|--------|----------------------------|-------|-------|-------|
|      1 | (0.125, 0.125, z )         | +     | +     | +     |
|      2 | (0.125, 0.625, 3 / 4 - z ) | -     | +     | +     |
|      3 | (0.375, 0.375, 1 - z )     | -     | -     | +     |
|      4 | (0.375, 0.875, 1 / 4 + z ) | +     | -     | +     |

The determined magnetic structure is described by the basis vector combination ( iA x , iC y , F z ) with magnitudes M x , M y and M z . The magnetic moment at position r belonging to site index n = 1 -4 is given by

$$M _ { n } ( r ) \ & = \ \hat { x } M _ { x } v _ { x } ( n ) \sin q \cdot r + \hat { y } M _ { y } v _ { y } ( n ) \sin q \cdot r \\ & + \ \hat { z } M _ { z } v _ { z } ( n ) \cos q \cdot r,$$

where ˆ x y , ˆ , ˆ z are unit vectors along the orthorhombic a , b , c axes, respectively. The pre-factors v x y z , , are obtained from the basis vectors along the corresponding axes, but without the displacement phase factors e -i q · ( r n -r 1 ) , i.e. v x = [+ --+] (from A x ), v y = [+ + --] (from C y ) and v z = [+ + ++] (from F z ). Eq. (B1) describes all iridium sites, including those related by F -centering translations, where r is the actual position of the ion and n is the site index at the equivalent position (1 -4) in the primitive unit cell.

The magnetic moments are expressed in terms of their Fourier components as M r n ( ) = ∑ k = ± q M k , n e -i k r · , where

$$\text{where} \\ M _ { q, n } \, & = \, \left \{ i \left [ \hat { x } \frac { M _ { x } } { 2 } v _ { x } ( n ) + \hat { y } \frac { M _ { y } } { 2 } v _ { y } ( n ) \right ] \\ & \quad + \, \hat { z } \frac { M _ { z } } { 2 } v _ { z } ( n ) \right \} e ^ { - i q \cdot ( r _ { n } - r _ { 1 } ) }, \quad \text{(B2)} \\ \dots \, \vec { \epsilon } \quad \vec { \epsilon } \quad \dots \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \,.$$

with M -q , n = M ∗ q , n as the magnetic moment distribution is real.

## Appendix C: Magnetic Structure Factors

This section gives the analytic expressions for the magnetic structure factors for all the magnetic basis vectors.

Starting from the general expression in eq. (2) we obtain the magnetic structure factors for the C A , , and G basis vectors as follows

$$\text{vectors as follows} \\ \mathcal { F } ^ { C } ( Q ) = \begin{cases} \ 1 6 M _ { \pm q, 1 } \cos \frac { \pi l } { 6 } e ^ { i \xi _ { \pm } }, \ h, k, l \text{ all even and} \\ h + k + l = 4 p + 2, p \text{ integer} \\ \ 1 6 M _ { \pm q, 1 } \cos \frac { \pi l } { 6 } \sin \frac { \pi ( h + l - k ) } { 4 } e ^ { i ( \chi _ { \pm } + \pi / 2 ) }, \\ h, k, l \text{ all odd} \\ 0, & \text{ otherwise} \end{cases}$$

$$\mathcal { F } ^ { A } ( Q ) = \begin{cases} \ 1 6 M _ { \pm q, 1 } \sin \frac { \pi l } { 6 } e ^ { i ( \xi _ { \pm } - \pi / 2 ) }, \ h, k, l \, \text{all even and} } & \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \end{cases} \\ \mathcal { F } ^ { A } ( Q ) = \begin{cases} \ 1 6 M _ { \pm q, 1 } \sin \frac { \pi l } { 6 } e ^ { i ( \xi _ { \pm } - \pi / 2 ) }, \ h, k, l \, \text{all even and} } & \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \quad \ 1 6 M _ { \pm q, 1 } \sin \frac { \pi l } { 6 } \cos \frac { \pi ( h + l - k ) } { 4 } e ^ { i ( \chi _ { \pm } - \pi / 2 ) }, \\ \ h, k, l \, \text{all odd} } & \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\, \ \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \quad \ 0, & \quad \text{otherwise} & \quad \quad \quad$$

$$\mathcal { F } ^ { G } ( Q ) = \begin{cases} 1 6 M _ { \pm q, 1 } \sin \frac { \pi l } { 6 } e ^ { i ( \xi _ { \pm } + \pi / 2 ) }, \ h, k, l \, \text{all even and } } \text{entation} \\ h + k + l = 4 p + 2, p \, \text{integer} \\ 1 6 M _ { \pm q, 1 } \sin \frac { \pi l } { 6 } \sin \frac { \pi ( h + l - k ) } { 4 } e ^ { i ( \chi _ { \pm } + \pi ) }, & \text{the a ve} \\ h, k, l \, \text{all odd} \\ 0, & \text{otherwise}, & \text{the sens} \\ \text{ifed} \,. \\ \text{ing plan} \\ \end{cases}$$

where to obtain an analytic closed-form we used the 'ideal' iridium positions z =5/8+1/12.

## Appendix D: Azimuth Dependence of the Magnetic Resonant X-ray Diffraction Intensity

Here we obtain the azimuth dependence of the magnetic x-ray diffraction intensity and show how it can be used to determine the relative phase between magnetic basis vectors. As discussed in Sec. III C the magnetic diffraction intensity is proportional to |F | ‖ 2 , where F ‖ = F · ˆ k ′ is the projection of the structure factor vector onto the scattered beam direction. To calculate this dot product it is convenient to express both vectors in terms of their components in the Laboratory reference frame 25 illustrated in Fig. 9 and defined by the orthonormal set of vectors ˆ U 2 = ˆ k × ˆ k ′ , ˆ U 3 = -ˆ Q and ˆ U 1 = ˆ U 2 × ˆ U 3 , i.e. ˆ U 2 is normal to the scattering plane and ˆ U 1 is in the scattering plane making an angle θ with both k and k ′ , where 2 θ is the total scattering angle. The scattered beam direction then can be written as ˆ k ′ = (cos θ , 0, -sin θ ) L , where the L subscript indicates that the components are with respect to the Laboratory basis of ˆ U 1,2,3 vectors, as opposed to the Crystal basis, defined by the orthonormal set of vectors ˆ x y , ˆ , ˆ z (along the orthorhombic axes a b , , c ).

The magnetic structure factor vector F defined in (2) is most naturally expressed in terms of the Crystal basis as ( F x , F y , F z ) C . The components in the Laboratory basis ( F 1 , F 2 , F 3 ) L are obtained via a linear transformation

FIG. 9: Schematic diagram of the x-ray diffraction experiment at zero azimuth (Ψ = 0). k and k ′ are the incident and scattered x-ray beam wavevectors, and Q = k ′ -k is the scattering wavevector. Main (bottom) panel shows the orientation of the Laboratory frame unit vectors ˆ U 1,2,3 and the orthorhombic crystal axes ( a b c , , ). Top left inset: view along b showing the angle ζ subtended from the scattering plane by the a vector. Right inset: view along Q . During an azimuth scan the crystal axes are rotated around Q by an angle Ψ (in the sense show in the figure), the reference position Ψ = 0 is defined as the azimuth angle for which b lies in the scattering plane making the smallest angle with the incident beam direction ˆ k .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0246v1-figure-873fc21f-2d84-4b73-90d1-54d9bbc44a14.jpg)

encoded in a matrix T (Ψ) via

$$\left [ \begin{array} { c } \mathcal { F } _ { 1 } \\ \mathcal { F } _ { 2 } \\ \mathcal { F } _ { 3 } \end{array} \right ] = T ( \Psi ) \left [ \begin{array} { c } \mathcal { F } _ { x } \\ \mathcal { F } _ { y } \\ \mathcal { F } _ { z } \end{array} \right ], \quad \quad \text{(D1)}$$

where we have explicitly introduced the dependence of the transformation matrix on the azimuth angle Ψ, which defines the rotation angle of the crystal axes around the scattering wavevector Q . To find the transformation matrix we first consider the reference position of zero azimuth, Ψ = 0, depicted in Fig. 9. In this case the orientation of the crystal axes with respect to the Laboratory axes is specified by two angles η and ζ . η is the angle between the b -axis and the scattering wavevector Q , given by

$$\eta = \cos ^ { - 1 } \hat { Q } \cdot \hat { b }$$

and ζ is the angle between the a -axis and the scattering plane, given by

$$\zeta = \tan ^ { - 1 } \frac { Q \cdot \hat { c } } { Q \cdot \hat { a } } \.$$

Using the diagram in Fig. 9 the transformation matrix is

obtained as

$$T ( 0 ) = \left [ \begin{array} { c c c } - \cos \zeta \cos \eta & \sin \eta & - \sin \zeta \cos \eta \\ \sin \zeta & 0 & - \cos \zeta \\ - \cos \zeta \sin \eta & - \cos \eta & - \sin \zeta \sin \eta \end{array} \right ]. \quad ( D 2 )$$

When moving away from the reference position to a finite azimuth angle Ψ the coordinates in the Laboratory frame need to be multiplied by the rotation matrix such that the resulting transformation matrix from the Crystal to the Laboratory bases in (D1) becomes

$$T ( \Psi ) = R ( \Psi ) \ T ( 0 ). \text{ \quad \ \ } \text{(D4)}$$

$$R ( \Psi ) = \left [ \begin{array} { c c c } \cos \Psi & \sin \Psi & 0 \\ - \sin \Psi & \cos \Psi & 0 \\ 0 & 0 & 1 \end{array} \right ], \quad \text{(D3)} \quad \text{Ex}$$

Explicitly, substituting (D2) and (D3) into (D4) yields

$$T ( \Psi ) = \left [ \begin{array} { c c c } \sin \Psi \sin \zeta - \cos \Psi \cos \zeta \cos \eta & \cos \Psi \sin \eta & - \cos \Psi \sin \zeta \cos \eta - \sin \Psi \cos \zeta \\ \sin \Psi \cos \zeta \cos \eta + \cos \Psi \sin \zeta & - \sin \Psi \sin \eta & \sin \Psi \sin \zeta \cos \eta - \cos \Psi \cos \zeta \\ - \cos \zeta \sin \eta & - \cos \eta & - \sin \zeta \sin \eta \end{array} \right ]. \quad \quad ( D 5 )$$

Once the structure factor components are obtained in the Laboratory basis, then the projection of the structure factor vector along the scattered beam direction is readily obtained as

$$\mathcal { F } _ { | | } = \cos \theta \ \mathcal { F } _ { 1 } - \sin \theta \ \mathcal { F } _ { 3 }. \quad \quad ( \text{D6} ) \quad \stackrel { \cdots } { \text{of t} }$$

To illustrate the azimuth dependence of the magnetic x-ray diffraction intensity we consider the magnetic peak at ( -6, 6, 8)+ q in Fig. 5a). Here only the A x and F z basis vectors contribute, therefore F = ( F x , 0, F z ) C where the structure factor components F x and F z are of the type F A and F F in (C2,3), respectively. Using (D6) combined with (D1) and (D5) we obtain

The cross-term is therefore directly sensitive to the phase difference between the two contributing basis vectors. If F z and A x are in-phase or π out-of-phase (denoted as F z ± A x ) the product F F z ∗ x is purely imaginary, so the cross-term cancels and the intensity is given by the sum of the contributions from the two basis vectors taken separately. This is clearly inconsistent with the data (purple curve in Fig. 5a)). A phase difference of ± π/ 2 between the two basis vectors (denoted as F z ± iA x ) gives a finite cross-term, which changes sign between those two cases, and the data in Fig. 5a) identifies the basis vector combination as F z + iA x (red solid line).

$$\mathcal { F } _ { | | } ( \Psi ) = g _ { x } ( \Psi ) \ \mathcal { F } _ { x } + g _ { z } ( \Psi ) \ \mathcal { F } _ { z },$$

where the geometric factors are

$$g _ { x } ( \Psi ) \ = \ & ( \sin \Psi \sin \zeta - \cos \Psi \cos \zeta \cos \eta ) \cos \theta + \\ & \cos \zeta \sin \eta \sin \theta \\ g _ { z } ( \Psi ) \ = \ & \sin \zeta \sin \eta \sin \theta - \\ & \quad ( \cos \Psi \sin \zeta \cos \eta + \sin \Psi \cos \zeta ) \cos \theta.$$

The magnetic scattering intensity is then proportional to

$$| \mathcal { F } _ { | | } ( \Psi ) | ^ { 2 } & = | g _ { x } ( \Psi ) \ \mathcal { F } _ { x } | ^ { 2 } + | g _ { z } ( \Psi ) \ \mathcal { F } _ { z } | ^ { 2 } \\ & + 2 \text{Re} ( \mathcal { F } _ { z } \mathcal { F } _ { x } ^ { * } ) \ g _ { x } ( \Psi ) g _ { z } ( \Psi ),$$

where Re() means the real part. The first two terms give the sum of the separate contributions from the two basis vectors, whereas the last (cross-term) is sensitive to the relative phases between the two basis vectors. Using eqs. (3,C2) the cross-term is proportional to

$$\mathcal { F } _ { z } \mathcal { F } _ { x } ^ { * } \in \begin{cases} \mp i \sin \frac { \pi l } { 6 } \cos \frac { \pi l } { 6 } &, & F _ { z } \pm A _ { x } \\ \pm \sin \frac { \pi l } { 6 } \cos \frac { \pi l } { 6 } &, & F _ { z } \pm i A _ { x }. \end{cases}$$

TABLE V: Equivalence between the iridium sites in γ -Li 2 IrO 3 ( Cccm space group) with those in β -Li 2 IrO 3 ( Fddd space group) (when both are projected onto the ac plane). Last column shows additional translations to be applied to primitive lattice sites in the Fddd structure (column 2) in addition to a common translation by (-1/8,0,1/8) to get the positions in the Cccm structure (column 1). For each structure the 'ideal' iridium positions have been used.

| Site Cccm   |   Site Fddd | Translation   |
|-------------|-------------|---------------|
| 1           |           4 | -(0,0,1)      |
| 2           |           3 | (0,0,0)       |
| 3           |           3 | (0,1/2,1/2)   |
| 4           |           4 | -(0,1/2,1/2)  |
| 1 ′         |           2 | (1/2,1/2,0)   |
| 2 ′         |           1 | (1/2,0,-1/2)  |
| 3 ′         |           1 | (1/2,1/2,0)   |
| 4 ′         |           2 | (1/2,1,1/2)   |

## Appendix E: Mapping of magnetic basis vectors between β - and γ -Li IrO 2 3

In this section we describe a formal mapping of the magnetic basis vectors between the notation used for the β - (space group Fddd ) and γ -polytypes of Li 2 IrO 3 (space group Cccm ). The iridium lattice arrangement in both crystal structures is made up of zig-zag chains verticallylinked along the c -axis, with the only difference being that in the β -polytype successive chains along c alternate in orientation between the a ± b directions, whereas in the γ -polytype a pair of successive chains are parallel and this direction then alternates for the next pair of parallel chains and so on, between the a ± b directions. When both structures are viewed in projection onto the ac plane [as shown in Fig. 6a) and Fig. 4 of Ref. 20] the two iridium lattice arrangements then look identical. This suggests that a formal mapping can be constructed between the notations used to describe the sites of the primitive cells of the two materials, and also between the magnetic basis vectors, such that one can then describe a magnetic structure equivalently in the notation of one or the other space group (of course the magnetic structures would only be 'equivalent' up to the projection of the sites onto the b -axis). Table V lists the equivalence between the sites in the primitive cell of the two structures (up to the y -coordinate of sites), where for the γ -polytype we follow the notation in Ref. 20 (two iridium sublattices with sites 1-4 and 1 -4 ). ′ ′ Using this site equivalence one can construct a mapping of the magnetic basis vectors as listed in Table VI. The magnetic structure described by

TABLE VI: Mapping of magnetic basis vectors between the Fddd and Cccm notations.

| Basis vector        | Basis vector                  |
|---------------------|-------------------------------|
| Fddd (4-site [1-4]) | Cccm (8-site,[1-4,1 ′ -4 ′ ]) |
| A [+ -- +]          | ( A , - A ) [+ -- + - ++ - ]  |
| C [+ + -- ]         | ( - F , F ) [ ---- ++++]      |
| F [+ +++]           | ( F , F ) [+ +++++++]         |

the basis vector combination ( iA x , iC y , F z ) in the Fddd setting then corresponds, in the notation of the Cccm setting, to i ( A , -A ) x , i ( -1) m ( F , -F ) y , ( F F , ) z with m = 1.
A. Biffin R. D. Johnson Sungkyun Choi F. Freund S. Manni A. Bombardi P. Manuel P. Gegenwart R. Coldea
Unconventional magnetic order on the hyperhoneycomb Kitaev lattice in $β$-Li2IrO3: full solution via magnetic resonant x-ray diffraction
The recently-synthesized iridate $\beta$-Li$_2$IrO$_3$ has been proposed as a
candidate to display novel magnetic behavior stabilized by frustration effects
from bond-dependent, anisotropic interactions (Kitaev model) on a
three-dimensional "hyperhoneycomb" lattice. Here we report a combined study
using neutron powder diffraction and magnetic resonant x-ray diffraction to
solve the complete magnetic structure. We find a complex, incommensurate
magnetic order with non-coplanar and counter-rotating Ir moments, which
surprisingly shares many of its features with the related structural polytype
"stripyhoneycomb" $\gamma$-Li$_2$IrO$_3$, where dominant Kitaev interactions
have been invoked to explain the stability of the observed magnetic structure.
The similarities of behavior between those two structural polytypes, which have
different global lattice topologies but the same local connectivity, is
strongly suggestive that the same magnetic interactions and the same underlying
mechanism governs the stability of the magnetic order in both materials,
indicating that both $\beta$- and $\gamma$-Li$_2$IrO$_3$ are strong candidates
to realize dominant Kitaev interactions in a solid state material.