Maximal distant entanglement in Kitaev tube

We study the Kitaev model on a finite-size square lattice with periodic boundary conditions in one direction and open boundary conditions in the other. Based on the fact that the Majorana representation of Kitaev model is equivalent to a brick wall model under the condition t = Δ = μ, this system is shown to support perfect Majorana bound states which is in strong localization limit. By introducing edge-mode fermionic operator and pseudo-spin representation, we find that such edge modes are always associated with maximal entanglement between two edges of the tube, which is independent of the size of the system.

We study the Kitaev model on a finite-size square lattice with periodic boundary conditions in one direction and open boundary conditions in the other. Based on the fact that the Majorana representation of Kitaev model is equivalent to a brick wall model under the condition t = Δ = μ, this system is shown to support perfect Majorana bound states which is in strong localization limit. By introducing edge-mode fermionic operator and pseudo-spin representation, we find that such edge modes are always associated with maximal entanglement between two edges of the tube, which is independent of the size of the system.
Topological materials have become the focus of intense research in the last years 1-4 , since they not only exhibit new physical phenomena with potential technological applications, but also provide a fertile ground for the discovery of fermionic particles and phenomena predicted in high-energy physics, including Majorana 5-10 , Dirac [11][12][13][14][15][16][17] and Weyl fermions [18][19][20][21][22][23][24][25][26] . These concepts relate to Majorana edge modes. A gapful phase can be topologically non-trivial, commonly referred to as topological insulators and superconductors, the band structure of which is characterized by nontrivial topology. The number of Majorana edge modes is determined by bulk topological invariant. In general, edge states are the eigenstates of Hamiltonian that are exponentially localized at the boundary of the system. A particularly important concept is the bulk-edge correspondence, which links the nontrivial topological invariant in the bulk to the localized edge modes. On the other hand, Majorana edge modes have been actively pursued in condensed matter physics 27-33 since spatially separated Majorana fermions lead to degenerate ground states, which encode qubits immune to local dechoerence 34 . There have been theoretical proposals for detecting Majorana fermions in 2D semiconductor heterostructures 35,36 , topological insulator-superconductor proximity 5,37-40 , 1D spin-orbit-coupled quantum wires 6,9,[41][42][43][44][45][46][47][48] and cold atom systems 33,[49][50][51][52][53] . Experimentally, it is claimed that indirect signatures of Majorana fermions in topological superconductors have been observed 7,8,[54][55][56][57][58][59][60][61][62] . So far, the thoretically predicted Majorana bound state in literatures requires the system in thermodynamic limit. An interesting question is whether there exists the Majorana bound state in a small sized system, or the topological feature is a prerequisite for Majorana bound state. The existence of such a type of edge mode would indicate that the bulk topology is not necessary to the spatially separated Majorana fermions and may provide an alternative way to detect and utilize Majorana fermions.
In this paper, we study the Majorana edge modes in the Kitaev model on a square lattice based on analytical solutions. In contrast to previous studies based on open boundary conditions in two directions, we focus on a finite-length cylindrical lattice. We show that the Majorana representation of Kitaev model is related to a brick wall model, based on which this model in a finite-length cylindrical geometry supports the perfect Majorana bound states under the condition μ = Δ = t . The perfect Majorana bound state is in the strong localization limit. This Majorana zero mode has two notable features: (i) The edge-mode states exhibit maximal entanglement between the two edges of the cylinder; (ii) By introducing edge-mode pesudospin operators, we find that the edge mode relates to a conserved observable. Remarkably, the expectation values of two types of pseudospins for eigenstates indicate the coexistence of both bosonic and fermionic excitations. And the eigenstates also possess maximal entanglement about the bosonic and fermionic modes. These results provide a way to detect the Majorana bound states in p-wave superconductors.

Model
We consider the Kitaev model on a square lattice which is employed to depict 2D p-wave superconductors. The Hamiltonian of the tight-binding model on a square lattice takes the following forma where r is the coordinates of lattice sites and c r is the fermion annihilation operators at site r. Vectors a = ai, aj, are the lattice vectors in the x and y directions with unitary vectors i and j. The hopping between (pair operator of) neighboring sites is described by the hopping amplitude t (the real order parameter Δ). The last term gives the chemical potential. Imposing boundary conditions on both directions, the Hamiltonian can be exactly diagonalized. The Kitaev model on a honeycomb lattice and chain provides well-known examples of systems with such a bulk-boundary correspondence [63][64][65][66][67][68][69] . It is well known that a sufficient long chain has Majorana modes at its two ends 70 . A number of experimental realizations of such models have found evidence for such Majorana modes 7,54,56,57,71 . In contrast to previous studies based on system in thermodynamic limit, we focus on the Kitaev model on a finite lattice system. This is motivated by the desire to get a clear physical picture of the egde mode via the investigation of a small system. We first study the present model from the description in terms of Majorana fermions.
We introduce Majorana fermion operators r r r r r r which satisfy the relations

r a rr
It represents a dimerized brick wall lattice (or honeycomb lattice) with extra hopping term b r a r+a .

Majorana Edge Modes
Let us consider a simple case to show that Majorana modes can appear on some edges. Taking t = Δ = μ the Hamiltonian reduces to

Consider the Fourier transformations of Majorana operators
which obeys i.e., H BW has been block diagonalized. We note that for K = 0, we have In general, all the eigenstates of H BW can be classified into two groups {|Φ + 〉} and {|Φ − 〉}, which are constructed as the forms We are interested in the feature of edge-mode operator d M . It is easy to check that  In short, a zero-energy mode is characterized by a conventional fermion operator, which is also referred as edge-mode operator. Any standard fermion operator has its own vacuum and particle states, or hole and particle states. We have shown that the corresponding hole and particle states for the edge-mode operators d M and † d M are both EPR pair states in the spin representation. It reveals the non-locality of edge mode, through the particle and hole states are not eigenstates of the system.
For the eigenstate, the long-range correlation still exists. In the section Method, we show that the eigenstate |Φ ± 〉 can be regarded as entangled states between boson and fermion. It is expected that such a framework can be applied to more general cases.

Summary
In this paper we have studied the edge modes of a finite size Kitaev model on a square lattice. The advantage of studying the finite system is that the obtained result can be demonstrated in synthetic lattice system. We studied the Majorana edge modes for the Kitaev model in a cylindrical geometry. The Majorana representation of the Hamiltonian turns out to be equivalent to a brick wall model under some conditions. The analytical solutions show that there exist perfect Majorana edge modes, which are in the strong localization limit. We provide a new way to analyze the excitation mechanisms in the framework of pseudospins for the edge modes. These modes, in contrast to the modes in Kitaev chain, can appear in small finite systems. This may provide a new venue for observing Majorana fermions in experiments.

Method
Pseudospin description. To get insight into the feature of the edge-mode related eigenstates in such a cylindrical Kitaev model, we introduce two types of pseudospin operators , [ , ] , (31) and where α β γ = x y z , , , , . Based on these relations the combination operator τ = + α α α J s obeys the standard angular momentum relation We note that J x is a conserved observable since = J H [ , ] 0 x BW . On the other hand, both operators J x and H BW are also invariant under a local particle-hole transformation , which is defined as , tells us operators J x ,  and H BW share a common eigen vectors |Φ ± 〉. In fact, one pseudospin can be transformed to the other (and vice versa) by applying the transformation . Direct derivation shows that which result in ± ± ± ± s 3/8 Together with x we find that state |Φ ± 〉 does not have definite values of s and τ. Unlike a standard spin operator which has its own vector space, operators {s α } and {τ β } share a common vector space. Actually, for two edge sub-system, there are total four possible states which can be written down as