
2
FIG. 1. Two-dimensional tight-binding system on a square
lattice. The geometry of the system is a rectangle, where the
lengths in the xdirection and in the ydirection are Lxand
Ly, respectively. The black circles express the unit cells at the
lattice point r= (nx, ny) (nx= 1,...,Lx, ny= 1,...,Ly).
This paper is organized as follows: In Sec. II, we in-
troduce our 2D non-Hermitian tight-binding model and
propose two classes of systems where we can apply the
non-Bloch band theory. Then, we can get the condition
for the generalized Brillouin zone. In Sec. III, we ana-
lyze the three 2D non-Hermitian models in terms of the
non-Bloch band theory. In Sec. IV, we discuss the 2D
non-Bloch band theory proposed in this work from the
viewpoint of the formation of standing wave and com-
ment on difficulty to construct the non-Bloch band the-
ory in general 2D non-Hermitian systems. Finally, in
Sec. V, we summarize our result.
II. MODEL
First of all, we introduce a 2D non-Hermitian tight-
binding system. Throughout this paper, the system lies
on a square lattice, and the geometry of the system is a
rectangle, as shown in Fig. 1. The real-space Hamiltonian
of our system is given by
H=X
r
q
X
µ,ν=1 Nx
X
i=−Nx
t(x)
i,µν c†
r+iˆ
x,µcr,ν
+
Ny
X
i=−Ny
t(y)
i,µν c†
r+iˆ
y,µcr,ν
,(1)
where c†
r,µ is a creation operator, r=
(nx, ny) (nx= 1, . . . , Lx, ny= 1, . . . , Ly) is a lat-
tice point, qis the number of internal degrees of freedom
in a unit cell, and ˆ
xand ˆ
yare unit vectors in the x
direction and in the ydirection, respectively. The parti-
cles hop up to the Nxth and Nyth nearest-neighbor unit
cells along the xand ydirections, respectively. We note
that when either t(x)
i,νµ 6=t(x)
−i,µν ∗
or t(y)
i,νµ 6=t(y)
−i,µν ∗
is satisfied, the system becomes non-Hermitian. Now,
the real-space eigen-equation is written as
H|ψi=E|ψi,(2)
where the eigenstates are given by
|ψi=...,ψ(nx,ny),1, . . . , ψ(nx,ny,1),q, . . . T.(3)
Next, we describe a way to construct the non-Bloch
band theory in our system. In the case of q= 1, it is
straightforward to construct the non-Bloch band theory
in the system. This is because the system can be de-
composed into two 1D non-Hermitian models. In Ap-
pendix A, we derive the condition for the generalized
Brillouin zone in this case. On the other hand, in the
case of q≥2, we find that when the system satisfies
some conditions, one can calculate the generalized Bril-
louin zone in terms of the conventional non-Bloch band
theory. In the following, we explain the two cases where
the description of the non-Bloch bands is possible.
A. Case A
In a 1D non-Hermitian system, when the energy of
the eigenstate with the Bloch wavenumber +kand that
of the eigenstate with −kare degenerate, the non-
Hermitian skin effect is suppressed because kbecomes
real [8, 66, 67]. The degeneracy leading to the suppres-
sion of the non-Hermitian skin effect comes from symme-
tries, e.g., a global symmetry, such as a PT symmetry,
and a crystalline symmetry, such as a mirror symmetry.
Importantly, we here propose that such suppression is re-
alized also in a 2D non-Hermitian system. For example,
when the Bloch Hamiltonian H(kx, ky) satisfies
U−1HT(kx,−ky)U=H(kx, ky),(4)
where Uis a unitary matrix satisfying U2= 1, we can
show that the non-Hermitian skin effect in the ydirection
is suppressed. In fact, when Eq. (4) is satisfied, the wave
with the Bloch wavevector (kx, ky) is reflected to that
with (kx,−ky) at the boundary of the system parallel to
the xdirection. Then, the condition for the formation
of standing wave means that the imaginary parts of the
wavevector for the incident and reflected waves are equal.
Therefore, we can obtain Im (ky) = Im (−ky), meaning
that kyis real, and the non-Hermitian skin effect in the
ydirection is suppressed.
In this case, the bulk eigenstates extend over the ydi-
rection because the Bloch wavenumber kybecomes real.
Hence, in the limit of a large system size, the asymptotic
behavior of the bulk eigenstates in the cylinder geometry,
where a periodic boundary condition in the ydirection
is imposed to the system, matches that in the full-open
geometry. Therefore, the energy bands in the cylinder
geometry asymptotically reproduces the energy levels in
a finite open plane. Importantly, since the cylinder ge-
ometry can be regarded as a pseudo-1D system, one can
calculate the energy bands of this system in terms of the
conventional non-Bloch band theory. Thus, we can con-
struct the non-Bloch band theory in the system.
Based on the above concept, we establish the non-
Bloch band theory in the 2D non-Hermitian systems
without the non-Hermitian skin effect in the ydirection.
With Eq. (4), by assuming that the Bloch wavenumber