Twisted Hubbard model for Sr2IrO4: magnetism and possible high temperature superconductivity Fa Wang and T. Senthil, PRL 2011
Iridium oxide materials: various kinds of exotic physics Na4Ir3O8: insulating quantum spin liquid (Okamoto, Takagi et al, 2007) Sr2IrO4: Jeff = 1/2 spin-orbit entangled Mott insulator (B.J. Kim et al, 2008) Other interesting proposals: A2IrO3: quantum spin Hall effect? (Shitade, Nagaosa et al, PRL 2009), realization of Kitaev spin liquid? (Jackeli, Khaliullin, PRL 2009) Ln2Ir2O7 pyrochlore iridates: (i) spontaneous anomalous Hall effect (Nakatsuji et al, 2010) for Ln = Pr; (ii) other Ln: 3d correlated topological insulator? (Pesin, Balents, Nat Phys 2010).
Sr2IrO4: preliminaries K2NiF4 structure similar to La2CuO4 -+(./,0)123)1(.(,'344 56-65371203.289,3 ""! "#! "#! Layers of IrO2 ; Ir4+ ions on a 2d square lattice. Each Ir surrounded by O octahedra. Similarity to cuprates noticed in 1990s: M.A.Subramanian, Gang Cao,...
Sr2IrO4: Mott insulator 24$ H'/0I H'/2I G4 H'/GI H'/ I!-#&88)6$(#&#-$6&'$3-$4-(6%)3-4$&($&$+ Cao et al, PR B 1998 H. Takagi et al, 2008 However Sr2CoO4 and Sr2RhO4 are metallic.
Why a Mott insulator? Ir 4+ is in a 5d 5 configuration. ¼ Cubic environment: big t2g - eg splitting. SO coupling splits t2g level into Jeff = 1/2 and Jeff = 3/2 bands. Jeff = 1/2 band is half-filled and narrow => moderate U can drive it into a Mott insulator B. J. Kim et al, PRL 08
Spin-orbit entangled Mott insulating magnet J eff,+1/2 = 1 3 ( xy, +1/2 + yz, 1/2 + i zx, 1/2) J eff, 1/2 = 1 3 ( xy, 1/2 yz, +1/2 + i zx, 1/2) Jeff moments order into a magnetic state below 240 K. Weak ferromagnetism (from canting of AF moments): 0.14 μb Spin-orbital content confirmed by resonant magnetic X-ray diffraction (B.J. Kim, Takagi et al, Science 2009)
How similar to cuprates?? Despite obvious similarities, some worries 1. Different active orbitals 2. Spin anisotropy and related strong spin orbit coupling Is there a useful effective model which can be compared with, eg, Hubbard model popular for cuprates?
LDA + Spin-Orbit + U bandstructure Jin, Jaejun Yu, et al, PR B 2009 (a) 1 Energy (ev) 0-1 -2-3 Fit to obtain tightbinding model (+ U) in t2g subspace: real orbital diagonal hoppings (i.e only xy -> xy, yz->yz, xz->xz) + on-site spin-orbit term Similar result: H. Watanabe et al, PRL 2010
Y X Y X y x Interpretation of tightbinding model Wang and TS, 2011 Rotation of oxygen octahedra => must carefully distinguish between local X,Y, Z basis and global x,y, z basis. Note 1. Crystal field splitting between t2g and eg is diagonal in local X,Y,Z basis => define t2g orbitals in local orbital basis. 2. If Ir - Ir hopping is mediated by O p-orbitals, symmetry => Ir-Ir hopping diagonal in local orbital basis 3. Global basis for spin S = 1/2 => Ir-Ir hopping real. Conclusion: LDA+SO+U tight-binding fit involves Wannier states with t2g orbitals in local X,Y,Z basis but spin in global basis.
Effective model for Jeff = 1/2 states: strategy Wang and TS, 2011 What are correct Jeff = 1/2 states? Must define both t2g orbitals and spin in local X,Y, Z basis. J z eff =+1/2 = 1 3 (+ XY, XZ, + YZ, ), J z eff = 1/2 = 1 3 ( XY, + XZ, + YZ, ). Effective model in Jeff = 1/2 subspace: Start with full tightbinding t2g tightbinding model, rotate spin from global to local basis, then project onto Jeff = 1/2 states.
Effective model for Jeff = 1/2 states: result H = (t + ɛ α ɛ j t) d j,α d k,α t d j,α d k,α <jk>,α <<<jk>>>,α t d j,α d k,α + U j <<jk>>,α d j, d j, d j, d j, (3) Effective t-t -t -U Hubbard model for Jeff = 1/2 band.
Consequence: pseudospin isotropy At half-filling, expect can describe Mott insulator by usual isotropic pseudospin Jeff = 1/2 model
Coupling to external Zeeman field is `twisted Wang, TS, 2011 Observable moment M related to J through M j,x M j,y = 2µ B cos θ ɛ j sin θ 0 ɛ j sin θ cos θ 0 J j,1 J j,2. 0 0 1 M j,z J j,3 Twisted magnetic field coupling => anisotropic magnetization even though intrinsic physics is isotropic.
Consequences for magnetism in Sr2IrO4 Isotropic pseudospin ordering => observed moment rotates with octahedra rotation. => weak ferromagnetism. Y X Y X y x FM moment = (AF moment)sin θ θ = octahedron tilt angle Known AF moment of J =1/2 Heisenberg model => FM moment 0.14µ B similar to experiment.
Consequences for magnetism (cont d) 1. Measured susceptibility near q = 0 is mixture of uniform and staggered susceptibilities 2. Anisotropy between susceptibility to field along ab-plane versus along c axis. 3. Despite apparent spin anisotropy, spin correlation length above TN should behave like that of isotropic magnet. All of these agree with experiment.
Spin correlation length above TN Shigeki Fujiyama, H. Takagi et al, unpublished B inplane correlation 228 K Sr2IrO4 2dim. Ising BKT 230 K 233 K 235 K 237 K 239 K 244 K 249 K Intensity (log. scale) 227 K 228 K 229 K 230 K 232 K 234 K 239 K 244 K 249 K 254 K 2dim. Heisenberg 2.0 23.0 0.96 0.98 1.00 1.02 1.04 22.5 22.0 21.5 D Diffuse magnetic X-ray scattering inplane correlation 0.92 1.00 1.10 (T-T N ) / T N T/TN 0.008 FIG. 2: Temperature evolutions of inplane (ξ a ) and interplane (ξ c ) spin correlation lengths. Solid 0.004 lines are fitted curves to theoretically proposed formula for S = 1/2 square lattice quantum Heisenberg model (Chakravarty et al. and Makivic) (red), to two dimensional Ising (blue) and xy (green) Thanks to Dr. 0.000 Fujiyama correlations. for permission to show this data ξ a exp(1.5(t T N ) 1/2 ) (BKT). The best fitted curves to these two power-law temperature -0.004 evolution show fatal contradictions to the observed ξ a as shown in Fig. 2. Therefore, in spite of the agreement that symmetry breaking of SU(2) in magnetic correlation causes -0.008.00 1.05 1.10 0.992 0.996 1.000 1.004 1.008 canted moments in the ordered state as argued in Ref. [10], the magnetic correlation in the paramagnetic state is concluded to be predominantly isotropic and the crucial source for three dimensional ordering is a small interplane magnetic correlation with 1.2 mk. This
Compare to cuprate Effective isotropic Hubbard model parameters: t 0.2eV, t t t 4, t 10, U 10t. Similar to cuprate? 1. Roughly same parameter range as cuprate Hubbard model except overall energy scale smaller by 2. 2. t /t positive for iridate - opposite to cuprate
Doped system - superconductivity? If the doped square lattice 2d Hubbard model is superconducting, can we expect doping the iridate will give a `high temperature superconductor? From comparison to cuprate (i) t /t > 0 => Electron doping better than hole doping for iridate (Example: Replace Sr by La) (ii) Overall energy scale smaller by 2 => Tc in range of 50 K???
Differences with cuprates: more coherent c-axis transport? Interlayer tunneling matrix element in cuprate: (cosk x cosk y ) 2 factor which blocks c-axis transport of nodal states. Iridate: different orbital content implies no suppression of nodal c-axis hopping. Better PHYSICAL c-axis coherence? REVIEW B 82, 115117 2010 0 50 100 150 200 250 300 350 400 log ρ (Ω cm) ρ (Ω cm) 10 5 10 3 10 1 10-1 1.6 10-1 1.2 10-1 810-2 410-2 (a) (b) T a Sr 2 IrO 4-δ δ =0 I=0.05mA T MI 3.5 10-2 010 0 0 50 100 150 200 250 300 350 400 T(K) Monday, April 300 25, 2011 ρ c ρ a ρ c ρ a ρ c (Ω cm) T a ρ c ρ a ~T ρ a (Ω cm) 2.5 10-2 010 0 0 5 10 15 20 T(K) δ 0.04 610-4 410-4 210-4 I=0.05mA Overall resistivity anisotropy of about 100 in undoped insulator
log ρ a (Ω cm) M(µ B /Ir) ρ c (Ω cm) 10-2 KORNETA et al. Experimental prospects? Recent report: Metal-insulator transition with 0 50 100 150 200 250 300 350 400 ELECTRON-DOPED Sr 2 IrO 4 0 0 50 100 150 200 250 300 350 10 2 (a) δ=0 10 0 FIG. 1. Color online A schematic of the Ir1-O2-Ir1 bond angle 0< δ <0.04 10-4 at T=90 K for a =0 andsr b IrO 0.04 the changes in are δ 0.04 2 4-δ exaggerated for clarity ; a schematic of the band structure corresponding to c =0 and d 0.04. Note that decreases for 10-6 =0 but slightly increases for 0.04 as temperature is lowered. 0.05 δ =0 0.03 0.04 III. EXPERIMENTAL δ 0.04 RESULTS M a ρ a 0.03 Sr 0.02 2 IrO 4 crystallizes in µ a H=0.2T o reduced tetragonal structure space-group I4 0.02 1 /acd due to a rotation of the IrO 6 octahedra about the c axis by an angle 11. 2 4 This 0.01rotation corresponds c to a distorted in-plane T Ir1-O2-Ir1 bond angle that c M 0.01 decreases with decreasing temperature. (b) Our single-crystal 0 0 x-ray diffraction data confirm this trend for stoichiometric 10-1 ρ (Ω cm) 1.4 10-1 crystals =0 and also reveal a that for =0.04 increases slightly 7Twith decreasing 10-2 temperature from 157.028 at 295 K to 157.072 at 90 T -1/4 (K -1/4 ) 110-1 10 K, -3 and the latter angle is significantly larger than that for 0.23 =0, i.e., 0.25 =156.280 0.27 at 90 K, as shown in Fig. 1 and Table I. The increment =0.792 is 610-2 large for such a small oxygen depletion. Moreover, the volume of µ theh=0t o unit cell V for 0.04 δ contracts by an astonishing 0.14% compared to that for =0 see (c) Table I. These data 210-2 0 indicate 50 that 100dilute 150 oxygen 200 250 vacancies 300 relax 350 and reverse its temperature dependence T (K) with increasing while significantly reducing the structural distortion at low T. No such changes in the lattice parameters would be observable in x-ray diffraction data should the oxygen depletion be confined to the crystal surface and not uniformly distributed within the bulk. Oxygen doping 10 3 ρ a (Ω cm) log ρ (Ω cm) ρ (Ω cm) S c (µv/k) 10 5 10 1 10-1 1.6 10-1 1.2 10-1 (a) (b) 10-9 Cornetta, Cao et al, PR B 2010 PHYSICAL REVIEW B 82, 115117 2010 Sr 2 IrO 4-δ PHYSICAL δ =0 REVIEW B 82, 115117 2010 I=0.05mA 3.5 10-2 T T(K) MI 810-2 ρ c T a 410-2 δ 0.04 I=0.05mA ρ a 010 0 0 50 100 150 200 250 300 350 400 T(K) 300 (c) 200 δ =0 Sr IrO 2 4-δ δ 0.04 100 FIG. 3. Color online Temperature dependence of a log a for several representative values of, b the a- and c-axis magnetization, M a and M c for The =0 Ir1-O2-Ir1 thin lines bond angle 0.04 is an thick important lines focus and of a this right scale forstudy, 0.04, as itand controls c c the at 0hopping and 7 T. of Inset the in 5d c : electrons log a and vs T 1/4. superexchange interactions between Ir atoms via the bridging O sites 10 and is therefore expected to influence physical properties. For example, the a-axis resistivity a c-axis retransition near T MI =105 K, resulting in a reduction of a c by a factor of 10 4 10 1 from just below T MI to T =1.8 K Fig. 2 b. The strong low-t anisotropy reflected in the values of a and c and their sensitivity to are consistent with a nearly two-dimensional 2D, strongly correlated electron system. Below 20 K, a has linear-t dependence without saturation to a residual resistivity limit; and although there is a plateau in c for 5 T 35 K, it is followed by a very rapid downturn near T a =5 K Fig. 2 b inset, which indicates a sudden, rapid decrease in inelastic scattering. Oxygen depletion also changes the electronic density of states g E F of Sr 2 IrO 4, as reflected in the thermoelectric power S T, as shown in Fig. 2 c. A peak in the c axis S c T for =0.04 is only 1/3 of the peak value observed for =0. Since S T measures the voltage induced by a temperature gradient which cannot be confined to the surface of the ρ c ρ a ρ c (Ω cm) ρ a ~T 610-4 410-4 210-4 2.5 10-2 010 0 0 5 10 15 20 0 0 100 200 300 400 500 600 T (K) FIG. 2. Color online The a- and c-axis resistivities, a and c, as a function of temperature for a =0 and b 0.04; c c-axis thermoelectric power, S c T for =0 and 0.04. Inset in b : a and c vs T for 1.7 T 20 K; note a downturn near 5 K in c and linear temperature dependence in a. Will it be possible to dope Sr2IrO4 to higher doping levels? May be try field-effect doping or making oxide interface???? T a ρ c ρ a (Ω cm)
Summary Iridium oxides: spin-orbit coupling enhances importance of correlation effects even with moderate U. Sr2IrO4: spin-orbit entangled moments in a Mott insulator Hubbard model similar to cuprates but smaller overall energy scale. Doping: Possible high-tc superconductor, interesting lessons by comparing with cuprates. Other interesting iridate: Na2IrO3 - Jeff = 1/2 frustrated AF Mott insulator on honeycomb lattice (Singh, Gegenwart, PR B 2010; X. Liu,...Gegenwart, Hill, arxiv:1104.4046) Nature yet to be clarified? effects of doping?