arxiv: v2 [cond-mat.str-el] 3 Jul 2014

Transcription

1 Boundary Fidelity and Entanglement in the symmetry rotected toological hase of the SSH model arxiv: v [cond-mat.str-el] 3 Jul 14 J. Sirker Deartment of Physics and Research Center OPTIMAS, TU Kaiserslautern, D Kaiserslautern, Germany; Deartment of Physics and Astronomy, University of Manitoba, Winnieg, Manitoba R3T N, Canada M. Maiti Deartment of Physics and Research Center OPTIMAS, TU Kaiserslautern, D Kaiserslautern, Germany; Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna, Moscow Region, Russia N. P. Konstantinidis Deartment of Physics and Research Center OPTIMAS, TU Kaiserslautern, D Kaiserslautern, Germany; Max-Planck-Institut für Physik komlexer Systeme, 1187 Dresden, Germany N. Sedlmayr Deartment of Physics and Research Center OPTIMAS, TU Kaiserslautern, D Kaiserslautern, Germany; Institute de Physique Théorique, CEA/Saclay, Orme des Merisiers, 9119 Gif-sur-Yvette Cedex, France Abstract. We resent a detailed study of the fidelity, the entanglement entroy, and the entanglement sectrum, for a dimerized chain of sinless fermions a simlified Su-Schrieffer-Heeger (SSH) model with oen boundary conditions which is a well-known examle for a model suorting a symmetry rotected toological (SPT) hase. In the non-interacting case the Hamiltonian matrix is tridiagonal and the eigenvalues and -vectors can be given exlicitly as a function of a single arameter which is known analytically for odd chain lengths and can be determined numerically in the even length case. From a scaling analysis of these data for essentially semi-infinite chains we obtain the fidelity suscetibility and show that it contains a boundary contribution which is different in the toologically ordered than in the toologically trivial hase. For the entanglement sectrum and entroy we confirm redictions from massive field theory for a block in the middle of an infinite chain but also consider blocks containing the edge of the chain. For the latter case we show that in the SPT hase additional entanglement as comared to the trivial hase is resent which is localized at the boundary. Finally, we extend our study to the dimerized chain with a nearest-neighbour interaction using exact diagonalization, Arnoldi, and densitymatrix renormalization grou methods and show that a hase transition into a toologically trivial charge-density wave hase occurs. Keywords: Dimerized chain, Su-Schrieffer-Heeger model, Peierls transition, symmetry rotected toological order, fidelity, entanglement sectrum

2 CONTENTS Contents 1 Introduction Diagonalization of the SSH model 4.1 Periodic boundary conditions and N even Oen boundary conditions and N odd Oen boundary conditions and N even Bulk and boundary fidelity suscetibility Periodic boundary conditions Oen boundary conditions and N odd Oen boundary conditions and N even Entanglement entroy General setu A block in an infinite chain A block at the end of a semi-infinite chain Entanglement sectra 5.1 A block in an infinite chain A block at the end of a semi-infinite chain The SSH model with nearest-neighbour reulsion Fidelity suscetibility Entanglement sectra LCRG results for infinite chains Summary and Conclusions 3 1. Introduction The Su-Schrieffer-Heeger (SSH) model has been introduced as a tight-binding model to describe conducting olymers such as olyacetylene [1, ]. A simlified version of the SSH Hamiltonian for sinless fermions and a static lattice dimerization is given by N 1 H = t [1 + ( 1) j δ](c j+1 c j + h.c.) + U j=1 N 1 (n j 1/)(n j+1 1/) (1.1) where t is the hoing amlitude, δ the dimerization arameter, U the nearestneighbour reulsion, and N the number of lattice sites. The creation oerator for a sinless fermion at site j is denoted by c j, and n j = c j c j is the occuation number oerator. The model ossesses time reversal symmetry T and article-hole symmetry C with T = 1 and C = 1 and thus, for U =, belongs into the BDI class in a classification scheme of single-article Hamiltonians [3]. We will concentrate first on the non-interacting case, U =, and ostone a discussion of the interacting case to section 6. For half-filling and any finite δ (dimerized case) the excitation sectrum is then gaed with a ga E δ in the thermodynamic limit. It has also been realized early on that the excitations in the j=1

3 CONTENTS 3 dimerized case are toological solitons and anti-solitons (domain walls between the two ossible dimerization atterns) which carry exactly half an electron charge [, 4]. In recent years, the SSH model has also attracted interest as one of the simlest examles for a model with a symmetry rotected toological (SPT) hase [5, 6]. The toologically non-trivial hase for the Hamiltonian (1.1) with N even is realized for δ >. It can only be transformed into a toological trivial hase by either breaking the symmetries which rotect it or by closing the excitation ga. Thus galess edge modes have to be resent at the boundary of an SSH chain with δ > and the toological trivial vacuum. An oen SSH chain with δ <, on the other hand, is toologically trivial and there are no edge modes. Instead of distinguishing the two hases by the resence or absence of edge modes, one can equivalently consider the Zak-Berry hase γ (bulk-boundary corresondence) which is a Z bulk invariant with γ = π in the SPT hase and γ = in the toologically trivial hase [7, 8, 9, 1]. The comlexity of the ground state of a quantum many-body system and its resonse to changes in the microscoic arameters of the Hamiltonian can be characterized by its entanglement roerties and by the fidelity, resectively. An interesting question is the relation between these quantum information measures and SPT order. For model (1.1) without interactions the von-neumann entanglement entroy as a function of dimerization δ in a eriodic chain, obtained after tracing out one half of the system, has been investigated in Ref. [9]. A detailed study of the entanglement entroy and sectrum as a function of the block size in the SPT and toological trivial hases has, however, not been erformed yet. While it is clear that the ground state in the SPT hase is still only short-ranged entangled so that the bulk entanglement roerties are the same as in the trivial hase, additional entanglement can arise at the boundaries. Studying these boundary contributions to the entanglement is one of the main goals of our study. The fidelity suscetibility χ near hase transitions has also been studied in various models with toological hases [14, 15], however, these studies concentrated on the bulk contribution χ. For an oen system there will in addition also be a boundary contribution, χ = χ + χ B /N, where N is the number of lattice sites, which to the best of our knowledge has not been investigated for models with SPT hases before. It is exactly this boundary contribution χ B which should be sensitive to the resence or absence of edge states. Our aer is organized as follows: In section we review the exact solution of the SSH chain for eriodic boundary conditions (PBC) and oen boundary conditions (OBC). In section 3 we calculate the fidelity suscetibility and analyze the boundary contribution χ B both in the SPT and in the toologically trivial hase. In section 4 we first briefly review known results for the entanglement entroy in critical and gaed quantum chains and the general formalism to calculate the entanglement sectrum for free fermion models. In the analysis of the obtained data for the SSH model we then show, in articular, that redictions from massive field theory for the scaling of with block size are not valid in the SPT hase. In section 5 we study the entanglement sectra and comare with known analytical results in the limit of large block sizes. By using exact diagonalization, Arnoldi, and density-matrix renormalization grou methods we then extend our investigations in section 6 to the SSH model with a Note, however, that tracing out art of the system creates virtual edges so that the entanglement sectra will have different degeneracies deending on whether a strong or a weak bond is cut. The degeneracies of the entanglement sectrum can be used to classify toological hases in onedimensional interacting systems [11, 1, 13].

4 CONTENTS 4 nearest-neighbour Coulomb reulsion. On the basis of the entanglement roerties we identify, in articular, a transition into a toologically trivial charge-density wave hase. In section 7 we summarize our main results and conclude.. Diagonalization of the SSH model We will briefly review the stes required to diagonalize the SSH model with eriodic and oen boundary conditions and give the eigenenergies and eigenvectors. For OBC we will consider both the case where the number of sites N is even and the case where it is odd. The calculations to obtain the fidelity suscetibility and the entanglement sectrum and entroy are then exlained in section 3 and sections 4-5, resectively..1. Periodic boundary conditions and N even This is the simlest case to handle. Taking the doubling of the unit cell due to the dimerization into account and erforming a Fourier transform, the Hamiltonian (1.1) can be written as H = t k ( c 1k c k )T ( cos k + iδ sin k cos k iδ sin k ) ( c1k c k ) (.1) where k = πn/n and n = 1,, N/. This is straightforwardly diagonalized. The eigenvalues are λ ± k = ±t cos k + δ sin k = ± t 1 + δ + (1 δ ) cos k (.) with a ga at the Fermi oints k F = ±π/ given by E = 4t δ. The new oerators in which the Hamiltonian is diagonal, H = ( ) λ + k α k α k β k β k, (.3) k are given by α k = 1 (Ac k1 + c k ), β k = 1 ( Ac k1 + c k ) (.4) with A = λ + k /[t(cos k + iδ sin k)]. The ground state of the Hamiltonian (.3) at half-filling is obtained by filling u the valence band, Ψ = k α k... Oen boundary conditions and N odd While in the PBC case only two momentum modes are couled so that we are left with the diagonalization of a matrix, see equation (.1), oen boundary conditions instead lead to a couling of all momentum modes so that a Fourier transform does not diagonalize the roblem. In real sace, on the other hand, the Hamiltonian is a symmetric tridiagonal matrix with zero diagonal elements and suerdiagonal and subdiagonal elements which alternate between 1 ± δ. The eigenvalues and -vectors are exlicitly known if the number of lattice sites is odd [16]. The eigenvalues are λ =, λ ± k = ± t 1 + δ + (1 δ ) cos θ k (.5)

5 CONTENTS 5 where θ k = πk/(n + 1) and k = 1,, (N 1)/. The sectrum thus contains a single zero energy mode. If we arametrize each eigenvector as ψ k = (x k 1, x k,..., x k N ) T, (.6) then the not-normalized elements of the eigenvector for the zero energy mode are ( x n 1 = 1 δ ) n 1, n = 1,, N + 1 ; x n =, n = 1,, N 1. (.7) 1 + δ For the other modes one finds x k n 1 = 1 + δ 1 δ sin x k n = λ± k 1 δ sin ( πk(n 1) N + 1 ) ( πkn N + 1 ) ( ) πkn + sin, n = 1,, N + 1 N + 1, n = 1,, N 1. (.8) Deending on the sign of the dimerization δ the edge mode is thus localized either at the right or the left end of the chain..3. Oen boundary conditions and N even For an even number of lattice sites we exect two edge states for δ > and no edge states for δ <. The eigenenergies λ ± k can still be written as in equation (.5), however, the arameters θ k are no longer known analytically and instead have to be determined by numerically solving the imlicit equation [16] with (1 δ)t (θ k, N/) + (1 + δ)t (θ k, N/ 1) =, (.9) T (θ k, n) = sin[(n + 1)θ k] sin θ k. (.1) The eigenvector corresonding to each of the N solutions θ k of equation (.9) is then given by x k n 1(δ) = 1 + δ 1 δ T (θ k, n ) + T (θ k, n 1) x k n(δ) = λ± k 1 δ T (θ k, n 1), (.11) where n takes values from 1,, N/. While for δ < all N states are extended, localized edge states are ossible for δ > corresonding to comlex solutions of equation (.9). More recisely, we find that the system has N real solutions and two comlex solutions if δ > δ c where δ c = N, (.1) while all solutions are real if δ < δ c. Physically, this can be understood as follows: The extended states in a finite system are rotected against a small dimerization by the finite size ga E of order 1/N. Only if the ga E δ induced by the dimerization overcomes the finite size ga can the sectrum change qualitatively. Alternatively, one can also think about this roblem in terms of the localization length ξ loc 1/δ of the edge state, see equation (.7), which has to become smaller than the system size in order to allow for localized states. In the thermodynamic limit, N, we have δ c so that the hase transition between the toological trivial and the SPT hase

6 CONTENTS 6 1 λ k +, λk δ Figure 1. Single-article sectrum λ ± k, equation (.5) with θ k determined by equation (.9), of an oen SSH chain as a function of dimerization δ for N = 5 sites. occurs, as exected, at δ =. The eigensectrum for a finite system as a function of dimerization shown in figure 1 makes it clear that there is always a small ga between the two edge modes which will only go to zero for all δ > in the thermodynamic limit. 3. Bulk and boundary fidelity suscetibility At a quantum hase transition between different states of matter we exect the many-body wave function to change. A measure for this change is the fidelity F (δ, ɛ) = Ψ(δ) Ψ(δ+ɛ) where ɛ arameterizes the difference in dimerization between the two states with F (δ, ) 1. Similar to the Anderson orthogonality catastrohe [17] we exect the overla between the two many-body wave functions to vanish exonentially with system size N for ɛ. It is therefore aroriate to consider the fidelity density defined as f(δ, ɛ) = 1 ln F (δ, ɛ). (3.1) N Since f(δ, ) is a minimum, the fidelity has the small ɛ exansion f(δ, ɛ) = χ(δ)ɛ + O(ɛ 3 ) where χ(δ) = 1 f ɛ (3.) ɛ= is the fidelity suscetibility. At a quantum hase transition, this quantity shows universal scaling behaviour [18, 19,, 1] and can thus be used to characterize the transition. Here we want to use the fidelity suscetibility to characterize the transition from the toologically trivial into the SPT hase in the SSH model (1.1). In an exansion in 1/N we can write χ = χ + χ B /N + O(N ). While the bulk contribution

7 CONTENTS 7 χ = lim N χ will be the same in both hases, we exect that the boundary fidelity suscetibility χ B is different. For a system of non-interacting electrons, the many-body wave functions are simle Slater determinants of the single article eigenstates. The fidelity F (δ, ɛ) can thus be exressed as F (δ, ɛ) = det A kl (δ, ɛ), (3.3) with the M M single article overla matrix A kl (δ, ɛ) = ψ k (δ) ψ l (δ + ɛ), (3.4) where M is the number of articles in the ground state. The fidelity suscetibility can now be obtained as χ(δ) = 1 N lim ln F (δ, ɛ) ɛ ɛ (3.5) and thus does not require taking numerical derivatives Periodic boundary conditions In this case there is no boundary so that χ B =. Furthermore, we exect χ(δ) = χ( δ) for N even because changing the sign of δ is in this case equivalent to a simle relabelling of the lattice sites. Using the results of section.1, the analytic calculation of the fidelity at half-filling is straightforward where α k and leads to Ψ(δ) Ψ(δ + ɛ) = k,k α k α k (3.6) = 1 (Ac 1k + c k )(Ã c 1k + c k ) k,k = k AÃ + 1 Ã are obtained by relacing δ δ + ɛ. For the fidelity density this f = 1 N k ln AÃ + 1 N 1 π π/ The bulk fidelity suscetibility can now be calculated in closed form χ = 1 3π π/ ln AÃ + 1 dk. (3.7) sin k dk (cos k + δ sin k) = 1 3 δ (1 + δ ). (3.8) The finite size corrections in the gaed hase are exonentially small, χ(n) χ N ex( N/ξ) for N 1, where ξ is the correlation length, see equation (4.14). At the critical oint δ = the fidelity suscetibility thus diverges as χ 1/ δ. This is consistent with the general scaling theory at a critical oint δ c which redicts χ δ δ c dν where d is the satial dimension and ν the critical exonent related to the divergence of the correlation length at the critical oint. In the case considered here we have δ c =, d = 1, and ν = 1.

8 CONTENTS δ 1 = 6, 4,, 1, χ /N Figure. Fidelity suscetibility for the SSH chain with OBC and an odd number of sites (symbols). The lines are linear fits for 1/N small. Note that the nextleading corrections become more imortant with decreasing δ with the exansion in 1/N breaking down at the critical oint δ =. 3.. Oen boundary conditions and N odd In the oen chain case, we exect a 1/N boundary contribution χ B (surface term) to the fidelity suscetibility. According to the bulk-boundary corresondence, a symmetry rotected toological hase as indicated by a bulk invariant Zak-Berry hase is related to the resence of edge modes. We thus exect that χ B is a measure which can distinguish between a toologically trivial and a non-trivial hase. We start by considering the case of an odd number of sites N where the single article eigenstates are known exactly, see section.. Note that in this case a single edge mode is resent both for δ < and for δ >. In order to obtain both the bulk and the boundary contribution, we use equations (3.3), (3.4), and (3.5) with the exact single article eigenstates given by equations (.7) and (.8). We are interested in the half-filled case, however, for N odd we can only ut either (N 1)/ or (N + 1)/ articles into the system. This means that the zero energy edge mode is either emty or occuied. Because the Hamiltonian is article-hole symmetric, the fidelity is, however, exactly the same in both cases. In figure we show exemlarily the scaling of the fidelity suscetibility for different dimerizations, demonstrating that the leading correction to the bulk suscetibility is indeed of order 1/N. The bulk and boundary suscetibilities extracted from the scaling are shown in figure 3. Note that both are symmetric in δ as they should be because δ δ is u to a relabelling of the lattice sites a symmetry of the Hamiltonian (1.1) if N is odd. The bulk contribution χ is indeendent of the boundary conditions and agrees with the analytical solution (3.8) obtained using eriodic boundary conditions. For the boundary contribution χ B we numerically find a behaviour close to the critical oint which is consistent with χ B 1/ δ η with η. Note, however, that we have two different edges for the N odd case which can give boundary contributions which scale differently. The analysis of χ B is therefore easier

9 CONTENTS (a) (b) 15 χ χ B δ δ Figure 3. Symbols denote (a) the bulk suscetibility, and (b) the boundary suscetibility for N odd extracted from the scaling in figure. The data in (a) are comared with the analytical solution (lines), equation (3.8). In (b) the lines reresent χ B.33 δ obtained by fitting the data for δ.. in the N even case considered next where either both edges have a localized state or both are trivial Oen boundary conditions and N even In the N even case we can again calculate the overla matrix (3.4) from the single article eigenstates, equation (.11). However, this time the arameters θ k have to be determined by the imlicit equation (.9). Esecially for the bound states where θ k becomes imaginary this requires a high accuracy making the calculations numerically challenging for large system sizes. The boundary suscetibility obtained in this case is shown in figure 4. Note that χ B in the toological trivial hase (δ < ) is qualitatively different from χ B in the SPT hase (δ > ). In articular, the sign is different. The data in figure 4(b), furthermore, indicate that also the exonents of, what seems to be a ower-law divergence at the critical oint, are different by about 5%. The resence or absence of edge states is thus clearly visible in χ B. For large chain lengths N the edges should become indeendent so that χ B is simly the sum of the contributions from each of the two edges. In articular, this imlies that lim [χ B(N even, δ) + χ B (N even, δ)]/ = χ B (N 1 odd, ±δ). (3.9) N This behaviour is confirmed by the fits shown in figure 4(b). The data for the N odd case are thus better fitted by the sum of two ower laws with different exonents rather than by the single ower law used in figure 3. Finally, we can also study the finite size scaling of the fidelity suscetibility right at the critical oint, see figure 5. While both data sets for N even and N odd show

10 CONTENTS 1 3 (a) 1 1 (3) () (b) 1 χ B 1 1 (1) χ B δ δ Figure 4. (a) χ B in the thermodynamic limit for an even number of sites (circles). Lines are a guide to the eye. (b) Fits of the data (circles) for (1) N even, δ < : χ (1) B.177/ δ.73, () N even, δ > : χ () B.98/δ1.979, and (3) N odd: data and function χ (3) B = (χ(1) B + χ() B )/, see equation (3.9). For the fits the data with δ.1 have been used. 1 8 χ N Figure 5. Finite size scaling of χ at the critical oint δ = for N even (circles) and N odd (squares) with oen boundary conditions. The lines are fits with χ.1791n (N even), and χ.313n N/3 (N odd).

11 CONTENTS 11 a scaling χ N consistent with the general scaling rediction χ N /dν 1 [], the refactors are different. To summarize, we have established that the boundary suscetibility χ B behaves qualitatively different in the SPT (δ > ) than in the toologically trivial hase (δ < ). In articular, χ B (δ < ) < and χ B (δ > ) >. 4. Entanglement entroy The entanglement entroy is the von-neumann entroy of a reduced density matrix S A = Tr(ρ A log ρ A ) (4.1) with ρ A = Tr B ρ where ρ = Ψ Ψ is the full density matrix with Ψ being the ground state of the many-body Hamiltonian. The system has been divided u into two blocks A and B. It is easy to rove that S A = S B. The entanglement entroy is an effective measure to characterize the ground states of many-body systems. A gain in entanglement entroy can, furthermore, be a driving force for a Peierls transition [] leading to a dimerization of the chain as considered here. For one-dimensional critical quantum systems it has been shown that shows a universal deendence on the length l of the block A. For a system of size N with PBC one finds S CF T ent = c 3 ln [ N π sin πl N ] + s 1 (4.) where c is the central charge of the underlying conformal field theory (CFT) in the scaling limit, and s 1 a non-universal constant [3]. In the thermodynamic limit N, the CFT formula (4.) reduces to S CF T ent = c 3 ln l + s 1. (4.3) By a conformal maing one can show that a similar relation also holds for a system at small but finite temeratures [4]. For the undimerized chain of sinless fermions, equation (1.1) with δ =, one finds c = 1 and s 1.76 [5]. These results are altered though if one changes the boundary conditions from eriodic to oen. In this case the entanglement entroy of a block which includes the end of the chain reads [3] S CF T ent = c 6 ln [ N π ] πl sin + ln g + s 1 / (4.4) N where ln g is the boundary entroy related to the ground state degeneracy of the considered system [6]. For the critical undimerized chain we have g = 1. In numerical studies the scaling law for OBC, equation (4.4), is, however, often obscured by corrections coming from leading oerators of the underlying CFT [7]. For critical systems described by Luttinger liquid theory one finds, for examle, S CF T ent + cos(k F l)l K (4.5) where k F is the Fermi momentum and K the Luttinger liquid arameter [8, 9, 3, 31, 3]. In a massive hase one exects, on the other hand, that the entanglement entroy saturates once the block size become larger than the correlation length ξ. If the system

12 CONTENTS 1 is in a massive hase but close to a critical oint then also the leading corrections can be calculated and one finds for a block in an infinite chain = c 3 ln(ξ 1/a) + U 1 K (l/ξ α ) (4.6) 8 where U and a are non-universal constants, ξ α the correlation lengths (with ξ 1 being the largest), and K the modified Bessel function [33]. This tye of scaling behaviour has recently been confirmed numerically for certain sin models [34] General setu In the following, we want to study the entanglement entroy of the SSH model (1.1). We will consider two cases: First, a block in an infinite chain. Here we can simly use PBC which makes it ossible to erform the thermodynamic limit exactly. This case will rovide a numerical check of the formula (4.6). Second, we will consider the case of a block at the end of a semi-infinite chain. Here we can use either odd chain lengths where the eigenstates are known exlicitly, or even chain lengths. The odd chain length has an edge state which is either filled or emty, which are equivalent by article-hole symmetry. In contrast at half-filling a finite oen system with an even number of sites in the SPT hase has two edge states, only one of which is occuied at half filling. However, the boundary states at the edges in the thermodynamic limit are a suerosition of states localized at the left and the right edge. Thus we can also consider a semi-infinite chain with a artially occuied edge state. Since the correlation length ξ δ 1 is finite for δ it suffices to consider chain lengths N ξ in order to be effectively in the thermodynamic limit. We are articularly interested in the difference in entanglement entroies with and without a localized state being resent at the boundary of the chain. For a free fermion model there is a general setu to obtain the eigenvalues of the reduced density matrix ρ A [35, 36]. The main observation is that ρ A can be reresented as the exonential of a free-fermion oerator H A the so-called entanglement Hamiltonian with ) ρ A = 1 Z ex( H A) = 1 l ( Z ex ɛ k a k a k (4.7) and Z = Tr ex( H A ). On the other hand, all roerties of a free-fermion system can be obtained from the matrix C of two-oint correlation functions because all higher correlation functions can be exressed through the two-oint functions by using Wick s theorem. The matrix elements of C for all two-oint functions within the considered block A are, in articular, determined by the reduced density matrix ρ A, i.e., matrix elements are given by α=1 k=1 C ij = Tr(ρ A c i c j). (4.8) The eigenvalues ζ k of the correlation matrix can then be related to the eigenvalues ɛ k of the entanglement Hamiltonian in equation (4.7) and one finds ζ k = (e ɛ k + 1) 1. (4.9) According to equation (4.7) the l -many eigenvalues z of ρ A itself are now obtained by considering all ossible fillings of the l-many levels ɛ k z = 1 ( ) ex ɛ k n () k (4.1) Z n ()

13 CONTENTS 13 where n () = (n () 1,, n() l ) is a vector of single-level occuation numbers, n () k {, 1}, and Z = z. The entanglement entroy is then given by = z ln z. (4.11) 4.. A block in an infinite chain Using the results from section (.1) for a chain with PBC at half-filling, the two-oint correlation function can be easily calculated and the thermodynamic limit erformed exactly. For even distances one finds C ij = δ ij /, i.e., only the onsite correlation function (occuation number) is nonzero. For odd distances the result is c n c m+1 = 1 π π/ π/ cos k cos k(m + 1 n) + δ sin k sin k(m + 1 n) dk (4.1) cos k + δ sin k and the correlation function c n+1 c m is obtained by relacing δ δ. For the non-dimerized case, this reduces to the well-known result c nc m = 1 sin [ π (m n)]. (4.13) π m n By methods of steeest descent we can also obtain the asymtotic behaviour of the correlation function (4.1) at large distances. In articular, we find that the dominant correlation length is given by ξ 1 ξ1 1 = 1 ln 1 + δ 1 δ. (4.14) There should now be three different regimes deending on the ratio of block size l and largest correlation length ξ ξ 1 : (1) For ξ l 1 the entanglement entroy will essentially follow the scaling law for a critical system, equation (4.), as obtained by CFT. Exemlarily, we resent in figure 6 results for for small dimerizations and odd block lengths l. Indeed, we see that the smaller the dimerization is the longer equation (4.) aroximately holds as a function of block size. () For ξ l, on the other hand, the entanglement entroy will settle to a constant value = c 3 ln(ξ/a) + U where both constants a and U are non-universal. This behaviour is shown exemlarily in figure 7(a). In figure 7(b) we show, furthermore, that the saturation value deends on whether two weak bonds (l even), a weak and a strong bond (l odd), or two strong bonds (l even) are cut. We have chosen the block for l even in such a way that δ < corresonds to cutting two weak and δ > to cutting two strong bonds. For block sizes l ξ each non-trivial edge caused by the cutting of a strong bond gives an indeendent contribution S edge to the entanglement entroy. More recisely, we can define the edge contribution as S edge = lim (l, l odd) (l + 1, l + 1 even), (4.15) l which is shown as a function of dimerization in figure 8. In articular, in the limit δ 1 the entanglement entroy will be dominated entirely by the edge contribution so that Sent δ 1 (l even) = S δ 1 edge, Sδ 1 ent (l odd) = S δ 1 edge, and Sδ 1 ent (l even) with Sedge δ 1 = ln. The insets of figure 8 are consistent with a ower law behaviour of S edge for small δ and for δ 1. (3) Finally, there is also the regime ξ l where the formula (4.6) should hold which gives the leading correction to the constant value obtained in the limit ξ l. In

14 CONTENTS δ = 1, 3, 5, 1,, l Figure 6. for a block in an infinite chain and small dimerizations (symbols). The solid line is the CFT result, equation (4.), using the value s 1 =.76 obtained for the non-dimerized chain [5]. The dashed lines are a guide to the eye. 1 (a).8 δ=.9,.7,.5,.3, (b) S edge S edge l Figure 7. (a) for a block with an odd number of sites in an infinite chain and large dimerizations. (b) for l even and δ =. (circles), l odd (squares), and l even and δ = +. (triangles). Each non-trivial edge gives a contribution S edge to. The dashed lines are a guide to the eye.

15 CONTENTS 15.7 S edge (a) (b) δ ln() - S edge. S edge δ δ Figure 8. Edge contribution (4.15) to the entanglement entroy of a block in an infinite chain. The line marks the limiting value Sedge δ 1 = ln. The dashed line is a guide to the eye. Inset (a): Fit for δ 1 with S edge δ.694. Inset (b): Fit for δ 1 with S edge ln.73(1 δ ) the case of the dimerized model considered here, the leading correlation length ξ will show u twice: once with wave vector k = π/ and once with wave vector k = π/ [37, 38]. Therefore, to leading order, we exect the entanglement entroy to scale as = (l ) 1 4 K (l/ξ) (4.16) with ξ as given in equation (4.14). Fixing (l ) by the numerical data for large block sizes, there is thus no free fitting arameter. The comarison in figure 9 with data for δ =., in which case ξ 5, confirms this scaling rediction. For large arguments, the Bessel function asymtotically scales as K (x) π x ex( x). The inset of figure 9 demonstrates that this is indeed the correct scaling of the entanglement entroy for large block sizes. The data for the three different cuts, in articular, collase onto a single line after subtracting the constant art (l ) A block at the end of a semi-infinite chain To calculate the entanglement entroy for a block which includes the end of a semiinfinite chain, we can use the analytic results for the single-article eigenstates of a chain with OBC and N odd given in section. and the numerical results for N even from section.3. In both cases, the matrix elements of the correlation matrix C ij, see equation (4.8), can be exressed as C ij = (Ψ k i ) Ψ k j (4.17) k occuied with the sum running over all occuied single-article states. Diagonalizing the correlation matrix then gives access to the single-article eigenenergies ɛ k of the entanglement Hamiltonian using the relation (4.9). Results for small dimerizations are shown in figure 1. While the entanglement entroy for δ <, where a strong bond terminates the semi-infinite chain, behaves

16 CONTENTS l l Figure 9. (circles) for a block in an infinite chain at δ =. with (from bottom to to) l even and δ <, l odd, and l even and δ > comared to the scaling rediction (4.16) (lines). The inset shows that all three curves for l/ξ 1 collase onto a single line given by the asymtotics of the Bessel function (see text). qualitatively similar to the case of a block in an infinite chain (see figure 6) we find a much more comlicated scaling for δ > and N odd. In this case a localized boundary state is resent leading to a cometition between an increase of the boundary contribution to with increasing δ and a decrease of the bulk contribution due to the decreasing correlation length ξ 1/δ. As a result, the curves for different dimerizations cross, see figure 1(b). The final ossibility is to consider the semi-infinite chain with a localized edge state which is not fully occuied, as realized in a chain N even at half-filling, see figure 1(c). As for the fully occuied edge state there is a cometition between the increasing boundary entroy and decreasing bulk entroy with δ. However, due to the higher degeneracy associated with the artially occuied edge state, see also the entanglement sectra in the following section, the curves for the same dimerizations will cross at a much larger sub-block size, than for a fully occuied edge state. Therefore we show in figure 1(c) data for dimerizations which are an order of magnitude larger than in the other two cases shown in figure 1(a,b). For large block size l and finite dimerization we again exect that saturates. Let us first concentrate on the case δ < shown in figure 11. In this case is monotonically increasing. We can thus try to fit the data for large block sizes by = (l ) 1 8 K (l/ ξ) (4.18) similar to the CFT result (4.6) but with ξ as a free arameter. Doing so we obtain very good fits with arameters as given in table 1. The correlation length ξ obtained from the fits seems to be about a factor smaller than the bulk correlation length ξ given in equation (4.14). Note, however, that the heuristic fitting formula (4.18) is too naive to be able to fully cature the behaviour of. A two-oint correlation function near a boundary does not only deend on the distance between the two oerators but

17 CONTENTS 17 CFT - CFT - CFT δ = 1, 4, (b) 1 3 δ = 1, 4, (a) (c) 1 δ = 1, 4, l 1 5 Figure 1. for a block at the end of a semi-infinite chain at small dimerizations with l even for (a) δ < where the boundary is trivial, (b) δ > where a localized boundary state is resent and filled, and (c) δ > for N even where a localized boundary state is resent but only half-filled. In all cases the CFT result (4.4) in the non-dimerized case with s 1 =.76 has been subtracted. The dashed lines are a guide to the eye..8.7 δ=. δ= δ= l l Figure 11. for a block at the end of a semi-infinite chain with l even, fitted by equation (4.18) with the fitting arameter ξ given in table 1. The inset shows that the scaling is still roughly consistent with the Bessel function asymtotics.

18 CONTENTS 18 Table 1. Parameters in equation (4.18) used to fit the data for large block sizes shown in figure 11. δ (l ) ξ ξ ξ/ ξ l δ =. δ = δ =.6.6 δ = l Figure 1. for a block at the end of a semi-infinite chain with δ > and a fully occuied boundary state. Shown are data for even (circles) and odd block sizes (squares). also on the distance of the oerators from the boundary so that a third length scale has to aear. Let us now return to the case δ > with the fully occuied boundary state where figure 1(b) already indicates that the scaling behaviour of the entanglement entroy with block size is comletely changed. Figure 1 shows that the cometition between boundary and bulk contributions leads to a maximum in which becomes sharer and whose osition shifts to smaller block sizes with increasing dimerization. Consequently, the constant value (l ) is now aroached from above. Clearly a scaling as in (4.18) is no longer valid. To verify that the additional entanglement entroy is indeed caused by the boundary state we can define S bound (l) = Sent δ> (l, l odd) S δ< ent (l + 1, l + 1 even). (4.19) For both entroies in equation (4.19) a weak bond is cut so that for large block sizes l any remaining difference is directly related to the localized boundary state resent for δ >. Figure 13 shows that S bound (l) is exonentially decaying with block size. According to equation (.7) a non-trivial boundary leads to an exonentially localized zero energy state with a localization length ξ 1 loc = ln ( 1 + δ 1 δ ). (4.) The exonential fits in figure 13 using this correlation length are indeed in excellent agreement with the data for S bound (l) confirming that the boundary state leads to

19 CONTENTS 19 S bound δ =. δ = S bound S bound l δ =.6 δ = S bound l Figure 13. S bound as defined in equation (4.19) for different δ > where the boundary state is fully occuied. The lines dislay S bound = a ex( l/ξ loc ) with the localization length (4.) and a used as the only fitting arameter δ =. 1 δ = l δ =.6 1. δ = l Figure 14. for δ > and a artially filled boundary state. Shown are data for even block sizes (circles). The lines are exonential fits, see equation (4.1). additional entanglement which is exonentially localized at the boundary. For δ >, N even, and half-filling the boundary state is only half-filled. As shown in figure 14 this qualitatively changes the behaviour of as comared to the case where the edge state is comletely filled. In articular, there is no longer a maximum yet the scaling with block size is also not described by a modified Bessel function as was the case for δ <. Instead, the scaling seems urely exonential (l) = (l ) a ex( l/ξ) (4.1)

20 CONTENTS 1 S bound S bound δ =. δ = l δ =.6 δ = S bound S bound 4 6 l Figure 15. S bound as defined in equation (4.19) for δ > with a artially filled boundary state. The lines dislay S bound = a ex( l/ξ loc ) with the localization length (4.) and a used as the only fitting arameter. with ξ being the leading correlation length, see equation (4.14), and a a fitting arameter. The boundary contribution of the artially filled edge state can again be calculated from equation (4.19) and is exonentially decaying with block size similar to the case of the comletely filled edge state, see figure 15. To summarize, we have shown that the resence of localized boundary states in a massive hase comletely changes the scaling of (l) with block size l for a block at the end of a semiinfinite chain. The boundary state leads, in articular, to additional entanglement which is exonentially loacalized. Further insight into the different entanglement roerties of a block with or without a localized boundary state can be gained by directly investigating the entanglement sectrum. 5. Entanglement sectra Instead of just considering the entanglement entroy, which is a ma from the l dimensional sace of eigenvalues of ρ A into the real numbers, it is often instructive to consider the sectrum of ρ A the so-called entanglement sectrum directly A block in an infinite chain As for the entanglement entroy we can use the analytical result for the elements of the correlation matrix in the thermodynamic limit to determine the sectrum of the reduced density matrix. In figure 16 the logarithm of the eigenvalues z of the reduced density matrix, see equation (4.1), for the three different cuts ossible are shown. The ictograms in the figure show whether the block size is even or odd (but not the actual block size!) and if strong (double lines) or weak (single dashed lines) are cut. While the structure is comlicated for finite block sizes l, the distribution of eigenvalues becomes very simle in the limit l. In this case, the single article energies ɛ k of the entanglement Hamiltonian (4.7) become equally saced and the level sacing ɛ can be obtained analytically using corner transfer methods [39, 4]. One

21 CONTENTS 1 -ln(z ) (a) (b) (c) (18) () (6) (8) (1) 1 5 (17) (6) (8) (8) (1) (16) ε ε ε Figure 16. Entanglement sectrum for a block in an infinite chain at δ =.1 with (a) odd block size, (b) even block size cutting two weak bonds, and (c) even block size cutting two strong bonds. The black stars are data for l = 19 (l = ), the red crosses for l = 69 (l = 7). The dashed lines are the single article levels in the thermodynamic limit, the numbers in bracket give the degeneracy d k of each level. finds ɛ = πi(x )/I(x) with x = (1 δ)/(1 + δ), (5.1) x = 1 x, and I(x) being the ellitic integral of the first kind. This is the correct level sacing for odd block size as well as for even block size if two weak bonds are cut. Cutting two strong bonds, on the other hand, leads to a level sacing ɛ, see figure 16. The eigenvalues z of the reduced density matrix in the thermodynamic limit can then be obtained by filling u the single article levels. The dominant eigenvalue of ρ A is obtained by filling u the Fermi sea of negative eigenvalues ɛ k. In the case shown in figure 16(b) where the block size is even and two weak bonds are cut this leads to a unique dominant eigenvalue. In the case where either one [figure 16(a)] or two [figure 16(c)] strong bonds are cut the lowest eigenvalue is, however, degenerate. This can be understood as follows: Cutting a strong bond leads to an edge sin. For the case l odd where one strong bond is cut this leads to an eigenenergy ɛ (l+1)/. This level can now either be emty or filled leading to a two-fold degeneracy of the leading eigenvalue. For the case l even where two

22 CONTENTS strong bonds are cut, eigenvalues ɛ l/ = ɛ l/+1 aear which aroach zero in the thermodynamic limit. Consequently, the dominant eigenvalue of ρ A becomes four-fold degenerate. Thus the low-energy entanglement sectrum allows to directly identify the number of edge states caused by cutting a strong bond. Equation (5.1) only gives the level sacing. In order to obtain the allowed values for the eigenvalues z of ρ A one has to use, in addition, the normalization condition 1 = Tr ρ A = ( ) z = d k e ɛ kɛ ɛ = ln d k e kɛ. (5.) k= Here d k is the number of degenerate eigenvalues {z } for which z = ex( ɛ kɛ) with k fixed. The levels ɛ + kɛ or ɛ + kɛ, resectively, with k =, 1,, 3, are shown as dashed lines in figure A block at the end of a semi-infinite chain While the case of a block in an infinite chain is well understood, the interesting case of the entanglement sectrum for a block at the end of a dimerized semi-infinite chain with or without a localized boundary state has not been considered so far. We have already shown that a localized boundary state causes additional entanglement which is also localized at the boundary, see figure 13. Here we want to further investigate this case by directly looking at the entanglement sectrum. For N odd, we can distinguish four cases now: δ < (toological trivial case) or δ > (SPT case) with l either even or odd. In addition, we have to consider the two cases l even or odd with N even and δ >, where the edge state is only artially filled, searately. We start with the toological trivial case δ <. Data for even and odd block sizes with δ =.1 are shown in figure 17. For even block sizes, see figure 17(a), the leading eigenvalue is unique and the searation of the single-article levels is given by equation (5.1) as in the case of an even block in an infinite chain where two weak bonds are cut, see figure 16(b). However, the degeneracy structure of the higher levels is different which exlains the differences in the entanglement entroy discussed in the revious subsection. For the case of an odd block size, shown in figure 17(b), the degeneracy structure of the sectrum is the same as for a block of odd length in an infinite chain, see figure 16(a), however, the level sacing in the semi-infinite chain is ɛ instead of ɛ for the infinite chain. Note that the finite size corrections for a block at the end of a semi-infinite chain in the toological trivial case are of similar magnitude as for the infinite chain, i.e., for the lowest single-article levels the data for the smaller block sizes l = 19, and the larger ones, l = 69, 7, hardly deviate from each other. This is comletely different in the SPT hase the case δ =.1 with N odd is shown as examle in figure 18 where finite size corrections are much larger. In the limit l, on the other hand, both the level sacing and the degeneracies agree between the cases δ <, l even [figure 17(a)] and δ >, l odd [figure 18(a)] as well as δ <, l odd [figure 17(b)] and δ >, l even [figure 18(b)]. This means that the localized boundary state resent for δ > does not lay any role in the limit l and that the sectra in the trivial and the SPT hase become identical in this limit. The dramatic change in the entanglement sectra for small block sizes if a localized state is resent which leads to the non-monotonic behaviour of the entanglement entroy shown in figure 1 can be understood as follows: The nontrivial edge leads to an exonentially localized zero energy state, equation (.7), with a localization length as given in equation (4.). k=

23 CONTENTS 3 5 (a) (1) (b) (18) -ln(z ) () () (1) (5) ε (6) (9) (1) (1) (6) ε () Figure 17. Entanglement sectrum for a block at the end of a semi-infinite chain at δ =.1 with (a) even block size, and (b) odd block size. The black stars are data for l = 19 (l = ), the red crosses for l = 69 (l = 7). The dashed lines are the single article levels in the thermodynamic limit, the numbers in bracket give the degeneracy of each level. For the case δ =.1 considered here the localization length is ξ loc 1/δ = 5 lattice sites. The low-energy art of the entanglement sectrum will thus become similar to the one without the localized state if l ξ loc while it will be substantially modified for l ξ loc. The modification of the entanglement sectrum and the entanglement entroy is thus an effect of finite block size which is directly tied to the aearance of localized states in the SPT hase. The final ossibility is to consider the case N even, δ > where the boundary state is only artially filled, see figure 19. In this case the degeneracy of the levels is doubled comared to the case with a fully occuied boundary state, shown in 18. Furthermore, corrections to the l limit for finite blocks are smaller, exlaining the different behaviour seen in the entanglement entroy. To summarize, we have shown that the resence of a localized boundary state significantly alters the entanglement sectra for small block sizes. In the limit l, on the other hand, a free fermion sectrum with an equidistant level sacing is recovered. While creating virtual edges by considering a block in an infinite chain gives the same degeneracy of the lowest level than a hysical edge for a corresonding block at the end of a semi-infinite chain, the degeneracies of higher levels as well as the level sacing are, in general, different. This exlains the different scaling roerties

24 CONTENTS 4 5 (a) (1) (b) (18) -ln(z ) () (5) () (1) (6) (9) ε (1) (1) (6) ε () Figure 18. Entanglement sectrum for a block at the end of a semi-infinite chain at δ =.1 with a fully occuied boundary state for (a) odd block size, and (b) even block size. The black stars are data for l = 19 (l = ), the red crosses for l = 69 (l = 7). The dashed lines are the single article levels in the limit l, the numbers in bracket give the degeneracy of each level. of (l) for a block at the end as comared to a block in the middle. 6. The SSH model with nearest-neighbour reulsion Finally, we will investigate the effects of interactions by adding a nearest-neighbour reulsion term, see the Hamiltonian in equation (1.1). It has been shown that toological hases of interacting one-dimensional quantum systems can be classified by their entanglement roerties [11, 1, 13]. Interaction effects in one-dimensional dimerized fermionic systems with SPT hases have already been studied in Ref. [4]. However, in this work a sinful dimerized Hubbard model was considered where for δ 1 (SPT hase), even in the limit of infinite onsite Hubbard reulsion, one can still remove an electron from the edge site without any energy cost. Thus the cases δ < and δ > remain toologically distinct in the resence of a Hubbard interaction. For the sinless model with nearest-neighbour reulsion considered here, the situation is different. For U the system will be forced into a charge-density-wave (CDW) state with every second site occuied so that, even for δ 1, the edges are no longer decouled from the rest of the chain. We therefore exect a hase transition in the

25 CONTENTS 5 (a) (1) (b) (1) (4) 15 (8) 15 -ln(z ) 1 5 () (8) ε 1 5 (8) (1) ε () Figure 19. Entanglement sectrum for a block at the end of a semi-infinite chain with N even at δ =.1 (artially filled boundary state) with (a) odd block size, and (b) even block size. The black stars are data for l = 19 (l = ), the red crosses for l = 69 (l = 7). The dashed lines are the single article levels in the thermodynamic limit, shifted by ln relative to the odd system size case due to the doubling in degeneracy of each level, see the numbers in arenthesis. model (1.1) as a function of interaction U where the fidelity and the entanglement measures will change. We focus first on oen chains with an even number of sites and use exact diagonalization and Arnoldi algorithms to calculate the ground state roerties for system sizes u to N = 4. In order to firmly establish the change of the entanglement roerties at the hase transition we, furthermore, resent light cone renormalization grou (LCRG) [41] results for infinite chains Fidelity suscetibility In figure we show the fidelity suscetibility at different interaction strengths U for δ =.5 [toological hase for U = ] and δ =.5 [trivial hase for U = ]. For U =.5 and U =.5 the fidelity suscetibility behaves qualitatively similar to the non-interacting case with the bulk suscetibility χ being reduced with increasing interaction strength. Furthermore, we still find χ B (δ < ) < and χ B (δ > ) > with χ B given by the sloe of the curves in figure for 1/N small. This is an indication that the two cases remain toologically distinct even for intermediate interaction strengths. For U = 5, however, a drastic change is observed. In this case χ B is negative both for δ < and for δ >. Aarently the boundaries now behave qualitatively similar suggesting that a hase transition has occurred and edge states no longer exist. We