Single particle Green's functions and interacting topological insulators

Show preview

Hide preview

Single particle Green’s functions and interacting topological insulators
V. Gurarie Department of Physics, University of Colorado, Boulder CO 80309
(Dated: November 11, 2010)
We study topological insulators characterized by the integer topological invariant Z, in even and odd spacial dimensions. These are well understood in case when there are no interactions. We extend the earlier work on this subject to construct their topological invariants in terms of their Green’s functions. In this form, they can be used even if there are interactions. Specializing to one and two spacial dimensions, we further show that if two topologically distinct topological insulators border each other, the difference of their topological invariants is equal to the difference between the number of zero energy boundary excitations and the number of zeroes of the Green’s function at the boundary. In the absence of interactions Green’s functions have no zeroes thus there are always edge states at the boundary, as is well known. In the presence of interactions, in principle Green’s functions could have zeroes. In that case, there could be no edge states at the boundary of two topological insulators with different topological invariants. This may provide an alternative explanation to the recent results on one dimensional interacting topological insulators.
PACS numbers: 73.43.-f, 73.20.-i
I. INTRODUCTION
The original proposals establishing existence of topo- logical insulators with the Z2 topological invariant (in the presence of spin-orbit coupling with time- reversal invariance)1,2 and their subsequent observations in experiment3,4 contributed to a start of a tremendous amount of activity on this subject.
Recently a classification table was worked out for all possible noninteracting topological insulators5,6. 10 sym- metry classes of topological insulators were identified. In each spacial dimensionality, representative systems of 5 classes are topologically trivial, while the remaining 5 classes can be topologically nontrivial.
The term “topologically nontrivial” in this context im- plies the following. First of all, a noninteracting topolog- ical insulator, when coupled to an external electromag- netic field, is typically (but not always) expected to have a response to this field given by a Chern-Simons term in even dimensional space or by a θ-term in odd dimensional space5,7. Second, a topologically nontrivial insulator has zero energy excitations at its boundary (or at the bound- ary of two topologically distinct insulators). These are usually referred to as edge states. To derive all of these properties, typically one takes advantage of the fact that these systems are noninteracting fermions and thus are essentially exactly solvable.
A natural question is then whether adding interactions to topological insulators preserves their observable prop- erties, such as the existence and robustness of the edge states.
One particular approach to this problem was developed over 20 years ago in a series of papers by G. Volovik and collaborators, summarized in his book8. This approach is based on the following simple idea. Noninteracting fermionic systems are described by a Hamiltonian which is equivalent to a matrix Hab(k), which is a function of momenta k, as explained in Ref. [5]. This matrix satisfies
certain conditions, which define the symmetry class of the topological insulator. For example, for insulators with the chiral (sublattice) symmetry, the condition reads
H(k) = −Σ†H(k)Σ, (1) where Σ is a certain unitary matrix. Whether a particu- lar system belongs to a non-trivial topological class can be investigated by studying the topological properties of the function H(k) which defines a map from the space of the momenta k to the space of the matrices H.
Interacting systems are no longer described by a matrix of this kind. However, one can define a Green’s function (single particle retarded or advanced Green’s function) for both interacting and noninteracting systems. For noninteracting systems, it takes especially simple form
G(Ω,k) = [Ω−H(k)]−1 . (2) The Green’s function computed at zero frequency is sim- ply equal to minus inverse of H and thus inherits all its symmetries, implying in case of Eq. (1), that
G(0,k) = −Σ†G(0,k)Σ. (3) At this point, we say that an interacting system belongs to a particular symmetry class if its Green’s function G(0,k) satisfies the relations such as Eq. (3). Subse- quently, we can study whether a particular interacting system is topological by studying the function G(0,k), instead of the earlier H(k).
All of this would stay in the realm of pure math, if it were not for a further observation in Ref. [8], on the example of symmetry class D systems in 2 dimensions (class D is a shorthand for a superconductor without any symmetries, such as 2D p+ ip superconductor), that two topologically distinct systems (at least in two dimen- sional space) must have zero energy excitations at their boundary. This observation was made by examining the Green’s function G(Ω,k) only, without any reference to
ar X
iv :1
01 1.
22 73
v1 [
co nd
-m at.
me s-h
all ]
10 N
ov 20
10

2 whether the system is interacting or noninteracting, al- though as we will see in this paper, an implicit assump- tion of the absence of interactions was indeed made.
The goal of this paper is to construct the framework which could in principle generalize these observations to arbitrary number of dimensions for all 10 symmetry classes. At the basis of this framework lies the construc- tion of the topological invariants in terms of the Green’s functions in arbitrary number of dimensions. One should note that simply rewriting the topological invariants in terms of G−1(0,k) instead ofH(k) is not sufficient. To be able to relate the topological invariant to the zero energy excitations one needs to have it written in space-time, as topological characteristics of the map from Ω,k space to G.
Furthermore, in this work we are able to eliminate the implicit assumption of the absence of interactions. We do it by observing that two topologically distinct systems, in the presence of interactions, can have zeroes of the Green’s functions at their boundary instead of the poles. The zero of the Green’s function is defined as the zero of one of its eigenvalues, or equivalently of its determi- nant. If this indeed happens, two interacting topological insulators can have no edge states at the boundary be- tween them and can be continuously deformed from one to another without ever closing the bulk gap. Of course, this picture does not explain when this scenario happens, rather it just says that it in principle can happen.
In some cases one can give arguments that zero en- ergy excitations cannot simply vanish as interaction is increased9. These arguments, when they are at work, lead us to conclude that the poles cannot simply get re- placed by the zeroes. On the other hand, it is also possi- ble that while the edge states disappear as single particle excitations, they remain as collective excitations and con- tinue to carry current. In this case, perhaps it is possible to reconcile the picture from Ref. [9] and the one from this paper.
At the same time, recent studies of interacting topolog- ical insulators in one dimensional space10–12 show that in some cases the edge states do disappear as interactions are increased. To make contact with this work, in this pa- per we give the details of the Green’s function construc- tion in one dimensional space, for systems in symmetry class AIII or BDI5.
The rest of the paper is organized as follows. In sec. II we explain how to generalize the symmetry classes of topological insulators to the case when interactions are present. In Sec. III we discuss the structure of poles and zeroes of the Green’s functions. In Sec. IV we ex- plain why the zeroes of the Green’s function can change the topological invariants without any zero energy ex- citations. In Sec. V we discuss the application of this formalism to two dimensions for systems of the integer quantum Hall type, which mostly follows the derivation from Ref. [8], but takes into account the possible zeroes of the Green’s functions. This construction can be gen- eralized to any symmetry class without chiral symmetry,
which are classes A, AI, AII, D and C, in even number of spacial dimensions5. Finally, Sec. VI we study chiral topological insulators in remaining classes AIII, BDI, CII, DIII, CI, in odd number of spacial dimensions. We espe- cially concentrate on one dimensional AIII or BDI topo- logical insulators. At the end of the paper in Sec. VII, we discuss possible directions of further development of this formalism.
II. GREEN’S FUNCTIONS OF INTERACTING TOPOLOGICAL INSULATORS
We begin by reviewing the symmetries of the noninter- acting topological insulators. Those are characterized by matrix Hamiltonians Hij . Here, and in the rest of this paper, the indices i, j denote sites of a lattice, spin, and other degrees of freedom present in our problem. The total number of values i takes is equal to the number of degrees of freedom in the problem, we will denote this Df . Thus we take for the Hamiltonian of the problem
Hˆ = aˆ†iHij aˆj , (4)
where aˆ†i , aˆi are fermionic creation and annihilation op- erators. In this formula, as well as in many others below, summation over the indices i, j is implied (in some, but not all, cases we will write the summation symbol ex- plicitly). We also note that sometimes in what follows i will simply mean imaginary unit rather than an index; we hope the distinction will be clear from context.
These may satisfy one of three constraints. The first is time reversal invariance
U†TH∗UT = H. (5)
The second is charge conjugation symmetry
U†CH∗UC = −H. (6)
Here UC and UT are unitary matrices, satisfying
U∗CUC = ±1, U∗TUT = ±1. (7)
In other words, each of these two symmetries, if it is present, can be implemented in two distinct ways, with either +1 in Eq. (7) or with −1.
In the presence of these two constraints, the Hamilto- nian also automatically satisfies the constrain
Σ†HΣ = −H, (8)
where Σ = U∗CUT is a unitary matrix. This is called chiral, or sublattice, symmetry. The sublattice symmetry can also be present in case when neither time reversal invariance or charge conjugation symmetry are present (although it cannot be present if only one of these two is missing). All together, taking into account Eq. (7), one can see that there are exactly 10 separate symmetry

3 classes of Hamiltonians, distinguished by the presence of absence of these constraints5,6,13,14.
Note that the Green’s functions, defined by
G(Ω) = [Ω−H]−1 , (9) satisfy the following constraints in case when their parent Hamiltonians satisfy either of Eqs. (5), (6), and (8). Time reversal invariance leads to
G(Ω) = U†TG T (Ω)UT , (10)
charge conjugation leads to
G(Ω) = −U†CGT (−Ω)UC , (11) (here the symbol GT denotes the transposed Green’s function) and sublattice symmetry leads to
G(Ω) = −Σ†G(−Ω)Σ. (12) We would like to generalize the symmetries given by
Eqs. (5), (6), and (8) to the case when interactions are present, in such a way that the interacting Green’s func- tion G(Ω) still satisfies either of the relations Eqs. (10), (11), and (12). In the rest of this section, we show that when properly generalized to the presence of interactions, the symmetries of the Hamiltonian lead to exactly the same constraints on the Green’s functions as in the ab- sence of interactions15.
We first need to clarify what we mean by the interact- ing Hamiltonian and the interacting Green’s function. The interacting Hamiltonian is constructed in terms of the creation and annihilation operators in the same way as the noninteracting one, except it is no longer quadratic in these operators. Interacting Green’s function is then constructed in the following way.
We introduce the retarded Green’s function16 Gij(Ω). It can be analytically continued onto the upper half plane of complex Ω, where, on the purely imaginary axis, it becomes the imaginary time ordered Green’s function Gij(iω), where Ω = iω and ω will be taken to be real throughout this text. Since the system is supposed to be insulating and has a gap in the spectrum, Gij(Ω) can also be continued onto the lower half plane through the gap where, just below the real axis, it becomes the advanced Green’s function. For the real values of Ω where the density of states of the problem is not zero, the Green’s function has cuts on the real axis through which it can- not be analytically continued (if it is, it will spill into another Riemannian sheet). In some cases we will be in- terested in topological insulators with boundary modes which will manifest themselves in excitations at zero and close to zero. In this case, to define the Green’s function analytic everywhere we can think of a finite system so that the singularities of G(Ω) are confined to a number of poles. We will elaborate on the analytic properties of the Green’s function in the next section.
In addition to the Green’s functions defined here, it is often customary to introduce time-ordered, anti-time
ordered and Keldysh Green’s functions17. To distinguish them from the Green’s functions defined in the previ- ous paragraph, we will always specify which function we deal with whenever we work with one of those functions, retaining the name “Green’s function” for the analytic retarded-advanced-imaginary-time-ordered functions.
Now let us generalize the symmetries of the Hamilto- nian to the case when interactions are present. We define unitary operators UˆT , UˆC and Σˆ, which act on creation
and annihilation operators aˆ†i and aˆi (where i labels the sites of a lattice as well as other degrees of freedom if present) in the following way
Uˆ†T aˆiUˆT = ∑ j
UT ij aˆj , U † T aˆ † i UˆT =
∑ j
aˆ†jUT † ji, (13)
Uˆ†C aˆiUˆC = ∑ j
U†Cij aˆ † j , Uˆ
† C aˆ † i UˆC =
∑ j
aˆjUCji, (14)
and
Σˆ†aˆiΣˆ = ∑ j
aˆ†jΣ † ji, Σˆ
†aˆ†i Σˆ = ∑ j
Σˆij aˆj . (15)
We now declare that a many body Hamiltonian Hˆ is time reversal invariant if
Uˆ†T Hˆ ∗UˆT = Hˆ, (16)
is particle-hole symmetric if
Uˆ†CHˆUˆC = Hˆ, (17)
and is chiral if
Σˆ†Hˆ∗Σˆ = Hˆ. (18)
Here Hˆ∗ is a complex conjugate many body Hamiltonian. The creation and annihilation operators are not affected by the operation of complex conjugation.
Now it is straightforward to check that if the Hamilto- nian is noninteracting, that is if
Hˆ = ∑ ij
aˆ†iHij aˆj , (19)
then Eqs. (16) and (13) are equivalent to Eq. (5), Eqs. (17) and (14) are equivalent to Eq. (6), while Eqs. (18) and (15) are equivalent to Eq. (8). For ex- ample,
Uˆ†T Hˆ ∗UˆT = Uˆ
† T aˆ † iH∗ij aˆjUˆT = aˆ†kU†TkiH∗ijUTjlaˆl =
aˆ†iHij aˆj = Hˆ. (20) Similarly,
Uˆ†CHˆUˆC = Uˆ † C aˆ † iHij aˆjUˆC = aˆkUCkiHijU†Cjkaˆ†k =
−aˆjH∗jiaˆk = aˆ†iHij aˆj (21)

4 (this assumes, without loss of generality, that trH = 0), and finally
Σˆ†Hˆ∗Σˆ = Σˆ†aˆ†iHTij aˆjΣˆ = ΣikaˆkHTij aˆ†lΣ†lj = −aˆ†lΣ†ljHjiΣikaˆk = aˆ†iHij aˆj = Hˆ. (22) The generalizations given here are not the only way to
generalize these symmetries to the case when interactions are present. One alternative way to generalize the chiral symmetry is given in appendix A. However, they seem to be most natural from the point of view of realistic applications.
Next, one can check that even in the presence of in- teractions Eqs. (10), (11) and (12) still hold. This is a straightforward yet tedious calculation which can be skipped on the first reading.
We work with the time ordered and anti-time ordered Green’s functions
Gtij(t) = −i 〈0|T aˆi(t) a†j(0) |0〉 , Gt˜ij(t) = −i 〈0| T˜ aˆi(t) a†j(0) |0〉 . (23)
Here the symbols T and T˜ denote time and anti-time ordering of subsequent time dependent operators. We
denote the corresponding Green’s functions Gt and Gt˜
(to avoid confusion with the notation GT which implies transposed Green’s function GTij = Gji). Given these time ordered Green’s functions, we can construct re- garded and advanced Green’s functions according to
GRij(t) = θ(t) ( Gtij(t)−Gt˜ij(t)
) ,
GAij(t) = θ(−t) ( Gtij(t)−Gt˜ij(t)
) . (24)
Fourier transform of those gives us either GR(Ω) with complex Ω in the upper half plane or GA(Ω) with com- plex Ω in the lower half plane, which are known to be analytic continuation of each other (at least in the case of the gapful spectrum), and define the same analytic func- tionG(Ω) in the entire complex plane. To prove that thus defined G(Ω) satisfies Eqs. (10), (11) or (12), depending on the symmetries present, it is sufficient to prove that the time ordered and anti-time ordered Green’s function satisfy one of the three relations,
Gt(t) = U†TG tT (t)UT , (25)
Gt(t) = −U†CGt T
(−t)UC , (26) and
Gt(t) = −Σ†Gt(−t)Σ (27) (same for the anti-time ordered Green’s function).
Start with the time reversal invariance. The presence of the symmetry Eq. (16) implies that for every many-
body state |n〉, there exists a state Uˆ†T |n∗〉 which has
exactly the same energy as the state |n〉 (here |n∗〉 de- notes the complex conjugate state). Indeed, assuming
that Hˆ |n〉 = En |n〉 gives
HˆUˆ†T |n∗〉 = Uˆ†T UˆT HˆUˆ†T |n∗〉 = Uˆ†T Hˆ∗ |n∗〉 = EnUˆ†T |n∗〉 . (28)
Without loss of generality we assume that the ground state is time reversal invariant, that is
Uˆ†T |0∗〉 = |0〉 , |0∗〉 = UˆT |0〉 . (29) If it were not time reversal invariant, we could still take
expressions such as |0〉± Uˆ†T |0∗〉 as a ground state, which are invariant.
Now assume that t > 0. Then the following sequence of transformations can be carried out
iGtij(t) = 〈0| eiHˆtaˆie−iHˆtaˆ†j |0〉 = 〈0| Uˆ†T eiHˆ
∗tUˆT aˆiU † T e −iHˆ∗tUˆT aˆ
† jUˆ † T UˆT |0〉 =
〈0∗| eiHˆ∗tUˆT aˆiU†T e−iHˆ ∗tUˆT aˆ
† jUˆ † T |0∗〉 =
U†Tik 〈0∗| eiHˆ ∗taˆke
−iHˆ∗taˆ†l |0∗〉Ulj = U†Tik 〈0| aˆle−iHˆtaˆ†keiHˆt |0〉Ulj = U†Tik iG
tT kl(t)UTlj . (30)
To go from the third line to the fourth line of Eq. (30) we took advantage of Eq. (29). To go from the third to the fourth line, we used that for general states |n〉 and |m〉 and for generic operator Aˆ, the following is obviously true
〈n| Aˆ |m〉 = 〈m∗| AˆT |n∗〉 . (31) We could also prove a similar relation for t < 0, but this is not really necessary since we can just use the retarded Green’s functions to construct G(Ω). At the same time, we can easily prove the same for the anti-time ordered
function Gt˜. Thus we have proved Eq. (25) and, by ex- tension, Eq. (10).
We now move on to charge conjugation symmetry. Just as before, each excited state is doubly degenerate, and the ground state satisfies
UˆC |0〉 = |0〉 . (32) We now do the following sequence of transformations for t > 0
iGtij(t) = 〈0| eiHˆtaˆie−iHˆtaˆ†j |0〉 = 〈0| Uˆ†CeiHˆtUˆC aˆiUˆ†Ce−iHˆtUˆC aˆ†jUˆCUˆ†C |0〉 = UCjl 〈0| eiHˆtaˆ†ke−iHˆtaˆl |0〉U†Cki = −UCjl 〈0|T aˆleiHˆtaˆ†ke−iHˆt |0〉U†Cki = −UCjl iGtlk(−t)U†Cki (33)
From this Eq. (26) easily follows. In addition, same re- lations hold for negative time t < 0 as well as for the
anti-time ordered Green’s functions Gt˜. Thus we have proved Eq. (11).

5 Finally, for the chiral (sublattice) symmetry, we repeat the same process. We take the ground state to be invari- ant under the chiral transformation, which leads to
Σˆ† |0∗〉 = |0〉 , |0∗〉 = Σˆ |0〉 . (34)
We find, for t > 0,
iGtij(t) = 〈0| eiHˆtaˆie−iHˆtaˆ†j |0〉 = 〈0|Σ†eiHˆ∗tΣˆaˆiΣˆ†e−iHˆ∗tΣˆaˆ†jΣˆ†Σˆ |0〉 = Σ†ik 〈0∗| eiHˆ
∗taˆ†ke −iHˆ∗taˆl |0∗〉Σlj =
Σ†ik 〈0| aˆ†l e−iHˆtaˆkeiHˆt |0〉Σlj = −Σ†ik 〈0|Te−iHˆtaˆkeiHˆtaˆ†l |0〉Σlj = −Σ†ik iGtkl(−t)Σlj . (35)
Eq. (27) follows immediately, and, together with the ex-
tension of this relation for t < 0 and for Gt˜, Eq. (12) follows as well.
Thus we have proven that Eqs. (10), (11), and (12) hold even in the presence of interactions.
Finally, note that regardless whether or not the system is interacting, its Green’s function satisfies
Gij(Ω) = G † ij(Ω
∗). (36)
Here G† is the Hermitian conjugate of the Green’s func- tion understood as a matrix.
III. ANALYTIC PROPERTIES OF THE GREEN’S FUNCTIONS
In this section, we discuss the analytic properties of the Green’s functions, important for subsequent applica- tions. The main result of this section is that the deter- minant of the Green’s function can always be written in the following form
detG =
∏Dh−Df s=1 (Ω− rs)∏Dh n=1 (Ω− �n)
. (37)
Here rs are real numbers representing positions of the zeroes of the determinant, �n are the excitation energies which are also obviously real and which represent the poles of the determinant, Dh is the combined dimension of the Hilbert spaces of the system with one extra par- ticle than the ground state and one less particle than the ground state, and Df is the number of degrees of freedom, or number of distinct creation or annihilation operators as discussed earlier. Generally
Df ≤ Dh, (38)
where the equality is achieved, in particular, if the system is not interacting. This can be seen if one notes that in the presence of interactions Dh consists of all the states obtained from the ground state by adding or removing
a particle (which can be done in Nf distinct ways) plus states which differ from those by particle-hole excitations.
If the system is not interacting, Dh = Df and the determinant of its Green’s function is particularly easy to calculate starting from G = [Ω−H]−1, to give
detG =
Df∏ n=1
1
Ω− �n , (39)
which is consistent with Eq. (37). In the presence of interactions, the determinant of the
Green’s functions, in addition to its poles, has zeroes as follows from Eq. (37). Unlike the poles which have clear physical meaning as excitation energies, the zeroes are more subtle and, as we will see below, their positions re- flect the matrix elements of the creation and annihilation operators. However, it should be clear that as interac- tions of the system are switched off, the zeroes move to coincide with the positions of some of the poles, resulting in disappearance of one zero and one pole at a time from the determinant, until in the absence of interactions it takes the form Eq. (39). Conversely, if the interactions are turned on, the zeroes and poles appear in pairs.
In the rest of this section we present the proof of Eq. (37). This can be skipped on the first reading. We begin by writing down the spectral decomposition (sometimes referred to as Lehmann decomposition) of the Green’s function
Gij(Ω) = ∑ n
〈0| aˆi |n〉 〈n| aˆ†j |0〉 Ω− ω+n
+ ∑ m
〈0| aˆ†j |m〉 〈m| aˆi |0〉 Ω− ω−m
.
(40) Here |n〉 are the states which have one more particle com- pared to the ground state, and
ω+n = En(N + 1)− E0(N) ≥ 0, (41) where E0(N) is the ground state of the system of N par- ticles (we take the chemical potential as included in the Hamiltonian, thus dE0(N)/dN = 0), and En(N + 1) is the excited state of a system with one more particle. Sim- ilarly, |m〉 have one less particle than the ground state, and
ω−m = E0(N)− Em(N − 1) ≤ 0. (42) It is convenient for the moment to work with combined notations, where all states are denoted by |n〉 regardless or not whether they have one more particle or one less particle than the ground state and all the energy is de- noted by �n which can be either positive or negative and coincide with either ω+n or ω
− m. Then we can write
Gij(Ω) = ∑ n
U†inUnj Ω− �n , (43)
where Unj = 〈n| aˆ†j |0〉 if |n〉 has one more particle com- pared to the ground state and Unj = 〈0| aˆ†j |n〉 if |n〉 has

6 one less particle. Note that∑ n
U†inUnj = δij . (44)
In the matrix Unj , n goes over exactly Dh values while j goes over Df values, where, as we recall, Df ≤ Dh.
Now we are in the position to study the determinant of the Green’s function detG, to derive Eq. (37).
If the system is non-interacting, then it has as many excited states as the number of its degrees of freedom or Df = Dh. Then Unj is a unitary matrix which can be diagonalized. It follows that the eigenvalues of the Green’s functions are
λµ = 1
Ω− �µ , (45)
so that we immediately reproduce Eq. (39). If, however, the system is interacting, then Dh can
be larger than Df and Unj is a rectangular matrix. In that case, we write down the formal expression for the determinant of the Green’s function, starting from the definition Eq. (43)
detG = ∑ σ
∑ n1,n2,...
∏ i
U†iniUniσ(i) Ω− �ni
(−1)P (σ) (46)
Here as usual σ(i) is a permutation of the numbers 1, 2, . . . , Nf , and P denotes its parity.
It is now easy to see that if any two n coincide, for example if n1 = n2, then the summation over permuta- tions gives zero. This means that the determinant of the Green’s function has simple poles at each point �n (and does not have any higher order poles). This justifies the denominator of Eq. (37).
To explain the origin of the numerator, we first note that the determinant of the Green’s function is clearly a rational function of Ω, thus the numerator must be a polynomial. Moreover, it is well known16 that at very large Ω the Green’s function asymptotically approaches a very simple form Gij ∼ δij/Ω. This means that its de- terminant must go as 1/ΩDf . This implies that the poly- nomial in the numerator should be of the order Dh−Df , with the highest order term ΩDh−Df having the coeffi- cient equal to unity.
Thus the numerator of the determinant must take pre- cisely the form shown in Eq. (37). The only thing which we have not yet proven is that all rs are real.
To prove that, let us prove that the eigenvalues of the Green’s function cannot have zeroes (cannot vanish) when Ω is in the upper half plane or in the lower half plane. In this proof we follow the arguments given in Refs. [16,18]. Let us examine the eigenvalue equation
∑ jn
U†inUnj Ω− �n ψ
(µ) j (Ω) = λµ(Ω)ψ
(µ) i (Ω). (47)
Then for Ω = E where E is a real number which makes the Green’s function hermitian
λµ(E) = ∑ n
∣∣∣∑j Unjψ(µ)j (E)∣∣∣2 E − �n =
∑ n
Z (µ) n
E − �n , (48)
where the notation
Z(µ)n (E) =
∣∣∣∣∣∣ ∑ j
Unjψ (µ) j (E)
∣∣∣∣∣∣ 2
≥ 0 (49)
was introduced. In addition, if Ω is taken to infinity, λµ(Ω) ∼ 1/Ω as we have already discussed.
Now let us take Ω around the contour shown on Fig. 1. The horizontal part of the contour goes infinitesimally above the real axis. On this contour, it follows from Eqs. (48) and (49) that the imaginary part of an eigen- value λµ(Ω) is always zero or negative (negative if Ω is close to one of �n). At the same time, when Ω goes over a large semicircle, λµ(Ω) = 1/Ω and its imaginary part is also negative. Thus Imλµ(Ω) ≤ 0 for all Ω belonging to the contour shown on Fig. 1.
However, there exist a well known theorem in complex analysis which states that if an analytic function λ(Ω) has poles or zeroes in a certain domain, then if one plots λ(Ω) as Ω encircles this domain counterclockwise, the number of times λ(Ω) encircles zero in the counterclock- wise direction is equal to the number of zeroes minus the number of poles of this analytic function inside this do- main. We see that as Ω follows the contour shown on Fig. 1, λµ(Ω) does not encircle zero at all. Since λµ(Ω) cannot have any poles in the upper half plane, it follows that it cannot have any zeroes in the upper half plane.
Similarly, λµ(Ω) cannot have any zeroes also in the lower half plane. This means they can have zeroes only on the real axis.
If the eigenvalues of a matrix have zeroes on the real axis only, this implies that the determinant of the matrix can have zeroes on the real axis only as well. Then the fact that all rs in Eq. (37) are real follows.
This concludes our proof of Eq. (37).
IV. MODE COUNTING VIA THE GREEN’S FUNCTIONS
At the heart of the approach to the topological insula- tors developed by Volovik8 lies the following topological invariant
N1 = Tr
∫ ∞ −∞
dω
pii K∂ωG. (50)
Here, first of all, ω is a real variable such that Ω = iω. Second, Kij is a matrix inverse to Gij∑
j
KijGjl = δil. (51)

7 Ω
poles
FIG. 1: The contour Ω is taken along in order to prove that the eigenvalues of the Green’s functions cannot vanish in the upper half plane.
The notation G−1, which would seem to be natural to use instead of K in this instance, is going to be reserved for a different quantity which we will define later.
The existence of this topological invariant reflects the fact that pi1(GL(Df ,C)) = Z. It is straightforward to see that if the system is perturbed and its Green’s function changes by small amount as G → G + δG, this quantity remains unchanged.
To see the meaning of this quantity, we express it in terms of the eigenvalues of the Green’s function
N1 = ∑ µ
∫ ∞ −∞
dω
pii
∂ ln (λµ)
∂ω =
∫ ∞ −∞
dω
pii
∂ ln detG
∂ω . (52)
It is easy to see by direct substitution of Eq. (37) that the invariant is equal to
N1 = ∑ n
sign �n − ∑ s
sign rs. (53)
that is the difference between the total number of pos- itive and negative poles, minus the difference between the total number of positive and negative zeroes of the determinant of the Green’s function.
For noninteracting systems, N1 is simply the difference between the number of positive and negative eigenvalues of the system8, as noninteracting systems have no zeroes. For interacting systems, N1 turns out to be more compli- cated. Notice, however, that N1 can change only if a pole crosses zero (it increases by one if the pole moves from negative to positive values), or if a zero crosses zero (it increases by one if a zero moves from positive to negative values).
This immediately leads to the following conclusion. Suppose we have an insulating system (it has a gap in
the spectrum). If the system is noninteracting, the only way it can increase the number of positive energy levels is by moving an energy level from negative to positive values (and thus closing the gap somewhere in the mid- dle of the process where this energy level crosses zero). If the system is interacting, it can change the number of positive energy levels by moving a zero from positive to negative values (where it eventually may “fuse” with a pole and remove it from the spectrum). This process happens without ever closing the gap. A possible sce- nario outlining the motion of zeroes and poles is shown on Fig. 2.
(d)
λ
E
λ
E
λ
E
λ
E (a) (b) (c)
FIG. 2: A possible scenario outlining the change in the po- sition of the poles (solid lines) and zeroes (dashed lines) as a function of some parameter λ, as the interaction strength is varied. (a) an energy level crossing zero as a function of λ for a noninteracing system. (b) Two zeroes and and two nearby poles appear in the determinant of the Green’s func- tion Eq. (37) as the interactions are increased (c) zeroes move closer to the positions of the original poles of the noninter- acting problem (d) zeroes “annihilate” the poles resulting in a system with exactly the same change of the topological in- variant, but without any energy levels crossing zero energy.
This observation lies at the core of the possible transi- tions between two topologically distinct insulating states without closing the gap, in the presence of interactions. We will further elaborate on it later in this paper.
V. TOPOLOGICAL INSULATORS IN TWO DIMENSIONAL SPACE
In two spacial dimensions there are three types of topo- logical insulators characterized by the topological invari- ant Z. These are class A (an example is integer quantum Hall effect), class D (an example is a p + ip supercon- ductor) and class C (an example is a 2D superconductor which is not time reversal invariant but with spin rotation invariance)5,13. More generally, representative systems of a total of 5 symmetry classes can possess a Z invariant in even spacial dimensions.
The topological invariant which characterizes these systems in the context of integer quantum Hall effect in the absence of interactions is referred to as TKNN invariant19. It is usually written as a flux of the Berry connection through the Brillouin zone. However it is also possible to write it in terms of the Green’s functions, as was first observed in Ref. [20]. The advantage of the ex- pression in terms of Green’s functions is that it is then straightforward to generalize it to the case when interac-

8 tions are present. To write down this invariant, we first introduce the
Green’s function in the momentum space, Gab(Ω,k) (this assumes that our system is translationally invariant). Here the indices a, b label the remaining degrees of free- dom once momentum is introduced (these can be the ba- sis of the lattice, spin, and other degrees of freedom). We also introduce the notations k0 = ω, where Ω = iω and ω is real, in addition to the momentum vector components k1 and k2.
Now the TKNN invariant takes the form20
N3 = (54)
�αβγ 6
tr
∫ ∞ −∞
dω
∫ d2k
(2pi)2 G−1∂kαGG
−1∂kβGG −1∂kγG,
where summation over α, β, γ taking values 0, 1, and 2 is implied, and tr denotes the trace over the matrix indices of G. G−1 is a matrix inverse of G, defined as in∑
b
G−1ab (Ω,k)Gbc(Ω,k) = δac. (55)
Following Ref. [8], we denote this invariant by N3, to emphasize its similarity with Eq. (50). The value of this invariant is an integer number which, in case of integer quantum Hall effect, is proportional to the Hall conduc- tivity σxy.
However, we can now imagine interactions to our sys- tem. The presence of interactions does not change the fact that N3 is a topological invariant. Indeed, its ex- istence reflects the homotopy group pi3(GL(M,C)) = Z for a matrix of a size M by M and is totally unrelated to whether the Green’s function describes interacting or noninteracting system.
However, a question still remains whether the existence of this topological invariant in a system with interactions would have any consequences for its physical properties. In Ref. [8] it was shown that in the absence of interac- tions, if there are two domains separated by a straight one dimensional domain wall where N3 takes on different values on either side of the domain wall, this implies that there are zero energy excitations on the domain wall, that is, the edge states.
Let us see what remains of this phenomenon in the presence of interactions. We will see that if two domains with two different values of the invariant share a domain wall which separates them, then even in the presence of interactions one can have zero energy excitations on the domain wall. However, the new effect which appears in the presence of interactions is that instead of a zero en- ergy excitation, the Green’s function can also develop a zero at zero energy on the domain wall. This would still lead to the change in topological invariant from the first to the second domain. Thus in the presence of in- teractions, sometimes it is possible to have to adjacent distinct topological insulators without any edge states in between.
To derive this result, we closely follow the discussion in Ref. [8]. Let us imagine that we have two distinct topological insulators in the domains x > 0 and x < 0 of the two dimensional plane (assuming the two dimen- sional cartesian coordinates x and y). x = 0 is the po- sition of the domain wall. Suppose the parameters of the Hamiltonian smoothly change as a function of x, so that asymptotically at large positive and large negative x the Hamiltonian looks like that of a topological insulator with the invariant NR3 and N
L 3 .
The Green’s function of this topological insulator is no longer just the function of the momenta, as the mo- mentum kx is no longer a good quantum number. How- ever, the translational invariance in the y-direction is still present. Generally, the Green’s function, a matrix in the position space, can be written as Gab(Ω;x, x
′; y−y′). One can perform a Fourier transform of the Green’s function in the y-direction, along the domain wall, and a Wigner transform17,21 of the Green’s function in the x-direction perpendicular to the domain wall, to arrive at the Green’s function Gab(Ω,k, x).
With the help of the Wigner transformed Green’s func- tion, we define a local x-dependent topological invariant (k0 = ω, k1 = kx, and k2 = ky as before and α, β, γ are summed over the values 0, 1, and 2)
N3(x) = (56)
�αβγ 6
tr
∫ ∞ −∞
dω
∫ d2k
(2pi)2 G−1∂kαGG
−1∂kβGG −1∂kγG.
Here G−1(Ω,k, x) is defined as satisfying the equation∑ b
G−1ab (Ω,k, x)Gbc(Ω,k, x) = δac. (57)
Here, and below, we use the notation tr to denote the trace over matrix indices which do not include coordi- nates or momenta, and the notation Tr whenever the trace includes the summation over coordinates and mo- menta, such as in Eq. (58).
Crucially, G−1 is no longer the inverse of the Green’s function and is distinct from the Wigner transform of the matrix K defined earlier in Eq. (51). K at the same time is the true inverse of the Green’s function. Far away from the domain wall where the Hamiltonian and the Green’s function stop depending on x explicitly, G−1 and the Wigner transform of K coincide.
Since the topological invariant can take integer values only, as x varies we expect that N3 charges abruptly from NL3 to N
R 3 at some value of x, say x = 0.
Using a technique developed in Ref. [8], we can show that the difference NR3 − NL3 is related the topological invariant N1 defined in Eq. (50). To do that, consider N1 calculated at a certain value of momentum ky = Λ > 0
N1(ky) = Tr
∫ dω
pii K(Ω, ky) ∂ωG(Ω, ky), (58)
where the symbol Tr denotes summation over all the in- dices not represented by py and the second line of this

9 equation is added for later convenience. Construct the difference
1
2
( N1|ky=Λ − N1|ky=−Λ
) , (59)
where Λ is some very large positive momentum. For the noninteracting system, this gives the excess of positive energy excitations at ky = Λ over those at ky = −Λ, as follows from Eq. (53). If this number is nonzero, this means that there is a branch of excitations which crossed zero as a function of ky. In other words, there are zero energy excitations (edge states) at the boundary between two topological insulators.
For interacting system, this is a more subtle quantity as we have seen in the previous section. It can be nonzero even if there are no edge states, by moving a zero of the Green’s function across zero energy as ky changes. The formalism we are developing cannot tell which of the two scenarios happen in general, but nevertheless it can relate this number to the difference NR3 −NL3 .
To derive this, we need a relation between K and G−1. The exact relationship between those is very complicated. If however the x-dependence of the Green’s function is weak (the Hamiltonian varies slowly as x is taken across the domain wall) the relationship between them can be derived by the so-called Moyal product expansion or a gradient expansion17. This is defined in the following way. Suppose we have two functions A(x, y) and B(x, y), and a third function C(x, y) which is a convolution of the first two,
C(x, z) =
∫ dy A(x, y)B(y, z). (60)
Then the Wigner transform of these functions is approx- imately related to the Wigner transform of the first two by the following expansion
C(k, x) = A(k, x)B(k, x) + (61)
1
2i (∂xA(k, x) ∂kB(k, x)− ∂kA(k, x) ∂xB(k, x)) + . . . ,
valid when x dependence of these functions is weak. Note that the trace of C can be calculated in a straight-
forward way∫ dxC(x, x) =
∫ dxdk
2pi C(k, x) =
∫ dxdk
2pi A(k, x)B(k, x).
(62) The inverse Green’s function K’s convolution with the
Green’s function G is a delta-function. This allows us to calculate K approximately via G−1 with the help of the gradient expansion
K = G−1 + 1
2i
( G−1∂xGG−1GkxGG
−1− G−1∂kxGG
−1∂xGG−1 )
+ . . . . (63)
We use this to calculate Eq. (58) using the gradient expansion. We substitute Eq. (63) into Eq. (58) (there is
no need to use Moyal product expansion in the product of K and ∂ωG, due to Eq. (62)). We find
N1(ky) = tr
∫ dxdkxdω
4pi2 ( G−1∂xGG−1∂kxGG
−1∂ωG− G−1∂kxGG
−1∂xGG−1∂ωG )
=∑ αβγ
�αβγtr
∫ dxdkxdω
12pi2 G−1∂αGG−1∂βGG−1∂γG. (64)
Here α, β, γ go over x, kx, ω. The next step relies on the construction of a vector
nα = �αβγδtrG −1∂βGG−1∂γGG−1∂δG, (65)
where α, β, γ, δ now span the four-dimensional space of x, kx, ky, and ω. This vector satisfies
∂αnα = 0, (66)
as can be checked by direct differentiation. The relevant difference of the invariant Eq. (59) can be thought of as the flux of this vector through a three dimensional surface formed by the planes py = Λ and py = −Λ, surrounding the singular point ω = 0, x = 0, while kx and ky are tuned to appropriate values.
That same flux can be written as a flux through the surfaces x = L and x = −L. Finally, this can be brought to the form
1
2 (N1(Λ)−N1(−Λ)) = N3(L)−N3(−L). (67)
This formula relates what happens at the boundary to the difference of topological invariants. In case of non- interacting system, it proves that for two topologically distinct systems bordering each other, there are always boundary zero energy excitations (chiral edge states). In the presence of interactions, it is possible for the left hand side of this formula to be nonzero without having any zero energy excitations at the boundary as we saw be- fore. Whether this happens in practice requires further study of the zeroes of the Green’s functions in interacting topological insulators.
VI. TOPOLOGICAL INSULATORS WITH CHIRAL SYMMETRY
We now proceed to generalize the previous discussion to the topological insulators with chiral or sublattice symmetry. These include the insulators in classes AIII (sublattice symmetry only), BDI (sublattice symmetry and time reversal invariance), CII (sublattice symmetry and spin orbit coupling), DIII (time reversal invariant superconductors with spin orbit coupling such as phase B of He3) and CI (spin singlet time reversal invariant superconductor), all of which can possess a Z topologi- cal invariant in odd spacial dimensions only, unlike the remaining insulators including those discussed in the pre- vious chapters which could possess a Z topological invari- ants only in even spacial dimensions. Expressions for the

10
topological invariants such as Eq. (54) work only in the absence of chiral symmetry (they are designed to work in even number of spacial dimensions). In what follows we propose equivalent expressions for chiral topological insu- lators and proceed to relate those to zero energy bound- ary excitations and zeroes of the Green’s function.
The simplest of these insulators are the ones in class AIII, since they do not have any symmetries besides the sublattice symmetry. Let us first review the noninteract- ing topological insulators of this type.
A. Noninteracting chiral topological insulators
Consider a noninteracting topological insulator in class AIII. A noninteracting fermionic problem is described by its Hamiltonian
Hˆ = ∑ ij
Hij aˆ†i aˆj , (68)
where the indices i, j refer to points on a lattice, spin, and/or species of fermions in the problem. Such a system belongs to class AIII if there exist a unitary matrix Σij ,∑
j
Σ†ijΣjk = δik, (69)
such that ∑ jk
Σ†ijHjkΣkl = −Hil. (70)
This Hamiltonian is said to be insulating if it has a single particle gap.
It is straightforward to see that ∑ j ΣijΣjk = δik as
a consequence of Eq. (70), thus the eigenvalues of Σ are either +1 or −1.
The existence of this “chiral” or “sublattice” symme- try given by Eq. (70) immediately leads to a number of consequences. For any single particle state ψi of energy �n, satisfying
Hijψj = �nψi, (71) there exists a state
∑ j Σijψj with energy −�n. Among
the states there may be zero energy single particle states satisfying
Hijψ(0)j = 0. (72) Those are also eigenstates of Σij with the eigenvalues ei- ther +1 or −1. The zero energy eigenstates whose eigen- value of Σij is +1 are called the right zero modes, and those with the eigenvalue −1 the left zero modes,
Σψ (0) L = −ψ(0)L , Σψ(0)R = ψ(0)R . (73)
The difference of the number of the right and left zero modes is a topological invariant (this difference cannot
change if Hij is deformed smoothly). In fact, that differ- ence is called the index of the operator Vαβ from Eq. (75) and the fact that it is a topological invariant is referred to as the “index theorem”. The existence of the index theorem for the chiral topological insulators, as well as its absence in non-chiral topological insulators, in particular in p+ ip superconductors, was emphasized in Ref. 22.
It is convenient to work in the basis of the eigenstates of Σ. In this basis, Σ looks like
Σ =
( δαβ 0 0 −δαβ
) . (74)
We often refer to the sites of the lattice which are +1 eigenstates of Σ as even and −1 eigenstates as odd sub- lattice. The Hamiltonian takes form
H = (
0 Vαβ V †αβ 0
) . (75)
Here α and β label sites of each of the sublattices. It is straightforward to check that Eq. (75) is consistent with, and follows from, Eq. (70).
For a translationally invariant system, we can intro- duce momentum (crystal momentum if on the lattice). In this basis, the Hamiltonian takes the form
H = (
0 Vab(k)
V †ab(k) 0
) . (76)
Here a and b label remaining states within the lattice basis after the momentum is introduced.
Such a Hamiltonian is known to be characterized by a topological invariant in any odd dimensional space. The invariant is constructed the following way. Vαβ is a ma- trix function of momentum k. We can construct an ex- pression
Id ∼ �α1α2...αdtr ∫ ddk
d∏ n=1
V −1∂αnV. (77)
Here ∂α = ∂kα and appropriate summation is implied. The integral goes over all the entire d-dimensional space, or, in case if the problem is defined on a lattice, over the Brillouin zone. This expression takes on values propor- tional to an integer, which, mathematically speaking, re- flects the fact that pid(GL(M,C)) = Z for odd d. This in- variant is the odd-dimensional counterpart of the TKNN invariant, which exists in the symmetry classes without chiral symmetry only, and is nonzero only in even d.
It is especially simple to see that this is an integer if the matrix V is a one by one matrix, or reduces to a number z(k), and if d = 1. Then
I1 =
∫ dk
4pi
( 1
z
∂z
∂k − 1 z∗ ∂z∗
∂k
) . (78)
This measures the number of times z(k) winds around the origin as k goes over the Brillouin zone.

11
One famous example of the topologically nontrivial in- sulator in one dimension with the symmetry Eq. (70) are the Su-Schrieffer-Heeger23 solitons. Indeed, these au- thors considered a problem of a particle hopping on a 1D lattice with the Hamiltonian
Hˆ = ∑ i
( t aˆ†2iaˆ2i+1 + t
′ aˆ†2i+1aˆ2i+2 )
+ h.c. (79)
It is straightforward to see that the matrix H in this case takes the form
H(k) = (
0 t+ t′eik t+ t′e−ik 0
) . (80)
We see that t+ t′eik winds around the origin if t′ > t and does not if t′ < t. It immediately follows that if for i > 0, t > t′, and for i < 0, t < t′, at the boundary between t > t and t < t′ insulators there is going to be a zero mode, the Su-Schrieffer-Heeger soliton (this can of course be established also in a direct fashion by constructing the solutions to the Schro¨dinger equation).
Technically speaking the Hamiltonian Eq. (79) pos- sesses not only chiral symmetry Eq. (70) but also time- reversal invariance, so it belongs not to the class AIII, but rather to the class BDI. However, in 1D the distinction between these two classes is not significant (although it becomes more significant in spaces of higher dimension- alities).
More generally, any problem where a particle hops on a bipartite lattice possesses the symmetry Eq. (70) and is said to belong to class AIII (in the absence of any other symmetries), with the matrix Σ acting by +1 on even and by −1 on odd sites, as long as a particle can only hope from an odd site to an even site or vice versa.
B. Interacting chiral topological insulators
We would like to see if adding interactions to the Hamiltonian (68) preserves its topological properties, es- pecially the edge states. To do that, we need to generalize Eq. (70) to interacting systems. We do it in the way dis- cussed earlier in Sec. II, in particular in Eqs. (15), (18) and (12).
As a result of the symmetry defined by Eq. (18) for any
state |n〉 there exists a state Σˆ† |n∗〉 which has exactly the same energy as the state |n〉,
Hˆ |n〉 = E |n〉 , HˆΣ† |n∗〉 = EΣˆ† |n∗〉 . (81) Notice also that if the state |n〉 has N particles, then its conjugate state Σˆ† |n∗〉 has Df−N particles. To see that, we write
Df∑ i=1
aˆ†i aˆiΣˆ † |n∗〉 = Σˆ†Σˆ
Df∑ i=1
aˆ†i aˆiΣˆ † |n∗〉 =
Σˆ† Df∑ i=1
aˆiaˆ † i |n∗〉 = Σˆ†
Df − Df∑ i=1
aˆ†i aˆi
 |n∗〉 . (82)
The lowest energy state of the system obeying this kind of symmetry has exactly N = Df/2 particles, that is, the system is at half filling. Indeed, suppose the number of particles in the lowest energy state is not Df/2. That means that the ground state energy at particle number fixed at N , E0(N), has a minimum at N 6= Df/2. At the same time, it satisfies
E0 (N) = E0 (Df −N) (83) as a consequence of Eqs. (81) and (82). Therefore, this function also has exactly the same minimum at N1 = Df −N particles. This implies that the compress- ibility of the system, d2E0/dN
2, is not positive definite. From now one we restrict our attention to systems with positive compressibility, which as we just saw is equiva- lent to systems at half filling.
An interacting system invariant under Eq. (18) pos- sesses a topological invariant, a generalization of Eq. (77) for the case when there are interactions and reducing to it in case where there are no interactions. To define it, we again introduce the Green’s function Gij(Ω). For a translationally invariant system, it is a function of mo- mentum k. Otherwise, it has the same matrix structure as the noninteracting Hamiltonian Eq. (76).
We construct the following matrix out of the Green’s functions (here, as before, ω is taken to be real)
Q(ω,k) = G−1(iω,k) ΣG(iω,k). (84)
Here matrix multiplication is implied, while G−1 is the inverse of the matrix G(Ω,k). It is clear that
Q2 = 1. (85)
We now construct the expression
ID ∼ �α1,α2,...,αD tr ∫ ∞
0
dω
∫ ddk Q
D∏ n=1
∂αnQ. (86)
Here D = d+ 1 is the dimension of space-time, ∂i = ∂ki and ∂0 implies ∂ω. This expression is a topological invariant, which mathematically reflects the fact that piD (U(2M)/U(M)× U(M)) = Z for even D. One can check that for the noninteracting systems where G = [Ω−H]−1, Eq. (86) reduces to Eq. (77).
Note that the integration over ω goes from 0 to infinity, not from −∞ to ∞. The reason for that is that the integral from −∞ to ∞ is zero, as can be verified using
Q(ω) = Σ†Q(−ω)Σ, (87) which in turn follows from the definition Eq. (84) as well as the chiral symmetry Eq. (12).
To understand why Eq. (86) is a topological invariant, we need to further elucidate the structure of the chiral Green’s function. Taking advantage of the explicit form of the matrix Σ given in Eq. (74), we arrive at the fol- lowing structure
G(Ω) =
( ∆
(1) αβ(Ω) Wαβ(Ω)
W †αβ(Ω) ∆ (2) αβ(Ω)
) . (88)

12
Here the matrices ∆, W satisfy, thanks to Eq. (12)
∆(1,2)(Ω) = −∆(1,2)(−Ω), W (−Ω) = W (Ω). (89) For completeness we give the explicit expressions for these matrices in the appendix B.
For a translationally invariant system, the Green’s function can be written in the momentum space as
G(Ω,k) =
( ∆
(1) ab (Ω,k) Wab(Ω,k)
W †ab(Ω,k) ∆ (2) ab (Ω,k)
) . (90)
In the absence of any excitations around zero energy, that is, in case when one deals with an insulator, the functions ∆, W , are regular at Ω = 0. In particular this means that ∆(Ω = 0) must vanish, while W (0) is a constant which for now we assume to be nonvanishing. This leads to
Q(ω → 0)→ (−1 0
0 1
) . (91)
At the same time, at Ω→∞, the Green’s function must behave as 1/Ω, so W vanishes in this limit and
Q(ω →∞)→ (
1 0 0 −1
) . (92)
Thus the topological invariant measures how Q interpo- lates between these two values.
Finally, one can check that Eq. (86) is indeed a topo- logical invariant by checking that if Q is varied slightly, to Q + δQ, Eq. (86) does not change. This turns out to be a direct consequence of Q2 = 1.
The topological invariant Eq. (86) is a good starting point to analyze the chiral topological insulators in 3D, including exotic singlet time-reversal invariant supercon- ductors (class CI), Helium III-B (class DIII) and a parti- cle hopping on a bipartite lattice (class AIII)24–26. How- ever, in what follows we will use it to analyze chiral prob- lems in one spacial dimension.
C. Low energy excitations at the boundaries of interacting one dimensional class AIII topological
insulators
We now have in our possession the topological invari- ant for interacting system in chiral symmetry classes. However, what we are really interested in is whether two insulators with two different values of I which are in spacial contact possess zero energy excitations at the boundary. More precisely, we would like to show that it is possible to add particles without any energy cost at the boundary, just like in our discussion of topological insulators in two spacial dimensions in Sec. V. Let us show that this is indeed the case.
From now on we specialize to one spacial dimension (although generalizations of what follows to higher di- mensions should be possible to construct). Suppose we
have a 1D interacting chiral topological insulator, whose parameters vary in space (so that it is not translationally invariant). However at far away positive spacial infinity it becomes translationally invariant, with the topologi- cal invariant IR, while at far away negative infinity its topological invariant is IL.
Just as we have already discussed in case of two spacial dimensions in Sec. V, as the system is no longer trans- lationally invariant, its Green’s function is no longer a function of momentum only. Nevertheless, we can intro- duce the Wigner transformed Green’s function, the one Fourier transformed with respect to the difference of co- ordinates. Such a Wigner transformed function depends on both the momentum k and the coordinate x. Yet at far away spacial infinity G becomes independent of x and becomes equal to its value for a translationally invariant system.
We can define Q in the same way as before,
Q(ω, k, x) = G−1(iω, k, x) ΣG(iω, k, x). (93)
Here G−1(Ω, k, x) is defined as the inverse of the ma- trix G(Ω, k, x), which is not the true matrix inverse of G. With its help, we can define a position dependent topological invariant I2 by reducing Eq. (86) to the case D = 2 (d = 1). We denote it by I as opposed to I2 as in this subsection we always work in one dimensional space so using I without specifying dimensionality should not be the cause of confusion.
I(x) = 1
16pii tr
∫ ∞ 0
dω
∫ dk Q (∂ωQ∂kQ− ∂kQ∂ωQ) .
(94) Then
IL = I|x=−L , IR = I|x=L (95) for some large L such that at x = L or x = −L we are far from the boundary of the two insulators.
Now we define a vector
nα = �αβγtrQ∂βQ∂γQ. (96)
Here α, β, and γ go over values 0, 1, 2, which refer to the three effective coordinates, ω, k and x. Notice that
∂αnα = 0. (97)
This important relation follows from the fact that
�αβγtr ∂αQ∂βQ∂γQ = 0. (98)
This in turn can be proving by using that Q2 = 1 as well as ∂QQ = −Q∂Q. Indeed,
�αβγtr ∂αQ∂βQ∂γQ =
�αβγtr ∂αQQ 2∂βQQ
2∂γQ =
−�αβγtrQ∂αQ∂βQQ2∂γQQ = −�αβγtr ∂αQ∂βQ∂γQ. (99)
We see that the expression Eq. (98) is equal to minus itself, so it is zero.

13
Similar manipulations can be used to prove that Eq. (86) is indeed a topological invariant, that is, does not change under small deformations in Q, as long as Q2 = 1.
Due to Eq. (97), the following integral is zero∫ dSαnα = 0, (100)
where the integral is taken over any closed surface in the ω, k, x space. In particular, consider a surface formed by the planes x = L, x = −L and ω = 0 (as well as ω =∞, where the integral vanishes), with dSα being the element of this surface. We note that the sum of the integrals over x = Λ and x = −Λ is equal precisely to IR − IL. This allows us to express this difference in terms of the integral over the surface ω = 0, or more precisely
IR − IL = 1 16pii
lim ω→0
tr
∫ dxdk Q (∂xQ∂kQ− ∂kQ∂xQ) .
(101) The limit is needed because at ω = 0 the expression to be integrated is singular as we will see below.
Now let us relate this difference to the number of zero energy excitations. Earlier when working in 2D we cal- culated the number of zero energy excitations by a clever application of Eq. (58). This is not possible to do in chi- ral systems however (formally the expression Eq. (58) is always zero due to Eq. (12)). Instead we take advantage of a similar, but distinct formula
N = − lim ω→0
ωTr [ΣK∂ωG] , (102)
Here K is given by Eq. (51). We recall that as before, K coincides with G−1 in translationally invariant systems only.
Let us check that Eq. (102) is indeed equal to the num- ber of zero modes. First of all, in the noninteracting case, G = [iω −H]−1. This allows us to write
N0 = lim ω→0
iωTr Σ [iω −H]−1 = NR −NL. (103)
Here NR is the number of the right zero modes and NL is the number of the left zero modes (recall that zero
modes satisfy Σψ (0) R,L = ±ψ(0)R,L). This is the correct way
of counting zero modes, since it is this difference that is insensitive to small changes in the Hamiltonian as dis- cussed after Eq. (73).
In case when there are interactions, the analysis of Eq. (102) is somewhat more involved. Define the eigen- values and eigenvectors of the Green’s function as in Eq. (47), or in other words∑
j
Gij(Ω)ψ (µ) j (Ω) = λµ(Ω)ψ
(µ) i (Ω). (104)
The trace in Eq. (102) can be understood as a sum over these eigenvectors
N = − lim ω→0
ω ∑ µ
ψ(µ) ΣK ∂ωGψ (µ), (105)
where matrix multiplication is implied. If ψ(Ω) is an eigenvector of the Green’s function G(Ω)
with the eigenvalue λ(Ω), then Σ†ψ(−Ω) is an eigenvector of the Green’s function G(Ω) with an eigenvalue −λ(−Ω), as a consequence of Eq. (12).
A special role is played by the eigenvectors ψ(0). In- deed, if it is an eigenvector of G(0), then Σ†ψ(0) is also its eigenvector with the eigenvalue −λ(0). These eigen- vectors do not contribute to Eq. (105) in the limit of ω → 0.
However, there are also the eigenvalues such that λ(Ω) ∼ Ω as Ω→ 0, or such that λ(Ω) ∼ 1/Ω as Ω→ 0. Then ψ(0) and Σ†ψ(0) are both eigenfunctions with the same eigenvalue. In that case, just like in the analysis in the previous section Eq. (73), ψ(0) are the eigenvectors not only of G but also of Σ, with the eigenvalues +1 or −1.
We call the eigenvectors with the eigenvalues λ(Ω) ∼ 1/Ω right and left poles (right poles satisfy ΣψR = ψR, left poles satisfy ΣψL = −ψL), while we call the eigenvec- tors with the eigenvalues λ(Ω) ∼ Ω right and left zeroes. Then by direct substitution we find, from Eq. (105),
N = NR −NL − ZR + ZL, (106) where ZR and ZL are the number of right and left zeroes of the Green’s function.
Now we can relate N to the difference between the values of the topological invariant by following the ap- proach described in Sec. V. This involves expressing the product of two Green’s functions in Eq. (102) in terms of their Wigner transform using gradient expansion, as well as expressing K in terms of G−1 using Eq. (63). This produces two terms. The first one, where we use the gradient expansion while substituting G−1 for K, looks like
− ω 2i
TrΣ [−G−1∂xGG−1∂ω,kG+G−1∂kGG−1∂ω,xG] .
(107) Note that this term should be zero after integration over k and x due to Eq. (62). Nevertheless we would like to keep it because it will make the final expression which will be integrated over k and x simpler.
In the second term we use Eq. (63) for K to find
− ω 2i
TrΣ [ G−1∂xGG−1GkG−1∂ωG−
G−1∂kGG−1GxG−1∂ωG ] . (108)
Combining these terms together gives
N = lim ω→0
ω tr Σ
∫ dxdk
4pii
[ G−1∂xG∂k
( G−1∂ωG
)− G−1∂kG∂x
( G−1∂ωG
)] . (109)
At the same time, Eq. (101) can be rewritten, upon sub- stituting Eq. (93), as
IR − IL = lim ω→0
tr Σ
∫ dxdk
4pii
( ∂xGG
−1∂kGG−1−

14
∂kGG −1∂xGG−1
) . (110)
The expressions Eqs. (109) and (110) are not identical. Nevertheless, they are related in the following way. The entries of the Green’s function, as introduced in Eq. (88), satisfy the relations given by Eq. (89). As a result, the Green’s function has to have one of the following possible asymptotic behaviors as ω is taken to zero. The Green’s function might have a pole as ω is taken to zero at this value of k and x leading to
G(ω, k, x) ∼ (
δ (1) ab (k,x)
iω Wab(k, x)
W †ab(k, x) δ (2) ab (k,x)
iω
) . (111)
Alternatively, it is regular
G(ω, k, x) ∼ ( iω δ
(1) ab (k, x) Wab(k, x)
W †ab(k, x) iω δ (2) ab (k, x)
) . (112)
Here δ(1) and δ(2) are the appropriate expansion coeffi- cients of ∆(1) and ∆(2) from Eq. (88).
An interesting feature of both of these expressions is that, in the limit ω → 0, Q is either Σ or −Σ, just as in Eqs. (91) and (92) and is k and x independent. As a result, naive substitution of either Eq. (111) or Eq. (112) into Eq. (110) gives zero in the limit of ω → 0. However, at the values of k and x where δ vanish in Eq. (111) or W vanish in Eq. (112), Q varies rapidly, giving rise to a delta-function-like contribution to Eq. (110). Notice that the values of k and x where W vanishes corresponds to the zero of the Green’s function. We mostly elabo- rate on this argument in two and one dimensional space, although it should be generalizable to higher dimensions.
One can check by a direct substitution that if Eq. (111) holds, then Eq. (109) is equal to Eq. (110). If, how- ever, Eq. (112) holds then Eq. (109) is equal to minus Eq. (110). This was checked for G being a two by two matrix and a four by four matrix and we conjecture (al- though this has not been checked) that this is valid for matrices of all sizes.
As a result, we find that
IR − IL = P + Z, N = P − Z, (113) where P is the contribution of the poles and Z is the contribution of the zeroes. Since the change in the topo- logical invariant must be an integer, P and Z are also integers.
Therefore, we conclude that the change in the topo- logical invariant is due to either the poles of the Wigner transformed Green’s functions or the zeroes of the Wigner transformed Green’s functions, which occur at some values of x and k, where the appropriate values of x are close to the domain wall separating the two topo- logical insulators. The poles’ contributions are exactly opposite to those of the zeroes.
In particular, the topological invariant will change even if there is no poles, P = 0, as long as there are zeroes
present, Z 6= 0. Then two topological insulators with distinct topological invariants will have no zero energy excitations at their boundary. It is tempting to conclude that this scenario is at work for the one dimensional BDI systems discussed in recent papers10–12. However, it is not currently clear why the topological invariant consid- ered there had to change by 4 (by 8 if one works with Majorana chains instead of Dirac chains considered here) and not by 1 as would have been natural here (see ap- pendix C to see that Z can be any integer).
VII. CONCLUSIONS
We have presented a scenario in this paper where an interacting topological insulator can lose its edge states as the interactions are turned on. The method developed here has not yet been developed to a point where it al- lows to understand when this actually happens. Instead, it explains the mechanism by which this can happen. One should also be careful in interpreting the results of this formalism. The edge states are understood here as poles in the single particle Green’s functions. If there is trans- port along the edge in such a way that the elementary particle carrying the current along the edge are collec- tive excitations, this formalism may interpret it as the absence of the edge states.
In this work we restricted our attention to the topo- logical insulators with the topological invariant Z. The relationship between the invariant and the edge states was explored in two and one spacial dimensions only. It is straightforward to generalize this to higher number of dimensions, in particular to three dimensional space. It is also important to generalize this formalism to the invari- ant Z2, especially in view of current interest to 2D and 3D AII topological insulators. The expression for the Z2 topological invariant in terms of the Green’s functions was found in Ref. [27]. What remains is to find what it tells us about the edge states (and zeroes of the Green’s functions). All of this will be the subject of future work.
The author is grateful to A. Altland and A. W. W. Ludwig for many stimulating discussions in the course of this work, and would like to thank S. Ryu, C. L. Kane, L. Radzihovsky, M. Hermele and E. Altman for useful comments. Some of this work was done during visits to the Institute for Theoretical Physics at the University of Cologne, as well as to the Kavli Institute for Theoretical Physics in Santa Barbara. This work was supported by the NSF grants DMR-0449521, PHY-0904017, and PHY- 0551164.
Appendix A: Alternative generalization of chiral symmetry in the presence of interactions
Instead of Eq. (18), we could use the following con- struction to generalize the chiral symmetry to the inter-

15
acting case. Consider an operator Σˆ such that
Σˆ†HˆΣˆ = −Hˆ. (A1)
This operator does not mix creation and annihilation op- erators, so that its action on those operators is simply
Σˆ†aˆiΣˆ = Σij aˆj , Σˆ†aˆ † i Σˆ = aˆ
† jΣ † ji. (A2)
We can check that noninteracting chiral Hamiltonians au- tomatically satisfy Eq. (A1) due to Eq. (70). An example of an interacting Hamiltonian in 1D invariant under this operation would be (here we have spinful fermions hop- ping on a 1D lattice)
Hˆ = ∑
i,σ=↑,↓ ti
( aˆ†σiaˆσ,i+1 + aˆ
† σ,i+1aˆσi
) + (A3)
U ∑ i
( aˆ†↑i+1aˆ
† ↓iaˆ↓iaˆ↑i + aˆ
† ↑iaˆ † ↓iaˆ↓iaˆ↑i+1
) .
Indeed, it changes sign under the sign change of every creation and annihilation operator on odd sites of the lattice.
As a result of the symmetry Eq. (A1), for every state
|n〉 there exist a state Σˆ |n〉 whose energy is exactly op- posite to that of |n〉. In particular, in addition to the ground state |0〉, the system should have a maximally excited state
∣∣0˜⟩, with the energy opposite to that of the ground state.
We can then define two Green’s functions. One, de- noted byG(Ω) is defined with respect to the ground state.
The other, denoted G˜(Ω), is defined with respect to the maximally excited state. One can now check that the sum of these two functions,
Gs(Ω) = 1
2
( G(Ω) + G˜(Ω)
) (A4)
satisfies the same relationship Eq. (12) which Green’s functions for the chiral systems are supposed to satisfy. It follows that this is a legitimate generalization of the chiral symmetry to the interacting case, distinct from Eq. (18). It is possible to prove that in the absence of interactions Eq. (A4) coincides with the usual Green’s function, but this does not have to be true in the pres- ence of interactions. We leave the question whether any interesting models of this type can be constructed and studied to further work.
Appendix B: Spectral decomposition of the chiral Green’s functions
The general spectral decomposition Eq. (40) valid for any Green’s function can be further specialized to the case where chiral symmetry is present. In the presence of chiral symmetry, any state |n〉 with one extra particle compared to the ground state acquires a partner, a state
Σˆ† |n∗〉 with one less particle than the ground state. As a result, the spectral decomposition becomes
Gij(Ω) = ∑ n
〈0| aˆi |n〉 〈n| aˆ†j |0〉 Ω− �n +
〈0| aˆ†jΣˆ† |n∗〉 〈n∗| Σˆ aˆi |0〉 Ω + �n
.
(B1) Here |n〉 are the states with one extra particle compared to the ground state, while �n = ω
+ n > 0, defined in
Eq. (41). The second term in this expression can be transformed in the following way
〈0| aˆ†jΣˆ† |n∗〉 〈n∗| Σˆ aˆi |0〉 = 〈0| Σˆ†Σˆaˆ†jΣˆ† |n∗〉 〈n∗| ΣˆaˆiΣˆ†Σˆ |0〉 = 〈n| Σˆ† aˆjΣˆΣˆ† |0∗〉 〈0∗| ΣˆΣˆ†aˆ†i Σˆ |n〉 =
Σil 〈0| aˆl |n〉 〈n| aˆ†k |0〉Σ†kj . (B2)
To go from the second to the third line, Eq. (31) was used, while in the final line takes advantage of Eqs. (15) and (34). Finally, this gives
Gij(Ω) = (B3)∑ n
[ 〈0| aˆi |n〉 〈n| aˆ†j |0〉
Ω− �n + Σil 〈0| aˆl |n〉 〈n| aˆ†k |0〉Σ†kj
Ω + �n
] .
Note that this is compatible with the constraint Eq. (12). With the help of this expression, we can write the
following explicit construction for matrices defined in Eq. (88). We introduce the matrix element
Uin = 〈0| aˆi |n〉 . (B4)
We also employ the notations U (1) αn and U
(2) αn , which are
equal to Uin as i belongs to the first or the second sub- lattice. Now Eq. (B3) gives
∆ (1) αβ =
∑ n
U (1)αnU (1)†
nβ
2Ω
Ω2 − �2n ,
∆ (2) αβ =
∑ n
U (2)αnU (2)†
nβ
2Ω
Ω2 − �2n ,
Wαβ = ∑ n
U (1)αnU (2)†
nβ
2�n Ω2 − �2n
. (B5)
Appendix C: Zeroes of the chiral Green’s functions and the topological invariant
Suppose the chiral Green’s function takes the form
G(ω) =
( iωδ(1) W W ∗ iωδ(2)
) , (C1)
where δ(1) and δ(2) are real positive functions of x and k while W is a complex function of x and k. This is a 2 by 2 version of the Green’s function Eq. (112). Let us calculate the change in the topological invariant due to

16
this form of the Green’s function. Substituting this into Eq. (110) we find
IR−IL = lim ω→0
∫ dkdx
2pii
ω2δ(1)δ(2) (∂pW∂xW ∗ − ∂xW∂pW ∗)(
δ(1)δ(2)ω2 +WW ∗ )2 .
(C2) Only the terms with nonzero ω → 0 limit are shown.
Suppose W vanishes at some special values of k and x. At this point, in the limit ω → 0, the Green’s function obviously vanishes, or as we say, has a zero.
We write W = X + iY where X and Y are real and they both vanish at those same values of k and x. Then
IR − IL = lim ω→0
∫ dkdx
pi
ω2δ(1)δ(2) (∂xX∂kY − ∂kX∂xY )( δ(1)δ(2)ω2 +X2 + Y 2
)2 . (C3)
As can be checked in a straightforward manner,
lim ω→0
ω2
(ω2 +X2 + Y 2) 2 = pi δ(X)δ(Y ). (C4)
This gives
IR − IL = ± ∫ dXdY δ(X)δ(Y ) = ±1. (C5)
One can check that the result of substitution of the form Eq. (C1) into Eq. (109) gives exactly the same expression as Eq. (C2) except with the opposite sign.
1 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005).
2 B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 106802 (2006).
3 M. Ko¨nig, S. Wiedmann, C. Bru¨ne, A. Roth, H. Buhmann, L. W. Molenkamp, X.-L. Qi, and Shou-Cheng, Science 318, 766 (2007).
4 D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nature 452, 970 (2008).
5 S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Lud- wig, New J. of Phys. 12, 065010 (2010).
6 A. Kitaev (2009), arxiv:0901.2686. 7 X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 78,
195424 (2008). 8 G. E. Volovik, The Universe in a Helium Droplet (Oxford
University Press, Oxford, 2003). 9 S. Ryu, J. E. Moore, and A. W. W. Ludwig (2010), arxiv:1010.0936.
10 L. Fidkowski and A. Kitaev, Phys. Rev. B 81, 134509 (2010).
11 L. Fidkowski and A. Kitaev (2010), arxiv:1008.4138. 12 A. M. Turner, F. Pollmann, and E. Berg (2010),
arxiv:1008.4346. 13 A. Altland and M. Zirnbauer, Phys. Rev. B 55, 1142
(1997). 14 M. Zirnbauer, J. Math. Phys. 37, 4986 (1996). 15 M. S. Foster and A. W. W. Ludwig, Phys. Rev. B 73,
155104 (2006). 16 A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshin-
skii, Methods of quantum field theory in statistical physics (Dover Publications, New York, 1975).
17 A. Kamenev, in Nanophysics : coherence and transport: Les Houches 2004 session LXXXI (Elsevier, Amsterdam, Boston, 2005), arXiv:cond-mat/0412296.
18 L. D. Landau and E. M. Lifshitz, Statistical Physics (Butterworth-Heinemann, Oxford, UK, 1980).
19 D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett. 49, 405 (1982).
20 Q. Niu, D. J. Thouless, and Y. S. Wu, Phys. Rev. B 31, 3372 (1985).
21 A Wigner transform of a function of two variables f(x, x′) is defined as f(k,R) =
∫ dr eikrf(R+ r/2, R− r/2).
22 V. Gurarie and L. Radzihovsky, Phys Rev B 75, 212509 (2007).
23 W. Su, J. Schrieffer, and A. Heeger, Phys. Rev. Lett. 22, 2099 (1980).
24 A. P. Schnyder, S. Ryu, and A. W. W. Ludwig, Phys. Rev. Lett. 102, 196804 (2009).
25 P. Hosur, S. Ryu, and A. Vishwanath, Phys. Rev. B. 81, 045120 (2010).
26 M. A. Silaev and G. E. Volovik (2010), arxiv:1005.4672. 27 Z. Wang, X.-L. Qi, and S.-C. Zhang (2010),
arxiv:1004.4229.