Spin-orbit physics of j=1/2 Mott insulators on the triangular lattice
Michael Becker,1 Maria Hermanns,1 Bela Bauer,2 Markus Garst,1 and Simon Trebst1 1Institute for Theoretical Physics, Cologne University, Zulpicher¨ Straße 77, 50937 Cologne, Germany 2Station Q, Microsoft Research, Santa Barbara, CA 93106-6105, USA (Dated: April 16, 2015) The physics of spin-orbital entanglement in effective j = 1/2 Mott insulators, which have been experi- mentally observed for various 5d transition metal oxides, has sparked an interest in Heisenberg-Kitaev (HK) models thought to capture their essential microscopic interactions. Here we argue that the recently synthesized Ba3IrTi2O9 is a prime candidate for a microscopic realization of the triangular HK model – a conceptually interesting model for its interplay of geometric and exchange frustration. We establish that an infinitesimal Ki- taev exchange destabilizes the 120◦ order of the quantum Heisenberg model. This results in the formation of an extended Z2-vortex crystal phase in the parameter regime most likely relevant to the real material, which can be experimentally identified with spherical neutron polarimetry. Moreover, using a combination of analytical and numerical techniques we map out the entire phase diagram of the model, which further includes various ordered phases as well as an extended nematic phase around the antiferromagnetic Kitaev point.
I. INTRODUCTION yond the hexagonal lattice, triggered mainly by the synthe- sis of novel Iridate compounds, which includes e.g. the sis- 16 17 ter compounds β-Li2IrO3 and γ-Li2IrO3 that form three- The physics of transition metal oxides with partially filled dimensional Ir lattice structures. In this manuscript, we turn to 5d shells is governed by a largely accidental balance of elec- 18 the recently synthesized Iridate Ba3IrTi2O9 and argue that tronic correlations, spin-orbit entanglement, and crystal field it realizes a Heisenberg-Kitaev (HK) model on a triangular effects, with all three components coming up roughly equal lattice36. This model is of deep conceptual interest as it ex- in strength. With different materials exhibiting slight tilts to- hibits a subtle interplay of the two elementary sources of frus- wards one of the three effects a remarkably broad variety of tration – geometric frustration arising from its non-bipartite quantum states has recently been suggested, which includes lattice structure as well as exchange frustration arising from exotic states such as Weyl semi-metals, axion insulators, or 1 the Kitaev couplings. The ground states of the classical HK topological Mott insulators . A particularly intriguing sce- model on the triangular lattice have been already addressed by nario is the formation of Mott insulators in which the lo- Rousochatzakis et al.20. With the help of Monte Carlo simula- cal moments are spin-orbit entangled Kramers doublets. An tions these authors demonstrated that a small but finite Kitaev example are the j = 1 Mott insulators observed for vari- 2 exchange in addition to an antiferromagnetic Heisenberg in- 2–4 ous Iridates . The Iridium valence in the latter typically is teraction stabilizes a -vortex crystal. The -vortices can 4+ 5 Z2 Z2 Ir corresponding to a 5d electronic configuration. With be viewed as defects of the SO(3) order parameter associated the crystal field of the octahedral IrO oxygen cage splitting 6 with the 120◦ ordering of the antiferromagnetic Heisenberg off the two eg levels, this puts 5 electrons with an effective model. It is important to note that this physics plays out near 1 s = 2 magnetic moment into the t2g orbitals, which entan- the Heisenberg limit of the HK model – the relevant micro- gled by strong spin-orbit coupling leaves the system with a scopic parameter regime for all Iridates synthesized so far. As fully filled j = 3 band and a half-filled j = 1 band2,3,5. The 2 2 such Ba3IrTi2O9 is a prime candidate to observe this exotic smaller bandwidth of the latter then allows for the opening of phase. a Mott gap even for the relatively moderate electronic corre- After a discussion of the material aspect of Ba IrTi O and 1 3 2 9 lations of the 5d compounds. Interest in such j = 2 Mott a motivation of the HK model in Sec.II, we examine in insulators has been sparked by the theoretical observation6–8 Sec. III the formation of a Z2-vortex crystal from the analyti- that the microscopic interaction between their spin-orbit en- cal perspective of an expanded Luttinger-Tisza approximation tangled local moments not only includes an isotropic Heisen- and quantitatively describe its experimental signatures in po- berg exchange but also highly anisotropic interactions whose larized neutron scattering experiments. In Sec.IV we address easy axis depends on the spatial orientation of the exchange the full phase diagram of the HK model and discuss the var- 9 arXiv:1409.6972v2 [cond-mat.str-el] 15 Apr 2015 path tracing back to the orbital contribution of the moments . ious phases with the help of analytical as well as numerical In a hexagonal lattice geometry, as it is found for the lay- methods. Finally, in Sec.V we close with a summary. ered Na2IrO3 and α-Li2IrO3 compounds, these anisotropic interactions provide an implementation of the celebrated Ki- 10 taev model known for its spin liquid ground states. A trove II. MATERIAL PHYSICS OF Ba3IrTi2O9 of experimental data11, ab initio calculations12, and model 13 simulations for these hexagonal systems has spurred an on- 4+ Ba3IrTi2O9 forms layers of Ir ions in a triangular geome- going discourse illuminating the actual spin-orbital ordering try, which are separated from each other by two layers of Ti4+ mechanism in these materials. ions. An important characteristic of the Ir layer geometry il- Much recent activity14,15 has been targeted towards the lustrated in Figs.1 a) and b) is that it exhibits the two neces- physics of j = 1/2 Mott insulators for lattice geometries be- sary ingredients for Kitaev-type exchange couplings. First, 2
a) crystal structure b) single Iridium layer c) Ir-O-O-Ir exchange path
z y Ir O zˆ yˆ x zˆ yˆ zˆ xˆ yˆ Ti Ba xˆ view along (111) xˆ view along (311)
Figure 1. (Color online) a) Crystal structure of Ba3IrTi2O9. b) View of single Iridium layers from two different perspectives. Within the plane the x, y, and z exchange paths are indicated by the grey planes. The planes labeled by x (y, z) are normal to the coordinate axis xˆ (yˆ, zˆ). c) The exchange between the Iridium moments (blue) is mediated by two coplanar exchange paths. every pair of Iridium ions is coupled via two separate ex- ground state of the antiferromagnetic Heisenberg Hamilto- change paths as indicated in Fig.1 c) leading to a destruc- nian on the triangular lattice, which corresponds to couplings tive interference and subsequent suppression of the isotropic JH > 0 and JK = 0 for Hamiltonian (1), is characterized by 6–8 21 Heisenberg exchange . In comparison to the tricoordinated a 120◦ ordering of spins . At the classical level this order- ˆ Iridates (Na,Li)2IrO3, which exhibit Ir-O-Ir exchange paths, ing is captured by a spin orientation Si = SΩ(ri) with the ˆ the triangular Ba3IrTi2O9 exhibits somewhat longer Ir-O-O-Ir unit vector Ω120◦ (r) = e1 cos (Q r) + e2 sin (Q r) where exchange paths resulting in an overall lessening of the mag- the commensurate wavevector Q connects· the center· with a netic exchange strength. Second, the three principal bond di- 4π corner of the Brillouin zone, Q = 3 (1, 0). The orthonormal rections of the triangular lattice structure cut through three frame ei with i = 1, 2, 3 and e3 = e1 e2 constitutes an different edges of the IrO6 oxygen cages resulting in a distinct SO(3) order parameter. The energy per site× for this classical 6–8 locking of the exchange easy axis along the three directions state is given by as illustrated Fig.1 a) and ultimately giving rise to the three components of the Kitaev exchange. Note that the Ir layer is normal to the (111) direction, hence the three directions are 2 1 ε ◦ = S 3J + J . (2) all equivalent. The description of the microscopic physics is 120 − 2 H K thus given in terms of a Heisenberg-Kitaev (HK) Hamiltonian
X X γ γ 22 = J S S + J S S , (1) Crucially, the 120◦ ordering possesses Z2 vortices as HHK H i · j K i j ij γ ij topologically stable point defects, which can be understood h i kh i by considering the first homotopy group of its order parame- where Si is a spin-operator located on site i of the trian- ter Π1(SO(3)) = Z2. T gular lattice spanned by the lattice vectors ax = (1, 0) , T ay = ( 1/2, √3/2) , and az = ax ay, see Fig.2 a). Here and− in the following, we measure− lengths− in units of the lattice constant a. The first term is the standard Heisenberg ◦ coupling, JH , that describes an SU(2) invariant interaction be- A. Kitaev interaction destabilizes 120 ordering tween the spin-orbit entangled j = 1/2 moments on nearest- neighbor lattice sites. The Kitaev interaction, JK , on the other hand, explicitly breaks spin-rotation invariance and acts only For any finite JK the 120◦ state becomes immediately un- between single components, Sγ , of adjacent spins. The pre- stable with respect to fluctuations, which we demonstrate in cise component depends on the link between the lattice sites, the following. We parametrize the fluctuations with the help T see Fig.2 a); for our particular choice here, the γ-components of two real fields π(r) = (π1(r), π2(r)) , of spins interact via JK if sites are connected by a lattice vec- tor aγ with γ = x, y, z. q 2 Ωˆ(r) = Ωˆ ◦ (r) 1 (π(r)) (3) 120 − r e Q r e Q r r e ◦ + π1( )( 1 sin ( ) + 2 cos ( )) + π2( ) 3, III. 120 ORDER AND Z2-VORTEX CRYSTAL − · ·
We will start our discussion of the ground states of Hamilto- so that Ωˆ 2(r) = 1 is mantained. Plugging this Ansatz in the nian (1) by first elucidating the magnetic structure around the Hamiltonian and expanding up to second order in the fluctua- (2) antiferromagnetic Heisenberg point, where an extended Z2 tion fields one obtains for the energy = Nε120◦ + with vortex crystal phase is found in agreement with Ref.20. The N denoting the number of lattice sites.E The fluctuationE part 3
a) b) It becomes maximally negative for a wavevector ay x x JK 4 X γ 2 Si Sj kinst = aγ sin(Q aγ )(e ) (7) az ax y y J 3 · 3 Si Sj Q H γ=x,y,z z z Si Sj T JK 1 h y 2 z 2 x 2 i z 2 y 2 = (e3) + (e3) 2(e3 ) ) , (e3) (e3) , JH √3 − − Figure 2. (Color online) a) The triangular lattice with the three lattice vectors aγ . Solid, dashed and dotted bonds carry the three distinct that can be expressed in terms of the normal e3. In the spe- Kitaev interactions, see text. b) First Brillouin zone of the triangular cial case where the spins of the 120◦ ordering are confined lattice. The position and size of the coloured dots indicate the posi- within the x-y plane and e3 = ˆz, this wavevector is just given T tion and weight of Bragg peaks, respectively, expected in the static by kinst = JK /JH (1/√3, 1) . So it is the tilting Goldstone spin structure factor for the vortex crystal. Each color corresponds Z2 modes that trigger the instability of the 120◦ antiferromag- to a different spin-component as listed in panel a). netic ordering in the presence of a finite Kitaev interaction JK . reads
2 B. Incommensurate antiferromagnet: Z2-vortex crystal (2) X 2 JH S X = ε ◦ (π(r )) (π π 2π π ) E − 120 i − 2 1i 1j − 2i 2j i ij Indeed, allowing for a slowly spatially varying orthogonal h i X h frame e (r) one finds in the limit J J the effective 2 γ γ γ γ Qr Qr i K H + JK S e3 e3 π2iπ2j + e1 e1 cos( i) cos( j) energy functional = R d2r with| | γ ij E L kh i 2 γ γ γ γ 3JH S X 2 + + e e sin(Qr ) sin(Qr ) e e sin(Q(r + r )) π π = e−(r) 2iqK aγ e (r), 2 2 i j − 1 2 i j 1i 1j L 4 γ −∇ − · ∇ γ h i i γ=x,y,z + eγ sin(Qr ) + eγ cos(Qr ) eγ π π + (i j) , (8) − 1 i 2 i 3 1i 2j ↔ (4) where e± = (e ie )/√2. The Kitaev interaction induces a 1 ± 2 coupling qK = 2JK /(√3JH ) to constant gauge fields given with the abbreviation πai = πa(ri) for a = 1, 2. The fluctu- by the triangular lattice vectors aγ , that can be identified as ation eigenmodes are determined with the help of the Fourier 20 1 P ikri Lifshitz invariants as previously pointed out in Ref. . The transform πa(ri) = √ k 1.BZ e πa(k). In the absence N ∈ magnetization can thus minimize its energy by allowing for of the Kitaev interaction, JK = 0, one obtains a spatial modulation of the SO(3) order parameter on large length scales proportional to 1/q J /J . 2 h K H K (2) JH S X ∝ = (1 cos(kaγ ))π1∗(k)π1(k) JK =0 2 E k∈1.BZ − γ=x,y,z 1. Luttinger-Tisza approximation i + (1 + 2 cos(kaγ ))π2∗(k)π2(k) (5) The character of this modulated classical ground state can be obtained by minimizing the Hamiltonian treating the or- with πa∗(k) = πa( k). Whereas the π1 mode becomes soft − ˆ 2 at the center of the Brillouin zone, i.e., at k = 0, the en- thonormal constraint, ei ej = δij, or equivalently, Ω = 1, · 24 ergy of the π2 mode vanishes at its edge, e.g. for momenta within an improved Luttinger-Tisza approximation . The lat- k = Q. The zero modes π (k = 0) and π ( Q) thus ter is a good approximation for large length scales qK r 1, ± 1 2 ± | | identify three Goldstone modes that correspond to a long- or, alternatively, for small momenta q qK . Note that this | | wavelength rotation and tilting of the local orthogonal frame, latter limit does not commute with JK 0, and, as a conse- → respectively. In particular, the energy dispersion of the tilting quence, does not smoothly connect with the Heisenberg point. tilt 2 P We start with the functional mode εk JK =0 = JH S γ=1,2,3(1+2 cos(k aγ )) close to | tilt · 2 3 2 momentum Q possesses the form εQ+k JK =0 JH S 8 k . = (9) | ≈ E Adding a finite Kitaev coupling JK immediately results in X X X J S2 Ωˆ Ωˆ + J S2 Ωˆ γ Ωˆ γ λ (Ωˆ 2 1). a negative energy eigenvalue and, therefore, destabilises the K i · j K i j − i i − ij γ ij i 120◦ ground state. We can still diagonalize for the eigenen- h i kh i ergies perturbatively in J . In lowest order and in the long- K The unit length of the vector Ωˆ is locally imposed with the wavelength limit the zero modes do not hybridize, and we ob- i help of the Lagrange multipliers λ . Upon spatial Fourier tain for the tilting mode a dispersion relation that is given in i transformation, Ωˆ(r) = P eiqrΩˆ , the functional takes the the long-wavelength limit, k Q , by q q | | | | form X X 3 X γ ˆ α αβ ˆ β ˆ α ˆ α εtilt J S2 k2 2J S2 k a sin(Q a )(e )2. /N = Ω q (q)Ωq + λ q pΩpΩq λ0, Q+k ≈ H 4 − K · γ · γ 3 E − J − − − γ=x,y,z q q,p (6) (10) 4
αβ where λ0 = λq q=0. The matrix J possesses only diagonal tmin(2aη aγ ) with η = γ. The resulting Bragg peaks in entries with | the− static structure− factor6 are visualized in Fig.2 b), which nicely agrees with previous numerical findings for the clas- αα 2 20 (q) = JH S cos(ax q) + cos(ay q) + cos(az q) sical model . The relative weight of secondary and primary J · · · Bragg peaks are predicted to be 1/42 = 1/16 within the above + J S2 cos(a q) (11) K α · approximation. We find that the corresponding energy is in- dependent of the choice of origin r = (x , y )T as well as and J αβ(q) = 0 for α = β. At the Heisenberg point 0 0 0 6 the phase φ. JK = 0, the diagonal components of the matrix are min- ˆ imal for momenta at the corner of the Brillouin zone, e.g. In the Luttinger-Tisza approximation the length of the Ω vector is compromised to differ from unity, P3 (Ωˆ α(r))2 = q = Q, thus leading to 120◦ ordering. A finite J , however, α=1 t K ˆ 6 favours in general incommensurate order with wave vectors 1. Whereas the length Ωt(r) varies in space it nevertheless αα | | ˆ away from Q as J (q) become minimal for momenta of the remains always finite so that the orientation of Ωt(r) is al- (1) ways well defined. Note that in the limit J 0 the dis- form qα = Q t aα with t R. On the other hand, Fourier K − ∈ 2 → components, Ωˆ α , of the spin with such incommensurate tance tmin 0 and εLT(0) = S 3JH /2 recovers the exact (1) → − qα ground state energy whereas the state itself, Ωˆ α (r), does not wave vectors induce finite Fourier components, λ (1) with t=0 2qα reproduce the 120 ordering as expected. α = 1, 2, 3, of the Lagrange multiplier. Finite Lagrange± mul- ◦ tipliers λ (1) , in turn induce two finite secondary Fourier 2qα ± ˆ α (2) 2. Vector chirality and vortices components Ω (2) with qα,β = Q t (2aβ aα) where Z2 qα,β − − β = α and so on. 6In the following, we discuss a Luttinger-Tisza approxima- It turns out that the approximate classical ground state (13) tion where we limit ourselves to the lowest finite Fourier com- corresponds to a triangular lattice of condensed Z2 vortices, ˆ α ˆ α thus confirming the numerical results of Ref. 20. This is best ponents Ω (1) and Ω (2) for the spin and λ0 and λ (1) q q 2qα α α,β ± seen by defining chirality vectors on upward pointing triangles for the Lagrange multiplier; all higher Fourier modes are ne- of the lattice glected. In principle, this approximation can be systematically improved by including higher order modes. Minimizing the κ(r) = 2 S S + S S + S S . 3√3 r × r+ax r+ax × r+ay r+ay × r functional (10) within this approximation we obtain for the (14) energy per site The length of κ(r) measures the rigidity of the local 120◦ or- S2 h π + 6t π 3t 22 ε (t) = J cos + 17 sin − (12) dering and it vanishes at the center of each Z2 vortex . The LT − 9 H 3 6 chirality vector profile, that derives from Eq. (13), is shown in π + 6t π + 15t Fig.3 and clearly reveals the vortex crystal. + 8 sin + sin Z2 6 6 Within our Luttinger-Tisza approximation we find three π + 6t π + 6ti zero modes for the Z2 vortex crystal represented by the phase + JK cos + 8 sin , 3 6 φ and the vector r0. The latter are expected as the vortex crys- tal spontaneously breaks translational symmetry so that a con- which still depends on the parameter t that quantifies the dis- stant shift of the origin r0 does not cost any energy. The corre- tance of the primary Bragg peak from the corner of the Bril- sponding low-energy excitations are just the effective acous- (1) louin zone, qα = Q t aα. The value of tmin identifying tic phonon excitations of the vortex crystal. If the coupling − the position of the minimum of the function (12) finally de- between the two-dimensional atomic triangular lattice planes termines the ground state energy εLT(tmin). This analytical of Ba3IrTi2O9 is sufficiently small, these low-energy modes estimate for the ground state energy is found to be in excel- will destroy true long-range order of the Z2 vortex crystal at lent agreement with numerical estimates obtained from Monte any finite temperature that will be reflected in a characteris- Carlo simulations discussed in Sec.IVC. tic broadening of the Bragg peaks in the structure factor of The corresponding state is given by Fig.2 b). 4S n Sγ (r) Re eiφ (13) ≈ 3√3 × C. Polarized neutron scattering o i(Q taγ ) (r r ) 1 X i(Q t(2aη aγ )) (r r ) e − · − 0 + e − − · − 0 , 4 η=γ The structure factor of the Z2 vortex crystal possesses as a 6 hallmark of the Kitaev interaction a characteristic correlation where Sγ is the γ-component of the spin and the ground state between the positions of the Bragg peaks and the associated is obtained by setting t = tmin. The first term in Eq. (13) spin-components, see Fig.2 b). We suggest to resolve this is the most important, primary Fourier component which also correlation with the help of spherical neutron polarimetry. possesses the smallest deviation of momentum from the cor- The probability that an incoming neutron with spin σin ner of the Brillouin zone, Q. The secondary Fourier compo- is scattered into a spin-state σout is given by the energy- nents have a smaller weight and are shifted further away by integrated scattering cross section σσout,σin (q), where q is the 5
dual 120 order nematic phase dual Z2 vortex crystal Z2 vortex crystal