arxiv: v1 [cond-mat.mes-hall] 31 Mar 2019

Transcription

1 Emergent localize states at the interface of a twofol PT -symmetric lattice Jung-Wan Ryu, 1 Nojoon Myoung, Sungjong Woo, 1 Ara Go, 1 Sang-Jun Choi, 1, an Hee Chul Park 1, 1 Center for Theoretical Physics of Complex Systems, Institute for Basic Science (IBS), Daejeon 3416, South Korea Department of Physics Eucation, Chosun University, Gwangju 6145, Republic of Korea (Date: April, 019) arxiv: v1 [con-mat.mes-hall] 31 Mar 019 We consier the role of non-triviality resulting from a non-hermitian Hamiltonian that conserves twofol PT -symmetry assemble by interconnections between a PT -symmetric lattice an its time reversal partner. Twofol PT -symmetry in the lattice prouces aitional surface exceptional points that play the role of new critical points, along with the bulk exceptional point. We show that there are two istinct regimes possessing symmetry-protecte localize states, of which localization lengths are robust against external gain an loss. The states are emonstrate by numerical calculation of a quasi-1d laer lattice an a D bilayere square lattice. PACS numbers: Parity time (PT ) symmetric systems exhibit a phase transition through spontaneous symmetry breaking, from an unbroken PT-symmetric phase that keeps the eigenenergy real in non-hermitian systems, to a broken phase that contains complex conjugate energies [1]. As the general characteristic of non-hermitian systems, non- Hermiticity evelops from imaginary phase accumulation an imaginary potential ue to an imbalance of particle/energy flow where H H [ 5]. PT -symmetry is protecte in non-hermitian systems with a balance of energy gain an loss represente by the commutation relation [H,PT ] = 0, where H is a Hamiltonian an PT is a combination of parity an time-reversal symmetry operators. A wie range of PT -symmetric systems have been explore over several fiels, incluing optics [4 10], electronic circuits [11], atomic physics [1], an magnetic metamaterials [13]. Phase transitions in these systems occur via exceptional points (EPs), which are egenerate points of eigenenergies in non-hermitian systems that generate a Möbius strip structure of eigenenergies in parametric space because of the square-root branching property of the singular point [14, 15]. Such a topological structure has been reporte in microwave experiments [16, 17], optical microcavities [18], an a chaotic exciton-polariton billiar [19]. In non-hermitian systems, states localize at the eges, interfaces, an efects have recently attracte consierable attention not only in funamental stuies such as topology an symmetry but also in applications to quantum technologies. Non-Hermitian flat bans generate localize zero moes analogous to Majorana zero moes in conense matter physics, of which real energies are zero but imaginary energies are non-zero [0, 1]. These non-hermitian zero moes (NHZMs) are protecte by non-hermitian particle hole (NHPH) symmetry, which is also calle charge-conjugate symmetry, Electronic aress: aietheia@ibs.re.kr Electronic aress: hcpark@ibs.re.kr where ǫ i = ǫ j with Re(ǫ i) = 0 if i = j [, 3]. It has been shown that NHPH symmetry generates topological efect moes at the interface between non-hermitian lattices base on topologcially trivial Hermitian lattices [4, 5]. Otherwise, topologically protecte NHZMs have been propose in non-hermitian lattices base on topologically non-trivial Hermitian lattices, such as the Su Schrieffer Heeger moel [6 31]. In aition to NHZMs, non-hermitian boun states (NHBSs) are protecte by PT -symmetry, where ǫ i = ǫ j with Im(ǫ i) = 0 if i = j [3, 33]. Anomalous localize ege states in non- Hermitian lattice moels have been attribute to a wining number aroun an EP in momentum space [34 36]; localize states also exist at the interface between two lattices with the same topological orer but with istinct quantum phases, such as unbroken an broken PTsymmetric phases [37, 38]. A PT -symmetric system can be realize by combining the even-parity of the real potential an the oparity of the imaginary potential with respect to the PT -symmetric axis. Pairing the systems to be timereversal partners inuces twofol PT -symmetry, which generates two ifferent symmetric axes an guarantees two EPs. Such multiple EPs by multifol PT -symmetry as well as single EPs by simple PT -symmetry have been stuie both theoretically an experimentally [39 41]. In this work, we introuce a lattice containing twofol PT -symmetry to stuy how interplay between non-hermiticity an bulk symmetry affects the wavefunctions of the lattice. We fin two pairs of interface states that ecay exponentially in the space istinguishe by the bulk states. Particularly, we focus on robust localize states (NHZMs an NHBSs) of which localization length is unaffecte by variation of external gain/loss in two istinct regimes between the bulk EP an surface EPs. Let us prepare the PT -symmetric laer lattice with a basis of two sites (see Fig. 1 (a)). Non-Hermiticity is aopte by respective gain an loss at the two sites in the unit cell, a balance amount of which preserves PT -symmetry. A series of such unit cells forms a laer

2 FIG. 1: (color online). (a) Sketch of a PT -symmetric laer lattice with gain (re) an loss (blue) on the upper an lower lattices, respectively. (b) Real(black soli line) an imaginary (re ashe line) parts of the complex energy bans of a PT - symmetric laer lattice for γ < γ b (left), γ = γ b (mile), an γ > γ b (right). The parameters are = 3 an t = 1. (c) Sketch of the intertwine PT -symmetric laer lattice consiere here. FIG. : (color online). (a) Real an (b) imaginary parts of the complex eigenenergies of Fig. 1 (a) as a function of γ. γ b is the exceptional point relate to the PT -symmetric phase transition. (c) Real an () imaginary parts of the complex eigenenergies of Fig. 1 (c) as a function of γ. The gray shae region represents the bulk state eigenenergies, an the re an blue lines inicate the eigenenergies of the interface states. γ b is the bulk exceptional point an γ ± are the surface exceptional points for the interface states. The parameters are = 3 an t = 1, an the number of unit cells is N = 400. lattice with inter-cell hopping. The lattice Hamiltonian for this system reas Ĥ0 = l= ĥl, where ĥ l = c l hc l +c l T c l+1 +c l+1 Tc l (1) an c l ( l,a, l,b ). The matrix h iγσ z σ x escribes the non-hermicity with γ an the intra-cell hopping with, while T tσ 0 escribes the inter-cell hopping where σ x,z are the Pauli matrices an σ 0 is the ientity matrix. The matrix representation of the Hamiltonian for a single laer lattice in momentum space is written as H 0 (q) = h tcosqσ 0 through the Bloch wavefunction ψ(q) n einq c n, where q is the momentum vector. The (non-normalize) eigenstates are ψ η (q) = ( iγ/ η (γ/),) T, which are inepenent of q, where η = ± assigns the upper/lower bans. The corresponing eigenvalues of the Hamiltonian are E η (q) = tcosq + η (γ/), which are complex energy bans as E(k) = E r (k) + ie i (k) accoring to ifferent γ in Fig. 1 (b). One may notice that there exists a phase transition between unbroken an broken PT -symmetric phases through the EP. At this EP, two complex energy bans merge into a single energy ban, as shown in mile panel of Fig. 1 (b); otherwise, the real an imaginary bans are separate. The real an imaginary parts of the complex spectrum for a finite-size PT -symmetric laer lattice with 400 unit cells are numerically calculate as a function of γ in Fig. (a) an (b), respectively. Results show a collective PT -symmetric phase transition at the bulk EP γ b =. For γ < γ b, the energy spectrum pair is real in spite of non-hermiticity, with each part attracting each other up to the bulk EP unlike the Hermitian cases, as shown in Fig. (a). For γ > γ b, meanwhile, the energy spectrum pairs are complex conjugates, an repel each other on the imaginary energy plane with increasing γ from the bulk EP, as shown in Fig. (b). We note that the two bans are separable at any value of γ except the bulk EP, because there is no ban crossing for any q [4]. The twofol PT -symmetric laer lattice can be forme by merging the eges of the prepare lattice an its time-reversal partner to create an interface in a process calle intertwinement (see Fig. 1 (c)). The Hamiltonian operator of this intertwine lattice consists of Ĥ = ĤL +ĤR +Ĥc, () whereĥl = 1 l= ĥl is ahamiltonianforthe left semiinfinite non-hermitian lattice, an ĤR = T l=0ĥlt 1 is a Hamiltonian for the right semi-infinite lattice with complex conjugate operator T as a time-reversal operator. AcouplingHamiltonianĤc = c 0 Tc 1+h.c.connects the two time-reversal partners. Intertwinement introuces aitional PT -symmetry with respect to the z-axis an aitional EPs accoring to the symmetry. The energy spectrum of the finite-size twofol PT -symmetric laer lattice is treate by using a triiagonal Hamiltonian. As shown in Fig. (c) an (), the four emergent states (re/blue lines) not only show the PT -symmetric phase transition but also separation from the bulk states, while the bulk states are coincient with the states of the simple PT-symmetric lattice. The emergent states are exponentially localize at the interface, of which localization lengths efine by the ex-

3 of interface states, incluing NHBSs, by comparing the phase iagrams in which the white areas inicate forbien regimes of localize states. We note that, while invisibleinthe energyspectrain Fig. (c)an(), phase transitions appear for the localize states at the bulk EP, as shown in Fig. 3 (a). The interface states can be obtaine through the Schröinger equation H Φ = E Φ. Let us consier a linear combination of the simple laer lattice eigenstates to solve the intertwine lattice [see Supplemental Materials]. The ansatz of the wavefunction in the system is Φ = Φ L + Φ R = η=± {α η ψ(q η ) L +β η ψ( q η ) R }. The matching conition of the counter-propagating wavefunctions at the interface is ( γ ) cos(q+ ( γ ), sinq + sinq + +q ) = (3) 3 FIG. 3: (color online). (a) Localization lengths of o (blue line) an even (re line) interface states using the same parameters as Fig.. The NHZMs an NHBSs are robust against γ in the blue an re shae regimes II an III, respectively. The orange an green lines represent the analytic solutions of localization lengths of the o an even interface states, respectively. (b, c) Phase iagrams for the localization lengths of o an even states, respectively. Blue soli, black soli, an blue ashe lines in (b) inicate γ, γ b, an γ b, respectively. Black soli, re soli, an re ashe lines in (c) inicate γ b, γ +, an γ +b, respectively. No localize states exist in the white areas. The black otte lines inicate the localization lengths shown in (a). ponential ecay of the wavefunctions ψ e x/λ as a function of γ are shown in Fig. 3 (a). Two regimes having constant localization lengths are inicate by the shae regions in Fig.3 (a). In the blue-shae regime γ < γ < γ b, the localize interface states conserve NHPH symmetry as NHZMs, an in the re-shae regime γ b < γ < γ +, the localize interface states conserve PT -symmetry as NHBSs. Otherwise, the localize states are regare as efect states in a ban gap. The localization lengths of the NHZMs an NHBSs are inepenent of the non-hermiticity control parameter γ, which is relate to the interplay between the two ifferent PT -symmetric phases of interface an bulk states on the lattice protecting twofol PT -symmetry. The NHZMs appear when the bulk an interface states have unbroken an broken phases, respectively, while conversely the NHBSs appear when the bulk an interface states have broken an unbroken phases, respectively. Otherwise, the interface an bulk states have the same unbroken or broken phases. The phase iagrams in Fig. 3 (b) an (c) show the localization lengths of the interface states of NHZMs an NHBSs as a function of intra-cell hopping an non-hermiticity parameter γ, respectively. When < t, we can recognize that there exist only a pair where cosq η = E+η (γ/) is the momentum vector corresponing to each energy ban, with all energy scales being imensionless by t from here. This conition provies exact solutions for the interface states as ( E η ± = η 1± )( ( ) γ ±, (4) ) which are eigenenergies of the two paire states separate from the bulk bans, where the inex ± inicates the even/o interface states as the re/blue curves in Fig. (c) an (), respectively. It is easy to see that the two surface EPs are γ ± = ± an that the interface states unergo phase transitions from unbroken to broken PT -symmetric phases. There is no interface state in the Hermitian limit, i.e., γ = 0, an the interface states with real eigenenergies separate from the bulk states as γ increases. The even an o interface states finally merge into the bulk states at γ +b = + an γ b =, respectively. To analyze the characteristics of the emergent interface states, we introuce complex momentum vectors of the states as q ± η = k ± η + iκ± η, where k± η an κ ± η are real an ± inicates the even/o state. The ispersion relation is generalize by complex momentum cos(k +iκ) = coskcoshκ isinksinhκ. By substituting the energies of the interface states as foun in Eq. (4) for the energy in the ispersion relation, cosq η ± = E± η +η (γ/), (5) we can formulate localization lengths as λ ± = 1/κ ± η at the interfaces of the intertwine lattice with respect to the istinct regimes of control parameter γ [see Supplemental Materials]. For the NHZMs, we can evaluate the inverse localization length an wave number for oscillation as follows, coshκ η =, (γ/) cosk η = η. (6)

4 4 FIG. 4: (color online). (a) Imaginary parts of complex eigenenergies as a function of γ in the D laer lattice. The upper inset plots the eigenenergies on the complex plane for γ = 7.6, an the lower inset illustrates the D lattice with a 1D straight interface along the y-axis. (b) Intensity istribution of a selecte lossless NHBS for γ = 7.6. (c) Top view of the D lattice schematic with an arbitrarily shape close interface, with a cut-view showing the interface. The upper an lower layers have opposite profiles of gain an loss. () Intensity istribution of a selecte lossless NHBS for γ = 6.5. The parameters in this figure are = 3 an t x = t y = 1, an the number of unit cells is In case of the NHBSs, the inverse localization length an wave number of the localize states are coshκ + η = +1, + (γ/) cosk+ η = η. (7) The localization lengths 1/κ ± η from Eq. (6) an Eq. (7) agree well with numerical results, as shown in Fig. 3 (a). It is meaningful to unerstan the surface EPs through the eigenenergies of the twofol PT -symmetric four sites containing self-energy, which have the same symmetries as our moel. The effective Hamiltonian reas H eff = τ 0 ( σ x +i γ σ z)+τ x ( σ x ), (8) where τ i is the Pauli matrix with ientity matrix τ 0 for inter-cell hopping. The Hamiltonian is block iagonalize by means of unitary matrix L = (τ x +τ z ) σ 0 /. One of the blocks is H + = iγ /σ z ( + 1)σ x with eigenvalues E η + = η (+1) (γ /), an the other block is H = iγ /σ z ( 1)σ x with eigenvalues Eη = η ( 1) (γ /). We can see two ifferent critical values of γ ± with respect to the separate blocks. Suppose that the on-site energy is renormalize by γ = γ +4 in terms of the self-energy of the semiinfinite leas; then the block iagonal matrix has two EPs, γ ± = ±, which are the same as the surface EPs in the intertwine lattice. It shoul be note that three symmetry-protecte EPs exist in the twofol PTsymmetric laer lattice: two surface EPs arising from the twofol PT -symmetry, an a bulk EP from the PTsymmetric lattice itself. We now propose a two-imensional (D) lattice in which a issipationless one-imensional (1D) waveguie is realizable by repeating 1D twofol PT -symmetric laer lattices. For example, a straight waveguie is consiere on the Dlattice as the interface along the y-axis, as shownin the lowerinset offig. 4 (a). This Dlattice satisfies both translational symmetry along the y-axis an twofol PT -symmetry of each 1D laer lattice. Figure 4 (a) plots the imaginary parts of the eigenenergies as a function of γ in the D bilayere square lattice with unit cells, where it is apparent that the imaginary eigenenergies are broaene when compare to those of the 1D lattice in Fig. () because of finite-size effects. As shown in the upper inset of Fig. 4 (a), the istribution of the complex eigenenergies is similar to that of the 1D lattice, except for the real energy istribution ue to the ispersion towar the y-axis. The zero-imensional localize NHBSs, which are robust across a wie range of non-hermiticity in the 1D lattice, exten into the 1D NHBSs in the D lattice. The corresponing eigenstates, which exten along the y irection with no issipation along the x-axis, are epicte in Fig. 4 (b), with wavelengths relate to the real parts of the eigenenergies. Finally, we show the robustness of the issipationless 1D NHBSs against geometrical eformations. We esign a D PT -symmetric layer compose of square lattices that contain a locally inverte area with an arbitrary shape, as shown in Fig. 4 (c). Notably, the geometry presents many abrupt corners where the interface bens but still conserves symmetry. We can see that the NHBSs survive without issipation at the interfaces with sharp corners, as shown in Fig. 4 (), while there exists bening loss aroun the sharply eforme region of a trivial waveguie. Therefore, NHBSs can be implemente, without bening loss, as the effective waveguie in a D lattice while conserving local twofol PT -symmetry. In conclusion, we propose a system in which localize states emerge solely through the non-triviality resulting from non-hermiticity. Our twofol PT -symmetric laer lattice contains symmetry-protecte interface states as NHZMs an NHBSs with corresponing phase transitions at the two surface EPs an bulk EP as a function of γ. The characteristics of these two states are as follows. The NHZMs protecte by NHPH symmetry are localize at the interface between the o surface EP an the bulk EP, while the NHBSs protecte by PT -symmetry are localize at the interface between the even surface EP an the bulk EP. Both have constant localization lengths unaffecte by changes in the non-hermiticity parameter. As an example, we showehere that the NHBSs can form an effective waveguie without issipation at the interface between non-hermitian time-reversal partners. We expect the characteristics of these two symmetry-protecte

5 5 interface states to open up a new fiel of synthetic non- Hermitian systems. grant (NRF016-R1D1A1B ) an Korea Institute for Avance Stuy(KIAS) fune by the Korea government (MSIT). Acknowlegments This work was supporte by Project coe (IBS-R04- D1), the National Research Founation of Korea (NRF) [1] C. M. Bener an S. Boettcher, Real Spectra in Non- Hermitian Hamiltonians Having PT Symmetry, Phys. Rev. Lett. 80, 543 (1998). [] N. Hatano, an D. R. Nelson, Localization Transitions in Non-Hermitian Quantum Mechanics, Phys. Rev. Lett. 77, (1996). [3] I. Klich, Close hierarchies an non-equilibrium steay states of riven systems, Ann. of Phys. 404, (019). [4] R. El-Ganainy, K. G. Makris, D. N. Christooulies, an Z. H. Musslimani, Theory of couple optical PTsymmetric structures, Opt. Lett. 3, 63 (007). [5] C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christooulies, M. Segev, an D. Kip, Observation of parity-time symmetry in optics, Nat. Phys. 6, (010). [6] Z. H. Musslimani, K. G. Makris, R. El-Ganainy, an D. N. Christooulies, Optical Solitons in PT Perioic Potentials, Phys. Rev. Lett. 100, (008). [7] K. G. Makris, R. El-Ganainy, D. N. Christooulies, an Z. H. Musslimani, Beam Dynamics in PT Symmetric Optical Lattices, Phys. Rev. Lett. 100, (008). [8] S. Klaiman, U. Günther, an N. Moiseyev, Visualization of Branch Points in PT-Symmetric Waveguies, Phys. Rev. Lett. 101, (008). [9] A. Guo, G. J. Salamo, D. Duchesne, R. Moranotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, an D. N. Christooulies, Observation of P T -Symmetry Breaking in Complex Optical Potentials, Phys. Rev. Lett. 103, (009). [10] A. Regensburger, C. Bersch, M. A. Miri, G. Onishchukov, D. N. Christooulies, an U. Peschel, Parity-time synthetic photonic lattices, Nature (Lonon) 488, (01). [11] J. Schinler, A. Li, M. C. Zheng, F. M. Ellis, an T. Kottos, Experimental stuy of active LRC circuits with PT symmetries, Phys. Rev. A 84, (011). [1] Y. N. Joglekar, D. Scott, M. Babbey, an A. Saxena, Robust an fragile PT-symmetric phases in a tight-bining chain, Phys. Rev. A 8, (R) (010). [13] N. Lazaries an G. P. Tsironis, Phys. Gain-Driven Discrete Breathers in PT-Symmetric Nonlinear Metamaterials, Rev. Lett. 110, (013). [14] T. Kato, Perturbation Theory of Linear Operators (Springer, Berlin, 1996). [15] W. D. Heiss, The physics of exceptional points. J. Phys. A 45, (01). [16] C. Dembowski, H.-D. Gräf, H. L. Harney, A. Heine, W. D. Heiss, H. Rehfel, an A. Richter, Experimental Observation of the Topological Structure of Exceptional Points, Phys. Rev. Lett. 86, 787 (001). [17] C. Dembowski, B. Dietz, H.-D. Gräf, H. L. Harney, A. Heine, W. D. Heiss, an A. Richter, Observation of a Chiral State in a Microwave Cavity, Phys. Rev. Lett. 90, (003). [18] S.-B. Lee, J. Yang, S. Moon, S.-Y. Lee, J.-B. Shim, S. W. Kim, J.-H. Lee, an K. An, Observation of an Exceptional Point in a Chaotic Optical Microcavity, Phys. Rev. Lett. 103, (009). [19] T. Gao, E. Estrecho, K. Y. Bliokh, T. C. H. Liew, M. D. Fraser, S. Brobeck, M. Kamp, C. Schneier, S. Höfling, Y. Yamamoto, F. Nori, Y. S. Kivshar, A. G. Truscott, R. G. Dall, an E. A. Ostrovskaya, Observation of non- Hermitian egeneracies in a chaotic exciton-polariton billiar, Nature 56, (015). [0] L. Ge, Parity-time symmetry in a flat-ban system, Phys. Rev. A 9, (015). [1] L. Ge, Non-Hermitian lattices with a flat ban an polynomial power increase, Photon. Res. 6, A10 (018). [] L. Ge, Symmetry-protecte zero-moe laser with a tunable spatial profile, Phys. Rev. A 95, 0381 (017). [3] B. Qi, L. Zhang, an L. Ge, Defect States Emerging from a Non-Hermitian Flatban of Photonic Zero Moes, Phys. Rev. Lett. 10, (018). [4] S. Malzar, C. Poli, an H. Schomerus, Topologically Protecte Defect States in Open Photonic Systems with Non-Hermitian Charge-Conjugation an Parity- Time Symmetry, Phys. Rev. Lett. 115, 0040 (015). [5] S. Malzar an H. Schomerus, Bulk an ege-state arcs in non-hermitian couple-resonator arrays, Phys. Rev. A 98, (018). [6] H. Schomerus, Topologically protecte migap states in complex photonic lattices, Opt. Lett. 38, (013). [7] C. Poli, M. Bellec, U. Kuhl, F. Mortessagne, an H. Schomerus, Selective enhancement of topologically inuce interface states in a ielectric resonator chain, Nat. Commun. 6, 6710 (015). [8] J. M. Zeuner, M. C. Rechtsman, Y. Plotnik, Y. Lumer, S. Nolte, M. S. Runer, M. Segev, an A. Szameit, Observation of a Topological Transition in the Bulk of a Non- Hermitian System, Phys. Rev. Lett. 115, (015). [9] H. Zhao, P. Miao, M. H. Teimourpour, S. Malzar, R. El- Ganainy, H. Schomerus, an L. Feng, Topological hybri silicon microlasers, Nat. Commun. 9, 981 (018). [30] M. Parto, S. Wittek, H. Hoaei, G. Harari, M. A. Banres, J. Ren, M. C. Rechtsman, M. Segev, D. N. Christooulies, an M. Khajavikhan, Ege-Moe Lasing in 1D Topological Active Arrays, Phys. Rev. Lett. 10, (018). [31] W. P. Su, J. R. Schrieffer, an A. J. Heeger, Solitons in polyacetylene, Phys. Rev. Lett. 4, 1698 (1979). [3] S. Weimann, M. Kremer, Y. Plotnik, Y. Lumer, S. Nolte,

6 K. G. Makris, M. Segev, M. C. Rechtsman, an A. Szameit, Topologically protecte boun states in photonic parity-time-symmetric crystals, Nat. Mater. 16, (017). [33] X. Ni, D. Smirnova, A. Poubny, D. Leykam, Y. Chong, A. B. Khanikaev, Exceptional points in topological ege spectrum of PT symmetric omain walls, Phys. Rev. B 98, (018). [34] T. E. Lee, Anomalous Ege State in a Non-Hermitian Lattice, Phys. Rev. Lett. 116, (016). [35] V. M. Martinez Alvarez, J. E. Barrios Vargas, an L. E. F. Foa Torres, Non-Hermitian robust ege states in one imension: Anomalous localization an eigenspace conensation at exceptional points, Phys. Rev. B 97, 11401(R) (018). [36] D. Leykam, K. Y. Bliokh, C. Huang, Y. D. Chong, an F. Nori, Ege Moes, Degeneracies, an Topological Numbers in Non-Hermitian Systems, Phys. Rev. Lett. 118, (017). [37] H. Zhao, S. Longhi, an L. Feng, Robust light state by quantum phase transition in non-hermitian optical materials, Sci. Rep. 5, 170 (015). [38] M. Pan, H, Zhao, P. Miao, S. Longhi, an L. Feng, Photonic zero moe in a non-hermitian photonic lattice, Nat. Commun. 9, 1308 (018). [39] L. Ge an A. D. Stone, Parity-Time Symmetry Breaking beyon One Dimension: The Role of Degeneracy, Phys. Rev. X 4, , (014). [40] Z. Liu, Q. Zhang, X. Liu, Y. Yao, an J.-J. Xiao, Absence of Exceptional Points in Square Waveguie Arrays with Apparently Balance Gain an Loss, Sci. Rep. 6, 711 (016). [41] K. Ding, G. Ma, M. Xiao, Z. Q. Zhang, an C. T. Chan, Emergence, Coalescence, an Topological Properties of Multiple Exceptional Points an Their Experimental Realization, Phys. Rev. X, 6, (016). [4] H. Shen, B. Zhen, an L. Fu, Topological Ban Theory for Non-Hermtian Hamiltonian, Phys. Rev. Lett. 10, (018). 6

7 7 Supplemental Materials: Emergent localize states through intertwine non-hermitian lattice I. HAMILTONIAN OF A LADDER LATTICE The Hamiltonian of the laer lattice is given by EΨ j = H 0 Ψ j +H 1 Ψ j+1 +H + 1 Ψ j 1, (9) where H 0 = ( ) ( ) ǫa t 0,H ǫ 1 =, (10) b 0 t an Ψ j = (φ a j,φb j )T. We can set Ψ j+1 = Ψ j e ik an Ψ j 1 = Ψ j e ik ue to the translational symmetry of the unit cells. Finally, ( ) ǫa tcosk H =. (11) ǫ b tcosk Solving the eigenproblem of H when ǫ a = ǫ b = iγ/, we obtain the ban structure for the PT-symmetric laer lattice shown in the inset of Fig. (a). Figure 5 shows the real an imaginary parts of the complex energy bans with = 3 an t = 1 when γ =.0, 6.0, an II. NUMERICAL METHOD FOR EIGENENERGIES AND EIGENSTATES IN THE TWOFOLD PT-SYMMETRIC LADDER LATTICE We numerically obtain the eigenenergies an eigenstates in a finite-size twofol PT-symmetric laer lattice with N unit cells given by.... H 0 H 1 H + 1 H 0 H 1 H 1 + H 0 H, (1) 1 H 1 + H 0... where H 0 = ( ) ǫa ǫ b ( ) H 0 = ǫb ǫ a for a laer lattice, an (13) for a twiste laer lattice. (14) Solving this N N matrix, we can obtain N eigenenergies. For instance, in the Hermitian case of ǫ a = ǫ b = δ/ (δ is real), the eigenenergies as functions of δ are shown in Fig. 6. III. ANALYTIC SOLUTION OF THE INTERFACE STATES The left an right lattices of Fig. 7 can be escribe by H ψ L = E ± (k) ψ L te ipk (a p+1,a +b p+1,b )+te i(p+1)k (a p,a +b p,b ), (15) H ψ R = E ± (k) ψ R te iqk (b q 1,A +a q 1,B )+te i(q 1)k (b q,a +a q,b ), (16)

8 8 FIG. 5: (color online). (a c) Real an ( f) imaginary parts of the complex energy bans of PT-symmetric laer lattices with (a, ) γ =.0, (b, e) γ = 6.0, an (c, f) γ = where H = Σ n [iγ( n,a n,a n,b n,b ) (17) ( n,a n,b + n,b n,a ) t( n+1,a n,a n+1,b n,b +h.c.)], (18) E ± (k) = tcosk ± (γ/), (19) ψ L Σ p n= e ink (a n,a +b n,b ), (0) ψ R Σ n=q e ink (b n,a +a n,b ), (1) a ± = iγ/ (γ/), () b ± =. (3) Here, cosk ± = E± +(γ/) t an sink ± = 1 cos k ± for L, while k ± k ± for R. The ansatz is Φ = Φ L + Φ R, where Φ L = α + ψ(k + ) L + α ψ(k ) L an Φ R = β + ψ( k + ) R + β ψ( k ) R. We fin the conition satisfying H Φ = E Φ as ( γ ) cos(k+ ( γ ). sin(k )sin(k + )+ +k ) = (4) This is erive by ( α + a + e ik+ α a e ik +β + b + e ik+ +β b e ik ) q,a +(α + a + +α a β + b + β b ) p,a (5) +( α + b + e ik+ α b e ik +β + a + e ik+ +β a e ik ) q,b +(α + b + +α b β + a + β a ) p,b = 0. Then, the square matrix below shoul vanish for non-trivial solutions: a + e ik+ a e ik b + e ik+ b e ik α + 0 a + a b + b α 0 b + e ik+ b e ik a + e ik+ a e ik β = + 0. (6) b + b a + a β 0

9 9 FIG. 6: (color online). (a) Real eigenenergies as a function of δ when ǫ a = ǫ b = δ, = 3, t = 1, an N = 100 in a laer lattice of the inset in Fig. (a). (b) Real eigenenergies as a function of δ when ǫ a = ǫ b = δ, = 3, t = 1, an N = 100 in a twiste laer lattice of the inset in Fig. (c). FIG. 7: (color online). An intertwine PT -symmetric laer lattice we consiere. By putting cosk ± = E± +(γ/) t an sink ± = 1 cos k ± into Eq. (4), then the energy (E) of the interface states can be evaluate as a function of γ by ( E± I = ± 1 t ( E± O = ± 1+ t )( ( ) γ t, (7) ) )( ( ) γ +t, (8) ) where E I ± an E O ± are the eigenenergies of the eigenstates with o an even parities about the interface, respectively. Putting Eq. (7) into cosk ± = E± (γ/) t, complex wave numbers k± I of the interface states are given by cosk I ± = (1 t ) ( t ( ) ) γ t ± ( γ ). (9) This k± I can be classifie into three phases by the two transition points γ I = t an γ b =. Setting k± I = k ± I + iκi ±, the imaginary parts κi ± of complex wave number ki ±, which are the reciprocals of the localization

10 lengths of the interface states, are while the real parts k I ± are κ I ± = (1 t cosh 1 ) ( t ( γ ) ) + ( γ ) t for γ < γ I, (30) ( ) κ I ± = cosh 1 for γ I < γ < γ b, (31) t κ I ± = (1 t ) ( t ( γ ) ) + ( γ sinh 1 ) ti for γ b < γ, (3) k ± I = 0 or ±π for γ < γ I, (33) ) k ± I = cos ( 1 ( γ ) for γ I < γ < γ b, (34) t k I ± = π for γ b < γ. (35) Putting Eq. (8) into cosk ± = E± (γ/) t, complex wave numbers k± O of the interface states are given by (1+ ) ( t cosk± O = +t ( ) ) γ ± ( γ ). (36) t Likewise, k± O can be classifie into three phases by the two transition points γ b = an γ O = +t. Setting k± O = k ± O + iκ O ±, the imaginary parts κ O ± of complex wave number k±, O which are the reciprocals of the localization lengths of the interface states, are while the real parts k O ± are κ O ± = cosh 1 (1+ t ) ( +t ( γ ) ) ( γ ) t for γ < γ b, (37) ( ) κ O ± = cosh 1 t +1 for γ b < γ < γ O, (38) κ O ± = sinh 1 (1+ t ) ( +t ( γ ) ) + ( γ ) ti for γ O < γ, (39) k ± O = 0 or ±π for γ < γ b, (40) ) k ± (± O = cos 1 +t ( γ ) for γ b < γ < γ O, (41) t k O ± = ±π for γ O < γ. (4) 10 The Hamiltonian of the D laer lattice is given by where IV. HAMILTONIAN OF A D LADDER LATTICE EΨ j,k = H 0 Ψ j,k +H 1 Ψ j+1,k +H + 1 Ψ j 1,k +H Ψ j,k+1 +H + Ψ j,k 1, (43) H 0 = ( ) ( ) ( ) ǫa tx 0 ty 0,H ǫ 1 =,H b 0 t =, (44) x 0 t y an Ψ j,k = (φ a j,k,φb j,k )T. We can set Ψ j+1,k = Ψ j,k e ikx, Ψ j 1,k = Ψ j,k e ikx, Ψ j+1,k = Ψ j,k+1 e iky, an Ψ j,k 1 = Ψ j,k e iky ue to the translational symmetry of the unit cells. Finally, ( ) ǫa t H = x cosk x t y cosk y. (45) ǫ b t x cosk x t y cosk y