Quantum Phase Transition and Entanglement in Topological Quantum Wires

TL;DR

This paper studies quantum phase transitions in topological quantum wires, revealing that local entanglement nonanalyticities signal topological changes and uncovering a new entanglement pattern transition independent of spectral gap closure.

Contribution

It demonstrates that local entanglement nonanalyticities universally characterize topological phase transitions in one-dimensional noninteracting fermion systems and identifies a novel entanglement pattern transition.

Findings

Local entanglement becomes discontinuous at topological phase transitions.

Universal nonanalyticity observed in SSH and Kitaev chain models.

A new phase transition identified by entanglement pattern change, independent of spectral gap closure.

Abstract

We investigate the quantum phase transition of the Su-Schrieffer-Heeger (SSH) model by inspecting the two-site entanglements in the ground state. It is shown that the topological phase transition of the SSH model is signified by a nonanalyticity of local entanglement, which becomes discontinuous for finite even system sizes, and that this nonanalyticity has a topological origin. Such a peculiar singularity has a universal nature in one-dimensional topological phase transitions of noninteracting fermions. We make this clearer by pointing out that an analogous quantity in the Kitaev chain exhibiting the identical nonanalyticity is the local electron density. As a byproduct, we show that there exists a different type of phase transition, whereby the pattern of the two-site entanglements undergoes a sudden change. This transition is characterised solely by quantum information theory and does not accompany the closure of the spectral gap. We analyse the scaling behaviours of the entanglement in the vicinities of the transition points.

Figures

Equations

$$H=\sum _{n=1}^N\left(t_1{a}_n^{†}{b}_n+t_2{b}_n^{†}{a}_{n+1}+\text{H.c.}\right),$$

$$t_1=1-\lambda ,t_2=1+\lambda$$

$$H=\sum _{k\in \text{BZ}}{\varphi }_k^{†}\left[\vec{h}(k)\cdot \vec{\sigma }\right]{\varphi }_k,$$

$$\vec{h}(k)=(t_1+t_2\cos k)\hat{x}+(t_2\sin k)\hat{y}.$$

$$\frac{i}{4}{c}^T\bar{M}c=\sum _{k\in \text{BZ}}{\varphi }_k^{†}\left[\hat{h}(k)\cdot \vec{\sigma }\right]{\varphi }_k,$$

$$\frac{i}{4}{c}_{{\text{a}}_n{\text{b}}_m}^T{\bar{M}}_{{\text{a}}_n{\text{b}}_m}{c}_{{\text{a}}_n{\text{b}}_m}=\frac{1}{N}\sum _{k\in \text{BZ}}{\left(\begin{array}{c}{a}_n^{†}{e}^{ikn}\\ b_m^{†}{e}^{ikm}\end{array}\right)}^T\hat{h}(k)\cdot \vec{\sigma }\left(\begin{array}{c}{a}_n{e}^{-ikn}\\ b_m{e}^{-ikm}\end{array}\right)={\eta }_{{\text{a}}_n{\text{b}}_m}(a_n^{†}{b}_m+b_m^{†}{a}_n),$$

$${\eta }_{{\text{a}}_n{\text{b}}_m}$$

$$e^{i\phi (k)}$$

$${\rho }_{{\text{a}}_n{\text{b}}_m}=\frac{1}{4}\left(\begin{array}{cccc}1-{\eta }^2& 0& 0& 0\\ 0& 1+{\eta }^2& -2\eta & 0\\ 0& -2\eta & 1+{\eta }^2& 0\\ 0& 0& 0& 1-{\eta }^2\end{array}\right)$$

$${\mathcal{C}}_{{\text{a}}_n{\text{b}}_m}=\max \left\{0,\frac{1}{2}({\eta }_{{\text{a}}_n{\text{b}}_m}^2+2{\eta }_{{\text{a}}_n{\text{b}}_m}-1)\right\},$$

$${\mathcal{C}}_1(\lambda )={\mathcal{C}}_2(-\lambda ).$$

$${\eta }_1(\lambda )=\frac{1}{N}\sum _{k\in \text{BZ}}{\hat{h}}_x(k)\stackrel{N\to \infty }{\to }\frac{1}{2\pi }{\int }_0^{2\pi }{\hat{h}}_x(k)dk.$$

$$\delta (N)\equiv \left|{\eta }_1(0_{+})-{\eta }_1(0_{-})\right|=\frac{2}{N}.$$

$$\begin{array}{rl}{\frac{\partial }{\partial \lambda }F(T,\lambda )|}_{T\to 0}& =N⟨-a_1^{†}{b}_1+b_1^{†}{a}_2+\text{H.c.}{⟩}_{T\to 0}\\ & =N\left[{\eta }_1(\lambda )-{\eta }_2(\lambda )\right],\end{array}$$

$$\frac{\partial }{\partial \lambda }{\eta }_1(\lambda )\simeq \frac{2}{\pi }\log \left[{\left(\frac{e}{2}\right)}^2|\lambda |\right],$$

$${\rho }_1=\frac{1}{N}\sum _{k\in \text{BZ}}{\hat{h}}_z(k)(a_1^{†}{a}_1-a_1{a}_1^{†}),$$

$${\mathcal{C}}_1(\lambda )\simeq -1.476(\lambda -{\lambda }_{+}).$$

$$H_{\text{spin}}=\sum _{n=1}^{N-1}(t_1{a}_n^{†}{b}_n+t_2{b}_n^{†}{a}_{n+1}+\text{H.c.})+t_1(a_N^{†}{b}_N+\text{H.c.})+(-1{)}^{N-1}{t}_2(b_N^{†}{a}_1+\text{H.c.}).$$

$$H=\frac{i}{4}{c}^{T}Mc,$$

$$WM{W}^T=\left(\begin{array}{ccccc}0& {ϵ}_1& & & \\ -{ϵ}_1& 0& & & \\ & & \ddots & & \\ & & & 0& {ϵ}_N\\ & & & -{ϵ}_N& 0\end{array}\right)$$

$$H=\frac{i}{2}\sum _{j=1}^N{ϵ}_j{c}_{2j-1}^{\prime }{c}_{2j}^{\prime }=\sum _{j=1}^N{ϵ}_j\left({a_j^{\prime }}^{†}{a}_j^{\prime }-\frac{1}{2}\right),$$

$${\rho }_0=\prod _{j=1}^N{a}_j^{\prime }{a_j^{\prime }}^{†}=\prod _{j=1}^N\left(\frac{1}{2}-i\frac{1}{2}{c}_{2j-1}^{\prime }{c}_{2j}^{\prime }\right).$$

$$\bar{H}=\frac{i}{4}{c}^T\bar{M}c=\sum _{j=1}^N\left({a_j^{\prime }}^{†}{a}_j^{\prime }-\frac{1}{2}\right),$$

$$\omega ({\rho }_0)=\prod _{j=1}^N\left(\frac{1}{2}-i\frac{1}{2}{\theta }_{2j-1}^{\prime }{\theta }_{2j}^{\prime }\right)=\frac{1}{2^N}\exp \left(-\frac{i}{2}{\theta }^T\bar{M}\theta \right).$$

$$\bar{M}_{jk}=\left\{\begin{array}[]{ll}-i\text{Tr}(\rho_{0}c_{j}c_{k})&\text{ for }j\not=k\\ 0&\text{ for }j=k\end{array}\right.$$

$$\text{Tr}({\rho }_0{c}_{j_1}{c}_{j_2}\cdots c_{j_{2n}})=i^n\text{Pf}(\bar{M}{|}_{j_1{j}_2\cdots j_{2n}})$$

$$\omega ({\rho }_A)=\frac{1}{2^{N_A}}\exp \left(-\frac{i}{2}{\theta }_A^T{\bar{M}}_A{\theta }_A\right),$$

$$\text{Tr}({\rho }_A{X}_A)=\text{Tr}({\rho }_0{X}_A).$$

$$W_A{\bar{M}}_A{W}_A^T=\left(\begin{array}{ccccc}0& {\eta }_1& & & \\ -{\eta }_1& 0& & & \\ & & \ddots & & \\ & & & 0& {\eta }_N\\ & & & -{\eta }_N& 0\end{array}\right),$$

$${\rho }_A=\prod _{j=1}^{N_A}\left(\frac{1}{2}-i\frac{{\eta }_j}{2}{c}_{A,2j-1}^{\prime }{c}_{A,2j}^{\prime }\right),$$

Full text

Introduction

Quantum phase transition is one of the pillars underpinning condensed matter physics [1]. Conventional wisdom states that different quantum phases are generally discriminated in terms of the symmetry carried by the ground state or other features that have an underlying topological interpretation [2]. The former is described by local order parameters associated with the symmetries and the latter by topological orders, which are nowadays classified into intrinsic and symmetry-protected ones [3, 4, 2, 5, 6, 7, 8, 9]. In both cases, a continuous transition between different phases is mediated by a spectrally gapless critical point, in the vicinity of which thermodynamic quantities exhibit scaling behaviours classified into universality classes [10].

Modern understanding of quantum phase transition has been significantly enriched by incorporating the concept of entanglement [11, 12, 13, 14, 15, 16, 7, 8, 17]. Quantum phases are determined after all by the way how different particles or different parties in the system are mutually related. In this sense, it is natural to expect that entanglement would bear the fingerprint of the quantum phase. This perspective is especially powerful in the study of topological orders, which are a purely quantum effect. For example, states with an intrinsic topological order have a long-range entanglement and a nonzero topological entanglement entropy [14, 15, 7]. Symmetry-protected topological orders are signified by a degenerate entanglement spectrum [8]. Topological quantum phase transition would then be thought of as a rearrangement of the pattern of entanglement.

While it is a common practice to study macroscopic bipartite entanglements in topological phases, they apparently reveal only a partial aspect of many-body entanglement and, on the practical side, are hardly accessible in experiments. It is thus worthwhile to carry out a more detailed inspection of the many-body entanglement for a deeper understanding of topological phases. In particular, when it comes to the aspect of phase transition, local entanglement may be enough to gain information on the critical singularities, as is suggested by earlier works on symmetry-breaking quantum phase transitions in Heisenberg spin chains [11, 12]. If then, an interesting question is how its singular nature differs from that of the symmetry-breaking transitions. Besides, from the viewpoint that different quantum phases are imprinted in different patterns of entanglement, to examine many-body entanglement in topological models is an interesting problem in its own right.

The aim of this paper is to investigate the quantum phase transition of one-dimensional topological models in terms of the two-site entanglements, namely, the concurrences, in the ground state [18]. As a prototypical model, we consider the Su-Schrieffer-Heeger (SSH) model on a one-dimensional lattice, as shown in Figure 1 [19]. The system has $N$ unit cells, each consisting of two sites $\{\text{a}_{n},\text{b}_{n}\}$. The model Hamiltonian is given by

[TABLE]

where $\{a_{n},b_{n}\}$ denote the fermion operators for the $n$-th unit cell and the periodic boundary condition $\{a_{N+1},b_{N+1}\}=\{a_{1},b_{1}\}$ is taken. We take the hopping rates

[TABLE]

to have a single control parameter $\lambda\in[-1,1]$. This model has been well studied in the context of the band topology [5, 9]. It is in a topological phase for $\lambda>0$ and in a trivial phase for $\lambda<0$ (this distinction of the phase, of course, depends on the choice of the unit cell).

We first derive an analytic formula for the concurrence between any two sites and show that in the thermodynamic limit, the first derivative of the concurrence between adjacent sites with respect to $\lambda$ diverges logarithmically at the critical point $\lambda_{0}=0$, the exact form of which is also derived. This result is similar to the case of the symmetry-breaking quantum phase transition in the Heisenberg spin chain [11, 12]. However, due to the topological origin, there exists an interesting difference: for finite even $N$, the concurrence is discontinuous at $\lambda=\lambda_{0}$ with a gap inversely proportional to $N$, while it remains analytic for odd $N$. This feature contrasts with the case of symmetry-breaking quantum phase transitions wherein the nonanalyticity appears only in the thermodynamic limit. We provide a geometric interpretation of that, directly relating the singularity with the change of the band topology. This phenomenon has a universal nature in one-dimensional topological phase transitions of noninteracting fermions. As an example, we show that for the Kitaev chain [4], the local electron density is an analogous quantity exhibiting the identical nonanalyticity at the critical point: its first derivative, i.e., the local compressibility, diverges logarithmically in the thermodynamic limit and it is discontinuous for finite even system sizes.

Another finding we present in this paper is that there exists a different type of phase transition in the SSH model whereby the many-body entanglement of the ground state undergoes a sudden change in the following sense. As a means to characterise the many-body entanglement, we represent the pairwise pattern of all the concurrences as a (simple) graph, where each edge means the existence of entanglement, i.e., a nonzero concurrence, between the two vertices (sites). We call this graph an “entangled graph”, following Ref. [20]. If two many-body states have different entangled graphs, we will regard them as having different classes of many-body entanglement and hence belonging to different phases. This kind of characterisation of many-body entanglement has a relatively long tradition [21, 22, 20]. In the present work, our particular motivation is coming from the fact that at two extreme phases $\lambda=\pm 1$, the ground state is dimerised in such a way that either two sites in every unit cell form a singlet ($\lambda=-1$) or every adjacent pair of sites across unit cells form a singlet ($\lambda=+1$), the entanglement of which can be naturally represented by the entangled graphs shown in Figure 2. The problem now is then to figure out how the entangled graph for $\lambda\in(-1,1)$ interpolates between the two. It turns out that at $\lambda_{\pm}\simeq\pm 0.138$, the entangled graph changes suddenly along with a nonanalyticity of the concurrence. This singularity is reminiscent of the phenomenon called the entanglement sudden death [23]. We emphasise, however, that this transition should not be confused with conventional quantum phase transitions because it has nothing to do with a nonanalyticity of the ground-state wavefunction itself: the nonanalytic behaviour comes from the way the entanglement is defined and quantified. A further remark on this transition will be made later. In what follows, “quantum phase transition” will refer only to the transition at $\lambda=\lambda_{0}$.

Results

Reduced density matrices in the Su-Schrieffer-Heeger model

It is convenient to switch to the momentum space, in which the Hamiltonian (1) takes the Bogoliubov-de Gennes form:

[TABLE]

where $\phi_{k}=(a_{k}\leavevmode\nobreak\ b_{k})^{T}$, $\vec{\sigma}=\sigma_{x}\hat{x}+\sigma_{y}\hat{y}+\sigma_{z}\hat{z}$ is the vector composed of Pauli matrices, and $\vec{h}(k)$ is given by

[TABLE]

Our aim is to obtain the concurrences between all pairs of lattice sites, which we denote by $\mathcal{C}(\text{a}_{n},\text{b}_{m})$, $\mathcal{C}(\text{a}_{n},\text{a}_{m})$, and $\mathcal{C}(\text{b}_{n},\text{b}_{m})$. For this, we need to obtain the corresponding reduced density matrices $\rho_{\text{a}_{n}\text{b}_{m}}$, $\rho_{\text{a}_{n}\text{a}_{m}}$, and $\rho_{\text{b}_{n}\text{b}_{m}}$. To this end, we follow the method presented in Ref. [8], which we also recast in Methods. The first step needed is to spectrally flatten the Hamiltonian, which means that all positive eigenvalues of the Hamiltonian are replaced by +1, while all negative ones by -1. A crucial observation is that the spectral flattening actually amounts to replacing $\vec{h}(k)$ with a normalised one $\hat{h}(k)=\hat{h}_{x}(k)\hat{x}+\hat{h}_{y}(k)\hat{y}=\vec{h}(k)/|\vec{h}(k)|$. We can write

[TABLE]

where it is understood that the lefthand side is the conventional Majorana representation of the righthand side with $\bar{M}$ being a real skew-symmetric matrix. The next step is to project $c$ and $\bar{M}$ onto the subspace under consideration. For example, in order to obtain $\rho_{\text{a}_{n}\text{b}_{m}}$, one takes

[TABLE]

where

[TABLE]

The quantity $\eta_{\text{a}_{n}\text{b}_{m}}$ will play a central role in our analysis. Once the $4\times 4$ matrix $\bar{M}_{\text{a}_{n}\text{b}_{m}}$ is obtained from Eq. (6), one can follow the prescription in Methods to have

[TABLE]

in the basis $\{|0\rangle,b_{m}^{\dagger}|0\rangle,a_{n}^{\dagger}|0\rangle,a_{n}^{\dagger}b_{m}^{\dagger}|0\rangle\}$, the subscript in $\eta_{\text{a}_{n}\text{b}_{m}}$ being omitted for brevity.

Phase diagram of many-body entanglement

The concurrence can be directly calculated from the reduced density matrix (9) as

[TABLE]

which is nonzero only for $\eta_{\text{a}_{n}\text{b}_{m}}>\sqrt{2}-1$ [18]. It turns out that $\eta_{\text{a}_{n}\text{a}_{m}}=\eta_{\text{b}_{n}\text{b}_{m}}=0$ for all $\{n,m\}$, which implies that $\rho_{\text{a}_{n}\text{a}_{m}}$ and $\rho_{\text{b}_{n}\text{b}_{m}}$ are fully mixed states and thus $\mathcal{C}_{\text{a}_{n}\text{a}_{m}}=\mathcal{C}_{\text{b}_{n}\text{b}_{m}}=0$. Furthermore, $\mathcal{C}(\text{a}_{n},\text{b}_{m})$ also all vanish except for $\mathcal{C}(\text{a}_{n},\text{b}_{n})$ and $\mathcal{C}(\text{b}_{n},\text{a}_{n+1})$. As a result, there are only three possible entangled graphs associated with the ground state, as shown in Figure 2. Hereafter, for a notational simplicity, we let $\eta_{1}\equiv\eta_{\text{a}_{n}\text{b}_{n}}$, $\mathcal{C}_{1}\equiv\mathcal{C}_{\text{a}_{n}\text{b}_{n}}$, $\eta_{2}\equiv\eta_{\text{b}_{n}\text{a}_{n+1}}$, and $\mathcal{C}_{2}\equiv\mathcal{C}_{\text{b}_{n}\text{a}_{n+1}}$. The two concurrences are related as

[TABLE]

In Figure 2, we plot $\mathcal{C}_{1}(\lambda)$, $\mathcal{C}_{2}(\lambda)$, and $d\mathcal{C}_{2}(\lambda)/d\lambda$ with respect to $\lambda$. One can distinguish four different phases $\mathbf{P}_{\{0,1\}}$ and $\mathbf{Q}_{\{0,1\}}$. The subscript represents the $Z_{2}$ index for the band topology with 0 (1) representing the trivial (topological) phase. In the trivial (topological) phase, $\mathcal{C}_{1}>\mathcal{C}_{2}$ ($\mathcal{C}_{1}<\mathcal{C}_{2}$). In the phase $\mathbf{Q}$, the entangled graph is connected, while in the phase $\mathbf{P}$, it is disconnected. The phase transition of the band topology occurs at $\lambda_{0}=0$ and that of the entangled graph occurs at $\lambda_{\pm}\simeq\pm 0.138$.

Quantum phase transition at $\lambda=\lambda_{0}$

The nonanalyticity of the concurrence here can be ascribed to the sudden change of the band topology. The vector $\vec{h}(k)$ in Eq. (3) traverses a circle on the $x-y$ plane as $k$ sweeps over the Brillouin zone $0\leq k<2\pi$, as in Figure 3. The $Z_{2}$ topological index is then determined by whether the circle encloses the origin or not [17]. In order to make this topological distinction more pronounced, one can take the unit vector $\hat{h}(k)$, as in Figure 3. In the trivial phase ($\lambda<0$), $\hat{h}(k)$ wanders on the half circle in $x>0$. To close the loop, it traverses the same path twice. In the topological phase ($t_{1}<t_{2}$), on the other hand, $\hat{h}(k)$ traverses the full circle once.

Such a topological difference is in fact captured by the quantity $\eta_{1}$:

[TABLE]

At the phase transition, the trajectory of $\hat{h}(k)$ should change its shape from one to the other in Figure 3. In the thermodynamic limit wherein the trajectory becomes continuous, such a change cannot be made continuously. In view of Eq. (12), this discontinuity should be reflected as a sudden jump of $\eta_{1}$ in the vicinity of the critical point, rendering the derivative of $\eta_{1}$ with respect to $\lambda$, and hence that of $\mathcal{C}_{1}$, diverging at the critical point.

This nonanalyticity, as shown in Figure 2, is reminiscent of earlier results for the cases of symmetry-breaking quantum phase transitions in spin chains [11, 12]. However, a stark difference is found for finite systems due to the different origins of the singularity. In the present case, the behaviour of $\eta_{1}(\lambda)$ depends on the parity of $N$, as shown in Figure 4. For even $N$, $\eta_{1}(\lambda)$ is discontinuous at $\lambda=0$ even though the system size is finite. This can be understood by re-examining Figure 3 and Eq. (12). One can realise that for even $N$, as the trajectory of $\vec{h}(k)$ is made of $N$ equally spaced points on the circle, the difference between $\eta_{1}(0_{+})$ and $\eta_{1}(0_{-})$ is made solely by the contribution of $\hat{h}_{x}(k=\pi)$, which is either $1$ or $-1$ depending on the phase. We thus find that

[TABLE]

For odd $N$, due to the exclusion of $k=\pi$, the discrete nature of $\eta_{1}(\lambda)$ disappears.

We perform additional calculations taking more realistic situations into account. Figure 4(d) shows $\mathcal{C}_{2}(\lambda)$ in the presence of disorder, which is introduced by adding $\sum_{n}(\epsilon_{n}^{a}a_{n}^{\dagger}a_{n}+\epsilon_{n}^{b}b_{n}^{\dagger}b_{n})$ to the Hamiltonian with $\epsilon_{n}^{\{a,b\}}$ being taken randomly and uniformly from the interval $[-0.1,0.1]$. It turns out that small disorder does not significantly alter the essential features of the entanglement including the distinction between the cases of even and odd $N$. Figure 4(e) shows $d\mathcal{C}_{2}/d\lambda$ obtained by taking the two sites at the centre of the chain in the case of the open boundary condition. In this case, the peak of $d\mathcal{C}_{2}/d\lambda$ does not coincide with $\lambda=0$ because the symmetry with respect to $\lambda\leftrightarrow-\lambda$, as in Eq. (11), is broken down. Also, the difference according to the parity of $N$ is absent because its origin—the inclusion or exclusion of $k=\pi$ in the Brillouin zone—now loses its meaning. Here, the chemical potential is placed in the gap either below or above the two zero-energy states that appear due to the open boundary condition. Either case produces the same results. As is expected, as $N$ increases, the singularity is more pronounced and the peak approaches $\lambda=0$.

It is worthwhile to interpret the nonanalyticity in terms of the conventional language of the phase transition. In fact, the derivative of the free energy with respect to $\lambda$ picks up $\eta$:

[TABLE]

where $F(T,\lambda)=-k_{B}T\ln\text{Tr}[e^{-H(\lambda)/(k_{B}T)}]$. As the phase transition is second order, the derivative of $\eta$, and hence the derivative of the concurrence, diverges at the critical point. In the thermodynamic limit, $\eta_{1}(\lambda)$ diverges logarithmically for $|\lambda|\ll 1$ as

[TABLE]

which can be derived by using the property of the elliptic integral: $-\int_{0}^{\pi/2}(\lambda^{2}+\sin^{2}\theta)^{-1/2}d\theta\simeq\log(\lambda/4)$.

The above arguments remain valid for arbitrary one-dimensional topological phase transitions governed by a two-band Hamiltonian as in Eq. (3). In the thermodynamic limit, there should be a local quantity corresponding to $\eta_{1}$, which becomes nonanalytic at the critical point, reflecting the change of the winding number. For finite systems, it is discontinuous if the band-touching point belongs to the (discrete) Brillouin zone, which would be naturally related to a commensurability problem. For example, in the Kitaev chain [4], the reduced density matrix for a single site, say, site 1, is obtained as

[TABLE]

where $a_{1}$ is the fermion operator for site 1. Note that the coefficient is identical to Eq. (12) except for the orientation of the vector being different. In this case, the local electron density $\langle a_{1}^{\dagger}a_{1}\rangle$ as a function of the chemical potential is discontinuous at the critical point for finite even $N$ and the local compressibility diverges logarithmically as Eq. (15) in the thermodynamic limit. The latter behaviour has been addressed in a recent work [24].

Phase transition of entanglement at $\lambda=\lambda_{\pm}$

The nonanalyticity of the concurrence here is originated from the property of entanglement, which is defined as “not being separable”. In the Hilbert space of bipartite density operators, the set of all separable states forms a compact convex set [25, 26]. As a result, if one traces a continuous path in the Hilbert space from an entangled to a separable state, the entanglement suddenly disappears when one crosses the hyperplane separating the sets of separable and entangled states. When this occurs in a dynamical problem, the phenomenon is called an entanglement sudden death (or sudden birth in the opposite way) [23]. The entanglement sudden death is observed occasionally when a state evolves in a dissipative environment. However, it is rare to see an analogous phenomenon in the course of a quantum phase transition (to our best efforts, we could not find out a prior example). In view of Eq. (10), the existence of this phase transition should be robust against a small perturbation to the system. As the system remains gapped, one can transform it into a dynamical problem by considering an adiabatic evolution varying $\lambda$, making the link to the entanglement sudden death clearer.

Near $\lambda=\lambda_{+}$, $\mathcal{C}_{1}(\lambda)$ changes linearly for $\lambda<\lambda_{+}$ and vanished for $\lambda>\lambda_{+}$. In the thermodynamic limit, the former behaviour is given by

[TABLE]

For finite $N$, $\lambda_{+}$ changes with $N$ as shown in Figure 5. For $N=2$ or $N=4$, $\mathcal{C}_{1}>0$ for $\lambda<0$ and $\mathcal{C}_{1}=0$ for $\lambda>0$, and hence we cannot find $\lambda_{+}$.

As we have addressed in Introduction, this singularity is not originated from the nonanalyticity of the ground-state wave function and hence the transition is not a quantum phase transition in the conventional sense. Nonetheless, the transition may find meaning in the context of quantum information theory. For example, supposing two parties, say, Alice and Bob, possess each of the two sites in a unit cell, the usefulness of the state as a resource for a quantum information processing depends on the phase [27]. Here, instead of considering entanglement of the fermion occupation number, which is somewhat impractical, one can make the scenario more practical by turning the SSH model into an equivalent $s=1/2$ spin-chain model with nearest-neighbour interaction. That is, regarding $a_{n}$ and $b_{n}$ in Eq. (1) as $s=1/2$ spin lowering operators, in the subspace with the half-filling $\sum_{n}(a_{n}^{\dagger}a_{n}+b_{n}^{\dagger}b_{n})=N$, one can realise that the equivalent spin-chain Hamiltonian is

[TABLE]

Note that the last term introduces a gauge field for even $N$, which would be responsible for the different behaviour at $\lambda=\lambda_{0}$ for even and odd $N$. A deeper operational meaning of the entanglement phase transition at $\lambda=\lambda_{\pm}$ is open to question due to the lack of a relevant quantum informational protocol in the multipartite setting.

Discussion

We have fully characterised two-site entanglements in the ground state of the SSH model and obtained the phase diagram of the entanglement. It was found out that there are two kinds of singularities in the entanglement: one due to the topological quantum phase transition and the other one due to the entanglement sudden death.

Several remarks are in order. (i) The presented singularities are expected to be verifiable in optical lattices. For example, one could prepare a fermionic Mott insulator in a one-dimensional optical lattice and adiabatically change the potential shape by superimposing an additional lattice potential so that the number of sites is twice the number of atoms in the end [28, 29]. More detailed and rigorous analysis of its feasibility is left as a future work. (ii) As can be seen in Figure 2, one can detect the topological order of the SSH model simply by comparing local quantities $\mathcal{C}_{1}$ and $\mathcal{C}_{2}$, instead of referring to the entanglement spectrum. This might greatly facilitate the experimental detection of the topological phase. However, the validity of this idea in the presence of disorder or interaction and its generalisation to other models are unclear. (iii) While our way of characterising the many-body entanglement seems reasonable for the SSH model, one can employ a different characterisation to see a different aspect of the many-body entanglement. The phases in Figure 2 would then be divided into subphases, enriching the phase diagram.

Methods

Reduced density matrix for noninteracting gapped fermion systems

We recast the formulae presented in Refs. [4, 30, 8] to derive a reduced density matrix of the ground state of a gapped quadratic fermion Hamiltonian. Consider a system described by $N$ fermion annihilation (creation) operators $a_{j}$ ($a_{j}^{\dagger}$) with $1\leq j\leq N$. It is convenient to introduce Majorana operators $c_{l}$ with $1\leq l\leq 2N$ such that $a_{j}=\frac{1}{2}(c_{2j-1}+ic_{2j})$. Let $c$ be the vector with $2N$ elements $c_{l}$. Hereafter, we will use a similar convention for other operators as well. The Hamiltonian can be generally written as

[TABLE]

where $M$ is a $2N\times 2N$ real skew-symmetric matrix with $M_{jk}=-M_{kj}$. The matrix $M$ is block diagonalised as

[TABLE]

by a $2N\times 2N$ real orthogonal matrix $W$ with $W^{T}W=WW^{T}=I$, where $\epsilon_{j}>0$. The Hamiltonian is then written as

[TABLE]

where $c^{\prime}=Wc$ and $a^{\prime}_{j}=\frac{1}{2}(c^{\prime}_{2j-1}+ic^{\prime}_{2j})$. The ground state in the density matrix form is thus

[TABLE]

For later use, let us define a spectral flattening of the Hamiltonian $H$ as

[TABLE]

which amounts to replacing all $\epsilon_{j}$ with 1. We also define a Grassmann representation $\omega(X)$ for a polynomial $X$ of Majorana operators, which is done by replacing all Majorana operators in $X$ with Grassmann variables. For example, a Grassmann representation of $\rho_{0}$ is obtained by replacing $c_{l}$ and $c^{\prime}_{l}$ with Grassmann variables $\theta_{l}$ and $\theta^{\prime}_{l}=\sum_{m}W_{lm}\theta_{m}$, respectively, as

[TABLE]

Note that $\bar{M}$ is obtained from the spectrally flattened Hamiltonian (23). A state that has a Gaussian form of a Grassmann representation as in Eq. (24) is called a Gaussian state. It can be checked that $\bar{M}$ is in fact a two point correlation Matrix

[TABLE]

and higher order correlations are given by

[TABLE]

for $j_{1}<j_{2}<\cdots<j_{2n}$, where $\bar{M}|_{j_{1}j_{2}\cdots j_{2n}}$ is a $2n\times 2n$ submatrix of $\bar{M}$ restricted to the designated indices and $\text{Pf}(\cdot)$ is the Pfaffian. All odd-order correlations vanish.

Now let us divide the system into two parts $A$ and $B$. $A$ is described by $N_{A}$ pairs of fermion operators and $B$ by the rest. Our aim is to calculate the reduced density matrix $\rho_{A}=\text{Tr}_{B}(\rho_{0})$. One can check that if the Grassmann representation of $\rho_{A}$ is given by

[TABLE]

where $\theta_{A}$ and $\bar{M}_{A}$ are the corresponding submatrices restricted to the indices for the subsystem $A$, then for an arbitrary polynomial $X_{A}$ of the Majorana operators supported on A,

[TABLE]

This implies that Eq. (27) is indeed correct. As $\bar{M}_{A}$ is also skew-symmetric, it is block diagonalised as

[TABLE]

where $W_{A}$ is real orthogonal and $0\leq\eta_{j}\leq 1$. Taking the submatrix $c_{A}$ of $c$ and letting $c^{\prime}_{A}=W_{A}c_{A}$, we finally obtain

[TABLE]

which has the Grassmann representation as in Eq. (27).

Acknowledgements

We thank Gil Young Cho, Isaac Kim, Ki-Seok Kim, and Heung-Sun Sim for fruitful discussions. This research was supported (in part) by the R&D Convergence Program of NST (National Research Council of Science and Technology) of Republic of Korea (Grant No. CAP-15-08-KRISS). K.W.K. acknowledges support from the “Overseas Research Program for Young Scientists” by Korea Institute for Advanced Study (KIAS).

Author contributions statement

J.C. conceived the preliminary idea and then J.C. and K.W.K. equally contributed to the analysis. J.C. wrote the manuscript.

Additional Information

Competing financial interests: The authors declare no competing financial interests.