Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Non-Hermitian Topology in Hermitian Topological Matter

Shu Hamanaka hamanaka.shu.45p@st.kyoto-u.ac.jp Department of Physics, Kyoto University, Kyoto 606-8502, Japan Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland    Tsuneya Yoshida yoshida.tsuneya.2z@kyoto-u.ac.jp Department of Physics, Kyoto University, Kyoto 606-8502, Japan Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland    Kohei Kawabata kawabata@issp.u-tokyo.ac.jp Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan
(May 16, 2024)
Abstract

Non-Hermiticity leads to distinctive topological phenomena absent in Hermitian systems. However, connection between such intrinsic non-Hermitian topology and Hermitian topology has remained largely elusive. Here, considering the bulk and boundary as an environment and system, we demonstrate that anomalous boundary states in Hermitian topological insulators exhibit non-Hermitian topology. We study the self-energy capturing the particle exchange between the bulk and boundary, and demonstrate that it detects Hermitian topology in the bulk and induces non-Hermitian topology at the boundary. As an illustrative example, we show the non-Hermitian topology and concomitant skin effect inherently embedded within chiral edge states of Chern insulators. We also find the emergence of hinge states within effective non-Hermitian Hamiltonians at surfaces of three-dimensional topological insulators. Furthermore, we comprehensively classify our correspondence across all the tenfold symmetry classes of topological insulators and superconductors. Our work uncovers a hidden connection between Hermitian and non-Hermitian topology, and provides an approach to identifying non-Hermitian topology in quantum matter.

Topological phases of matter are a central topic in modern condensed matter physics [1, 2]. Gapped phases of noninteracting fermions are systematically classified by the tenfold fundamental symmetry [3], culminating in the periodic table of topological insulators and superconductors [4, *Ryu-10, 6, 7]. A hallmark of topological insulators is the bulk-boundary correspondence: the nontrivial bulk topology yields anomalous gapless states at boundaries.

Beyond the Hermitian regime, topological characterization of non-Hermitian systems has recently attracted growing interest [8, 9, *Esaki-11, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, *Kawabata-19, 22, *YSW-18-Chern, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, *Denner-23JPhysMater, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52]. Non-Hermiticity arises from the exchange of particles and energy with the environment [53, 54]. Even within closed systems, non-Hermiticity of self-energy characterizes finite-lifetime quasiparticles [55, *Papaj-19, *Shen-18QO, 58, 59, 60, 61, 62, 63, 64, 65]. Importantly, non-Hermiticity enables a unique gap structure of complex spectra—point gap—and concomitant topological phases that have no analogs in Hermitian systems [20, 29]. As the bulk-boundary correspondence in non-Hermitian systems, nontrivial point-gap topology leads to the non-Hermitian skin effect [37, 38] and anomalous boundary states [38, 39, *Denner-23JPhysMater, 44, 48]. Such intrinsic non-Hermitian topological phenomena, including the skin effect, have been realized in various experiments involving open classical and quantum systems [66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, *Bandres-18, 78, 79, 80, *Ghatak-19-skin-exp, 82, *Hofmann-19-skin-exp, 84, 85, 86, *Liang-22, 88, 89, 90, 91].

Notably, the topological classification indicates the correspondence of Hermitian topology in d𝑑ditalic_d dimensions and non-Hermitian topology in d1𝑑1d-1italic_d - 1 dimensions [29, 92, 93]. In line with this correspondence, non-Hermitian perturbations can induce point gaps for anomalous boundary states in topological insulators [94, 95, *Zhu-23, 97, 98, 99, 100], as observed in recent photonic [101] and phononic [102] experiments. However, the direct connection between Hermitian and non-Hermitian topology has remained elusive, given that non-Hermiticity is merely added as an external perturbation. Consequently, the fundamental mechanism underlying the correspondence between the Hermitian bulk and non-Hermitian boundary has been unclear.

In this Letter, we reveal non-Hermitian topology inherently embedded within Hermitian topological matter. Regarding the bulk and boundary in Hermitian topological insulators as an environment and system, respectively, we show that effective boundary Hamiltonians exhibit non-Hermitian topology (Fig. 1). As an illustrative example, we demonstrate the non-Hermitian skin effect in chiral edge states of Chern insulators. Furthermore, we systematically classify our correspondence across all the tenfold symmetry classes of topological insulators and superconductors (Table 1). Our work uncovers a hidden connection between Hermitian and non-Hermitian topology.

Refer to caption
Figure 1: Effective non-Hermitian Hamiltonian Hedge+Σsubscript𝐻edgeΣH_{\rm edge}+\Sigmaitalic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + roman_Σ at the edge of a Hermitian Hamiltonian. The self-energy ΣΣ\Sigmaroman_Σ describes the particle exchange between the bulk Hbulksubscript𝐻bulkH_{\rm bulk}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT and edge Hedgesubscript𝐻edgeH_{\rm edge}italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT, yielding non-Hermitian topology.

Non-Hermitian Hamiltonians in closed systems.—Our central idea involves conceptually dividing the bulk and boundary of a closed system into an environment and system, respectively (Fig. 1). The single-particle Hamiltonian of the entire system reads

H=(HbulkTTHedge),𝐻matrixsubscript𝐻bulk𝑇superscript𝑇subscript𝐻edgeH=\begin{pmatrix}H_{\rm bulk}&T\\ T^{{\dagger}}&H_{\rm edge}\end{pmatrix},italic_H = ( start_ARG start_ROW start_CELL italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT end_CELL start_CELL italic_T end_CELL end_ROW start_ROW start_CELL italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_CELL start_CELL italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) , (1)

where Hbulksubscript𝐻bulkH_{\rm bulk}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT (Hedgesubscript𝐻edgeH_{\rm edge}italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT) is the Hamiltonian in the bulk (at the boundary), and T𝑇Titalic_T denotes the coupling between the bulk and boundary. From the original Schrödinger equation H(|ψbulk|ψedge)T=(E+iη)(|ψbulk|ψedge)T𝐻superscriptsubscriptket𝜓bulksubscriptket𝜓edge𝑇𝐸i𝜂superscriptsubscriptket𝜓bulksubscriptket𝜓edge𝑇H\,(\ket{\psi}_{\rm bulk}\,\ket{\psi}_{\rm edge})^{T}=\left(E+\text{i}\eta% \right)(\ket{\psi}_{\rm bulk}\,\ket{\psi}_{\rm edge})^{T}italic_H ( | start_ARG italic_ψ end_ARG ⟩ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT | start_ARG italic_ψ end_ARG ⟩ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT = ( italic_E + i italic_η ) ( | start_ARG italic_ψ end_ARG ⟩ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT | start_ARG italic_ψ end_ARG ⟩ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT with an infinitesimal number η>0𝜂0\eta>0italic_η > 0 reflecting causality, we project the boundary degree of freedom and derive the effective Hamiltonian Heff(E)=Hedge+Σ(E)subscript𝐻eff𝐸subscript𝐻edgeΣ𝐸H_{\mathrm{eff}}\left(E\right)=H_{\mathrm{edge}}+\Sigma\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) = italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + roman_Σ ( italic_E ) with the self-energy [103, 104, 105]

Σ(E)T(E+iηHbulk)1T.Σ𝐸superscript𝑇superscript𝐸i𝜂subscript𝐻bulk1𝑇\displaystyle\Sigma\left(E\right)\coloneqq T^{{\dagger}}\left(E+\text{i}\eta-H% _{\rm bulk}\right)^{-1}T.roman_Σ ( italic_E ) ≔ italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T . (2)

The self-energy Σ(E)Σ𝐸\Sigma\left(E\right)roman_Σ ( italic_E ) captures the continuous exchange of particles between the bulk and boundary for given energy E𝐸E\in\mathbb{R}italic_E ∈ blackboard_R. Consequently, Σ(E)Σ𝐸\Sigma\left(E\right)roman_Σ ( italic_E ) acquires non-Hermiticity, describing finite lifetimes of quasiparticles that escape from the boundary to the bulk.

Crucially, topology of the original Hermitian Hamiltonian H𝐻Hitalic_H should leave an imprint on that of the effective non-Hermitian Hamiltonian Heff(E)subscript𝐻eff𝐸H_{\rm eff}\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ). For example, when H𝐻Hitalic_H is a quantum Hall (Chern) insulator, Heff(E)subscript𝐻eff𝐸H_{\rm eff}\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) effectively describes chiral edge states for E𝐸Eitalic_E within an energy gap. Intuitively, chirality or nonreciprocity of the anomalous boundary states should yield non-Hermitian topology, as can be seen in the celebrated Hatano-Nelson model [106, *Hatano-Nelson-97]. In this Letter, we substantiate this intuition and demonstrate that the self-energy bridges Hermitian bulk topology and non-Hermitian boundary topology. Below, for fixed E𝐸Eitalic_E (i.e., Markov approximation [62, 108]), we study non-Hermitian topology of Heff(E)subscript𝐻eff𝐸H_{\rm eff}\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) in prototypical topological insulators.

Non-Hermitian topology in the Su-Schrieffer-Heeger model.—We begin with a one-dimensional topological insulator, Su-Schrieffer-Heeger (SSH) model [109]. The Bloch Hamiltonian reads

HSSH(k)=(v+tcos(k))σx+(tsin(k))σy,subscript𝐻SSH𝑘𝑣𝑡𝑘subscript𝜎𝑥𝑡𝑘subscript𝜎𝑦\displaystyle H_{\mathrm{SSH}}\left(k\right)=\left(v+t\cos{k}\right)\sigma_{x}% +\left(t\sin{k}\right)\sigma_{y},italic_H start_POSTSUBSCRIPT roman_SSH end_POSTSUBSCRIPT ( italic_k ) = ( italic_v + italic_t roman_cos ( start_ARG italic_k end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin ( start_ARG italic_k end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , (3)

where σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s (i=x,y,z𝑖𝑥𝑦𝑧i=x,y,zitalic_i = italic_x , italic_y , italic_z) denote Pauli matrices, and v>0𝑣0v>0italic_v > 0 and t>0𝑡0t>0italic_t > 0 are the intracell and intercell hopping amplitudes, respectively.

Refer to caption
Figure 2: Non-Hermitian topology in the Su-Schrieffer-Heeger model. (a, b) Correspondence of Hermitian and non-Hermitian topology for the (a) nontrivial and (b) trivial phases. (c, d) Complex eigenvalues of the effective non-Hermitian Hamiltonian at the edge for the (c) nontrivial and (d) trivial phases.

The SSH model respects chiral symmetry σzHSSH(k)σz=HSSH(k)subscript𝜎𝑧subscript𝐻SSH𝑘subscript𝜎𝑧subscript𝐻SSH𝑘\sigma_{z}H_{\rm SSH}\left(k\right)\sigma_{z}=-H_{\rm SSH}\left(k\right)italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT roman_SSH end_POSTSUBSCRIPT ( italic_k ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = - italic_H start_POSTSUBSCRIPT roman_SSH end_POSTSUBSCRIPT ( italic_k ) and belongs to class AIII, characterized by the integer topological invariant [7]. In the topological phase (v/t<1𝑣𝑡1v/t<1italic_v / italic_t < 1), the bulk topology leads to the emergence of an edge state |ψ0x(v/t)x|x(10)Tproportional-toketsubscript𝜓0subscript𝑥tensor-productsuperscript𝑣𝑡𝑥ket𝑥superscript10𝑇\ket{\psi_{0}}\propto\sum_{x}\left(-v/t\right)^{x}\ket{x}\otimes\left(1~{}0% \right)^{T}| start_ARG italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ⟩ ∝ ∑ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( - italic_v / italic_t ) start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT | start_ARG italic_x end_ARG ⟩ ⊗ ( 1 0 ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT with zero energy E=0𝐸0E=0italic_E = 0 under the open boundary conditions. When we regard the unit cell at the left edge as a system and the remaining portion an environment [Fig. 2 (a, b)], the edge degree of freedom is coupled with the bulk and can escape into the bulk. This decaying property is quantified by the self-energy in Eq. (2), calculated as [105]

Σ(E)=(000iπ(t2v2)δ(E)θ(tv)).Σ𝐸matrix000i𝜋superscript𝑡2superscript𝑣2𝛿𝐸𝜃𝑡𝑣\displaystyle\Sigma\left(E\right)=\begin{pmatrix}0&0\\ 0&-\text{i}\pi\left(t^{2}-v^{2}\right)\delta\left(E\right)\theta\left(t-v% \right)\end{pmatrix}.roman_Σ ( italic_E ) = ( start_ARG start_ROW start_CELL 0 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL - i italic_π ( italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_δ ( italic_E ) italic_θ ( italic_t - italic_v ) end_CELL end_ROW end_ARG ) . (4)

The complex eigenvalues of the effective non-Hermitian Hamiltonian Heff(0)=vσx+Σ(0)subscript𝐻eff0𝑣subscript𝜎𝑥Σ0H_{\rm eff}\left(0\right)=v\sigma_{x}+\Sigma\left(0\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 0 ) = italic_v italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + roman_Σ ( 0 ) are obtained as 00 and ii-\text{i}\infty- i ∞, the former of which corresponds to a topologically stable zero state with an infinite lifetime and the latter of which an unstable zero state with a vanishing lifetime. In the trivial phase (v/t>1𝑣𝑡1v/t>1italic_v / italic_t > 1), by contrast, no eigenstates appear at E=0𝐸0E=0italic_E = 0, and hence the self-energy vanishes Σ(0)=0Σ00\Sigma\left(0\right)=0roman_Σ ( 0 ) = 0.

Notably, the topological nature of the zero-energy edge state causes non-Hermitian topology of the effective Hamiltonian Heff(0)subscript𝐻eff0H_{\rm eff}\left(0\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 0 ). In the topological phase, the two single-particle eigenenergies are gapped with a reference energy on the imaginary axis [Fig. 2 (c)], which is a non-Hermitian extension of energy gap called a point gap [20, 29]. Inheriting from chiral symmetry of the original SSH model HSSH(k)subscript𝐻SSH𝑘H_{\rm SSH}\left(k\right)italic_H start_POSTSUBSCRIPT roman_SSH end_POSTSUBSCRIPT ( italic_k ) in Eq. (3), the edge non-Hermitian Hamiltonian Heff(0)subscript𝐻eff0H_{\rm eff}\left(0\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 0 ) also respects chiral symmetry σzHeff(0)σz=Heff(0)subscript𝜎𝑧superscriptsubscript𝐻eff0subscript𝜎𝑧subscript𝐻eff0\sigma_{z}H_{\rm eff}^{{\dagger}}\left(0\right)\sigma_{z}=-H_{\rm eff}\left(0\right)italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( 0 ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = - italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 0 ). Thanks to chiral symmetry, topologically-protected zero-energy states persist as long as the point gap is open. Such persistence is ensured by point-gap topology, given by the zeroth Chern number of the Hermitian matrix iHeff(0)σzisubscript𝐻eff0subscript𝜎𝑧\text{i}H_{\rm eff}\left(0\right)\sigma_{z}i italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 0 ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT [29]. In fact, Heff(0)=vσx+Σ(0)subscript𝐻eff0𝑣subscript𝜎𝑥Σ0H_{\rm eff}\left(0\right)=v\sigma_{x}+\Sigma\left(0\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 0 ) = italic_v italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + roman_Σ ( 0 ) in Eq. (4) exhibits the zeroth Chern number 1111 (00) in the topological (trivial) phase. Thus, the self-energy Σ(0)Σ0\Sigma\left(0\right)roman_Σ ( 0 ) detects nontrivial Hermitian topology in the bulk and induces non-Hermitian topology at the edge.

Refer to caption
Figure 3: Non-Hermitian topology in the Chern insulator (t=1.0𝑡1.0t=1.0italic_t = 1.0, m=1.3𝑚1.3m=-1.3italic_m = - 1.3, E=0𝐸0E=0italic_E = 0). The open boundary conditions (OBC) are imposed along the x𝑥xitalic_x direction. (a, b) Chiral edge state as a non-Hermitian subsystem of the Chern insulator. (c) Complex spectrum of the effective non-Hermitian Hamiltonian under the periodic boundary conditions (PBC; red; Lx=3000subscript𝐿𝑥3000L_{x}=3000italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 3000, Ly=30subscript𝐿𝑦30L_{y}=30italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 30) and OBC (blue; Lx=Ly=30subscript𝐿𝑥subscript𝐿𝑦30L_{x}=L_{y}=30italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 30) along the y𝑦yitalic_y direction. (d) Collection of all right (orange) and left (green) skin states. The infinitesimal number η𝜂\etaitalic_η is chosen as η=1/Ly𝜂1subscript𝐿𝑦\eta=1/\sqrt{L_{y}}italic_η = 1 / square-root start_ARG italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG [110].

Non-Hermitian topology in a Chern insulator.—Next, we consider a Chern insulator on the square lattice described by [111, *QHZ-08]

HChern(𝒌)subscript𝐻Chern𝒌\displaystyle H_{\mathrm{Chern}}\left(\bm{k}\right)italic_H start_POSTSUBSCRIPT roman_Chern end_POSTSUBSCRIPT ( bold_italic_k ) =(tsin(kx))σx+(tsin(ky))σyabsent𝑡subscript𝑘𝑥subscript𝜎𝑥𝑡subscript𝑘𝑦subscript𝜎𝑦\displaystyle=\left(t\sin{k_{x}}\right)\sigma_{x}+\left(t\sin{k_{y}}\right)% \sigma_{y}= ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT
+(m+tcos(kx)+tcos(ky))σz𝑚𝑡subscript𝑘𝑥𝑡subscript𝑘𝑦subscript𝜎𝑧\displaystyle\qquad+\left(m+t\cos{k_{x}}+t\cos{k_{y}}\right)\sigma_{z}+ ( italic_m + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT (5)

with t,m𝑡𝑚t,m\in\mathbb{R}italic_t , italic_m ∈ blackboard_R. We regard the one-dimensional edge at x=1𝑥1x=1italic_x = 1 as a system and the remaining bulk an environment, and show that the effective edge Hamiltonian yields non-Hermitian topology and concomitant skin effect (Fig. 3).

This model HChern(𝒌)subscript𝐻Chern𝒌H_{\mathrm{Chern}}\left(\bm{k}\right)italic_H start_POSTSUBSCRIPT roman_Chern end_POSTSUBSCRIPT ( bold_italic_k ) is characterized by the first Chern number, obtained as C1=sgn(m/t)subscript𝐶1sgn𝑚𝑡C_{1}=\mathrm{sgn}\left(m/t\right)italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = roman_sgn ( italic_m / italic_t ) (C1=0subscript𝐶10C_{1}=0italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0) for |m/t|<2𝑚𝑡2\left|m/t\right|<2| italic_m / italic_t | < 2 (|m/t|>2𝑚𝑡2\left|m/t\right|>2| italic_m / italic_t | > 2). Imposing the open boundary conditions along the x𝑥xitalic_x direction and periodic boundary conditions along the y𝑦yitalic_y directions, we analytically obtain the self-energy [105]

Σ(E,ky)=t2(m+tcos(ky))22(E+iηtsinky)(σ0σy).Σ𝐸subscript𝑘𝑦superscript𝑡2superscript𝑚𝑡subscript𝑘𝑦22𝐸i𝜂𝑡subscript𝑘𝑦subscript𝜎0subscript𝜎𝑦\displaystyle\Sigma\left(E,k_{y}\right)=\frac{t^{2}-\left(m+t\cos{k_{y}}\right% )^{2}}{2\left(E+\text{i}\eta-t\sin k_{y}\right)}\left(\sigma_{0}-\sigma_{y}% \right).roman_Σ ( italic_E , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_m + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( italic_E + i italic_η - italic_t roman_sin italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ( italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) . (6)

We calculate the complex eigenvalues of the effective non-Hermitian Hamiltonian Heff(E,ky)=(tsinky)σy+(m+tcosky)σz+Σ(E,ky)subscript𝐻eff𝐸subscript𝑘𝑦𝑡subscript𝑘𝑦subscript𝜎𝑦𝑚𝑡subscript𝑘𝑦subscript𝜎𝑧Σ𝐸subscript𝑘𝑦H_{\rm eff}\left(E,k_{y}\right)=\left(t\sin k_{y}\right)\sigma_{y}+\left(m+t% \cos k_{y}\right)\sigma_{z}+\Sigma\left(E,k_{y}\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = ( italic_t roman_sin italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + ( italic_m + italic_t roman_cos italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT + roman_Σ ( italic_E , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) [Fig. 3 (c)]. These eigenvalues form a loop in the complex plane and host a point gap for reference energy E0subscript𝐸0E_{0}italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT inside the loop. This loop structure intuitively originates from chirality of the anomalous boundary states, in a similar manner to the Hatano-Nelson model [106, *Hatano-Nelson-97]. The nontrivial point-gap topology is captured by the winding of the complex spectrum,

W1dk2πi(ddklogdet[H(k)E0])subscript𝑊1contour-integral𝑑𝑘2𝜋i𝑑𝑑𝑘𝐻𝑘subscript𝐸0W_{1}\coloneqq\oint\frac{dk}{2\pi\text{i}}\left(\frac{d}{dk}\log\det\left[H% \left(k\right)-E_{0}\right]\right)italic_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≔ ∮ divide start_ARG italic_d italic_k end_ARG start_ARG 2 italic_π i end_ARG ( divide start_ARG italic_d end_ARG start_ARG italic_d italic_k end_ARG roman_log roman_det [ italic_H ( italic_k ) - italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] ) (7)

for a Bloch Hamiltonian H(k)𝐻𝑘H\left(k\right)italic_H ( italic_k ) in one dimension [20, 29]. In fact, Heff(E,ky)subscript𝐻eff𝐸subscript𝑘𝑦H_{\rm eff}\left(E,k_{y}\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) with fixed E𝐸Eitalic_E exhibits W1=sgn(m/t)subscript𝑊1sgn𝑚𝑡W_{1}=\mathrm{sgn}\left(m/t\right)italic_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = roman_sgn ( italic_m / italic_t ), leading to the correspondence C1=W1subscript𝐶1subscript𝑊1C_{1}=W_{1}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT between the Hermitian bulk and non-Hermitian boundary [99, 100].

The bulk-boundary correspondence for the nontrivial point-gap topology W10subscript𝑊10W_{1}\neq 0italic_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≠ 0 manifests itself as the non-Hermitian skin effect [37, 38]. We calculate the complex spectrum of Heff(E,ky)subscript𝐻eff𝐸subscript𝑘𝑦H_{\rm eff}\left(E,k_{y}\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) under the open boundary conditions along both x𝑥xitalic_x and y𝑦yitalic_y directions [Fig. 3 (c)]. Clearly, the complex spectrum does not form a loop but an arc; such extreme sensitivity to the boundary conditions is a signature of the non-Hermitian skin effect [14, 22, 24]. Consistently, most of the right and left eigenstates of Heffsubscript𝐻effH_{\rm eff}italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT are localized at either corner of the square lattice [Fig. 3 (d)].

Refer to caption
Figure 4: Localized current distribution due to the non-Hermitian skin effect. The open boundary conditions (OBC) are imposed along the x𝑥xitalic_x direction. (a, b) While the inflow and outflow balance (a) under the periodic boundary conditions (PBC) along the y𝑦yitalic_y direction, they do not (b) under OBC. (c, d) Terminal current for various energy E𝐸Eitalic_E under (c) PBC and (d) OBC (t=1.0𝑡1.0t=1.0italic_t = 1.0, m=1.3𝑚1.3m=-1.3italic_m = - 1.3, Lx=Ly=30subscript𝐿𝑥subscript𝐿𝑦30L_{x}=L_{y}=30italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 30). (e) Inflow Jinflowsubscript𝐽inflowJ_{\rm inflow}italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT (green), net current J=JinflowJoutflow𝐽subscript𝐽inflowsubscript𝐽outflowJ=J_{\rm inflow}-J_{\rm outflow}italic_J = italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT - italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT (red), and outflow Joutflowsubscript𝐽outflow-J_{\rm outflow}- italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT (blue) for E=0𝐸0E=0italic_E = 0. (f) Localization length of the current J𝐽Jitalic_J (red) and the most localized skin state (green) as functions of m/t𝑚𝑡m/titalic_m / italic_t, obtained from the scaling for the five sites near the boundary. The infinitesimal number η𝜂\etaitalic_η is chosen as η=1/Ly𝜂1subscript𝐿𝑦\eta=1/\sqrt{L_{y}}italic_η = 1 / square-root start_ARG italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG [110].

Skin current.—We show that the non-Hermitian skin effect of the chiral edge states physically results in the localized current distribution (Fig. 4). Imposing the open boundary conditions along the x𝑥xitalic_x direction, we investigate the E𝐸Eitalic_E-resolved local current [105, 103, 113]

J(E)=[Hedge,Gedge(E)][Hedge,Gedge(E)],𝐽𝐸subscript𝐻edgesubscript𝐺edge𝐸superscriptsubscript𝐻edgesubscript𝐺edge𝐸J\left(E\right)=-\big{[}H_{\mathrm{edge}},G_{\mathrm{edge}}\left(E\right)\big{% ]}-\big{[}H_{\mathrm{edge}},G_{\mathrm{edge}}\left(E\right)\big{]}^{\dagger},italic_J ( italic_E ) = - [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ] - [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT , (8)

with the edge Green’s function Gedge(E)(E+iηHeff(E))1subscript𝐺edge𝐸superscript𝐸i𝜂subscript𝐻eff𝐸1G_{\rm edge}\left(E\right)\coloneqq\left(E+\text{i}\eta-H_{\mathrm{eff}}\left(% E\right)\right)^{-1}italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ≔ ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. Under the periodic boundary conditions along the y𝑦yitalic_y direction, the inflow Jinflowsubscript𝐽inflowJ_{\rm inflow}italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT and outflow Joutflowsubscript𝐽outflowJ_{\rm outflow}italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT between the bulk and boundary balance, leading to the absence of the net current J=JinflowJoutflow=0𝐽subscript𝐽inflowsubscript𝐽outflow0J=J_{\rm inflow}-J_{\rm outflow}=0italic_J = italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT - italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT = 0. Under the open boundary conditions along the y𝑦yitalic_y direction, by contrast, the nonzero net current J𝐽Jitalic_J arises around the corners. This unique current distribution results from the chiral edge states. In fact, chirality of the edge states leads to a large inflow at one corner and a large outflow at the opposite corner.

Notably, the localized current distribution is also a direct consequence of the corner skin effect. We find that the localization length of the current density shows behavior consistent with that of the skin states [Fig. 4 (f)]. A defining feature of the skin effect is the localization of right and left eigenstates at the opposite edges, as shown in Fig. 3 (d). From Eq. (8), a skin state with the localization length ξ𝜉\xiitalic_ξ contributes to the the current density as Jyj1ey/ξj2e(Ly)/ξsimilar-to-or-equalssubscript𝐽𝑦subscript𝑗1superscript𝑒𝑦𝜉subscript𝑗2superscript𝑒𝐿𝑦𝜉J_{y}\simeq j_{1}e^{-y/\xi}-j_{2}e^{-(L-y)/\xi}italic_J start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ≃ italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_y / italic_ξ end_POSTSUPERSCRIPT - italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - ( italic_L - italic_y ) / italic_ξ end_POSTSUPERSCRIPT with constants j1,j2subscript𝑗1subscript𝑗2j_{1},j_{2}\in\mathbb{R}italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ∈ blackboard_R, consistent with Fig. 4 (e, f). No local current arises if both right and left eigenstates of Heff(E)subscript𝐻eff𝐸H_{\rm eff}\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) are localized at the same boundary, as opposed to the skin states.

We also demonstrate the 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect in a time-reversal-invariant topological insulator [114, 105]. It is also notable that similar E𝐸Eitalic_E-dependent skin effect has been discussed in a recent work [115].

Refer to caption
Figure 5: Non-Hermitian topology in a three-dimensional topological insulator. (a) Hinge state in an effective non-Hermitian Hamiltonian at a surface. (b) Complex spectrum of the effective non-Hermitian Hamiltonian under the periodic boundary conditions along the y𝑦yitalic_y direction (y𝑦yitalic_y-PBC; orange; Lx=Ly=200subscript𝐿𝑥subscript𝐿𝑦200L_{x}=L_{y}=200italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 200, Lz=30subscript𝐿𝑧30L_{z}=30italic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 30) and open boundary conditions (y𝑦yitalic_y-OBC; blue; Lx=200subscript𝐿𝑥200L_{x}=200italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 200, Ly=Lz=30subscript𝐿𝑦subscript𝐿𝑧30L_{y}=L_{z}=30italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 30) (t=1.0𝑡1.0t=1.0italic_t = 1.0, m=1.6𝑚1.6m=-1.6italic_m = - 1.6, δ=0.2𝛿0.2\delta=0.2italic_δ = 0.2). The periodic boundary conditions are imposed along the x𝑥xitalic_x direction. The infinitesimal number η𝜂\etaitalic_η is chosen as η=1/Lz𝜂1subscript𝐿𝑧\eta=1/\sqrt{L_{z}}italic_η = 1 / square-root start_ARG italic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG [110].

Three dimensions.—We further study a three-dimensional topological insulator described by [4]

H3DTI(𝒌)=(m+tcos(kx)+tcos(ky)+tcos(kz))τysubscript𝐻3DTI𝒌𝑚𝑡subscript𝑘𝑥𝑡subscript𝑘𝑦𝑡subscript𝑘𝑧subscript𝜏𝑦\displaystyle H_{\mathrm{3DTI}}\left(\bm{k}\right)=\left(m+t\cos{k_{x}}+t\cos{% k_{y}}+t\cos{k_{z}}\right)\tau_{y}italic_H start_POSTSUBSCRIPT 3 roman_D roman_T roman_I end_POSTSUBSCRIPT ( bold_italic_k ) = ( italic_m + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG ) ) italic_τ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT
+(tsin(kx))σxτx+(tsin(ky))σyτx𝑡subscript𝑘𝑥subscript𝜎𝑥subscript𝜏𝑥𝑡subscript𝑘𝑦subscript𝜎𝑦subscript𝜏𝑥\displaystyle\qquad+\left(t\sin{k_{x}}\right)\sigma_{x}\tau_{x}+\left(t\sin{k_% {y}}\right)\sigma_{y}\tau_{x}+ ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT
+(tsin(kz))σzτx+δ(cos(kx)+cos(ky))σyτy𝑡subscript𝑘𝑧subscript𝜎𝑧subscript𝜏𝑥𝛿subscript𝑘𝑥subscript𝑘𝑦subscript𝜎𝑦subscript𝜏𝑦\displaystyle\qquad+\left(t\sin{k_{z}}\right)\sigma_{z}\tau_{x}+\delta\left(% \cos{k_{x}}+\cos{k_{y}}\right)\sigma_{y}\tau_{y}+ ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + italic_δ ( roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) + roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT (9)

with Pauli matrices σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s and τisubscript𝜏𝑖\tau_{i}italic_τ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s (i=x,y,z𝑖𝑥𝑦𝑧i=x,y,zitalic_i = italic_x , italic_y , italic_z), and real parameters m,t,δ𝑚𝑡𝛿m,t,\delta\in\mathbb{R}italic_m , italic_t , italic_δ ∈ blackboard_R. Similar to the SSH model, this Hamiltonian respects chiral symmetry τzH3DTI(𝒌)τz=H3DTI(𝒌)subscript𝜏𝑧subscript𝐻3DTI𝒌subscript𝜏𝑧subscript𝐻3DTI𝒌\tau_{z}H_{\mathrm{3DTI}}\left(\bm{k}\right)\tau_{z}=-H_{\mathrm{3DTI}}\left(% \bm{k}\right)italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT 3 roman_D roman_T roman_I end_POSTSUBSCRIPT ( bold_italic_k ) italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = - italic_H start_POSTSUBSCRIPT 3 roman_D roman_T roman_I end_POSTSUBSCRIPT ( bold_italic_k ) and is characterized by the three-dimensional winding number W3subscript𝑊3W_{3}\in\mathbb{Z}italic_W start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ∈ blackboard_Z [7]. In the topological phase, the nontrivial bulk topology W30subscript𝑊30W_{3}\neq 0italic_W start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ≠ 0 leads to the emergence of Dirac surface states. Applying the open boundary conditions along the z𝑧zitalic_z direction, we regard the two-dimensional surface at z=1𝑧1z=1italic_z = 1 as a system and the remaining bulk an environment (Fig. 5). The effective surface Hamiltonian Heff(E,kx,ky)subscript𝐻eff𝐸subscript𝑘𝑥subscript𝑘𝑦H_{\rm eff}\left(E,k_{x},k_{y}\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E , italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) hosts point-gap topology, characterized by the first Chern number of iHeff(E,kx,ky)τzisubscript𝐻eff𝐸subscript𝑘𝑥subscript𝑘𝑦subscript𝜏𝑧\text{i}H_{\rm eff}\left(E,k_{x},k_{y}\right)\tau_{z}i italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E , italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT. In contrast to the two-dimensional case, this point-gap topology does not result in the skin effect [38, 48]. Instead, anomalous hinge states with the imaginary dispersion ikxproportional-toabsentisubscript𝑘𝑥\propto\text{i}k_{x}∝ i italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT appear. The distinct boundary physics should reflect the different nature of the coupling between the bulk and boundary.

Table 1: Periodic table of Hermitian topological insulators and superconductors in spatial dimensions d=1,2,3𝑑123d=1,2,3italic_d = 1 , 2 , 3. The tenfold Altland-Zirnbauer (AZ) symmetry classes consist of time-reversal symmetry (TRS), particle-hole symmetry (PHS), and chiral symmetry (CS). In the entries specified by , effective non-Hermitian Hamiltonians at boundaries exhibit the skin effect; in 3D class DIII, specified by ∗∗, the skin effect arises only for the odd number of the topological invariant [48].
 AZ class  TRS  PHS  CS  d=1𝑑1d=1italic_d = 1  d=2𝑑2d=2italic_d = 2  d=3𝑑3d=3italic_d = 3
A 00 00 00 00 superscript\mathbb{Z}^{*}blackboard_Z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 00
AIII 00 00 1111 \mathbb{Z}blackboard_Z 00 \mathbb{Z}blackboard_Z
AI +11+1+ 1 00 00 0 0 0
BDI +11+1+ 1 +11+1+ 1 1111 \mathbb{Z}blackboard_Z 0 0
D 00 +11+1+ 1 00 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT superscript\mathbb{Z}^{*}blackboard_Z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 0
DIII 11-1- 1 +11+1+ 1 1111 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 2superscriptsubscript2\mathbb{Z}_{2}^{*}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT superscriptabsent\mathbb{Z}^{**}blackboard_Z start_POSTSUPERSCRIPT ∗ ∗ end_POSTSUPERSCRIPT
AII 11-1- 1 00 00 0 2superscriptsubscript2\mathbb{Z}_{2}^{*}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 2superscriptsubscript2\mathbb{Z}_{2}^{*}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT
CII 11-1- 1 11-1- 1 1111 222\mathbb{Z}2 blackboard_Z 0 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT
C 00 11-1- 1 00 0 22superscript2\mathbb{Z}^{*}2 blackboard_Z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT 0
CI +11+1+ 1 11-1- 1 1111 0 0 222\mathbb{Z}2 blackboard_Z

Classification.—While we have hitherto focused on the prototypical models, anomalous boundary states in topological insulators generally host non-Hermitian topology, as summarized in Table 1. In fact, when the original Hermitian Hamiltonians in Eq. (1) respect the Altland-Zirnbauer (AZ) symmetry [3, 7], the effective non-Hermitian Hamiltonians Heff(E)=Hedge+Σ(E)subscript𝐻eff𝐸subscript𝐻edgeΣ𝐸H_{\rm eff}\left(E\right)=H_{\rm edge}+\Sigma\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) = italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + roman_Σ ( italic_E ) respect the corresponding symmetry known as the AZ symmetry [105, 29]. Consequently, Heff(E)subscript𝐻eff𝐸H_{\rm eff}\left(E\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) shows d𝑑ditalic_d-dimensional non-Hermitian topology inheriting from (d+1)𝑑1\left(d+1\right)( italic_d + 1 )-dimensional Hermitian topology. In Table 1, various boundary phenomena are classified based on symmetry and dimensions [38, 48]. For example, helical Dirac surface states of time-reversal-invariant topological insulators host the 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect, contrasting with the chiral-symmetric topological insulator in Eq. (Non-Hermitian Topology in Hermitian Topological Matter).

Discussions.—In this Letter, we reveal intrinsic non-Hermitian topology of anomalous boundary states in Hermitian topological matter. We demonstrate that the self-energy quantifying the particle exchange between the bulk and boundary detects Hermitian topology in the bulk and non-Hermitian topology at the boundary. In Refs. [98, 99, 100], point-gap topology of anomalous boundary states was investigated in the presence of non-Hermitian external perturbations. In this Letter, by contrast, we show non-Hermitian topology inherently embedded within topological boundary states. Non-Hermitian topology was identified also in scattering matrices of Hermitian topological insulators [116, 117]. Conversely, we microscopically construct effective Hamiltonians and demonstrate their intrinsic non-Hermitian topology.

While we focus on clean noninteracting systems in this Letter, our formalism should be extended to topological matter with disorder or many-body interactions, which we investigate in future work. Additionally, it merits further study to develop a field-theoretical understanding of our correspondence. In this respect, our finding implies that anomaly of topological boundary states has a close connection with a different type of anomaly accompanying non-Hermitian topology [44]. It is also worthwhile to revisit fermion doubling from our perspective [118, 119, 120].

Acknowledgements.
We thank Akito Daido and Daichi Nakamura for helpful discussions. S.H. thanks Hosho Katsura for letting him know Ref. [121]. K.K. thanks Taylor L. Hughes for helpful discussions and appreciates the program “Topological Quantum Matter: Concepts and Realizations” held at Kavli Institute for Theoretical Physics (KITP), University of California Santa Barbara (National Science Foundation under Grant No. NSF PHY-1748958). We appreciate the long-term workshop “Quantum Information, Quantum Matter and Quantum Gravity” (YITP-T-23-01) held at Yukawa Institute for Theoretical Physics (YITP), Kyoto University. S.H. acknowledges travel support from MEXT KAKENHI Grant-in-Aid for Transformative Research Areas A “Extreme Universe” collaboration (Grant Nos. 21H05182 and 22H05247). S.H. is supported by JSPS Research Fellow, JSPS Overseas Challenge Program for Young Researchers, and MEXT WISE Program. S.H. and T.Y. are supported by JSPS KAKENHI Grant Nos. 21K13850 and 23KK0247. S.H. and T.Y are grateful for the support and hospitality of the Pauli Center for Theoretical Studies. K.K. is supported by MEXT KAKENHI Grant-in-Aid for Transformative Research Areas A “Extreme Universe” No. 24H00945.

References

  • Hasan and Kane [2010] M. Z. Hasan and C. L. Kane, Colloquium: Topological insulators, Rev. Mod. Phys. 82, 3045 (2010).
  • Qi and Zhang [2011] X.-L. Qi and S.-C. Zhang, Topological insulators and superconductors, Rev. Mod. Phys. 83, 1057 (2011).
  • Altland and Zirnbauer [1997] A. Altland and M. R. Zirnbauer, Nonstandard symmetry classes in mesoscopic normal-superconducting hybrid structures, Phys. Rev. B 55, 1142 (1997).
  • Schnyder et al. [2008] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Classification of topological insulators and superconductors in three spatial dimensions, Phys. Rev. B 78, 195125 (2008).
  • Ryu et al. [2010] S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Ludwig, Topological insulators and superconductors: tenfold way and dimensional hierarchy, New J. Phys. 12, 065010 (2010).
  • Kitaev [2009] A. Kitaev, Periodic table for topological insulators and superconductors, AIP Conf. Proc. 1134, 22 (2009).
  • Chiu et al. [2016] C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Classification of topological quantum matter with symmetries, Rev. Mod. Phys. 88, 035005 (2016).
  • Rudner and Levitov [2009] M. S. Rudner and L. S. Levitov, Topological Transition in a Non-Hermitian Quantum Walk, Phys. Rev. Lett. 102, 065703 (2009).
  • Sato et al. [2012] M. Sato, K. Hasebe, K. Esaki, and M. Kohmoto, Time-Reversal Symmetry in Non-Hermitian Systems, Prog. Theor. Phys. 127, 937 (2012).
  • Esaki et al. [2011] K. Esaki, M. Sato, K. Hasebe, and M. Kohmoto, Edge states and topological phases in non-Hermitian systems, Phys. Rev. B 84, 205128 (2011).
  • Hu and Hughes [2011] Y. C. Hu and T. L. Hughes, Absence of topological insulator phases in non-Hermitian PT-symmetric Hamiltonians, Phys. Rev. B 84, 153101 (2011).
  • Schomerus [2013] H. Schomerus, Topologically protected midgap states in complex photonic lattices, Opt. Lett. 38, 1912 (2013).
  • Longhi et al. [2015] S. Longhi, D. Gatti, and G. D. Valle, Robust light transport in non-Hermitian photonic lattices, Sci. Rep. 5, 13376 (2015).
  • Lee [2016] T. E. Lee, Anomalous Edge State in a Non-Hermitian Lattice, Phys. Rev. Lett. 116, 133903 (2016).
  • Leykam et al. [2017] D. Leykam, K. Y. Bliokh, C. Huang, Y. D. Chong, and F. Nori, Edge Modes, Degeneracies, and Topological Numbers in Non-Hermitian Systems, Phys. Rev. Lett. 118, 040401 (2017).
  • Xu et al. [2017] Y. Xu, S.-T. Wang, and L.-M. Duan, Weyl Exceptional Rings in a Three-Dimensional Dissipative Cold Atomic Gas, Phys. Rev. Lett. 118, 045701 (2017).
  • Shen et al. [2018] H. Shen, B. Zhen, and L. Fu, Topological Band Theory for Non-Hermitian Hamiltonians, Phys. Rev. Lett. 120, 146402 (2018).
  • Takata and Notomi [2018] K. Takata and M. Notomi, Photonic Topological Insulating Phase Induced Solely by Gain and Loss, Phys. Rev. Lett. 121, 213902 (2018).
  • Martinez Alvarez et al. [2018] V. M. Martinez Alvarez, J. E. Barrios Vargas, and L. E. F. Foa Torres, Non-Hermitian robust edge states in one dimension: Anomalous localization and eigenspace condensation at exceptional points, Phys. Rev. B 97, 121401(R) (2018).
  • Gong et al. [2018] Z. Gong, Y. Ashida, K. Kawabata, K. Takasan, S. Higashikawa, and M. Ueda, Topological Phases of Non-Hermitian Systems, Phys. Rev. X 8, 031079 (2018).
  • Kawabata et al. [2019a] K. Kawabata, S. Higashikawa, Z. Gong, Y. Ashida, and M. Ueda, Topological unification of time-reversal and particle-hole symmetries in non-Hermitian physics, Nat. Commun. 10, 297 (2019a).
  • Yao and Wang [2018] S. Yao and Z. Wang, Edge States and Topological Invariants of Non-Hermitian Systems, Phys. Rev. Lett. 121, 086803 (2018).
  • Yao et al. [2018] S. Yao, F. Song, and Z. Wang, Non-Hermitian Chern Bands, Phys. Rev. Lett. 121, 136802 (2018).
  • Kunst et al. [2018] F. K. Kunst, E. Edvardsson, J. C. Budich, and E. J. Bergholtz, Biorthogonal Bulk-Boundary Correspondence in Non-Hermitian Systems, Phys. Rev. Lett. 121, 026808 (2018).
  • McDonald et al. [2018] A. McDonald, T. Pereg-Barnea, and A. A. Clerk, Phase-Dependent Chiral Transport and Effective Non-Hermitian Dynamics in a Bosonic Kitaev-Majorana Chain, Phys. Rev. X 8, 041031 (2018).
  • Lee and Thomale [2019] C. H. Lee and R. Thomale, Anatomy of skin modes and topology in non-Hermitian systems, Phys. Rev. B 99, 201103(R) (2019).
  • Liu et al. [2019] T. Liu, Y.-R. Zhang, Q. Ai, Z. Gong, K. Kawabata, M. Ueda, and F. Nori, Second-Order Topological Phases in Non-Hermitian Systems, Phys. Rev. Lett. 122, 076801 (2019).
  • Lee et al. [2019a] C. H. Lee, L. Li, and J. Gong, Hybrid Higher-Order Skin-Topological Modes in Nonreciprocal Systems, Phys. Rev. Lett. 123, 016805 (2019a).
  • Kawabata et al. [2019b] K. Kawabata, K. Shiozaki, M. Ueda, and M. Sato, Symmetry and Topology in Non-Hermitian Physics, Phys. Rev. X 9, 041015 (2019b).
  • Zhou and Lee [2019] H. Zhou and J. Y. Lee, Periodic table for topological bands with non-Hermitian symmetries, Phys. Rev. B 99, 235112 (2019).
  • Herviou et al. [2019] L. Herviou, J. H. Bardarson, and N. Regnault, Defining a bulk-edge correspondence for non-Hermitian Hamiltonians via singular-value decomposition, Phys. Rev. A 99, 052118 (2019).
  • Zirnstein et al. [2021] H.-G. Zirnstein, G. Refael, and B. Rosenow, Bulk-Boundary Correspondence for Non-Hermitian Hamiltonians via Green Functions, Phys. Rev. Lett. 126, 216407 (2021).
  • Borgnia et al. [2020] D. S. Borgnia, A. J. Kruchkov, and R.-J. Slager, Non-Hermitian Boundary Modes and Topology, Phys. Rev. Lett. 124, 056802 (2020).
  • Yokomizo and Murakami [2019] K. Yokomizo and S. Murakami, Non-Bloch Band Theory of Non-Hermitian Systems, Phys. Rev. Lett. 123, 066404 (2019).
  • Chang et al. [2020] P.-Y. Chang, J.-S. You, X. Wen, and S. Ryu, Entanglement spectrum and entropy in topological non-Hermitian systems and nonunitary conformal field theory, Phys. Rev. Research 2, 033069 (2020).
  • Wanjura et al. [2020] C. C. Wanjura, M. Brunelli, and A. Nunnenkamp, Topological framework for directional amplification in driven-dissipative cavity arrays, Nat. Commun. 11, 3149 (2020).
  • Zhang et al. [2020] K. Zhang, Z. Yang, and C. Fang, Correspondence between Winding Numbers and Skin Modes in Non-Hermitian Systems, Phys. Rev. Lett. 125, 126402 (2020).
  • Okuma et al. [2020] N. Okuma, K. Kawabata, K. Shiozaki, and M. Sato, Topological Origin of Non-Hermitian Skin Effects, Phys. Rev. Lett. 124, 086801 (2020).
  • Denner et al. [2021] M. M. Denner, A. Skurativska, F. Schindler, M. H. Fischer, R. Thomale, T. Bzdušek, and T. Neupert, Exceptional topological insulators, Nat. Commun. 12, 5681 (2021).
  • Denner et al. [2023] M. M. Denner, T. Neupert, and F. Schindler, Infernal and exceptional edge modes: non-Hermitian topology beyond the skin effect, J. Phys. Mater. 6, 045006 (2023).
  • Okugawa et al. [2020] R. Okugawa, R. Takahashi, and K. Yokomizo, Second-order topological non-Hermitian skin effects, Phys. Rev. B 102, 241202(R) (2020).
  • Kawabata et al. [2020] K. Kawabata, M. Sato, and K. Shiozaki, Higher-order non-Hermitian skin effect, Phys. Rev. B 102, 205118 (2020).
  • Fu et al. [2021] Y. Fu, J. Hu, and S. Wan, Non-Hermitian second-order skin and topological modes, Phys. Rev. B 103, 045420 (2021).
  • Kawabata et al. [2021] K. Kawabata, K. Shiozaki, and S. Ryu, Topological Field Theory of Non-Hermitian Systems, Phys. Rev. Lett. 126, 216405 (2021).
  • Zhang et al. [2022] K. Zhang, Z. Yang, and C. Fang, Universal non-Hermitian skin effect in two and higher dimensions, Nat. Commun. 13, 2496 (2022).
  • Sun et al. [2021] X.-Q. Sun, P. Zhu, and T. L. Hughes, Geometric Response and Disclination-Induced Skin Effects in Non-Hermitian Systems, Phys. Rev. Lett. 127, 066401 (2021).
  • Kim and Park [2021] K.-M. Kim and M. J. Park, Disorder-driven phase transition in the second-order non-Hermitian skin effect, Phys. Rev. B 104, L121101 (2021).
  • Nakamura et al. [2024] D. Nakamura, T. Bessho, and M. Sato, Bulk-Boundary Correspondence in Point-Gap Topological Phases, Phys. Rev. Lett. 132, 136401 (2024).
  • Wang et al. [2024] H.-Y. Wang, F. Song, and Z. Wang, Amoeba Formulation of Non-Bloch Band Theory in Arbitrary Dimensions, Phys. Rev. X 14, 021011 (2024).
  • Nakai et al. [2024] Y. O. Nakai, N. Okuma, D. Nakamura, K. Shimomura, and M. Sato, Topological enhancement of nonnormality in non-Hermitian skin effects, Phys. Rev. B 109, 144203 (2024).
  • Bergholtz et al. [2021] E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Exceptional topology of non-Hermitian systems, Rev. Mod. Phys. 93, 015005 (2021).
  • Okuma and Sato [2023] N. Okuma and M. Sato, Non-Hermitian Topological Phenomena: A Review, Annu. Rev. Condens. Matter Phys. 14, 83 (2023).
  • Konotop et al. [2016] V. V. Konotop, J. Yang, and D. A. Zezyulin, Nonlinear waves in 𝒫𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetric systems, Rev. Mod. Phys. 88, 035002 (2016).
  • El-Ganainy et al. [2018] R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, and D. N. Christodoulides, Non-Hermitian physics and PT symmetry, Nat. Phys. 14, 11 (2018).
  • [55] V. Kozii and L. Fu, Non-Hermitian Topological Theory of Finite-Lifetime Quasiparticles: Prediction of Bulk Fermi Arc Due to Exceptional Point, arXiv:1708.05841 .
  • Papaj et al. [2019] M. Papaj, H. Isobe, and L. Fu, Nodal arc of disordered Dirac fermions and non-Hermitian band theory, Phys. Rev. B 99, 201107 (2019).
  • Shen and Fu [2018] H. Shen and L. Fu, Quantum Oscillation from In-Gap States and a Non-Hermitian Landau Level Problem, Phys. Rev. Lett. 121, 026403 (2018).
  • Yoshida et al. [2018] T. Yoshida, R. Peters, and N. Kawakami, Non-Hermitian perspective of the band structure in heavy-fermion systems, Phys. Rev. B 98, 035141 (2018).
  • Hirsbrunner et al. [2019] M. R. Hirsbrunner, T. M. Philip, and M. J. Gilbert, Topology and observables of the non-Hermitian Chern insulator, Phys. Rev. B 100, 081104(R) (2019).
  • McClarty and Rau [2019] P. A. McClarty and J. G. Rau, Non-Hermitian topology of spontaneous magnon decay, Phys. Rev. B 100, 100405(R) (2019).
  • Bergholtz and Budich [2019] E. J. Bergholtz and J. C. Budich, Non-Hermitian Weyl physics in topological insulator ferromagnet junctions, Phys. Rev. Research 1, 012003(R) (2019).
  • Michishita and Peters [2020] Y. Michishita and R. Peters, Equivalence of Effective Non-Hermitian Hamiltonians in the Context of Open Quantum Systems and Strongly Correlated Electron Systems, Phys. Rev. Lett. 124, 196401 (2020).
  • Okuma and Sato [2021] N. Okuma and M. Sato, Non-Hermitian Skin Effects in Hermitian Correlated or Disordered Systems: Quantities Sensitive or Insensitive to Boundary Effects and Pseudo-Quantum-Number, Phys. Rev. Lett. 126, 176601 (2021).
  • Lehmann et al. [2021] C. Lehmann, M. Schüler, and J. C. Budich, Dynamically Induced Exceptional Phases in Quenched Interacting Semimetals, Phys. Rev. Lett. 127, 106601 (2021).
  • [65] J. Cayao and M. Sato, Non-Hermitian phase-biased Josephson junctions, arXiv:2307.15472 .
  • Poli et al. [2015] C. Poli, M. Bellec, U. Kuhl, F. Mortessagne, and H. Schomerus, Selective enhancement of topologically induced interface states in a dielectric resonator chain, Nat. Commun. 6, 6710 (2015).
  • Zeuner et al. [2015] J. M. Zeuner, M. C. Rechtsman, Y. Plotnik, Y. Lumer, S. Nolte, M. S. Rudner, M. Segev, and A. Szameit, Observation of a Topological Transition in the Bulk of a Non-Hermitian System, Phys. Rev. Lett. 115, 040402 (2015).
  • Zhen et al. [2015] B. Zhen, C. W. Hsu, Y. Igarashi, L. Lu, I. Kaminer, A. Pick, S.-L. Chua, J. D. Joannopoulos, and M. Soljac̆ić, Spawning rings of exceptional points out of Dirac cones, Nature 525, 354 (2015).
  • Weimann et al. [2017] S. Weimann, M. Kremer, Y. Plotnik, Y. Lumer, S. Nolte, K. G. Makris, M. Segev, M. C. Rechtsman, and A. Szameit, Topologically protected bound states in photonic parity-time-symmetric crystals, Nat. Mater. 16, 433 (2017).
  • Xiao et al. [2017] L. Xiao, X. Zhan, Z. H. Bian, K. K. Wang, X. Zhang, X. P. Wang, J. Li, K. Mochizuki, D. Kim, N. Kawakami, W. Yi, H. Obuse, B. C. Sanders, and P. Xue, Observation of topological edge states in parity-time-symmetric quantum walks, Nat. Phys. 13, 1117 (2017).
  • St-Jean et al. [2017] P. St-Jean, V. Goblot, E. Galopin, A. Lemaître, T. Ozawa, L. L. Gratiet, I. Sagnes, J. Bloch, and A. Amo, Lasing in topological edge states of a one-dimensional lattice, Nat. Photon. 11, 651 (2017).
  • Parto et al. [2018] M. Parto, S. Wittek, H. Hodaei, G. Harari, M. A. Bandres, J. Ren, M. C. Rechtsman, M. Segev, D. N. Christodoulides, and M. Khajavikhan, Edge-Mode Lasing in 1D Topological Active Arrays, Phys. Rev. Lett. 120, 113901 (2018).
  • Bahari et al. [2017] B. Bahari, A. Ndao, F. Vallini, A. E. Amili, Y. Fainman, and B. Kanté, Nonreciprocal lasing in topological cavities of arbitrary geometries, Science 358, 636 (2017).
  • Zhao et al. [2018] H. Zhao, P. Miao, M. H. Teimourpour, S. Malzard, R. El-Ganainy, H. Schomerus, and L. Feng, Topological hybrid silicon microlasers, Nat. Commun. 9, 981 (2018).
  • Zhou et al. [2018] H. Zhou, C. Peng, Y. Yoon, C. W. Hsu, K. A. Nelson, L. Fu, J. D. Joannopoulos, M. Soljac̆ić, and B. Zhen, Observation of bulk Fermi arc and polarization half charge from paired exceptional points, Science 359, 1009 (2018).
  • Harari et al. [2018] G. Harari, M. A. Bandres, Y. Lumer, M. C. Rechtsman, Y. D. Chong, M. Khajavikhan, D. N. Christodoulides, and M. Segev, Topological insulator laser: Theory, Science 359, eaar4003 (2018).
  • Bandres et al. [2018] M. A. Bandres, S. Wittek, G. Harari, M. Parto, J. Ren, M. Segev, D. Christodoulides, and M. Khajavikhan, Topological insulator laser: Experiments, Science 359, eaar4005 (2018).
  • Cerjan et al. [2019] A. Cerjan, S. Huang, K. P. Chen, Y. Chong, and M. C. Rechtsman, Experimental realization of a Weyl exceptional ring, Nat. Photon. 13, 623 (2019).
  • Zhao et al. [2019] H. Zhao, X. Qiao, T. Wu, B. Midya, S. Longhi, and L. Feng, Non-Hermitian topological light steering, Science 365, 1163 (2019).
  • Brandenbourger et al. [2019] M. Brandenbourger, X. Locsin, E. Lerner, and C. Coulais, Non-reciprocal robotic metamaterials, Nat. Commun. 10, 4608 (2019).
  • Ghatak et al. [2020] A. Ghatak, M. Brandenbourger, J. van Wezel, and C. Coulais, Observation of non-Hermitian topology and its bulk-edge correspondence in an active mechanical metamaterial, Proc. Natl. Acad. Sci. USA 117, 29561 (2020).
  • Helbig et al. [2020] T. Helbig, T. Hofmann, S. Imhof, M. Abdelghany, T. Kiessling, L. W. Molenkamp, C. H. Lee, A. Szameit, M. Greiter, and R. Thomale, Generalized bulk-boundary correspondence in non-Hermitian topolectrical circuits, Nat. Phys. 16, 747 (2020).
  • Hofmann et al. [2020] T. Hofmann, T. Helbig, F. Schindler, N. Salgo, M. Brzezińska, M. Greiter, T. Kiessling, D. Wolf, A. Vollhardt, A. Kabaši, C. H. Lee, A. Bilušić, R. Thomale, and T. Neupert, Reciprocal skin effect and its realization in a topolectrical circuit, Phys. Rev. Research 2, 023265 (2020).
  • Xiao et al. [2020] L. Xiao, T. Deng, K. Wang, G. Zhu, Z. Wang, W. Yi, and P. Xue, Non-Hermitian bulk-boundary correspondence in quantum dynamics, Nat. Phys. 16, 761 (2020).
  • Weidemann et al. [2020] S. Weidemann, M. Kremer, T. Helbig, T. Hofmann, A. Stegmaier, M. Greiter, R. Thomale, and A. Szameit, Topological funneling of light, Science 368, 311 (2020).
  • Gou et al. [2020] W. Gou, T. Chen, D. Xie, T. Xiao, T.-S. Deng, B. Gadway, W. Yi, and B. Yan, Tunable Nonreciprocal Quantum Transport through a Dissipative Aharonov-Bohm Ring in Ultracold Atoms, Phys. Rev. Lett. 124, 070402 (2020).
  • Liang et al. [2022] Q. Liang, D. Xie, Z. Dong, H. Li, H. Li, B. Gadway, W. Yi, and B. Yan, Dynamic Signatures of Non-Hermitian Skin Effect and Topology in Ultracold Atoms, Phys. Rev. Lett. 129, 070401 (2022).
  • Palacios et al. [2020] L. S. Palacios, S. Tchoumakov, M. Guix, I. P. S. Sánchez, and A. G. Grushin, Guided accumulation of active particles by topological design of a second-order skin effect, Nat. Commun. 12, 4691 (2020).
  • Zhang et al. [2021] X. Zhang, Y. Tian, J.-H. Jiang, M.-H. Lu, and Y.-F. Chen, Observation of higher-order non-Hermitian skin effect, Nat. Commun. 12, 5377 (2021).
  • Wang et al. [2023] A. Wang, Z. Meng, and C. Q. Chen, Non-Hermitian topology in static mechanical metamaterials, Sci. Adv. 9, eadf7299 (2023).
  • [91] R. Shen, T. Chen, B. Yang, and C. H. Lee, Observation of the non-Hermitian skin effect and Fermi skin on a digital quantum computer, arXiv:2311.10143 .
  • Lee et al. [2019b] J. Y. Lee, J. Ahn, H. Zhou, and A. Vishwanath, Topological Correspondence between Hermitian and Non-Hermitian Systems: Anomalous Dynamics, Phys. Rev. Lett. 123, 206404 (2019b).
  • Bessho and Sato [2021] T. Bessho and M. Sato, Nielsen-Ninomiya Theorem with Bulk Topology: Duality in Floquet and Non-Hermitian Systems, Phys. Rev. Lett. 127, 196404 (2021).
  • Li et al. [2022] Y. Li, C. Liang, C. Wang, C. Lu, and Y.-C. Liu, Gain-Loss-Induced Hybrid Skin-Topological Effect, Phys. Rev. Lett. 128, 223903 (2022).
  • Zhu and Gong [2022] W. Zhu and J. Gong, Hybrid skin-topological modes without asymmetric couplings, Phys. Rev. B 106, 035425 (2022).
  • Zhu and Gong [2023] W. Zhu and J. Gong, Photonic corner skin modes in non-Hermitian photonic crystals, Phys. Rev. B 108, 035406 (2023).
  • Ou et al. [2023] Z. Ou, Y. Wang, and L. Li, Non-Hermitian boundary spectral winding, Phys. Rev. B 107, L161404 (2023).
  • Ma et al. [2024] X.-R. Ma, K. Cao, X.-R. Wang, Z. Wei, Q. Du, and S.-P. Kou, Non-Hermitian chiral skin effect, Phys. Rev. Research 6, 013213 (2024).
  • Schindler et al. [2023] F. Schindler, K. Gu, B. Lian, and K. Kawabata, Hermitian Bulk – Non-Hermitian Boundary Correspondence, PRX Quantum 4, 030315 (2023).
  • Nakamura et al. [2023] D. Nakamura, K. Inaka, N. Okuma, and M. Sato, Universal Platform of Point-Gap Topological Phases from Topological Materials, Phys. Rev. Lett. 131, 256602 (2023).
  • Liu et al. [2024] G.-G. Liu, S. Mandal, P. Zhou, X. Xi, R. Banerjee, Y.-H. Hu, M. Wei, M. Wang, Q. Wang, Z. Gao, H. Chen, Y. Yang, Y. Chong, and B. Zhang, Localization of Chiral Edge States by the Non-Hermitian Skin Effect, Phys. Rev. Lett. 132, 113802 (2024).
  • [102] J. Wu, R. Zheng, J. Liang, M. Ke, J. Lu, W. Deng, X. Huang, and Z. Liu, Spin-dependent localization of helical edge states in a non-Hermitian phononic crystal, arXiv:2312.12060 .
  • Datta [1995] S. Datta, Electronic Transport in Mesoscopic Systems (Cambridge University Press, Cambridge, England, 1995).
  • Moiseyev [2011] N. Moiseyev, Non-Hermitian Quantum Mechanics (Cambridge University Press, Cambridge, England, 2011).
  • [105] See the Supplemental Material, which includes Refs. [122, 123, 124, 121, 125], for details on effective non-Hermitian Hamiltonians, symmetry classification of non-Hermitian self-energy, analytical derivations of self-energy, current formulas, and 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect in time-reversal-invariant topological insulators.
  • Hatano and Nelson [1996] N. Hatano and D. R. Nelson, Localization Transitions in Non-Hermitian Quantum Mechanics, Phys. Rev. Lett. 77, 570 (1996).
  • Hatano and Nelson [1997] N. Hatano and D. R. Nelson, Vortex pinning and non-Hermitian quantum mechanics, Phys. Rev. B 56, 8651 (1997).
  • Hatano [2021] N. Hatano, What is the resonant state in open quantum systems?, J. Phys.: Conf. Ser. 2038, 012013 (2021).
  • Su et al. [1979] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons in Polyacetylene, Phys. Rev. Lett. 42, 1698 (1979).
  • [110] η=1/L𝜂1𝐿\eta=1/\sqrt{L}italic_η = 1 / square-root start_ARG italic_L end_ARG is chosen to be much larger than the level spacing 1/Lsimilar-toabsent1𝐿\sim 1/L∼ 1 / italic_L but much smaller than the macroscopic energy scale 1similar-toabsent1\sim 1∼ 1.
  • Qi et al. [2006] X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Topological quantization of the spin Hall effect in two-dimensional paramagnetic semiconductors, Phys. Rev. B 74, 085308 (2006).
  • Qi et al. [2008] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Topological field theory of time-reversal invariant insulators, Phys. Rev. B 78, 195424 (2008).
  • Datta [2005] S. Datta, Quantum Transport (Cambridge University Press, Cambridge, England, 2005).
  • Bernevig et al. [2006] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Quantum Spin Hall Effect and Topological Phase Transition in HgTe Quantum Wells, Science 314, 1757 (2006).
  • Eek et al. [2024] L. Eek, A. Moustaj, M. Röntgen, V. Pagneux, V. Achilleos, and C. M. Smith, Emergent non-Hermitian models, Phys. Rev. B 109, 045122 (2024).
  • Franca et al. [2022] S. Franca, V. Könye, F. Hassler, J. van den Brink, and C. Fulga, Non-Hermitian Physics without Gain or Loss: The Skin Effect of Reflected Waves, Phys. Rev. Lett. 129, 086601 (2022).
  • Ochkan et al. [2024] K. Ochkan, R. Chaturvedi, V. Könye, L. Veyrat, R. Giraud, D. Mailly, A. Cavanna, U. Gennser, E. M. Hankiewicz, B. Büchner, J. van den Brink, J. Dufouleur, and I. C. Fulga, Non-Hermitian topology in a multi-terminal quantum Hall device, Nat. Phys. 20, 395 (2024).
  • Fujikawa and Suzuki [2004] K. Fujikawa and H. Suzuki, Path Integrals and Quantum Anomalies (Oxford University Press, Oxford, 2004).
  • Karsten [1981] L. H. Karsten, Lattice fermions in euclidean space-time, Phys. Lett. B 104, 315 (1981).
  • Chen et al. [2023] W.-Q. Chen, Y.-S. Wu, W. Xi, W.-Z. Yi, and G. Yue, Fate of Quantum Anomalies for 1d lattice chiral fermion with a simple non-Hermitian Hamiltonian, J. High Energ. Phys. 2023 (5), 90.
  • Sirker et al. [2014] J. Sirker, M. Maiti, N. P. Konstantinidis, and N. Sedlmayr, Boundary fidelity and entanglement in the symmetry protected topological phase of the SSH model, J. Stat. Mech. 2014, P10032 (2014).
  • Ryndyk et al. [2009] D. A. Ryndyk, R. Gutiérrez, B. Song, and G. Cuniberti, Green Function Techniques in the Treatment of Quantum Transport at the Molecular Scale, in Energy Transfer Dynamics in Biomaterial Systems, edited by I. Burghardt, V. May, D. A. Micha, and E. R. Bittner (Springer, Berlin, Heidelberg, 2009) pp. 213–335.
  • Lieu et al. [2020] S. Lieu, M. McGinley, and N. R. Cooper, Tenfold Way for Quadratic Lindbladians, Phys. Rev. Lett. 124, 040401 (2020).
  • Kawabata et al. [2023] K. Kawabata, Z. Xiao, T. Ohtsuki, and R. Shindou, Singular-Value Statistics of Non-Hermitian Random Matrices and Open Quantum Systems, PRX Quantum 4, 040312 (2023).
  • Shin [1997] B. C. Shin, A formula for Eigenpairs of certain symmetric tridiagonal matrices, Bull. Aust. Math. Soc. 55, 249 (1997).

Supplemental Material for
“Non-Hermitian Topology in Hermitian Topological Matter”

I Effective non-Hermitian Hamiltonian

We consider an edge Hamiltonian Hedgesubscript𝐻edgeH_{\rm edge}italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT that is attached to a bulk Hamiltonian Hbulksubscript𝐻bulkH_{\rm bulk}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT. The whole system is described as

H=(HbulkTTHedge),𝐻matrixsubscript𝐻bulk𝑇superscript𝑇subscript𝐻edgeH=\begin{pmatrix}H_{\rm bulk}&T\\ T^{{\dagger}}&H_{\rm edge}\end{pmatrix},italic_H = ( start_ARG start_ROW start_CELL italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT end_CELL start_CELL italic_T end_CELL end_ROW start_ROW start_CELL italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_CELL start_CELL italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) , (I.1)

where T𝑇Titalic_T denotes the coupling between the bulk and the edge. The Green’s function G𝐺Gitalic_G of H𝐻Hitalic_H is defined by [103]

[(E+iη)H]G=1,delimited-[]𝐸i𝜂𝐻𝐺1\left[\left(E+\text{i}\eta\right)-H\right]G=1,[ ( italic_E + i italic_η ) - italic_H ] italic_G = 1 , (I.2)

or

(E+iηHbulkTTE+iηHedge)(GbulkGbulkedgeGedgebulkGedge)=1,matrix𝐸i𝜂subscript𝐻bulk𝑇superscript𝑇𝐸i𝜂subscript𝐻edgematrixsubscript𝐺bulksubscript𝐺bulkedgesubscript𝐺edgebulksubscript𝐺edge1\begin{pmatrix}E+\text{i}\eta-H_{\rm bulk}&-T\\ -T^{{\dagger}}&E+\text{i}\eta-H_{\rm edge}\end{pmatrix}\begin{pmatrix}G_{\rm bulk% }&G_{\rm bulk-edge}\\ G_{\rm edge-bulk}&G_{\rm edge}\end{pmatrix}=1,( start_ARG start_ROW start_CELL italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT end_CELL start_CELL - italic_T end_CELL end_ROW start_ROW start_CELL - italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_CELL start_CELL italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) ( start_ARG start_ROW start_CELL italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT end_CELL start_CELL italic_G start_POSTSUBSCRIPT roman_bulk - roman_edge end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_G start_POSTSUBSCRIPT roman_edge - roman_bulk end_POSTSUBSCRIPT end_CELL start_CELL italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) = 1 , (I.3)

with an infinitesimal number η>0𝜂0\eta>0italic_η > 0. This equation reads

(E+iηHbulk)GbulkedgeTGedge𝐸i𝜂subscript𝐻bulksubscript𝐺bulkedge𝑇subscript𝐺edge\displaystyle\left(E+\text{i}\eta-H_{\rm bulk}\right)G_{\rm bulk-edge}-TG_{\rm edge}( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) italic_G start_POSTSUBSCRIPT roman_bulk - roman_edge end_POSTSUBSCRIPT - italic_T italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT =0,absent0\displaystyle=0,= 0 , (I.4)
TGbulkedge+(E+iηHedge)Gedgesuperscript𝑇subscript𝐺bulkedge𝐸i𝜂subscript𝐻edgesubscript𝐺edge\displaystyle-T^{{\dagger}}G_{\rm bulk-edge}+\left(E+\text{i}\eta-H_{\rm edge}% \right)G_{\rm edge}- italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_G start_POSTSUBSCRIPT roman_bulk - roman_edge end_POSTSUBSCRIPT + ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ) italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT =1,absent1\displaystyle=1,= 1 , (I.5)

leading to

Gbulkedge=(E+iηHbulk)1TGedgeGbulk(0)TGedge,subscript𝐺bulkedgesuperscript𝐸i𝜂subscript𝐻bulk1𝑇subscript𝐺edgesuperscriptsubscript𝐺bulk0𝑇subscript𝐺edgeG_{\rm bulk-edge}=\left(E+\text{i}\eta-H_{\rm bulk}\right)^{-1}TG_{\rm edge}% \eqqcolon G_{\rm bulk}^{(0)}TG_{\rm edge},italic_G start_POSTSUBSCRIPT roman_bulk - roman_edge end_POSTSUBSCRIPT = ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ≕ italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT italic_T italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , (I.6)

with the Green’s function Gbulk(0)(E+iηHbulk)1superscriptsubscript𝐺bulk0superscript𝐸i𝜂subscript𝐻bulk1G_{\rm bulk}^{(0)}\coloneqq\left(E+\text{i}\eta-H_{\rm bulk}\right)^{-1}italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ≔ ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT of the isolated bulk. Moreover, we have

(E+iηHedgeΣ)Gedge=1,𝐸i𝜂subscript𝐻edgeΣsubscript𝐺edge1\left(E+\text{i}\eta-H_{\rm edge}-\Sigma\right)G_{\rm edge}=1,( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT - roman_Σ ) italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT = 1 , (I.7)

where we introduce the self-energy ΣΣ\Sigmaroman_Σ of the bulk as

ΣT(E+iηHbulk)1T=TGbulk(0)T.Σsuperscript𝑇superscript𝐸i𝜂subscript𝐻bulk1𝑇superscript𝑇superscriptsubscript𝐺bulk0𝑇\Sigma\coloneqq T^{{\dagger}}\left(E+\text{i}\eta-H_{\rm bulk}\right)^{-1}T=T^% {{\dagger}}G_{\rm bulk}^{(0)}T.roman_Σ ≔ italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T = italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT italic_T . (I.8)

The effective non-Hermitian Hamiltonian at the edge reads

HeffHedge+Σ=Hedge+T(E+iηHbulk)1T.subscript𝐻effsubscript𝐻edgeΣsubscript𝐻edgesuperscript𝑇superscript𝐸i𝜂subscript𝐻bulk1𝑇H_{\rm eff}\coloneqq H_{\rm edge}+\Sigma=H_{\rm edge}+T^{{\dagger}}\left(E+% \text{i}\eta-H_{\rm bulk}\right)^{-1}T.italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ≔ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + roman_Σ = italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T . (I.9)

Let Ensubscript𝐸𝑛E_{n}italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT be an eigenenergy of Hbulksubscript𝐻bulkH_{\rm bulk}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT and |ψnketsubscript𝜓𝑛\ket{\psi_{n}}| start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⟩ be the corresponding eigenstate. Then, we have the spectral decomposition

Σ=T(n|ψnψn|E+iηEn)T.Σsuperscript𝑇subscript𝑛ketsubscript𝜓𝑛brasubscript𝜓𝑛𝐸i𝜂subscript𝐸𝑛𝑇\Sigma=T^{{\dagger}}\left(\sum_{n}\frac{\ket{\psi_{n}}\bra{\psi_{n}}}{E+\text{% i}\eta-E_{n}}\right)T.roman_Σ = italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( ∑ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT divide start_ARG | start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⟩ ⟨ start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG | end_ARG start_ARG italic_E + i italic_η - italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) italic_T . (I.10)

Example

As the simplest example, we consider a single-band model [122]

Hij=t(δi,j+1+δi,j1)(t>0).subscript𝐻𝑖𝑗𝑡subscript𝛿𝑖𝑗1subscript𝛿𝑖𝑗1𝑡0H_{ij}=t\left(\delta_{i,j+1}+\delta_{i,j-1}\right)\quad\left(t>0\right).italic_H start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = italic_t ( italic_δ start_POSTSUBSCRIPT italic_i , italic_j + 1 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_i , italic_j - 1 end_POSTSUBSCRIPT ) ( italic_t > 0 ) . (I.11)

Let |ψn=(ψ1(n),ψ2(n),,ψL(n))ketsubscript𝜓𝑛superscriptsubscript𝜓1𝑛superscriptsubscript𝜓2𝑛superscriptsubscript𝜓𝐿𝑛\ket{\psi_{n}}=(\psi_{1}^{(n)},\psi_{2}^{(n)},\cdots,\psi_{L}^{(n)})| start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⟩ = ( italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT , italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT , ⋯ , italic_ψ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ) be an eigenstate of Hbulksubscript𝐻bulkH_{\rm bulk}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT. Then, the Schrödinger equation H|ψn=En|ψn𝐻ketsubscript𝜓𝑛subscript𝐸𝑛ketsubscript𝜓𝑛H\ket{\psi_{n}}=E_{n}\ket{\psi_{n}}italic_H | start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⟩ = italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT | start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⟩ reads

tψj1(n)+tψj+1(n)=Enψj(n)(j=2,3,,L1),ψ0(n)=ψL+1(n)=0.formulae-sequence𝑡superscriptsubscript𝜓𝑗1𝑛𝑡superscriptsubscript𝜓𝑗1𝑛subscript𝐸𝑛superscriptsubscript𝜓𝑗𝑛𝑗23𝐿1superscriptsubscript𝜓0𝑛superscriptsubscript𝜓𝐿1𝑛0t\psi_{j-1}^{(n)}+t\psi_{j+1}^{(n)}=E_{n}\psi_{j}^{(n)}\quad\left(j=2,3,\cdots% ,L-1\right),\quad\psi_{0}^{(n)}=\psi_{L+1}^{(n)}=0.italic_t italic_ψ start_POSTSUBSCRIPT italic_j - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + italic_t italic_ψ start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_ψ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ( italic_j = 2 , 3 , ⋯ , italic_L - 1 ) , italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = italic_ψ start_POSTSUBSCRIPT italic_L + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = 0 . (I.12)

The eigenenergies and eigenstates are given as

En=2tcoskn,ψj(n)=2L+1sin(knj);knnπL+1(n=1,2,,L).formulae-sequencesubscript𝐸𝑛2𝑡subscript𝑘𝑛formulae-sequencesuperscriptsubscript𝜓𝑗𝑛2𝐿1subscript𝑘𝑛𝑗subscript𝑘𝑛𝑛𝜋𝐿1𝑛12𝐿E_{n}=2t\cos k_{n},\quad\psi_{j}^{(n)}=\sqrt{\frac{2}{L+1}}\sin\left(k_{n}j% \right);\quad k_{n}\coloneqq\frac{n\pi}{L+1}\quad\left(n=1,2,\cdots,L\right).italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 2 italic_t roman_cos italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_ψ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = square-root start_ARG divide start_ARG 2 end_ARG start_ARG italic_L + 1 end_ARG end_ARG roman_sin ( italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_j ) ; italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ≔ divide start_ARG italic_n italic_π end_ARG start_ARG italic_L + 1 end_ARG ( italic_n = 1 , 2 , ⋯ , italic_L ) . (I.13)

Thus, the self-energy of the bulk is obtained as

Σ=t2n=1L|ψj(n)|2E+iηEn=2t2L+1n=1Lsin2knE+iη2tcoskn.Σsuperscript𝑡2superscriptsubscript𝑛1𝐿superscriptsuperscriptsubscript𝜓𝑗𝑛2𝐸i𝜂subscript𝐸𝑛2superscript𝑡2𝐿1superscriptsubscript𝑛1𝐿superscript2subscript𝑘𝑛𝐸i𝜂2𝑡subscript𝑘𝑛\Sigma=t^{2}\sum_{n=1}^{L}\frac{|\psi_{j}^{(n)}|^{2}}{E+\text{i}\eta-E_{n}}=% \frac{2t^{2}}{L+1}\sum_{n=1}^{L}\frac{\sin^{2}k_{n}}{E+\text{i}\eta-2t\cos k_{% n}}.roman_Σ = italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT divide start_ARG | italic_ψ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_E + i italic_η - italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG = divide start_ARG 2 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_L + 1 end_ARG ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG italic_E + i italic_η - 2 italic_t roman_cos italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG . (I.14)

In the semi-infinite limit L𝐿L\to\inftyitalic_L → ∞, we have

Σ=2t2n=1LΔknπsin2knE+iη2tcoskn2t20πdkπsin2kE+iη2tcosk=E2(114t2E2).Σ2superscript𝑡2superscriptsubscript𝑛1𝐿Δsubscript𝑘𝑛𝜋superscript2subscript𝑘𝑛𝐸i𝜂2𝑡subscript𝑘𝑛2superscript𝑡2superscriptsubscript0𝜋𝑑𝑘𝜋superscript2𝑘𝐸i𝜂2𝑡𝑘𝐸2114superscript𝑡2superscript𝐸2\Sigma=2t^{2}\sum_{n=1}^{L}\frac{\Delta k_{n}}{\pi}\frac{\sin^{2}k_{n}}{E+% \text{i}\eta-2t\cos k_{n}}\to 2t^{2}\int_{0}^{\pi}\frac{dk}{\pi}\frac{\sin^{2}% k}{E+\text{i}\eta-2t\cos k}=\frac{E}{2}\left(1-\sqrt{1-\frac{4t^{2}}{E^{2}}}% \right).roman_Σ = 2 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT divide start_ARG roman_Δ italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG italic_π end_ARG divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG italic_E + i italic_η - 2 italic_t roman_cos italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG → 2 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT divide start_ARG italic_d italic_k end_ARG start_ARG italic_π end_ARG divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k end_ARG start_ARG italic_E + i italic_η - 2 italic_t roman_cos italic_k end_ARG = divide start_ARG italic_E end_ARG start_ARG 2 end_ARG ( 1 - square-root start_ARG 1 - divide start_ARG 4 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ) . (I.15)

The imaginary part of ΣΣ\Sigmaroman_Σ appears for E𝐸Eitalic_E within the bulk bandwidth |E|<2t𝐸2𝑡\left|E\right|<2t| italic_E | < 2 italic_t, which reflects the particle exchange between the bulk and the edge.

II Symmetry classification of non-Hermitian self-energy

We study symmetry of effective non-Hermitian Hamiltonians introduced in Sec. I. Specifically, we show that when original Hermitian Hamiltonians belong to the Altland-Zirnbauer symmetry class [3], the corresponding effective non-Hermitian Hamiltonians belong to the Altland-Zirnbauer symmetry class [29]. The similar correspondence was derived, for example, for quadratic Lindbladians [123] and reflection matrices in scattering processes [124]. Because of this correspondence, the topological classification of d𝑑ditalic_d-dimensional Hermitian Hamiltonians coincides with that of (d1)𝑑1\left(d-1\right)( italic_d - 1 )-dimensional effective non-Hermitian Hamiltonians [29, 92, 93, 99, 100].

II.1 Time-reversal symmetry

Suppose that the original Hermitian Hamiltonian H𝐻Hitalic_H respects time-reversal symmetry:

𝒯H𝒯1=H𝒯superscript𝐻superscript𝒯1𝐻\mathcal{T}H^{*}\mathcal{T}^{-1}=Hcaligraphic_T italic_H start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_H (II.1)

with a unitary matrix 𝒯𝒯\mathcal{T}caligraphic_T. Since time reversal acts only on the internal degrees of freedom, we have

𝒯Hbulk𝒯1=Hbulk,𝒯Hedge𝒯1=Hedge,𝒯T𝒯1=T.formulae-sequence𝒯subscriptsuperscript𝐻bulksuperscript𝒯1subscript𝐻bulkformulae-sequence𝒯subscriptsuperscript𝐻edgesuperscript𝒯1subscript𝐻edge𝒯superscript𝑇superscript𝒯1𝑇\mathcal{T}H^{*}_{\rm bulk}\mathcal{T}^{-1}=H_{\rm bulk},\quad\mathcal{T}H^{*}% _{\rm edge}\mathcal{T}^{-1}=H_{\rm edge},\quad\mathcal{T}T^{*}\mathcal{T}^{-1}% =T.caligraphic_T italic_H start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT , caligraphic_T italic_H start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , caligraphic_T italic_T start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_T . (II.2)

Consequently, we have

𝒯ΣT(E)𝒯1=T(E+iηHbulk)1T=Σ(E),𝒯superscriptΣ𝑇𝐸superscript𝒯1superscript𝑇superscript𝐸i𝜂subscript𝐻bulk1𝑇Σ𝐸\mathcal{T}\Sigma^{T}\left(E\right)\mathcal{T}^{-1}=T^{{\dagger}}\left(E+\text% {i}\eta-H_{\rm bulk}\right)^{-1}T=\Sigma\left(E\right),caligraphic_T roman_Σ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_E ) caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T = roman_Σ ( italic_E ) , (II.3)

and hence

𝒯HeffT(E)𝒯1=Heff(E).𝒯superscriptsubscript𝐻eff𝑇𝐸superscript𝒯1subscript𝐻eff𝐸\mathcal{T}H_{\rm eff}^{T}\left(E\right)\mathcal{T}^{-1}=H_{\rm eff}\left(E% \right).caligraphic_T italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_E ) caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) . (II.4)

Thus, for arbitrary E𝐸Eitalic_E, the effective non-Hermitian Hamiltonian Heffsubscript𝐻effH_{\rm eff}italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT respects time-reversal symmetry [29]. Notably, Heffsubscript𝐻effH_{\rm eff}italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT is not invariant under time reversal (i.e., 𝒯Heff𝒯1Heff𝒯superscriptsubscript𝐻effsuperscript𝒯1subscript𝐻eff\mathcal{T}H_{\rm eff}^{*}\mathcal{T}^{-1}\neq H_{\rm eff}caligraphic_T italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ≠ italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT).

II.2 Particle-hole symmetry

Suppose that the original Hermitian Hamiltonian H𝐻Hitalic_H respects particle-hole symmetry:

𝒞H𝒞1=H𝒞superscript𝐻superscript𝒞1𝐻\mathcal{C}H^{*}\mathcal{C}^{-1}=-Hcaligraphic_C italic_H start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H (II.5)

with a unitary matrix 𝒞𝒞\mathcal{C}caligraphic_C. Since particle-hole transformation acts only on the internal degrees of freedom, we have

𝒞Hbulk𝒞1=Hbulk,𝒞Hedge𝒞1=Hedge,𝒞T𝒞1=T.formulae-sequence𝒞subscriptsuperscript𝐻bulksuperscript𝒞1subscript𝐻bulkformulae-sequence𝒞subscriptsuperscript𝐻edgesuperscript𝒞1subscript𝐻edge𝒞superscript𝑇superscript𝒞1𝑇\mathcal{C}H^{*}_{\rm bulk}\mathcal{C}^{-1}=-H_{\rm bulk},\quad\mathcal{C}H^{*% }_{\rm edge}\mathcal{C}^{-1}=-H_{\rm edge},\quad\mathcal{C}T^{*}\mathcal{C}^{-% 1}=-T.caligraphic_C italic_H start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT , caligraphic_C italic_H start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , caligraphic_C italic_T start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_T . (II.6)

Consequently, we have

𝒞Σ(E)𝒞1=T(Eiη+Hbulk)1T=Σ(E),𝒞superscriptΣ𝐸superscript𝒞1superscript𝑇superscript𝐸i𝜂subscript𝐻bulk1𝑇Σ𝐸\mathcal{C}\Sigma^{*}\left(E\right)\mathcal{C}^{-1}=T^{{\dagger}}\left(E-\text% {i}\eta+H_{\rm bulk}\right)^{-1}T=-\Sigma\left(-E\right),caligraphic_C roman_Σ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_E ) caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E - i italic_η + italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T = - roman_Σ ( - italic_E ) , (II.7)

and hence

𝒞Heff(E)𝒞1=Heff(E).𝒞superscriptsubscript𝐻eff𝐸superscript𝒞1subscript𝐻eff𝐸\mathcal{C}H_{\rm eff}^{*}\left(E\right)\mathcal{C}^{-1}=-H_{\rm eff}\left(-E% \right).caligraphic_C italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_E ) caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( - italic_E ) . (II.8)

Thus, at the particle-hole-symmetric point E=0𝐸0E=0italic_E = 0, the effective non-Hermitian Hamiltonian Heff(E=0)subscript𝐻eff𝐸0H_{\rm eff}\left(E=0\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E = 0 ) respects particle-hole symmetry [29]. Notably, Heffsubscript𝐻effH_{\rm eff}italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT does not respect particle-hole symmetry even at E=0𝐸0E=0italic_E = 0 [i.e., 𝒞HeffT(E=0)𝒞1Heff(E=0)𝒞superscriptsubscript𝐻eff𝑇𝐸0superscript𝒞1subscript𝐻eff𝐸0\mathcal{C}H_{\rm eff}^{T}\left(E=0\right)\mathcal{C}^{-1}\neq-H_{\rm eff}% \left(E=0\right)caligraphic_C italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_E = 0 ) caligraphic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ≠ - italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E = 0 )].

II.3 Chiral symmetry

Suppose that the original Hermitian Hamiltonian H𝐻Hitalic_H respects chiral symmetry:

ΓHΓ1=HΓ𝐻superscriptΓ1𝐻\Gamma H\Gamma^{-1}=-Hroman_Γ italic_H roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H (II.9)

with a unitary matrix ΓΓ\Gammaroman_Γ. Since chiral transformation acts only on the internal degrees of freedom, we have

ΓHbulkΓ1=Hbulk,ΓHedgeΓ1=Hedge,ΓTΓ1=T.formulae-sequenceΓsubscript𝐻bulksuperscriptΓ1subscript𝐻bulkformulae-sequenceΓsubscript𝐻edgesuperscriptΓ1subscript𝐻edgeΓ𝑇superscriptΓ1𝑇\Gamma H_{\rm bulk}\Gamma^{-1}=-H_{\rm bulk},\quad\Gamma H_{\rm edge}\Gamma^{-% 1}=-H_{\rm edge},\quad\Gamma T\Gamma^{-1}=-T.roman_Γ italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT , roman_Γ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , roman_Γ italic_T roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_T . (II.10)

Consequently, we have

ΓΣ(E)Γ1=T(Eiη+Hbulk)1T=Σ(E),ΓsuperscriptΣ𝐸superscriptΓ1superscript𝑇superscript𝐸i𝜂subscript𝐻bulk1𝑇Σ𝐸\Gamma\Sigma^{{\dagger}}\left(E\right)\Gamma^{-1}=T^{{\dagger}}\left(E-\text{i% }\eta+H_{\rm bulk}\right)^{-1}T=-\Sigma\left(-E\right),roman_Γ roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E - i italic_η + italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_T = - roman_Σ ( - italic_E ) , (II.11)

and hence

ΓHeff(E)Γ1=Heff(E).Γsuperscriptsubscript𝐻eff𝐸superscriptΓ1subscript𝐻eff𝐸\Gamma H_{\rm eff}^{{\dagger}}\left(E\right)\Gamma^{-1}=-H_{\rm eff}\left(-E% \right).roman_Γ italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( - italic_E ) . (II.12)

Thus, at the chiral-symmetric point E=0𝐸0E=0italic_E = 0, the effective non-Hermitian Hamiltonian Heff(E=0)subscript𝐻eff𝐸0H_{\rm eff}\left(E=0\right)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E = 0 ) respects chiral symmetry [29]. Notably, Heffsubscript𝐻effH_{\rm eff}italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT does not respect sublattice symmetry even at E=0𝐸0E=0italic_E = 0 [i.e., ΓHeff(E=0)Γ1Heff(E=0)Γsubscript𝐻eff𝐸0superscriptΓ1subscript𝐻eff𝐸0\Gamma H_{\rm eff}\left(E=0\right)\Gamma^{-1}\neq-H_{\rm eff}\left(E=0\right)roman_Γ italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E = 0 ) roman_Γ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ≠ - italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E = 0 )].

III Analytical derivation of the self-energy

We analytically derive the self-energy between the bulk and boundaries for the Su-Schrieffer-Heeger (SSH) model and a Chern insulator.

III.1 Su-Schrieffer-Heeger model

We consider the SSH model [109]

Hbulk(k)=(v+tcos(k))σx+(tsin(k))σy,subscript𝐻bulk𝑘𝑣𝑡𝑘subscript𝜎𝑥𝑡𝑘subscript𝜎𝑦\displaystyle H_{\mathrm{bulk}}\left(k\right)=\left(v+t\cos{k}\right)\sigma_{x% }+\left(t\sin{k}\right)\sigma_{y},italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_k ) = ( italic_v + italic_t roman_cos ( start_ARG italic_k end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin ( start_ARG italic_k end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , (III.1)

where σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s (i=x,y,z𝑖𝑥𝑦𝑧i=x,y,zitalic_i = italic_x , italic_y , italic_z) are Pauli matrices, and v>0𝑣0v>0italic_v > 0 and t>0𝑡0t>0italic_t > 0 are the intracell and intercell hopping amplitudes, respectively. Below, we assume the open boundary conditions. In the topologically nontrivial phase r<1𝑟1r<1italic_r < 1 (rv/t𝑟𝑣𝑡r\coloneqq v/titalic_r ≔ italic_v / italic_t), a pair of zero-energy states appears at the boundaries. We assume that the total number of sites is odd, 2L+12𝐿12L+12 italic_L + 1, for the sake of simplicity, for which eigenstates are exactly obtained under the open boundary conditions (see, for example, Refs. [121, 125]). Let |ψn=(ψ1(n),ψ2(n),,ψ2L+1(n))ketsubscript𝜓𝑛superscriptsubscript𝜓1𝑛superscriptsubscript𝜓2𝑛superscriptsubscript𝜓2𝐿1𝑛\ket{\psi_{n}}=(\psi_{1}^{(n)},\psi_{2}^{(n)},\cdots,\psi_{2L+1}^{(n)})| start_ARG italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ⟩ = ( italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT , italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT , ⋯ , italic_ψ start_POSTSUBSCRIPT 2 italic_L + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ) be an eigenstate. The bulk states host the eigenenergy

En=±tr2+1+2rcos(kn)(knnπL+1;n=1,,L),subscript𝐸𝑛plus-or-minus𝑡superscript𝑟212𝑟subscript𝑘𝑛formulae-sequencesubscript𝑘𝑛𝑛𝜋𝐿1𝑛1𝐿\displaystyle E_{n}=\pm t\sqrt{r^{2}+1+2r\cos{k_{n}}}\quad\left(k_{n}\coloneqq% \frac{n\pi}{L+1};~{}n=1,\cdots,L\right),italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = ± italic_t square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) end_ARG ( italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ≔ divide start_ARG italic_n italic_π end_ARG start_ARG italic_L + 1 end_ARG ; italic_n = 1 , ⋯ , italic_L ) , (III.2)

and the corresponding bulk eigenstates are obtained as

ψ2j1(n)=rsin(knj)+sin(kn(j1))(L+1)(1+2rcos(kn)+r2)(j=1,,L+1),ψ2j(n)=±sin((knj))L+1(j=1,,L).formulae-sequencesuperscriptsubscript𝜓2𝑗1𝑛𝑟subscript𝑘𝑛𝑗subscript𝑘𝑛𝑗1𝐿112𝑟subscript𝑘𝑛superscript𝑟2𝑗1𝐿1superscriptsubscript𝜓2𝑗𝑛plus-or-minussubscript𝑘𝑛𝑗𝐿1𝑗1𝐿\displaystyle\psi_{2j-1}^{(n)}=\frac{r\sin\left(k_{n}j\right)+\sin\left(k_{n}(% j-1)\right)}{\sqrt{(L+1)\left(1+2r\cos{k_{n}}+r^{2}\right)}}\hskip 14.22636pt% \left(j=1,\cdots,L+1\right),\quad\quad~{}\psi_{2j}^{(n)}=\pm\frac{\sin{\left(k% _{n}j\right)}}{L+1}\hskip 14.22636pt\left(j=1,\cdots,L\right).italic_ψ start_POSTSUBSCRIPT 2 italic_j - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = divide start_ARG italic_r roman_sin ( italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_j ) + roman_sin ( italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_j - 1 ) ) end_ARG start_ARG square-root start_ARG ( italic_L + 1 ) ( 1 + 2 italic_r roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG end_ARG ( italic_j = 1 , ⋯ , italic_L + 1 ) , italic_ψ start_POSTSUBSCRIPT 2 italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = ± divide start_ARG roman_sin ( start_ARG ( italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_j ) end_ARG ) end_ARG start_ARG italic_L + 1 end_ARG ( italic_j = 1 , ⋯ , italic_L ) . (III.3)

On the other hand, the zero-energy eigenstate localized at the left edge is given as

ψ2j=0(j=1,,L),ψ2j+1subscript𝜓2𝑗0𝑗1𝐿subscript𝜓2𝑗1\displaystyle\psi_{2j}=0\hskip 14.22636pt\left(j=1,\cdots,L\right),\quad\quad% \psi_{2j+1}italic_ψ start_POSTSUBSCRIPT 2 italic_j end_POSTSUBSCRIPT = 0 ( italic_j = 1 , ⋯ , italic_L ) , italic_ψ start_POSTSUBSCRIPT 2 italic_j + 1 end_POSTSUBSCRIPT =(r)j1r21r2(L+1)(j=0,,L)absentsuperscript𝑟𝑗1superscript𝑟21superscript𝑟2𝐿1𝑗0𝐿\displaystyle=(-r)^{j}\sqrt{\frac{1-r^{2}}{1-r^{2(L+1)}}}\hskip 14.22636pt% \left(j=0,\cdots,L\right)= ( - italic_r ) start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 1 - italic_r start_POSTSUPERSCRIPT 2 ( italic_L + 1 ) end_POSTSUPERSCRIPT end_ARG end_ARG ( italic_j = 0 , ⋯ , italic_L ) (III.4)

for r<1𝑟1r<1italic_r < 1.

Now, we compute the self-energy Σ=Σedge+ΣbulkΣsubscriptΣedgesubscriptΣbulk\Sigma=\Sigma_{\mathrm{edge}}+\Sigma_{\mathrm{bulk}}roman_Σ = roman_Σ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT. The contribution ΣedgesubscriptΣedge\Sigma_{\mathrm{edge}}roman_Σ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT from the edge state is obtained as

Σedge=t2(1r2)1r2(L+1)1E+iησ0σz2t2(1r2)E+iησ0σz2subscriptΣedgesuperscript𝑡21superscript𝑟21superscript𝑟2𝐿11𝐸i𝜂subscript𝜎0subscript𝜎𝑧2superscript𝑡21superscript𝑟2𝐸i𝜂subscript𝜎0subscript𝜎𝑧2\displaystyle\Sigma_{\mathrm{edge}}=\frac{t^{2}(1-r^{2})}{1-{r}^{2(L+1)}}\frac% {1}{E+\text{i}\eta}\frac{\sigma_{0}-\sigma_{z}}{2}\rightarrow\frac{t^{2}(1-r^{% 2})}{E+\text{i}\eta}\frac{\sigma_{0}-\sigma_{z}}{2}roman_Σ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT = divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 1 - italic_r start_POSTSUPERSCRIPT 2 ( italic_L + 1 ) end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η end_ARG divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG → divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_E + i italic_η end_ARG divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG (III.5)

in the limit L𝐿L\rightarrow\inftyitalic_L → ∞ for r<1𝑟1r<1italic_r < 1. Here, σ0subscript𝜎0\sigma_{0}italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the 2×2222\times 22 × 2 identity matrix. Next, from Eqs. (III.2) and (III.3), the self-energy ΣbulksubscriptΣbulk\Sigma_{\mathrm{bulk}}roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT from the bulk states is obtained as

Σbulk=subscriptΣbulkabsent\displaystyle\Sigma_{\mathrm{bulk}}=roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT = t2L+1n=1Lr2sin2knr2+1+2rcos(kn)(1E+iηtr2+1+2rcos(kn)+1E+iη+tr2+1+2rcos(kn))σ0σz2superscript𝑡2𝐿1superscriptsubscript𝑛1𝐿superscript𝑟2superscript2subscript𝑘𝑛superscript𝑟212𝑟subscript𝑘𝑛1𝐸i𝜂𝑡superscript𝑟212𝑟subscript𝑘𝑛1𝐸i𝜂𝑡superscript𝑟212𝑟subscript𝑘𝑛subscript𝜎0subscript𝜎𝑧2\displaystyle\frac{t^{2}}{L+1}\sum\limits_{n=1}^{L}\frac{r^{2}\sin^{2}{k_{n}}}% {r^{2}+1+2r\cos{k_{n}}}\Biggl{(}\frac{1}{E+\text{i}\eta-t\sqrt{r^{2}+1+2r\cos{% k_{n}}}}+\frac{1}{E+\text{i}\eta+t\sqrt{r^{2}+1+2r\cos{k_{n}}}}\Biggr{)}\frac{% \sigma_{0}-\sigma_{z}}{2}divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_L + 1 end_ARG ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT divide start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) end_ARG ( divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η - italic_t square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) end_ARG end_ARG + divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η + italic_t square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG ) end_ARG end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG
\displaystyle\rightarrow\, t20πdkπr2sin2kr2+1+2rcos(k)(1E+iηtr2+1+2rcos(k)+1E+iη+tr2+1+2rcos(k))σ0σz2superscript𝑡2superscriptsubscript0𝜋𝑑𝑘𝜋superscript𝑟2superscript2𝑘superscript𝑟212𝑟𝑘1𝐸i𝜂𝑡superscript𝑟212𝑟𝑘1𝐸i𝜂𝑡superscript𝑟212𝑟𝑘subscript𝜎0subscript𝜎𝑧2\displaystyle t^{2}\int_{0}^{\pi}\frac{dk}{\pi}\frac{r^{2}\sin^{2}{k}}{r^{2}+1% +2r\cos{k}}\Biggl{(}\frac{1}{E+\text{i}\eta-t\sqrt{r^{2}+1+2r\cos{k}}}+\frac{1% }{E+\text{i}\eta+t\sqrt{r^{2}+1+2r\cos{k}}}\Biggr{)}\frac{\sigma_{0}-\sigma_{z% }}{2}italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_π end_POSTSUPERSCRIPT divide start_ARG italic_d italic_k end_ARG start_ARG italic_π end_ARG divide start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k end_ARG ) end_ARG ( divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η - italic_t square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k end_ARG ) end_ARG end_ARG + divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η + italic_t square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 + 2 italic_r roman_cos ( start_ARG italic_k end_ARG ) end_ARG end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG (III.6)

for L𝐿L\rightarrow\inftyitalic_L → ∞. By introducing the complex variable z=eik𝑧superscript𝑒i𝑘z=e^{\text{i}k}italic_z = italic_e start_POSTSUPERSCRIPT i italic_k end_POSTSUPERSCRIPT, the self-energy further leads to the contour integral

ΣbulkE+iη2σ0σz212πiC𝑑zf(z),subscriptΣbulk𝐸i𝜂2subscript𝜎0subscript𝜎𝑧212𝜋isubscriptcontour-integral𝐶differential-d𝑧𝑓𝑧\displaystyle\Sigma_{\mathrm{bulk}}\rightarrow\frac{E+\text{i}\eta}{2}\frac{% \sigma_{0}-\sigma_{z}}{2}\frac{1}{2\pi\text{i}}\oint_{C}dz~{}f(z),roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT → divide start_ARG italic_E + i italic_η end_ARG start_ARG 2 end_ARG divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG divide start_ARG 1 end_ARG start_ARG 2 italic_π i end_ARG ∮ start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT italic_d italic_z italic_f ( italic_z ) , (III.7)

where C𝐶Citalic_C is the unit circle in the complex plane, and f(z)𝑓𝑧f(z)italic_f ( italic_z ) is defined as

f(z)1z(z21)2z2+(r+r1)z+11z2+2uz+1,2u(r+1r)1r(E+iηt)2.formulae-sequence𝑓𝑧1𝑧superscriptsuperscript𝑧212superscript𝑧2𝑟superscript𝑟1𝑧11superscript𝑧22𝑢𝑧12𝑢𝑟1𝑟1𝑟superscript𝐸i𝜂𝑡2\displaystyle f(z)\coloneqq\frac{1}{z}\frac{(z^{2}-1)^{2}}{z^{2}+\left(r+r^{-1% }\right)z+1}\frac{1}{z^{2}+2uz+1},\quad 2u\coloneqq\left(r+\frac{1}{r}\right)-% \frac{1}{r}\left(\frac{E+\text{i}\eta}{t}\right)^{2}.italic_f ( italic_z ) ≔ divide start_ARG 1 end_ARG start_ARG italic_z end_ARG divide start_ARG ( italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_r + italic_r start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) italic_z + 1 end_ARG divide start_ARG 1 end_ARG start_ARG italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_u italic_z + 1 end_ARG , 2 italic_u ≔ ( italic_r + divide start_ARG 1 end_ARG start_ARG italic_r end_ARG ) - divide start_ARG 1 end_ARG start_ARG italic_r end_ARG ( divide start_ARG italic_E + i italic_η end_ARG start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (III.8)

The poles of f(z)𝑓𝑧f(z)italic_f ( italic_z ) are

z=0,r,1r,z±,(z±u±u21),𝑧0𝑟1𝑟subscript𝑧plus-or-minussubscript𝑧plus-or-minusplus-or-minus𝑢superscript𝑢21z=0,-r,-\frac{1}{r},z_{\pm},\quad\left(z_{\pm}\coloneqq-u\pm\sqrt{u^{2}-1}% \right),italic_z = 0 , - italic_r , - divide start_ARG 1 end_ARG start_ARG italic_r end_ARG , italic_z start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT , ( italic_z start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ≔ - italic_u ± square-root start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 end_ARG ) , (III.9)

whose residues are calculated as

Res[z=0]=1,Res[z=r]=t2(1r2)(E+iη)2,Res[z=1r]=t2(1r2)(E+iη)2,Res[z=z±]=±2t2ru21(E+iη)2.formulae-sequenceresidue𝑧01formulae-sequenceresidue𝑧𝑟superscript𝑡21superscript𝑟2superscript𝐸i𝜂2formulae-sequenceresidue𝑧1𝑟superscript𝑡21superscript𝑟2superscript𝐸i𝜂2residue𝑧subscript𝑧plus-or-minusplus-or-minus2superscript𝑡2𝑟superscript𝑢21superscript𝐸i𝜂2\displaystyle\Res[z=0]=1,~{}\Res[z=-r]=-\frac{t^{2}(1-r^{2})}{\left({E+\text{i% }\eta}\right)^{2}},~{}\Res[z=-\frac{1}{r}]=\frac{t^{2}(1-r^{2})}{\left({E+% \text{i}\eta}\right)^{2}},~{}\Res[z=z_{\pm}]=\pm\frac{2t^{2}r\sqrt{u^{2}-1}}{% \left({E+\text{i}\eta}\right)^{2}}.roman_Res [ italic_z = 0 ] = 1 , roman_Res [ italic_z = - italic_r ] = - divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , roman_Res [ italic_z = - divide start_ARG 1 end_ARG start_ARG italic_r end_ARG ] = divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , roman_Res [ italic_z = italic_z start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ] = ± divide start_ARG 2 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r square-root start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 end_ARG end_ARG start_ARG ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (III.10)

Depending on r𝑟ritalic_r, each pole is either located inside or outside the unit circle C𝐶Citalic_C, as summarized in the following:

  1. (i)

    For r>1𝑟1r>1italic_r > 1 and E2<t2(r2+1)superscript𝐸2superscript𝑡2superscript𝑟21E^{2}<t^{2}(r^{2}+1)italic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ), only the three poles z=0,r1,z+𝑧0superscript𝑟1subscript𝑧z=0,-r^{-1},z_{+}italic_z = 0 , - italic_r start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , italic_z start_POSTSUBSCRIPT + end_POSTSUBSCRIPT are located inside the unit circle C𝐶Citalic_C, leading to

    Σbulk=(E+iη2+12t2(1r2)E+iη+12[t2(r2+1)(E+iη)2]24t4r2E+iη)σ0σz2.subscriptΣbulk𝐸i𝜂212superscript𝑡21superscript𝑟2𝐸i𝜂12superscriptdelimited-[]superscript𝑡2superscript𝑟21superscript𝐸i𝜂224superscript𝑡4superscript𝑟2𝐸i𝜂subscript𝜎0subscript𝜎𝑧2\displaystyle\Sigma_{\mathrm{bulk}}=\left(\frac{E+\text{i}\eta}{2}+\frac{1}{2}% \frac{t^{2}(1-r^{2})}{E+\text{i}\eta}+\frac{1}{2}\frac{\sqrt{\left[t^{2}(r^{2}% +1)-\left(E+\text{i}\eta\right)^{2}\right]^{2}-4t^{4}r^{2}}}{E+\text{i}\eta}% \right)\frac{\sigma_{0}-\sigma_{z}}{2}.roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT = ( divide start_ARG italic_E + i italic_η end_ARG start_ARG 2 end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_E + i italic_η end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG square-root start_ARG [ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) - ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_t start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_E + i italic_η end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG . (III.11)
  2. (ii)

    For r>1𝑟1r>1italic_r > 1 and E2>t2(r2+1)superscript𝐸2superscript𝑡2superscript𝑟21E^{2}>t^{2}(r^{2}+1)italic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT > italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ), only the three poles z=0,r1,z𝑧0superscript𝑟1subscript𝑧z=0,-r^{-1},z_{-}italic_z = 0 , - italic_r start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , italic_z start_POSTSUBSCRIPT - end_POSTSUBSCRIPT are located inside the unit circle C𝐶Citalic_C, leading to

    Σbulk=(E+iη2+12t2(1r2)E+iη12[t2(r2+1)(E+iη)2]24t4r2E+iη)σ0σz2.subscriptΣbulk𝐸i𝜂212superscript𝑡21superscript𝑟2𝐸i𝜂12superscriptdelimited-[]superscript𝑡2superscript𝑟21superscript𝐸i𝜂224superscript𝑡4superscript𝑟2𝐸i𝜂subscript𝜎0subscript𝜎𝑧2\displaystyle\Sigma_{\mathrm{bulk}}=\left(\frac{E+\text{i}\eta}{2}+\frac{1}{2}% \frac{t^{2}(1-r^{2})}{E+\text{i}\eta}-\frac{1}{2}\frac{\sqrt{\left[t^{2}(r^{2}% +1)-\left(E+\text{i}\eta\right)^{2}\right]^{2}-4t^{4}r^{2}}}{E+\text{i}\eta}% \right)\frac{\sigma_{0}-\sigma_{z}}{2}.roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT = ( divide start_ARG italic_E + i italic_η end_ARG start_ARG 2 end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_E + i italic_η end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG square-root start_ARG [ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) - ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_t start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_E + i italic_η end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG . (III.12)
  3. (iii)

    For r<1𝑟1r<1italic_r < 1 and E2<t2(r2+1)superscript𝐸2superscript𝑡2superscript𝑟21E^{2}<t^{2}(r^{2}+1)italic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ), only the three poles z=0,r,z+𝑧0𝑟subscript𝑧z=0,-r,z_{+}italic_z = 0 , - italic_r , italic_z start_POSTSUBSCRIPT + end_POSTSUBSCRIPT are located inside the unit circle C𝐶Citalic_C, leading to

    Σbulk=(E+iη212t2(1r2)E+iη+12[t2(r2+1)(E+iη)2]24t4r2E+iη)σ0σz2.subscriptΣbulk𝐸i𝜂212superscript𝑡21superscript𝑟2𝐸i𝜂12superscriptdelimited-[]superscript𝑡2superscript𝑟21superscript𝐸i𝜂224superscript𝑡4superscript𝑟2𝐸i𝜂subscript𝜎0subscript𝜎𝑧2\displaystyle\Sigma_{\mathrm{bulk}}=\left(\frac{E+\text{i}\eta}{2}-\frac{1}{2}% \frac{t^{2}(1-r^{2})}{E+\text{i}\eta}+\frac{1}{2}\frac{\sqrt{\left[t^{2}(r^{2}% +1)-\left(E+\text{i}\eta\right)^{2}\right]^{2}-4t^{4}r^{2}}}{E+\text{i}\eta}% \right)\frac{\sigma_{0}-\sigma_{z}}{2}.roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT = ( divide start_ARG italic_E + i italic_η end_ARG start_ARG 2 end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_E + i italic_η end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG square-root start_ARG [ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) - ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_t start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_E + i italic_η end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG . (III.13)
  4. (iv)

    For r<1𝑟1r<1italic_r < 1 and E2>t2(r2+1)superscript𝐸2superscript𝑡2superscript𝑟21E^{2}>t^{2}(r^{2}+1)italic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT > italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ), only the three poles z=0,r,z𝑧0𝑟subscript𝑧z=0,-r,z_{-}italic_z = 0 , - italic_r , italic_z start_POSTSUBSCRIPT - end_POSTSUBSCRIPT are located inside the unit circle C𝐶Citalic_C, leading to

    Σbulk=(E+iη212t2(1r2)E+iη12[t2(r2+1)(E+iη)2]24t4r2E+iη)σ0σz2.subscriptΣbulk𝐸i𝜂212superscript𝑡21superscript𝑟2𝐸i𝜂12superscriptdelimited-[]superscript𝑡2superscript𝑟21superscript𝐸i𝜂224superscript𝑡4superscript𝑟2𝐸i𝜂subscript𝜎0subscript𝜎𝑧2\displaystyle\Sigma_{\mathrm{bulk}}=\left(\frac{E+\text{i}\eta}{2}-\frac{1}{2}% \frac{t^{2}(1-r^{2})}{E+\text{i}\eta}-\frac{1}{2}\frac{\sqrt{\left[t^{2}(r^{2}% +1)-\left(E+\text{i}\eta\right)^{2}\right]^{2}-4t^{4}r^{2}}}{E+\text{i}\eta}% \right)\frac{\sigma_{0}-\sigma_{z}}{2}.roman_Σ start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT = ( divide start_ARG italic_E + i italic_η end_ARG start_ARG 2 end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_E + i italic_η end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG square-root start_ARG [ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) - ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_t start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_E + i italic_η end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG . (III.14)

Combining both contributions from the bulk and edge, we obtain the total self-energy as

Σ(E)=(E+iη2+12t2(1r2)E+iη+12sgn[t2(r2+1)E2][t2(r2+1)(E+iη)2]24t4r2E+iη)σ0σz2.Σ𝐸𝐸i𝜂212superscript𝑡21superscript𝑟2𝐸i𝜂12sgndelimited-[]superscript𝑡2superscript𝑟21superscript𝐸2superscriptdelimited-[]superscript𝑡2superscript𝑟21superscript𝐸i𝜂224superscript𝑡4superscript𝑟2𝐸i𝜂subscript𝜎0subscript𝜎𝑧2\displaystyle\Sigma(E)=\left(\frac{E+\text{i}\eta}{2}+\frac{1}{2}\frac{t^{2}(1% -r^{2})}{E+\text{i}\eta}+\frac{1}{2}\mathrm{sgn}\left[t^{2}(r^{2}+1)-E^{2}% \right]\frac{\sqrt{\left[t^{2}(r^{2}+1)-\left(E+\text{i}\eta\right)^{2}\right]% ^{2}-4t^{4}r^{2}}}{E+\text{i}\eta}\right)\frac{\sigma_{0}-\sigma_{z}}{2}.roman_Σ ( italic_E ) = ( divide start_ARG italic_E + i italic_η end_ARG start_ARG 2 end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_E + i italic_η end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_sgn [ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) - italic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] divide start_ARG square-root start_ARG [ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 1 ) - ( italic_E + i italic_η ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_t start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_E + i italic_η end_ARG ) divide start_ARG italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG . (III.15)

The imaginary part of the self-energy appears for E𝐸Eitalic_E within the bulk bandwidth |vt|<|E|<|v+t|𝑣𝑡𝐸𝑣𝑡\left|v-t\right|<\left|E\right|<\left|v+t\right|| italic_v - italic_t | < | italic_E | < | italic_v + italic_t | or for zero energy E=0𝐸0E=0italic_E = 0. At the critical point r=1𝑟1r=1italic_r = 1, Eq. (III.15) reduces to Eq. (I.15). Additionally, at zero energy E=0𝐸0E=0italic_E = 0, the self-energy reduces to

Σ(0)={0(r>1);iπt2(1r2)δ(0)(σ0σz)/2(r<1).Σ0cases0𝑟1i𝜋superscript𝑡21superscript𝑟2𝛿0subscript𝜎0subscript𝜎𝑧2𝑟1\displaystyle\Sigma(0)=\begin{cases}0\hskip 14.22636pt&(r>1);\\ -\text{i}\pi t^{2}(1-r^{2})\delta(0)\left(\sigma_{0}-\sigma_{z}\right)/2\hskip 1% 4.22636pt&(r<1).\end{cases}roman_Σ ( 0 ) = { start_ROW start_CELL 0 end_CELL start_CELL ( italic_r > 1 ) ; end_CELL end_ROW start_ROW start_CELL - i italic_π italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_δ ( 0 ) ( italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) / 2 end_CELL start_CELL ( italic_r < 1 ) . end_CELL end_ROW (III.16)

Thus, in the topologically nontrivial phase (r<1𝑟1r<1italic_r < 1), the imaginary part of the self-energy diverges.

III.2 Chern insulator

We consider the Chern insulator [111]

Hbulk(kx,ky)=(tsin(kx))σx+(tsin(ky))σy+(m+tcos(kx)+tcos(ky))σzsubscript𝐻bulksubscript𝑘𝑥subscript𝑘𝑦𝑡subscript𝑘𝑥subscript𝜎𝑥𝑡subscript𝑘𝑦subscript𝜎𝑦𝑚𝑡subscript𝑘𝑥𝑡subscript𝑘𝑦subscript𝜎𝑧\displaystyle H_{\mathrm{bulk}}\left(k_{x},k_{y}\right)=\left(t\sin{k_{x}}% \right)\sigma_{x}+\left(t\sin{k_{y}}\right)\sigma_{y}+\left(m+t\cos{k_{x}}+t% \cos{k_{y}}\right)\sigma_{z}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + ( italic_m + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT (III.17)

with Pauli matrices σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s (i=x,y,z𝑖𝑥𝑦𝑧i=x,y,zitalic_i = italic_x , italic_y , italic_z). Here, t>0𝑡0t>0italic_t > 0, denotes the hopping amplitude, and m𝑚m\in\mathbb{R}italic_m ∈ blackboard_R denotes the onsite potential. We assume the square lattice under the open boundary conditions along the x𝑥xitalic_x direction and the periodic boundary conditions along the y𝑦yitalic_y direction. Defining R(ky)𝑅subscript𝑘𝑦R(k_{y})italic_R ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) and S𝑆Sitalic_S as

R(ky)(tsin(ky))σy+(m+tcos(ky))σz,St2(σziσx),formulae-sequence𝑅subscript𝑘𝑦𝑡subscript𝑘𝑦subscript𝜎𝑦𝑚𝑡subscript𝑘𝑦subscript𝜎𝑧𝑆𝑡2subscript𝜎𝑧isubscript𝜎𝑥\displaystyle R(k_{y})\coloneqq\left(t\sin{k_{y}}\right)\sigma_{y}+\left(m+t% \cos{k_{y}}\right)\sigma_{z},\quad S\coloneqq\frac{t}{2}\left(\sigma_{z}-\text% {i}\sigma_{x}\right),italic_R ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ≔ ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + ( italic_m + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_S ≔ divide start_ARG italic_t end_ARG start_ARG 2 end_ARG ( italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT - i italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ) , (III.18)

we express the Hamiltonian as

Hbulk(ky)=x=1L|xx|R(ky)+x=1L1[|xx+1|S+|x+1x|S].subscript𝐻bulksubscript𝑘𝑦superscriptsubscript𝑥1𝐿tensor-productket𝑥bra𝑥𝑅subscript𝑘𝑦superscriptsubscript𝑥1𝐿1delimited-[]tensor-productket𝑥bra𝑥1𝑆tensor-productket𝑥1bra𝑥superscript𝑆\displaystyle H_{\mathrm{bulk}}(k_{y})=\sum\limits_{x=1}^{L}\ket{x}\bra{x}% \otimes R(k_{y})+\sum\limits_{x=1}^{L-1}\Bigg{[}\ket{x}\bra{x+1}\otimes S+\ket% {x+1}\bra{x}\otimes S^{\dagger}\Bigg{]}.italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_x = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT | start_ARG italic_x end_ARG ⟩ ⟨ start_ARG italic_x end_ARG | ⊗ italic_R ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) + ∑ start_POSTSUBSCRIPT italic_x = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L - 1 end_POSTSUPERSCRIPT [ | start_ARG italic_x end_ARG ⟩ ⟨ start_ARG italic_x + 1 end_ARG | ⊗ italic_S + | start_ARG italic_x + 1 end_ARG ⟩ ⟨ start_ARG italic_x end_ARG | ⊗ italic_S start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ] . (III.19)

In the topologically nontrivial phase |m|<2t𝑚2𝑡\absolutevalue{m}<2t| start_ARG italic_m end_ARG | < 2 italic_t, a chiral edge state appears around the boundaries x=1𝑥1x=1italic_x = 1 and x=L𝑥𝐿x=Litalic_x = italic_L. We use the ansatz

|ψedge(ky)=x=1Lβx(ky)|x|u(ky)ketsubscript𝜓edgesubscript𝑘𝑦superscriptsubscript𝑥1𝐿tensor-productsuperscript𝛽𝑥subscript𝑘𝑦ket𝑥ket𝑢subscript𝑘𝑦\displaystyle\ket{\psi_{\mathrm{edge}}\left(k_{y}\right)}=\sum\limits_{x=1}^{L% }\beta^{x}\left(k_{y}\right)\ket{x}\otimes\ket{u\left(k_{y}\right)}| start_ARG italic_ψ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ = ∑ start_POSTSUBSCRIPT italic_x = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT italic_β start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | start_ARG italic_x end_ARG ⟩ ⊗ | start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ (III.20)

for the chiral edge state localized around x=1𝑥1x=1italic_x = 1 with a kysubscript𝑘𝑦k_{y}italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT-dependent two-component vector |u(ky)ket𝑢subscript𝑘𝑦\ket{u\left(k_{y}\right)}| start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩. Here, 1/log|β(ky)|1𝛽subscript𝑘𝑦-1/\log\left|\beta\left(k_{y}\right)\right|- 1 / roman_log | italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | gives a kysubscript𝑘𝑦k_{y}italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT-dependent localization length of the chiral edge state along the x𝑥xitalic_x direction. Then, we reduce the eigenvalue equation Hbulk(ky)|ψedge(ky)=λedge(ky)|ψedge(ky)subscript𝐻bulksubscript𝑘𝑦ketsubscript𝜓edgesubscript𝑘𝑦subscript𝜆edgesubscript𝑘𝑦ketsubscript𝜓edgesubscript𝑘𝑦H_{\mathrm{bulk}}(k_{y})\ket{\psi_{\mathrm{edge}}(k_{y})}=\lambda_{\mathrm{% edge}}(k_{y})\ket{\psi_{\mathrm{edge}}(k_{y})}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | start_ARG italic_ψ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ = italic_λ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | start_ARG italic_ψ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ to

(R(ky)+Sβ(ky)+Sβ1(ky))|u(ky)𝑅subscript𝑘𝑦𝑆𝛽subscript𝑘𝑦superscript𝑆superscript𝛽1subscript𝑘𝑦ket𝑢subscript𝑘𝑦\displaystyle\left(R\left(k_{y}\right)+S\beta\left(k_{y}\right)+S^{{\dagger}}% \beta^{-1}\left(k_{y}\right)\right)\ket{u\left(k_{y}\right)}( italic_R ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) + italic_S italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) + italic_S start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ) | start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ =λedge(ky)|u(ky)(2xL1)absentsubscript𝜆edgesubscript𝑘𝑦ket𝑢subscript𝑘𝑦2𝑥𝐿1\displaystyle=\lambda_{\mathrm{edge}}(k_{y})\ket{u\left(k_{y}\right)}\hskip 28% .45274pt(2\leq x\leq L-1)= italic_λ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ ( 2 ≤ italic_x ≤ italic_L - 1 ) (III.21)

in the bulk and

(R(ky)+Sβ(ky))|u(ky)𝑅subscript𝑘𝑦𝑆𝛽subscript𝑘𝑦ket𝑢subscript𝑘𝑦\displaystyle\left(R\left(k_{y}\right)+S\beta\left(k_{y}\right)\right)\ket{u% \left(k_{y}\right)}( italic_R ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) + italic_S italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ) | start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ =λedge(ky)|u(ky)(x=1)absentsubscript𝜆edgesubscript𝑘𝑦ket𝑢subscript𝑘𝑦𝑥1\displaystyle=\lambda_{\mathrm{edge}}(k_{y})\ket{u\left(k_{y}\right)}\hskip 28% .45274pt(x=1)= italic_λ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ ( italic_x = 1 ) (III.22)

at the left edge. From these equations, we have

S|u(ky)=0,superscript𝑆ket𝑢subscript𝑘𝑦0\displaystyle S^{\dagger}\ket{u\left(k_{y}\right)}=0,italic_S start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT | start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ = 0 , (III.23)

further leading to

|u(ky)=(1i).ket𝑢subscript𝑘𝑦1i\displaystyle\ket{u\left(k_{y}\right)}=\left(\begin{array}[]{c}1\\ \text{i}\end{array}\right).| start_ARG italic_u ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ⟩ = ( start_ARRAY start_ROW start_CELL 1 end_CELL end_ROW start_ROW start_CELL i end_CELL end_ROW end_ARRAY ) . (III.26)

Combining Eqs. (III.22) and (III.26), we obtain

[it2β(ky)σx+(tsinky)σy+(m+cosky+t2β(ky))σz]|u=λedge(ky)|u.delimited-[]i𝑡2𝛽subscript𝑘𝑦subscript𝜎𝑥𝑡subscript𝑘𝑦subscript𝜎𝑦𝑚subscript𝑘𝑦𝑡2𝛽subscript𝑘𝑦subscript𝜎𝑧ket𝑢subscript𝜆edgesubscript𝑘𝑦ket𝑢\displaystyle\left[-\frac{\text{i}t}{2}\beta\left(k_{y}\right)\sigma_{x}+\left% (t\sin k_{y}\right)\sigma_{y}+\left(m+\cos k_{y}+\frac{t}{2}\beta\left(k_{y}% \right)\right)\sigma_{z}\right]\ket{u}=\lambda_{\mathrm{edge}}(k_{y})\ket{u}.[ - divide start_ARG i italic_t end_ARG start_ARG 2 end_ARG italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + ( italic_m + roman_cos italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + divide start_ARG italic_t end_ARG start_ARG 2 end_ARG italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ] | start_ARG italic_u end_ARG ⟩ = italic_λ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) | start_ARG italic_u end_ARG ⟩ . (III.27)

From Eq. (III.27), λedge(ky)subscript𝜆edgesubscript𝑘𝑦\lambda_{\mathrm{edge}}\left(k_{y}\right)italic_λ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) and β(ky)𝛽subscript𝑘𝑦\beta\left(k_{y}\right)italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) are given as

λedge(ky)=tsin(ky),β(ky)=(mt+cosky).formulae-sequencesubscript𝜆edgesubscript𝑘𝑦𝑡subscript𝑘𝑦𝛽subscript𝑘𝑦𝑚𝑡subscript𝑘𝑦\displaystyle\lambda_{\mathrm{edge}}(k_{y})=t\sin{k_{y}},\quad\beta(k_{y})=-% \left(\frac{m}{t}+\cos k_{y}\right).italic_λ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) , italic_β ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = - ( divide start_ARG italic_m end_ARG start_ARG italic_t end_ARG + roman_cos italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) . (III.28)

Thus, the self-energy is calculated as

Σedge(ky)=t2(1β2)2(E+iηtsinky)(σ0σy).subscriptΣedgesubscript𝑘𝑦superscript𝑡21superscript𝛽22𝐸i𝜂𝑡subscript𝑘𝑦subscript𝜎0subscript𝜎𝑦\displaystyle\Sigma_{\mathrm{edge}}(k_{y})=\frac{t^{2}(1-\beta^{2})}{2(E+\text% {i}\eta-t\sin k_{y})}(\sigma_{0}-\sigma_{y}).roman_Σ start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_β start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG 2 ( italic_E + i italic_η - italic_t roman_sin italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) end_ARG ( italic_σ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) . (III.29)

For ky=0subscript𝑘𝑦0k_{y}=0italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 0 or ky=πsubscript𝑘𝑦𝜋k_{y}=\piitalic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = italic_π, this self-energy essentially reduces to that of the SSH model [see Eq. (III.5)]. Notably, these expressions are valid only for

1>|β|=|mt+cos(ky)|1𝛽𝑚𝑡subscript𝑘𝑦\displaystyle 1>\absolutevalue{\beta}=\absolutevalue{\frac{m}{t}+\cos{k_{y}}}1 > | start_ARG italic_β end_ARG | = | start_ARG divide start_ARG italic_m end_ARG start_ARG italic_t end_ARG + roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) end_ARG | (III.30)

so that the above chiral edge state can be normalized.

IV Current formula

We derive the current formula in Eq. (8) of the main text (see also “Chapter 8: Non-equilibrium Green’s function formalism” in Ref. [103], as well as “Chapter 9: Coherent transport” and “Appendix: advanced formalism” in Ref. [113]). We consider a noninteracting fermionic Hamiltonian described by

H^=i,ja^iHija^j,^𝐻subscript𝑖𝑗superscriptsubscript^𝑎𝑖subscript𝐻𝑖𝑗subscript^𝑎𝑗\displaystyle\hat{H}=\sum_{i,j}\hat{a}_{i}^{\dagger}H_{ij}\hat{a}_{j},over^ start_ARG italic_H end_ARG = ∑ start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_H start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , (IV.1)

where a^isubscript^𝑎𝑖\hat{a}_{i}over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT (a^isubscriptsuperscript^𝑎𝑖\hat{a}^{\dagger}_{i}over^ start_ARG italic_a end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT) is a fermionic annihilation (creation) operator at site i𝑖iitalic_i, and Hijsubscript𝐻𝑖𝑗H_{ij}italic_H start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT is a single-particle Hermitian Hamiltonian. The Heisenberg equation for the annihilation operator reads

ida^i(t)dt=[a^i(t),H^]=jHija^j(t).iPlanck-constant-over-2-pi𝑑subscript^𝑎𝑖𝑡𝑑𝑡subscript^𝑎𝑖𝑡^𝐻subscript𝑗subscript𝐻𝑖𝑗subscript^𝑎𝑗𝑡\displaystyle\text{i}\hbar\frac{d\hat{a}_{i}(t)}{dt}=[\hat{a}_{i}(t),\hat{H}]=% \sum_{j}H_{ij}\hat{a}_{j}(t).i roman_ℏ divide start_ARG italic_d over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_d italic_t end_ARG = [ over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_t ) , over^ start_ARG italic_H end_ARG ] = ∑ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_t ) . (IV.2)

We rewrite this equation as

iddt(𝒂^bulk(t)𝒂^edge(t))=(HbulkTTHedge)(𝒂^bulk(t)𝒂^edge(t)),iPlanck-constant-over-2-pi𝑑𝑑𝑡matrixsubscript^𝒂bulk𝑡subscript^𝒂edge𝑡matrixsubscript𝐻bulk𝑇superscript𝑇subscript𝐻edgematrixsubscript^𝒂bulk𝑡subscript^𝒂edge𝑡\displaystyle\text{i}\hbar\frac{d}{dt}\begin{pmatrix}\hat{\boldsymbol{a}}_{% \mathrm{bulk}}(t)\\ \hat{\boldsymbol{a}}_{\mathrm{edge}}(t)\end{pmatrix}=\begin{pmatrix}H_{\rm bulk% }&T\\ T^{{\dagger}}&H_{\rm edge}\end{pmatrix}\begin{pmatrix}\hat{\boldsymbol{a}}_{% \mathrm{bulk}}(t)\\ \hat{\boldsymbol{a}}_{\mathrm{edge}}(t)\end{pmatrix},i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG ( start_ARG start_ROW start_CELL over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t ) end_CELL end_ROW start_ROW start_CELL over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t ) end_CELL end_ROW end_ARG ) = ( start_ARG start_ROW start_CELL italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT end_CELL start_CELL italic_T end_CELL end_ROW start_ROW start_CELL italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_CELL start_CELL italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) ( start_ARG start_ROW start_CELL over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t ) end_CELL end_ROW start_ROW start_CELL over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t ) end_CELL end_ROW end_ARG ) , (IV.3)

where 𝒂^bulk(t)subscript^𝒂bulk𝑡\hat{\boldsymbol{a}}_{\mathrm{bulk}}(t)over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t ) [𝒂^edge(t)subscript^𝒂edge𝑡\hat{\boldsymbol{a}}_{\mathrm{edge}}(t)over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t )] is a row vector of annihilation operators acting on the bulk (edge). We expand the annihilation operator for the bulk as

𝒂^bulk(t)=𝒂^bulk(0)(t)+𝝌^(t),subscript^𝒂bulk𝑡superscriptsubscript^𝒂bulk0𝑡^𝝌𝑡\displaystyle\hat{\boldsymbol{a}}_{\mathrm{bulk}}(t)=\hat{\boldsymbol{a}}_{% \mathrm{bulk}}^{(0)}(t)+\hat{\boldsymbol{\chi}}(t),over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t ) = over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t ) + over^ start_ARG bold_italic_χ end_ARG ( italic_t ) , (IV.4)

where 𝒂^bulk(0)(t)superscriptsubscript^𝒂bulk0𝑡\hat{\boldsymbol{a}}_{\mathrm{bulk}}^{(0)}(t)over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t ) denotes the unperturbed annihilation operator satisfying i(d/dt)𝒂^bulk(0)(t)=Hbulk𝒂^bulk(0)(t)iPlanck-constant-over-2-pi𝑑𝑑𝑡superscriptsubscript^𝒂bulk0𝑡subscript𝐻bulksuperscriptsubscript^𝒂bulk0𝑡\text{i}\hbar\left(d/dt\right)\hat{\boldsymbol{a}}_{\mathrm{bulk}}^{(0)}(t)=H_% {\rm bulk}\hat{\boldsymbol{a}}_{\mathrm{bulk}}^{(0)}(t)i roman_ℏ ( italic_d / italic_d italic_t ) over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t ) = italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t ). From the original Heisenberg equation in Eq. (IV.3), we have

(iddtHbulk)𝝌^(t)iPlanck-constant-over-2-pi𝑑𝑑𝑡subscript𝐻bulk^𝝌𝑡\displaystyle\left(\text{i}\hbar\frac{d}{dt}-H_{\rm bulk}\right)\hat{% \boldsymbol{\chi}}(t)( i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) over^ start_ARG bold_italic_χ end_ARG ( italic_t ) =T𝒂^edge(t),absent𝑇subscript^𝒂edge𝑡\displaystyle=T\hat{\boldsymbol{a}}_{\mathrm{edge}}(t),= italic_T over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t ) , (IV.5)
(iddtHedge)𝒂^edge(t)iPlanck-constant-over-2-pi𝑑𝑑𝑡subscript𝐻edgesubscript^𝒂edge𝑡\displaystyle\left(\text{i}\hbar\frac{d}{dt}-H_{\rm edge}\right)\hat{% \boldsymbol{a}}_{\mathrm{edge}}(t)( i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ) over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t ) =T(𝒂^bulk(0)(t)+𝝌^(t)).absentsuperscript𝑇superscriptsubscript^𝒂bulk0𝑡^𝝌𝑡\displaystyle=T^{\dagger}\left(\hat{\boldsymbol{a}}_{\mathrm{bulk}}^{(0)}(t)+% \hat{\boldsymbol{\chi}}(t)\right).= italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t ) + over^ start_ARG bold_italic_χ end_ARG ( italic_t ) ) . (IV.6)

We introduce the retarded Green’s function gbulk(0)(tt)subscriptsuperscript𝑔0bulk𝑡superscript𝑡g^{(0)}_{\mathrm{bulk}}(t-t^{\prime})italic_g start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) for the bulk

(iddtHbulk)gbulk(0)(tt)=δ(tt).iPlanck-constant-over-2-pi𝑑𝑑𝑡subscript𝐻bulksubscriptsuperscript𝑔0bulk𝑡superscript𝑡𝛿𝑡superscript𝑡\displaystyle\left(\text{i}\hbar\frac{d}{dt}-H_{\rm bulk}\right)g^{(0)}_{% \mathrm{bulk}}(t-t^{\prime})=\delta(t-t^{\prime}).( i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ) italic_g start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = italic_δ ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) . (IV.7)

Here, the argument of the Green’s function depends only on the time difference tt𝑡superscript𝑡t-t^{\prime}italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT since the Hamiltonian Hbulksubscript𝐻bulkH_{\mathrm{bulk}}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT is independent of time. Using the Green’s function, we obtain the formal solution of Eq. (IV.5) by the convolution integral,

𝝌^(t)=𝑑tgbulk(0)(tt)T𝒂^edge(t),^𝝌𝑡differential-dsuperscript𝑡subscriptsuperscript𝑔0bulk𝑡superscript𝑡𝑇subscript^𝒂edgesuperscript𝑡\displaystyle\hat{\boldsymbol{\chi}}(t)=\int dt^{\prime}g^{(0)}_{\mathrm{bulk}% }(t-t^{\prime})T\hat{\boldsymbol{a}}_{\mathrm{edge}}(t^{\prime}),over^ start_ARG bold_italic_χ end_ARG ( italic_t ) = ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_T over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (IV.8)

and hence

(iddtHedge)𝒂^edge(t)=T𝒂^bulk(0)(t)+𝑑tΣ(tt)𝒂^edge(t)iPlanck-constant-over-2-pi𝑑𝑑𝑡subscript𝐻edgesubscript^𝒂edge𝑡superscript𝑇superscriptsubscript^𝒂bulk0𝑡differential-dsuperscript𝑡Σ𝑡superscript𝑡subscript^𝒂edgesuperscript𝑡\displaystyle\left(\text{i}\hbar\frac{d}{dt}-H_{\rm edge}\right)\hat{% \boldsymbol{a}}_{\mathrm{edge}}(t)=T^{\dagger}\hat{\boldsymbol{a}}_{\mathrm{% bulk}}^{(0)}(t)+\int dt^{\prime}\Sigma(t-t^{\prime})\hat{\boldsymbol{a}}_{% \mathrm{edge}}(t^{\prime})( i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ) over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t ) = italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t ) + ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT roman_Σ ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) (IV.9)

with the self-energy defined as

Σ(tt)Tgbulk(0)(tt)T.Σ𝑡superscript𝑡superscript𝑇subscriptsuperscript𝑔0bulk𝑡superscript𝑡𝑇\Sigma(t-t^{\prime})\coloneqq T^{\dagger}g^{(0)}_{\mathrm{bulk}}(t-t^{\prime})T.roman_Σ ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ≔ italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_g start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_T . (IV.10)

Similarly, introducing the retarded Green’s function gedge(tt)subscript𝑔edge𝑡superscript𝑡g_{\mathrm{edge}}(t-t^{\prime})italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) for the edge by

(iddtHedge)gedge(tt)=𝑑t′′Σ(tt′′)gedge(t′′t)+δ(tt),iPlanck-constant-over-2-pi𝑑𝑑𝑡subscript𝐻edgesubscript𝑔edge𝑡superscript𝑡differential-dsuperscript𝑡′′Σ𝑡superscript𝑡′′subscript𝑔edgesuperscript𝑡′′superscript𝑡𝛿𝑡superscript𝑡\displaystyle\left(\text{i}\hbar\frac{d}{dt}-H_{\rm edge}\right)g_{\mathrm{% edge}}(t-t^{\prime})=\int dt^{\prime\prime}\Sigma(t-t^{\prime\prime})g_{% \mathrm{edge}}(t^{\prime\prime}-t^{\prime})+\delta(t-t^{\prime}),( i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ) italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT roman_Σ ( italic_t - italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) + italic_δ ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (IV.11)

we obtain the formal solution of Eq. (IV.9) as

𝒂^edge(t)=𝑑tgedge(tt)T𝒂^bulk(0)(t).subscript^𝒂edge𝑡differential-dsuperscript𝑡subscript𝑔edge𝑡superscript𝑡superscript𝑇superscriptsubscript^𝒂bulk0superscript𝑡\displaystyle\hat{\boldsymbol{a}}_{\mathrm{edge}}(t)=\int dt^{\prime}g_{% \mathrm{edge}}(t-t^{\prime})T^{\dagger}\hat{\boldsymbol{a}}_{\mathrm{bulk}}^{(% 0)}(t^{\prime}).over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t ) = ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) . (IV.12)

Similar calculations for the creation operator lead to

(iddt+HedgeT)𝒂^edge(t)=𝑑tΣ(tt)𝒂^edge(t)TT𝒂^bulk(0)(t)iPlanck-constant-over-2-pi𝑑𝑑𝑡superscriptsubscript𝐻edge𝑇superscriptsubscript^𝒂edge𝑡differential-dsuperscript𝑡superscriptΣ𝑡superscript𝑡superscriptsubscript^𝒂edgesuperscript𝑡superscript𝑇𝑇superscriptsubscript^𝒂bulksuperscript0𝑡\displaystyle\left(\text{i}\hbar\frac{d}{dt}+H_{\rm edge}^{T}\right)\hat{% \boldsymbol{a}}_{\mathrm{edge}}^{\dagger}(t)=-\int dt^{\prime}\Sigma^{*}(t-t^{% \prime})\hat{\boldsymbol{a}}_{\mathrm{edge}}^{\dagger}(t^{\prime})-T^{T}\hat{% \boldsymbol{a}}_{\mathrm{bulk}}^{(0)^{{\dagger}}}(t)( i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG + italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t ) = - ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT roman_Σ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) - italic_T start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ( italic_t ) (IV.13)

and hence

𝒂^edge(t)=𝑑tgedge(tt)TT𝒂^bulk(0)(t).superscriptsubscript^𝒂edge𝑡differential-dsuperscript𝑡superscriptsubscript𝑔edge𝑡superscript𝑡superscript𝑇𝑇superscriptsubscript^𝒂bulksuperscript0superscript𝑡\displaystyle\hat{\boldsymbol{a}}_{\mathrm{edge}}^{\dagger}(t)=\int dt^{\prime% }g_{\mathrm{edge}}^{*}(t-t^{\prime})T^{T}\hat{\boldsymbol{a}}_{\mathrm{bulk}}^% {(0)^{{\dagger}}}(t^{\prime}).over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t ) = ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_T start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT over^ start_ARG bold_italic_a end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) . (IV.14)

We define the Fourier transform of f(tt)𝑓𝑡superscript𝑡f(t-t^{\prime})italic_f ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) as

[f(tt)]d(tt)eiE(tt)f(tt).delimited-[]𝑓𝑡superscript𝑡𝑑𝑡superscript𝑡superscript𝑒i𝐸Planck-constant-over-2-pi𝑡superscript𝑡𝑓𝑡superscript𝑡\displaystyle\mathcal{F}\left[f(t-t^{\prime})\right]\coloneqq\int d(t-t^{% \prime})e^{\text{i}\frac{E}{\hbar}(t-t^{\prime})}f(t-t^{\prime}).caligraphic_F [ italic_f ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ] ≔ ∫ italic_d ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT i divide start_ARG italic_E end_ARG start_ARG roman_ℏ end_ARG ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_f ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) . (IV.15)

The Fourier transforms of the Green’s functions gbulk(0)(tt)subscriptsuperscript𝑔0bulk𝑡superscript𝑡g^{(0)}_{\mathrm{bulk}}(t-t^{\prime})italic_g start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) in Eq. (IV.7) and gedge(tt)subscript𝑔edge𝑡superscript𝑡g_{\mathrm{edge}}(t-t^{\prime})italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) in Eq. (IV.11) are

Gbulk(0)(E)[gbulk(0)(tt)]=1E+iηHbulk,Gedge(E)[gedge(tt)]=1E+iηHedgeΣ(E),formulae-sequencesuperscriptsubscript𝐺bulk0𝐸delimited-[]superscriptsubscript𝑔bulk0𝑡superscript𝑡1𝐸i𝜂subscript𝐻bulksubscript𝐺edge𝐸delimited-[]subscript𝑔edge𝑡superscript𝑡1𝐸i𝜂subscript𝐻edgeΣ𝐸\displaystyle G_{\mathrm{bulk}}^{(0)}(E)\coloneqq\mathcal{F}\left[g_{\mathrm{% bulk}}^{(0)}(t-t^{\prime})\right]=\frac{1}{E+\text{i}\eta-H_{\mathrm{bulk}}},% \quad G_{\mathrm{edge}}(E)\coloneqq\mathcal{F}\left[g_{\mathrm{edge}}(t-t^{% \prime})\right]=\frac{1}{E+\text{i}\eta-H_{\mathrm{edge}}-\Sigma(E)},italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) ≔ caligraphic_F [ italic_g start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ] = divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT end_ARG , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ≔ caligraphic_F [ italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ] = divide start_ARG 1 end_ARG start_ARG italic_E + i italic_η - italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT - roman_Σ ( italic_E ) end_ARG , (IV.16)

with Σ(E)[Σ(tt)]Σ𝐸delimited-[]Σ𝑡superscript𝑡\Sigma(E)\coloneqq\mathcal{F}\left[\Sigma(t-t^{\prime})\right]roman_Σ ( italic_E ) ≔ caligraphic_F [ roman_Σ ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ]. Here, an infinitesimal number η>0𝜂0\eta>0italic_η > 0 is introduced so that the initial condition can be satisfied appropriately.

Now, we introduce the two-time current operator [103, 113]

J(t,t)(ddt+ddt)C(t,t),𝐽𝑡superscript𝑡𝑑𝑑𝑡𝑑𝑑superscript𝑡𝐶𝑡superscript𝑡\displaystyle J(t,t^{\prime})\coloneqq\left(\frac{d}{dt}+\frac{d}{dt^{\prime}}% \right)C(t,t^{\prime}),italic_J ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ≔ ( divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG + divide start_ARG italic_d end_ARG start_ARG italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG ) italic_C ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (IV.17)

with the two-time correlation function

Cij(t,t)a^j(t)a^i(t),subscript𝐶𝑖𝑗𝑡superscript𝑡delimited-⟨⟩superscriptsubscript^𝑎𝑗superscript𝑡subscript^𝑎𝑖𝑡\displaystyle C_{ij}(t,t^{\prime})\coloneqq\langle\hat{a}_{j}^{\dagger}(t^{% \prime})\hat{a}_{i}(t)\rangle,italic_C start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ≔ ⟨ over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_t ) ⟩ , (IV.18)

where the bracket denotes the average over a thermal equilibrium state. In the limit ttsuperscript𝑡𝑡t^{\prime}\rightarrow titalic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT → italic_t, the diagonal element of J(t,t)𝐽𝑡superscript𝑡J(t,t^{\prime})italic_J ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) gives the local current. From Eqs. (IV.12) and (IV.13), the correlation function of the edge, Cedge;ij(t,t)a^edge;j(t)a^edge;i(t)subscript𝐶edge𝑖𝑗𝑡superscript𝑡delimited-⟨⟩superscriptsubscript^𝑎edge𝑗superscript𝑡subscript^𝑎edge𝑖𝑡C_{\mathrm{edge};ij}(t,t^{\prime})\coloneqq\langle\hat{a}_{{\mathrm{edge}};j}^% {\dagger}(t^{\prime})\hat{a}_{{\mathrm{edge}};i}(t)\rangleitalic_C start_POSTSUBSCRIPT roman_edge ; italic_i italic_j end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ≔ ⟨ over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT roman_edge ; italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) over^ start_ARG italic_a end_ARG start_POSTSUBSCRIPT roman_edge ; italic_i end_POSTSUBSCRIPT ( italic_t ) ⟩, satisfies

iddtCedge(t,t)=Cedge(t,t)Hedge𝑑t′′gedge(t,t′′)TCbulk(0)(t,t′′)T𝑑t′′Cedge(t′′,t)Σ(t,t′′),iPlanck-constant-over-2-pi𝑑𝑑𝑡subscript𝐶edge𝑡superscript𝑡subscript𝐶edge𝑡superscript𝑡subscript𝐻edgedifferential-dsuperscript𝑡′′subscript𝑔edgesuperscript𝑡superscript𝑡′′superscript𝑇superscriptsubscript𝐶bulk0𝑡superscript𝑡′′𝑇differential-dsuperscript𝑡′′subscript𝐶edgesuperscript𝑡′′superscript𝑡superscriptΣ𝑡superscript𝑡′′\displaystyle\text{i}\hbar\frac{d}{dt}C_{\mathrm{edge}}(t,t^{\prime})=-C_{% \mathrm{edge}}(t,t^{\prime})H_{\mathrm{edge}}-\int dt^{\prime\prime}g_{\mathrm% {edge}}(t^{\prime},t^{\prime\prime})T^{\dagger}C_{\mathrm{bulk}}^{(0)}(t,t^{% \prime\prime})T-\int dt^{\prime\prime}C_{\mathrm{edge}}(t^{\prime\prime},t^{% \prime})\Sigma^{\dagger}(t,t^{\prime\prime}),i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t end_ARG italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = - italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT - ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_T - ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) , (IV.19)

with the isolated bulk correlation function Cbulk;ij(0)(t,t)a^bulk;j(0)(t)a^bulk;i(0)(t)superscriptsubscript𝐶bulk𝑖𝑗0𝑡superscript𝑡delimited-⟨⟩subscriptsuperscript^𝑎superscript0bulk𝑗superscript𝑡subscriptsuperscript^𝑎0bulk𝑖𝑡C_{\mathrm{bulk};ij}^{(0)}(t,t^{\prime})\coloneqq\langle\hat{a}^{(0)^{{\dagger% }}}_{\mathrm{bulk};j}(t^{\prime})\hat{a}^{(0)}_{\mathrm{bulk};i}(t)\rangleitalic_C start_POSTSUBSCRIPT roman_bulk ; italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ≔ ⟨ over^ start_ARG italic_a end_ARG start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_j end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) over^ start_ARG italic_a end_ARG start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_i end_POSTSUBSCRIPT ( italic_t ) ⟩. Here, we assume that the bulk maintains the local equilibrium even after attaching the edge [103, 113]. Similarly, from Eqs. (IV.9) and (IV.14), we have

iddtCedge(t,t)=HedgeCedge(t,t)+𝑑t′′TCbulk(0)(t′′,t)Tgedge(t,t′′)+𝑑t′′Σ(t,t′′)Cedge(t,t′′).iPlanck-constant-over-2-pi𝑑𝑑superscript𝑡subscript𝐶edge𝑡superscript𝑡subscript𝐻edgesubscript𝐶edge𝑡superscript𝑡differential-dsuperscript𝑡′′superscript𝑇superscriptsubscript𝐶bulk0superscript𝑡′′superscript𝑡𝑇superscriptsubscript𝑔edge𝑡superscript𝑡′′differential-dsuperscript𝑡′′Σsuperscript𝑡superscript𝑡′′subscript𝐶edge𝑡superscript𝑡′′\displaystyle\text{i}\hbar\frac{d}{dt^{\prime}}C_{\mathrm{edge}}(t,t^{\prime})% =H_{\mathrm{edge}}C_{\mathrm{edge}}(t,t^{\prime})+\int dt^{\prime\prime}T^{% \dagger}C_{\mathrm{bulk}}^{(0)}(t^{\prime\prime},t^{\prime})Tg_{\mathrm{edge}}% ^{\dagger}(t,t^{\prime\prime})+\int dt^{\prime\prime}\Sigma(t^{\prime},t^{% \prime\prime})C_{\mathrm{edge}}(t,t^{\prime\prime}).i roman_ℏ divide start_ARG italic_d end_ARG start_ARG italic_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) + ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_T italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) + ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT roman_Σ ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) . (IV.20)

The first term of the right-hand side in Eqs. (IV.19) and (IV.20) represents the particle current inside the edge; we below omit it since it is irrelevant to the current between the bulk and edge, and hence non-Hermitian topology. Additionally, the second (third) term in Eqs. (IV.19) and (IV.20) describes the particle current from the bulk to the edge (from the edge to the bulk). Then, the net current is

J(t,t)=Jinflow(t,t)Joutflow(t,t)𝐽𝑡superscript𝑡subscript𝐽inflow𝑡superscript𝑡subscript𝐽outflow𝑡superscript𝑡\displaystyle J(t,t^{\prime})=J_{\mathrm{inflow}}(t,t^{\prime})-J_{\mathrm{% outflow}}(t,t^{\prime})italic_J ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) - italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) (IV.21)

with

iJinflow(t,t)iPlanck-constant-over-2-pisubscript𝐽inflow𝑡superscript𝑡\displaystyle\text{i}\hbar J_{\mathrm{inflow}}(t,t^{\prime})i roman_ℏ italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) 𝑑t′′gedge(t,t′′)TCbulk(0)(t,t′′)T+𝑑t′′TCbulk(0)(t′′,t)Tgedge(t,t′′),absentdifferential-dsuperscript𝑡′′subscript𝑔edgesuperscript𝑡superscript𝑡′′superscript𝑇superscriptsubscript𝐶bulk0𝑡superscript𝑡′′𝑇differential-dsuperscript𝑡′′superscript𝑇superscriptsubscript𝐶bulk0superscript𝑡′′superscript𝑡𝑇superscriptsubscript𝑔edge𝑡superscript𝑡′′\displaystyle\coloneqq-\int dt^{\prime\prime}g_{\mathrm{edge}}(t^{\prime},t^{% \prime\prime})T^{\dagger}C_{\mathrm{bulk}}^{(0)}(t,t^{\prime\prime})T+\int dt^% {\prime\prime}T^{\dagger}C_{\mathrm{bulk}}^{(0)}(t^{\prime\prime},t^{\prime})% Tg_{\mathrm{edge}}^{\dagger}(t,t^{\prime\prime}),≔ - ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_T + ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_T italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) , (IV.22)
iJoutflow(t,t)iPlanck-constant-over-2-pisubscript𝐽outflow𝑡superscript𝑡\displaystyle\text{i}\hbar J_{\mathrm{outflow}}(t,t^{\prime})i roman_ℏ italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) 𝑑t′′Cedge(t′′,t)Σ(t,t′′)𝑑t′′Σ(t,t′′)Cedge(t,t′′).absentdifferential-dsuperscript𝑡′′subscript𝐶edgesuperscript𝑡′′superscript𝑡superscriptΣ𝑡superscript𝑡′′differential-dsuperscript𝑡′′Σsuperscript𝑡superscript𝑡′′subscript𝐶edge𝑡superscript𝑡′′\displaystyle\coloneqq\int dt^{\prime\prime}C_{\mathrm{edge}}(t^{\prime\prime}% ,t^{\prime})\Sigma^{\dagger}(t,t^{\prime\prime})-\int dt^{\prime\prime}\Sigma(% t^{\prime},t^{\prime\prime})C_{\mathrm{edge}}(t,t^{\prime\prime}).≔ ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) - ∫ italic_d italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT roman_Σ ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) . (IV.23)

At thermal equilibrium, from the cyclicity of the trace, the correlation functions Cbulk(0)(t,t)superscriptsubscript𝐶bulk0𝑡superscript𝑡C_{\mathrm{bulk}}^{(0)}(t,t^{\prime})italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) and Cedge(t,t)subscript𝐶edge𝑡superscript𝑡C_{\mathrm{edge}}(t,t^{\prime})italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) depend only on the argument tt𝑡superscript𝑡t-t^{\prime}italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. After the Fourier transformation, the E𝐸Eitalic_E-resolved currents read,

iJinflow(E)iPlanck-constant-over-2-pisubscript𝐽inflow𝐸\displaystyle\text{i}\hbar J_{\mathrm{inflow}}(E)i roman_ℏ italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT ( italic_E ) [Jinflow(tt)]=Gedge(E)TCbulk(0)(E)T+TCbulk(0)(E)TGedge(E),absentdelimited-[]subscript𝐽inflow𝑡superscript𝑡subscript𝐺edge𝐸superscript𝑇superscriptsubscript𝐶bulk0𝐸𝑇superscript𝑇superscriptsubscript𝐶bulk0𝐸𝑇superscriptsubscript𝐺edge𝐸\displaystyle\coloneqq\mathcal{F}\left[J_{\mathrm{inflow}}(t-t^{\prime})\right% ]=-G_{\mathrm{edge}}(E)T^{\dagger}C_{\mathrm{bulk}}^{(0)}(E)T+T^{\dagger}C_{% \mathrm{bulk}}^{(0)}(E)TG_{\mathrm{edge}}^{\dagger}(E),≔ caligraphic_F [ italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ] = - italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) italic_T + italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) italic_T italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) , (IV.24)
iJoutflow(E)iPlanck-constant-over-2-pisubscript𝐽outflow𝐸\displaystyle\text{i}\hbar J_{\mathrm{outflow}}(E)i roman_ℏ italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT ( italic_E ) [Joutflow(tt)]=Cedge(E)Σ(E)Σ(E)Cedge(E).absentdelimited-[]subscript𝐽outflow𝑡superscript𝑡subscript𝐶edge𝐸superscriptΣ𝐸Σ𝐸subscript𝐶edge𝐸\displaystyle\coloneqq\mathcal{F}\left[J_{\mathrm{outflow}}(t-t^{\prime})% \right]=C_{\mathrm{edge}}(E)\Sigma^{\dagger}(E)-\Sigma(E)C_{\mathrm{edge}}(E).≔ caligraphic_F [ italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT ( italic_t - italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ] = italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) - roman_Σ ( italic_E ) italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) . (IV.25)

To calculate the correlation function Cbulk(0)(E)superscriptsubscript𝐶bulk0𝐸C_{\mathrm{bulk}}^{(0)}(E)italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) of the isolated bulk, we diagonalize the Hamiltonian,

H^bulk=nEnγ^bulk;n(0)γ^bulk;n(0),γ^bulk;n(0)ma^bulk;m(0)Umn,formulae-sequencesubscript^𝐻bulksubscript𝑛subscript𝐸𝑛subscriptsuperscript^𝛾superscript0bulk𝑛subscriptsuperscript^𝛾0bulk𝑛subscriptsuperscript^𝛾superscript0bulk𝑛subscript𝑚subscriptsuperscript^𝑎superscript0bulk𝑚subscript𝑈𝑚𝑛\hat{H}_{\mathrm{bulk}}=\sum_{n}E_{n}\hat{\gamma}^{{(0)}^{\dagger}}_{\mathrm{% bulk};n}\hat{\gamma}^{(0)}_{\mathrm{bulk};n},\quad\hat{\gamma}^{(0)^{\dagger}}% _{\mathrm{bulk};n}\coloneqq\sum_{m}\hat{a}^{(0)^{\dagger}}_{\mathrm{bulk};m}U_% {mn},over^ start_ARG italic_H end_ARG start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_n end_POSTSUBSCRIPT over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_n end_POSTSUBSCRIPT , over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_n end_POSTSUBSCRIPT ≔ ∑ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT over^ start_ARG italic_a end_ARG start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_m end_POSTSUBSCRIPT italic_U start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT , (IV.26)

where Ensubscript𝐸𝑛E_{n}italic_E start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is eigenenergy of the single-particle Hamiltonian Hbulksubscript𝐻bulkH_{\rm bulk}italic_H start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT, and Umnsubscript𝑈𝑚𝑛U_{mn}italic_U start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT is a unitary matrix of the single-particle eigenstates. Then, we have

Cbulk;ij(0)(E)=kUikUkjδ(EEk)γ^bulk;k(0)γ^bulk;k(0)=f0(E)kUikUkjδ(EEk),superscriptsubscript𝐶bulk𝑖𝑗0𝐸subscript𝑘subscript𝑈𝑖𝑘subscriptsuperscript𝑈𝑘𝑗𝛿𝐸subscript𝐸𝑘delimited-⟨⟩subscriptsuperscript^𝛾superscript0bulk𝑘subscriptsuperscript^𝛾0bulk𝑘subscript𝑓0𝐸subscript𝑘subscript𝑈𝑖𝑘superscriptsubscript𝑈𝑘𝑗𝛿𝐸subscript𝐸𝑘\displaystyle C_{\mathrm{bulk};ij}^{(0)}(E)=\sum\limits_{k}U_{ik}U^{\dagger}_{% kj}\delta(E-E_{k})\langle\hat{\gamma}^{(0)^{\dagger}}_{\mathrm{bulk};k}\hat{% \gamma}^{(0)}_{\mathrm{bulk};k}\rangle=f_{0}\left(E\right)\sum_{k}U_{ik}U_{kj}% ^{{\dagger}}\delta\left(E-E_{k}\right),italic_C start_POSTSUBSCRIPT roman_bulk ; italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) = ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_U start_POSTSUBSCRIPT italic_i italic_k end_POSTSUBSCRIPT italic_U start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k italic_j end_POSTSUBSCRIPT italic_δ ( italic_E - italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) ⟨ over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_k end_POSTSUBSCRIPT over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_k end_POSTSUBSCRIPT ⟩ = italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_U start_POSTSUBSCRIPT italic_i italic_k end_POSTSUBSCRIPT italic_U start_POSTSUBSCRIPT italic_k italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_δ ( italic_E - italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) , (IV.27)

with the Fermi distribution function f0(Ek)=γ^bulk;k(0)γ^bulk;k(0)subscript𝑓0subscript𝐸𝑘delimited-⟨⟩subscriptsuperscript^𝛾superscript0bulk𝑘subscriptsuperscript^𝛾0bulk𝑘f_{0}\left(E_{k}\right)=\langle\hat{\gamma}^{(0)^{\dagger}}_{\mathrm{bulk};k}% \hat{\gamma}^{(0)}_{\mathrm{bulk};k}\rangleitalic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) = ⟨ over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_k end_POSTSUBSCRIPT over^ start_ARG italic_γ end_ARG start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_bulk ; italic_k end_POSTSUBSCRIPT ⟩. Using the Green’s function in Eq. (IV.16), i.e.,

Gbulk;ij(0)(E)=kUikUkjEEk+iη=kUikUkj[𝒫(1EEk)iπδ(EEk)],superscriptsubscript𝐺bulk𝑖𝑗0𝐸subscript𝑘subscript𝑈𝑖𝑘subscriptsuperscript𝑈𝑘𝑗𝐸subscript𝐸𝑘i𝜂subscript𝑘subscript𝑈𝑖𝑘subscriptsuperscript𝑈𝑘𝑗delimited-[]𝒫1𝐸subscript𝐸𝑘i𝜋𝛿𝐸subscript𝐸𝑘G_{\mathrm{bulk};ij}^{(0)}(E)=\sum_{k}\frac{U_{ik}U^{{\dagger}}_{kj}}{E-E_{k}+% \text{i}\eta}=\sum_{k}U_{ik}U^{{\dagger}}_{kj}\left[\mathcal{P}\left(\frac{1}{% E-E_{k}}\right)-\text{i}\pi\delta\left(E-E_{k}\right)\right],italic_G start_POSTSUBSCRIPT roman_bulk ; italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) = ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT divide start_ARG italic_U start_POSTSUBSCRIPT italic_i italic_k end_POSTSUBSCRIPT italic_U start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k italic_j end_POSTSUBSCRIPT end_ARG start_ARG italic_E - italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT + i italic_η end_ARG = ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_U start_POSTSUBSCRIPT italic_i italic_k end_POSTSUBSCRIPT italic_U start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k italic_j end_POSTSUBSCRIPT [ caligraphic_P ( divide start_ARG 1 end_ARG start_ARG italic_E - italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG ) - i italic_π italic_δ ( italic_E - italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) ] , (IV.28)

we have

Cbulk(0)(E)=if0(E)2π[Gbulk(0)(E)Gbulk(0)(E)],superscriptsubscript𝐶bulk0𝐸isubscript𝑓0𝐸2𝜋delimited-[]superscriptsubscript𝐺bulk0𝐸superscriptsubscript𝐺bulksuperscript0𝐸\displaystyle C_{\mathrm{bulk}}^{(0)}(E)=\frac{\text{i}f_{0}(E)}{2\pi}\left[G_% {\mathrm{bulk}}^{(0)}(E)-G_{\mathrm{bulk}}^{(0)^{\dagger}}(E)\right],italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) = divide start_ARG i italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π end_ARG [ italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) - italic_G start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ( italic_E ) ] , (IV.29)

and hence

TCbulk(0)(E)T=if0(E)2π[Σ(E)Σ(E)]=f0(E)2πΓ(E),superscript𝑇superscriptsubscript𝐶bulk0𝐸𝑇isubscript𝑓0𝐸2𝜋delimited-[]Σ𝐸superscriptΣ𝐸subscript𝑓0𝐸2𝜋Γ𝐸\displaystyle T^{\dagger}C_{\mathrm{bulk}}^{(0)}(E)T=\frac{\text{i}f_{0}(E)}{2% \pi}\left[\Sigma(E)-\Sigma^{\dagger}(E)\right]=\frac{f_{0}(E)}{2\pi}\Gamma(E),italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) italic_T = divide start_ARG i italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π end_ARG [ roman_Σ ( italic_E ) - roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) ] = divide start_ARG italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π end_ARG roman_Γ ( italic_E ) , (IV.30)

with the broadening matrix

Γ(E)i[Σ(E)Σ(E)].Γ𝐸idelimited-[]Σ𝐸superscriptΣ𝐸\displaystyle\Gamma(E)\coloneqq\text{i}\left[\Sigma(E)-\Sigma^{\dagger}(E)% \right].roman_Γ ( italic_E ) ≔ i [ roman_Σ ( italic_E ) - roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) ] . (IV.31)

Next, from Eqs. (IV.12) and (IV.14), we have

Cedge(t,t)=𝑑t1𝑑t2gedge(t,t2)TCbulk(0)(t2,t1)Tgedge(t,t1)subscript𝐶edge𝑡superscript𝑡differential-dsubscript𝑡1differential-dsubscript𝑡2subscript𝑔edgesuperscript𝑡subscript𝑡2superscript𝑇superscriptsubscript𝐶bulk0subscript𝑡2subscript𝑡1𝑇superscriptsubscript𝑔edge𝑡subscript𝑡1\displaystyle C_{\mathrm{edge}}(t,t^{\prime})=\int dt_{1}\int dt_{2}g_{\mathrm% {edge}}(t^{\prime},t_{2})T^{\dagger}C_{\mathrm{bulk}}^{(0)}(t_{2},t_{1})Tg_{% \mathrm{edge}}^{\dagger}(t,t_{1})italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t , italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = ∫ italic_d italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∫ italic_d italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_T italic_g start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_t , italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) (IV.32)

and hence

Cedge(E)=Gedge(E)TCbulk(0)(E)TGedge(E),subscript𝐶edge𝐸subscript𝐺edge𝐸superscript𝑇superscriptsubscript𝐶bulk0𝐸𝑇superscriptsubscript𝐺edge𝐸\displaystyle C_{\mathrm{edge}}(E)=G_{\mathrm{edge}}(E)T^{\dagger}C_{\mathrm{% bulk}}^{(0)}(E)TG_{\mathrm{edge}}^{\dagger}(E),italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) = italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) italic_T start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_bulk end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_E ) italic_T italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) , (IV.33)

further leading to

Cedge(E)=f0(E)2πGedge(E)Γ(E)Gedge(E)=f0(E)2πAedge(E),subscript𝐶edge𝐸subscript𝑓0𝐸2𝜋subscript𝐺edge𝐸Γ𝐸superscriptsubscript𝐺edge𝐸subscript𝑓0𝐸2𝜋subscript𝐴edge𝐸\displaystyle C_{\mathrm{edge}}(E)=\frac{f_{0}(E)}{2\pi}G_{\mathrm{edge}}(E)% \Gamma(E)G_{\mathrm{edge}}^{\dagger}(E)=\frac{f_{0}(E)}{2\pi}A_{\mathrm{edge}}% (E),italic_C start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) = divide start_ARG italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π end_ARG italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) roman_Γ ( italic_E ) italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) = divide start_ARG italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π end_ARG italic_A start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) , (IV.34)

where Aedge(E)subscript𝐴edge𝐸A_{\mathrm{edge}}(E)italic_A start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) is the spectral function

Aedge(E)i[Gedge(E)Gedge(E)]=Gedge(E)Γ(E)Gedge(E).subscript𝐴edge𝐸idelimited-[]subscript𝐺edge𝐸superscriptsubscript𝐺edge𝐸subscript𝐺edge𝐸Γ𝐸superscriptsubscript𝐺edge𝐸\displaystyle A_{\mathrm{edge}}(E)\coloneqq\text{i}\left[G_{\mathrm{edge}}(E)-% G_{\mathrm{edge}}^{\dagger}(E)\right]=G_{\mathrm{edge}}(E)\Gamma(E)G_{\mathrm{% edge}}^{\dagger}(E).italic_A start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ≔ i [ italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) - italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) ] = italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) roman_Γ ( italic_E ) italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) . (IV.35)

Then, from Eqs. (IV.30) and (IV.34), we obtain

Jinflow(E)subscript𝐽inflow𝐸\displaystyle J_{\mathrm{inflow}}(E)italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT ( italic_E ) =f0(E)2πi[Γ(E)Gedge(E)Gedge(E)Γ(E)],absentsubscript𝑓0𝐸2𝜋iPlanck-constant-over-2-pidelimited-[]Γ𝐸superscriptsubscript𝐺edge𝐸subscript𝐺edge𝐸Γ𝐸\displaystyle=\frac{f_{0}(E)}{2\pi\text{i}\hbar}\left[\Gamma(E)G_{\mathrm{edge% }}^{\dagger}(E)-G_{\mathrm{edge}}(E)\Gamma(E)\right],= divide start_ARG italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π i roman_ℏ end_ARG [ roman_Γ ( italic_E ) italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) - italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) roman_Γ ( italic_E ) ] , (IV.36)
Joutflow(E)subscript𝐽outflow𝐸\displaystyle J_{\mathrm{outflow}}(E)italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT ( italic_E ) =f0(E)2πi[Aedge(E)Σ(E)Σ(E)Aedge(E)],absentsubscript𝑓0𝐸2𝜋iPlanck-constant-over-2-pidelimited-[]subscript𝐴edge𝐸superscriptΣ𝐸Σ𝐸subscript𝐴edge𝐸\displaystyle=\frac{f_{0}(E)}{2\pi\text{i}\hbar}\left[A_{\mathrm{edge}}(E)% \Sigma^{\dagger}(E)-\Sigma(E)A_{\mathrm{edge}}(E)\right],= divide start_ARG italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π i roman_ℏ end_ARG [ italic_A start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) roman_Σ start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_E ) - roman_Σ ( italic_E ) italic_A start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ] , (IV.37)

and hence

J(E)=Jinflow(E)Joutflow(E)=f0(E)2π([Σ(E),Gedge(E)]+[Σ(E),Gedge(E)])𝐽𝐸subscript𝐽inflow𝐸subscript𝐽outflow𝐸subscript𝑓0𝐸2𝜋Planck-constant-over-2-piΣ𝐸subscript𝐺edge𝐸superscriptΣ𝐸subscript𝐺edge𝐸\displaystyle J(E)=J_{\mathrm{inflow}}(E)-J_{\mathrm{outflow}}(E)=\frac{f_{0}(% E)}{2\pi\hbar}\left(\left[\Sigma(E),G_{\mathrm{edge}}(E)\right]+\left[\Sigma(E% ),G_{\mathrm{edge}}(E)\right]^{\dagger}\right)italic_J ( italic_E ) = italic_J start_POSTSUBSCRIPT roman_inflow end_POSTSUBSCRIPT ( italic_E ) - italic_J start_POSTSUBSCRIPT roman_outflow end_POSTSUBSCRIPT ( italic_E ) = divide start_ARG italic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_E ) end_ARG start_ARG 2 italic_π roman_ℏ end_ARG ( [ roman_Σ ( italic_E ) , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ] + [ roman_Σ ( italic_E ) , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_E ) ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) (IV.38)

with the effective non-Hermitian Hamiltonian Heff(E)Hedge+Σ(E)subscript𝐻eff𝐸subscript𝐻edgeΣ𝐸H_{\mathrm{eff}}(E)\coloneqq H_{\mathrm{edge}}+\Sigma(E)italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_E ) ≔ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT + roman_Σ ( italic_E ). From the relation

[Σ,Gedge]+[Σ,Gedge]=[Hedge,Gedge][Hedge,Gedge],Σsubscript𝐺edgesuperscriptΣsubscript𝐺edgesubscript𝐻edgesubscript𝐺edgesuperscriptsubscript𝐻edgesubscript𝐺edge\displaystyle\big{[}\Sigma,G_{\mathrm{edge}}\big{]}+\big{[}\Sigma,G_{\mathrm{% edge}}\big{]}^{\dagger}=-\big{[}H_{\mathrm{edge}},G_{\mathrm{edge}}\big{]}-% \big{[}H_{\mathrm{edge}},G_{\mathrm{edge}}\big{]}^{\dagger},[ roman_Σ , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] + [ roman_Σ , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT = - [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] - [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT , (IV.39)

Eq. (IV.38) reduces to Eq. (8) in the main text. In passing, “f0/2πsubscript𝑓02𝜋Planck-constant-over-2-pif_{0}/2\pi\hbaritalic_f start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / 2 italic_π roman_ℏ” is dropped in Eq. (8).

We also introduce the E𝐸Eitalic_E-resolved current density j(y,E)𝑗𝑦𝐸j(y,E)italic_j ( italic_y , italic_E ) at site y𝑦yitalic_y, summing up the internal degree of freedom α𝛼\alphaitalic_α as

j(y,E)αy,α|([Σ,Gedge]+[Σ,Gedge])|y,α.𝑗𝑦𝐸subscript𝛼bra𝑦𝛼Σsubscript𝐺edgesuperscriptΣsubscript𝐺edgeket𝑦𝛼\displaystyle j(y,E)\coloneqq\sum\limits_{\alpha}\bra{y,\alpha}\left(\big{[}% \Sigma,G_{\mathrm{edge}}\big{]}+\big{[}\Sigma,G_{\mathrm{edge}}\big{]}^{% \dagger}\right)\ket{y,\alpha}.italic_j ( italic_y , italic_E ) ≔ ∑ start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ⟨ start_ARG italic_y , italic_α end_ARG | ( [ roman_Σ , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] + [ roman_Σ , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) | start_ARG italic_y , italic_α end_ARG ⟩ . (IV.40)

Notably, while j(y,E)𝑗𝑦𝐸j\left(y,E\right)italic_j ( italic_y , italic_E ) can be nonzero, the net current y=1Lj(y,E)superscriptsubscript𝑦1𝐿𝑗𝑦𝐸\sum_{y=1}^{L}j\left(y,E\right)∑ start_POSTSUBSCRIPT italic_y = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT italic_j ( italic_y , italic_E ) flowing between the bulk and edge always vanishes owing to the cyclicity of the trace:

y=1Lj(y,E)=Tr([Hedge,Gedge][Hedge,Gedge])=0,superscriptsubscript𝑦1𝐿𝑗𝑦𝐸tracesubscript𝐻edgesubscript𝐺edgesuperscriptsubscript𝐻edgesubscript𝐺edge0\displaystyle\sum_{y=1}^{L}j\left(y,E\right)=\Tr\left(-\big{[}H_{\mathrm{edge}% },G_{\mathrm{edge}}\big{]}-\big{[}H_{\mathrm{edge}},G_{\mathrm{edge}}\big{]}^{% \dagger}\right)=0,∑ start_POSTSUBSCRIPT italic_y = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT italic_j ( italic_y , italic_E ) = roman_Tr ( - [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] - [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) = 0 , (IV.41)

where the trace is taken over both sites y𝑦yitalic_y and the internal degrees α𝛼\alphaitalic_α of freedom. This is consistent with the absence of net current at equilibrium [113].

Additionally, the current density j(y,E)𝑗𝑦𝐸j(y,E)italic_j ( italic_y , italic_E ) vanishes under the periodic boundary conditions along the y𝑦yitalic_y direction. Owing to translation invariance, the Hamiltonians can be written as

Hedge=k;α,βhedgeαβ(k)|k,αk,β|,Heff=k;α,βheffαβ(k)|k,αk,β|,formulae-sequencesubscript𝐻edgesubscript𝑘𝛼𝛽superscriptsubscriptedge𝛼𝛽𝑘ket𝑘𝛼bra𝑘𝛽subscript𝐻effsubscript𝑘𝛼𝛽superscriptsubscripteff𝛼𝛽𝑘ket𝑘𝛼bra𝑘𝛽\displaystyle H_{\mathrm{edge}}=\sum_{k;\alpha,\beta}h_{\mathrm{edge}}^{\alpha% \beta}(k)\ket{k,\alpha}\bra{k,\beta},\quad H_{\mathrm{eff}}=\sum_{k;\alpha,% \beta}h_{\mathrm{eff}}^{\alpha\beta}(k)\ket{k,\alpha}\bra{k,\beta},italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_k ; italic_α , italic_β end_POSTSUBSCRIPT italic_h start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_β end_POSTSUPERSCRIPT ( italic_k ) | start_ARG italic_k , italic_α end_ARG ⟩ ⟨ start_ARG italic_k , italic_β end_ARG | , italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_k ; italic_α , italic_β end_POSTSUBSCRIPT italic_h start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_β end_POSTSUPERSCRIPT ( italic_k ) | start_ARG italic_k , italic_α end_ARG ⟩ ⟨ start_ARG italic_k , italic_β end_ARG | , (IV.42)

where k𝑘kitalic_k is momentum, and α𝛼\alphaitalic_α and β𝛽\betaitalic_β specify the internal degree of freedom. Then, from the cyclicity of the trace, we have

Trα[hedge(k),heff(k)]=0,subscripttrace𝛼subscriptedge𝑘subscripteff𝑘0\displaystyle\Tr_{\alpha}[h_{\mathrm{edge}}(k),h_{\mathrm{eff}}(k)]=0,roman_Tr start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT [ italic_h start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ( italic_k ) , italic_h start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( italic_k ) ] = 0 , (IV.43)

where the trace is taken only for the internal degree of freedom. Hence, we have Trα[Hedge,Heff]=0subscripttrace𝛼subscript𝐻edgesubscript𝐻eff0\Tr_{\alpha}[H_{\mathrm{edge}},H_{\mathrm{eff}}]=0roman_Tr start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_H start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ] = 0 and

j(y,E)=y|Trα([Hedge,Gedge]+[Hedge,Gedge])|y=0.𝑗𝑦𝐸bra𝑦subscripttrace𝛼subscript𝐻edgesubscript𝐺edgesuperscriptsubscript𝐻edgesubscript𝐺edgeket𝑦0\displaystyle j(y,E)=-\bra{y}\Tr_{\alpha}\left([H_{\mathrm{edge}},G_{\mathrm{% edge}}]+[H_{\mathrm{edge}},G_{\mathrm{edge}}]^{{\dagger}}\right)\ket{y}=0.italic_j ( italic_y , italic_E ) = - ⟨ start_ARG italic_y end_ARG | roman_Tr start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT ( [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] + [ italic_H start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT roman_edge end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) | start_ARG italic_y end_ARG ⟩ = 0 . (IV.44)

This is consistent with the absence of the local current under the periodic boundary conditions, as shown in Fig. 4 (c) of the main text.

V 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect in time-reversal-invariant topological insulators

We demonstrate the 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect in a time-reversal-invariant topological insulator. As a prototypical model, we study the Bernevig-Hughes-Zhang (BHZ) model [114] with the spin-orbit coupling,

HBHZ(kx,ky)=(tsin(kx))σzτx+(tsin(ky))τy+(m+tcos(kx)+tcos(ky))τz+Δσxτy,subscript𝐻BHZsubscript𝑘𝑥subscript𝑘𝑦𝑡subscript𝑘𝑥subscript𝜎𝑧subscript𝜏𝑥𝑡subscript𝑘𝑦subscript𝜏𝑦𝑚𝑡subscript𝑘𝑥𝑡subscript𝑘𝑦subscript𝜏𝑧Δsubscript𝜎𝑥subscript𝜏𝑦\displaystyle H_{\mathrm{BHZ}}\left(k_{x},k_{y}\right)=\left(t\sin{k_{x}}% \right)\sigma_{z}\tau_{x}+\left(t\sin{k_{y}}\right)\tau_{y}+\left(m+t\cos{k_{x% }}+t\cos{k_{y}}\right)\tau_{z}+\Delta\sigma_{x}\tau_{y},italic_H start_POSTSUBSCRIPT roman_BHZ end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) = ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) ) italic_σ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + ( italic_t roman_sin ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_τ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + ( italic_m + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG ) + italic_t roman_cos ( start_ARG italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG ) ) italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT + roman_Δ italic_σ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_τ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , (V.1)

where σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s and τisubscript𝜏𝑖\tau_{i}italic_τ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT’s (i=x,y,z𝑖𝑥𝑦𝑧i=x,y,zitalic_i = italic_x , italic_y , italic_z) are Pauli matrices, and t,m,Δ𝑡𝑚Δt,m,\Delta\in\mathbb{R}italic_t , italic_m , roman_Δ ∈ blackboard_R are real parameters. The BHZ model preserves time-reversal symmetry 𝒯H(𝒌)𝒯1=H(𝒌)𝒯𝐻𝒌superscript𝒯1𝐻𝒌\mathcal{T}H\left(\bm{k}\right)\mathcal{T}^{-1}=H\left(-\bm{k}\right)caligraphic_T italic_H ( bold_italic_k ) caligraphic_T start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = italic_H ( - bold_italic_k ) with 𝒯=iσy𝒦𝒯isubscript𝜎𝑦𝒦\mathcal{T}=\text{i}\sigma_{y}\mathcal{K}caligraphic_T = i italic_σ start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT caligraphic_K (𝒯2=1superscript𝒯21\mathcal{T}^{2}=-1caligraphic_T start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - 1) and belongs to class AII, to which the 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT topological invariant is assigned. We apply the open boundary conditions along the x𝑥xitalic_x direction and regard the edge at x=1𝑥1x=1italic_x = 1 as a system and the remaining bulk as an environment. In Fig. S1, we show the right eigenstates of the effective non-Hermitian Hamiltonian under the open boundary conditions along the y𝑦yitalic_y direction. We find that the Kramers pairs of the right eigenstates are localized at the opposite edges, confirming the 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect.

Refer to caption
Figure S1: 2subscript2\mathbb{Z}_{2}blackboard_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT skin effect in a time-reversal-invariant topological insulator (t=1.0𝑡1.0t=1.0italic_t = 1.0, m=1.6𝑚1.6m=-1.6italic_m = - 1.6, Δ=0.3Δ0.3\Delta=0.3roman_Δ = 0.3). The amplitude of the right eigenstates of the effective non-Hermitian Hamiltonian are shown.