Observation of possible topological in-gap surface states in the Kondo insulator SmB6 by photoemission

Abstract

SmB₆ , a well-known Kondo insulator, exhibits a transport anomaly at low temperature. This anomaly is usually attributed to states within the hybridization gap. Recent theoretical work and transport measurements suggest that these in-gap states could be ascribed to topological surface states, which would make SmB₆ the first realization of topological Kondo insulator. Here by performing angle-resolved photoemission spectroscopy experiments, we directly observe several dispersive states within the hybridization gap of SmB₆ . These states show negligible k_z dependence, which indicates their surface origin. Furthermore, we perform photoemission circular dichroism experiments, which suggest that the in-gap states possess chirality of the orbital angular momentum. These states vanish simultaneously with the hybridization gap at around 150 K. Together, these observations suggest the possible topological origin of the in-gap states.

ARTICLE

Received 7 Jul 2013 | Accepted 26 Nov 2013 | Published 18 Dec 2013

J. Jiang1,*, S. Li2,*, T. Zhang1,*, Z. Sun3,4, F. Chen5, Z.R. Ye1, M. Xu1, Q.Q. Ge1, S.Y. Tan1, X.H. Niu1,M. Xia1, B.P. Xie1, Y.F. Li2, X.H. Chen3,5, H.H. Wen2 & D.L. Feng1

SmB6, a well-known Kondo insulator, exhibits a transport anomaly at low temperature. This anomaly is usually attributed to states within the hybridization gap. Recent theoretical work and transport measurements suggest that these in-gap states could be ascribed to topological surface states, which would make SmB6 the rst realization of topological Kondo insulator. Here by performing angle-resolved photoemission spectroscopy experiments, we directly observe several dispersive states within the hybridization gap of SmB6. These states show negligible kz dependence, which indicates their surface origin. Furthermore, we perform photoemission circular dichroism experiments, which suggest that the in-gap states possess chirality of the orbital angular momentum. These states vanish simultaneously with the hybridization gap at around 150 K. Together, these observations suggest the possible topological origin of the in-gap states.

DOI: 10.1038/ncomms4010 OPEN

Observation of possible topological in-gap surface states in the Kondo insulator SmB6 by photoemission

1 State Key Laboratory of Surface Physics, Department of Physics, and Advanced Materials Laboratory, Fudan University, Shanghai 200433, China. 2 National Laboratory of Solid State Microstructures and Department of Physics, Nanjing University, Nanjing 210093, China. 3 Department of Physics, University of Science and Technology of China, Hefei 230026, China. 4 National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230029, China. 5 Hefei National Laboratory for Physical Science at Microscale, University of Science and Technology of China, Hefei 230026, China. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to D.L.F. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010

Kondo insulators (KIs) or heavy-fermion semiconductors13 are exotic materials with strong electron correlations, in which the localized 4f electrons give rise to novel ground

states. SmB6 is a well-known KI. At high temperature, the many-body interactions between the local spins of Sm ions and conduction electrons (the Kondo screening) are weak, and the system is a correlated metal. With decreasing temperature, an energy gap develops due to the hybridization of 4f bands and conduction 5d bands13, featured by a rapidly rising resistance. However, there is a long-standing puzzle: its resistivity does not diverge but saturates at temperatures below 4 K (refs 46). The residual conductivity was attributed to some in-gap states, evidenced by optical, neutron scattering, specic heat and angle-resolved photoemission spectroscopy (ARPES) experiments710. Although various models have been proposed to address the in-gap states5,11, their exact nature remains illusive after decades of studies.

Recently, the novel concept of topological KI (TKI) may shed new light on the in-gap states in SmB6 (refs 12,13). In the so-called topological insulator (TI) such as Bi2Se3 or Bi2Te3 (refs 14,15), there will be surface states that are protected by band topology from impurities that do not break time-reversal symmetry. Although the usual TIs are dened in weakly interacting systems, KIs with strong correlations can also be topologically classied. A topologically non-trivial KI could be realized by strong spin-orbit coupling associated with a hybridization band gap12,13. If SmB6 is a TKI as predicted12,13,16,17, there will be surface states that naturally explains the origin of in-gap states and the residual conductivity. In particular, unlike traditional TI materials in which surface transport was usually concealed by the defect-induced bulk carriers, SmB6 seems to possess a true insulating bulk that will make it a promising candidate material for future spintronics applications18. Several recent experiments have observed pronounced surface-dominating transport in SmB6 (refs 10,1821) suggesting the TKI scenario. However, direct electronic structure evidences for these potential topological surface states are highly demanded for the understanding of their exact nature.

In this letter we present an ARPES study of SmB6 single crystals. We observe several bands within the hybridization gap of SmB6, which show negligible kz dependence, indicative of their surface origin. Furthermore, the photoemission circular dichroism (CD) of these in-gap states suggests the chirality of the orbital angular momentum (OAM), and these states vanish simultaneously with the hybridization gap around 150 K. Our results thus strongly suggest their possible topological origin.

ResultsSample characterization and valence band structure. Sample preparation and experimental details are described in the Methods section. Figure 1a shows the powder X-ray diffraction pattern of SmB6 taken after grinding the single crystals into powder. No evident secondary phase is observed, showing the high quality of the crystals. The resistivity of our SmB6 sample exhibits a sharp increase below 50 K and starts to saturate below5 K (Fig. 1b), similar to previous reports46. As shown in Fig. 1b,c, magnetoresistance is clearly present, which is negative below B20 K and positive at higher temperatures. The negative magnetoresistance can be understood in the context of Kondo effect, where the Kondo singlet will be partially broken by the magnetic eld, which releases the localized d-band electrons and lowers the resistivity. In contrast, the positive magnetoresistance at high temperatures might be induced by the orbital scattering effect. Figure 1d shows the low-energy electron diffraction (LEED) pattern of the cleaved surface. Clean 1 1 pattern of

the (001) surface is observed. The bulk and projected Brillouin zones of the (001) surface are shown in Fig. 1e.

The photoemission intensities measured along GX and XM directions over a large energy scale are presented in Fig. 2a,b. Three at bands can be observed at 18, 150 and 950 meV below the Fermi energy (EF), as marked in the integrated energy distribution curves (EDCs). They are the 6H5/2, 6H7/2 and 6F

multiplets of the Sm2 4f6-4f5 nal state, respectively10,22. In addition, there is a highly dispersive band centred at X (dashed line). In previous ARPES studies22, it was ascribed to the bulk 5d band based on the band calculations23.

T (K)

14

0

Experiment Calculation Devitation

10

Intensity (countse+3)

12

001

011

0 T14 T

10

SmB6

8

111

002

021

R ()

6

211

1

0.1

4

022

003

031

311

222

032

321

2

20 30

40 50 60 70 80

90

1 10 100

2[afii9835]

k(001)

0

4

8

12

16

2 K

X

M

X

MR (%)

M

R

X

M

k(100)

M

X

k(010)

0 20 40 60 80 100 120

140

LEED pattern 120.2 eV @300 K

Figure 1 | SmB6 sample characterization. (a) Powder X-ray diffraction pattern of an SmB6 crystal. Solid line is the Rietveld tting using the TOPAS software. The s.d. is negligible. (b) The temperature dependence of resistance at zero eld and 14 T. A negative magnetoresistance is observed in the low temperature region, while it becomes positive in the high temperature region. (c) The magnetoresistance of SmB6 taken at 2 K. A negative magneto-resistance is obvious here. (d) LEED pattern of the cleaved SmB6 (001) surface. Bright spots in square lattice reect the pristine 1 1 surface, with a lattice

constant of 4.13 . (e) Bulk and surface projected Brillouin zone of SmB6 and the high symmetry points.

H (kOe)

2 NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010 ARTICLE

a

0.0

X X

0.5

4f flat bands

1.0

5d

EE F(eV)

1.5

2.0

7 K

7 K

21 eV

High

Low

b

M

X M

0.0

0.5

4f flat bands

Intensity (a.u.)

5d

EE F(eV)

1.0

1.5

2.0

27 eV

0.5 0.5

k// (1)

X

Figure 2 | Valence band structure of SmB6. (a,b) Left panel: large energy scale photoemission intensity plot along the GX and XM directions, respectively. Right panel: integrated EDCs to highlight the positions of the at 4f bands as marked by the arrows. Dashed lines in panels a and b indicate the Sm 5d band that crosses EF at high temperatures. Data were taken at 7 K with 21 and 27 eV photons, respectively, at Stanford Synchrotron Radiation Lightsource.

Dispersive in-gap states around C and X. To search for the ingap states in SmB6, the photoemission intensity at EF is collected over the surface projected Brillouin zone (Fig. 3a). Data were taken at 8 K, well into the KI phase. Indeed, four large oval-shaped Fermi surfaces are clearly observed around X, whose total area covers about 32.7% of the projected two-dimensional Brillouin zone, together with some spectral weight around G. In the second Brillouin zone, the Fermi surface contour around X is weak while the spectral weight around G is much enhanced probably because of the matrix element effects of the photo-emission process. The nature of these states can be further revealed in their dispersions. In Fig. 3b, we show the photo-emission intensity plot over a larger energy scale along cut #1 across X. The aforementioned 5d-like band is marked as d, and the band between the two at 4f bands is referred to as b. More details of the near-EF states and the hybridization gap are further illustrated in Fig. 3c. Here the hybridization gap is manifested by the suppressed spectral weight between EF and the 6H5/2 4f band at 18 meV (referred to as f). However, there is an in-gap

dispersive band centred at X (referred to as a), which gives the large oval-shaped Fermi surface as shown in Fig. 3a. The dispersion of a is clearly visible in the momentum distribution curves (MDCs) in Fig. 3d as well, with Fermi crossings at

0.29 1. The band becomes obscure when it crosses f and disappears below 30 meV. The photoemission intensity along a

cut across G is shown in Fig. 3e, where the weak feature located at G exhibits an electron-like band dispersion as further illustrated by their MDCs in Fig. 3f. This is more readily visible in the second Brillouin zone as shown in Fig. 3g,h. This band (referred to as g) gives the tiny central electron pocket with a Fermi momentum of B0.09 1 in Fig. 3a. Its area is just about 1.1% of the projected two-dimensional Brillouin zone. The bright spot on band f at G could be attributed to its crossing with g, and the crossing of the two g branches is estimated to occur at 233 meV. Furthermore, an oval-shaped Fermi surface a0 can be observed around G in the 2nd Brillouin zone, which shows almost the same size as the a-Fermi surface around X, although it appears weaker in the 1st Brillouin zone. Its twofold symmetry and size indicate that it is most likely to be the shadow Fermi surface from the umklapp scattering of the a-Fermi surface around X owing to the existence of a 1 2 surface reconstruction

in this sample as reported in ref. 10 before. However, the superstructure spot of such a reconstruction is missing in the LEED pattern of our samples (Fig. 1d). Probably, it is too weak for the sensitivity of our LEED apparatus.

kz dependence of the electronic structure. Since the bulk SmB6 is an insulator based on the transport data and calculations, the ingap states a and g should be some metallic surface states. To further illustrate this, we varied photon energies to reveal the kz dependence of the electronic structure. The photoemission intensity distribution across X taken with various photon energies are presented in Fig. 4a. The photon energies we used cover a full kz period in the extended Brillouin zone (Fig. 4e). The intensity distributions do not show any noticeable kz dependence. In the corresponding MDCs shown in Fig. 4b, we tracked the dispersion of the a band taken with 25 eV photons and overlaid it on the data taken with other photon energies, which clearly show that the dispersions and Fermi crossings of a are kz independent. This behaviour indicates the surface origin of the a pockets. However, in addition, one could observe a shoulder feature that does show some photon energy dependence, which thus might be contributed by a bulk state. Furthermore, Fig. 4c shows the photo-emission intensity distributions taken with various photon energies over a large energy scale across X. Since kx, ky and kz are equivalent in a cubic compound, the bulk 5d bands are expected to have similar fast dispersions along kz as those in kx and ky directions. However, at the rst sight, the 5d-like fast dispersive bands together with others do not show strong kz dependence, for one can notice that the envelope of the photoemission intensity in Fig. 4c appears to be independent of the photon energies. Nevertheless, the bright interior of the intensity envelope for 29 and 31 eV data is a typical consequence of the projection of intensities from different kz due to the poor kz resolution of

ARPES24. As is further illustrated by the broad MDCs taken at

500 meV around X in Fig. 4d, the nite kz resolution of ARPES spectral weight smears out the intensity distribution for the 5d band with a strong kz dispersion. Combining these observations, we could conclude that the a band probably contains both surface and bulk contributions, while the b and d bands are most probably originated from the bulk 5d band.

CD of the electronic structure. To study whether the surface states are topologically non-trivial or not, the chirality of their spin (S) and OAM shall be examined. This is another remarkable hallmark of the topological surface state: both the spin and OAM are interlocked and rotate with the electron momentum (k) and, their sum, the total angular momentum J, is a good quantum

NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010

a

c

e

g

25 eV @ 8 K

2

21.2 eV

21.2 eV

25 eV

[afii9825]

0

2nd BZ

[afii9825]

[afii9825]

Cut #3

[afii9825]

[afii9825]

X

1st BZ

Cut #2

k y(/a)

1

1

0

EE F(meV)

[afii9828]

[afii9828]

[afii9850]

20

40

X

Cut #4

M

Cut #1

[afii9828]

[afii9826]

High Low

Cut #1

12 K

Cut #2

25 K

Cut #3

8 K

1

0

1

kx(/a)

b d

f h

21 eV

0

100

300

EF

EF

EF

[afii9825]

[afii9825]

[afii9826]

[afii9825]

Intensity (a.u.)

20 meV

40 meV

[afii9828]

[afii9828]

EE F(meV)

200

20 meV

20 meV

[afii9826]

High

Low

[afii9829]

400

Cut #1

9 K

0.4 X

0.4

0.4 X

0.4

0.4

0.4

0.4

0.4

k// (1)

k// (1)

Figure 3 | Dispersive in-gap states around C and X. (a) Photoemission intensity map at the Fermi energy for SmB6 taken with 25 eV photons at Hiroshima Synchrotron Radiation Center. The intensity was integrated over a window of (EF 5 meV, EF 5 meV). The left panel is a sketch of the projected

two-dimensional Brillouin zone. (b) Photoemission intensity plot over a large energy scale along cut #1 taken with 21 eV photons at Stanford synchrotron Radiation Lightsource. (c,d) The photoemission intensity plot and momentum distribution curves (MDCs) along cut #1 in the projected two-dimensional Brillouin zone as indicated in panel a, respectively. (e,f) The photoemission intensity plot and MDCs along cut #2, respectively. (g,h) The photo-emission intensity plot and MDCs along cut #3, respectively. The dashed lines in panels bh indicate the dispersions of various bands, ignoring band warping at crossings due to hybridization. The photon energies and temperatures for measurements are labelled in individual panels. The a and g Fermi surfaces cover about 32.7 and 1.1% of projected two-dimensional Brillouin zone, respectively. The Fermi velocities of a and g are B0.24 and B0.22 eV , respectively; the linear extrapolations of the a and g band dispersions cross at 654 and 233 meV below EF, respectively.

number. Spin-resolved ARPES requires long acquisition time because of its low count rate and, at this stage, it is not feasible for the weak surface state signal discussed above. Alternatively, instead of examining the chirality of the spin, one could examine the chirality of the OAM with the so-called CD of ARPES, namely, the difference of photoemission intensities taken with right circularly polarized (RCP) light and left circularly polarized (LCP) light. For Bi2Se3, it has been shown that photoemission CD signal is proportional to the inner product between OAM and light propagation vector25. This technique has been demonstrated to be a powerful tool to investigate the OAM texture of the surface states in TIs25,26.

Figures 5a and b show two Fermi surface maps taken with 25 eV RCP light and LCP light, respectively. Comparing the two data sets, we can clearly see the difference between the two intensity maps taken with different circular polarizations. The intensity of the upper electron pocket on the positive ky side is higher than that in the lower electron pocket on the negative ky side for the RCP data, while the LCP data exhibit an opposite behaviour. The electron pocket at the right side also shows the switching of high intensity regions between the positive ky side

and negative ky side with different circularly polarized light. However, the electron pocket at the left side does not show strong CD, which is probably because of the asymmetric matrix element effects. Since the incident light is at 50 angle from kz in the kxkz plane, the emission of the photoelectrons on the left and right Fermi surface sheets are asymmetric with respect to the light propagation direction. The CD can be more clearly represented by the normalized difference (RCP LCP)/(RCP LCP) in Fig. 5c.

Figure 5df shows the similar results for another SmB6 x sample

measured with 35 eV photons. More specically, Fig. 5g plots the CD values in Fig. 5c,f along the right a-Fermi pocket, which could be tted by individual sinusoidal functions. The amplitude of CD is thus found to be photon energy dependent.

As for the g band at G, its weak intensity makes the CD effects not so obvious. However, for the data taken at 25 eV, Fig. 5h shows that there apparently exists certain intensity inversion between the RCP and LCP data in the MDCs across G along the ky axis. We can see that higher intensities appear on the negative ky side for the RCP data but on the positive ky side for the LCP data. Therefore, both a and g bands show similar CD that is antisymmetric with respect to the ky 0 axis.

4 NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010 ARTICLE

a

25 eV

30 eV

35 eV

40 eV

45 eV

50 eV

75 eV

0

d

[afii9825]

[afii9825]

[afii9826]

[afii9825]

[afii9826]

[afii9825]

[afii9825]

[afii9825]

EF500 meV

34 eV

20

[afii9825]

[afii9826]

EE F(eV) E E F (meV)

40

31 eV

Intensity (a.u.)

[afii9826]

[afii9826]

[afii9826]

[afii9826]

29 eV

27 eV

25 eV

23 eV

b

21 eV

19 eV

Intensity (a.u.)

17 eV

0.4 0.4

X

20 meV

k// (1)

e

6

5

4

25 eV

30 eV

35 eV

40 eV

45 eV

50 eV

75 eV

X

75 eV

50 eV

X 45 eV

40 eV

35 eV 25 eV

30 eV

0.3 0.3

X 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3

25 eV

X X X X X X

0.4 0.4

X 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4

X X X X X X

EF

c

0

k z(/c)

M

0.5

X

1

17 eV

19 eV

1.5

3

2

High

Low

0

1

19 eV

23 eV

27 eV

29 eV

31 eV

34 eV

ky(/b)

k// (1)

Figure 4 | Photon energy dependence of the in-gap state around X. (a) Photon energy dependence of the photoemission intensity measured along cut #4 shown in Fig. 2a. The a and b band dispersions are indicated by the dashed lines. The bright regions correspond to the location where the 4f and 5d bands hybridize, which are particularly strong at some photon energies due to matrix element effects. Data were taken at 1 K at BESSY II. (b) The corresponding MDC plots of panel a. The red dashed lines are dispersions tracked from the MDC peak in the 25-eV data, and they are overlaid on data taken at other photon energies, indicating the lack of photon energy dependence. The blue circles refer to the shoulders in the MDCs that might be contributed by the bulk band. (c) The photoemission intensity measured along cut #1 in Fig. 2a for various photon energies. Data were taken at 7 K at Stanford Synchrotron Radiation Lightsource. (d) MDCs at 500 meV for data in panel c. (e) The sampled momentum cuts for various photon energies

calculated with the inner potential 14 eV following ref. 10, indicating that the photon range that we used cover an extended Brillouin zone in kz.

Many CD behaviours observed here resemble those of the nontrivial surface state found in Bi2Se3 before25,26, which indicates that both a and g bands are likely to be topologically non-trivial. However, we notice that the CD behaviour of the a band is much more complicated than that of Bi2Se3. For example, the CD for the left a pocket is much weaker than the others, and the upper and lower a pockets exhibit different CD behaviours from that of the right a pocket. Particularly, the CD for both upper and lower a pockets does not change sign across the Brillouin zone boundary in Fig. 5c,f. In addition to certain nal-state effect27, the complicated CD might in part be due to the coexistence of bulk and surface components in the a pocket. For example, if the bulk 5d state and the surface state on the a pocket have different CD and their photoemission intensities vary differently with the emission angle and photon energy, one may expect asymmetric CD on left and right a pockets, and the photon energy-dependent

CD amplitude shown in Fig. 5g. To fully understand these CDARPES data, such as the different functional forms of CD in the upper and right pockets, detailed density functional theory calculations of the OAM texture are required for future further studies, especially when there might be out-of-plane OAM component. Nevertheless, the pattern in our CD-ARPES data does suggest a non-trivial OAM texture (thus, also spin texture through spin-orbit coupling) for the a and g bands, which is a

necessary manifestation in photoemission if they are topologically non-trivial.

Temperature dependence of the electronic structure. If a and g are topological surface states of a topological KI, they should vanish when the hybridization gap closes at high temperatures. Figure 6a presents temperature-dependent ARPES intensity measured across X. a stays almost unchanged at To80 K, but its band velocity starts to increase above 80 K. Eventually, it appears that a and b merge into one highly dispersive band at 150 K, which ts the large-scale band dispersion shown in Fig. 2b. Therefore, a, b and d could be the broken sections of the 5d band induced by the 4f5d hybridization. However, one also sees that a goes straight through f (Fig. 3c,d), which suggests that the a-Fermi surface probably contains contributions both from the not-fully gapped bulk 5d band and from the topological surface state, although the latter disappears before ending at a Dirac point. This is also consistent with the conclusion that we draw from the photon energy dependence data.

The temperature evolution of the hybridization can be further demonstrated by the integrated EDCs of the images in Fig. 6a. After some processing described in the caption of Fig. 6, the resulting EDCs are presented in Fig. 6b. As temperature increases,

NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010

a b

RCP LCP

c

(RCPLCP)/(RCP+LCP)

g

40

25 eV @ Hisor

M

X

25 eV

M

X

25 eV

M

X

25 eV

0

35 eV @ SLS

1 )

k y(

40

0 2 (rad)

/2

3/2

d e

RCP LCP (RCPLCP)/(RCP+LCP)

RCPLCP

RCP+LCP(%)

f

h

RCP

LCP

M

X

35 eV

M

X

35 eV

M

X

35 eV

Intensity (a.u.)

kx=01

1 )

+

[0.081, 0.08

k y(

1] RCP

LCP

0

0.8 0.4 0 0.4 0.8 ky (1)

kx (1)

Figure 5 | CD of the surface states at X and C. (a,b) Fermi surface maps of SmB6 taken with right circularly polarized (RCP) and left circularly polarized (LCP) light, respectively. The intensity was integrated over a window of ( 5 meV, 5 meV). The data were taken at 8 K with 25 eV photons at Hiroshima

Synchrotron Radiation Center (HiSOR). (c) The differential map of the RCP and LCP photoemission intensities in panels a and b that are normalizedto their sum intensity. (df) The same as in panels ac, but for another SmB6 x sample with a slight boron deciency. The intensity was integrated over a

window of ( 10 meV, 10 meV) and the data were taken at 27 K with 35 eV photons at Swiss Light Source. Note the g band is not so visible here,

probably because of its negligible matrix element at this photon energy. (g) The CD values taken along the right pocket in panel c (red dots) and in panel f (blue dots), where the polar angle O is dened within panel c with respect to X. The dashed curves are sine function ts of the experimental data. (h) Top: the two MDCs along ky taken at EF and kx 0 with RCP and LCP light, respectively, where the CD of the g band are indicated by the arrows. Bottom:

analogous to the top except that the data are integrated over a kx window of ( 0.08 1, 0.08 1) to cover the g pocket. The data were taken at 8 K with

25 eV photons at HiSOR.

a

b

12 K

150 K

20 K

12 K

[afii9825]

[afii9826]

[afii9850]

1.4

1

0.6

40 K

60 K

80 K

100 K

150 K

0

Intensity (a.u.)

150 K

20

100 K

80 K

5d

60 K

40

40 K

20 K

12 K

EE F(meV)

40 20 0

0.4 X

0.4

0.4

X 0.4

0.4

0.4

X 0.4

0.4

X 0.4

0.4

0.4 X

0.4

X 0.4

0.4

X

EEF (meV)

c

d e

E 10meV

1.0

0.8

0.6

0.4

0.2

0.0

12 K 25 K 45 K 80 K 150 K

0.4 0.4

0.4 0.4

0.4 0.4

0.4 0.4

0.4 0.4

0

[afii9828]

20

Intensity (a.u.)

150 K80 K45 K25 K12 K

Intensity (a.u.)

[afii9825] [afii9828]

[afii9850]

High

Low

5d

[afii9828]

40

0.4 0.4

0

50

100

150

Temperature (K)

k// (1) k// (1)

Figure 6 | Temperature dependence of the photoemission data of the surface state around X. (a) Photoemission intensity of SmB6 along cut #1in Fig. 2a at 12, 20, 40, 60, 80, 100 and 150 K. (b) Temperature dependence of the angle-integrated spectra of panel a, the Fermi cutoff effects were removed by dividing the angle-integrated spectra by the energy-resolution convoluted FermiDirac distribution functions at individual temperatures.(c) Photoemission intensity of SmB6 along cut #2 in Fig. 2a at 12, 25, 45, 80 and 150 K. (d) Temperature dependence of the MDCs at 10 meV below EF for data in panel c, and the shaded regions represent the difference from the 150-K data. (e) Area of the shaded region normalized by the 12-K area. The error bar is the relative uncertainty of the integrated photoelectron count. All data were taken with 21.2 eV photons.

6 NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010 ARTICLE

the peak near 20 meV (derived from f) is gradually

suppressed, and the dip near EF (a measure of the hybridization strength) is gradually lled. This indicates that the 4f states become localized and completely lose coherence at TE150 K.

Consequently, they are decoupled with the 5d bands and the hybridization disappears28,29. Similarly, in the photoemission intensity near G plotted in Fig. 6c, the g band also vanishes with the hybridization gap at 150 K, which can be further visualized by the MDCs in Fig. 6d, and the integrated intensity in Fig. 6e as well. The thermal broadening effects from 12 to 80 K are not obvious in Fig. 6d; thus, the disappearance of the intensity at 150 K should be intrinsic. Therefore, the temperature dependencies of the surface states are consistent with the behaviours expected for a TKI.

DiscussionSince the (001) surface Brillouin zone contains one G and two X points on the zone boundary, we have two a bands enclosing two surface X points but only one g band around G. Regardless of whether or not a is a topologically trivial surface state, there will be totally an odd number of branches of surface states, which indicates that the whole system is topologically non-trivial. In fact, the basic electronic structure of these surface states qualitatively agrees with those topological surface states predicted in refs 16,17 remarkably well, although there are some complications in the CD that are still to be understood. In contrast, we note that unlike the case of weakly correlated systems, where theoretical calculations work quantitatively well for describing TIs14,15, there are some quantitative differences between our data and the calculations of SmB6. For example, the

Fermi crossings of a and g are different from the calculations, and the observed a Fermi surface has ve to six times area of the calculated ones. The discrepancies might be partially because the chemical potentials used in theories are different from that in the real surface, and/or the actual cleaved surface is different from the ideal surface used in the calculations. Furthermore, these discrepancies also demonstrate that correlation might affect the topological surface states, since different calculations could vary somewhat depending on how the correlations were treated16,17. Finally, it is also intriguing to note that the Fermi momenta of the two surface states on the (101) surface are 0.0383 and0.0955 1, respectively as found in a recent quantum oscillation measurements of SmB6 (ref. 30), which are rather different from our measurements of the (001) surface (0.09 1 for the small pocket near G and 0.29 1 along the short axis of the large oval Fermi surface).

In summary, we present a direct observation of the in-gap surface states in SmB6. The CD, and photon energy and temperature dependences of these in-gap states were systematically studied, which exhibit characteristic behaviours of a topological surface state for a KI, and thus support that SmB6 is a possible TKI. Our ndings lay the foundation for understanding the anomalous transport properties of SmB6 and further exploration of the interplay between strong correlations and topological effects in TKIs in general.

Note added. After nishing this manuscript, we note that there are other three pieces of independent ARPES work on SmB6 (refs 3133). Two of them report similar results31,32, except that our temperature dependence data agree with those reported in ref. 31, while ref. 32 suggests the surface states vanish at a much lower temperature of 30 K. The size of the a pocket measured here is similar to that reported in ref. 31 but much smaller than that reported in ref. 32. These discrepancies may be caused by different samples and requires further investigation. On the basis of data

taken at 38 K, ref. 33 claims that the in-gap states around X must be bulk states. This partially agrees with our ndings. However, the detailed data here taken at lower temperatures do indicate that the a band contains both surface and bulk components.

Methods

Sample synthesis. High-quality single crystals of SmB6 were grown either with the arc-melting technique or the Al-ux method. The arc melting technique was conducted by establishing a temperature gradient in the slow cooling process. First, the metal Sm were cut into small pieces and mixed with the boron powder in the ratio of Sm:B 1:6. The mixture was pressed into a pellet with the diameter of 1 cm

and thickness of about 0.5 cm. Then the pellet was heated up in the water-cooled clean oven for arc melting. Welding was repeated for ve times to ensure the sample uniformity. In the Al-ux method, a chunk of Sm (99.9%) together with the powders of Boron (99.99%) and Al (99.99%) were mixed with a ratio of 1:6:400. Then the mixture was loaded into an alumina crucible. The entire mixture was heated to 1,773 K and then maintained at 1,773 K for about 2 days before slowly cooling it down to 873 K at 5 K h 1. During all the heating progress, the mixture was kept in the circumstance with owing Argon gas. The samples with Al ux were soaked in dense NaOH solution, and then washed with dilute HNO3 solution.

Shiny crystals in the sizes of millimetre or sub-millimetre were obtained.

ARPES measurement. The samples were cleaved along the (001) plane and measured under ultra-high vacuum of 3 10 11 torr. The in-house ARPES

measurements were performed with SPECS UVLS discharge lamp (21.2 eV He-Ia light). The synchrotron ARPES experiments were performed at the Beamline 5-4 of Stanford Synchrotron Radiation Lightsource and the one-cubed ARPES station of BESSY II. The CD experiments were performed at Beamline 9 A of Hiroshima Synchrotron Radiation Center and the Surface and Interface Spectroscopy (SIS) beamline of Swiss Light Source. All data were taken with Scienta electron analysers, the overall energy resolution was better than 7 meV and the angular resolution was0.3.

References

1. Aeppli, G. & Fisk, Z. Kondo insulators. Comm. Condens. Matter Phys. 16, 155 (1992).

2. Riseborough, P. Heavy fermion semiconductors. Adv. Phys. 49, 257320 (2000).3. Coleman, P. Heavy Fermions: Electrons at the Edge of Magnetism, Handbook of Magnetism and Advanced Magnetic Materials, Vol 1, 95148 (Wiley, 2007).

4. Menth, A., Buehler, E. & Geballe, T. H. Magnetic and semiconducting properties of SmB6. Phys. Rev. Lett. 22, 295 (1969).

5. Allen, J. W., Batlogg, B. & Wachter, P. Large low-temperature Hall effect and resistivity in mixed-valent SmB6. Phys. Rev. B 20, 48074813 (1979).

6. Cooley, J. C., Aronson, M. C., Fisk, Z. & Caneld, P. C. SmB6: Kondo insulator or exotic metal? Phys. Rev. Lett. 74, 1629 (1995).

7. Nanba, T. et al. Gap state of SmB6. Physica B 186188, 440443 (1993).8. Nyhus, P., Cooper, S. L., Fisk, Z. & Sarrao, J. Low-energy excitations of the correlation-gap insulator SmB6: a light-scattering study. Phys. Rev. B 55, 12488 (1997).

9. Alekseev, P. A. et al. Magnetic excitations in SmB6 single crystals. Physica B 186188, 384386 (1993).

10. Miyazaki, H. et al. Momentum-dependent hybridization gap and dispersive ingap state of the Kondo semiconductor SmB6. Phys. Rev. B 86, 075105 (2012).

11. Curnoe, S. & Kikoin, K. A. Electron self-trapping in intermediate-valent SmB6.

Phys. Rev. B 61, 15714 (2000).12. Dzero, M. et al. Topological Kondo insulators. Phys. Rev. Lett. 104, 106408 (2010).

13. Dzero, M. et al. Theory of topological Kondo insulators. Phys. Rev. B 85, 045130 (2012).

14. Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045 (2010).

15. Qi, X. L. & Zhang, S. C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057 (2011).

16. Takimoto, T. SmB6: a promising candidate for a topological insulator. J. Phys. Soc. Jpn 80, 123710 (2011).

17. Lu, F. et al. Correlated topological insulators with mixed valence. Phys. Rev. Lett. 110, 096401 (2013).

18. Kim, D. J., Grant, T. & Fisk, Z. Limit cycle and anomalous capacitance in the Kondo insulator SmB6. Phys. Rev. Lett. 109, 096601 (2012).

19. Wolgast, S. et al. Low temperature surface conduction in the Kondo insulator SmB6. Preprint at http://arxiv.org/abs/1211.5104

Web End =http://arxiv.org/abs/1211.5104 (2012).

20. Zhang, X. et al. Hybridization, inter-ion correlation, and surface states in the Kondo insulator SmB6. Phys. Rev. X 3, 011011 (2013).

21. Botimer, J. et al. Robust surface Hall effect and nonlocal transport in SmB6:

indication for an ideal topological insulator. Preprint at http://arxiv.org/abs/ 1211.6769

Web End =http://arxiv.org/abs/ http://arxiv.org/abs/ 1211.6769

Web End =1211.6769 (2012).

NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4010

22. Denlinger, J. D. et al. Advances in photoemission spectroscopy of f-electron materials. Physica B 281&282, 716722 (2000).

23. Massidda, S. et al. Electronic structure of divalent hexaborides. Z. Phys. B 102, 8389 (1997).

24. Ye, Z. R. et al. Orbital selective correlations between nesting/scattering/Lifshitz transition and the superconductivity in AFe1-xCoxAs (A Li, Na). Preprint at

http://arxiv.org/abs/1303.0682

Web End =http://arxiv.org/abs/1303.0682 (2013).25. Park, S. R. et al. Chiral orbital-angular momentum in the surface state of Bi2Se3. Phys. Rev. Lett. 108, 046805 (2012).

26. Wang, Y. H. et al. Observation of a warped helical spin texture in Bi2Se3 from circular dichroism angle-resolved photoemission spectroscopy. Phys. Rev. Lett. 107, 207602 (2011).

27. Scholz, M. R. et al. Reversal of the circular dichroism in angle-resolved photoemission from Bi2Te3. Phys. Rev. Lett. 110, 216801 (2013).

28. Nozawa, S. et al. Ultrahigh-resolution and angle-resolved photoemission study of SmB6. J. Phys. Chem. Solid 63, 12231226 (2002).

29. Souma, S. et al. Direct observation of pseudogap of SmB6 using ultrahigh resolution photoemission spectroscopy. Physica B 312-313, 329330 (2002).

30. Li, G. et al. Quantum oscillations in Kondo Insulator SmB6. Preprint at http://arxiv.org/abs/1306.5221

Web End =http://arxiv.org/abs/1306.5221 (2013).

31. Xu, N. et al. Surface and bulk electronic structure of the strongly correlated system SmB6 and implications for a topological Kondo insulator. Phys. Rev. B 88, 121102(R) (2013).

32. Neupane, M. et al. Surface electronic structure of the topological Kondo Insulator candidate correlated electron system SmB6. Nat. Commun 4, 2991 (2013).

33. Frantzeskakis, E. et al. Kondo hybridization and the origin of metallic states at the (001) surface of SmB6. Preprint at http://arxiv.org/abs/1308.0151

Web End =http://arxiv.org/abs/1308.0151 (2013).

Acknowledgements

We are grateful for the discussions with Professor R.B. Tao and Professor G.M. Zhang, and the experimental support by Dr D.H. Lu, Dr Y. Zhang and Dr M. Hashimoto at

Stanford Synchrotron Radiation Lightsource (SSRL), Dr E. Rienks at BESSY II,Dr M. Arita and Professor Shimada at Hiroshima Synchrotron Radiation Center, and Dr N. Plumb and Professor M. Shi at Swiss Light Source. This work is supported in part by the National Science Foundation of China and National Basic Research Program of China (973 Program) under the Grant numbers 2012CB921400, 2011CB921802, 2011CBA00100, 2012CB821403 and 2012CB21400. SSRL is operated by the US DOE, Ofce of Basic Energy Science, Divisions of Chemical Sciences and Material Sciences. This work is conducted within the collaborative network of National Center of Microstructures and Quantum Manipulation.

Author contributions

J.J., Z.R.Y., Z.S., S.Y.T., M.Xu, Q.Q.G., X.H.N., M.Xia and B.P.X. performed ARPES measurements. S.L. and Y.F.L. grew most of the samples and conducted sample characterization measurements, F. C. and X. H. C. grew the samples for the CD-ARPES measurement, J.J., T.Z. and D.L.F. analysed the ARPES data, T.Z., D.L.F., J.J., Z. S. and H.H.W wrote the paper. D.L.F. and H.H.W. are responsible for the infrastructure, project direction and planning.

Additional information

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsand permissions

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsand permissions

Web End =reprintsandpermissions /

How to cite this article: Jiang, J. et al. Observation of possible topological in-gap surface states in the Kondo insulator SmB6 by photoemission. Nat. Commun. 4:3010 doi: 10.1038/ncomms4010 (2013).

This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. To view a copy of this licence visit http://creativecommons.org/licenses/by/3.0/

Web End =http:// http://creativecommons.org/licenses/by/3.0/

Web End =creativecommons.org/licenses/by/3.0/ .

8 NATURE COMMUNICATIONS | 4:3010 | DOI: 10.1038/ncomms4010 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.