arxiv: v2 [cond-mat.mes-hall] 15 Sep 2011

Transcription

1 Topologicl phses nd surfce flt bnds in superconductors without inversion symmetry rxiv:.438v [cond-mt.mes-hll] 5 Sep Andres P. Schnyder nd Shinsei Ryu Mx-Plnck-Institut für Festkörperforschung, Heisenbergstrsse, D-7569 Stuttgrt, Germny Deprtment of Physics, University of Cliforni, Berkeley, CA 947, USA (Dted: October, 8 We exmine different topologicl phses in three-dimensionl non-centrosymmetric superconductors with time-reversl symmetry by using three different types of topologicl invrints. Due to the bulk boundry correspondence, non-zero vlue of ny of these topologicl numbers indictes the ppernce of zero-energy Andreev surfce sttes. We find tht some of these boundry modes in nodl superconducting phses re dispersionless, i.e., they form topologiclly protected flt bnd. The region where the zero-energy flt bnd ppers in the surfce Brillouin zone is determined by the projection of the nodl lines in the bulk onto the surfce. These dispersionless Andreev surfce bound sttes hve mny observble consequences, in prticulr, zerobis conductnce pek in tunneling mesurements. We lso find tht in the gpless phse there pper Mjorn surfce modes t time-reversl invrint moment which re protected by Z topologicl invrint. PACS numbers: f, 73..-r, 3.65.Vf, v, c, 73..Fz The hllmrk of topologicl insultors nd superconductors (SCs is the existence of topologiclly protected conducting boundry modes. The recent experimentl observtion of these edge nd surfce sttes in spin-orbit inducedz topologicl insultors in two nd three dimensions [], respectively, hs led to surge of interest nd excitement []. An exhustive clssifiction of topologiclly protected boundry modes occurring in gpped free fermion systems in terms of symmetry nd sptil dimension ws given in Refs. [3 5]. Interestingly, this clssifiction scheme, which is known s the periodic tble of topologicl insultors nd SCs, predicts three dimensionl (3D topologicl SC which stisfies time-reversl symmetry, but breks spin-rottion symmetry. Indeed, the B phse of 3 He is one exmple of this so-clled clss DIII topologicl superfluid, whose different topologicl sectors cn be distinguished by n integer topologicl invrint. Recent trnsverse coustic impednce mesurements in 3 He-B confirmed the existence of the predicted surfce Mjorn bound stte [6]. However, finding n electronic nlog of the superfluid B phse of 3 He remins n outstnding chllenge. In this pper we rgue tht some of the 3D non-centrosymmetric SCs might be exmples of electronic topologicl SCs in symmetry clss DIII [7]. We nlyze the topologicl phse digrm of these systems nd demonstrte quite generlly tht djcent to fully gpped topologicl phses there exist non-trivil gpless superconducting phses with topologiclly protected nodl lines (rings. To chrcterize these gpless lines we introduce set of topologicl invrints nd show tht, due to the bulkboundry correspondence, the presence of topologiclly stble nodl rings implies the ppernce of dispersionless zeroenergy Andreev surfce sttes. These topologiclly protected surfce flt bnds mnifest themselves in scnning tunneling spectroscopy (STS s zero bis conductnce pek, feture which could be used s n experimentl signture of the topologicl non-trivility. In non-centrosymmetric SCs the bsence of inversion in the crystl structure genertes ntisymmetric spin-orbit couplings (SOC nd leds to mixing of spin-singlet nd spintriplet piring sttes. These re the properties tht give rise to topologiclly non-trivil qusi-prticle bnd structures in these systems. Strting with CePt 3 Si [8], multitude of noncentrosymmetric SCs hs recently been discovered, including, mong others, Li Pd x Pt 3 x B [9]. Model Hmiltonin As generic phenomenologicl description pplicble to ny of the forementioned mterils we employ single bnd model with ntisymmetric SOC nd tret superconductivity t the men field level. Thus, let us consider H = k Ψ k H(kΨ k with Ψ k = (c k,c k,c k,c k T, where c σk is the electron cretion opertor with spin σ nd momentum k nd the Bogoliubov-de Gennes (BdG Hmiltonin is given by ( h(k (k H(k = (k h T. ( ( k The norml stte Hmiltoninh(k describes non-intercting electrons in crystl without inversion centerh(k = ε k σ + γ k σ, where ε k = ε k is the spin-independent prt of the spectrum,σ,,3 stnd for the three Puli mtrices, ndσ denotes the unit mtrix. The second term in h(k represents n ntisymmetric SOC with coupling constntγ k. Due to the presence of the prity breking SOC γ k the order prmeter in q. ( is in generl n dmixture of spin-singlet ψ k nd spin-triplet d k piring sttes (k = (ψ k σ +d k σ(iσ, where ψ k nd d k re even nd odd functions of k, respectively. The direction of the spin-triplet component d k is ssumed to be prllel to γ k, s for this choice the ntisymmetric SOC is not destructive for triplet piring []. Hence, we prmetrize the d-vector nd the SOC s d k = t l k nd γ k = αl k, respectively. For the spinsinglet component we ssume s-wve piring ψ k = s nd choose the mplitudes t,s to be rel nd positive. In order to exemplify the topologicl properties of the BdG Hmiltonin (, we consider norml stte tight-binding bnd structure on the cubic lttice ε k = t (cosk x +cosk y +cosk z µ, with hopping mplitude t nd chemicl potentil µ. We will set (t,µ,α, t = (4.,4.8,.,. henceforth. The specific form of the SOC γ k depends on the crystl structure [], i.e., gγ g k = γ k,

2 g ( 3 ν = + ν = - ν = O g s (b - ν = s C 4v FIG.. (color online. Phse digrm s function of singlet piring mplitude s nd SOC g [see qs. ( nd (3] for the point group ( O nd (b C 4v. The gpped phses re chrcterized by the winding number ν with ν = (white, ν = ± (drk/light blue htched, ν = 5 (light brown htched, nd ν = +7 (drk brown htched. Grey shded nd dotted regions re nodl superconducting phses with N C = ± (red/blck dotted, N C = ± (light/drk grey, N C3 = + (white dotted, nd N C4 = + (grey. where g is ny symmetry opertion in the point group G of the crystl. Hving in mind Li Pd x Pt 3 x B, we ssume for the pseudovector l k the following form comptible with the symmetry requirements of the cubic point group O l k = sink x sink y g sink x(cosk y +cosk z sink y (cosk x +cosk z, ( sink z sink z (cosk x +cosk y with the constntg, nd where we neglect higher order terms. Furthermore, we lso consider the point group C 4v, relevnt for CePt 3 Si, in which cse l k reds l k = (sink y ê x sink x ê y (3 +g sink x sink y sink z (cosk x cosk y ê z. It is importnt to note tht the qusi-prticle bnd topology of H(k, s defined by q. (, is minly determined by the momentum dependence of l k long the Fermi surfce sheets. Hence, the results we obtin re expected to remin qulittively unchnged upon inclusion of further-neighbor hopping terms in the bnd structureε k. Topologicl Invrints To chrcterize the topologicl properties of H(k we introduce three different topologicl invrints. But before doing so, we observe tht H(k stisfies both time-reversl symmetry (TRS, with T =, nd prticle-hole symmetry (PHS, with C = +, which re the defining symmetry properties of symmetry clss DIII in the terminology of Ref. [3]. Combining TRS nd PHS yields third discrete symmetry, the chirl symmetrys = T C, i.e., there is unitry mtrixs which nticommutes withh(k. It is importnt to note tht while both TRS nd PHS relteh(k to H T ( k, S is symmetry which is stisfied by H(k t ny given pointk in the Brillouin zone (. As shown in Ref. [3] topologicl sectors in the fully gpped phses ofh(k re distinguished by the winding number ν = d 3 k 4π εµνρ Tr [ (q µ q(q ν q(q ρ q ],(4.5 where the integrl is over the st nd q(k is the offdigonl block of the flt-bnd mtrix ofh(k []. In the nodl superconducting phses the winding number ν is no longer quntized. However, we cn consider H(k restricted to D loops in reciprocl spce nd define topologicl number in terms of D momentum spce loop integrl to chrcterize the topology of the gpless phses. We observe tht H(k confined to generic momentum spce loop no longer stisfies TRS nor PHS, but it still obeys chirl symmetry S. Hence, H(k restricted to loop in the belongs to symmetry clss AIII [3] nd its topologicl chrcteristics re described by the D winding number N L = πi L dltr[q (k l q(k], (5 where the integrl is evluted long the loop L in the. Observe tht for ny closed loopltht does not intersect with gpless regions in the, N L is quntized to integer vlues. If L is chosen such tht it encircles line node, then N L determines the topologicl stbility (i.e., the topologicl chrge of the gpless line [3, 4]. Finlly, we lso consider H(k restricted to time-reversl invrint (TRI loop L, which is mpped onto itself under k k. In tht cse we obtin D Hmiltonin stisfying both TRS nd PHS (i.e., belonging to symmetry clss DIII. The topologicl properties of such D system re chrcterized by the following invrint [] W L = K Pf [ q T (K ] / det[q(k], (6 wherek denotes the two TRI moment on the looplndpf is the Pfffin. Note tht W L is either + or - for ny TRI loop tht does not cross gpless regions in the. Topologicl Phse digrm Numericl evlution of the topologicl numbers (4 nd (5 yields the topologicl phse digrm of H(k, which is shown in Fig. s function of second order SOC g nd reltive strength of singlet nd triplet piring components. Fully gpped phses with different topologicl properties (i.e., the phses lbeled by ν = ±,, 5,+7 re seprted in the phse digrm by regions of nodl superconducting phses (grey shded nd dotted res. The fully gpped phses with ν = ± re electronic nlogs of 3 He-B. The nodl superconducting phses exhibit topologiclly stble nodl rings, which re centered round high symmetry xes of the (see Figs. nd 3. In order to determine the topologicl chrcter of these nodl lines (nd hence of the corresponding gpless phses it is sufficient to consider the topologicl invrint N L only for loops L tht run long high symmetry xes. Thus, for the cubic point groupo we choose the loops C : Γ M X Γ ndc : Γ M R Γ, wheres for the tetrgonl point group C 4v we consider C 3 : Γ Z R X Γ nd C 4 : Γ Z A M Γ. For the cubic point group we find tht whenever (N C,N C = (±, there re topologiclly stble nodl rings centered round the ( xis (nd symmetry relted directions. When (N C,N C = (,± the gpless lines re oriented long the ( xis, wheres when(n C,N C = (±,± the rings re locted round the

3 3 R M ( (c X k y / π Γ k x / π (b (d A M ( (c R X k y / π Z Γ k x / π (b (d FIG.. (color online Nodl rings ( nd ( surfce sttes (c,d for the point group O with (g, s = (.3,.5. This prmeter choice corresponds to the red dotted region in Fig.. (b Topologicl invrint N (, q. (7, s function of surfce momentum k. Grey nd drk blue indicte N ( = ±, while light blue is N ( =. (c Bnd structure for slb with ( fce s function of surfce momentum k with k =.75π. (d nergy dispersion of the lowest lying stte with positive energy. The color scle is such tht blck corresponds to zero energy. The sttes t zero energy in (c nd (d re loclized t the surfce. The flt bnds in (c nd (d re singly degenerte (i.e., one brnch per surfce, wheres the linerly dispersing zero mode t the center of the in (d is doubly degenerte. ( direction. (A similr nlysis lso holds for the group C 4v. Andreev surfce sttes A non-zero quntized vlue of ny of the three topologicl numbers (4, (5 nd (6 implies the existence of zero-energy Andreev surfce sttes. First of ll, in fully gpped phses with topologiclly non-trivil chrcter there pper linerly dispersing Mjorn surfce modes [3, 5 7]. In order to understnd the ppernce of zeroenergy Andreev surfce sttes in the gpless phses, we now mke use of the topologicl invrintn L with cleverly chosen loop L. Let us consider q. ( in slb configurtion with(lmn fce. In this geometry the HmiltoninH (lmn retins trnsltionl invrince long the two independent directions prllel to the(lmn surfce. Hence,H (lmn (k cn be viewed s fmily of D systems prmetrized by the two surfce momentk = (k,k. Since H (lmn(k obeys chirl symmetry (but breks in generl TRS nd PHS, its topologicl properties re given by the D winding number of clss AIII N (lmn (k = dk Tr [ q (k q(k ], (7 πi wherek is the bulk momentum perpendiculr to the surfce, nd = / k. Note tht N (lmn is the sme s N L, FIG. 3. (color online Sme s Fig. but for the point group C 4v, for slb with ( fce, nd with (g, s = (.,.5. This prmeter choice corresponds to the white dotted re in Fig. b. q. (5, withlchosen longk, following non-contrctible cycle of the torust 3. Now, the key observtion is tht the bove line integrl is closely relted the loop integrl N L, with L = C i, tht determines the topologicl chrge of the superconducting nodl lines. Tht is, for those surfce moment k for which the loop long k in q. (5 psses through just one non-trivil nodl ring, N (lmn (k is equl to the topologicl chrge of this given nodl ring. Hence, if we plotn (lmn (k s function of surfce moment (see Figs. b,3b, we find tht the boundries seprting regions with different winding number re identicl to the projection of the nodl lines onto the (lmn plne. Furthermore, since non-zero quntized vlue of N (lmn implies the existence of zero energy sttes t the end points of the D HmiltoninH (lmn (k [3, 8], we find tht there re zero-energy Andreev bound sttes on the(lmn surfce locted within the projected nodl rings. This conclusion is corroborted by numericl computtions of the zeroenergy surfce sttes both for the point groupo ndc 4v (see Figs. nd 3. When two nodl rings overlp in the (lmn projection of the, then the quntized vlue of N (lmn in the overlpping region is determined by the dditive contribution of the topologicl chrges of the two rings. In prticulr, one cn hve sitution where the two contributions cncel, in which cse there is no zero-energy surfce stte within the overlpping region. Finlly, using n nlogous rgument s in the previous prgrph, we cn lso employ the Z number (6 to deduce the presence of zero energy modes t TRI moment of the surfce []. One exmple of this is the Krmers pir of surfce zero modes locted t the center of the surfce in Fig. d (cf. Refs. [6, 7]. Remrkbly, this is surfce M-

4 4 ρ(w (rb. units 6 4 ( O ( SDOS bulk DOS - - w (b C 4v ( SDOS bulk DOS - - w FIG. 4. (color online. Surfce nd bulk density of sttes for the point group ( O nd (b C 4v. The surfces re oriented perpendiculr to the ( nd ( xes, respectively. The employed prmeter vlues re the sme s in Figs. nd 3. jorn mode in gpless (nodl superconducting phse [9]. xperimentl signtures One of the most direct signtures of the topologicl spects of non-centrosymmetric SCs re the surfce Andreev bound sttes. These cn be probed by ngleresolved photoemission mesurements, or by STS of the surfce density of sttes (SDOS. STS hs proved to be n effective tool to explore surfce sttes of two-dimensionl unconventionl superconductors, see, e.g., [ 4]. The bulk density of sttes of 3D gpless SCs with nodl lines vnishes linerly t zero energy. In contrst, the surfce flt bnds led to diverging zero-energy pek in the SDOS (see Fig. 4. The zero-bis pek in the SDOS is strongly dependent on the surfce orienttion. From this dependence it is in principle possible to (prtilly mp out the loction of the topologiclly stble nodl lines in the bulk. In ddition, one cn tke dvntge of the fct tht n pplied mgnetic field leds to splitting of the zero-energy pek. Agin, this splitting is strongly dependent on the orienttion of the mgnetic field xis with respect to the nodl lines. Another possibility is to use sptilly resolved STS to investigte the SDOS in the presence of impurities on the surfce. It is expected tht surfce impurities will led to strong sptil modultions of the SDOS, which might give some informtion bout the topologicl chrcteristics of the nodl lines in the bulk. Stbility of surfce modes The zero-energy surfce flt bnds in time-reversl symmetric non-centrosymmetric SCs re topologiclly protected ginst the opening of gp nd re therefore stble ginst wek symmetry preserving deformtions. Conversely, ny perturbtion tht leds to gp opening of the surfce sttes is expected to be ccompnied by the breking of the symmetries of the time-reversl symmetric SC, i.e., TRS or certin types of trnsltionl invrince. One possible scenrio, for exmple, is tht interctions might led to spontneous TRS breking t the boundry of the SC, such s to the coexistence of TRS breking nd TRS preserving order prmeters ner the surfce. This would be observble in experiments, for instnce, s splitting of the zero-bis conductnce pek. In conclusion, using three different topologicl invrints, we exmined the topologicl properties of generl 3D noncentrosymmetric superconductors with TRS. We showed tht in nodl superconducting phses there lwys pper dispersionless Andreev surfce bnds. We estblished correspondence between these zero-energy surfce flt bnds nd the topologiclly protected nodl lines in the bulk, thereby reveling the topologicl origin of the surfce flt bnd. In prticulr, we demonstrted tht the projection of the nodl lines on the surfce coincides with the boundry of the surfce flt bnd. We emphsize tht the presented formlism (or generliztion thereof cn be pplied to ny 3D unconventionl SC tht preserves TRS. One prticulrly interesting fmily of compounds is Li Pd x Pt 3 x B. In these SCs the substitution of Pd by Pt seems to be relted to the reltive strength of singlet nd triplet piring sttes [5]. Hence, it might be possible to observe in Li Pd x Pt 3 x B the trnsition between two topologiclly distinct quntum phses s function of Pt concentrtion. Acknowledgments The uthors thnk B. Béri, A. Furuski, L. Klm, A. Ludwig, R. Nki, P. Horsch, nd M. Sigrist for discussions. A.P.S. is grteful to the Aspen Center for Physics for hospitlity during the preprtion of this work. S.R. is supported by Center for Condensed Mtter Theory t UCB. [] M. König, S. Wiedmnn, C. Brüne, A. Roth, H. Buhmnn, L. W. Molenkmp, X-L. Qi, S-C. Zhng, Science 38, 766 (7; D. Hsieh, D. Qin, L. Wry, Y. Xi, Y. Hor, R. Cv, nd M. Hsn, Nture 45, 97 (8. [] M. Z. Hsn nd C. L. Kne, Rev. Mod. Phys. 8, 345 (; X.-L. Qi nd S.-C. Zhng, rxiv:8.6 (unpublished. [3] A. P. Schnyder, S. Ryu, A. Furuski, A. W. W. Ludwig, Phys. Rev. B (8; AIP Conf. Proc. 34 (9. [4] A. Y. Kitev, AIP Conf. Proc. 34 (9. [5] S. Ryu, A. P. Schnyder, A. Furuski, A. W. W. Ludwig, New Journl of Physics 65 (. [6] S. Murkw, Y. Tmur, Y. Wd, M. Wsi, M. Sitoh, Y. Aoki, R. Nomur, Y. Okud, Y. Ngto, M. Ymmoto, S. Higshitni, nd K. Ngi, Phys. Rev. Lett. 3, 553 (9. [7] Topologicl SCs in clss DIII my lso be found mong centrosymmetric SCs with triplet piring, see, e.g., G.. Volovik nd L. P. Gorkov, Sov. Phys. JTP 6, 843 (985. [8]. Buer, G. Hilscher, H. Michor, Ch. Pul,. W. Scheidt, A. Gribnov, Yu. Seropegin, H. Noël, M. Sigrist, nd P. Rogl, Phys. Rev. Lett. 9, 73 (4. [9] K. Togno, P. Bdic, Y. Nkmori, S. Orimo, H. Tkey, nd K. Hirt, Phys. Rev. Lett. 93, 474 (4; P. Bdic, T. Kondo, nd K. Togno, J. Phys. Soc. Jpn. 74, 4 (5. [] P. A. Frigeri, D. F. Agterberg, A. Kog, nd M. Sigrist, Phys. Rev. Lett. 9, 97 (4. [] K. V. Smokhin, Annls of Physics (9. [] See supplementry mteril for detils of the clcultion. [3] B. Béri, Phys. Rev. B 8, 3455 (. [4] M. Sto, Phys. Rev. B (6. [5] X.-L. Qi, T. L. Hughes, S. Rghu, nd S-C. Zhng, Phys. Rev. Lett. 87 (9. [6] A. B. Vorontsov, I. Vekhter, nd M. schrig, Phys. Rev. Lett., 73 (8. [7] M. schrig, C. Iniotkis, nd Y. Tnk, rxiv:.486.

5 5 [8] S. Ryu nd Y. Htsugi, Phys. Rev. Lett. 89, 77 (. [9] For similr Mjorn Andreev stte t edges of D nodl superconductors, see M. Sto nd S. Fujimoto, Phys. Rev. Lett. 5, 7 (. [] Y. Tnk nd S. Kshiwy, Phys. Rev. Lett. 74, 345 (995. [] S. Kshiwy nd Y. Tnk, Rep. Prog. Phys. 63, 64 (. [] Y. Tnk, Y. Mizuno, T. Yokoym, K. Yd, nd M. Sto Phys. Rev. Lett. 5, 97 (. [3] K. Yd, M. Sto, Y. Tnk, nd T. Yokoym Phys. Rev. B 83, 6455 (. [4] M. Sto, Y. Tnk, K. Yd, nd T. Yokoym, Phys. Rev. B 83, 45 (. [5] H. Q. Yun, D. F. Agterberg, N. Hyshi, P. Bdic, D. Vndervelde, K. Togno, M. Sigrist, nd M. B. Slmon, Phys. Rev. Lett. 97, 76 (6. Supplementry Mterils We first discuss bsic symmetry properties of superconductors with time-reversl invrince nd then go on to derive the topologicl numbers (4, (5, nd (6 from the min text. We shll keep the nlysis s generl s possible, such tht it my be pplied to rbitrry superconducting systems. In Section D we will then specilize to the Bogoliubov-de Gennes Hmiltonin ( describing single-bnd non-centrosymmetric superconductor. Appendix A: Symmetries of the Bogoliubov-de Gennes Hmiltonin Let us consider generl time-reversl invrint superconductor belonging to symmetry clss DIII in the terminology of Refs. [3, 6, 7] ( h(k (k H(k = (k h T, (A ( k with the N-bnd norml stte Hmiltonin h(k nd the superconducting gp mtrix (k, which obeyes (k = T ( k becuse of Fermi sttistics. With N bnds (orbitls nd two spin degrees of freedom, the totl dimension of the Bogoliubov-de Gennes Hmiltonin t momentum k is 4N 4N. A clss DIII superconductor stisfies two independent nti-unitry symmetries: time-reversl symmetry T = KU T, with T =, nd prticle-hole symmetry C = KU C, with C = +. Here, K stnds for the complex conjugtion opertor. Time-reversl symmetry constrins H(k s U T H ( ku T = +H(k, (A withu T = dig(u T,u T ndu T is N N unitry mtrix tht implements time-reversl invrince of the norml stte Hmiltonin, i.e., u T h ( ku T = h(k nd ut T = u T. Observe, q. (A implies u T (k = (ku T. Prticlehole symmetry cts on the Bogoliubov-de Gennes HmiltoninH(k s U C H ( ku C = H(k, (A3 whereu C = σ ½ N,σ,,3 stnd for the three Puli mtrices, nd ½ N is the N N unit mtrix. Combining timereversl nd prticle-hole symmetry we obtin third discrete symmetry, which is given by U S H(kU S = H(k, (A4 with U S = iu T U C. In other words, there is unitry mtrixu S tht nticommutes withh(k nd thereby endows the Hmiltonin with chirl structure. Nmely, in the bsis in whichu S is digonl,h(k tkes block off-digonl form ( D(k H(k = VH(kV = D, (A5 (k where the unitry trnsformtion V is given by V = ( ½N +iu T, (A6 ½ N iu T nd the block off-digonl component redsd(k = h(k+ i (ku T. In the off-digonl bsis the unitry mtrix U T is given s ŨT = VU T V T = σ u T. Thus, time-reversl symmetry cts on D(k s follows u T D T ( ku T = D(k. (A7 To compute the invrint, q. (6, it is dvntgeous to perform second bsis trnsformtion which brings U T into the simple form ŪT = WU T W T = iσ ½ N. This cn be chieved with the help of the unitry mtrix W = dig(½ N,u T. In this new bsis the Bogoliubov-de Gennes Hmiltonin reds ( H(k = W H(kW D(k = D, (A8 (k with the block off-digonl component D(k = D(ku T = h(ku T i (k. (A9 We note tht time-reversl symmetry opertes on D(k s D T ( k = D(k. Appendix B: Flt Bnd Hmiltonin nd Winding Number For the derivtion of the topologicl invrints it is convenient to dibticlly deformh(k, q. (A, into flt bnd Hmiltonin Q(k. The only ssumptions tht we need for computing Q(k re: (i the Hmiltonin hs full spectrl gp nd (ii there is unitry mtrix U S nticommuting with H(k. Thus, the following derivtion of Q(k is pplicble to ny chirl symmetric Hmiltonin with full bulk gp, in prticulr lso to the three-dimensionl topologicl superconductors in symmetry clss AIII, DIII, nd CI [3, 5, 8 3]. In wht follows, we work in bsis in which H(k tkes block off-digonl form. The flt bnd Hmiltonin is defined in terms of the projection opertor P(k which projects onto filled Bloch eigensttes ofh(k t given momentumk. The

6 6 projector P(k, in turn, is defined in terms of the eigenfunctions of H(k ( ( ( D(k χ ± (k χ ± D (k η ± = ±λ (k (k (k η ±, (B (k where =,...,N is the combined bnd nd spin index. We ssume there is spectrl gp round zero energy with λ (k >, nd for definitiveness we choose λ (k > for ll. Multiplying eqution (B from the left by H(k yields DD χ ± (k = λ χ±, D Dη ± = λ η±. (B Hence, the eigenfunctions(χ ±,η ± cn be obtined from the eigenvectors ofdd ord D DD u = λ u, D Dv = λ v. (B3 The eigenvectorsu,v re tken to be normlized to one, i.e., u u = v v =, for ll (here, the index is not summed over. The eigensttes of D D follow from the eigensttes of DD vi v = N D u, (B4 with the normliztion fctor N. Using q. (B3 one cn check tht v is indeed n eigenvector ofd D, D Dv = D D(N D u = N λ D u = λ v, (B5 for ll. The normliztion fctorn is given by xmples of topologicl insultors nd superconductors with completely degenerte bnds re the Dirc representtives of Ref. [5]. The integer-vlued topologicl invrint chrcterizing topologicl superconductors is now simply given by the winding number of q(k. It cn be defined in ny odd sptil dimension. In three dimensions we hve ν 3 = d 3 k 4π εµνρ Tr [ (q µ q(q ν q(q ρ q ], nd in one sptil dimension it reds ν = dktr [ q k q ]. πi (B (B3 Alterntively, it is lso possible to define the winding number in terms of the unflttened off-digonl block D(k of the Hmiltonin. For exmple, for the winding number in one sptil dimension this reds ν = 4πi = π Im dktr [ D k D {D } k D ] dktr[ k lnd]. (B4 u DD u = λ u u = λ N = λ, (B6 Appendix C: Z Invrint for Symmetry Clss DIII for ll. It follows tht the eigenfunctions of H(k re ( χ ± = ( u = ( u ±v ±D. (B7 u /λ η ± With this, the projector P(k onto the filled Bloch sttes becomes P = ( u (u v v (B8 = ( ½N ( u v ½ N v u. Finlly, we obtin for the flt bnd Hmiltonin Q, which is defined sq = ½ 4N P [3], Q = ( u v v u = ( u u D λ D λ u u. (B9 In other words, the off-digonl block ofq(k reds q(k = λ (k u (ku (kd(k, (B where u (k denotes the eigenvectors of DD. For system with completely degenerte bnds, λ = λ, for ll, the bove formul simplifies to q(k = λ(k u (ku (kd(k = D(k. (B λ(k In this section we compute the Z topologicl invrint for symmetry clss DIII ind = ndd = sptil dimensions. It is most convenient to perform this derivtion in the bsis (A8, in which the 4N 4N Bogoliubov-de Gennes Hmiltonin tkes the form H(k = ( D(k D (k, D(k = D T ( k. (C In this representtion, the time-reversl symmetry opertor is given by T = KU T = Kiσ ½ N nd the flt bnd Hmiltonin reds Q(k = ( q(k q (k, q(k = q T ( k. (C The presence of time-reversl symmetry llows us to define the Kne-MeleZ invrint [5, 33 37], W = K Pf [w(k] det[w(k], (C3 with K time-reversl invrint momentum ndpf the Pfffin of n nti-symmetric mtrix. Here, w(k denotes the sewing mtrix w b (k = u + ( k T u + b (k, (C4 where,b =,...,N nd u ± (k is the -th eigenvector of Q(k with eigenvlue ±. The Pfffin is n nlog of

7 7 the determinnt tht cn be defined only for n n ntisymmetric mtrices A. It is given in terms of sum over ll elements of the permuttion groups n Pf(A = n n sgn(σ A σ(i,σ(i. n! σ S n i= Due to the block off-digonl structure of q. (C set of eigen Bloch functions ofq(k cn be constructed s [5] or, lterntively, s u ± (k N = ( n ±q, (C5 (kn u ± (k S = ( ±q(kn, (C6 n where n re N momentum independent orthonorml vectors. For simplicity we choose (n b = δ b. In pssing, we note tht both u ± (k N nd u ± (k S re well-defined globlly over the entire Brillouin zone. To compute the topologicl number we choose the bsis u ± (k N. Combining qs. (C4 nd (C5 yields w b (k = n (, n q( k ( nb iσ ½ N K q (kn b = ( n, n q( k ( q T (kn b n b = ( n q T (kn b n q( kn b = qb T (k. (C7 In the second lst line we used q. (C, i.e., q( k = q T (k. In conclusion, the Z topologicl number in sptil dimensionsd = ndd = is given by W = Pf [ q T (K ], (C8 K det[q(k] where K denotes the four (two time-reversl invrint moment of the two-dimensionl (one-dimensionl Brillouin zone. Appendix D: Flt Bnd Hmiltonin for the Non-centrosymmetric Superconductor, q. ( Let us now pply the formlism developed in the preceding sections to the Bogoliubov-de Gennes Hmiltonin ( from the min text, describing single bnd noncentrosymmetric superconductor [38, 39]. First, we note tht time-reversl symmetry for H(k, q. (, is implemented by U T H ( ku T = +H(k with U T = σ iσ. Hence, we need to set u T = iσ in q. (A. It then follows from q. (A6 tht H(k cn be brought into block off-digonl form by the unitry trnsformtion V = ( ½ σ. (D ½ +σ ( (b M M + FIG. 5. ( Bnd structure of Hmiltonin ( for the point group O in slb geometry with ( fce s function of surfce momentum k with k =, i.e., long the dshed line in pnel (b. Here, we set (g, s = (.3,.5. (b Brillouin zone of the ( surfce with the vlues of W ( (K, q. (D, t the four time reversl invrint moment K {Γ,M,M,M 3}. The trnsformed Hmiltonin is given by ( D(k VH(kV = D (k + Γ + M 3, (D with the off-digonl block ( Bk +Al D(k = k z A(lk x ily k A(lk x +ily k B k Alk z = B k σ +Al k σ, (D3 nd where we hve introduced the short-hnd nottion A = α+i t, B k = ε k +i s. (D4 Alterntively, we cn lso choose to work in the bsis (A8, in which cse the off-digonl component reds ( A(l x D(k = k il y k B k Alk z B k Alk z A(lk x +ily k. (D5 For the computtion of the flt bnd Hmiltonin Q(k it is, however, more convenient to use q. (D3. Repeting the steps of section B, we clculte the eigenvectorsu (k of D k D k = A l k + B k +(AB k +B k A l k σ, (D6 where l k = l k. The eigenfunctionsu (k of D k D k cn be obtined by digonlizing l k σ. Hence, when (lk x,ly k (,, we find tht the eigenvectorsu (k re given by ( l x u / (k = k il y k lk (l k lk z ±l k lk z. (D7 According to q. (B, the off-digonl block of the flt bnd Hmiltonin Q(k is defined in terms of the eigenvectorsu (k. Thus, we need to compute u λ (ku (k (D8 =, k [ = (λ k +λ k σ +(λ k λ k l ] k σ, λ k λ k l k

8 8 with the two positive eigenvlues λ k = B k Al k nd λ k = B k + Al k. Note tht the lst term in the second line of q. (D8 contins removble singulrities t the points k where l k =. For those points in the Brillouin zone one needs to crefully tke the limit k k to obtin the correct vlue of q. (D8. Finlly, by use of q. (B together with qs. (D3 nd (D8 we find for the off-digonl block of the flt bnd Hmiltonin [ q(k = {Al k (λ k λ k +B k (λ k +λ k }σ λ k λ k +{Al k (λ k +λ k +B k (λ k λ k } l k l k σ ]. (D9 Now, for the invrint we need to bring q(k into the bsis in whichu T = iσ ½. This is chieved by letting Using q. (C8 we get W = K q(k iq(kσ. Pf [ iσ q T (K ] det[iσ q T (K] = K (D B K B K, (D where we hve mde use of the fct thtl k is n ntisymmetric function, i.e.,l k = l k.. surfce stte As discussed in the min text, the number (D cn be used to deduce the presence of Andreev surfce sttes t timereversl invrint moment of the surfce. To exemplify this, let us consider Hmiltonin ( in slb geometry with (lmn fce. At the four time-reversl invrint moment K of the(lmn surfce the invrint is defined by W (lmn (K = Pf [ iσ q T (K,K ] K det [ ]. (D iσ q T (K,K q. (D is quntized to+ or, withw (lmn (K = indicting the presence of Krmers degenerte surfce modes t the surfce momentumk (see Fig. 5. [6] M. R. Zirnbuer, J. Mth. Phys. 37, 4986 (996. [7] A. Altlnd nd M. R. Zirnbuer, Phys. Rev. B 55, 4 (997. [8] G.. Volovik, in The Universe in Helium Droplet, The Interntionl Series of Monogrphs on Physics Vol. 7 (Oxford University Press, New York, 3; G.. Volovik, in xotic Properties of Superfluid 3He, Series I N Modern Condensed Mtter Physics Vol. (World Scientific, Singpore, 99. [9] R. Roy, rxiv: (unpublished. [3] X.-L. Qi, T. Hughes, nd S.-C. Zhng, Phys. Rev. B 78, 9544 (8. [3] A. P. Schnyder, S. Ryu, nd A. W. W. Ludwig, Phys. Rev. Lett., 9684 (9. [3] A. P. Schnyder, P. M. R. Brydon, D. Mnske, nd C. Timm, Phys. Rev. B 8, 8458 (. [33] C. Kne nd. Mele, Phys. Rev. Lett. 95, 68 (5. [34] C. Kne nd. Mele, Phys. Rev. Lett. 95, 468 (5. [35] L. Fu nd C. L. Kne, Phys. Rev. B (7. [36] J.. Moore nd L. Blents, Phys. Rev. B (7. [37] R. Roy, Phys. Rev. B 79, 953 (9. [38] P. A. Frigeri, D. F. Agterberg, A. Kog, nd M. Sigrist, Phys. Rev. Lett. 9, 97 (4. [39] K. V. Smokhin nd V. P. Mineev, Phys. Rev. B 77, 45 (8.