Möbius Kondo Insulators

Möbius Kondo Insulators

Po-Yao Chang pychang@physics.rutgers.edu Center for Materials Theory, Rutgers University, Piscataway, New Jersey, 08854, USA    Onur Erten Center for Materials Theory, Rutgers University, Piscataway, New Jersey, 08854, USA Max Planck Institute for the Physics of Complex Systems, Dresden, 01187, Germany    Piers Coleman coleman@physics.rutgers.edu Center for Materials Theory, Rutgers University, Piscataway, New Jersey, 08854, USA Department of Physics, Royal Holloway, University of London, Egham, Surrey TW20 0EX, UK

Heavy fermion materials have recently attracted attention for their potential to combine topological protection with strongly correlated electron physics. To date, the ideas of topological protection have been restricted to the heavy fermion or “Kondo” insulators with the simplest point-group symmetries. Here we argue that the presence of nonsymmorphic crystal symmetries in many heavy fermion materials opens up a new family of topologically protected heavy electron systems. Re-examination of archival resistivity measurements in nonsymorphic heavy fermion insulators CeBiPt and CeNiSn reveals the presence of low temperature conductivity plateau, making them candidate members of the new class of material. We illustrate our ideas with a specific model for CeNiSn, showing how glide symmetries generate surface states with a novel Möbius braiding that can be detected by ARPES or non-local conductivity measurements. One of the interesting effects of strong correlation, is the development of partially localization or “Kondo breakdown” on the surfaces, which transforms Möbius surface states into quasi-one dimensional conductors, with the potential for novel electronic phase transitions.

The discovery of topological phases of matter, initiated by the pioneering works on quantum Hall states in 1980sLaughlin (); Thouless (); Haldane () has now evolved into the broad notion of symmetry-protected topological states of matter. Heavy fermion or “Kondo” insulators have recently emerged as a particularly promising platform to study the interplay between topological phases and strong electron interactionsDzero (); Dzero2012 (). In topological Kondo insulators (TKIs), the strong interaction between conduction electrons and local magnetic moments leads to the formation of a narrow gap associated with the development of Kondo singlets which screen the local momentsAeppli (); Fisk (); Tsunetsugu (); Riseborough (). In SmB, the oldest known Kondo insulator, the existence of metallic surface states has been demonstrated by transport experimentsWolgast (); Kim (); Kim2014 (); Thomas () and angle-resolved photoemission spectroscopy (ARPES)Neupane (); Xu2013 (); Jiang2013 (). These results have identified SmB as a promising candidate for a TKI.

One of the important aspects of these materials is the interplay between crystalline symmetry and topological order. To date, the main focus of interest in heavy fermion materials has been limited to the simplest crystalline symmetries. In this work, we expand this notion to a wider class of heavy fermion materials in which unique topological features can arise from the combination of fractional translations and by point group transformations known as nonsymmorphic symmetries. Examples of Kondo insulators with non-symorphic symmetries include CeBiPtHundley (); Cooley_PRB97 (), and CeNiSn, CeRhSb, CeIrSb Hiess94 (); Sera97 (); Kawasaki_PRB2007 (). The observation of resistivity saturation at low temperatures in CeBiPt under pressureCooley_PRB97 () and CeNiSnSlebarski () with Sb doping, closely resemble the conductivity plateau of topological SmB, strongly suggesting that these nonsymmorphic Kondo insulators are topologically nontrivial.

Figure 1: Crystal structure of CeNiSn: showing a the location of the Ce and Ni atoms, denoted by red and blue spheres. There are four equivalent Ce sites related by glide reflections and screw rotations. The Ce and Ni atoms form zigzag chains which reside on two planes labeled (light gray) and (dark gray), stacked along the -axis. b, showing the glide reflection which connects the Ce sites on different chains . c showing the screw rotation which connects the chains in different planes .

To illustrate the topological effects of nonsymmorphic symmetries in heavy fermion systems, here we study CeNiSn as a representative member of this new new class of TKIs. A central element of our theory is a tight-binding model of CeNiSn with all its salient symmetries. The structure of CeNiSn belongs to nonsymmorphic space group No. 62 (), containing three glide reflections, three screw rotations and inversion symmetry, as shown in Fig. 1. Glide reflections and screw rotations are nonsymmorphic symmetries, which combine a point group operation (mirror or rotation) and a fractional lattice translation. We find that these symmetries permit nonsymmorphic Kondo insulators to develop a protected surface state, composed of two Dirac cones. Unlike regular topological insulators, in which scattering between two surface Dirac cones can open a gap without breaking time-reversal symmetry, the autonomy of these surface states is stabilized by glide reflection and time-reversal symmetry Shiozaki2015 (); Wang_arxiv2015 (). Moreover, nonsymmorphic symmetries give rise to a momentum dependent twist that enables the surface states to be detached from the bulk on the glide plane. Following recent studiesLu_2016 (); Fang2015 (); Shiozaki (); Shiozaki2015 (); Wang_arxiv2015 (), we refer to these states as Möbius-twisted surface states. From the bulk-boundary correspondence, we are able to define a topological invariant, and discuss the experimental signatures of such a phase. One of the important effects that sets these topological insulators apart from their weakly interacting counterparts, is the possibility of breakdown of the Kondo effect at the surfaceAlexandrov_PRL2015 (). We find that this breakdown has a particularly dramatic effect on the Möbius-twisted surface states, giving rise to quasi-one dimensional Fermi surfaces.

I Tight-binding Hamiltonian and nonsymmorphic symmetries

We begin by constructing a tight-binding Hamiltonian for CeNiSn. CeNiSn has an orthorhombic -TiNiSi structure belonging to nonsymmorphic space group No. 62, , which contains an inversion , a screw rotation , and a glide reflection , where denotes a rotation about the -axis, refers to the mirror operation in the plane perpendicular to the to the -axis, and is the translation operator along . There are four equivalent Ce sites in the unit cell which we label as as shown in Fig 1. The Cerium sites form zig-zag chains in the ac plane which are stacked along the c direction. Fig. 1b shows how glide reflection connects inter-chain sites . These layers are then arranged in an alternating fashion along the b direction; the alternating layers are related related by the screw rotation as shown in Fig. 1c. In the following discussion, we re-scale the dimensions , and of the unit cell to be unity. When applying the glide reflection and the screw rotation symmetries twice, the system is shifted by a lattice translation, but the process also involves a double reflection or rotation. The half-integer character of the electrons means that reflections or -rotations square to , and this additional factor means appears in the square the square of glide reflection and screw rotations as follows:


Band structure calculationsYanase92 (); Hammond95 (); Hiess97 () indicate that the relevant orbitals near the chemical potential derive from the Ce -electrons and Ni -electrons. We now construct a simplified model involving these two sets of orbitals. A key ingredient of our model is the hybridization between the and states which involves the transfer of one unit of angular momentum from spin, to orbital angular momentum. As a result, the hybridization develops a p-wave form-factorAlexandrov_PRL2015 (), and can be modelled by a simpler model of spin-orbit orbitals hybridizing with -wave conduction electrons. We project the Wannier states of these two sets of orbitals onto the common sites of the Cerium atoms. The resulting tight-binding Hamiltonian has the structure


where is the hybridization matrix, and are the nearest hopping matrices for the conduction and -electrons respectively. The detailed structure of this Hamiltonian, which respects the full nonsymmorphic symmetries of the lattice, is provided in the methods.

Figure 2: A pair of Möbius twisted surface states on the surface. a. Three dimensional plot of the dispersion on the surface. The two Möbius twisted surface states are detached from the bulk along the glide line (). b. Schematic showing dispersion along the glide line , where red and blue lines correspond to states with positive and negative glide eigenvalues , respectively. Along the surface states are doubly degenerate due to a combination of time-reversal and glide symmetry, . The imaginary glide eigenvalues at the point force the members of each Kramers pair to belong to different glide sectors, but at the points the real glide eigenvalues mean that the members of a Kramers pair are in the same glide sector. The connectivity between Kramers pairs at the and points gives rise to the Möbius character, for if one starts at and follows the red loop to , one has to pass a second time around the loop on the blue line, before one returns to the origin.
Figure 3: topological invariants and the corresponding surface states: a, Four phases as the function of the chemical potential of conduction electrons , with and the parameters used in equation (4) are , , and . b-d, The energy spectra are computed in a slab geometry for , respectively. Red and blue lines are the energy dispersion of surface states with positive and negative glide eigenvalues , respectively. A Möbius twisted surface is located at the glide plane () shown in c and d. This Möbius twisted surface is composed by two Dirac cones. Due to Kramers pairs at point have the same glide eigenvalue and , this Möbius twisted surface is detached from the bulk at the glide plane.

Ii Topological surface states with a Möbius twist

One symmetry-preserving surface which respects to the glide reflection is the surface. This surface is perpendicular to the glide plane ( plane) and is also invariant under lattice translations parallel to the surface. The surface energy dispersion as a function of is computed by diagonalizing the Hamiltonian in a slab geometry. The corresponding surface Brillouin zone (BZ) is shown in Fig. 2a. The glide lines on the surface BZ are the set of glide reflection invariant momenta, which are at (path ) and (path ). Along these lines, the Hamiltonian from equation (4) commutes with and can be block diagonalized into two sectors with two eigenvalues for , along the glide lines.

On the glide lines along or , a pair of surface Dirac cones is stabilized by the glide reflection and time-reversal symmetry. To demonstrate this state (see Fig. 2b), we focus on path . At the point, the glide eigenvalues are real (), which implies that the members of each Kramers pair derive from the same glide sector, i.e., the glide eigenvalues for two Kramers pairs are and . By contrast, at the point, the glide eigenvalues are imaginary (), so time reversal inverts the glide eigenvalue, which indicates that the members of each Kramers pair come from opposite glide sectors, i.e., the glide eigenvalues for two Kramers pairs are both . When we connect two Kramers pairs at point to two Kramers pairs at point we obtain the hourglass structure of this surface stateWang_arxiv2015 (), which contains two Dirac cones at the point (Fig. 2b). This surface state contains a Möbius twist, for if we follow the arrow from Fig. 2b along the loop , we need go around the loop twice: once on a red and once on a blue branch, before returning to the origin. Due to this unusual connectivity, the surface state can be detached from the bulk along the loop . Fig. 2a displays the result of a band-calculation on a strip, showing the Möbius-twisted character.

The presence or absence of a Möbius-twisted surface state on the (010) surface defines a variable. When we combine this with the additional variable associated with the possibility of forming a strong topological insulator, by introducing an additional odd number of Dirac cones on every surface, we see that the combination of time-reversal and non-symmorphic symmetries gives rise to a topological invariant , for which we can construct a corresponding index (see Supplementary Information and also Ref. Shiozaki2015 ()), as shown in Fig. 3a. corresponds to a trivial insulator with no gapless surface states. corresponds to a strong topological insulating phase with one single surface Dirac cone (Fig. 3b). corresponds to a nonsymmorphic topological insulator with a Möbius twisted surface state along path (Fig. 3c), while corresponds to a strong topological insulating phase with three surface Dirac cones (Fig. 3d).

In our model calculations, we also observe a double Dirac cone like surface state on surface, where the crossings are located at . However at the mirror plane , this surface state is gapped and is not protected by mirror symmetry and time-reversal symmetry , so this state will likely be absent in the real material.

Figure 4: a, Showing the effect of Kondo breakdown on the dispersion of the Möbius-twisted surface state, causing the Dirac points to sink into the valence band. b, Three dimensional plot showing the quasi-one dimensional dispersion of Möbius-twisted surface state in the presence of Kondo breakdown.

Iii Discussion

We have shown that CeNiSn and Kondo insulators with nonsymmorphic symmetries have the potential to form a new class of topological Kondo insulators with unusual surface states. CeNiSn is of course a low-carrier density metal, with a small Fermi surface derived from an indirect band-gap closure or a lightly doped conduction bandTerashima02 (); Paschen:2016 (), but such small bulk Fermi surfaces are readily localized by disorder or substitution. This is the likely explanation of the observation of a resistivity plateau below 10K in antimony-doped CeNiSnSbSlebarski (), where the observation of a resistivity plateau most likely derives from a metallic topological surface states, as in the case of SmB. Moreover the V-shape density of states deduced from NMR experiments and point-contact spectroscopyEkino (); Nakamura96 (); Nakamura_JPSJ1994 () can be accounted for as a signatures of Dirac cone surface states. The large magnetoresistance for fields perpendicular to the axisTerashima02 () may be a consequence of metallic surface in the plane, combined with an insulating surface on the and planes. Although the above discussion ahs focussed on CeNiSn, we note that the nonsymmorphic Kondo insulator CeBiPt also displays a reistivity plateau. Remarkably, in the presence of pressure, the resistivity plateau persists up to 100KCooley_PRB97 () [See supplementary material]. These early experimental results provide strong circumstantial support for the topological nature of these nonsymmorphic Kondo insulators and provide a strong motivation for further detailed investigation.

The confirmation of these ideas requires a direct probe of Möbius conducting surface states, either by ARPES measurements, or by non-localWolgast () or sample thickness dependent transport measurementsSyers (). One of the interesting challenges is to delineate Möbius surface states from single Dirac surface states. According to Ref. Nakajima (), one-dimensional chiral edge modes from ferromagnetic domain walls on the surface of SmB have been observed with quantized conductance. If an analogous ferromagnetic order emerges on a nonsymmorphic surface (induced by field or magnetic impurities), any Möbius state present will become gapped. The Landau levels of a Möbius surface state will contribute to quantized conductance , where depends on the surface chemical potential. When the Fermi energy is in the gap, the conductance only comes from the lowest Landau level with . The transport measurement presumably comes from both top and bottom surface that give rise to the total conductance , twice that observed in SmB. Another way to detect twisted surface states is from Hall bar measurements. In contrast with the usual topological insulators, where the quantum Hall conductivity switches from to on gating, the twisted surface state will generate three values of quantized Hall conductivity. When the gate voltage is below(above) two cones, we have (). Whereas the gate voltage is in between two cones, the Hall conductance will vanish, .

The strong electron correlations in Kondo insulators opens up the possibility of many interesting phenomena, absent in their weakly interacting counterparts. Here, a particularly important phenomenon is the possibility of a surface breakdown of the Kondo effectAlexandrov_PRL2015 (). Surface Kondo breakdown is based on the observation that the reduced coordination of the rare earth ions at the surface causes a reduction of in the surface Kondo temperature. In principle, competing instabilities such as magnetism can now be activated on the surface. The breakdown of the the Kondo singlets at the surface liberates a large number of conduction electrons which then dope the topological surface states. We have examined the effects of surface Kondo breakdown in Möbius surface states by recomputing the surface spectrum in the absence of the surface f-states. These calculations show that the Kondo breakdown causes the Dirac points in Möbius surface states to sink into the valence band as shown in Fig. 4a which then generates large Fermi surfaces [see Fig. 4b]. The detached nature of the Möbius surface state causes the resulting Fermi surfaces to become quasi one-dimensional along direction. The interaction of the partially unscreened surface local moments with these quasi one-dimensional Fermi surfaces is expected to lead to a wide variety of surface electronic instabilities, including unconventional superconductivityBennemann_Book () and charge or spin density wave instabilities.

Another interesting future direction is the possibility of nonsymmorphic topological superconductors. Promising candidates are UCoGe and URhGe, which share the same space group as CeNiSnHuy_PRL2007 (); Aoki_Nature2001 (). These materials exhibit spin-triplet superconductivity in coexistence with ferromagnetism. The topological classification of such superconductors is an intriguing future prospect.

Iv Methods

In the momentum space, the tight-binding Hamiltonian is , where is a sixteen component spinor, with



where is the hybridization matrix, and are the nearest hopping matrices for conduction and -electrons, respectively. In order to simplify our calculation, we introduce four sets of Pauli matrices: acts on the spin basis; acts on the basis of conduction electrons and -electrons; acts on the basis of the atom labels and ; acts on the basis of the layer labels and .

From equation (4), the hybridization matrix has the form


with , , , and .

The nearest hopping matrices for conduction electrons and -electrons are

where , are the hopping amplitudes along -direction, and are the bare energies of the isolated conduction electrons and -electrons. In the Supplementary Information we perform the construction of this tight-binding Hamiltonian in detail.

We write down the matrix representations of symmetries as follows:

  1. Time-reversal symmetry, , where with being the complex conjugation operator.

  2. Inversion symmetry, , where .

  3. Glide reflection symmetry , , where

  4. Screw rotation symmetry , , where

  5. Mirror symmetry , , where

In the spin-orbit coupled systems, reflection and rotation square to . We have , , and .

V Acknowledgments

The authors would like to thank Silke Paschen and Tyrel McQueen for discussions about CeNiSn. This work was supported by the Rutgers Center for Materials Theory group postdoc grant (Po-Yao Chang), US National Science Foundation grant grant DMR-1309929 (Onur Erten) and US Department of Energy grant DE-FG02-99ER45790 (Piers Coleman).

Vi Author contributions

All authors performed the calculations, discussed the results and prepared the manuscript.

Vii Competing financial interests

The authors declare no competing financial interests.


Appendix A Möbius Kondo Insulators: Supplementary Material

a.1 A tight-binding model based on nonsymmorphic symmetries

Here we derive the tight-binding Hamiltonian for CeNiSn using the following procedure:

  1. We define the symmetry operations on the fermionic operators.

  2. We then construct the Hamiltonian by writing down the all the nearest neighbor hopping and hybridization terms which respect the nonsymmorphic symmetries.

In order to simplify our calculation, we introduce four sets of Pauli matrices: acts on the spin basis; acts on the basis of conduction electrons and -electrons; acts on the basis of the atom labels and ; acts on the basis of the layer labels and .

In CeNiSn, -electrons of Ce atoms hybridize with -electrons of Ni atoms. Since the total angular momentum difference between these two states is one, we orbitally “downfold” the tight-binding model, replacing it by an equivalent model, reducing the total angular momentum of each band by two units. The resulting model involves the hybridization of spin-orbit coupled p-electrons with s-electrons The highly localized states are then modelled as -wave Kramers doublets,


while the mobile conduction electrons are -wave,


For further simplicity, we project the conduction electron Wannier states onto the sites of -electrons.

a.2 Symmetries:

Here, we describe the actions of the various symmetry operations in space group No. 62 () on the Fermi operators. This group contains three screw rotations , three glide reflections and an inversion . In fact, this set of six operators can all be generated from , in combination with inversion and translation operators, so it is sufficient for us to focus on these two non-symmorphic operators.

a.2.1 Glide reflection

The glide reflection transforms the -electrons and conduction electrons as follows


where is the co-ordinate of the unit cell, and and are the site and layer indices of atoms within one unit cell, and are unit vectors in the and directions, respectively, while is the spin index. The glide reflection acts in a similar way on the conduction electrons


a.2.2 Screw rotation

The screw rotation transforms the -electrons and conduction electrons as follows:




where is the unit vector in direction.

a.2.3 Inversion

The inversion transforms the -electrons and conduction electrons in the following way,




a.2.4 Other symmetries:

By combining inversion, translations, and the above two nonsymmorphic symmetries we can obtain all the remaining symmetries. Thus , and , since


and likewise,


In a similar fashion, we obtain and , since




a.3 Construction of the Hamitonian

a.3.1 Hopping terms

Figure 5: Hopping Hamiltonian: a, The hopping elements and with on the plane. b, The hopping elements with on the plane.

Now we consider the nearest neighbor hopping terms in the Hamiltonian [see Fig. 5].


where is the hopping amplitude for electrons along direction, and is the on-site chemical potential for electrons. The hopping Hamiltonian is explicitly invariant under , and , and since all the other non-symmorphic symmetries can be expanded in terms of these three (plus translation), the hopping is invariant under the full non-symmorphic group.

In momentum space, the hopping Hamiltonian is , where where is a sixteen component creation operator, with





a.3.2 Hybridization terms

Figure 6: Hybridization Hamiltonian on layer: a, The hybridization matrix between and is . b, The hybridization elements , and on the plane.

We now consider the hybridization terms between -electrons and conduction electrons. Let us start with layer . The hybridization Hamiltonian on layer is (see Fig. 6b)


where in the second line is the glide-reflection of the first four terms. This guarantees that . The hybridization matrices have the structure , and (see Fig. 6a). In momentum space, the hybridization Hamiltonian on layer is , where




The hybridization Hamiltonian on the layer can be obtained by performing the screw rotation on , .

Figure 7: Interlayer hybridization Hamiltonian: a, The interlayer hybridization matrix between and is . b, The interlayer hybridization elements and on the plane.

Next, we write down the interlayer hybridization terms. The interlayer hybridization Hamiltonian (see Fig. 7b) is


where denote a repeat of the first two terms, and (see Fig. 7a). By construction, interlayer hybridization is invariant under and , and is thus invariant under the full non-symmorphic group In momentum space, the interlayer hybridization Hamiltonian is , where