## Non-Hermitian Topological Anderson Insulators

Dan-Wei Zhang , 1 ∗ Ling-Zhi Tang , Li-Jun Lang , Hui Yan , and Shi-Liang Zhu 1 1 1 2 1 , † 1 Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, GPETR Center for Quantum Precision Measurement and SPTE,

South China Normal University, Guangzhou 510006, China and

2 National Laboratory of Solid State Microstructures and School of Physics, Nanjing University, Nanjing 210093, China (Dated: April 15, 2020)

Non-Hermitian systems can exhibit exotic topological and localization properties. Here we elucidate the non-Hermitian effects on disordered topological systems using a nonreciprocal disordered Su-Schrieffer-Heeger model. We show that the non-Hermiticity can enhance the topological phase against disorders by increasing bulk gaps. Moreover, we uncover a topological phase which emerges under both moderate non-Hermiticity and disorders, and is characterized by localized insulating bulk states with a disorder-averaged winding number and zero-energy edge modes. Such topological phases induced by the combination of non-Hermiticity and disorders are dubbed non-Hermitian topological Anderson insulators . We reveal that the system has unique non-monotonous localization behavior and the topological transition is accompanied by an Anderson transition. These properties are general in other non-Hermitian models.

Topological states of matter have been widely explored in condensed-matter materials [1-5] and various engineered systems, which include ultracold atoms [6-8], photonic lattices [9, 10], mechanic systems [11], classic electronic circuits [12-15], and superconducting circuits [1620]. One hallmark of topological insulators is the robustness of nontrivial boundary states against certain types of weak disorders, since the topological band gap (topological invariants) preserves under these perturbations [13]. However, the band gap closes for sufficiently strong disorders and the system becomes trivial as all states are localized according to the Anderson localization [21]. Surprisingly, there is a topological phase driven from a trivial phase by disorders, known as topological Anderson insulator (TAI) [22]. The TAI was first predicted in two-dimensional (2D) quantum wells and then was shown to exhibit in a wide range of systems [22-30], such as SuSchrieffer-Heeger (SSH) chains [31]. Recently, the TAI has been observed in one-dimensional (1D) cold atomic wires and 2D photonic waveguide arrays [32, 33].

On the other hand, recent advances in non-Hermitian physics show that non-Hermitian systems have many intriguing features and applications [34-37]. Growing efforts are made to study topological properties of nonHermitian systems [38-73], which include new topological invariants [66], the non-Hermitian skin effect [48], the revised bulk-edge correspondence [47-53], and gainand-loss induced topological phases [57]. Non-Hermitian systems can exhibit unique localization properties in the presence of disorders [74-81]. Notably, the topological phases have been studied in 1D non-Hermitian AubryAndr´ e-Harper model [82-87], which describes topological quasicrystals and can be mapped to a 2D quantum Hall system [88-90]. The topological non-Hermitian quasicrystals were predicted [86, 87] and the topological nature of Anderson transitions in the systems was revealed [84, 85]. However, the interplay among topology, disorder

FIG. 1. (Color online) Sketch of the non-Hermitian disordered SSH model. The dotted box is the unit cell, and m j and t ( l,r ) j are Hermitian intracell and nonreciprocal intercell hoppings.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-c10cbe79-be21-4410-b1ee-423c03201273.jpg)

and non-Hermiticity can induce rich physical phenomena that have been rarely explored, in particular, the TAI phase has not been revealed in non-Hermitian systems.

In this Letter, we elucidate the non-Hermitian effects on disordered topological systems and uncover an exotic TAI phase. We construct an SSH model with nonreciprocal and disordered hopping terms and propose a realspace winding number to characterize its topology. Our main results are as follows: (i) We show that the nonHermiticity can enhance the topological phase against disorders under certain circumstances due to the increase of bulk gaps. (ii) We uncover a topological phase induced by combination of non-Hermiticity and disorders, which is dubbed non-Hermitian topological Anderson insulator (NHTAI) and characterized by localized insulating bulk states with a nontrivial winding number and zero-energy edge modes. (iii) We reveal the unique non-monotonous localization behavior owing to its interplay with the skin effect and the topological transition accompanied by an Anderson transition in the disordered system. (iv) We illustrate that the non-Hermitian enhancement effect and the NHTAI are general in other non-Hermitian models.

Model and topological invariant.Let us begin by considering the SSH model with nonreciprocal and disordered hoppings, which is depicted in Fig. 1. The tightbinding model with two-site unit cell reads

$$\underset { \text{der} } { \ n a } ^ { \ } _ { \text{led} } H = \sum _ { j } ( m _ { j } a _ { j } ^ { \dagger } b _ { j } + h. c. ) + t _ { j } ^ { ( r ) } a _ { j + 1 } ^ { \dagger } b _ { j } + t _ { j } ^ { ( l ) } b _ { j } ^ { \dagger } a _ { j + 1 }, \ \ ( 1 )$$

where a † j ( b † j ) creates a particle on the sublattice site A (B) in the j -th lattice cell, and a j ( b j ) is the corresponding annihilation operator. Here m j denotes the j -dependent (Hermitian) intracell hopping energy, and t ( r,l ) j characterize the non-Hermitian intercell hoppings. This Hamiltonian has the chiral symmetry as H satisfies CHC -1 = -H , where the chiral operator is C = σ z ⊗ I with the Pauli matrix σ z referring to the sublattice and the identity matrix I acting on the lattice cell.

In contrast to the site-potential disorder, the pure tunneling disorder is crucial for preserving the chiral symmetry. In particular, we consider the hopping terms as

$$m _ { j } = t + W _ { 1 } \omega _ { j }, \ t _ { j } ^ { ( l ) } = t ^ { \prime } + W _ { 2 } \omega _ { j } ^ { \prime }, \ t _ { j } ^ { ( r ) } = t _ { j } ^ { ( l ) } + f ( \gamma ). \ ( 2 ) \quad \ \. _ { 1 } |$$

Here t and t ′ are intracell and intercell tunneling energies, ω j and ω ′ j are independent random numbers chosen uniformly in the range [ -1 1], , W 1 and W 2 are the disorder strengths. Unless mentioned otherwise, we hereafter assume the nonreciprocal term f ( γ ) = t γ ′ with the nonHermiticity parameter γ . We set t = 1 as the energy unit and focus on the case of W 2 = 0.

Similar as the topological phase in the SSH model with t ′ &gt; t , this non-Hermitian disordered SSH model is topologically characterized by zero-energy edge modes and the corresponding winding number [see Eq. (4)]. Note that the system recovers to the Hermitian disordered SSH chain when γ = 0 [26, 27, 32] and to the non-Hermitian clean chain when W 1 = W 2 = 0 [47-50], respectively. In the clean limit, the topological invariant of the nonreciprocal SSH model can be a non-Bloch winding number in complex momentum space [48] or a dual open-bulk winding number in real space [49].

We now generalize the open-bulk winding number to our non-Hermitian disordered SSH model. Given a disorder configuration denoted by s , we diagonalize the Hamiltonian (1) under open boundary conditions (OBCs) with two chiral-symmetric parts: H nR s | ± 〉 s = ± E n,s | nR ± 〉 s with | nR -〉 s = C nR | + 〉 s . In the biorthonormal basis, the corresponding left eigenstates | nL ± 〉 orthonormal to the right eigenstates can be taken from the columns of ( T -1 s ) † by writing H s = T s Λ s T -1 s with Λ s diagonal. The homotopically equivalent flat band version of the Hamiltonian H s under OBCs is the open-boundary Q s matrix, which is given by Q s = ∑ n ( | nR + 〉 ss 〈 nL + | - | nR -〉 ss 〈 nL -| ). Here the summing takes over the eigenstates in the bulk spectrum without edge modes. The open-bulk winding number in real space is then defined as

$$\nu _ { s } = \frac { 1 } { 2 L ^ { \prime } } \text{Tr} ^ { \prime } ( C Q _ { s } [ Q _ { s }, X ] ), \text{ \quad \ \ } ( 3 ) \text{ \quad \ \ } \text{the}$$

where X is the coordinate operator, L = L ′ + 2 l is the chain length with three intervals of lengthes l, L , l ′ , and Tr ′ denotes the trace over the middle interval of length L ′ [49]. Here ν s serves well for disordered systems as it does not require the translation invariance and is quantized to

FIG. 2. (Color online) (a) Disorder-averaged winding number ν as a function of non-Hermiticity γ and disorder strength W . (b,c) Four middle disorder-averaged energy eigenvalues E n (real and imaginary parts) under OBCs, as a function of W for γ = 1. From bottom to up are n = L/ 2 -1 , L/ 2 , L/ 2+ 1 , L/ 2 + 2, respectively. (d) The bulk gap E g as a function of γ from simulations (hollow dots) and Eq. (5) (sold line). (e) ν as a function of γ and W for a nonreciprocal term f ( γ ) = t ′ (1 -γ + γ 2 ) with t ′ = . Other parameters are t t = 1, t ′ = 1 2 . t [in (a-d)], L = 5 l = 100, W 1 = W W , 2 = 0 and N s = 200.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-b92e5c82-9796-4ee0-8bab-245055cc6876.jpg)

an integer in the limit of L →∞ [27, 91]. We can define the disorder-averaged winding number

$$\nu = \frac { 1 } { N _ { s } } \sum _ { s = 1 } ^ { N _ { s } } \nu _ { s },$$

where a modest configuration number N s suffices in practice. For γ = 0, the topological invariants in Eqs. (3) and (4) reduce to those in Hermitian systems [27, 32, 91-94], where the boundary condition is irrelevant. However, Eq. (3) correctly corresponds to the topological edge modes only under OBCs, due to the unconventional bulk-edge correspondence in the non-Hermitian cases [47-52].

Enhancing topological phase by non-Hermiticity.We first consider the non-Hermitian effect on the topological phase of the open SSH chain with t ′ &gt; t under disorders. Figure 2(a) shows the disorder-averaged winding number ν as a function of γ and disorder strengths W 1 = W and W 2 = 0. We can find that when γ varies from 0 to 3 5, the topological phase with . ν glyph[similarequal] 1 preserves from a region W/t glyph[lessorsimilar] 2 to a larger one W/t glyph[lessorsimilar] 5. We rewrite the eigenequation of the open chain H ψ s | n s 〉 = E ( s ) n | ψ n s 〉 with wave functions | ψ n s 〉 = [ ψ ( s ) n, 1 , ψ ( s ) n, 2 , ...ψ ( s ) n,x , ...ψ ( s ) n,L ] T and eigenenergies E ( s ) n , where x is the lattice site index. Then the disorder-averaged eigenenergies is given by E n = 1 N s ∑ N s s =1 E ( s ) n . Figures 2(b,c) show the real and imaginary parts of E n for four center eigenstates as a function of W , respectively. The energy spectrum of this

open chain is purely real and two zero-energy edge modes exhibit in the topological phase.

To understand the numerical results, we take a similarity transformation [48]: ˜ H s = S -1 H S s with the diagonal matrix S = diag(1 1 , , r, r, r 2 , r 2 , ..., r L/ 2 -1 , r L/ 2 -1 ), the eigenequation is equivalent to ˜ H ψ s | ˜ n s 〉 = E ( s ) n | ˜ ψ n s 〉 with | ˜ ψ n s 〉 = S -1 | ψ n s 〉 . Let r = √ 1 + γ for γ &gt; -1, then ˜ H becomes the Hermitian disordered SSH model with intracell and intercell hoppings ˜ m j = m j and ˜ = t ′ t ′ √ 1 + γ for W 2 = 0. This transformation indicates that all eigenenergies of the open-chain Hamiltonian are real [Figs. 2(b,c)]. The transformation also accumulates the wave functions of bulk states to one boundary, which is the non-Hermitian skin effect [48] and will be discussed later. We can numerically calculate the bulk gap of the open clean system E g = | E L/ 2+2 -E L/ 2 -1 | for W = 0. Furthermore, the bulk gap for the open non-Hermitian chain can be derived from the Hermitian SSH chain after the similarity transformation, i.e.,

$$\tilde { E } _ { g } = 2 | \tilde { t } ^ { \prime } - \tilde { t } | = 2 | t ^ { \prime } \sqrt { 1 + \gamma } - t |. \quad \ \ ( 5 ) \quad \underset { \text{aver} } { \nu \nu }$$

The numerical results of E g for L = 100 are consistent with the values of ˜ E g given by Eq. (5), as shown in Fig. 2(d). Thus, the enhancement of the topological phase in this case can be interpreted as the increase of the effective bulk gap by non-Hermiticity.

The enhancement effect would be broken for negative γ or large W 2 , however, it can preserve without a similarity transformation, such as in the cases with small W 2 or an additional intercell hopping [95]. Surprisingly, we find that the nonreciprocal hopping (the gain and loss) can induce a topological phase when the corresponding Hermitian limit of γ = 0 is a trivial or critical phase for t ′ ≤ t . The results for the critical case of t = t ′ are shown in Fig. 2(e) for f ( γ ) = t (1 -γ + γ 2 ) (in Fig. 9 for the disordered SSH model with non-Hermitian gain and loss [95]). From Fig. 2(e), we also find the combination of non-Hermiticity and disorders can give rise to the topological phase.

NHTAI from non-Hermiticity and disorder.We further study the effects of combined non-Hermiticity and disorder on an initially trivial phase in the SSH chain with t &gt; t ′ . We numerically calculate ν as a function of γ and disorder strengths W 1 = W and W 2 = 0, with the results for t ′ = 0 7 . t and L = 400 shown in Fig. 3(a). In the Hermitian and clean limit γ = W = 0, the system is in the trivial phase with ν glyph[similarequal] 0. Interestingly, we find ν ≈ 1 in a region with moderate non-Hermiticity and disorder strength in Fig. 3(a). Actually, ν can approach to unit in this region by increasing the lattice size, with an example as a function of W with γ = 0 6 shown in . Fig. 3(c). The four middle disorder-averaged eigenenergies of the open chain of L = 100 are plotted in Fig. 3(d), which shows that two zero-energy edge modes inside a small gap between upper and lower eigenstates corresponding to the topological phase. Note that the small

𝛾

FIG. 3. (Color online) The winding number ν as a function of W and γ (a) or ˜ t ′ (b) for L = 400. (c) ν as a function of W for fixed γ = 0 6 and . L = 100 200 400. (d) Four disorder-, , averaged eigenenergies in the center of the energy spectrum E n = Re[ E n ] under OBCs, as a function of W for γ = 0 6 . and L = 100. From bottom to up are n = L/ 2 -1 , L/ 2 , L/ 2+ 1 , L/ 2 + 2, respectively. Other parameters are t ′ = 0 7 , . t l = 0 2 . L W , 1 = W W , 2 = 0, and N s = 50.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-5af26156-2eec-4de5-9b0c-630a30cc7fe5.jpg)

gap essentially vanish when L →∞ for moderate disorders. It is excepted that the energy gap is replaced by a mobility gap and the band insulator becomes an Anderson insulator with the topology carried by localized bulk states, similar as those for the Hermitian disordered chiral wires [27]. Our results indicate that a topological insulator can be induced from a trivial phase in the Hermitian and clean limit by the combination of moderate non-Hermiticity and disorders, which is dubbed NHTAI as a non-Hermitian extension of the TAI [22-27, 32].

To show the connection between NHTAIs and TAIs, we perform the similarity transformation and map the nonHermitian open SSH chain to the Hermitian one with ˜ t ′ = t ′ √ 1 + γ . We calculate the corresponding winding number for the Hermitian Hamiltonian (where the right eigenstates are used) [27, 32]. The results for an open chain of L = 400 shown in Fig. 3(b) indicates the ˜ t ′ W region of the TAI for 0 7 . &lt; t /t &lt; ˜ ′ 0 97, which corresponds . well to the γ W -region of the NHTAI for 0 &lt; γ &lt; 0 95. . Thus, NHTAIs can be topologically connected to TAIs through the similarity transformation with the same energy spectrum under OBCs. Notably, the NHTAI has unique properties without Hermitian counterparts, such as unconventional bulk-edge correspondence and nonmonotonous localization behavior under OBCs as we will discuss below. Moreover, the NHTAI can exist without the similarity transformation when the energy spectrum is generally complex [95].

Localization properties.We now investigate the localization properties of the disordered system. The den-

sity distributions of the L/ 2-th eigenstate ( L = 100) for γ = 0 6 and . W/t = 0 1 2 3 in a disorder configuration , . , are shown in Figs. 4(a-c), respectively. In the clean limit [Fig. 4(a)], this eigenstate is a bulk state, which is extended under PBCs but pinned to the right edge (because t ( r ) j &gt; t ( ) l j ) of the lattice under OBCs, as a manifestation of the non-Hermitian skin effect [48]. For moderate disorder strength, this eigenstate becomes a zero-energy mode localized at the right edge of the system that is in the NHTAI phase. The disorder-averaged density distribution of the edge state is shown in the inset of Fig. 4(b). For large disorder strength [Fig. 4(c)], it becomes a localized state in the bulk and the skin effect is broken. The same results are for the ( L/ 2 + 1)-th eigenstate due to the chiral symmetry. We calculate the disorder-averaged inverse participation ratio for the n -th eigenstate I n and its averaging over all eigenstates ¯ , I which are given by

$$\mathcal { I } _ { n } = \frac { 1 } { N _ { s } } \sum _ { s = 1 } ^ { N _ { s } } \sum _ { x = 1 } ^ { L } | \psi _ { n, x } ^ { ( s ) } | ^ { 4 }, \ \bar { \mathcal { I } } = \frac { 1 } { L } \sum _ { n = 1 } ^ { L } \mathcal { I } _ { n }, \quad ( 6 ) \quad \text{fig.\ } \text{th.e} \quad \\ \text{der c}.$$

glyph[negationslash]

respectively. We find that I n &gt; 1 /L for all bulk states with E n = 0 when W &gt; 0, which implies that the entire energy spectrum (excluded E = 0) of the chiral chain are localized immediately after the disorder is turned on [27]. The results for ¯ and I I L/ 2 with L = 400 are shown in Fig. 4(d). Under OBCs, I L/ 2 glyph[greatermuch] 1 /L for all W indicates that the L/ 2-th eigenstate always localized in this case; and I L/ 2 with a maximum value in the NHTAI phase implies that the zero-energy edge modes are more localized. The global localization index ¯ increasing rapidly as a funcI tion of W under PBCs also shows that the bulk states are localized in the presence of disorders. Remarkably, ¯ varies non-monotonously when I W/t glyph[lessorsimilar] 1 under OBCs, owing to the interplay of Anderson localization and the skin effect. When W/t glyph[greaterorsimilar] 1, ¯ takes nearly the same valI ues for OBCs and PBCs, which indicates that the skin effect of bulk states is now destroyed by disorders.

Following Ref. [27], we derive an analytical formula for the localization length Λ of our model (with W 2 = 0) at energy E = 0 under OBCs. We obtain [95]

$$\Lambda ^ { - 1 } = \left | \ln \left [ \frac { 2 e t ^ { \prime } \sqrt { 1 + \gamma } | 2 - 2 W | ^ { \frac { 1 } { 2 W } - \frac { 1 } { 2 } } } { | 2 + 2 W | ^ { \frac { 1 } { 2 W } + \frac { 1 } { 2 } } } \right ] \right |, \quad ( 7 ) \quad \text{the } \begin{array} { c } \text{the } \colon \\ \text{topo} \\ \text{inter} \end{array}$$

where e is the natural constant. Figure. 4(e) shows the critical points with Λ -1 → 0 in the γ W -plane, where the localization length of zero-energy states diverges. The delocalized critical points match with the topological transition points in Fig. 3(a), indicating that the topological transition is accompanied by an Anderson localizationdelocalization transition. Our results reveal that the nonHermitian topological numbers can be carried by localized bulk states and disorders can drive a localized nonHermitian topological phase through a delocalized point, similar as those in Hermitian chiral chains [27].

ҧ

FIG. 4. (Color online) (a,b,c) Density distribution of the L/ 2th eigenstate with L = 100 for W/t = 0 1 2 3 in a disor-, . , der configuration, respectively. The inset figure in (b) is the disorder-averaged density distribution. (d) ¯ and I I L/ 2 (inset) with L = 400, respectively. (e) Λ -1 as a function of γ and W . The same parameters as in Fig. 3 and γ = 0 6 in (a-d). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-5ddaeb1c-ac36-4356-87f7-923974f26e7b.jpg)

Discussion and conclusion.We discuss the realization of the non-Hermitian disordered SSH model in some artificial systems. The first feasible system is ultracold atoms [6-8]. The TAI has been realized in an atomic SSH wire with controllable disordered hoppings [32], where the effective nonreciprocal hopping can be engineered by a collective one-body loss [54, 96]. Another system is photonic crystals [9, 10], where the TAI phase and tunable non-Hermiticity were achieved [33, 97-100]. The realization of nonreciprocal hoppings in optics was suggested [101, 102]. The third system is topological electronic circuits [12-15], where the Hermitian and nonreciprocal SSH chains were realized [103-105] and tunable hopping disorders can be added. In view that non-Hermiticity and disorders have been engineered in these artificial systems, the studied model with the NHTAI phase is realizable in current experiments. The robust or disorder-induced topological edge modes can be detected, and it would be interesting to measure the topological numbers.

glyph[negationslash]

We note that our findings are applicable for other nonHermitian models. The non-Hermitian enhancement effect and the NHTAI phase can preserve without the similarity transformation, in which cases the energy spectra are generally complex. For instance, when W 2 = 0, the similarity transformation is inapplicable but the main results obtained for W 2 = 0 preserve for W 2 glyph[lessorsimilar] W/ 4, as shown in Supplemental Material [95]. Another example is a modified model with an additional intercell hopping term [95]. The enhancement effect and the NHTAI phase also exhibit in a disordered SSH model with non-

Hermitian gain and loss [95]. The NHTAI phase even exists when the non-Hermiticity becomes random, although the skin effect is generally broken in this case [95].

In summary, we have explored the topological and localization properties of the SSH model with nonreciprocal and disordered hoppings. We have revealed the enhancement of the topological phase by non-Hermiticity and the NHTAI induced by combination of non-Hermiticity and disorders. The non-monotonous localization behavior and the topological nature of the Anderson transition have been elucidated. Moreover, the predicted NHTAI can be experimentally realized in some artificial systems.

This work was supported by the NKRDP of China (Grant No. 2016YFA0301800), the NSFC (Grants No. 11604103, No. 11704367, No. 11904109, and No. 91636218), the NSAF (Grant No. U1830111 and No. U1801661), the Key-Area Research and Development Program of GuangDong Province (Grant No. 2019B030330001), and the Key Program of Science and Technology of Guangzhou (Grant No. 201804020055).

- ∗ danweizhang@m.scnu.edu.cn † slzhu@nju.edu.cn
- [1] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82 , 3045 (2010).
- [2] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83 , 1057 (2011).
- [3] A. Bansil, H. Lin, and T. Das, Rev. Mod. Phys. 88 , 021004 (2016).
- [4] N. P. Armitage, E. J. Mele, and A. Vishwanath, Rev. Mod. Phys. 90 , 015001 (2018).
- [5] C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Rev. Mod. Phys. 88 , 035005 (2016).
- [6] D.-W. Zhang, Y.-Q. Zhu, Y. X. Zhao, H. Yan, and S.-L. Zhu, Advances in Physics 67 , 253 (2018).
- [7] N. R. Cooper, J. Dalibard, and I. B. Spielman, Rev. Mod. Phys. 91 , 015005 (2019).
- [8] N. Goldman, J. C. Budich, and P. Zoller, Nat. Phys. 12 , 639 (2016).
- [9] L. Lu, J. D. Joannopoulos, and M. Soljaˇ ci´, c Nature Photonics 8 , 821 (2014).
- [10] T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zilberberg, and I. Carusotto, Rev. Mod. Phys. 91 , 015006 (2019).
- [11] S. D. Huber, Nature Physics 12 , 621 (2016).
- [12] J. Ningyuan, C. Owens, A. Sommer, D. Schuster, and J. Simon, Phys. Rev. X 5 , 021031 (2015).
- [13] V. V. Albert, L. I. Glazman, and L. Jiang, Phys. Rev. Lett. 114 , 173902 (2015).
- [14] S. Imhof, C. Berger, F. Bayer, J. Brehm, L. W. Molenkamp, T. Kiessling, F. Schindler, C. H. Lee, M. Greiter, T. Neupert, and R. Thomale, Nature Physics 14 , 925 (2018).
- [15] C. H. Lee, S. Imhof, C. Berger, F. Bayer, J. Brehm, L. W. Molenkamp, T. Kiessling, and R. Thomale, Communications Physics 1 , 39 (2018).
- [16] M. D. Schroer, M. H. Kolodrubetz, W. F. Kindel,
- M. Sandberg, J. Gao, M. R. Vissers, D. P. Pappas, A. Polkovnikov, and K. W. Lehnert, Phys. Rev. Lett. 113 , 050402 (2014).
- [17] P. Roushan, C. Neill, Y. Chen, M. Kolodrubetz, C. Quintana, N. Leung, M. Fang, R. Barends, B. Campbell, Z. Chen, B. Chiaro, A. Dunsworth, E. Jeffrey, J. Kelly, A. Megrant, J. Mutus, P. J. J. O'Malley, D. Sank, A. Vainsencher, J. Wenner, T. White, A. Polkovnikov, A. N. Cleland, and J. M. Martinis, Nature 515 , 241 (2014).
- [18] X. Tan, D.-W. Zhang, Q. Liu, G. Xue, H.-F. Yu, Y.-Q. Zhu, H. Yan, S.-L. Zhu, and Y. Yu, Phys. Rev. Lett. 120 , 130503 (2018).
- [19] X. Tan, Y. X. Zhao, Q. Liu, G. Xue, H.-F. Yu, Z. D. Wang, and Y. Yu, Phys. Rev. Lett. 122 , 010501 (2019).
- [20] X. Tan, D.-W. Zhang, Z. Yang, J. Chu, Y.-Q. Zhu, D. Li, X. Yang, S. Song, Z. Han, Z. Li, Y. Dong, H.-F. Yu, H. Yan, S.-L. Zhu, and Y. Yu, Phys. Rev. Lett. 122 , 210401 (2019).
- [21] P. W. Anderson, Physical Review 109 , 1492 (1958).
- [22] J. Li, R.-L. Chu, J. K. Jain, and S.-Q. Shen, Phys. Rev. Lett. 102 , 136806 (2009).
- [23] C. W. Groth, M. Wimmer, A. R. Akhmerov, J. Tworzydglyph[suppress]o, l and C. W. J. Beenakker, Phys. Rev. Lett. 103 , 196805 (2009).
- [24] H. Jiang, L. Wang, Q.-F. Sun, and X. C. Xie, Phys. Rev. B 80 , 165316 (2009).
- [25] H.-M. Guo, G. Rosenberg, G. Refael, and M. Franz, Phys. Rev. Lett. 105 , 216601 (2010).
- [26] A. Altland, D. Bagrets, L. Fritz, A. Kamenev, and H. Schmiedt, Phys. Rev. Lett. 112 , 206602 (2014).
- [27] I. Mondragon-Shem, T. L. Hughes, J. Song, and E. Prodan, Phys. Rev. Lett. 113 , 046802 (2014).
- [28] P. Titum, N. H. Lindner, M. C. Rechtsman, and G. Refael, Phys. Rev. Lett. 114 , 056801 (2015).
- [29] P. V. Sriluckshmy, K. Saha, and R. Moessner, Phys. Rev. B 97 , 024204 (2018).
- [30] Z.-Q. Zhang, B.-L. Wu, J. Song, and H. Jiang, arXiv:1906.04064v1.
- [31] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 42 , 1698 (1979).
- [32] E. J. Meier, F. A. An, A. Dauphin, M. Maffei, P. Massignan, T. L. Hughes, and B. Gadway, Science 362 , 929 (2018).
- [33] S. Sttzer, Y. Plotnik, Y. Lumer, P. Titum, N. H. Lindner, M. Segev, M. C. Rechtsman, and A. Szameit, Nature 560 , 461 (2018).
- [34] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80 , 5243 (1998).
- [35] L. Feng, R. El-Ganainy, and L. Ge, Nature Photonics 11 , 752 (2017).
- [36] R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, and D. N. Christodoulides, Nature Physics 14 , 11 (2018).
- [37] M.-A. Miri and A. Al`, Science u 363 , eaar7709 (2019).
- [38] M. S. Rudner and L. S. Levitov, Phys. Rev. Lett. 102 , 065703 (2009).
- [39] J. M. Zeuner, M. C. Rechtsman, Y. Plotnik, Y. Lumer, S. Nolte, M. S. Rudner, M. Segev, and A. Szameit, Phys. Rev. Lett. 115 , 040402 (2015).
- [40] K. Esaki, M. Sato, K. Hasebe, and M. Kohmoto, Phys. Rev. B 84 , 205128 (2011).
- [41] Y. C. Hu and T. L. Hughes, Phys. Rev. B 84 , 153101 (2011).

- [42] B. Zhu, R. L¨, u and S. Chen, Phys. Rev. A 89 , 062102 (2014).
- [43] S. Malzard, C. Poli, and H. Schomerus, Phys. Rev. Lett. 115 , 200402 (2015).
- [44] D. Leykam, K. Y. Bliokh, C. Huang, Y. D. Chong, and F. Nori, Phys. Rev. Lett. 118 , 040401 (2017).
- [45] Y. Xu, S.-T. Wang, and L.-M. Duan, Phys. Rev. Lett. 118 , 045701 (2017).
- [46] X. Zhan, L. Xiao, Z. Bian, K. Wang, X. Qiu, B. C. Sanders, W. Yi, and P. Xue, Phys. Rev. Lett. 119 , 130501 (2017).
- [47] T. E. Lee, Phys. Rev. Lett. 116 , 133903 (2016).
- [48] S. Yao and Z. Wang, Phys. Rev. Lett. 121 , 086803 (2018).
- [49] F. Song, S. Yao, and Z. Wang, arXiv:1905.02211v1.
- [50] F. K. Kunst, E. Edvardsson, J. C. Budich, and E. J. Bergholtz, Phys. Rev. Lett. 121 , 026808 (2018).
- [51] Y. Xiong, Journal of Physics Communications 2 , 035043 (2018).
- [52] L. Jin and Z. Song, Phys. Rev. B 99 , 081103 (2019).
- [53] D. S. Borgnia, A. J. Kruchkov, and R.-J. Slager, arXiv: 1902.07217v2.
- [54] Z. Gong, Y. Ashida, K. Kawabata, K. Takasan, S. Higashikawa, and M. Ueda, Phys. Rev. X 8 , 031079 (2018).
- [55] S. Yao, F. Song, and Z. Wang, Phys. Rev. Lett. 121 , 136802 (2018).
- [56] H. Shen, B. Zhen, and L. Fu, Phys. Rev. Lett. 120 , 146402 (2018).
- [57] K. Takata and M. Notomi, Phys. Rev. Lett. 121 , 213902 (2018).
- [58] Y. Chen and H. Zhai, Phys. Rev. B 98 , 245130 (2018).
- [59] L.-J. Lang, Y. Wang, H. Wang, and Y. D. Chong, Phys. Rev. B 98 , 094307 (2018).
- [60] G. Harari, M. A. Bandres, Y. Lumer, M. C. Rechtsman, Y. D. Chong, M. Khajavikhan, D. N. Christodoulides, and M. Segev, Science 359 , eaar4003 (2018).
- [61] M. A. Bandres, S. Wittek, G. Harari, M. Parto, J. Ren, M. Segev, D. N. Christodoulides, and M. Khajavikhan, Science 359 , eaar4005 (2018).
- [62] H. Zhou, C. Peng, Y. Yoon, C. W. Hsu, K. A. Nelson, L. Fu, J. D. Joannopoulos, M. Soljaˇi´, c c and B. Zhen, Science 359 , 1009 (2018).

[63]

T.-S. Deng and W. Yi, Phys. Rev. B

100

, 035102 (2019).

- [64] M. Ezawa, Phys. Rev. B 99 , 201411 (2019).
- [65] K. Kawabata, S. Higashikawa, Z. Gong, Y. Ashida, and M. Ueda, Nature Communications 10 , 297 (2019).
- [66] A. Ghatak and T. Das, Journal of Physics: Condensed Matter 31 , 263001 (2019).
- [67] K. Kawabata, Y. Ashida, H. Katsura, and M. Ueda, Phys. Rev. B 98 , 085116 (2018).
- [68] J.-Q. Cai, Q.-Y. Yang, Z.-Y. Xue, M. Gong, G.-C. Guo, and Y. Hu, arXiv:1812.02610v2.
- [69] T. Liu, Y.-R. Zhang, Q. Ai, Z. Gong, K. Kawabata, M. Ueda, and F. Nori, Phys. Rev. Lett. 122 , 076801 (2019).
- [70] C. H. Lee, L. Li, and J. Gong, Phys. Rev. Lett. 123 , 016805 (2019).
- [71] T. Yoshida, K. Kudo, and Y. Hatsugai, arXiv: 1907.07596v2.
- [72] D. J. Luitz and F. Piazza, arXiv:1906.02224v1.
- [73] K. Yamamoto, M. Nakagawa, K. Adachi, K. Takasan, M. Ueda, and N. Kawakami, Phys. Rev. Lett. 123 , 123601 (2019).
- [74] N. Hatano and D. R. Nelson, Phys. Rev. Lett. 77 , 570 (1996).
- [75] N. Hatano and D. R. Nelson, Phys. Rev. B 56 , 8651 (1997).
- [76] A. Carmele, M. Heyl, C. Kraus, and M. Dalmonte, Phys. Rev. B 92 , 195107 (2015).
- [77] C. Mej´ ıa-Cort´s e and M. I. Molina, Phys. Rev. A 91 , 033815 (2015).
- [78] E. Levi, M. Heyl, I. Lesanovsky, and J. P. Garrahan, Phys. Rev. Lett. 116 , 237203 (2016).
- [79] Q.-B. Zeng, S. Chen, and R. L¨, Phys. Rev. A u 95 , 062118 (2017).
- [80] R. Hamazaki, K. Kawabata, and M. Ueda, arXiv:1811.11319v2.
- [81] V. M. Martinez Alvarez, J. E. Barrios Vargas, and L. E. F. Foa Torres, Phys. Rev. B 97 , 121401 (2018).
- [82] P. G. Harper, Proc. Phys. Soc. Sect. A 68 , 874 (1955).
- [83] S. Aubry and G. Andre, Ann. Israel Phys. Soc. 3 , 133 (1980).
- [84] S. Longhi, Phys. Rev. Lett. 122 , 237601 (2019).
- [85] H. Jiang, L.-J. Lang, C. Yang, S.-L. Zhu, and S. Chen, Phys. Rev. B 100 , 054301 (2019).
- [86] Q.-B. Zeng, Y.-B. Yang, and Y. Xu, arXiv:1901.08060v2.
- [87] J. Hou, Y.-J. Wu, and C. Zhang, arXiv:1906.03988v1.
- [88] Y. E. Kraus, Y. Lahini, Z. Ringel, M. Verbin, and O. Zilberberg, Phys. Rev. Lett. 109 , 106402 (2012).
- [89] L.-J. Lang and S. Chen, Phys. Rev. B 86 , 205135 (2012).
- [90] S. Ganeshan, K. Sun, and S. Das Sarma, Phys. Rev. Lett. 110 , 180403 (2013).
- [91] A. Kitaev, Annals of Physics 321 , 2 (2006).
- [92] E. Prodan, T. L. Hughes, and B. A. Bernevig, Phys. Rev. Lett. 105 , 115501 (2010).
- [93] R. Bianco and R. Resta, Phys. Rev. B 84 , 241106 (2011).
- [94] Y. X. Zhao, Y. Lu, and H.-Z. Lu, arXiv:1706.09783v1.
- [95] See Supplemental Material for more details.
- [96] J. Li, A. K. Harter, J. Liu, L. de Melo, Y. N. Joglekar, and L. Luo, Nature Communications 10 , 855 (2019).
- [97] C. Poli, M. Bellec, U. Kuhl, F. Mortessagne, and H. Schomerus, Nature Communications 6 , 6710 (2015).
- [98] S. Weimann, M. Kremer, Y. Plotnik, Y. Lumer, S. Nolte, K. G. Makris, M. Segev, M. C. Rechtsman, and A. Szameit, Nature Materials 16 , 433 (2016).
- [99] M. Pan, H. Zhao, P. Miao, S. Longhi, and L. Feng, Nature Communications 9 , 1308 (2018).
- [100] S ¸ . K. zdemir, S. Rotter, F. Nori, and L. Yang, Nature Materials 18 , 783 (2019).
- [101] S. Longhi, EPL (Europhysics Letters) 106 , 34001 (2014).
- [102] S. Longhi, D. Gatti, and G. Della Valle, Phys. Rev. B 92 , 094204 (2015).
- [103] Y. Hadad, J. C. Soric, A. B. Khanikaev, and A. Al`, u Nature Electronics 1 , 178 (2018).
- [104] Y. Wang, L.-J. Lang, C. H. Lee, B. Zhang, and Y. D. Chong, Nature Communications 10 , 1102 (2019).
- [105] T. Helbig, T. Hofmann, S. Imhof, M. Abdelghany, T. Kiessling, L. W. Molenkamp, C. H. Lee, A. Szameit, M. Greiter, and R. Thomale, arXiv:1907.11562v1.

## Supplemental Material for Non-Hermitian Topological Anderson Insulators

## 1. Localization length at zero energy

For the Hermitian disordered chiral SSH chain, the localization length Λ at energy E = 0 can be analytically obtained [27]. Indeed, the Schr¨dinger equation of o the Hermitian SSH model ˜ ˜ Hψ = 0 reads: ˜ ˜ t j ψ j -α,α + ˜ m ψ j ˜ j,α = 0, where α = ± 1 represents A and B sublattice, respectively. The solution is given by

$$\tilde { \psi } _ { N, \alpha } = \prod _ { j = 1 } ^ { N } \left ( \frac { \tilde { t } _ { j } } { \tilde { m } _ { j } } \right ) \tilde { \psi } _ { 0, \alpha }, \quad \quad \ ( 8 ) \quad \text{FIG}$$

where the unit cell is labeled by j = 0 1 , , ..., N . The inverse of the localization length is given by (in the thermodynamic limit disregarding the boundaries) [27]

$$\Lambda ^ { - 1 } & = \max _ { \alpha = \pm 1 } \left [ - \lim _ { N \to \infty } \frac { 1 } { N } \ln | \tilde { \psi } _ { N, \alpha } | \right ] \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad = | \lim _ { N \to \infty } \frac { 1 } { N } \sum _ { j = 1 } ^ { N } ( \ln | \tilde { t } _ { j } | - \ln | \tilde { m } _ { j } | ) |. \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

We consider the disordered hopping parameters: ˜ m j = ˜ + t W ω 1 j and ˜ t j = ˜ + t ′ W ω 2 ′ j , where t and t ′ are the characteristic intracell and intercell tunneling energies, ω j and ω ′ j are independent random real numbers chosen uniformly from the range [ -1 1], , W 1 and W 2 are the corresponding disorder strengths, and ˜ = 1 is t set as the energy unit. Note that the notations used here are slightly different from those in Ref. [27].

According to Birkhoff's ergodic theorem, one can use the ensemble average to evaluate the last expression in Eq. (9), which is then given by

$$\Lambda ^ { - 1 } & = \left | \int _ { - 1 } ^ { 1 } d \omega \int _ { - 1 } ^ { 1 } d \omega ^ { \prime } \ ( \ln | \tilde { t } ^ { \prime } + W _ { 2 } \omega | - \ln | 1 + W _ { 1 } \omega ^ { \prime } | ) \right |. \\ \tau _ { \infty } \ \cdot \. \ \. \ \. \ \. \ \. \. \. \.. \.. \........................................................................................................................................................................................................$$

The integrations can be calculated explicitly, and the arguments of the logarithms can become negative in the regime of large W . One can obtain

$$\Lambda ^ { - 1 } = \left | \ln \left [ \frac { | 2 \tilde { t } ^ { \prime } + 2 W _ { 2 } | ^ { \frac { \tilde { t } ^ { \prime } } { 2 W _ { 2 } } + \frac { 1 } { 2 } } } { | 2 \tilde { t } ^ { \prime } - 2 W _ { 2 } | ^ { \frac { \tilde { t } ^ { \prime } } { 2 W _ { 2 } } - \frac { 1 } { 2 } } } \frac { | 2 - 2 W _ { 1 } | ^ { \frac { 1 } { 2 W _ { 1 } } - \frac { 1 } { 2 } } } { | 2 + 2 W _ { 1 } | ^ { \frac { 1 } { 2 W _ { 1 } } + \frac { 1 } { 2 } } } \right ] \right |. \ ( 1 1 ) \quad \text{(used),} \\ \text{$\text{$\text{$\text{$comp.}$} } }.$$

By using the similarity transformation ˜ = H S -1 HS on our non-Hermitian disordered SSH Hamiltonian H with W 2 = 0 under OBCs, one can the Hermitian disordered SSH Hamiltonian ˜ . H After the transformation, we can obtain the solution of the Schr¨dinger equation of o Hψ = 0 with ψ = Sψ ˜ , which replaces the solution in Eq. (8) with the form:

$$\psi _ { N, \alpha } = \prod _ { j = 1 } ^ { N } \left ( \frac { t _ { j } } { m _ { j } } \right ) \psi _ { 0, \alpha }. \quad \quad ( 1 2 ) \quad \text{top} \\ \text{and}$$

glyph[negationslash]

FIG. 5. (Color online) Λ -1 as a function of γ and W . The parameters are t = 1, t ′ = 1 2 , . t W 1 = W W , 2 = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-deade0a9-1463-433e-b471-02783d5c9d04.jpg)

With the mapping, one has the hopping parameters t j = t r ′ + W ω 1 j = t ′ √ 1 + γ + W ω 1 j and m j = ˜ m j = 1 for ˜ = t t = 1 and W 2 = 0. Thus, by substituting ˜ t ′ = t ′ √ 1 + γ , W 1 = W and W 2 = 0 into Eq. (11), we obtain the inverse of the localization length of zero-energy states in our non-Hermitian disordered SSH chain under OBCs,

$$\Lambda ^ { - 1 } = \left | \ln \left [ \frac { 2 e t ^ { \prime } \sqrt { 1 + \gamma } | 2 - 2 W | \frac { 1 } { 2 W } - \frac { 1 } { 2 } } } { | 2 + 2 W | \frac { 1 } { 2 W } + \frac { 1 } { 2 } } \right ] \right |, \quad ( 1 3 )$$

where e is the natural constant. Figure. 5 shows the critical points with Λ -1 → 0 in the γ W -plane, where the delocalized critical points match with the topological transition points in Fig. 2(a). This result, as well as Fig. 4(e) and Fig. 3(a) in the main text, demonstrates that the topological transition is accompanied by an Anderson localization-delocalization transition in our non-Hermitian disordered SSH model.

glyph[negationslash]

## 2. Results for the case of W 2 = 0

glyph[negationslash]

When the intercell hopping disorder strength W 2 = 0, one can not find a similarity transformation of the nonHermitian disordered SSH chain under OBCs for every disorder configuration (especially in the strong disorder case). In this case, the (disorder-averaged) energy spectrum can be generally complex. Figures. 6(a,b) show the complex energy spectrum for W 1 = 4 W 2 = W as a function of W in a disorder configuration. The corresponding disorder-averaged winding number ν as a function of W and γ is shown in Fig. 6(c), which is close to the result of W 2 = 0 shown in Fig. 2(a). Thus, the non-Hermitian enhancement of the topological phase can still exhibit when W 2 = 0. For comparisons, we plot the winding number ν as a function of W and ˜ t ′ in Fig. 6(d), corresponding to the Hermitian open SSH chain with ˜ t ′ = t ′ √ 1 + γ after the similarity transformation when W 2 = 0, where the topological regime matches well with those in Figs. 2(a) and 6(c). We note that the enhancement effect will be

FIG. 6. (Color online) (a,b) The energy spectrum E n (real and imaginary parts) under OBCs as a function of W with fixed γ = 1 in a disorder configuration with strengths W 1 = 4 W 2 = W . (c,e) Disorder-averaged winding number ν as a function of W W ( 2 = W/ 4 and W 2 = W/ 2) and γ ; (d) ν as a function of W W ( 2 = 0) and ˜ . t ′ Other parameters are t = 1, t ′ = 1 2 , . t L = 5 l = 100, and N s = 200.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-658183cf-0893-4cee-a52f-f7a6de0201da.jpg)

FIG. 7. (Color online) (a,b) ν as a function of W and γ . Other parameters are t = 1, t ′ = 0 7 , . t l = 0 2 . L , and N s = 200.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-03f85324-6ca0-4b18-96cb-f838c8c05d8f.jpg)

𝑊1 = 𝑊 2 = 𝑊 𝑊1 = 0, 𝑊 2 = 𝑊 broken by increasing W 2 , such as when W 2 glyph[greaterorsimilar] W/ 2 in this case, as shown in Fig. 6(e).

In Figs. 7(a,b), we show ν as a function of W and γ for W 1 = 4 W 2 = W and W 1 = 2 W 2 = W , respectively. Although we here adopt a lattice of length L = 100 smaller than that of L = 400 in Fig. 3, the topological regime with the NHTAI can be roughly seen. Tuning on the disorder of W 1 up to W 1 = W/ 4, the topological regime almost preserves, and the NHTAI phase remains. However, the topological region will be reduced when W 1 becomes larger, such as the case of W 1 = W/ 2 in Fig. 7(b) with a narrow parameter region for the NHTAI phase.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-02ea84f2-a2f2-4596-89ef-c265ec7986af.jpg)

## 3. Modified Hamiltonian without a similarity transformation

glyph[negationslash]

We consider a modified SSH Hamiltonian with nonreciprocal and disordered hoppings that cannot be mapped to a Hermitian Hamiltonian through the similarity transformation under OBCs. The modified model contains an additional intercell hopping term with the strength t ′′ , as depicted in Fig. 8(a). This hopping term preserves the chiral symmetry of the system, however, the similarity transformation is inapplicable when t ′′ = 0 (see a similar example without disorders in Ref. [48]) and thus the modified model in this case does not has a direct Hermitian counterpart.

glyph[negationslash]

Since the chiral symmetry remains when t ′′ = 0, we can still calculate the disordered-averaged winding num-

FIG. 8. (Color online) (a) Sketch of the disordered nonreciprocal SSH model with an additional hopping term t ′′ . (b,d) ν as a function of W and γ for t ′′ = 0 2 . t and t ′′ = 0 1 , re-. t spectively. (c) The energy spectrum E n (real and imaginary parts) under OBCs as a function of W for γ = 1 in a disorder configuration. (e) ν as a function of W for fixed γ = 0 5, . t ′′ = 0 1 , and . t L = 100 200 400. Other parameters are , , t = 1, L = 100, l = 0 2 . L , N s = 200, and t ′ = 1 2 . t in (b,c) and t ′ = 0 7 . t in (d,e).

ber ν . In Fig. 8(b), we show ν as a function of W ( W 1 = W and W 2 = 0) and γ for t ′′ = 0 2 . t ( t ′ = 1 2 . t and L = 100). We find that the non-Hermitian enhancement effect on the topological phase exhibits in this case, although the energy spectrum is generally complex since the similarity transformation is inapplicable. An example of the complex energy spectrum for a disorder configuration and γ = 1 is shown in Fig. 8(c). Moreover, we find that the NHTAI phase can exist without the simi-

FIG. 9. (Color online) (a,b,c) The disorder-averaged winding number ν as a function of the disorder strength W and the standard deviation of random non-Hermiticity σ γ for t ′ = 2 , t t ′ = , and t t ′ = 0 7 , respectively. . t (d) The averaged density distribution of the L/ 2-th eigenstate for σ γ = 0 1 1 4 with . , , fixed W 1 = W 2 = 0 and t ′ = 0 7 . . t Other parameters are t = 1, L = 100, l = 0 2 . L W , 1 = W W , 2 = 0, and N s = 1000.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-6fa3b84a-2315-4247-9057-f4297964726c.jpg)

glyph[negationslash]

larity transformation when t ′′ = 0, with an example for t ′′ = 0 1 . t ( t ′ = 0 7 . t and L = 100) shown in Fig. 8(c). To show the parameter regime for the NHTAI phase more clearly, we plot ν as a function of W for γ = 0 5 . t and L = 100 200 400 in Fig. , , 8(e). We note that the enhancement effect and the NHTAI phase will be broken for large t ′′ .

## 4. Results for the case of random non-Hermiticity

In this section, we consider the effect of random nonHermiticity on the topological phenomena described in the main text. The hopping terms can be rewritten as

$$m _ { j } = t + W _ { 1 } \omega _ { j }, \ t _ { j } ^ { ( l ) } = t ^ { \prime } + W _ { 2 } \omega _ { j } ^ { \prime }, \ t _ { j } ^ { ( r ) } = t _ { j } ^ { ( l ) } + t \gamma _ { j }. \ ( 1 4 ) \quad \underset { \ell \ell \ell } { \ m o m e }$$

Here t = 1 is set as the energy unit, ω j and ω ′ j are random numbers chosen uniformly in the range [ -1 1] with , the disorder strengths W 1 and W 2 , and the nonreciprocal parameters for the j -th lattice cell γ j are independent random numbers satisfying the Gaussian distribution N (0 , σ 2 γ ) with the standard deviation σ γ . For simplicity, we focus on the case of W 2 = 0 and set W 1 = W .

In the presence of hopping disorder and random nonreciprocal parameter, we numerically find that this kind of non-Hermiticity cannot enhance the topological phase in the general case when t ′ &gt; t , with an typical example shown in Fig. 9(a) for t ′ = 2 . t However, the random non-Hermiticity still can induce and enhance the topological phase from the critical case of t ′ = t , as shown in Fig. 9(b), similar as the results shown in Fig. 2(e) in the main text. Moreover, the NHTAI phase still exists in under this condition by the combination of hopping disorders and random non-Hermiticity from the trivial phase in the clean and Hermitian limit ( t ′ &lt; t ), as shown in Fig. 9(c) for t ′ = 0 7 . . t Notably, under the random non-Hermiticity γ j with the zero mean value, the nonHermitian skin effect of the bulk states in the open chain is generally broken after averaging over many configurations. As shown in Fig. 9(d), we plot the averaged density distribution of the L/ 2-th eigenstate for σ γ = 0 1 1 4 . , , with fixed W 1 = W 2 = 0 and t ′ = 0 7 . . t It is clear that when σ γ = 0 1 1, this bulk state is extended in the lat-. , tice under OBCs, which indicates the breaking of the non-Hermitian skin effect. Other extended bulk states exhibit similar density distributions. When σ γ = 4, the L/ 2-th ( L/ 2 + 1-th) eigenstate becomes a zero-energy mode localized near the left (right) edge of the system that is in the topological phase.

## 5. Disordered SSH model with gain and loss

We further consider another disordered SSH-type model with non-Hermitian gain and loss, as depicted in Fig. 10(a). The model Hamiltonian is given by

$$\infty \quad \text{$-\circ$} \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \\ \text{ore} \quad H _ { 2 } & = \sum _ { j } ( m _ { j } a _ { j } ^ { \dagger } b _ { j } + h. c. ) + \sum _ { j } \frac { i \Gamma } { 2 } ( a _ { j } ^ { \dagger } a _ { j } - b _ { j + 1 } ^ { \dagger } b _ { j } ) + \frac { t ^ { \prime } } { 2 } a _ { j + 1 } ^ { \dagger } b _ { j } + h. c. ], \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

with the gain-and-loss parameter Γ and the random hopping strengths

$$m _ { j } = t + W _ { 1 } \omega _ { j }, \ i t ^ { \prime } _ { j } = i t ^ { \prime } + W _ { 2 } \omega ^ { \prime } _ { j }. \quad \ ( 1 6 )$$

Here ω j and ω ′ j are independent random real numbers chosen uniformly in the range [ -1 1], , W 1 and W 2 are the corresponding disorder strengths. In the clean limit W 1 = W 2 = 0, the corresponding Bloch Hamiltonian in momentum ( ) space is given by k H 2 ( k ) = ( + t t ′ cos k σ ) x + ( t ′ sin k -i Γ 2) / σ z , which has been proposed in Ref. [47]. This model is also studied in Ref. [48] under a (pseudo-)spin rotation σ z → σ y , which corresponds to the unitary transformation UH 2 ( k ) U -1 = H ′ 2 ( k ), with U = ( iσ x + I 2 ) / √ 2 and I 2 as the 2 × 2 identity matrix. The Bloch Hamiltonian H 2 ( k ) has the chiral (sublattice) symmetry CH 2 ( k ) C -1 = -H 2 ( k ) and the parity-time ( PT ) symmetry PT H 2 ( k )( PT ) -1 = H 2 ( k ), where the chiral, parity and time reversal operators are respectively defined by C = σ y , P = σ x and T = K with the complex conjugation K .

These symmetries can also be introduced for the system under OBCs. The real-space Hamiltonian H 2 under

FIG. 10. (Color online) (a) Sketch of the disordered SSH model with non-Hermitian gain and loss, denoted by ± i Γ 2. / The dotted box is the unit cell, and m j , it ′ j / 2, and t / ′ 2 denote the corresponding hopping strengths. (b) The winding number ν as a function of t ′ and Γ for W 1 = W 2 = 0 and L = 24. The red dashed lines denote the PT symmetry breaking threshold. (c) ν as a function of W and Γ for t ′ = 1 1 , . t W 2 = W and W 1 = 0. (d,e) ν as a function of W and Γ for t ′ = , t W 1 = W (with W 2 = 0) and W 2 = W (with W 1 = 0), respectively. Other parameters are t = 1, L = 100, l = 0 2 . L and N s = 200.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01172v3-figure-908db496-13b9-43fa-b116-bceb85fecd25.jpg)

OBCs has the chiral symmetry CH C 2 -1 = -H 2 , where the chiral operator is C = σ y ⊗ I N with I N the N × N identity matrix for the system of N lattice cells. The parity-time ( PT ) symmetry for H 2 under OBCs is defined by PTH 2 ( PT ) -1 = H 2 , with the parity operator P = σ x ⊗ I N and the time-reversal operator T = K⊗ I N . In the presence of disorders, the chiral symmetry preserves and thus the (disorder-averaged) open-bulk winding number in real space can be calculated in this model (by rewriting the Q s matrix with two chiral-symmetric parts defined by C = σ y ⊗ I N here).

For the clean case of W 1 = W 2 = 0, the PT symmetry remains unbroken in the system under OBCs before | Γ | reaches the threshold | Γ = 2 . In this case, we find that | t the increase of the non-Hermitian gain-and-loss parameter | Γ can enlarge the parameter regime of the topologi-| cal phase with ν = 1, as shown in Fig. 10(b). In particular, the topological phase can be induced from an initially Hermitian trivial phase with t /t &lt; ′ 1 and Γ = 0 by turning on and increasing the gain-and-loss non-Hermiticity. The induced topological phase will return to the trivial phase for large | Γ | in the PT broken region. Note that the results preserve under the transformation Γ → -Γ. The enhancement of the topological phase by the gain and loss in the presence of the PT symmetry can also be understood by a similarity transformation as the energy spectrum is real under this symmetry. We apply a unitary transformation on the Hamiltonian H 2 under OBCs: UH U 2 -1 = H ′ 2 , where the transformation operator

$$U = \frac { 1 } { \sqrt { 2 } } \left ( \begin{array} { c c } 1 & i \\ i & 1 \end{array} \right ) \otimes \mathbb { I } _ { N }.$$

Then the real-space Hamiltonian H ′ 2 describes an open SSH chain with nonreciprocal intracell hopping terms t ± Γ 2 [48]. / When | Γ | &lt; 2 , t there exists a similarity transformation mapping the non-Hermitian Hamiltonian H ′ 2 onto a Hermitian counterpart ˜ H ′ 2 with intracell hopping strength ˜ = t √ t 2 -Γ 2 / 4 and intercell hopping strength ˜ = t ′ t ′ . Thus, the increase of | Γ | effectively increases the ratio ˜ t /t ′ ˜ up to | Γ | = 2 , which leads to the t enhancement of the topological phase with ˜ t /t &gt; ′ ˜ 1. Note that here t = 1 as the energy unit and the topological transition point is ˜ t /t ′ ˜ = 1.

glyph[negationslash]

In the presence of hopping disorders with W 2 = W and W 1 = 0 under the PT symmetry, the non-Hermitian gain-and-loss still enhances the topological phase in the open chain with t ′ &gt; t against disorders, as shown in Fig. 10(c). In addition, we find that the gain-and-loss can induce and enhance the topological phase from the critical case of t = t ′ , as shown in Fig. 10(d). When W 1 = 0, the PT symmetry can be spontaneously broken in certain disorder configurations with large W 1 , and then the non-Hermitian enhancement effect in the region of t ′ &gt; t may disappear (not shown here). However, the NHTAI phase can still be induced from the critical case in this SSH model by the combination of disorder ( W 1 = W and W 2 = 0) and gain-and-loss non-Hermiticity, as shown in Fig. 10(e). We notice that the proposed disordered SSH model with non-Hermitian gain and loss is realizable in some artificial systems, such as ultracold atoms, photonic crystals, and electronic circuits.
Dan-Wei Zhang Ling-Zhi Tang Li-Jun Lang Hui Yan Shi-Liang Zhu
Non-Hermitian Topological Anderson Insulators
Non-Hermitian systems can exhibit unique topological and localization
properties. Here we elucidate the non-Hermitian effects on disordered
topological systems by studying a non-Hermitian disordered Su-Schrieffer-Heeger
model with nonreciprocal hoppings. We show that the non-Hermiticity can enhance
the topological phase against disorders by increasing energy gaps. Moreover, we
uncover a topological phase which emerges only under both moderate
non-Hermiticity and disorders, and is characterized by localized insulating
bulk states with a disorder-averaged winding number and zero-energy edge modes.
Such topological phases induced by the combination of non-Hermiticity and
disorders are dubbed non-Hermitian topological Anderson insulators. We also
find that the system has non-monotonous localization behaviour and the
topological transition is accompanied by an Anderson transition. These
properties are general in other non-Hermitian models.