Topological non-Hermitian skin effect

This article reviews recent developments in the non-Hermitian skin effect (NHSE), particularly on its rich interplay with topology. The review starts off with a pedagogical introduction on the modified bulk-boundary correspondence, the synergy and hybridization of NHSE and band topology in higher dimensions, as well as, the associated topology on the complex energy plane such as spectral winding topology and spectral graph topology. Following which, emerging topics are introduced such as non-Hermitian criticality, dynamical NHSE phenomena, and the manifestation of NHSE beyond the traditional linear non-interacting crystal lattices, particularly its interplay with quantum many-body interactions. Finally, we survey the recent demonstrations and experimental proposals of NHSE.

The past two decades have witnessed a burgeoning interest in the intriguing properties of non-Hermitian Hamiltonians [1,2].Compared with Hermitian Hamiltonians with real eigenenergies, which represent unitary time evolution of isolated quantum systems, non-Hermitian Hamiltonians provide an effective physical description of many non-conservative systems, including open quantum systems [3], solid-state systems with finite lifetime induced by interactions [4][5][6], and acoustic/photonic systems with gain and loss [7][8][9][10][11][12][13]. Most early investigations of non-Hermitian Hamiltonians have been devoted to systems with parity-time (PT) symmetry [1, 10,11,14].As a specific form of pseudo-Hermiticity [15], PT symmetry allows a non-Hermitian Hamiltonian with balanced gain and loss to have a real spectrum, enabling a stable unitary time-evolution for eigenstates of the system.On the other hand, a PT symmetric Hamiltonian does not guarantee that the symmetry is also possessed by each of the eigenstates.In the so-called PT-broken phase, a PT symmetric Hamiltonian can have pairs of complex-conjugate eigenenergies, whose eigenstates can be transformed into each other through a PT symmetry operation.The transition between PT-symmetric and PT-broken phases is accompanied by the emergence of exceptional points (EP) for the system, a type of spectral degeneracies of non-Hermitian systems where two or more eigenstates coalesce into one and the Hamiltonian matrix becomes rankdeficient [3,[16][17][18][19][20][21].Besides providing a signal of PT transition, EPs are physically interesting in their own good, as they can lead to various exotic features of non-Hermitian systems, such as unidirectional invisibility [22,23] and enhanced sensitivity [24][25][26][27] in photonics, and unusual topological [28][29][30][31][32][33][34][35] and dynamical properties [36,37] for non-Hermitian Hamiltonians encircling an EP.
Among the many fascinating aspects of non-Hermitian Hamiltonians, one noteworthy route that triggers an enormous number of extensive studies over the years is how non-Hermiticity affects the topological properties of lattice systems.In the study of topological phases of matter, a most fundamental principle is the bulk-boundary correspondence (BBC), i.e. the number of edge states under the open boundary conditions (OBCs) has an one-on-one correspondence to a topological invariant defined for the bulk states under the periodic boundary conditions (PBCs) [38][39][40].In early studies, it was commonly believed that the BBC still holds for non-Hermitian systems, despite that non-Hermitian topological phases may behave differently from their Hermitian counterparts [41][42][43][44][45][46][47].However, in 2016, T. E. Lee noticed that non-Hermiticity can break conventional BBC for topological phases, evidenced by an inconsistency between the spectra of a system under the periodic and open boundary conditions (PBCs and OBCs) [31].This phenomenon was firstly understood as a consequence of a half-integer winding number for non-Hermitian Hamiltonians encircling an EP in momentum space [31][32][33], but this winding number, whose definition involves only PBC Hamiltonians, does not predict the number of topological edge states under the OBCs accurately when the system is closed to a topological phase transition.Another interpretation is to consider the evolution of the system from PBCs to OBCs, which changes the topological structure of the spectrum by passing through a series of EPs [48], but it does not provides a solution to restore the BBC either.This problem is eventually solved in 2018 by the sensational formulation of the non-Hermitian skin effect (NHSE), under which all eigenstates are spatially localized at the boundary of a system under the OBCs [49,50].Independently discovered by two research groups [49,50], the NHSE has motivated the socalled non-Bloch band theory, recovering a generalized BBC for non-Hermitian Hamiltonians [50][51][52][53].A parallel formulation of the non-Bloch Hamiltonian is to consider the PBC-OBC spectral evolution with a tunable hopping strength across the boundary, effectively giving rise to an imaginary flux to the system [48,54].Ref. [54], which pre-dated the detailed study of the generalized Brillouin zone (GBZ) of 1D models with arbitrarily far hoppings, also developed a criterion for the presence of topological zero modes in generic 1D 2band PH-symmetric models with arbitrarily far hoppings, requiring only information on the PBC dispersion poles/zeros, and not the GBZ.It is also worth noticing that the eigenstates of non-Hermitian Hamiltonians are not necessarily mutually orthogonal, and the biorthogonal condition needs to be taken into account to obtain a proper topological invariant from bulk states [55].A recent work demonstrates a numerically efficient relationship between EPs and biorthogonality [56].
In this review, after introducing the generalized BBC of the non-Bloch band theory for systems with NHSE, we will further discuss how NHSE affects localization of topological edge states.That is, beyond modifying the BBC of topological phases, NHSE can also change the direction of localization of topological edge states as it provides another localization mechanism for eigenstates [94][95][96][97].This phenomenon is also associated with the single (defective) edge state reported in Ref. [31], which can be viewed as two exceptionally degenerate edge states in the thermodynamic limit [33,98].A more intriguing development of the interplay between NHSE and conventional topological properties is their simultaneous effect of NHSE in 2D or higher dimensions, which leads to a new class of hybrid skin-topological higher-order boundary modes, living in boundaries with a co-dimension of two or higher [99][100][101][102].It is found later that this hybrid skintopological effect can be regarded as a special kind of a new class of higher-order NHSE, where conventional topological protection is not a necessary condition for these higher-order boundary modes [103][104][105].Besides affecting conventional topological properties, NHSE itself has also been found to be a signature of a spectral winding topology unique in non-Hermitian systems with complex spectrum [106][107][108], which is the second main theme of this review.A quantized response of the spectral winding is later found in the steady-state response of such systems, without involving the linear response theory based on many-body ground states [109,110].On the other hand, nontrivial spectral winding is also found during the investigation of non-Hermitian quasicrystal, the spectral winding topology of which is found to not necessarily correspond to NHSE in real space, but in reciprocal space for systems without translational symmetry [111][112][113][114][115][116][117][118][119][120][121][122].In [123], interesting mathematical relations are further found between the NHSE, knot theory and spectral winding topology in the context of 3D exceptional metals.In the later half of this article, we also review how the NHSE can non-perturbatively modify state dynamics as well as signatures of criticality, in both single-body and interacting many-body systems.We also note that while this review will mainly focus on theoretical aspects, we will devote the last section to reviewing experimental demonstrations and proposals of NHSE-related phenomena.

II. MODIFIED BAND TOPOLOGY FROM THE NHSE
Topological BBC plays a central role in describing topological phases of matter, as it associates topological boundary modes with bulk topological invariants, analogous to the order parameters for conventional quantum phases and phase transitions.Without a doubt, the breakdown of BBC in non-Hermitian systems has attracted great attention since its discovery, and many efforts have been made in recovering it with different methods.For examples, in 1D non-Hermitian systems, BBC have been restored between topological winding numbers and different properties of the systems, including zero-energy edge states in an semi-infinite system [33], a singular spectrum obtained from a singular-value decomposition [124], and spatial growth of the bulk Green function [125], to mention a few.In this section, we will introduce the method of non-Bloch band theory and NHSE, which recovers a modified BBC between boundary states under OBC and topological invariants defined in a so-called generalized Brillouin zone (GBZ) -a complex continuation of the conventional Brillouin Zone.
A. NHSE and non-Bloch band topology: minimal non-Hermitian SSH model As a starting point, we shall first give a brief review of NHSE in a non-Hermitian Su-Schrieffer-Heeger (SSH) model [126] [see Fig. 1(a)], which is a minimal model for demonstrating the generalized bulk-boundary correspondence and non-Bloch band theory for non-Hermitian band topology, by reproducing the results in Ref. [50].A comprehensive study of non-Hermitian extensions to another representative 1D topological model, namely the Creutz model [127], can be found in Ref. [128].The real-space and Bloch Hamiltonians of the non-Hermitian SSH model are respectively given by with α † n the creation operator acting on α sublattice of the nth unit cell, and where h x = t 1 + t 2 cos k, h y = t 2 sin k,and σ x,y are the Pauli matrices.Non-Hermiticity is introduced through nonzero γ, which is responsible for non-reciprocity of this system.Under the OBC, the bulk eigenstates of this non-Hermitian SSH model are localized near the left boundary, known as the "non-Hermitian skin effect" (NHSE).The difference of this model between OBC and PBC can be understood with a similarity transformation of the Hamiltonian, which maps the model under OBC into a Hermitian one while keeping its eigenvalues unchanged.To see this, we first write down the eigenequation as with real-space eigenvectors taking the form of |ψ⟩ = Under OBC, we can apply a similarity transformation with and L the total number of unit cells of the 1D chain, so that HOBC becomes the Hermitian SSH model with intracell and intercell hoppings given by t1 = (t 1 − γ/2)(t 1 + γ/2), t2 = t 2 respectively.Therefore, HOBC In the non-Hermitian SSH model, a generalized BBC is restored between Hermitian Hamiltonians H OBC and H PBC through the similarity transformation.However, in general non-Hermitian Hamiltonians have complex OBC spectrum and cannot be transformed into Hermitian ones with a similarity transformation, thus we need a more universal approach to establish a generalized BBC of non-Hermitian Hamiltonians.We first explicitly write down the bulk eigenequation Eq.(3) under OBC, where E is eigenenergy and the index of unit cell n = 1, 2, • • • , L − 1, with boundary conditions Due to the translation symmetry of the bulk equation, we can take an ansatz of the eigenvectors, |ψ⟩ = j |ϕ ( j) ⟩ with (ϕ ( j) n,A , ϕ ( j) n,B ) = z n (ϕ ( j) A , ϕ ( j) B ), with j indexes different solutions of ϕ.Therefore Eq. ( 5) becomes or On the other hand, Eq. ( 7) also leads to which has two solutions z 1,2 satisfying Since z 1,2 correspond to the same eigenenergy E, a linear combination of (ϕ ( j) A ϕ ( j) B ) T with ( j = 1, 2) is also an eigenvector of the Hamiltonian, given by Combining with the boundary conditions of Eq. ( 6), we shall arrive at The solution of this equation is nontrival only when in the thermodynamic limit (L → ∞).It is equivalent to z 1 = 0 or t 1 − γ/2 + t 2 z 2 = 0, and in this situation there is only one pair of (z 1 , z 2 ) in Eq. (11), which is independent of the system size L and gives only a single solution of eigenenergy E. But obviously the number of the solutions of Eq. ( 12) should be proportional to the system size L, otherwise the system under OBC would not have continuum bands.Therefore we obtain the condition of continuum bands, Eq. ( 14) is also the kernel of non-Bloch band theory [51].That is, since it only restricts the absolute values of z 1,2 , the continuum bands can be given by Eq. ( 9) with with k varying from 0 to 2π.In other words, the phase factor e ik in usual Bloch waves of Hermitian systems is now replaced by z GBZ .Here GBZ stands for the "generalized Brillouin zone", a concept analogous to the Brillouin zone (BZ) with z BZ = e ik , but describing OBC spectrum of non-Hermitian systems (see Fig. 2 for a comparison of BZ and GBZ).In the literature, z GBZ is sometimes denoted as β.It can also be viewed as a complex deformation of the momentum, k → k + iκ, with κ = ln r describing the inverse localization length of skin modes [51,54,106].For this model, κ is constant, but generically κ would be a function of k, which we will discuss in the next subsection.
Recently, the topological phases of the non-Hermitian SSH model and the non-Bloch band theory has been studied in terms of the electronic polarization [129][130][131].Beyond OBC systems, the GBZ has also be used to describe systems with other boundary conditions such as a domain wall [132], onsite disorder [133] or a strong local impurity [65].

B. Non-Bloch band theory for general cases
For more general cases, the generalized Brillouin zone may not be a unit circle in the complex plane, i.e. the imaginary flux κ may be k-dependent, meaning that different PBC eigenstates correspond to OBC skin states with different decay rates.In such cases, it is impossible to find a similarity transformation to remove the NHSE of the whole system, as each k-eigenstate corresponds to a different localization length.Nevertheless, the GBZ z GBZ can still be determined through an analyzing the characteristic polynomial of the Hamiltonian.Here we shall briefly summarize the method for a general 1D single-band non-Hermitian system with nonreciprocal hoppings, described by the Hamiltonian with ĉx the annihilation operator acting on site x, and N the total number of lattice sites.By solving the eigenvalue equation HΨ s (E) = EΨ s (E) with an ansatz for an eigenstate with eigenenergy E as we obtain a characteristic polynomial equation for H with a given E, where H(z s ) is the same as the Hamiltonian in momentum space for z s = e ik , and there shall be m + n solutions of z s marked by its s index.Note that so far we only assume x starts from 1, but have not considered boundary conditions of the system.As a matter of fact, these solutions of z s correspond to eigenstates under the semi-infinite boundary conditions with x ranging from 1 to L and L tends to infinity, provided |z s | < 1 [108].
To obtain the GBZ of the system, one needs to take into account OBC for the system, namely setting t j = 0 in Eq. ( 16) for x + j > L or x + j < 1.An eigensolution of the OBC Hamiltonian can be constructed from a linear combination of Ψ s (E), Ψ(E) = s c s Ψ s (E), which also needs to satisfy the OBC.Since L can be arbitrarily large, in order for there to always be at least two surviving eigensolutions, we arrive at the conclusion that a nontrivial solution of {c s } exists only when |z m (E)| = |z m+1 (E)| for some m, which limits the allowed values of eigenenergies E [51,54,106].These allowed values are precisely those that form the OBC spectrum.
Although the above conclusion does not directly give us an analytical relation between the imaginary flux κ(k) and the momentum k, but it already provides an efficient scheme for obtaining the GBZ numerically.That is, for any complex value E, the solutions z m , z m+1 of the characteristic polynomial equation f (z s , E) = 0 form a pair of GBZ solutions if they have the same absolute value.OBC spectrum of the system is then reproduced by collecting all values of E satisfying this condition.Alternatively, an analytical solution of the GBZ can be obtained by introducing a concept of auxiliary GBZ and solving the corresponding algebraic equation [134].The Hamiltonian H(k + iκ(k)), which is defined on the GBZ and which thus no longer experiences the NHSE [52,135], is known as the surrogate Hamiltonian.We also mention in passing that the difficult problem for solving for the GBZ and surrogate Hamiltonian may be circumvented by mapping the NHSE problem onto an equivalent electrostatics problem, which we will disucss in Sect.IV F. A modified GBZ theory has also been introduced to describe NHSE and bulkboundary correspondence in disordered systems [136].

C. Transition of topological localization direction and half quantized winding numbers
In addition to modifying the BBC, NHSE provides another mechanism to induce exponential localization of eigenmodes, and hence can alter the localization properties of topological edge modes [94-96, 137, 138], as already experimentally demonstrated in a mechanical setting [139,140].To see this, let us revisit the non-Hermitian SSH model in Eq. (1), which is topologically nontrivial when |t 1 | < (t 2 ) 2 + (γ/2) 2 .Within this parameter regime, a pair of zero-energy edge modes appear and localize at the two ends of the 1D chain [see Fig. 3(a)], analogous to the Hermitian counterpart.When increasing t 1 , non-Hermiticity becomes increasingly dominant, and one of the topological edge modes becomes extended at t 1 = 1.5 [Fig.3(b)], and localized at the same edge as the other for t 1 > 1. 5 [Fig.3(b)].
This transition of localization direction can be understood from a competition between NHSE localization and topological localization [95,96].This is because here NHSE induces a unidirectional exponential localization of eigenmodes, and topological localization induces a bidirectional one.Therefore one topological edge mode will be pumped to the other end of the 1D chain when the NHSE dominates.A transition point thus emerges where the two localization mechanisms are balanced, and this topological edge mode becomes delocalized [141].Generally speaking, this transition occurs if the eigenenergy of this edge mode coincides with the PBC spectrum of the same system [95,142].In the non-Hermitian SSH model, these topological edge states are pinned at zero energy as required by sublattice symmetry, As a result, the localization transition accompanies an EP of the two bands, which also represents a transition between separable and inseparable bands, as shown in Fig. 3(c) to (e).
From the topological aspect, this transition can also be related to a half-quantized winding number [33,95].In 1D Hermitian systems, a sublattice symmetry ensures a Z-class topology characterized by an integer winding number ν, which counts the number of times that the Hamiltonian vector h(k) in the momentum space encircling the origin of h(k) = 0. Explicitly, this winding number is defined as with tan ϕ = h y /h x for the Hermitian SSH model [i.e.Eq. ( 2) with γ = 0].Here the origin of h(k) = 0 represents a degenerate point of the two energy bands, which splits into two EPs in the pseudospin space when γ 0 in the non-Hermitian SSH model [pink star and red triangle in Fig. 3(g) to (i)].Therefore two winding numbers can be defined regarding each of the two EPs, with h x = h xr + ih xi and h y = h yr + ih yi .The total winding number is shown to be given by hence the winding number ν takes a half-quantized value when the trajectory of the real part of h(k), denoted as h r (k), encircles only one of the two EPs.[33].At the transition point between integer and half-quantized winding numbers, h r (k) passes one EP in the pseudospin space, resulting in an EP between the two energy bands.Thus the half-quantized winding number also characterizes the localization transition and band-structure transition in the non-Hermitian SSH model.In generic two-band models in the absence of symmetries, h r (k) of generic two-band models traces a curve in a 3D pseudospin space, where the two 0D EPs evolve into a 1D singularity ring [35]].The winding of h r (k) around an EP now becomes a linkage between two 1D trajectories.

D. Topological characterization without the GBZ
Even though one can in principle always (numerically) evaluate the GBZ and compute the correct topological numbers on the GBZ, for generic models this is in practice often fraught with subtleties due to GBZ singularities, as well as the amplification of floating-point errors due to the NHSE.Fortunately, for particle-hole symmetric systems of the momentumspace form where z = e ik , it is possible to determine whether topological zero modes exists just from the polynomial structure of a(z) and b(z) since only off-diagonal elements of the 2-by-2 matrix Hamiltonian exist, without computing the GBZ at all.This was developed in Ref. [54], which pre-dated the formal development of the GBZ.Here coefficients a i , b i are the complex roots of Laurent polynomials a(z), b(z), both of which can be rescaled without changing the topology.
Ref. [54] provided three equivalent formulations of the criterion for having topological zero modes.
1.The first formulation is a combinatorial condition on the coefficients a i and b i of the polynomials a(z) and b(z), which can take generic forms: An isolated topological zero mode exists if and only if the r a + r b largest coefficients do not contain r a members from {a 1 , ..., a p a } and r b members from {b 1 , ..., b p b }.
Even though this formulation is very explicit in terms of the roots a i , b i , as well as the pole orders r a , r b , its geometric and topological interpretation is not explicit.
2. The second formulation below recasts the above formulation in terms of topological windings.An isolated topological zero mode exists if and only if where the windings W a (R) and W b (R) are given by with g = a, b.Hence W g counts the number of zeros #Z g (R) minus the number of poles #P g encircled by a circle |z| = R of radius R ∈ R. Due to the skin localization in NHSE, we do not perform the winding integral on the typical |z| = 1 contour, i.e. like in the usual BZ.The criterion adroitly does away the need to evaluate the GBZ directly, which can have a radius that depends on arg(z) in an arbitrarily complicated manner.
3. The previous formulation can also be expressed in terms of the energy surface vorticity and eigenmode winding [34,50], such that the winding condition becomes From the above topological criterion, which is applicable to PH-symmetric models with arbitrarily complicated hoppings (the coefficients of a(z), b(z)), we see that while evaluating the full PBC spectrum requires knowledge of the GBZ, finding out whether they are topological only requires the above winding conditions.

III. NHSE IN HIGHER DIMENSIONS A. Chern lattices with NHSE
Lattice systems with Chern topology underlie the highly sought-after quantum anomalous Hall effect [143][144][145][146][147], and have thus been thoroughly studied and reviewed.As such, we will just briefly mention some recent advances.
As in other lattice NHSE systems, the effective band structure and its properties have to be evaluated in the GBZ, not the usual BZ.In the GBZ, the NHSE is effectively "gauged away", and the Chern number of a non-Hermitian Chern insulator [52,[148][149][150] can be computed as usual via the formula where the ket |ψ⟩ ranges over the right occupied states and ⟨ψ| is its corresponding biorthogonal bra state.What can be different from Hermitian cases is that the occupied band can have its own spectral singularities and topological configurations.
As such, the Berry curvature and band metric can be discontinuous in the GBZ, leading to semiclassical response kinks, as demonstrated in a model where the Chern bands exhibit singularities with 3-fold rotational symmetry in the complex spectral plane [52].
Even though non-Hermitian Chern bands can be mathematically described in the same way in their GBZ as Hermitian Chern bands in their BZ, their phenomenology exhibits important differences.In particular, because the biorthogonal bras and kets are not really Hermitian conjugates of the "same" states, under NHSE, a quantized topological invariant may not give rise to quantized transport: [151,152] gives a non-quantized contact effect from the NHSE.That said, it is still possible to formulate chiral anomalies rigorously for non-Hermitian systems [153], although a distinction must be made between directed state amplification by the NHSE, and chiral pumping by the chiral anomaly.The interplay between these two effects lead to the non-Hermitian chiral magnetic skin effect predicted by Ref. [154], which also extends the known Nielsen-Ninomiya theorem for nonchiral gapless fermions protected by symmetry.Furthermore, since the chiral boundary states of a topological insulator can possess their own GBZ, it is possible to design the chiral edge/hinge states of a non-Hermitian topological insulator passive without gain/loss, and hence immune to the NHSE [155].
Closely related to non-Hermitian Chern lattices are non-Hermitian lattices with magnetic flux.The key mathematical difference is only that of the size of the magnetic unit cell.NHSE lattices with flux exhibit exotic properties illustrating the interplay of time-reversal breaking and the NHSE.At the intuitive semi-classical level, the cyclotron trajectories of the wave packets in a 2D lattice always form closed orbits in four-dimensional (4D) phase space [156,157], and the semi-classical quantization rules remain valid despite the nonreciprocity, with the propensity to preserve real Landau levels [158].See [159] for a generic phenomenological approach to the interplay of flux and the NHSE, and [160] for a comprehensive treatment of the flux response under nonhermiticity via localized flux intersections on a single plaquette.Interestingly, skin states localized at the boundary can be pushed back into the bulk by an applied magnetic field, thereby leading to flux-suppressed NHSE [161].
Since wavepackets at the edge of a Chern lattice are steered by chiral topological pumping, through the Chern polarization, the extent of gain/loss experienced by an edge wavepacket can be controlled.Additionally, by weakly coupling two Chern topological layers with oppositely directed gain, it is possible to devise a topological guiding mechanism where there is asymptotic gain when a topological island becomes sufficiently wide.This has been shown to give rise to a percolation-induced PT transition in a disordered setting [162].

B. Non-Hermitian Weyl and exceptional metals
In a higher dimensional system, the gap between non-Hermitian bands can close along trajectories known as exceptional rings [163][164][165][166][167], in analogy to the nodal rings of Hermitian systems, some with non-abelian structure .Likewise, topological surface or edge states can appear at surface/edge terminations, since for a Hamiltonian H(k) = H(k ∥ , k ⊥ ), gap closures at certain k ∥ can behave like topological phase boundaries in the parameter space of k ⊥ .
When the NHSE is present under OBCs in the k ⊥ direction, the topological invariant in k ⊥ must be computed in the GBZ, and not the ordinary BZ.As such, the projection (shadow) of the exceptional trajectory onto the boundary may not necessarily correspond to the topological phase boundary on the surface, demarcated by the bulk nodal lines [123,193].This is schematically illustrated in Fig. 4a,d for the non-Hermitian Hopf link and trefoil knots.
Nodal or exceptional structures with intricate momentum space structures possess interesting topological linkages or knots.However, their models necessarily contain higher harmonics in momentum space, which correspond to longranged real space couplings.A common approach for realistic implementation of such Hamiltonians is through Floquet engineering, where quenching between different simple nodal/exceptional configurations can give rise to much more intricate Dirac/nodal/exceptional structures [123,[194][195][196][197][198][199][200][201][202][203].Floquet engineering can also generate gap closures on demand, such as non-Hermitian Weyl semimetals [204].Providentially, the inclusion of non-Hermiticity favourably enlarges the parameter space for nodal knot and facilitates their implementation via with Hamiltonians H(k) = h(k) • σ with local hoppings.This is schematically illustrated in Fig. 4c.The advantage conferred by non-Hermicity can be understood by considering the criteria for gap closures, which is more easily satisfied via not one but two conditions |Re h| = |Im h| and Re h • Im h = 0. Now, the aforementioned disagreement (if present) between drumhead surface states and their non-Hermitian analog fundamentally stems from the NHSE.In higher-dimensional set-tings, the NHSE manifests as a 'inward' compression of the drumhead surface states into a smaller area contained within the drumhead region.This smaller region is coined the 'tidal state' region as it can be derived from a marine analogy from Ref. [123], as described below.Crucially, the boundary of these novel 'tidal states' corresponds to the gap closures of the skin states.
To understand the mathematical and geometric foundation of this tidal analogy, we examine its derivation from Ref. [123].There, a formal criterion to determine the 'tidal state' boundaries was established for particle-hole symmetric Hamiltonians (of the form given by Eqn.24), which encompasses a large class of systems.Since k ∥ coordinates are just spectators when taking the OBCs along the k ⊥ direction, we may directly invoke the formulation (Eqn.25) from Ref. [54] (described in Sect.II D) to find topological zero modes [54].We restate this below: To illustrate this procedure, we graphically work out the topological phase boundary for the simplest non-Hermitian nodal knot metal -the Hopf link, given by the Hamiltonian h = 2(sin The winding numbers involved in the criterion for topological modes are where we have the 'tidal boundary' demarcated by the in- This region is plotted in Fig. 4b and is clearly distinct from the corresponding Hermitian drumhead surface state (Fig. 4a).
To illuminate the nomenclature of this tidal construction, we consider the trefoil knot.By graphically plotting the inverse skin depth κ = log |z| (Fig. 4e-f), we see that these surfaces intersect precisely along the 'tidal boundaries', i.e. the skin mode solutions experience gap closure, which facilitates topological phase transitions in momentum space.To reconcile with our PBC-OBC interpolation, we identify the 'sea level' in our marine landscape analogy as the Bloch states characterized by log |z| = 0, i.e. purely real momenta.This 'sea level' is non-unique due to the rescaling given by the non-reciprocal similarity transformation, i.e. fluctuations of the sea level analogous to tides.The well-known 'drumhead states' in the corresponding Hermitian system are what we call the 'beaches' and this coincides with the intersection trenches between the energy bands and the sea level, i.e. 'tidal' boundaries (Fig. 4e).This completes the marine landscape analogy from [123].
To gain a broader perspective of nodal knot metals, we highlight an elegant relation between vorticity, the complex energy bands, and the Seifert surfaces in non-Hermitian nodal knots.Also uncovered in Ref. [123], this stems from the  bulk-boundary correspondence between the OBC 'tidal' region shape and the bulk point-gap topology (measured by the vorticity).As illustrated earlier, the former is constrained within the interior of PBC bulk nodal loop projections.Along this surface-projected nodal line, we can construct a director and count the number of times any point is encircled anticlockwise by u.This is the half-integer and is thus a hallmark of skin gap closure.However, we remark that non-trivial vorticity is insufficient to guarantee skin gap closure, and one has to inspect the complex energy band crossings.
By considering any closed path in the trefoil knot's surface BZ, we uncover the PBC loci to be a cobordism of one or more conjoined tube/s along this chosen path, flanked by an interior skeleton of skin states (Fig. 4g).These tubes will join at their ends to form a Riemann surface indicative of the vorticity structure (Fig. 4h).Crucially, tidal regions are topologically constrained to contain islands of vanishing vorticity.
Lastly, to further highlight the topological significance of the 'tidal states', we can appropriately reverse the directors u at each crossing in the knot diagram, so as to construct a 'dual' knot that bounds a Seifert surface.The dual Seifert surface has an intimate relation to the knot topology.Remarkably, this resultant layer structure of the Seifert surface resembles our celebrated 'tidal islands'.All in all, non-Hermitian nodal knot metals are indeed rich in various forms of topology, and have been simulated on various experimental platforms [123,205].

C. Higher-order NHSE and hybrid skin-topological systems
Higher dimensions also offer a fertile setting for topological and NHSE states to hybridize.In purely Hermitian settings, higher-order topological insulators (HOTIs) [206] result from the nontrivial interplay of topological localization in 2 or more directions, arranged such that the first-order (1D) topological polarization cancels.In non-Hermitian settings, the NHSE provides another distinct avenue for localization [207].Hybrid skin-topological systems, first proposed in Ref. [99], result from the non-trivial interplay of topological and NHSE localization in different directions.Finally, one can also define higher-order NHSE systems [105] characterized by NHSE localizations in more than one direction.Such higher-order NHSE have been proposed in suitably designed non-Hermitian superconductors [208], and also lead to new avenues for disorder induced phase transitions in 2D [209].

Hybrid skin-topological modes -Original proposal in a non-reciprocal lattice
In two or higher dimensions, the non-Hermitian skin effect and Hermitian topological boundary localization can be treated on equal footing, leading to the hybrid skintopological effect in net reciprocal lattices [99].Simply put, such systems support opposite non-reciprocity for different sublattices, which cancel each other in each unit cell, resulting in extended bulk states free from the NHSE.On the other hand, topological boundary modes are usually sublatticepolarized, hence they experience a spontaneous breaking of reciprocity [210] and are further pushed to lower dimensional boundaries (e.g.1D edge states are pushed to 0D corners for a 2D system).At a more sophisticated level, the interplay of topological and skin localizations can be exploited for the engineering of "PT-activated" states [211].
This higher-order (lower-dimensional) boundary localization was first unveiled in a non-Hermitian extension of the 2D Benalcazar-Bernevig-Hughes model [206], which can also be viewed as a mesh of two different nonreciprocal Su-Schrieffer-Heeger model along each of x and y directions, as shown in Fig. 5(a) [99].In this model, non-Hermiticity is introduced by two non-reciprocal real hopping parameters δ 1 and δ 2 (δ 3 and δ 4 ) along x (y) direction, which balance each other and lead to a net reciprocity when they are chosen to be the same.For properly chosen parameters, topological boundary modes appear and distribute evenly along the two 1D boundaries of a cylinder geometry (with PBC taken along y direction) of the lattice [Fig.5(b)], which further evolve into 0D corner modes when OBC is taken along both directions (i.e. a double OBC) [Fig.5(c)], representing a NHSE acting on these boundary modes.On the other hand, bulk modes in this system remain extended and are free of NHSE.This can be seen from the spectral properties in 5(d) and (e), where bulk bands are almost identical for the system under different boundary conditions.Similar to 1D systems, it indicates that no imaginary flux shall be introduced when changing the boundary conditions, a sign of the absence of NHSE for the bulk modes.In contrast, as seen in 5(e), eigenenergies of boundary modes show a typical feature of the emergence of NHSE, i.e. they form a loop-like spectrum for x-OBC/y-PBC (cylinder geometry), which shrinks into some open lines under double OBC.
Inspired by this explicit 2D construction, hybrid skintopological phases can be straightforwardly generalized to higher dimensions, which support much richer classes of skintopological boundary modes, as each dimension can contribute skin (S) or topological (T) boundary modes, or neither (0).For instance, as shown in Fig. 5(f), a 3D system realizing the hybrid skin-topological effect can be obtained by stacking a series of the 2D model in Fig. 5(a).In different parameter regimes, this model supports all the three possible classes of hybrid skin-topological modes in 3D, namely the STT corner modes, SST corner modes, and ST0 hinge modes, as demonstrated in Fig. 5

(g) to (i).
Soon after its discovery, experiments of hybrid skintopological effect have been proposed using cold atoms loaded in optical lattices [100], and physically realized in circuit lattices [78].In the non-Hermitian BBH model, non-reciprocity and its resultant hybrid skin-topological effect are induced by asymmetric hopping amplitudes, which are difficult to realize in quantum platforms such as cold atoms.Alternatively, in cold atoms, non-Hermiticity can be experimentally implemented through an extra resonant optical beam transferring the atoms to an excited state, representing an atom loss for the concerned system [212].In particular, Ref. [100] considers a two-orbital two-sublattice optical lattice with an orbital-dependent atom loss, which interplay with other interand intra-orbital hoppings and induces an effective nonreciprocity along x direction.By considering a specific twotone shaking of the lattice, the non-reciprocity is induced toward opposite directions for the two sublattices, leading to a cancellation of non-reciprocity in the bulk.Thus hybrid skin-topological effect emerges in this system when topological boundary modes appear, which manifests as a topological switch to control NHSE along the system's boundary.On the other hand, in circuit lattices, non-reciprocal hoppings can be induced and precisely tuned by negative impedance converter through current inversion, allowing for a direct realization of the non-Hermitian BBH model and its 3D extension supporting different types of hybrid skin-topological modes [78].We note in passing that 2D higher-order lattice can also acquire their requisite coupling values through Floquet engineering [213,214], based on an approach first developed in 1D [215].

Hybrid skin-topological effect in a honeycomb lattice without asymmetric couplings
Recently, a class of hybrid skin-topological modes has been discovered in a non-Hermitian Haldane model [143] with gain and loss, without the need for asymmetric couplings [101,102].
As shown in Fig. 6(a), the model is given by a honeycomb lattice with real nearest neighbor hopping (t 1 ) and complex next-nearest-neighbor hopping with an amplitude t 2 and a phase ϕ.Complex on-site mass terms m ± iγ are added to the two sublattices, which interplay with the chiral edge currents induced by the flux of ϕ and lead to a non-reciprocal pumping along the system's 1D boundary.An interesting property of this model is that the hybrid skin-topological skin effect emerges only along the zigzag boundary, but not the armchair boundary.These two types of boundaries are demonstrated in Fig. 6(b) and (c) respectively.and the only difference between them is the direction of flux.For the armchair boundary, the next-nearest-neighbor couplings have phases with alternate signs, thus the fluxes of neighbor triangular plaquettes cancel each other, and the NHSE is absent.In contrast, the zigzag boundary supports a non-vanishing flux, which induces the hybrid skin-topological modes in the presence of sublatticedependent gain and loss, as shown in Fig. 6(d) to (f).Such a gain-loss-induced hybrid skin-topological effect is discovered independently in Refs.[101,102].The former focuses mainly on different types of boundary and a PT phase transition of skin-topological modes, and the latter extends their results to a Floquet realization of the hybrid skin-topological effect induced by gain and loss, and proposes an auxiliary Hermitian Hamiltonian of higher-order topological phases to understand this effect from a different angle.

Higher-order NHSE
When the NHSE localization occurs in more than one direction (i.e.SS or SST modes), we obtain a higher-order NHSE system [103][104][105]209], as experimentally reported in Ref. [79,80,216,217].In conventional (Hermitian) higherorder topological phases, e.g.d-dimensional nth-order topological insulators, topologically protected boundary modes appear at the (d − n)-dimensional boundaries, with their number scaling as O(L d−n ), as sketched in Fig. 7(a)(b).In contrast, nth-order NHSE d dimensions induces O(L d−n+1 ) number of skin modes localized at the (d − n)-dimensional boundaries, while bulk modes remains delocalized in the system.In this sense, the hybrid skin-topological effect discussed above also represents a type of higher-order NHSE, where higher-order skin modes are also protected by conventional topological properties.Note that the terminology "Higher-order NHSE" has also been used to describe different phenomena in different literatures.For example, in Ref. [99], it refers to corner NHSE where all bulk modes localize at a corner of a 2D lattice.
Conventionally, first-order NHSE (i.e.bulk NHSE) originates from from intrinsic non-Hermitian spectral winding topology of the Bloch Hamiltonian H(k) [106][107][108][109]. Alterna- tively, this non-Hermitian topology can also be mapped to the band topology of a Hermitian Hamiltonian in extended space: This extended Hermitian Hamiltonian H(k, E) respects chiral symmetry with a Pauli matrix If the spectral winding of H(k) for a given reference energy E r is topologically nontrivial, the skin effect occurs for the non-Hermitian system, and the extended Hermitian Hamiltonian H(k, E) also possesses nontrivial (Hermitian) topology which supports zero-energy edge modes under the open boundary conditions.It associates the zeros modes of the chiral-symmetric Hermitian Hamiltonian with a nonzero spectral winding number of H(k) [108].With this mapping between Hermitian and non-Hermitian topology, a higherorder NHSE can be mapped to a chiral-symmetric Hermitian Hamiltonian with higher-order Hermitian topology.As higher-order topological phases are usually protected by spatial symmetries (e.g., inversion, mirror, and rotation symmetries) [108,[218][219][220], higher-order NHSE originated from intrinsic non-Hermitian topology can also be associated with these symmetries.Below, we briefly introduce another model proposed in Ref. [103], which exhibits the second-order NHSE.Its Bloch Hamiltonian reads where γ and λ are real parameters, and σ i 's(i = x, y, z) are Pauil matrices.Its extended Hermitian Hamiltonian is given by where τ i 's(i = x, y, z) are Pauli matrices that describe an additional pseudospin-1/2 degree of freedom.This Hamiltonian is known as the 2D Benalcazar-Bernevig-Hughes model and describes a second-order topological insulator [206]).The connection between the second-order NHSE in Eq.( 33) and second-order Hermitian topology in 2D BBH model can help us understand the origin of the second-order non-Hermitian skin effect.For example, in 2D BBH model, under the open boundary conditions along both directions, zero-energy modes appear at the corners for |γ/λ| < 1. Correspondingly, in the model of Eq. ( 33), eigenmodes isolated from bulk bands emerge when |γ/λ| < 1 [Fig.8(a)], which localize at two corners of the 2D lattice [Fig.8(b)], manifesting the second-order NHSE.Meanwhile, eigenmodes in the bulk bands remain delocalized in the bulk, as shown in Fig. 8(c).According to Ref. [103], the second-order NHSE in this model is protected by a rotation-type symmetry, represented by Based on this model, Ref. [221] shows that a perpendicular magnetic field can enhance this second-order NHSE, or even induce it in the otherwise trivial parameter regime.

D. Other NHSE phenomena unique to higher dimensions
Here, we review some other interesting phenomena resulting from the interplay of the NHSE with the enlarged spatial degrees of freedom in higher dimensions.

Translation-invariant bulk
We first discuss scenarios where the skin modes accumulate at the boundaries of a clean D-dimensional lattice.Implementing the D different OBCs sequentially, one can generalize the known result that the NHSE occurs if and only if the PBC spectrum covers nonzero area (i.e.possess nontrivial spectral winding) in the complex energy plane [106][107][108][109].
However, the generalization in Ref. [72] may fail to apply when the operation of taking OBCs in the various directions do not commute.This can happen, for instance, when the lattice cannot be "disentangled" or broken down into the sum of 1D NHSE lattices [222].In such scenarios, the spectrum obtained by first taking the GBZ in, say, the x-direction, and then followed by the iterated GBZ in the y-direction, is different from that obtained by taking the y-and then x-direction GBZ.This indicates a breakdown in the conventional paradigm of taking 1D GBZs, forcing us to take a different approach that considers a fundamentally multi-dimensional GBZ, as introduced by Ref. [222].In Ref. [222], it was further shown that a multi-dimensional nontrivially "entangled" lattice can in general possess a "transmutated" or dimensionally-reduced GBZ which consists of a union of low-dimensional subspaces of the GBZ under Hermitian settings.
Further violations of the nonvanishing spectral area criterion [72] occur when there are GBZ singularities, such that certain pairs of biorthonormal eigenvectors become orthogonal.This can allow an open boundary to change the spectrum from real to complex, even though it is usually the case that an open boundary makes a complex spectrum real due to arrested NHSE pumping.As a specific example, the 2D model with "rank-2" chirality [223] exhibits a real spectrum when open boundaries are opened in one direction, but the spectrum becomes complex upon opening boundaries in both lattice directions.A qualitatively similar observation was also observed at sufficiently high dimensions in Ref. [224].

Non translation-invariant bulk
We now discuss the NHSE arising right in the "bulk" due to translation breaking.In general, disorder such as dislocations can act like localized defect boundaries, such that the NHSE causes skin mode accumulation against them [70,71], albeit not necessarily in the same way as open boundaries, which possess translation invariance in the directions within the boundaries [73].
Suitably designed disorder can even admit skin modes of all decay lengths, such that they are described by a large range of Im(k).This could allow for the disordered spectrum to "fill up" the interior of the PBC spectrum loop [225], recapitulating the scenario of semi-infinite boundaries which also admit skin modes of all decay lengths [226].
Lattices with fractal structures poses the interesting scenario where the dimensionality may not even be welldefined [74,[227][228][229][230]. in Ref. [74], 2D fractal lattices are shown to host "inner skin effects", where details at different scales lead to different levels of skin mode accumulation and thus different levels of skin localization.Interestingly, mirror crystalline symmetry interplays with NHSE to result in sensitivity of the energy spectrum on the boundary condition only along mirror invariant lines -the mirror skin effect [231].
We will revisit the interplay of NHSE with breaking of translation invariance in Sec.VII A.

IV. TOPOLOGY ON THE COMPLEX ENERGY PLANE A. Point-gap and line-gap topology
Since gap closing is typically associated with a topological phase transition, it is important to provide a precise definition of energy gaps in different systems.Here we shall briefly review the definition of point gaps and line gaps for non-Hermitian Hamiltonians in Ref. [226].It is clear that for a Hermitian Hamiltonian, its real spectrum is defined to have an energy gap if and only if its energy bands do not cross a zero-dimensional point E = E F , which is called the Fermi energy [see Fig. 9(a)].In contrast, the spectrum of a non-Hermitian Hamiltonian spans the 2D complex energy plane, hence a complex-energy gap can be either a zero-dimensional point or a one-dimensional line.Similarly, a non-Hermitian Hamiltonian can be defined to have a zero-dimensional point gap (one-dimensional line gap) if and only if its complexenergy bands do not cross a reference point E = E P (reference line) in the complex plane, as shown in Fig. 9(b) and (c).More precisely, taking into account the restrictions on eigenenergies from different symmetry (e.g., eigenenergies come in (E, E * ) pairs due to the time-reversal symmetry), it is convenient to choose E P = 0, and a non-Hermitian Hamiltonian H(k) is defined to have a point gap if and only if it is invertible (i.e.,∀k det H(k) 0) and all the eigenenergies are nonzero (i.e.,∀k E(k) 0 ).Thus, a non-Hermitian Hamiltonian H(k) is defined to have a line gap in the real (imaginary) part of its complex spectrum, denoted as a real (imaginary) gap, if and only if it is invertible (i.e.,∀k det H(k) 0) and the real (imaginary) part of all the eigenenergies is nonzero These enriched definitions of energy gaps have led to several unusual topological properties for non-Hermitian systems.Intuitively, if a non-Hermitian Hamiltonian has a point (line) gap, it can be continuously deformed into a unitary (Hermitian/anti-Hermitian) matrix while keeping the point (line) gap and its symmetries, which means that the two Hamiltonians before and after the deformation are topological equivalent [see Fig. 9(b)(c)].These properties have been proven in Ref. [226], which play a crucial role in the complete topological classification of 38 symmetry classes for point gaps, and 54 symmetry classes for non-Hermitian line gaps [232] and Floquet systems [233].The studies of point-gap topology in non-Hermitian systems have also be generalized to systems with crystal symmetries such as an inversion symmetry, where topological invariants can be determined from high-symmetric momenta [234,235].
As another example, an immediate topological description of a point gap is the spectral winding number defined as in Eq. (35), which vanishes together with the point gap when Anderson localization is induced to the system by spatial disorder [236].A nontrival spectrum winding number associated with a point gap under PBCs also result in NHSE under OBC [106,108], as discussed in Sec.IV B.
In 3D topological phases, it has also been suggested that topological surface states originates from a 3D winding num- ber under PBC, which is the point-gap topological number for general 3D systems [235,237].On the other hand, a more recent study unveils that a system with a nontrivial 3D winding number may have distinguished behaviors for boundary states with OBC along different directions, as shown in Fig. 10 [238].Following the idea of real space topological invariants for non-Hermitian topological systems [239][240][241][242], Ref. [238] defines real space winding numbers for different OBCs, which establishes the bulk-boundary correspondence in the point-gap topology (with certain symmetry protection) of 3D non-Hermitian systems.

B. Topological origin of NHSE
We next describe an important connection between nontrivial spectral winding, which leads to spectral loops of nonzero area, with the presence of the NHSE in 1D non-Hermitian Hamiltonians.We define a spectral winding number with H(k) a non-Hermitian Hamiltonian, and E r a chosen complex value as a reference energy.In contrast to the winding number of h r (k) in the pseudospin space discussed in the last subsection, a spectral winding number describes the number of times that the complex spectrum of H(k) encircles E r anti-clockwisely throughout the Brillouin zone.Mathematically, w(E r ) can be related to the total number of zeros and poles of H(z) − E r with z = e ik , enclosed by the path of z on the complex plane with k varying from 0 to 2π, i.e. the BZ.Explicitly, where N zeros and N pole are the counting of zeros and poles weighted by their respective orders.Replacing H(k) with H(z GBZ ) = H(k + iκ(k)) in Eq. ( 35), a winding number can be defined for the GBZ, which has been proven to always enclose m zeros and a pole of order m, resulting in a zero winding number for the GBZ spectrum (equivalent to the OBC one exact for topological edge states) [106].This means that if a system supports nonzero spectral winding under the PBC, its GBZ must have κ(k) = 0, so that OBC eigenstates have nonzero decaying rates and exhibit NHSE.Alternatively, this topological bulk-boundary correspondence between OBC NHSE and PBC spectral winding number has also been unveiled with a doubled Green's function approach [107], PBC to OBC interpolations [243] and detailed studies of Toeplitz matrices [108].Such spectral properties of systems with NHSE have also been studied by mathematicians decades ago [244].
Further insight into the NHSE in terms of the detailed structure of the eigenstate trajectories, not just the eigenvalue trajectories, can be obtained through the use of Majorana stars [245,246].A N-component state vector can be mapped onto N −1 Bloch spin vectors (stars), such that its complex degrees of freedom can be visualized in terms of the real angles and real correlations between these N − 1 vector directions.These visualizations provide additional geometric and topological pictures into generic multicomponent non-Hermitian models.

C. Quantized response from spectral winding topology
As discussed in the previous Sec.IV B, a nontrivial spectral winding is responsible for the emergence of NHSE under the OBCs, representing a new type of bulk-boundary correspondence of spectral winding topology.A one-on-one correspondence has been established between the spectral winding and number edge states under the semi-infinite boundary conditions [108], which naturally arises for the Fock space mapped onto a real-space lattice [247].Experimental observation of arbitrary spectral winding has also been carried out by visualizing the frequency band structure of optical frequency modes [248].
In the framework of conventional band topology, Hermitian topological systems exhibit quantized physical phenomena that hinges on their respective topological invariants.Celebrated examples include the quantized charge transport in Thouless pumps for 1D systems, and the quantum Hall effect for 2D systems.It is thus natural to ask if non-Hermitian systems exhibit quantized response?Recently, a topological quantized response with a one-on-one correspondence to non-Hermitian spectral winding topology has been proposed by considering the Green's function of the Hamiltonian with PBC-OBC interpolations, with E r a chosen reference energy, H(β) the Hamiltonian, and β a parameter controlling the boundary conditions [109].
Namely, for a general single-band Hamiltonian of Eq. ( 16), the hopping parameters are set to t j → t j e −β for x + j > L or x + j < 1. Therefore the system is under the PBC when β = 0, and OBC when β tends to infinity.Note that a finitesize non-Hermitian system behaves as under the OBC when β exceeds a finite value associated with the non-Hermitian parameters [63,65].Finally, for a spectral winding number w(E r ) = m, the quantized response is defined as , where clear quantized plateaus from 2 to 0 are seen in the latter.The correspondence betwen ν 2 and the spectral winding regarding the chosen reference energy (red star) is further verified by comparing with the complex spectra in Fig. 11(c).Another way to demonstrate the topological quantized response is to consider ν m at β ≳ 0, which reflects the spectral winding topology of the PBC system, as shown in Fig. 11(d).
Physically, the quantized response quantity ν m can be associated with the directional signal amplification for non-Hermitian systems [58][59][60].Namely, for a local steady-state driving field (input signal) ϵ x (ω) with a frequency ω at a given location x, a response field (output signal) ψ x ′ (ω) can be obtained as [60], with G xx ′ an element of the Green's function.In other words, G xx ′ describes the amplification ratio between sites x and x ′ , thus the quantized response quantity ν m describes the changing rate of this ratio between the first and last m sites of a 1D chain, during a PBC-OBC transition of the system.Alternatively, ν m may also be detected by measuring the two-point impedance Z xx ′ between the two sites in circuit lattices, which is related to the Green's function via We end the section with a few remarks.First, instead of turning off the boundary coupling, a local on-site potential can also act as a boundary of the system, and leads to the same quantized response [250].
Secondly, in a system with two weakly coupled non-Hermitian chains, it is found that spectral winding topology of one chain can be detected by the quantized response solely of the other chain, reflecting an anomalous hybridization of spectral winding topology of the two chains [110].
Thirdly, nontrivial spectral winding has also been found to emerge for boundary states in 2D lattices, which can also be captured by the quantized response defined for the boundary of the systems [251].Physically, the boundary spectral winding originates from the interplay between non-Hermitian nonreciprocal pumping and conventional topological localization, analogous to the mechanism behind hybrid skin-topological effect introduced in Sec.III C 1.
Fourthly, while the Green's functions defined above are designed to capture signatures of the NHSE, bulk Green's functions in general i.e. [252] do not necessarily capture the onset of NHSE pumping [125], at least for a sufficient large system [253].

D. Complex band evolution as braiding processes
Another topological feature that can arise from the complex eigenenergies of non-Hermitian systems is the braiding between different energy bands i.e. braiding between the trajectories of complex eigenenergies as a parameter [254,255].Braids can be closed to form knots [256], which exhibit extremely rich topology since there is an infinite number of knot configurations that cannot be deformed into each other [194,200,[256][257][258].
Unlike spectral winding or point-gap topology, nontrivial braiding emerges only in non-Hermitian multiband systems with separable bands, i.e.E i (k) E j (k) for all band indices i j and momentum k.Conceptually, as a 1D momentum k varies from 0 to 2π, eigenenergy trajectories of different bands may wind around each other and form a "braid".The simplest example is that the trajectories of two bands exchange once and connect to each other, giving rise to a energy vorticity unique in non-Hermitian systems [34].
Different braids of these trajectories cannot be continuously deformed into each other without touching between different bands (usually some EPs), representing a topological feature unique in non-Hermitian multiband systems.Due to the periodicity of Brillouin zone, eigenenergy trajectories of separable bands are closed loops, hence their braids can be mapped to knots in the 3D energy-momentum space.Several examples of nontrivial energy braids and knots are demonstrated in Fig. 12 adopted from Ref. [254].In addition, Ref. [254] also develops an algorithm to construct tight-binding Hamiltonian for any desired knot, and propose a scheme to probe the knot structure via quantum quench.Such nontrivial braiding of non-Hermitian Bloch bands has been experimentally observed in coupled ring resonators with phase and amplitude modulation, by extracting complex band structure from measured transmission signals measured from the resonators, and reproducing the complex spectrum of a non-Hermitian lattice FIG.12. Several examples of topologically distinct braids, their knot closures and their realizations in 1D non-Hermitian Bloch bands.τ i and τ −1 i denote braid operators of the ith string crossing over and under the (i + 1)th string from the left, respectively.Braid diagrams in the second row are mapped to the knots in the third row when connected top to bottom.The fourth row demonstrate typical structures of energy bands for these braids and knots in the 3D space of momentum and complex energy.The figure is adopted from Ref. [254].
system in a frequency synthetic dimension [259].Braiding of exceptional arcs have also been demonstrated [260,261].
In a more recent study, it has been shown that nontrivial braiding also exists for non-Bloch bands describing non-Hermitian systems under OBC [255].The braiding of non-Bloch bands can be topologically different from that of Bloch bands of the same system, since the presence of the NHSE leads to different band structures under PBC and OBC.It also further reveals that the spectral winding and braiding of separable bands are essentially different types of topology, as the former is always trivial for non-Bloch bands, despite that they both originate from the paths of eigenenergies changing with momentum k in the complex energy plane.

E. Emergent spectral graph topology
While the complex spectrum under PBCs is characterized by winding numbers as the momentum k is traversed over a period from 0 to 2π, it generically collapses into lines, curves, closed loops and branches as OBCs are introduced, as shown in Fig. 13a.As such, OBC NHSE spectra are generically characterized by spectral graph topology, which remain invariant under conformal transformations in the energy spectrum as first systematically studied in Ref. [52], and further characterized and classified in Ref. [134,135,262].This spectral graph topology on the complex energy plane is an emergent feature of NHSE with no Hermitian analog.
OBC NHSE spectra generically assume graph-like structures because they are the "shrunken" versions of PBC spectral loops.This can be understood through complex flux pumping arguments, as initially suggested by Ref. [48], and subsequently expanded on in Ref. [54].As explained in Ref. [54], let us first interpolate between PBCs and OBCs by multiplying the boundary couplings by e iϕ , ϕ ∈ C. For real ϕ, this is equivalent to threading a real flux through the lattice, as can be seen via a gauge transform; if we let Im(ϕ) → ∞, we obtain the OBC limit.Intuitively, the effect of threading Re(ϕ) → Re(ϕ) + 2π diminishes as Im(ϕ) increases, since the boundary coupling would be exponentially suppressed in magnitude and its phase should therefore exhibit diminishing influence on the whole system.This is illustrated in Fig. 13(b) for the illustrative model with next-nearest-neighbor (NNN) hopping (Eqn.39), which hosts a richer graph structure compared to the prototypical models with only one non-reciprocal length scale like the non-Hermitian SSH model.
where σ ± = (σ x ± iσ y )/2 and z = e ik as usual.When there is no suppression to the boundary coupling i.e.Im(ϕ) = 0, we have the PBC loop, and threading the real flux over one cycle maps one eigenvalue to the next.This cyclic permutation of the eigenvalues still occurs when Im(ϕ) > 0, but because the boundary couplings are now much weaker, the eigenvalue flow should also be correspondingly diminished: as seen in Fig. 13(b), this is indeed achieved by having smaller, "shrunken" loops in the interior of the original PBC loops.Extrapolating, we expect that in the OBC limit, the effect of pumping real flux Re(ϕ) should entirely vanishand this can only be possible if the spectral loop have somehow "shrunken" until they become degenerate i.e. enclosing zero area, such that for every eigenvalue, there exists another one infinitesimally close in the complex plane.This spectral flow can also be represented as a spatial flow along the eigenspectra surfaces (of generalized boundary conditions) as we interpolate from PBC to OBC, as shown in Fig. 13c.Crucially, at the OBC limit, the spectra is the intersection of two |z| surfaces in the complex plane (Fig. 13d) since in order to simultaneously satisfy the OBCs at both ends where the wavefunctions vanish, the OBC eigenstate should be a superposition of two or more degenerate generalized Bloch solutions that decay equally fast.
1.The characteristic polynomial and corresponding allowed spectral graphs Spectral graph topology was first briefly studied in Ref. [52], where the branching pattern of simple non-Hermitian models is related to the number of coexisting nonreciprocal length scales.By starting from the simplest cases containing one or two non-reciprocal length scales, the OBC spectra are analytically worked out via a complex momentum deformation k → k + iκ(k).In this representation, the spatial eigenmode accumulation due to NHSE is nullified via a spatial basis rescaling from the original Hamiltonian H(k) to the surrogate Hamiltonian H(k) = H(k + iκ(k)).By doing this, the PBC spectra of H(k) will recover the desired OBC spectrum of the original Hamiltonian H(k), at the expense of inducing an emergent non-locality in real space as a trade-off.Concretely, for the simplest non-trivial case with two nonreciprocal length scales, we have the bivariate characteristic polynomial to be which we can solve for the eigenspectra ϵ(k) of the surrogate Hamiltonian H to satisfy where the graph topology adopts the shape of three straight lines radiating from the origin of the complex F(E) plane aligned along the cube roots of unity ω j .With this approach, the smallest complex deformation |κ(k)| needed to recover the surrogate Hamiltonian H(k) may also be worked out: where j is chosen to give the branch the smallest complex deformation |κ(k)|.The key takeaway is that the OBC loci ϵ(k) is only dependent on the structure of the couplings (e.g.number of non-reciprocal length scales, coupling parameters) and independent of the form of F(E) (which tells us e.g.number of bands).By extending this to generic non-reciprocal couplings of the characteristic polynomial form with p, q > 0, the OBC spectrum ϵ(k) can be shown numerically to take the shape of a N(p + q)-pointed star, i.e. the Brillouin Zone is folded p + q number of times.A systematic analysis was subsequently done by Ref. [135], where a very rich graph topology of generic bounded non-Hermitian spectra is uncovered, distinct from the topology of conventional band invariants and spectral winding familiar in both Hermitian and non-Hermitian settings.Here, the goal was to uncover the deep mathematical relationships between spectral graph topology of non-Hermitian systems and the algebrogeometric properties of the energy-momentum dispersion, i.e. the bivariate Laurent polynomial P(E, z), which also control Table I.Different forms of the canonical dispersion P (E, z) of Eq. 5 correspond to rich and diverse OBC spectral graphs Ē.
For each P (E, z), we can 1.associate a non-unique minimal Hamiltonian H(z), 2. identify emergent global symmetries of Ē not necessary present in H(z), and 3. characterize its spectral graph topology with its number of branches NS, NL and loops N , as well as its adjacency matrix.See the Supplement [65] for the latter as well as more examples with varying r.
1.At small |E| and large |E|, the OBC spectral graph of P (E, z) have respectively N S and N L rotationally symmetric branches centered at the origin.The rest of the spectral graph interpolates between them.

Generally, the origin is the only branching point if
).An exception might occur if the equal number of spectral branches at small and large |E| are displaced by a small rotation angle.Additional branches appear to connect these two regimes, typically leading to a flower-like shape (Table I(ii)).
3. When N S < N L , additional disjointed or isolated branches can appear at larger |E| (Table I(iii)).However, this condition alone does not guarantee the existence of isolated branches (Table I(iv) is connected).
4. Depending on the symmetry of P (E, z) under r → −r, the spectral graph at sufficiently large ±r are either identical, or mirror-reflected (r-reflection).
5. As we tune |r|, the number of rotationally symmetric branches interpolates between N L at large |r| and f (q + + q − ) at small |r|, as exemplified in [65].
6. Sometimes, branches emanating from the origin may join up into loops.This is especially common when FIG. 14. Classification table to illustrate the complexity of spectral graph topology, as adopted from [135].For each form of the canonical dispersion P(E, z), one can 1.associate a non-unique minimal Hamiltonian H(z), 2. identify emergent global symmetries of the spectral graph Ē not necessarily present in H(z), and 3. characterize its spectral graph topology with its number of branches N S , N L and loops N ℓ , as well as its adjacency matrix.
the localization of Wannier functions [263][264][265][266][267]. A more generic and sophisticated energy dispersion of the form was considered, containing a term involving the product of an energy-band related term E and the momenta z = e ik .Eqn. 41 is sufficiently generic and encompasses a wide class of non-Hermitian Hamiltonians since a vast group of previously unrelated Hamiltonians are now tied together via a conformal transformation in the complex energy E → f (E) for some analytical function f .Here, a kaleidoscope of interestingly shaped spectral graphs resembling stars, flowers or insects was uncovered.The simplest examples are tabulated in Fig. 14, together with the minimal non-unique Hamiltonian, its adjacency matrix representation, as well as, the emergent global symmetries of the eigenspectra.
Similar to how conventional eigenstate topology manifests as linear response quantization, the topological transition between different spectral graphs physically manifests as linear response kinks [52,268], with the different parts of the eigenstates mixing abruptly, leading to enigmatic gapped marginal transitions with no Hermitian analog, giving rise to emergent Berry curvature discontinuities [52,268] with physically measurable response signatures, as elaborated in [268] for an ultracold atomic system.

F. Electrostatics approach to solving the NHSE problem
Ref. [93] established a correspondence between the NHSE problem and the age-old problem of the electrostatic field of charge configurations.By mapping a NHSE problem onto an electrostatics problem, one can circumvent direct numerical evaluation of the NHSE spectrum, which for sophisticated spectral graphs may quickly require too many real-space sites, leading to rapid accumulation of floating point numerical errors.This also circumvents the difficulty associated with complicated and perhaps unsolvable algebraic equations to determine the GBZ by completely doing away with them, and instead only requires solving a simpler boundary-valued Poisson equation (Fig. 13e-f).
Specifically, Ref. [93] demonstrated the correspondence between (i) PBC spectral loops with equipotential conductors, (ii) OBC spectral eigenvalues with electric charges, (iii) PBC-OBC spectral flow as electric field lines and (iv) density of states in the complex E plane with charge density.These relations follow from identifying the inverse skin depth κ with the electrical potential V,as illustrated in Fig. 13.
The duality between charge density and spectral density can be understood by considering the DOS along an arbitrary curve ϵ in the complex E plane of a lattice with L sites, obtained via the Cauchy-Riemann relations: This bears mathematical resemblance to the induced charge density on a plane with discontinuous field strength This elegant analogy allows us to tackle the difficult inverse problem -engineer a non-Hermitian Hamiltonian H(k) with desired OBC spectral properties and desired spatial profile κ(k).A toolbox of familiar methods can be employed here -superposition of point charges and the method of images, further enriching the utility of this approach by including scenarios with non-Bloch band collapse [269], etc.The electrostatic analogy also provides intuition to the phenomenon of non-Hermitian pseudo-gaps and pseudobands [270]: Under the NHSE, the PBC and OBC bands may not be in one-to-one correspondence.To understand how, consider the spectral flow as PBCs are deformed into OBCs.In the electrostatics picture, the spectral flow corresponds to electric field lines, and for conductor geometries that are "too sharp", it in conceiveable that neighboring field lines may diverge and group themselves into different pseudo sub-bands.Physically in finite-size OBC systems, this may result in the appearance of topological edge states and bands that ostensibly do not correspond to the topological indices [270].

V. EMERGENT CRITICALITY FROM THE NHSE
Beyond significantly modifying the band structure and topology of a system, the NHSE also introduces a new length scale, the skin decay length κ −1 .This extra degree of freedom nontrivially affects the behavior of critical systems, as we will review below.

A. Critical non-Hermitian skin effect
Intriguingly, when two non-Hermitian systems with different inverse skin lengths (or more generally GBZs) are coupled together, a novel critical behaviour is observed -the critical non-Hermitian skin effect (CNHSE).First introduced in Ref. [66], the notion of CNHSE arose after the celebrated GBZ formalism (which we previously established in Section IIB to restore the conventional bulk-boundary correspondence) which only holds in the thermodynamic limit, was rigorously challenged.Fundamentally, the non-Hermiticity effects contribute their share of long-ranged influences which is crucial in critical phenomena.Consequently, critical skin states can even exhibit scale-free behavior while decaying exponentially in space, contrary to conventional critical states which are almost synonymous with power-law spatial decay.They also possess unusual size-dependent entanglement entropy behavior, which challenges the usual approaches for characterizing critical states through their entanglement entropy scaling [66].
Although the GBZ formalism holds in the thermodynamic limit, the spectra agreement with that of finite systems only holds far away from the critical point.At these critical points, the characteristic polynomial f (E, z) cannot be reduced to two systems because the two very different subsystems are coupled [66].This is concretely illustrated with a minimal model of two coupled non-Hermitian Hatano-Nelson chains with only non-reciprocal nearest-neighbor hoppings, described by the Hamiltonian and t ± a/b = t 1 ± δ a/b , as illustrated in the schematic Fig 15a .When the chains are decoupled, i.e. t 0 = 0, the characteristic polynomial f (z, E) = (g a (z) − E)(g b (z) − E) is reducible such that each factor of f (z, E) determines the skin eigensolutions of its respective chain.Yet, when the coupling is switched on, t 0 0, f (z, E) is no longer irreducible.For the simplest case where t + a = t − b = 1 and t − a = t + b = 0, we have the irreducible bivariant characteristic polynomial The resulting eigenenergy roots E = cos k ± t 2 0 + (V + i sin k) 2 are no longer Laurent polynomials in z = e ik that can be separately interpreted as de facto subsystems with local hoppings.To obtain the OBC spectrum (in the thermodynamic limit), we set FIG. 15.Emergent criticality in finite non-Hermitian lattices, as adopted from [66].(a-c) Critical behaviour in asymmetrically coupled systems involving two non-Hermitian systems with different inverse skin lengths.As the size of the lattice increases, the open boundary spectrum transitions from the decoupled thermodynamic limit (green spectra) to the coupled thermodynamic limit (red spectra).(d-g) Similarly, a critical topological phase transition can be obtained with cross inter-chain non-reciprocal couplings as the size of the lattice increases.This is robust with exponentially weak, but non-zero inter-chain coupling.
which does not reduce to the above-mentioned OBC spectrum of the two decoupled chains in the limit t 0 → 0 (Fig. 15b).
Likewise, the t 0 → 0 limit of the coupled GBZ loci are qualitatively different with the collapsed GBZs of the decoupled case (Fig. 15c).The corresponding OBC E ∞ spectrum and the GBZ for t 0 0 are qualitatively different.
In critical systems such as the above example, the eigenstates are formed from the superpositions of eigenstates from dissimilar subsystems.This can be understood from a more intuitive perspective.In the GBZ picture, the physical local hoppings are replaced with effectively non-local ones so as to 'unravel' the real-space eigenstate accumulation due to NHSE [52].In other words, the NHSE "renormalizes" the hopping strengths such that they increase dramatically with system size, such that the same bare physical couplings can be tuned into the strong or weak coupling regimes just by changing the system size.While mathematically, the CNHSE arises when the energy eigenequation f (E, z) exhibits an algebraic singularity involving dissimilar auxiliary GBZ across the transition, physically, this manifests as a discontinuity in eigenenergies and eigenstates in the thermodynamic limit.In physical finite systems, this discontinuity would have to be manifested as some type of finite-size scaling behavior.
As a result, the spectrum exhibits a strong finite-size scaling and the simplistic GBZ picture no longer holds.The scaling rule was analytically worked out in Ref. [67] with a minimal model and its universality was demonstrated for multiband models [67].Exact solutions exhibiting boundary scaling behaviour were also worked out in Ref. [64].Here, they transcended the difficulties and ambiguities presented in conventional numerical methods and via their analytical results, they uncovered the origin of size-dependent NHSE and quantitatively demonstrated the interplay effect of boundary hopping terms and lattice size [64].Recently, the scaling rule was shown to apply much more generically [271], with the GBZ shown to explicitly depend on the system size N according to |β| ∼ |β| N→∞ + b 1/(N+1) − 1, b a function of the model parameters.Finite-size spectral properties have also been noticed earlier in Ref. [137], where the energy gap of a non-Hermitian SSH chain exhibits an oscillating exponential decay or a realimaginary transition as the system's size grows, depending on how Hermiticity is introduced to the system The paradigmatic example of two coupled dissimilar non-Hermitian chains is further studied in detail in Ref. [272] by considering the interplay with inter-chain coupling and different types of skin mode localization.Moreover, topologicallyprotected zero modes arise, even when the individual chains do not harbor such zero modes, and exhibit critical phenomena as well.This was also analytically studied in detail, with a proposal to realizing it in topolectrical circuit lattices [273].A proposal to realize the CNHSE in open quantum systems was given in Ref. [68], by explicitly considering the Lindblad master equation.This work unraveled the Z 2 skin effect from the CNHSE, showcasing how both the dynamical CSE and the anomalous CSE arise from the modified GBZ equation.
CNHSE is most saliently revealed in size-dependent topological phase crossovers, where the system only exhibits topological modes at certain system sizes.To study that, one can build upon the prototypical example -the two-chain model (Eqn.44), but instead of non-reciprocal intra-chain couplings, we have non-reciprocal inter-chain couplings between adjacent unit cells, which is illustrated schematically in Fig. 15d, and described by the Hamiltonian: where . Similar to before, the OBC spectra transitions discontinuously from the decoupled spectra to the coupled spectra as N increases (Fig. 15e,f).This manifests as a gap closure before the emergence of a point gap with two zero-energy degenerate modes lying in its centre -a paradigmatic example of a topological phase transition, but with an intriguing size-induced effect.Notably, this phenomenon is robust for exponentially weak inter-chain coupling for sufficiently large N, as illustrated in the phase diagram in Fig. 15g.

B. Exotic non-Hermitian critical behavior
In general, criticality occurs whenever the bands become gapless.In non-Hermitian systems, gapless points can be more interesting either because of the richer variety of gaps (i.e. point and line gaps), or due to the defectiveness of exceptional gapless points.
The critical properties of non-Hermitian gapless points have been extensively studied, for example, their fidelity susceptibility [274], disorder effects [275] and their thermodynamic scaling [276].The criticality of two paradigmatic models, the extended Non-Hermitian SSH [277] and Kitaev models [278], were studied concretely by tracking the evolution of the gapless zero energy edge states.In turn, this unravels the relation between EPs and criticality as gap closing points are associated with the appearance of EPs, which in interacting contexts are deeply related to non-unitary conformal field theories (Yang-Lee singularities) [279,280].In Ref. [270], it was also noticed that certain non-Hermitian gaps may appear to host topological in-gap modes even when the topological index is trivial, due to the phenomenon of non-Hermitian pseudo-gaps.
Other than explicitly working out the critical boundaries of the phase diagrams of various paradigmatic models, a natural approach to describe criticality would be to use renormalization group (RG).Yet, the presence of non-Hermiticity poses considerable challenge to applying conventional RG theory directly since conventional RG flow may drive a critical state towards a non-critical state.To transcend this difficulty, the work [281] proposes a novel real-space block decimation RG scheme which is much more natural given that the critical hypersurfaces of non-Hermitian systems are obtained in real space under open boundary conditions.With their distinctive properties compared to conventional criticality, these rich non-Hermitian critical behaviors also inspire further investigations, such as driven dynamics associated with the Kibble-Zurek mechanism [282].

non-Hermitian quantum entanglement
Remarkably, the effect of non-Hermiticity on quantum criticality and entanglement phase transition is profound, as first comprehensively studied in Ref. [283].The subsequent work [284] studied the impact of NHSE on the entanglement dynamics and non-equilibrium phase transitions in open quantum systems.Firstly, they showed the NHSE suppresses the entanglement propagation, leading to a non-equilibrium steady state characterized by the area law of entanglement entropy, in contrast with the volume law for thermal equilibrium states.Secondly, they revealed a new type of entanglement phase transition induced by the NHSE, arising from the competition between coherent coupling and nonreciprocal dissipation; the non-equilibrium steady state exhibits the volume law for small dissipation but the area law for large dissipation, between which the entanglement entropy grows subextensively (i.e., logarithmically with respect to the subsystem size).Anomalously, this non-equilibrium quantum criticality is characterized by a nonunitary conformal field theory whose effective central charge is extremely sensitive to boundary conditions.This originates from an EP in the non-Hermitian Hamiltonian, as also previously put forth in Ref. [285], and the concomitant scale invariance of the skin modes localized according to a power law instead of exponential localization.Moreover, the NHSE leads to the purification and the reduction of von Neumann entropy even in Markovian open quantum systems described by the Lindblad master equation.
The effect of non-Hermiticity on entanglement entropy can be completely understood at the single-particle level for free fermion systems.Through Peschel's formula [286], one expresses the entanglement entropy as S = −Tr[ P log P + (I − P) log(I− P)], where P = µ∈occ.|ψ R µ ⟩⟨ψ L µ | is the single-particle projector onto the occupied bands and P is the truncation of P onto a region demarcated by real-space entanglement cuts.In the Hermitian context, this formula has enabled the explicit identification of topological spectral flow with the flow of P eigenvalues [287-290].In the recent years, similar results have been extended to non-Hermitian models [291, 292], with an enigmatic discovery of negative entanglement entropy at a phase transition within the non-Hermitian SSH model that can be linked to the bc-ghost non-unitary conformal field theory [293].In this work, the entanglement spectra are concretely studied with the PT-symmetry SSH model and the non-Hermitian Chern insulator model, and some attempt was made in linking the effective c = −2 central charge from entanglement entropy scaling with a bc-ghost non-unitary conformal field theory.

Exceptional bound states
It was subsequently found that negative entanglement entropy is a generic feature of lattice models containing EPs [285], because the defectiveness at an EP leads to singularities of the two-point function such that P exhibits special eigenvectors known as "exceptional bound states".The effective central charge c in S ∼ c 3 log N depends linearly on the order of the EP.Interestingly, these exceptional bound states corresponding to P eigenvalues that typically lie way outside the interval [0, 1], and thus contribute strongly to negative entanglement entropy.Due to this spectral gap, they are very robustly protected by the existence of the exceptional band crossing, and constitute a new class of robust bound states distinct from topological and NH skin states.
That said, exceptional bound states exist solely within topological bands, such that they are topologically-protected by construction [294].While Hermitian topologically bands are typically orthogonal to each other even when they cross energetically, being localized at opposite boundaries, they can still experience considerable overlap under the non-Hermitian skin effect.Indeed, in Ref [294], they overlap so well such that they become geometrically defective, and even exhibits super-volume law negative entanglement scaling in the presence of the macroscopic degeneracy of topological flat bands.
While exceptional bound states may seem rather physically elusive, being initially defined as the eigenstates of non-Hermitian Fermi gas propagators, they can also exist as the eigenstates of any physical Hamiltonian or graph Laplacian that are mathematically equivalent to a P operator of another parent system.Ref. [295] reports the first experimental detection of such exceptional bound states in a classical electrical circuit setup.
In the dynamical skin effect, the wave packet acceleration and inelastic scattering are explained by the interplay of the NHSE and the Hermitian wave packet spreading (Fig. 16a) [296].Fundamentally, this acceleration is transiently induced by the non-reciprocal hoppings in the lattice as the wavepacket traverses through the non-Hermitian lattice.Intriguingly, a localized stationary wavepacket can be accelerated by the inherent non-Hermiticity of the lattice and reach the boundary without being reflected (Fig. 16b).The selfacceleration of the wavepacket, in the early time dynamics of a system that exhibits NHSE, is further studied by Ref. [297] and was shown to be proportional to the area enclosed by the energy spectrum of the Bloch Hamiltonian under periodic boundary conditions (Fig. 16c).Intriguingly, non-Hermitian skin modes in semi-infinite lattices can self-reconstruct their shape after being scattered off by a space-time potential, via a phenomenon dubbed as 'self-healing' [61], as illustrated by the schematic Fig. 16d.The work further proves that in a non-Hermitian semi-infinite lattice with a left boundary, any topological edge skin mode at energy E with winding W(E) < 0 and Im(E) larger than the largest imaginary part of the OBC energies, is a self-healing wavefunction (as illustrated in Fig. 16e).
Another novel dynamic phenomenon associated with the NHSE is the non-Hermitian edge burst, which arises from the interplay between NHSE and imaginary (dissipative) gap closure.This is first demonstrated on a lossy lattice with a Bloch Hamiltonian with loss occurring only on one sublattice.Unlike typical quantum-walk models, this features the NHSE.This is true since Eqn.48 is simply the non-Hermitian SSH model with left-right asymmetric hopping.The model can also effectively arise from an open system governed by the quantum master equation: where the effective non-Hermitian Hamiltonian is x , where h is the Hermitian part of Eqn.48, while L x = 2γc B x is the dissipator.The edge burst always occur whenever the imaginary (dissipative) gap closes, i.e. |t 1 | ≤ |t 2 |, except at t 1 = 0.This correspondence has been previously worked out in Refs.[34,226,236]: the existence of NHSE is related to the non-zero area enclosed by the complex energy spectrum.All in all, the edge burst manifests as a substantial loss on the boundary.More fundamentally, its origin can be identified as a universal bulk-edge scaling relation derived via Green's functions [62], with algebraic decay stemming from the closure of the imaginary gap, and the small decay exponent at the boundary stemming from the NHSE.
The above example highlights the remarkable fact that the intimate relationship between the bulk and the edge continues to hold in the state dynamics.Concretely, it was shown that the Lyapunov exponent in the long-time behavior of bulk wave dynamics (far from the edges) can generally reveal non-Bloch symmetry-breaking phase transitions and the existence of the non-Hermitian skin effect [303].
Apart from topological toy models, there are surprises associated with NHSE found in open quantum systems.In [304], the Lindblad master equation was exactly solved for a dissipative topological SSH chains of fermions.The sensitivity on the boundary conditions is reflected in the rapidities governing the time evolution of the density matrix giving rise to a Liouvillian skin effect, which leads to several intriguing phenomena including boundary sensitive damping behavior, steady state currents in finite periodic systems, and diverging relaxation times in the limit of large systems.In a system with quantum jumps and stochasticity, both the short-and longtime relaxation dynamics provide a hidden signature of the skin effect found in the semiclassical limit [305].Even more remarkably, the directed funneling of light at an interface was shown to be possible purely from stochastic fluctuations, even though the hoppings are reciprocal on average [306].A A comprehensive understanding of the effects of NHSE on state dynamics will allow us to better design sensors and devices.For example, it was shown that trapping light at a topological interface with NHSE depends significantly on the initial state [307].Moreover, the implications of non-Bloch band theory to particle Bloch dynamics are profound and lead to new physics.For instance, at the collapse of non-Bloch bands, electrons irreversibly tunnel between Bloch bands in a chiral fashion, contrary to Hermitian systems where Zener tunnelling is oscillatory [269].Tangentially, the Hartman effect -the independence of the phase tunneling time on the barrier width, can exist without any PT symmetry requirements whenever the barrier itself exhibits NHSE [308].Finally, the NHSE can be selectively turned on and off under the presence of static or time-dependent electric fields; this is a consequence of the interplay between Stark localization and dynamic localization, and the NHSE [301].
The presence of the NHSE renders the study of quench dynamics of these non-Hermitian topological models challenging.NHSE dictates the collective localization of states on the boundary under OBCs, therefore coupling the dynamics of different momentum states.Furthermore, the eigenenergy spectra of both the initial and final Hamiltonian (related by the quantum quench) will have generically complex energy, particularly loops in the complex plane.Ref. [298] circumvents this issue by projecting the quench dynamics onto the generalized momentum sectors of the GBZ, in turn revealing the dynamic skyrmions in the generalized momentum-time domain, which are intimately related to the non-Bloch topological invariants of the pre-and post-quench Hamiltonians.This formalism would facilitate the direct detection of non-Bloch topological invariants in experiments.

A. Real spectra and asymptotic non-divergent states from the NHSE
The directed amplification from unbalanced couplings causes a generic initial state to evolve and spread out such that it is amplified more in one direction than the other.As such, under OBCs, it would eventually encounter a boundary and be unable to propagate further.Since directed propagation and amplification are tied to each other in such a NHSE lattice, the state's amplification would also be significantly suppressed by the boundary or even completedly stalled.In the latter case, this would correspond to an energy spectrum that is entirely real due to the NHSE.The reality of the eigenspectra thus has significant implications on the dynamical behaviour, i.e. the eigenstates do not blow up after long time-evolution (but also see [62]).
Non-Hermitian models with real eigenenergies are highly sought-after for their stability.There are many proposals and methods to engineer such systems.The most common way to guarantee real spectra is to enforce parity-time (PT) symmetry on the Hamiltonian, such that the gains and losses conspire to lead to eigenstates with conserved total amplitude [1, 10,11,309].Yet, having a PT-symmetric Hamiltonian is neither a necessary nor a sufficient condition because the PT symmetry itself has to remain unbroken [15].Alternatively, one could obtain real spectra from pseudo-Hermitian systems, which exist only in the a priori unknown engineered real spectrum system [15].There are also specific attempts at engineering real spectra by tuning model parameters, to induce real-complex transitions [310, 311], which will not be the focus of this review.We will like to highlight more general approaches to engineering real spectra via the NHSE itself, without necessarily invoking symmetries of the lattice couplings or their corresponding momentum-space Hamiltonians.
In Ref.
[312], it was shown that real OBC eigenenergies correspond to intersections of the inverse skin depth κ curves at purely real energies (called κ crossings), which can exist even if the PBC eigenenergies are already complex, i.e. if the κ curves cross κ = 0 at Im [E] 0. Hence, as long as the symmetry of the κ curves themselves are preserved, the OBC spectra will be real.Through a systematic investigation, the work [312] gives examples of simple ansatz models with purely real OBC spectrum, that are also as local as possible.More generally, via an electrostatics approach [93] (as outlined in Section IV), one can reverse engineer a parent Hamiltonians for any desired real OBC spectrum and skin localization.All in all, these works open the door to a plethora of models with stable eigenenergies, beyond the use of conventional symmetries.

B. NHSE in Floquet systems
Floquet topological phases are extensively studied in timeperiodic Hermitian Hamiltonians [313,314].The richness of Floquet topological phases in non-Hermitian systems was first explored in the work [315].Here, non-Hermiticityinduced Floquet topological phases with unlimited winding numbers with arbitrarily many real zero-and π-quasienergy edge states are engineered.This is achieved by subjecting a one-dimensional ladder-geometry lattice with a piecewise time-periodic quench.Effectively, periodic driving can induce long-range hoppings, thereby giving rise to emergent physics.In one period of the proposed driving, the quench alternates between a non-Hermitian Hamiltonian H 1 (with asymmetric couplings between sublattices in the same unit cell) and a Hermitian Hamiltonian H 2 , given by A systematic study was done later on the interplay of NHSE and periodic driving [302].Here, a new phenomenon dubbed the Floquet non-Hermitian skin effect (FNHSE) was discovered.The non-Hermiticity not only splits each spectral degenerate point of the parent Hermitian Floquet system into two EPs, but also induce many other EPs.Moreover, the quenched Hamiltonians comprising the periodic quench protocol do not need to exhibit NHSE in order for the periodically quenched system to exhibit FNHSE.The FNHSE does break the BBC but under certain parameter regimes where low-order truncation of the characteristic polynomial can be done with negligible error.Finally, the existence of two different types of Floquet edge modes can still be predicted exactly by introducing the generalized Brillouin zone (GBZ) in two time-symmetric frames.As a remark, the unsupervised identification of Floquet topological phase boundaries has been successful for Hermitian models, and holds the promise to be extended to non-Hermitian models [316].
The study of non-Hermitian Floquet systems was extended to disordered systems possessing non-Hermitian Floquet topological Anderson insulator phases [317], as well as, with spatial modulated on-site potential [318].The NHSE, in the presence of the NHSE, also results in exponentially enhanced Rabi frequencies due to the exponentially large amplified states [299].
Typically, the characterization of topological phases of the time-periodic open quantum systems is done via the use of a dynamical winding number [319,320] or the frequency-space Floquet Hamiltonian.For the latter, the eigen-energy forms the celebrated Wannier-Stark ladder, and at each frequency space lattice site, repeated Floquet bands are observed [321].Floquet topological phase transitions are thus identified by the collapsing of repeated Floquet band gap.Yet, this is not observed in non-Hermitian systems.In the work [322], the authors proposed the non-Floquet theory, which features a temporally non-unitary transformation on the Floquet state, to restore the Wannier-Stark ladders in the Floquet spectrum.Another work proposed a dual topology characterization scheme, so as to circumvent the need to construct the GBZ [323].
Floquet effects in systems with NHSE can be engineered for applications.In a lattice of coupled ring resonators [82], for example, by unconventionally fixing the on-site gain or loss in each ring and by identifying the lattice modes as Floquet eigenstates, the anomalous Floquet NHSE can be realized.Here, skin modes exist at every Floquet quasienergy, allowing for broadband asymmetric transmission, akin to anomalous Floquet insulators [324].In another work [204], it was shown that the presence of NHSE can transform Weyl semimetals into Weyl-exceptional-ring semimetals, as well as a zoo of exotic non-Hermitian topological phases.Finally, it was also proposed that Floquet driving a non-Hermitian topological superconductor can yield multiple Majorana edge modes, useful for realizing environmentally robust Floquet topological quantum computations [325].
Periodic driving can also generate nontrivial "mixed" higher-order topology [326].In Ref. [327], Floquet driving in anomalous Floquet topological insulators are show to generate intrinsic topologically non-trivial non-Hermitian boundary states which are furthermore scale-invariant.The entire system is thus characterized by an unprecedented 'mixed' higherorder topology, where a bulk system with Floquet topology induces a non-Hermitian topology on the boundary.

VII. NHSE BEYOND LINEAR NON-INTERACTING CRYSTAL LATTICES
Typically, the NHSE is formulated on a translationinvariant lattice described by a linear non-interacting tightbinding model.But in fact, the NHSE requires just nonreciprocal breaking of Hermiticity.More interesting incarnations abound when we relax these conditions of translation invariance, linearity, single-particle physics, etc, as we will review below.

A. Breaking of translation invariance -NHSE interplaying with disorder or impurities
The prototypical way to break translational invariance is to introduce disorder, which acts as partial "boundaries" that acquire non-local and non-perturbative influences due to the NHSE [275].Surprisingly, the overall decay scaling is independent of the system's size, in spite of the NHSE which exponentially localizes steady states [65,119].
To understand this, consider a local impurity represented as a modified coupling between the first and last sites, i.e.
where µ ± = µe ±α , with µ controlling the local impurity, α > 0, and x = 0, 1, . . ., L being the lattice site index.This goes beyond the typical interpolation between PBCs and OBCs, with µ ∈ [0, ∞) instead of just µ ∈ [0, 1].When one analytically solves the eigenstates for Eqn.52, we yield two isolated strongly localizing eigenstates at both ends, with the other eigenstates exponentially decaying with a common decay constant When µ = e 2α , i.e. quasi-PBC, it is possible to 'gauge' away to recover the original Hamiltonian under PBC.For µ > e 2α , the leftwards hopping from x = L to x = 0 is further enhanced whereas the opposite is further suppressed.This accumulation is |ψ L,n |/ψ 1,n | = e 2α/µ , i.e. scale-free decay profile.The reverse occurs for 0 < µ < e 2α .
Going beyond hopping with one length scale, we can consider the following which yields a three-fold symmetric spectrum for all values of µ.The impurity-free version of such models have been studied in Refs.[52,135,262,312].While the inverse decay lengths κ L (µ) of SFA states, as induced by the impurity, are insensitive to the exact configuration of non-reciprocal hoppings in the bulk, the OBC skin modes have k-dependent inverse decay lengths.The two phenomena compete and may coexist for some critical k-values, as shown in Fig. 17a, leading to a range of qualitatively different edge-localized regimes beyond usual NHSE states, and also dualities between strong and weak boundary couplings.Beyond this simple impurity toy model, the seminal phenomenon of Anderson localization is modified in the presence of non-Hermiticity, giving rise to the coexistence of localized and extended states even in one-and two-dimensional lattices [118,[332][333][334].More profoundly, chiral hinge states of higher-order non-Hermitian topological insulators remain robust once disorder is switched on, and transitions to surface states as the disorder strength is increased [335].Moreover, disorder can drive the system into a HOTI phase [209].One could also construct a real-space topological invariant for strongly-disordered non-Hermitian systems and in turn, predicting the non-Hermitian Anderson skin effect where the skin effect solely arises from the presence of disorder [121].The evolution of the mobility edge (the energy band boundary between localized states and extended states) and the competition between NHSE and localization effects is exemplified concretely using the Hatano-Nelson model with unidirectional hopping under on-site potential uncorrelated disorder [336].To fully encapsulate the effects of (particularly strong) on-site disorder in an analytical fashion, a "modified GBZ theory" has been suggested -essentially involving the search of the minimum of the polynomial [133].This yields an interval instead of a single point, in turn desirably restoring the bulk-boundary correspondence for disordered samples.The additional restrictions presented by the "modified GBZ theory" also correctly describes the interplay of NHSE and the magnetic field, where the latter similarly breaks translational invariance [133].Other examples of disordered systems with NHSE that have been recently studied include the many-body coupled Hatano-Nelson chains in the presence of a random disorder potential [311], and a quasiperiodic lattice (the disorder is emulated by the incommensurate quasi-periodic on-site potential) [337] with Rashba Spin-Orbit interaction [338].
Disorder may also be used to adroitly realize an effective semi-infinite 1D lattice system, with a complex eigenspectrum that completely fills up the interior of a PBC loop.This is done by concatenating Hatano-Nelson chain segments with random couplings.The result is equivalent to an ensemble of Hatano-Nelson chains with different inverse skin lengths κ [225].
The interplay of the point gap topology (responsible for giving closed PBC loops) and topological defects give rise to more interesting NHSE phenomena.In Ref. [70], two dislocations were introduced into a two-dimensional weak Hatano-Nelson lattice.Concurrently, a skin and anti-skin effect is realized on each defect -a macroscopic localization of states towards one dislocation and a concomitant depletion of states away from the other.A topological invariant for dislocation vector, with state accumulation from the dislocation.A topological invariant is identified, which takes the form of a Z 2 Hopf index that depends on the Burgers vector characterizing the dislocations.Crucially, the anti-skin effect uncovers an additional knob for tailoring the positions of eigenstates in non-Hermitian systems, which is pertinent for applications.
In higher dimensions, mismatch between the macroscopic symmetry and the lattice symmetry itself could also be a means to realize NHSE, even on a reciprocal system [339].This skin effect is solely dependent on the geometry of the system, hence facilitating new routes for wave structuring.
Random impurities can also manifest as fluctuating hopping amplitudes in a spatially-ordered lattice.The NHSE can still be realized in a stochastic system, dubbed the stochastic NHSE, even if the couplings are symmetric on the average [306].The stochastic skin effect stems from the point-gap topology of the Lyapunov exponents under PBC.

B. NHSE without a real-space lattice
Although almost always formulated as arising from asymmetric lattice couplings, the NHSE is just a result of the interplay between non-reciprocity and non-Hermitian gain/loss, and does not require a lattice.Indeed, Ref. [328] suggested that the appearance of the NHSE by an imaginary gauge field β in a finite crystal with OBC is a very general feature, that holds beyond the usual tight-binding models.As discussed, β has the interpretation of an inverse decay length and is responsible for the localization of wavefunctions under OBCs.
Staying with PBC, one can introduce non-Hermiticity via an imaginary vector potential β in the Schrödinger's equation.This is equivalent to complexifying the Bloch wave number, resulting in the energy spectrum to transition from purely real intervals to arbitrarily closed loops in the complex plane.As β increases, the individual closed curves will increase in area and eventually merge with adjacent curves, leading to an open curve in the complex energy plane, approaching the freeparticle dispersion curve in the large β limit.This is illustrated in Fig. 17b.
The imaginary vector potential is crucial in formulating an intriguing duality between non-Hermiticity and curved spacetime [340][341][342].Specifically, by mapping the continuum limit of non-Hermitian lattice models to the Schrödinger equation on a Poincaré half-plane, the inverse localization length κ manifests as an imaginary vector potential which curves the space.The significance of this result is profound.Theorists can study curved spaces on easily accessible experimental non-Hermitian systems such as electrical circuits.Correspondingly, experimentalists can employ readily accessible curved spaces, such as hyperbolic surfaces [343][344][345], to realize experimentally challenging non-Hermitian models with more easily implementable non-Hermitian building blocks.

C. Nonlinear NHSE systems
The interplay with the NHSE and the classical non-linearity leads to intriguing new phenomena, such as trapping effects [346,347], breathers, and solitons [329].
We describe the topological end breather -a novel oscillatory soliton in a nonlinear, non-reciprocal, non-Hermitian lattice that exhibits the NHSE [329].The end breather is strongly localized to a self-induced topological domain near the end of the lattice, in sharp contrast to the extended topological solitons in linear lattices [348].Fundamentally, this is aided by the NHSE which suppresses topologically trivial bulk states, leading to a domain wall near the lattice edge.To understand this, consider the non-reciprocal SSH-like model with Hamiltonian ) where the intercell hoppings are reciprocal, while the intracell hoppings are generally non-linear, of the form where I s is a saturation intensity state.An example of a topological end breather (Fig. 17c) is demonstrated with parameters γ 0 = 0 and γ s = √ 7ν/2.Such non-linear models, as opposed to interacting quantum models, can be realized in classical platforms such as electrical circuits equipped with non-linear elements [349][350][351][352] i.e. diodes or non-linear capacitors [349,[353][354][355], with linear non-reciprocity provided by operational amplifiers [356].

D. NHSE in interacting many-body systems
A plethora of condensed matter phenomena hinges on the many-body nature of quantum interactions giving rise to novel emergent phenomena with no single-particle analog.The NHSE is no exception, where new phenomena can manifest, such as an emergent real-space Fermi surface [357], clustering of the eigenspectrum [358], multifractality of many-body skin effect [359] and many-body distributions caused by NHSE interplaying with Fermionic repulsion [311,357,360], including [361] which was observed on a digital quantum processor.Fundamentally, the competition between repulsion and boundary state localization reshapes (Fig. 17d from Ref. [330]) how the NHSE manifests in non-Hermitian systems, as manifested in the full-counting statistics [362], non-Hermitian Laughlin states [363],Kitaev-Hubbard bosons [364], the Lieb-Liniger Bose gas with imaginary vector potential [365] or spin chain excitations [366].Interestingly, the NHSE itself can also be suppressed by certain correlation/interaction effects [362,367].
A special class of interacting few-body quantum systems can be exactly mapped onto one-body problems in a higherdimensional lattice.Specifically, Refs.[330,331] proposed a framework for mapping an N-body system in d dimensions to a single-body problem in Nd dimensions.In particular, it is possible to achieve robust cluster states where particles are localized next to each other, in the absence of boundaries, solely using particle statistics and appropriately engineered interactions (Fig. 17e from Ref. [330]).By identifying the original interacting chain with the singleparticle Hilbert space of a higher-dimensional configuration lattice, one can consequently observe skin states aggregating at the effective "boundaries" in the many-body configuration space [330,331,368].Fock space skin localization can even be engineered to enhance quantum "scarring" [369], where the ergodicity breaking is strengthened by the accumulation of states in the direction of decreasing Hamming distance from the initial Neel state.This Fock skin effect, however, is complicated by the macroscopically large Hilbert space of states at the same Hamming distance, and do not possess a simple correspondence with the asymmetry in the physical hoppings.Beyond studying the state dynamics and phase diagrams of particular archetypal models, Ref. [370] generalized the topological invariants for many-body non-Hermitian systems.
The NSHE can also emerge in the effective descriptions of various interacting models, even though the interactions themselves are not asymmetric hoppings.For instance, realspace dynamical mean-field theory reveals the NHSE in the pseudospectrum of some strongly-correlated systems [371].NHSE can also occur in the synthetic field moments space of zero-dimensional bosonic quantum dimers [372].An new, interacting form of the NHSE also occurs in the presence of interacting impurities, as manifested by so-called squeezed polarons which are impurity-localized dipole-like density profiles that are impervious to the lattice boundaries [373].More interesting multi-polar NHSE signatures have been suggested in [374], and it remains to be seen if such geometric anisotropy may be generalized to give rise to quantum Halllike states [375] with NHSE interplay.Interestingly, quasi-particle excitations of a closed manybody Hermitian system may exhibit effective dissipation due to their scattering off other degrees of freedom within the considered system, leading to experimentally measurable responses for solid-state systems [376].Finally, coupling an interacting Kitaev honeycomb spin model with the environment leads to an emergent non-equilibrium phase called the exceptional spin liquid [377].
Beyond many-body Hamiltonians, the NHSE -extensive exponential localization of "skin" modes -also arises in the dynamics of open systems.For instance, [378] found an unexpected sudden transition in the purity relaxation rate in the many-body unitary dynamics of qudits, a behavior which can be traced to the asymmetric matrix elements of the Toeplitz matrix underlying the purity dynamics.The NHSE was also revealed in the anomalous behavior of quantum emitters in non-Hermitian heat baths [379][380][381], and interpreted as a Maxwell pressure demon in the many-body context [382].Interesting, under a space-time duality mapping, the sensitivity to initial conditions in a quantum chaotic system can reinterpreted as the sensitivity to boundary conditions in a NHSE system [383].

A. Recent experimental demonstrations
Compared to merely realizing gain/loss, experimentally realizing the NHSE is more challenging, requiring the simultaneous presence of non-reciprocity and non-Hermitian gain/loss.It is not only until 2020 that the NHSE was first demonstrated in electrical [75], quantum optics [14] and mechanical platforms [88]."Topolectrical" circuit realizations implement the non-reciprocity primarily through operational amplifiers [75,77,356,[384][385][386][387], and their versatility [345,388,389] has enabled the experimental simulation of higher-order states [78,80,309] in two or more dimensions, potentially aided by machine learning techniques [80].Photonic setups [83] and single photon quantum walks have also fruitfully demonstrated the transition between the NHSE and other condensed matter phenomena of interest, such as non-Hermitian quasicrystals [390] and topological Anderson insulators [391].Mechanical [88,139] and acoustic setups [79,81] rely on the intrinsic non-reciprocity of the medium, and have also been highly successful in demonstrating the 1D and higher-order NHSE.In Ref. [392], second-order NHSE has been observed in active matter systems of Janus particles, manifesting as spontaneous particle edge guidance and corner accumulation of self-propelled particles.Although NHSE is fundamentally a single-particle phenomenon, experimental demonstrations in solid state NV-center platforms [393] and ultracold atomic lattices [90] constitute significant steps towards physical investigations of the interplay of the NHSE with many-body effects.

B. Experimental proposals and simulations for the NHSE
While the NHSE has aleady been observed in a select set of experiments, there are many other proposals for future experimental demonstration in a variety of physical platforms.Below, we briefly review some existing proposals, each suited for realizing different NHSE-related phenomena.We would omit an explicit discussion of proposals based on classical circuits, since they are already well discussed in the literature [394,395].

Photonics/Optics proposals
The NHSE can be implemented in mature photonic systems such as lasers.in Ref. [396], it was proposed that the interplay between nonlinear gain saturation and the non-Hermitian skin effect gives a laser with the opposite behavior from usual: multimode lasing occurs at low output powers, but pumping beyond a certain value produces a single lasing mode, with all other candidate modes experiencing negative effective gain, thereby giving rise to the NHSE.
Lasers can also indirectly induce the NHSE, such as in exciton-polaritons condensates excited by circularly polarized laser driving in 1D [397], 2D corners [398], with implications on multistability [399].The strong non-linearity in lasers The NHSE can also exist in waveguides or photonic crystals with the requisite loss and reciprocity breaking [400][401][402].In Ref. [400], it was shown that the NHSE can be induced and controlled by varying the positions of the atoms in a waveguide, with the atoms experiencing long-ranged effective couplings and non-Hermitian loss.Further, the NHSE can also be steerable via the gauge field in a coupled optical ring resonator array [403].The NHSE can also be induced in photonic crystals [404,405] i.e. in chiral photonic crystals with anomalous PT symmetry [406].The anomalous Floquet NHSE [82] as well as antihelical topological edge states [407] have also been proposed for a photonic ring resonator lattice.
Finally, single-photon circuits can also simulate the NHSE.Beyond the pioneering NHSE experiment based on photonic quantum walks [14], there has also been other proposals in 2D [408,409], which are potentially useful for also studying the interplay of photon-photonic interactions in a NSHE background, which boasts of a variety of interesting physics [410].

Quantum circuit proposals
Moving beyond few-photon quantum walks, universal quantum simulators for simulating a wide range of many-body phenomena are a rapidly developing technology.Here, we focus only on works relevant for physically demonstrating the interplay of many-body effects with the NHSE.
Implementing the NHSE requires a mechanism for loss, which is not naturally existing in unitary quantum circuits.Hence any such proposal for the NHSE must involve nonunitary evolution implemented through measurements or postselection. in Ref. [411], monitored quantum circuits emulate the non-Hermitian SSH model.It consists of rapidly alternating unitary evolution and measurement steups.The unitary stage provides the unitary evolution of effectively spinless electrons due to the Hermitian part of the Hamiltonian.The periodic measurements that are stochastically invoked correspond to the non-Hermitian σ y term.Under the Floquet-Magnus expansion, effective non-unitary evolution under the non-Hermitian SSH model can be obtained.
In popular state-of-the-art quantum simulators such as the IBM quantum computer, monitored measurements corresponding to particular Kraus operators may not be straightforwardly implemented.An alternative approach is to implement non-unitary evolution by embedding the non-unitary operator within a larger unitary operator, which is an important technique for realizing imaginary time evolution [412][413][414].It can be shown that this can be done with just one additional qubit, known as the ancilla qubit [412].To demonstrate the NHSE on a lattice, there will be a need to for implement a tight-binding model on the chain of spin qubits; this has been implemented for the Hermitian Heisenberg model [415], various 1D topological lattices [416] and Floquet time crystals [417,418], the 2D Chern lattice [419] and even higherorder topological lattices in up to 4 dimensions [420].More recently, non-unitary evolution on a lattice has been implemented by this post-selection approach on a lattice [421]; note that for each qubit experiencing loss, an ancilla qubit is required.Future quantum computer implementations of the NHSE would likely involve lossy qubits as well as a flux ladder.
Going beyond previous experiments that are restricted to preparing ground states via non-unitary imaginary-time evolution [412][413][414], the NHSE and its many-fermion "Fermi skin" profile are realized on current noisy intermediate-scale quantum processors [361], where the NHSE is observed at the many-body level for the first time.

Mechanical/Acoustic proposals
Mechanical systems allow for intuitive observation of non-Hermitian effects, particularly NHSE dynamics, as well as, the relationship between biorthogonality and the system's physical response.[422], The NHSE may be induced through piezoelectric sensors and actuators with Floquet feedback [423] and through flexural phonon modes [424].Mechanical demonstrations of the NHSE may also result in useful applications, such as optomechanically induced transparency [425].

Ultracold atomic proposals
Compared with most other proposals, ultracold atoms provide a promising platform for demonstrating and investigating intriguing quantum many-body physics.In such systems, non-Hermiticity is usually introduced through atom loss, which may be induced by a resonant beam that couples atoms to an excited state [212,426].Realization of NHSE has been proposed in ultracold atoms for both continuous models [427] and optical lattices [100,428], from the interplay of (pseudo)spindependent atom loss and a synthetic flux induced by periodic driving [100] or spin-orbit couplings [280,427,428].It has also been proposed and reported that NHSE can be realized in the momentum space of a two-component Bose-Einstein condensate of ultracold atoms [90,429].

NHSE in complex networks and active media
NHSE is a phenomenon on generic directed graphs, not just crystal lattices.As such, it can also be manifested in graph networks representing real-world processes, for instance in the non-linear dynamics of a rock-paper-scissors game [430].Such models are based on Lokta-Volterra population evolution models, which have been shown to give rise to interesting unexpected non-reciprocal and topological signatures [431][432][433][434][435].
Since the NHSE ultimately stems from non-reciprocity and not the lattice structure per se, it will also manifest in systems devoid of any crystal structure, i.e. continuous media.In [436] and [437], non-reciprocal effects are theoretically and experimentally investigated in active continuous media where the non-reciprocity arises from static deformations, and conserves linear momentum.In [437], a general framework was presented encompassing three archetypal classes of selforganization out of equilibrium: synchronization, flocking and pattern formation.These systems exhibit collective phenomena not lying at a configuration energy minimum, such as active time-(quasi)crystals, exceptional-point-enforced pattern formation and hysteresis.Realization of NHSE has been proposed with magnetic materials, induced by chiral coupling between dipolar-coupledmagnets [438,439].
C. Further discussion on physical NHSE signatures.
Finally, we note that pairs of oppositely localized NHSE state can be observed in reciprocal systems that are mathematically equivalent to appropriately coupled equal and opposite NHSE chains [76,440], or [102] in the case of 2D hybrid ST system.Note that the breaking of bulk-boundary correspondence in the spectrum does not necessarily imply the same in the impedance of a circuit -for an RLC circuit with very different PBC vs. OBC Laplacian spectra, the two-point impedance between most pairs of points can still be approximately the same, whether under PBCs or OBCs [441].
Recently, an implementation of NHSE and its corresponding spectral winding topology has been proposed in electronic mesoscopic systems, with asymmetric coupling between electrons of the concerned system and a reservoir [442].NHSE engineered in this way can be either charge-or spin-resolved in different setups, which can be probed by different transport measurements.On more general platforms, the various transport signatures of NHSE is thoroughly studied in Ref. [443].
Fundamental constraints on the observability of non-Hermitian effects in passive systems, including the NHSE, are derived in Ref. [444], which also discussed about the prospects for observing symmetry-protected edge states and EP signatures.Some observable signatures are embedded in the density of states, particularly the signatures of drastic mode nonorthogonality, which can be effectively exploited and detected by active elements in devices.
band topology from the NHSE 3 A. NHSE and non-Bloch band topology: minimal non-Hermitian SSH model 3 B. Non-Bloch band theory for general cases 5 C. Transition of topological localization direction and half quantized winding numbers 5 D. Topological characterization without the GBZ 6 III.NHSE in higher dimensions 7 A. Chern lattices with NHSE 7 B. Non-Hermitian Weyl and exceptional metals 8 C. Higher-order NHSE and hybrid skin-topological systems 10 1. Hybrid skin-topological modes -Original proposal in a non-reciprocal lattice 10 2. Hybrid skin-topological effect in a honeycomb lattice without asymmetric couplings on the complex energy plane 13 A. Point-gap and line-gap topology 13 B. Topological origin of NHSE 14 C. Quantized response from spectral winding topology 15 D. Complex band evolution as braiding processes 16 E. Emergent spectral graph topology 17 I. INTRODUCTION [see Fig.1(b)].These eigenstates are different from Bloch waves which are eigenstates under PBC and extended in lattice.OBC and PBC spectra are also qualitatively different, forming lines and loops on the complex energy planes respectively, as shown in Fig.1(c).Thus the gap-closing point of H nonH−SSH (k) (touching of the two complex energy bands) in parameter space may not correspond to that of the open-boundary hamiltonian H OBC , indicating the breakdown of bulk-boundary correspondence (BBC).
FIG. 1.(a) Schematic of the non-Hermitian SSH model.(b) Squared wavefunction amplitude |ψ(x)| 2 plots, showing that the bulk modes are localized near the left boundary.(c) Energy spectra of the non-Hermitian SSH model with open(purple) and periodic(red and blue solid lines denote different bands) boundary conditions.Parameters: t 2 = 2, γ = 1, t 1 = 1; L = 160.Figures reproduced from the model in Ref. [50].

FIG. 2 .
FIG. 2. The Brillouin zone (blue loop, a unit circle) and generalized Brillouin zone(pink loop) of the non-Hermitian SSH model, which lies within the Brillouin zone.In general, the generalized Brillouin zone is usually not circular for more general settings.Parameters are t 2 = 2, γ = 1, t 1 = 1. Figure reproduced from the non-Bloch band theory introduced in Refs.[50, 51].

FIG. 4 .
FIG. 4. Topological surface states of non-Hermitian nodal knot metals.(a) The disagreement between the surface state projections with (yellow) and without non-Hermiticity (blue) for the Hopf link.(b) Analytically solved surface state of the Hopf link with k ⊥ = k 2 .(c) Schematic comparing the tight-binding implementations of the Hermitian and non-Hermitian nodal knots.(d) Similar to (a) but for the trefoil knot.(e) Non-reciprocal similarity transforms can rescale z = e ik⊥ , leading to fluctuations of log |z| = 0 level, identified as the 'sea level'.The skin gap closures correspond to the surface states and manifest as 'tidal boundary', whereas the band intersections with the sea level recover the Hermitian drumhead surface states, with a close analogy to 'shore'.(f) More detailed construction of the tidal phase boundary of the trefoil knot.(g) Evolution of the PBC and OBC spectra with k 1 , revealing intricate relations between the complex energy bands, vorticity and the 'tidal surface states'.(h) The 'tidal region' is demarcated by the half-integer vorticity boundary.(i) A direct relationship between the tidal islands and the Seifert surface of the dual nodal knot metal.Figures reproduced from [123].
FIG. 4. Topological surface states of non-Hermitian nodal knot metals.(a) The disagreement between the surface state projections with (yellow) and without non-Hermiticity (blue) for the Hopf link.(b) Analytically solved surface state of the Hopf link with k ⊥ = k 2 .(c) Schematic comparing the tight-binding implementations of the Hermitian and non-Hermitian nodal knots.(d) Similar to (a) but for the trefoil knot.(e) Non-reciprocal similarity transforms can rescale z = e ik⊥ , leading to fluctuations of log |z| = 0 level, identified as the 'sea level'.The skin gap closures correspond to the surface states and manifest as 'tidal boundary', whereas the band intersections with the sea level recover the Hermitian drumhead surface states, with a close analogy to 'shore'.(f) More detailed construction of the tidal phase boundary of the trefoil knot.(g) Evolution of the PBC and OBC spectra with k 1 , revealing intricate relations between the complex energy bands, vorticity and the 'tidal surface states'.(h) The 'tidal region' is demarcated by the half-integer vorticity boundary.(i) A direct relationship between the tidal islands and the Seifert surface of the dual nodal knot metal.Figures reproduced from [123].

FIG. 5 .
FIG. 5. Hybrid skin-topological effect in a non-Hermitian Benalcazar-Bernevig-Hughes model and its 3D extension.(a) The model contains four sublattices in each unit cell, with non-reciprocal real hopping parameters δ 1,2,3,4 along x and y directions.(b) Emergence of 1D boundary modes under x-OBC/y-PBC (cylinder geometry).(c) Boundary modes accumulate towards the corners under double OBC due to the hybrid skin-topological effect.(d) and (e) Spectra of the system with PBC (brown), x-OBC/y-PBC (gray), and double OBC (black).(f) A 3D lattice from stacks of the 2D model in (a), which are coupled through extra non-reciprocal hoppings along z direction.(g) STT corner modes, (h) SST corner modes, and (i) ST0 hinge modes in the 3D model, originated from different hybridizations of NHSE and topological localization along different directions.The figures are reproduced from Ref. [99].

FIG. 6 .
FIG. 6. Hybrid skin-topological effect in a non-Hermitian Haldane model.(a) The model contains real nearest neighbor hopping t 1 , and complex next-nearest neighbor hopping t 2 e ±ϕ , with arrows indicating the direction of hopping with a positive phase.Non-Hermiticity is introduced by the complex on-site potential m ± iγ on the two sublattices respectively.(b) Zigzag boundary of the system.Long solid and dashed arrows [also in (d) to (f)] represent the chiral edge current for the two sublattices, which are toward and opposite to the localization direction of NHSE respectively.(c) Armchair boundary of the system.Arrows connecting different sites in (b) and (c) represent the same property as that in (a).(d) to (f) Distribution of edge modes in the model with negative, zero, positive γ respectively.The figures are reproduced from Ref. [101].

FIG. 7 .
FIG. 7. Schematics contrasting higher-order NHSE and higher-order Hermitian topological phases for 2D systems with L×L sites.(a) Hermitian first-order topological insulator.O(L) chiral or helical modes appear at the edges.(b) Hermitian second-order topological insulator.O(1) zero modes appear at the corners.(c) First-order non-Hermitian skin effect.O(L 2 ) skin modes appear at the edges.(d) Second-order non-Hermitian skin effect.O(L) skin modes appear at the corners.Figures are adopted from Ref. [103].

FIG. 9 .
FIG. 9. Qualitatively different types of energy gaps in Hermitian and non-Hermitian systems.(a) Two gapped Hermitian bands can be flattened to two points at E = ±1 along the energy axis, with a flattened Hamiltonian H 2 = 1.(b) Complex energy bands with a point gap can be flattened to a unit circle in the complex energy plane, with the system's Hamiltonian becoming a unitary one, H † H = 1.(c) Hermitian and anti-Hermitian flattening of complex energy bands with a line gap.A non-Hermitian Hamiltonian with a real (an imaginary) line gap can be flattened to a Hermitian (an anti-Hermitian) Hamiltonian with H 2 = +1 (H 2 = −1).Figures adopted from Ref. [226].

FIG. 10
FIG. 10.A 3D system with different bulk and surface behaviors under different boundary conditions.Gray color indicate the PBC spectrum, and blue and orange colors correspond to the spectra with boundary conditions specified in the figures.(a) The PBC spectrum has a point gap with a nontrivial 3D winding number +1.(b) When OBC is taken along z direction, surface states cover the point-gapped region with a nontrivial 3D winding number.(c) When OBC is taken along y direction, a sharp change of the spectrum indicates the occurrence of NHSE, and in-gap skin modes (orange) appear.Figures are adopted from Ref. [238].
with G m×m the m × m block of the top-right (bottom-left) corner of G(β) for positive (negative) m.A typical example for a model with w(E r ) = 0, 1, 2 for different values of E r is shown in Fig.11.The quantities of ln |G 2×2 | and ν 2 for a reference energy E r with w(E r ) = 2 as a function of β are displayed in Fig. 11(a) and (b)

FIG. 11 .
FIG. 11.Quantized response of non-Hermitian spectral winding topology.Results are obtained from the model of Eq. (16), with t 1 = 1, t −1 = 0.5, t 2 = 2, and t j = 0 for other values of j.(a) and (b) the two quantities extracted from the Green's function as functions of β, the parameters determining the boundary conditions.(c) Spectra at different values of β, corresponding to the five dashed lines in (b) respectively.Red star indicates the chosen reference energy E r = −0.96+ i for (a) and (b).(d) The quantized response quantity ν = Max[ν 1 , ν 2 ], as spectral winding number ranges from 0 to 2 for the chosen system.Figures reproduced from [109].

FIG. 13 .
FIG. 13.Emergent spectral graph topology.(a) An assortment of OBC spectra displaying a rich variety of spectral graph topology.(b)Eigenvalues of Eq. 39 flow into each other as a real boundary flux Re(ϕ) is threaded over a period, for the model given in Eq. 39[54].The outermost loop describes the PBC loop; successively smaller inner loops represent the spectrum as the boundary coupling is decreased through increasing attenuation factors e −Im(ϕ) .(c) Illustration of PBC-OBC spectral flow of Eq. 39, with PBC bulk eigenvalues (red) flow along the blue magenta curves upon threading κ ∝ Im[ϕ], which eventually converge to the OBC eigenvalues (black)[54].(d) An equivalent representation would be plotting the intersections of the inverse skin depth solution surfaces |κ i |[135].(e) The PBC-OBC spectral flow can be generalized to all non-Hermitian lattices, where the intimate relationship between the complex energies E and the inverse spatial decay lengths κ draw a parallel with the electrostatic potential landscape, with PBC and OBC spectral loci corresponding to grounded conductors (V = 0) and lines of induced charges respectively[93].Figures (a,d), (b,c) and (e,f) are taken from Refs.[135],[54] and[93] respectively.

FIG. 16 .
FIG. 16.Dynamical effects associated with the NHSE.(a-b) Wavepacket spreading and acceleration of an initially stationary localized wavepacket in a non-Hermitian lattice [296].Upon traversing to the boundary, the wavepacket does not get reflected, suggesting an inelastic scattering.(c) More generally, the self-acceleration plotted is proportional to the area enclosed by the PBC spectrum [297].(d) Self-healing of a wavepacket upon striking an arbitrary space-time potential.(e) Topological criterion to self-healing edge skin modes [61].(f) Example of Floquet edge states in a periodically driven non-Hermitian lattice [302].Figures (a,b), (c), (d,e) and (f) are reproduced from Refs.[296], [297], [61] and [302] respectively.

FIG. 17 .
FIG. 17. Unconventional manifestations of the NHSE.(a) Lattices with non-reciprocal impurities can result in the coexistence of scale-free accumulation and NHSE [65].(b) Introducing non-Hermiticity, the imaginary vector potential can result in the NHSE even under PBC [328].(c)With both nonlinearity and nonreciprocal non-Hermiticity, a novel oscillatory soliton, a topological end breather, is formed and is strongly localized to a self-induced topological domain near the end of the lattice[329].(d) A combination of particle statistics and suitably engineered many-body interactions can result in boundaries (sites of disallowed occupation) in the many-body configuration space.Topological and skin states can thus form without physical boundaries[330].(e) An example of a 1D chain with well-designed two-body interactions that result in a chiral propagating state along the diagonal boundary in the two-body configuration space[330].(f) This principle can be generalized to the strongly interacting limit, which results in the formation of localized non-Hermitian skin clusters[331], shaped by the connectivity structure of the many-body Hilbert space instead of the real-space lattice.Figures (a), (b), (c), (d,e) and (f) are reproduced from Refs.[65], [328],[329],[330] and[331] respectively.