# Topological Kondo Insulators

###### Abstract

Kondo insulators are particularly simple type of heavy electron material, where a filled band of heavy quasiparticles gives rise to a narrow band insulator. Starting with the Anderson lattice Hamiltonian, we develop a topological classification of emergent band structures for Kondo insulators and show that these materials may host three-dimensional topological insulating phases. We propose a general and practical prescription of calculating the strong and weak topological insulator indices for various lattice structures. Experimental implications of the topological Kondo insulating behavior are discussed.

###### pacs:

71.27.+a, 75.20.Hr, 74.50.+rKondo insulators are a particularly simple type of heavy fermion material, first discovered forty years agoGeballe69 (), in which highly renormalized f-electrons, hybridized with conduction electrons, form a completely filled band of quasiparticles with excitation gaps in the millivolt range Coleman2007 (); KIReviews (). While these materials are strongly interacting electron systems, their excitations and their ground-states can be regarded as adiabatically connected to non-interacting band-insulatorsallen ().

It was recently shown that time-reversal invariant band insulators can be classified by the topological structure of their ground-state wavefunctions Fu2007 (); Moore2007 (); Roy2009 (); Qi2008 (). One of the dramatic consequences of this discovery, is the existence of a new class of “topological” band-insulator in which strong spin-orbit coupling leads to a ground-state that is topologically distinct from the vacuum, giving rise to gapless surface excitations.

In this paper, we show that Kondo insulators, as adiabatic descendents of band insulators, can also be topologically classified. The strong spin-orbit coupling characteristic of these materials leads us to predict that a subset of Kondo insulators are topologically non-trivial, with anomalous surface excitations. In current models of topological insulators, the spin-orbit coupling is encoded in a spin-dependent hopping amplitudes between different unit cells. By contrast, in a topological Kondo insulator (TKI), we show that the topologically nontrivial insulating state is produced by the spin-orbit coupling associated with the hybridization between conduction and -electrons.

Below, we develop a simple model for topological KIs. The physics we study is motivated by the canonical Kondo insulating behavior of SmBGeballe69 () and CeBiPt Hundley1990 (). The realization of a particular topologically nontrivial insulating state depends on the position of renormalized -level relative to the bottom of the conduction band (see Fig. 1). To analyze the topology of the bands in these materials, we use a periodic Anderson lattice model.

In a KI, the insulating state arises due to hybridization between the conduction and -electrons, provided that the chemical potential lies inside the hybridization gap separating the quasiparticle bands. The spatial symmetry of the hybridization amplitude is determined by the symmetry of the underlying crystal-field Kramers doublets of the rare-earth ion and it is precisely this symmetry that is responsible for non-trivial topological structures in a KI. To analyze this topology, we first employ a tight-binding model on a simple cubic lattice, which is adiabatically connected to the Hamiltonian of the KI material. We show that in this case, band topology is uniquely determined by the noninteracting band structure of the system in the absence of hybridization. Secondly, we consider a more general KI in a lattice with a body centered cubic lattice, and show that regardless of microscopic details, there always exists a parameter range in which a KI is a strong topological insulator (STI).

We begin with the periodic Anderson Kondo lattice Hamiltonian, written in terms of the fermion operators associated with the crystal field symmetry of the underlying lattice

(1) |

where is the dispersion of a tight-binding band of conduction electrons. We assume that the ground state of the isolated magnetic ion is a Kramers doublet , where labels a particular representation of the crystal symmetry group. For instance, in a cerium-based Kondo system, the Ce ion has a Ce valence configuration, hence we have one -electron in a atomic shell. In Eq. (1) above, the operator creates an electron on site in a Wannier state described by the quantum number , is the bare energy of the -level, is the bare hybridization, and describes Hubbard repulsion between -electrons. One can relate the Wannier states at site as follows Coqblin1969 (); Flint2008 (): , where the form factors are two dimensional matrices

(2) |

where

(3) |

is a tight-binding generalization of the spherical Harmonics that preserves the translational symmetry of the hybridization, , where is reciprocal lattice vector. Here, are the positions of the Z nearest neighbor sites around the magnetic ion.

The low-energy properties of the model (1) are described in terms of renormalized quasiparticles formed via strong hybridization between the and states and on-site repulsion . In the regime where the states are predominantly localized, we can neglect the momentum dependence of the -electron self-energy so that the effective low-energy Hamiltonian reads Ikeda1996 ():

(4) |

where is the bare spectrum of conduction electrons taken relative to the chemical potential, is the renormalized -level, , , and denotes the unit matrix. The KI is formed if the chemical potential of the quasiparticles lies inside the hybridization gap, separating the two bands with the spectra , with .

From Eq. (2), we see that the form factors are momentum-dependent unitary matrices that relate the spin quantization axes of the Bloch states and the spin-orbit coupled Wannier states. Our choice of hybridization ensures that the mean-field Hamiltonian (Eq. 4) is a periodic function satisfying .

These form factors are uniquely determined by the wave-functions of a magnetic ion, . For example, in a tetragonal crystal field environment, the Ce multiplet is split into three doublets: , , and , where the mixing angle defines orientation of the corresponding states. In an orthorhombic environment, the Kramer’s doublets are generally described by a linear superposition of all three wave-functions Kikoin1994 (); Moreno2000 (). To discuss the topological properties of these KI (4), we need to consider separately the form factors for different ’s. We distinguish these states according to their orbital symmetry parameterized by the index and the pseudo-spin quantum number () Moreno2000 (). Hence, we have , , and .

The momentum-dependence of the hybridization gap follows from Eq. (2). At small momenta , , , and , where and define the direction of the unit vector , associated with the point on the Fermi surface. Note that the hybridization gap has a line of nodes along the -axis for the shapes , but generic combinations of all three form-factors characteristic of contain no nodes. The key results of this paper are most simply illustrated using the nodeless Kramers doublet as the ground-state of the magnetic ion.

To analyze the topology of the bands we use the fact that topology is invariant under any adiabatic deformation of the Hamiltonian. We begin our study with a tight-binding model for a KI on a simple cubic lattice. The technical analysis is readily generalized to more complicated cases as discussed below. The most important element of the analysis is the odd parity form factor of the electrons, . This parity property is the only essential input as far as the topological structure is concerned.

In Ref. FuKane2007 (), Fu and Kane demonstrate that in an insulator with time-reversal and space-inversion symmetry, the topological structure is determined by parity properties at the eight high-symmetry points, , in the 3D BZ which are invariant under time-reversal, up to a reciprocal lattice vector: (see insets in Fig. 1). In our case, these symmetries require that and , where the parity matrix and the unitary part of the time-reversal operator are given by

(5) |

where is the second Pauli matrix. For any space-inversion-odd form factor, it follows immediately that at a high-symmetry point. Hence, the Hamiltonian at this high symmetry point is simply , where is the four-dimensional identity matrix.

The parity at a high symmetry point is thus determined by . Four independent topological indices Kitaev () (one strong and three weak indices) can be constructed from : (i) The strong topological index is the product of all eight ’s: ; (ii) by setting (where ), three high-symmetry planes, , are formed that contain four high-symmetry points each. The product of the parities at these four points defines the corresponding weak-topological index, , . The existence of the three weak topological indices in 3D is related to a topological index for 2D systems (a weak 3D TI is similar to a stack of 2D topological insulators). Because there are three independent ways to stack 2D layers to form a 3D system, the number of independent weak topological indices is also three.

A conventional band insulator has all of the four indices , while an index indicates a topological state with the odd number of surface Dirac modes. For a KI with and , we find , and hence the Kondo insulating state is a strong-topological insulator which will be robust with respect to disorder. Weak-topological insulators and topologically trivial insulators can in principle be found for different band structures and different values of , Fig. 1. A particularly interesting possibility is to tune topological phase transitions between different types of insulators (e.g., by applying a pressure), Fig. 1. Although we have been specifically considering a tight-binding model with a primitive unit cell, all our conclusions apply directly to systems adiabatically connected to this model. In order to prove it explicitly and to investigate more general cases, we develop a different and more general technique similar to that proposed in Ref. Roy2009 ().

We shall now study an example of a KI for the specific shape and explore the parameter range in which it remains a STI. Here, the form factor is universally determined in the small-momentum limit by the -wave symmetry of the electron orbitals and the point group symmetry of the lattice. Expression for the form factor of a Kramer’s doublet is

(6) |

where and are the spherical harmonics. At larger momenta, the form-factor depends on the microscopic details of the lattice and the Kondo coupling. In general form factor can be written as .

We show below that momenta where are crucial in calculating the topological indices, so that our results are generic for a linear combination of all three shapes. The most obvious zero-point is located in the origin as shown in Eq. (6). We now prove that this zero-point is topologically protected and that its existence necessarily yields the existence of other zero-points with (this conclusion is similar to fermion doubling of relativistic fermions, which requires the presence of an even number of Dirac points on a compact manifold). To prove this, we draw a sphere at the origin with radius as shown in Fig. 2(b) and require on the sphere, which enables the definition of a Chern number , where and is the component of as defined above. The Chern number is a topological index quantized to an integer value. For small , Eq. (6) is asymptotically accurate, which gives . Notice that this topological index is invariant as changes adiabatically. The non-zero value of indicates that cannot be decreased to zero smoothly, and hence a point with must exist inside the sphere, which ensures . Since the BZ has a periodic structure and is a compact manifold, the same argument requires zero points outside the sphere with (). This conclusion can be verified explicitly in the simple cubic lattice model studied above, where at each of the eight high symmetry points.

Now we relax the assumption about the simple cubic lattice and allow for a more general structure. Due to time-reversal and space-inversion symmetries, both bands are doubly-generated in a KI. Therefore, the corresponding Bloch wavefunctions and (which are four-component vectors) can be chosen arbitrarily up to a “local” transformation in momentum space. For fermions, , and hence, we require that under time reversal with being the Levi-Civita symbol. For concreteness, we focus below on the case with at and at all other zero points of , which gave us a strong topological insulator in the model with a simple cubic lattice discussed above. For such a band structure, it can be easily checked that the wave-functions can not be defined globally in the entire BZ with the constraint Fu2006 (). However, the BZ can be separated into two regions (c.f., Ref. KSVG ()) () and () with being an arbitrarily small momentum such that only one zero point of , is enclosed inside the sphere, Fig. 2(b). In each of the two regions, singularity-free Bloch wave-functions can be constructed. For example, the valence band has following wave-functions in region

(7) |

and in region

(8) |

Here are the -components of the form-factor, and were defined above Eq. (4), and . These two sets of Bloch wave-functions are connected by a “gauge” transformation at the boundary, , between the regions and , with the matrix .

The topological structure of a 3D time-reversal invariant insulator is determined by the wavefunctions (7,8) on the six high-symmetry planes and FuKane2007 (); Moore2007 (); Roy2009 (). As shown in Fig. 2(b), the boundary between the two regions intersects with and the intersections are circles, . On such a circle, the matrix takes the form of , with being a fixed 3D unit vector, corresponding to a specific “gauge.” As a result, a winding number can be defined on as . However, it is defined modulo only, because a gauge transformation can change and hence change by an even number , where is the azimuth angle of in a high symmetry plane (note that a transformation would violate the symmetry constraint ). For , the corresponding winding numbers are zero for the case we studied here, because they do not intersect with the boundary. The topological indices can be computed from these winding numbers as follows multi ():

(9) |

where is equivalently defined for , or . For the types of KI band structures considered here, the topological indices can be universally determined by choosing a small enough and using the asymptotic form-factor of Eq. (6). As a result, we find that the generic Kondo system is a STI in full agreement with arguments above based on adiabatic deformation of the Hamiltonian onto a simple cubic lattice.

Let us briefly discuss the implications of our results for existing Kondo insulators. From our theory’s we expect that materials in which -electrons are close to integral valence are likely to be weak topological Kondo insulators and thus are unstable with respect to disorder. An interesting example is SmB for which recent LSDA+ band structure calculations Harmon2002 () show the position of the -level equals approximately one sixth of the bandwidth consistent with the core-level spectroscopy measurements of the -level occupation Chazalviel1976 (), placing the quasiparticle f-level of SmB close to the border separating STI and WTI phases. Another promising candidate for the manifestation of topologically nontrivial insulating state is CeNiSn. Recent transport data in CeNiSn CeNiSn2002 () shows suppression of semiconducting behavior in resistivity with increase in sample’s quality, although there is an evidence for the gap formation at K. Given that the -electrons in these systems are predominantly localized, it is tempting to speculate that the CeNiSn is a weak topological Kondo insulator ascribing the semi-metallic transport properties to metallic surface states. These are issues that we hope can be resolved in the near future through more accurate modelling and the use of high precision ARPES and STEM spectroscopy.

To summarize, we have developed a theory of topological 3D Kondo insulators. Within our model topologically nontrivial insulating states are realized over a wide parameter range. In particular, we have shown that strong topological insulating state occurs when the position of the re-normalized level is near the top, or the bottom of the conduction band. This suggests the most likely candidates for this kind of behavior are heavy fermion materials which are more mixed valent or have narrow conduction bands.

This work was supported by DARPA (M. D. and V. G.), JQI-NSF-PFC (K. S.), NSF-CAREER (V. G.), and DOE grant DE-FG02-99ER45790 (P. C.). We would like to thank Joel Moore and David Vanderbilt for discussions related to this project.

## References

- (1) A. Menth, E. Buehler & T. H. Geballe, Phys Rev Lett., 22, 295 (1969).
- (2) For a review, see e.g P. Coleman, Handbook of Magnetism and Advanced Magnetic Materials, Vol 1, 95-148 (Wiley, 2007).
- (3) G. Aeppli & Z. Fisk, Comm. Condens. Matter Phys. 16, 155 (1992); H. Tsunetsugu, M. Sigrist & K. Ueda, Rev. Mod. Phys. 69, 809 (1997); P. Riseborough, Adv. Phys. 49, 257 (2000);
- (4) R. Martin & J. Allen, J. Appl. Phys 50, 7561 (1979).
- (5) L. Fu, C.L. Kane & E.J. Mele, Phys. Rev. Lett. 98, 106803 (2007).
- (6) J.E. Moore & L. Balents, Phys. Rev. B 75, 121306(R)(2007).
- (7) Rahul Roy, Phys. Rev. B 79, 195321 (2009).
- (8) X.L. Qi, T.L. Hughes, & S.C. Zhang, Phys. Rev. B 78, 195424 (2008).
- (9) M. F. Hundley et al., Phys. Rev. B 42, 6842 (1990).
- (10) B. Coqblin & J. R. Schriefer, Phys. Rev. 185, 847 (1969).
- (11) Rebecca Flint, M. Dzero, & P. Coleman, Nature Physics 4, 643 (2008).
- (12) H. Ikeda & K. Miyake, Jour. of Phys. Soc. of Japan 65, 1769 (1996).
- (13) K. A. Kikoin, A. de Visser, K. Bakker, T. Takabatake, Z. Phys. B 94, 79 (1994).
- (14) Juana Moreno & P. Coleman, Phys. Rev. Lett. 84, 342 (2000).
- (15) Liang Fu & C. L. Kane, Phys. Rev. B 76, 045302 (2007).
- (16) A. Kitaev, pre-print arXiv:0901.2686v2 (2009).
- (17) Liang Fu & C. L. Kane, Phys. Rev. B 74, 195312 (2006).
- (18) M. Cheng, K. Sun, V. Galitski, & S. Das Sarma, pre-print arXiv:0908.2805v2 (2009).
- (19) Multiple regions and more complicated boundaries may be required for more complicated band structures and lattices. The intersections between a high-symmetry plane and the boundaries may contain several contours Roy2009 () and the topological index is determined by the sum of all winding numbers.
- (20) V. N. Antonov, B. Harmon, & A. N. Yaresko, Phys. Rev. B 66 165209 (2002).
- (21) J. N. Chazalviel et al., Phys. Rev. B 14, 4586 (1976).
- (22) T. Terashima et al., Phys. Rev. B 66, 075127 (2002).