Cubic Topological Kondo Insulators
Abstract
Current theories of Kondo insulators employ the interaction of conduction electrons with localized Kramers doublets originating from a tetragonal crystalline environment, yet all Kondo insulators are cubic. Here we develop a theory of cubic topological Kondo insulators involving the interaction of spin quartets with a conduction sea. The spin quartets greatly increase the potential for strong topological insulators, entirely eliminating the weaktopological phases from the diagram. We show that the relevant topological behavior in cubic Kondo insulators can only reside at the lower symmetry X or M points in the Brillouin zone, leading to three Dirac cones with heavy quasiparticles.
pacs:
72.15.Qm, 73.23.b, 73.63.Kv, 75.20.HrOur classical understanding of order in matter is built around Landau’s concept of an order parameter. The past few years have seen a profound growth of interest in topological phases of matter, epitomized by the quantum Hall effect and topological band insulators, in which the underlying order derives from the nontrivial connectedness of the quantum wavefunction, often driven by the presence of strong spinorbit coupling (1); (2); (3); (4); (5); (6); (7); (8); (9).
One of the interesting new entries to the world of topological insulators, is the class of heavy fermion, or “Kondo insulators” (10); (11); (12); (13); (14); (15); (16). The strongspin orbit coupling and highly renormalized narrow bands in these intermetallic materials inspired the prediction (12) that a subset of the family of Kondo insulators will be Z topological insulators. In particular, the oldest known Kondo insulator SmB(17) with marked mixed valence character, was identified as a particularly promising candidate for a strong topological insulator (STI): a conclusion that has since also been supported by bandtheory calculations (13); (16). Recent experiments (18); (19); (20) on SmB have confirmed the presence of robust conducting surfaces, large bulk resistivity and a chemical potential that clearly lies in the gap providing strong support for the initial prediction.
However, despite these developments, there are still many aspects of the physics in these materials that are poorly understood. One of the simplifying assumptions of the original theory (12) was to treat the states as Kramer’s doublets in a tetragonal environment. In fact, the tetragonal theory predicts that strong topological insulating behavior requires large deviations from integral valence, while in practice Kondo insulators are much closer to integral valence (11). Moreover, all known Kondo insulators have cubic symmetry, and this higher symmetry appears to play a vital role, for all apparent “Kondo insulators” of lower symmetry, such as CeNiSn(21) or CeRuSn(22) have proven, on improving sample quality, to be semimetals. One of the important effects of high symmetry is the stablization of magnetic fquartets. Moreover, Raman(23) experiments and various bandtheory studies(24); (25) that it is the Kondo screening of the magnetic quartets that gives rise to the emergence of the insulating state.
Motivated by this observation, here we formulate a theory of cubic topological Kondo insulators, based on a lattice of magnetic quartets. We show that the presence of a spinquartet greatly increases the possibility of strong topological insulators while eliminating the weaktopological insulators from the phase diagram, Fig. 1. We predict that the relevant topological behavior in simple cubic Kondo insulators can only reside at the lower point group symmetry X and M points in the Brillouin zone (BZ), leading to a three heavy Dirac cones at the surface. One of the additional consequences of the underlying Kondo physics, is that the coherence length of the surface states is expected to be very small, of order a lattice spacing.
While we outline our model of cubic Kondo insulators with a particular focus on SmB, the methodology generalizes to other cubic Kondo insulators. SmB has a simple cubic structure, with the B clusters located at the center of the unit cell, acting as spacers which mediate electron hopping between Sm sites. Bandtheory (25) and XPS studies (27) show that the 4f orbitals hybridize with bands which form electron pockets around the X points. In a cubic environment, the orbitals split into a doublet and a quartet, while the fivefold degenerate orbitals are split into double degenerate and triply degenerate orbitals. Band theory and Raman spectroscopy studies (23) indicate that the physics of the orbitals is governed by valence fluctuations involving electrons of the quartet and the conduction states, . The () quartet consists of the following combination of orbitals: . This then leads to a simple physical picture in which the quartet of states hybridizes with an quartet (Kramers plus orbital degeneracy) of states to form a Kondo insulator.
To gain insight into how the cubic topological Kondo insulator emerges it is instructive to consider a simplified onedimensional model consisting of a quartet of conduction bands hybridized with a quartet of bands (Fig. 2a). In one dimension there are two high symmetry points: () and X (), where the hybridization vanishes (12); (14); (26)). Away from the zone center , the and quartets split into Kramers doublets. The topological invariant is then determined by the product of the parities of the occupied states at the and X points. However, the quartet at the point is equivalent to two Kramers doublets, which means that is always positive, so that and a onedimensional topological insulator only develops when the and bands invert at the X point.
Generalizing this argument to three dimensions we see that there are now four high symmetry points , X, M and R. The bands are fourfold degenerate at both and R points which guarantees that (Fig. 2b). Therefore, we see that the 3D topological invariant is determined by band inversions at X or M points only, . If there is a band inversion at the X point, we get . In this way the cubic character of the Kondo insulator and, specifically, the fourfold degeneracy of the orbital multiplet protects the formation of a strong topological insulator.
We now formulate our model for cubic topological Kondo insulators. At each site, the quartet of f and d holes is described by an orbital and spin index, denoted by the combination (, ). The fields are then given by the eight component spinor
(1) 
where destroys an dhole at site j, while is the the Hubbard operator that destroys an hole at site j. The tightbinding Hamiltonian describing the hybridized  system is then
(2) 
in which the nearest hopping matrix has the structure
(3) 
where the diagonal elements describe hopping within the  and  quartets while the offdiagonal parts describe the hybridization between them, while is the vector linking nearest neighbors. The various matrix elements simplify for hopping along the zaxis, where they become orbitally and spin diagonal:
(4) 
where and is the ratio of orbital hopping elements. In the above, the overlap between the orbitals, which extend perpendicular to the zaxis is neglected, since the hybridization is dominated by the overlap of the the orbitals, which extend out along the zaxis. The hopping matrix elements in the and directions are then obtained by rotations in orbital/spin space. so that and where and denote 90 rotations about the y and negative x axes, respectively.
The Fourier transformed hopping matrices can then be written in the compact form
(5) 
where . Here are the bare energies of the isolated d and fquartets, while , and (). The hybridization is given by
(6) 
where we denote . Note how the hybridization between the even parity dstates and oddparity fstates is an odd parity function of momentum .
To analyze the properties of the Kondo insulator, we use a slave boson formulation of the Hubbard operators, writing , where creates an fhole in the quartet while denotes the singlet filled shell, subject to the constraint at each site.
We now analyze the properties of the cubic Kondo insulator, using a meanfield treatment of the slave boson field , replacing the slaveboson operator at each site by its expectation value: so that the f hopping and hybridization amplitude are renormalized: and . The meanfield theory is carried out, enforcing the constraint on the average. In addition, the chemical potentials and for both electrons and holes are adjusted selfconsistently to produce a band insulator, . This condition guarantees that four out of eight doubly degenerate bands will be fully occupied. The details of our meanfield calculation are given in the Supplementary Materials section. Here we provide the final results of our calculations.
In Fig. 3 we show that the magnitude reduces with temperature, corresponding to a gradual rise in the Sm valence, due to the weaker renormalization of the electron level. The degree of mixed valence of Sm is given then by . In our simplified meanfield calculation, the smooth temperature crossover from Kondo insulating behavior to local moment metal at high temperatures is crudely approximated by an abrupt secondorder phase transition.
Fig. 4 shows the computed band structure for the cubic Kondo insulator obtained from meanfield theory, showing the band inversion between the  and bands at the X points that generates the strong topological insulator. Moreover, as the value of the bare hybridization increases, there is a maximum value beyond which the bands no longer invert and the Kondo insulator becomes a conventional band insulator.
One of the interesting questions raised by this work concerns the many body character of the Dirac electrons on the surface. Like the lowlying excitations in the valence and conduction band, the surface states of a TKI involve heavy quasiparticles of predominantly character. The characteristic Fermi velocity of these excitations is renormalized with respect to the conduction electron band group velocities, where is the mass renormalization of the felectrons. In a band topological insulator, the penetration depth of the surface states , where is the bandgap, scale that is significantly larger than a unitcell size. Paradoxically, even though the Fermi velocity of the Dirac cones in a TKI is very low, we expect the characteristic penetration depth of the heavy wavefunctions into the bulk to be of order the lattice spacing . To see this, we note that , where the indirect gap of the Kondo insulator is of order the Kondo temperature . But since , where is the width of the conduction electron band, this implies that the penetration depth of the surface excitations is given by the size of the unit cell. Physically, we can interpret the surface Dirac cones as a result of broken Kondo singlets, whose spatial extent is of order a lattice spacing. This feature is likely to make the surface states rather robust against the purity of the bulk.
Various interesting questions are raised by our study. Conventional Kondo insulators are most naturally understood as a strongcoupling limit of the Kondo lattice, where local singlets form between a commensurate number of conduction electrons and localized moments. What then is the appropriate strong coupling description of topological Kondo insulators, and can we understand the surface states in terms of broken Kondo singlets? A second question concerns the temperature dependence of the hybridization gap. Experimentally, the hybridization gap observed in Raman studies(23) is seen to develop in a fashion strongly reminiscent of the meanfield theory. Could this indicate that fluctuations about meanfield theory are weaker in a fully gapped Kondo lattice than in its metallic counterpart?
We end with a few comments on the experimental consequences of the above picture. One of the most dramatic consequences of the quartet model is the prediction of three Dirac cones of surface excitation of substantially enhanced effective mass. By contrast, were the underlying groundstate of the fstate a Kramers doublet, then we would expect a single Dirac cone excitation. These surface modes should be observable via various low energy spectroscopies. For instance, the heavy mass of the quasi particles should appear in Shubnikovde Haas or cyclotron resonance measurements. The quasiparticle mass will depend on the Fermi momentum of the Dirac cones. Scanning tunneling spectroscopy measurements of the quasiparticle interference created by the Dirac conese, and high resolution ARPES measurements may provide a direct way to observe the predicted three Dirac cones.
To summarize, we have studied the cubic topological Kondo insulator, incorporating the effect of a fourfold degenerate multiplet. There are two main effects of the quartet states: first, they allow the low fractional filling of the band required for strong topological insulating behavior to occur in the almost integral valent environment of the Kondo insulator; second, they double the degeneracy of the bandstates at the highsymmetry and points in the Brillouin zone, removing these points from the calculation of the Z topological invariant so that the only important band crossings occur at the the X or M points. This leads to a prediction that three heavy Dirac cones will form on the surfaces (13); (16).
We would like to thank A. Ramires, V. Galitski, K. Sun, S. Artyukhin and J. P. Paglione for stimulating discussions related to this work. This work was supported by the Ohio Board of Regents Research Incentive Program grant OBRRIP220573 (M.D.), DOE grant DEFG0299ER45790 (V. A & P. C.), the U.S. National Science Foundation I2CAM International Materials Institute Award, Grant DMR0844115.
Supplementary Materials for Cubic Topological Kondo Insulators
Victor Alexandrov, Maxim Dzero and Piers Coleman
March 6, 2018
Appendix A Supplementary materials for Cubic Topological Kondo insulators
These notes provide:

details of the derivation of the tightbinding Hamiltonian for a cubic Kondo insulator.

derivation of the meanfield theory for the infinite limit

derivation of the meanfield equations.
a.1 Rotation matrices
To construct the Hamiltonian, we evaluate the hopping matrices along the zaxis, and then carry out a unitary transformation to evaluate the corresponding quantities for hopping along the x and y axes.
Consider a general rotation operator
(7) 
where and describe rotations about a principle axis of the crystal in real and spin space respectively. The Hamiltonian in a cubic environment is invariant under these transformations: . We now write the directional dependence of the Hamiltonian explicitly.
(8) 
We can always choose the rotation to transform in the cyclic manner: , then substituting ,
(9) 
hence we can assume the 3D Hamiltonian of the following form
(10) 
Where to rotate in the opposite direction we use . Using the Wigner Dfunctions we can construct the rotation in angular momentum space
(11) 
Here, denotes the rotation operator in the ”ZYZ” convention, corresponding to a rotation around the axis, followed by a rotation around the axis, and then the new  axis. The matrix elements of this operator are then
so that upon rotating the state using we obtain , and . One can now obtain the transformation matrices for a given multiplet ( etc). For the doublet we are to read off the matrices from the transformation equation
(12) 
where is an doublet. Similarly for the quartet,
(13) 
where the quartet is denoted by , where . The result is,
(14) 
Note that the rotation cyclically exchanges , thus applying it three times gives overall ”” due to fermionic statistics.
Finally, we can use (14) to derive the hamiltonian defined in (5) and (6) of the main paper.
a.2 Details of the meanfield theory
The full Hamiltonian for the problem also contains a term describing the local Hubbard interaction between the electrons:
(15) 
We consider the infinite , where we can project out all states with occupation number larger than one by replacing the bare felectron fields by Hubbard operators. We represent the Hubbard operators using a slave boson representation, as follows:
(16) 
supplemented by a constraint of no more than one electron per each site ():
(17) 
The partition function corresponding to the model Hamiltonian above and with constraint condition (17) reads :
(18) 
where the Lagrangian is
(19) 
Meanfield (saddlepoint) approximation corresponds to the following values of the bosonic fields:
(20) 
where both and are independent. In the meanfield theory we choose a “radial” gauge where the the phase of the field has been absorbed into the felectron fields. Also, for an insulator, we need a filled quartet of states at each site, so that
(21) 
Note, the parameter also renormalizes the hopping elements. Indeed, it follows that when , . However, the onsite occupancy is unrenormalized by the slave boson fields, since in the infinite limit, the onsite occupancy .
The first three terms together with the last term in the Lagrangian (19) can be written using the new fermionic basis and (20). It follows:
(22) 
where the renormalized electron dispersion is
(23) 
Correspondingly, the hybridization matrix in the new fermionic basis can be obtained via unitary transformation with the matrix
(24) 
Thus we need to calculate the elements of the matrix
(25) 
We can use the following relations
(26) 
In our subsequent discussion it will be convenient to write in a more compact form. To derive the corresponding expression we first recall that can be written as follows
(27) 
Again, note the change from minus sign to plus sign in front of term. This yields the agreement with the hybridization Hamiltonian obtained from the rotations method. We also have
(28) 
After some algebra we find
(29) 
These results become much more transparent if we express and in terms of the angle :
(30) 
Then we see that
(31) 
where now
(32) 
Similarly, we find
(33) 
Thus, we can rewrite (25) as follows
(34) 
Next we derive the meanfield equations.
a.3 Derivation of the meanfield equations
To derive the mean field equations we first integrate out electrons by making the following change of variables in the path integral:
(35) 
Then the resulting action is Gaussian and the electrons can be integrated out. This yields the effective action of the form
(36) 
The renormalized electron correlation function has a block diagonal form:
(37) 
where the diagonal elements are given by
(38) 
and we have introduced the following functions
(39) 
Then the determinant of the matrix is