Spin-Orbit-Free Topological Insulators

Spin-Orbit-Free Topological Insulators

Abstract

In this lecture for the Nobel symposium, we review previous research on a class of translational-invariant insulators without spin-orbit coupling. These may be realized in intrinsically spinless systems such as photonic crystals and ultra-cold atoms. Some of these insulators have no time-reversal symmetry as well, i.e., the relevant symmetries are purely crystalline. Nevertheless, topological phases exist which are distinguished by their robust surface modes. To describe these phases, we introduce the notions of (a) a halved-mirror chirality: an integer invariant which characterizes half-mirror planes in the 3D Brillouin zone, and (b) a bent Chern number: the traditional Thouless-Kohmoto-Nightingale-Nijs invariant generalized to bent 2D manifolds. Like other well-known topological phases, their band topology is unveiled by the crystalline analog of Berry phases, i.e., parallel transport across certain non-contractible loops in the Brillouin zone. We also identify certain topological phases without any robust surface modes – they are uniquely distinguished by parallel transport along bent loops, whose shapes are determined by the symmetry group. Finally, we describe the Weyl semimetallic phase that intermediates two distinct, gapped phases.

I Introduction

Insulating phases are deemed distinct if they cannot be connected by continuous changes of the Hamiltonian that preserve both the energy gap and the symmetries of the phase; in this sense we say that the symmetry protects the phase. Distinct phases have strikingly different properties – some nontrivial phases are characterized by robust boundary modes which are particularly accessible through experimental techniques, e.g., STM and ARPES. For many well-known topological insulators (TI’s), the existence of boundary bands is in one-to-one correspondence with the topology of the bulk wavefunctions.fidkowski2011 () Examples include the Chern insulator and the quantum spin Hall insulator; these non-interacting insulators fall under ten well-known symmetry classes which are distinguished by time-reversal, particle-hole and chiral symmetries.schnyder2009 () Given the completeness of this classification, attention has shifted to identifying topological phases which rely on other symmetries.chenfang2012 (); chenfang2014a (); shiozaki () The symmetries which are ubiquitous in condensed matter are the crystal space groups; among them the point groups involve transformations that preserve a spatial point. For example, the reflection symmetry plays an important role in SnTe,Hsieh2012 (); Xu2012 (); tanaka2012 () where the mirror Chern numberteo2008 () characterizes planes in the 3D Brillouin zone which are invariant under reflection; in short, we call these mirror planes. The Bloch wavefunctions in each mirror plane may be decomposed according to their representations under reflection, and each subspace may exhibit a quantum anomalous Hall effect.Haldane1988 ()

Including SnTe, all experimentally-realized TI’s have thus far been strongly spin-orbit-coupled, resulting in a variety of exotic phenomenon, e.g., Rashba spin-momentum locking on the surface of a TI.winklerbook () Considerably less attention has been addressed to spin-orbit-free systems, i.e., electronic insulators and semimetals in which spin-orbit coupling is negligibly weak, as well as intrinsically spinless systems such as photonic crystals and ultra-cold atoms. The topological classifications of spin-orbit-free and spin-orbit-coupled systems generically differ. For example, the mirror Chern numbers are trivially zero for any spin-orbit-free system with the same crystalline symmetries as SnTe.aris2013a () The first proposal of a spin-orbit-free TI with robust surface modes relies on a combination of -fold rotational and time-reversal symmetries, for and .fu2011 () Henceforth, we refer to this combined group as , where . Topological phases in this symmetry class are sometimes distinguished by their robust surface modes, but not always. A more sensitive probe of their band topology is through the crystalline analog of Berry phases, i.e., parallel transport across certain non-contractible loops in the Brillouin zone. As we elaborate in Sec. III, the topological phases without surface modes are uniquely distinguished by parallel transport along bent loops, whose shapes are determined by the symmetry group.aris2014a () Besides parallel transport,aa2012 () the entanglement spectrum is also known to distinguish topological phases without boundary modes, in various symmetry classes.turner2010A (); hughes2011 (); turner2012 (); chenfang2013a ()

As with many other well-known TI’s,fu2006 (); roy2009 () time-reversal symmetry (TRS) plays an essential role in protecting the phase. Can one do away with TRS entirely? Indeed, it is now known that spin-orbit-free topological phases can exist without TRS; the relevant symmetries are purely crystalline. These phases include the first TI protected by symmorphic point groups,aris2013a () as well as the first TI protected by non-symmorphic space groups.chaoxingliu2013 (). In this short review, we establish sufficient criteria for a symmetry group to protect robust boundary modes. By applying this criteria to the 32 crystallographic point groups, we exhaustively identify the groups, for and , as being able to support robust surface modes without additional symmetries. For example, the symmetries of are described in Fig. 1(a). To describe these systems, we introduce the notions of (a) a halved-mirror chirality: an integer invariant which characterizes half-mirror planes in the 3D Brillouin zone (BZ), and (b) a bent Chern number: the traditional TKNN invariantthouless1982 () generalized to bent 2D manifolds (illustrated in Fig. 1(b)). We find that a Weyl semimetallic phase intermediates two gapped phases with distinct halved chiralities. In Sec. II of this review, we characterize the systems, leaving a more comprehensive treatment of the other groups to Ref. aris2013a, .

Figure 1: (a) Top-down view of tetragonal BZ with symmetry; our line of sight is parallel to the rotational axis. Reflection-invariant planes are indicated by solid lines. Except the line through , all non-equivalent -invariant lines are indicated by circles for , and squares for . Bottom of (b): Half-mirror planes (HMP’s) in the 3D Brillouin zone (BZ) of a tetragonal lattice with symmetry. Red: HMP. Blue: HMP. We define a bent Chern number on the triangular pipe with its ends identified. Top of (b): Non-black lines are half-mirror lines (HML’s) in the corresponding 2D BZ of the 001 surface; each HML connects two distinct -invariant points.

Ii insulators and semimetals

Given that the mirror Chern number vanishes in spin-orbit-free systems with symmetry, how else can surface modes be made robust? Progress is made by identifying two -invariant points ( and ) in the (001) surface BZ; each -invariant point is mapped to itself under a four-fold rotation, up to translations by a reciprocal lattice vector, as we illustrate in Fig. 1(a). In addition, the black lines that intersect (resp. ) in Fig. 1(a) are invariant under reflection, indicating that there are additional symmetries at this point, e.g., the reflection maps . The combined symmetries of (resp. ) form the little grouptinkhambook () . At these points, the symmetry enforces that and orbitals are doubly-degenerate. That is, each pair of orbitals transforms in the two-dimensional irreducible representation (irrep) of . Indeed, if we choose a basis () that diagonalizes with eigenvalue , clearly transforms and vice versa. We then propose that surface bands of systems assume topologically distinct structures on the red line which connects and , as illustrated in Fig. 1(b). We call the red line a half-mirror line (HML) because it does not close into a circle (modulo reciprocal lattice vectors), yet all Bloch wavefunctions on the red line may be diagonalized by a single reflection operator (). Suppose we parametrize HML with , where (resp. ) at (resp. ). At and , we have just established that (001) surface bands (with character) form doubly-degenerate pairs with opposite mirror eigenvalues, irrespective of whether the system has TRS. For representations without spin, we may label bands with mirror eigenvalue as mirror-even (mirror-odd). Given these constraints at and , there are ways to connect mirror-even bands to mirror-odd bands, as illustrated schematically in Fig. 2; it should be noticed that band intersections do not split as they correspond to bands with opposite mirror eigenvalues. We define the halved mirror chirality as the difference in number of mirror-even chiral modes with mirror-odd chiral modes; if , the surface bands robustly interpolate across the energy gap. may be easily extracted by inspection of the surface energy-momentum dispersion: first draw a constant-energy line within the bulk energy gap and parallel to the HML, e.g., the blue line in Fig. 2. At each intersection with a surface band, we calculate the sign of the group velocity , and multiply it with the eigenvalue under reflection . Finally, we sum this quantity over all intersections along HML to obtain . In Fig. 2, we find two intersections as indicated by red squares, and . The -classification of (001) surface bands relies on two-dimensional irreps in the surface BZ; on surfaces which break symmetry, the surface bands transform in one-dimensional irreps, and cannot assume topologically distinct structures.

Figure 2: Distinct connectivities of the (001) surface bands along the half-mirror-lines. Black solid (dotted) lines indicate surface bands with eigenvalue () under reflection ; crossings between solid and dotted lines are robust due to reflection symmetry. For simplicity, we have depicted all degeneracies at momenta and as dispersing linearly with momentum. For , such crossings are in reality quadratic.aris2013a ()

In addition to HML, the blue line in Fig. 1(b) is a second half-mirror line (HML), and states along HML are characterized by an independent invariant . HML is the union of two mirror lines, and , which are related by a rotation. States along are invariant under the reflection , while along the relevant reflection is . The product of these orthogonal reflections is a rotation () about , thus . Since both orbitals are odd under a -rotation, an eigenstate of with eigenvalue also diagonalizes with eigenvalue , e.g., is even under but odd under . In effect, each Bloch wavefunction along HML can be labelled by the quantum number of a single reflection operator . The second invariant is then the difference in number of mirror-even () chiral modes with mirror-odd () chiral modes, in a straightforward generalization of .

Thus far we have described the halved chirality as a topological property of surface bands along HML, but we have not addressed how is encoded in the bulk wavefunctions. Taking to lie along the rotational axis, each HML in the surface BZ is the -projection of a half-mirror-plane (HMP) in the 3D BZ, as illustrated in Fig. 1(b). Each HMP connects two distinct -invariant lines, and all Bloch wavefunctions in a HMP may be diagonalized by a single reflection operator. HMP is parametrized by and , where along the first (second) -invariant line. Then the halved mirror chirality has the following gauge-invariant expression by bulk wavefunctions:aris2013a ()

(1)

() is defined as the Berry curvatureberry1984 (); zak1982 (); zak1989 () of occupied bands as contributed by the mirror-even (-odd) subspace.

Figure 3: (a) to (d) illustrate how the halved chirality of a HMP may change. The direction of arrows indicate whether Weyl nodes are created or annihilated. + (-) labels a Berry monopole with positive (negative) charge; () labels a crossing in the HMP between mirror-even (-odd) bands. (e) We provide a top-down view of the trajectories of Weyl nodes, in the transition between two distinct gapped phases. Our line of sight is parallel to . Black lines indicate mirror faces, while colored lines specially indicate HMP’s, with the same color legend as in Fig. 1(b).

How do transitions occur between gapped phases with different ? The invariants are well-defined so long as bulk states in the HMP’s are gapped, which is true of insulators. These invariants may also be used to characterize semimetals, so long as the gaps close away from the HMP’s. When gaps close, the band touchings are generically Weyl nodes,Murakami2007B (); wan2010 () though exceptions exist with a conjunction of time-reversal and inversion symmetries.burkov2011 () The chirality of each Weyl node is its Berry charge, which is positive (negative) if the node is a source (sink) of Berry flux. By the Nielsen-Ninomiya theorem, the net chirality of all Weyl nodes in the BZ is zero.NIELSEN1981 () To make progress, we divide the BZ into ‘unit cells’, such that the properties of one ‘unit cell’ determine all others by symmetry. As seen in Fig. 1(b), each ‘unit cell’ resembles the interior of a triangular pipe; it is known as the orbifold . The net chirality of an orbifold can be nonzero, and is determined by the Chern number on the 2D boundary of the orbifold. As the boundary resembles the surface of a triangular pipe, we call the invariant a bent Chern number (). systems are thus described by three invariants , which are related by If is odd, must be nonzero due to an odd number of Weyl nodes within the orbifold, which implies the system is gapless. If the system is gapped, must be even. However, the converse is not implied. This parity constraint may be understood in light of a Weyl semimetallic phase that intermediates two gapped phases with distinct halved chiralities. There are four types of events that alter the halved chirality of HMP; we explain how Weyl nodes naturally emerge in the process. (i) Suppose the gap closes between two mirror-even bands in HMP. Around this band touching, bands disperse linearly within the mirror plane, and quadratically in the direction orthogonal to the plane. Within HMP, the linearized Hamiltonian around the band crossing describes a massless Dirac fermion in the even representation of reflection. If the mass of the fermion inverts sign, changes by , implying that also changes by through Eq. (1). This quantized addition of Berry flux is explained by a splitting of the band-touching into two Weyl nodes of opposite chirality, and on opposite sides of the mirror face (Fig. 3(a)). In analogy with magnetostatics, the initial band touching describes the nucleation of a dipole, which eventually splits into two opposite-charge monopoles; the flux through a plane separating two monopoles is unity. (ii) The same argument applies to the splitting of dipoles in the mirror-odd subspace, which alters by , and by . For (iii) and (iv), consider two opposite-charge monopoles which converge on HMP and annihilate, causing to change by unity. The sign of this change is determined by whether the annihilation occurs in the mirror-even or odd subspace. The transition between two distinct gapped phases is then characterized by a transfer of Berry charge between two distinct HMP’s (Fig. 3(e)). In the intermediate semimetallic phase, the experimental implications include Fermi arcs on the (001) surface.wan2010 (); halasz2012 (); chen2012 ()

Iii insulators and semimetals

The topological connectivity of surface bands in systems relied on a criterion which may be generalized: all we need are two high-symmetry points on the boundary BZ, which admit multi-dimensional irreps at each point. Suppose we relaxed the reflection symmetries in the group , and now added time-reversal symmetry. The resultant group () also has two-dimensional irreps at -invariant momenta and . For example, the orbitals are eigenstates of four-fold rotation, while time-reversal maps . This two-fold degeneracy in a spin-orbit-free system is analogous to the Kramer’s degeneracy in spin-orbit coupled systems. Similarly, there exists two distinct gapped phases which are distinguished by robust surface modes;fu2011 () we refer to one as the trivial phase and the other a strong topological phase. Our case study of the insulator is motivated by two issues: (i) It is unclear if these two phases are physically distinguishable if we experimentally probe the bulk instead of the surface. (ii) In this review, we highlight the existence of a third topological phase which does not manifest robust surface modes. Does this ‘weak phase’ have any physical consequence? One answer to (i) and (ii) may be found through holonomy, i.e., parallel transport along certain non-contractible loops in the BZ.leone2011 (); ekert2000 (); recati2002 () An electron transported around a loop acquires a Berry-Zak phase,berry1984 (); zak1982 (); zak1989 () which has recently been measured by interference in cold atom experiments.atala2013 () For the purpose of unveiling the bulk topology of all three phases, we find that not all non-contractible loops work. Straight loops are commonly studied in the geometric theory of polarization, due to their relation with Wannier functions;kingsmith1993 (); Taherinejad2013 () however, these loops cannot identify our weak phase. Instead, we propose that all three phases are distinguished by parallel transport along bent loops, whose shapes are determined by the symmetry group – they are illustrated in Fig. 4(a) for the group, while the case is discussed in Ref. aris2014a, .

Figure 4: Bottom of (a) illustrates certain bent loops in the 3D Brillouin zone of a simple tetragonal lattice. Top of (a): (001)-surface BZ of the same lattice. (b) and (c): Berry-phase spectrum of trivial and strong phases respectively, for the model of Eq. (2).

It is known from Ref. fu2011, that the insulator is characterized by two indices , for ; describes the Bloch wavefunctions in the plane . In analogy with the spin-orbit-coupled insulator,fu2011 (); fu2007b (); roy2009a () shall be referred to as weak indices. (resp. ) corresponds to the trivial (resp. weak) phase, and describes a strong phase. The product is a strong index that determines the absence or presence of robust surface modes on the 001 surface; we take to lie along the principal axis. We give these weak indices a physically transparent interpretation from the perspective of holonomy. For illustration, we consider a model on a tetragonal lattice that is composed of two interpenetrating cubic sublattices. Our tight-binding basis comprises of orbitals, which transform in the two-dimensional irrep of . The Bloch Hamiltonian is

(2)

where , , , , and . In , and are Pauli matrices for , while and are identities in each 2D subspace. label orbitals on one sublattice, and label orbitals on the other. This Hamiltonian is four-fold symmetric: , and time-reversal symmetric: . The ground state of (2) comprises its two lowest-lying bands. The phase diagram of this model is plotted in Fig. 5(a) for different parametrizations of Eq. (2).

Figure 5: (a) Phase diagram of model (2), as a function of parameters and . Blue (uncolored) regions correspond to gapped (gapless Weyl) phases.aris2013a () The weak indices in each gapped phase are indicated by . The blue square in the center is approximately bound by and . The 001-surface spectrum is plotted for two representative points on the phase diagram: (b) for , and (c) for . , and are high-symmetry momenta defined in Fig. 4(a).

To probe the bulk topology, we perform parallel transport along a bent loop that connects two -invariant points. We define as the loop connecting in the plane, and connects in the plane. They are respectively depicted by red and brown lines in Fig. 4(a). Let us denote the eigenstates of the Bloch Hamiltonian by , where is a band index. The matrix representation of holonomy is known as the Wilson loop, and it is the path-ordered exponential of the Berry-Wilczek-Zee connection : wilczek1984 (); berry1984 () Here, denotes a loop and is a matrix with dimension equal to the number () of occupied bands. The gauge-invariant spectrum of is also known as the Berry-phase factors ().zak1982 (); zak1989 (); aa2012 () We show that the spectrum of encodes the weak index . If we define is the number of eigenvalues in the spectrum of , we find that this number is always even. Furthermore, the weak indices are related to by

(3)

This weak index is equivalent to an alternative formulation in Ref. fu2011, , where it is expressed as an invariant involving the Pfaffian of a matrix.aris2014a () Our formulation through holonomy reveals a geometric connection with the theory of matrices in : as we elaborate in Ref. aris2014a, , the parity of specifies one of two classes of a special rotation, which is in one-to-one correspondence with two sectors of ground states in the plane. Furthermore, we are able to derive the spectrum of for any number of occupied bands.aris2014a () In a minimal model with , there are only two possibilities (), which are distinguished by the parity of .

Let us demonstrate how to realize both classes of in the model (2). We consider a family of loops in planes of constant , such that is the red line in Fig. 4(a), and all other loops project to in . In the trivial phase (parametrized by ), we find , or equivalently . The eigenvalues of interpolate between (at ) to (at ), as illustrated in Fig. 4(b). The absence of surface modes on the line is demonstrated in Fig. 5(b). In comparison, the strong phase () is characterized by , or ; its surface modes are illustrated in Fig. 5(c). As is varied from to in Fig. 4(c), the Berry phases interpolate across the maximal range . We say that the Berry phases exhibit spectral flow if they interpolate across their maximally-allowed range as we tune a BZ parameter. Spectral flow is a unifying trait shared by many TI’s, including the Chern insulator,Haldane1988 (); thouless1982 (); qi2006 (); ProdanJMP2009 () the quantum spin Hall insulator,yu2011 (); soluyanov2011 () and the inversion-symmetric TI.aa2012 () In comparison, the weak phase is distinguished by , which clearly does not exhibit spectral flow. Finally, we point out that a different choice of Wilson loop can identify the strong index ,Taherinejad2013 () but cannot individually distinguish the weak indices: and .

Iv Conclusion

We propose sufficient criteria for surface modes whose robustness rely on a symmetry group. Minimally, (i) the symmetry must be unbroken by the presence of the surface. Additionally, either (ii) a reflection symmetry exists so that mirror-subspaces can display a quantum anomalous Hall effect, or (iii) there exist at least two high-symmetry points in the surface BZ, which admit higher-than-one dimensional irreps of the symmetry group. In addition to predicting new topological materials, these criteria are also satisfied by the well-known SnTe class,Hsieh2012 () and also the insulators.fu2006 (); kane2005A (); kane2005B (); bernevig2006a (); bernevig2006c (); koenig2007 (); fu2007b (); moore2007 (); roy2009 (); roy2009a (); hsieh2008 () Among the crystallographic point groups, only the and groups are preserved for a surface that is orthogonal to the rotational axis.tinkhambook () Though all groups satisfy (ii), the lack of spin-orbit coupling implies only systems can have nonvanishing mirror Chern numbers.aris2014a () While all groups by themselves only have one-dimensional irreps, the or group satisfies (iii) in combination with TRS, as is known for the topological crystalline insulators introduced in Ref. fu2011, . Finally, only a subset of the groups possess two-dimensional irreps which satisfy (iii): and .aris2013a ()

Acknowledgments

Chen Fang and Matthew Gilbert played an essential role in formulating the theory of insulators, whose short review we have presented here. AA and BAB were supported by Packard Foundation, 339-4063–Keck Foundation, NSF CAREER DMR-095242, ONR - N00014-11-1-0635 and NSF-MRSEC DMR-0819860 at Princeton University. This work was also supported by DARPA under SPAWAR Grant No.: N66001-11-1-4110. CF is supported by ONR-N0014-11-1-0728 for salary and DARPA-N66001-11-1-4110 for travel. MJG acknowledges support from the ONR under grant N0014-11-1-0728 and N0014-14-1-0123, a fellowship from the Center for Advanced Study (CAS) at the University of Illinois.

References

  1. L. Fidkowski, T. S. Jackson,  and I. Klich, Phys. Rev. Lett., 107, 036601 (2011).
  2. A. P. Schnyder, S. Ryu, A. Furusaki,  and A. W. W. Ludwig, AIP Conf. Proc., 1134, 10 (2009).
  3. Chen Fang, Matthew J. Gilbert, and B. Andrei Bernevig, Phys. Rev. B 86, 115112 (2012).
  4. Chen Fang, Matthew J. Gilbert, and B. Andrei Bernevig, Phys. Rev. Lett. 112, 106401 (2014).
  5. K. Shiozaki and M. Sato, Phys. Rev. B 90, 165114 (2014).
  6. T. H. Hsieh, et al., Nat. Commun. 3:982 doi: 10.1038/ncomms1969 (2012).
  7. S.-Y. Xu, et al., Nat. Commun. 3:1192 doi: 10.1038/ncomms2191 (2012).
  8. Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, T. Takahashi, K. Segawa,  and Y. Ando, Nat. Phys., 8, 800 (2012).
  9. J. C. Y. Teo, L. Fu,  and C. L. Kane, Phys. Rev. B, 78, 045426 (2008).
  10. F. D. M. Haldane, Phys. Rev. Lett., 61, 2015 (1988).
  11. R. Winkler, Spin-Orbit Coupling Effects in Two-Dimensional Electron and Hole Systems (Springer Tracts in Modern Physics, 2003).
  12. A. Alexandradinata, C. Fang, M. J. Gilbert,  and B. A. Bernevig, Phys. Rev. Lett., 113, 116403 (2014).
  13. L. Fu, Phys. Rev. Lett., 106, 106802 (2011).
  14. A. Alexandradinata and B. A. Bernevig, arXiv:1409.3236.
  15. A. Alexandradinata, X. Dai,  and B. A. Bernevig, Phys. Rev. B, 89, 155114 (2014).
  16. A. M. Turner, Y. Zhang,  and A. Vishwanath, Phys. Rev. B, 82, 241102 (2010).
  17. T. L. Hughes, E. Prodan,  and B. A. Bernevig, Phys. Rev. B, 83, 245132 (2011).
  18. Ari M. Turner, Yi Zhang, Roger S. K. Mong, and Ashvin Vishwanath, Phys. Rev. B 85, 165120 (2012).
  19. Chen Fang, Matthew J. Gilbert, and B. Andrei Bernevig, Phys. Rev. B 87, 035119 (2013).
  20. L. Fu and C. L. Kane, Phys. Rev. B, 74, 195312 (2006).
  21. R. Roy, Phys. Rev. B, 79, 195321 (2009a).
  22. Chao-Xing Liu, Rui-Xing Zhang, and Brian K. VanLeeuwen, Phys. Rev. B 90, 085304 (2014).
  23. D. J. Thouless, M. Kohmoto, M. P. Nightingale,  and M. den Nijs, Phys. Rev. Lett., 49, 405 (1982).
  24. M. Tinkham, Group Theory and Quantum Mechanics (Dover, 2003).
  25. M. V. Berry, Proc. R. Soc. Lond A, 392, 45 (1984).
  26. J. Zak, Phys. Rev. Lett., 48, 359 (1982).
  27. J. Zak, Phys. Rev. Lett., 62, 2747 (1989).
  28. S. Murakami, New J. Phys., 9, 356 (2007).
  29. X. Wan, A. Turner, A. Vishwanath,  and S. Y. Savrasov, Phys. Rev. B, 83, 205101 (2011).
  30. A. A. Burkov, M. D. Hook,  and L. Balents, Phys. Rev. B, 84, 235126 (2011).
  31. H. B. Nielsen and M. Ninomiya, Nucl. Phys. B, 185, 20 (1981).
  32. G. B. Halasz and L. Balents, Phys. Rev. B, 85, 035103 (2012).
  33. C. Fang, M. J. Gilbert, X. Dai,  and B. A. Bernevig, Phys. Rev. Lett., 108, 266802 (2012).
  34. R. Leone, J. Phys. A: Math. Theor., 44, 295301 (2011).
  35. A. Ekert, M. Ericsson, P. Hayden, H. Inamori, J. A. Jones, D. K. L. Oi,  and V. Vedral, J. Mod. Opt., 47, 2501 (2000).
  36. A. Recati, T. Calarco, P. Zanardi, J. I. Cirac,  and P. Zoller, Phys. Rev. A, 66, 0302309 (2002).
  37. M. Atala, et al., Nature Physics, 9, 795 (2013).
  38. R. D. King-Smith and D. Vanderbilt, Phys. Rev. B, 47, 1651 (1993).
  39. M. Taherinejad, K. F. Garrity,  and D. Vanderbilt, Phys. Rev. B, 89, 115102 (2014).
  40. L. Fu, C. L. Kane,  and E. J. Mele, Phys. Rev. Lett., 98, 106803 (2007).
  41. R. Roy, Phys. Rev. B, 79, 195322 (2009b).
  42. F. Wilczek and A. Zee, Phys. Rev. Lett., 52, 2111 (1984).
  43. X.L. Qi, Y.S. Wu,  and S.C. Zhang, Phys. Rev. B, 74, 045125 (2006).
  44. E. Prodan, J. Math. Phys., 50, 083517 (2009).
  45. R. Yu, X. L. Qi, A. Bernevig, Z. Fang,  and X. Dai, Phys. Rev. B, 84, 075119 (2011).
  46. A. A. Soluyanov and D. Vanderbilt, Phys. Rev. B, 83, 235401 (2011).
  47. C. L. Kane and E. J. Mele, Phys. Rev. Lett., 95, 226801 (2005a).
  48. C. L. Kane and E. J. Mele, Phys. Rev. Lett., 95, 146802 (2005b).
  49. B.A. Bernevig and S.C. Zhang, Phys. Rev. Lett., 96, 106802 (2006).
  50. B. A. Bernevig, T. L. Hughes,  and S.C. Zhang, Science, 314, 1757 (2006).
  51. M. König, S. Wiedmann, C. Brüne, A. Roth, H. Buhmann, L. Molenkamp, X.-L. Qi,  and S.-C. Zhang, Science, 318, 766 (2007).
  52. J. E. Moore and L. Balents, Phys. Rev. B, 75, 121306 (2007).
  53. D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava,  and M. Z. Hasan, Nature, 452, 970 (2008).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
204658
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description