Curving the space by non-Hermiticity.

Quantum systems are often classified into Hermitian and non-Hermitian ones. Extraordinary non-Hermitian phenomena, ranging from the non-Hermitian skin effect to the supersensitivity to boundary conditions, have been widely explored. Whereas these intriguing phenomena have been considered peculiar to non-Hermitian systems, we show that they can be naturally explained by a duality between non-Hermitian models in flat spaces and their counterparts, which could be Hermitian, in curved spaces. For instance, prototypical one-dimensional (1D) chains with uniform chiral tunnelings are equivalent to their duals in two-dimensional (2D) hyperbolic spaces with or without magnetic fields, and non-uniform tunnelings could further tailor local curvatures. Such a duality unfolds deep geometric roots of non-Hermitian phenomena, delivers an unprecedented routine connecting Hermitian and non-Hermitian physics, and gives rise to a theoretical perspective reformulating our understandings of curvatures and distance. In practice, it provides experimentalists with a powerful two-fold application, using non-Hermiticity to engineer curvatures or implementing synthetic curved spaces to explore non-Hermitian quantum physics.

Introduction

System-environment couplings lead to a plethora of intriguing non-Hermitian phenomena[1-9], such as non-orthogonal eigenstates, the non-Hermitian skin effect[10-16], real energy spectra in certain parameter regimes[17,18], collapsed energy spectra, and coalesced eigenstates at an exceptional point (EP)[12,19,20], and drastic responses to boundary conditions[21,22]. While these phenomena have been extensively explored in quantum sciences and technologies[2-5,8,23-27], peculiar theoretical tools are often required to study non-Hermitian physics. Though bi-orthogonal vectors and metric operators are introduced to restore orthogonality[11,19,28-30], the underlying physics of these mathematical tools is not clear yet. Moreover, it remains challenging to prove the real energy spectra of certain non-Hermitian systems, as the existence of the symmetry does not guarantee a real energy spectrum and sophisticated mathematical techniques are required[17,18,31].

In this work, we show a duality between non-Hermitian Hamiltonians in flat spaces and their counterparts in curved spaces. On the theoretical side, this duality leads to a geometric framework providing a unified explanation of several non-Hermitian phenomena. For instance, it is the finite curvature that requires an orthonormal condition distinct from that in flat spaces, enforces eigenstates to localize at edges, and gives rise to the supersensitivity to boundary conditions. Dual models in curved spaces could be Hermitian, providing a simple proof of the existence of real energy spectra in certain non-Hermitian systems. Moreover, in sharp contrast to existing schemes of studying curved spaces[32-35], which were built on the conventional wisdom that a flat space needs to be physically distorted to become curved, our results show that non-Hermiticity is a controllable knob for tuning curvatures even when the space appears to be flat, for instance, in lattices with fixed lattice spacing. This duality therefore may reform our understanding of distance and curvatures.

In practice, our duality has a two-fold implication. On the one hand, it establishes non-Hermiticity as a unique tool to simulate intriguing quantum systems in curved spaces. For instance, it offers an approach of using non-Hermitian systems in flat spaces to solve the grand challenge of accessing gravitational responses of quantum Hall states (QHS) in curved spaces[34,36-38]. On the other hand, the duality allows experimentalists to use curved spaces to explore non-Hermitian physics. Whereas a variety of non-Hermitian phenomena have been addressed in experiments, delicate designs of dissipations are often required[2,4,7-9,39-42]. Our results show that curved spaces can serve as an unprecedented means to explore non-Hermitian Hamiltonians without resorting to dissipations.

Results

Hatano-Nelson model and hyperbolic surfaces

Our duality can be demonstrated using the celebrated Hatano-Nelson (HN) model[43], which reads,where n = 0, 1, . . . N − 1 is the lattice index of a one-dimensional (1D) chain, ψ is the eigenstate, E is the corresponding eigenenergy, and t and t are the nearest-neighbor tunneling amplitudes towards the left and the right, respectively. Under the open boundary condition (OBC), ψ0 = ψ = 0, , where characterizes the strength of non-Hermiticity, k = mπ/((N − 1)d), m = 1, 2, . . . N − 2, and d is the lattice constant. The eigenenergy reads .

Similar to Hermitian lattice models, the effective theory of Eq. (1) in the continuum limit describes the motion of a non-relativistic (relativistic) particle at (away from) the band bottom and top, with a quadratic (linear) dispersion relation, as shown in Fig. 1a. At the band bottom, the effective theory is written aswhere is the effective mass, and . Solutions to Eq. (2), , have the same energy, . An eigenstate under OBC is their superposition, with and ψ(y0) = ψ(y) = 0. y0 and y specify the positions of the two edges. At the band top, we have M → − M. Eq. (2) is a dimension reduction of the Schrödinger equation on a Poincaré half-plane,where , , and − κ is the curvature (Supplementary Material). The metric tensor is , . Since k is a good quantum number, Eq. (3) reduces to Eq. (2) when k = 0. A finite k adds an onsite potential to the HN model,where . The dimensionless quantity characterizes the strength of V.

Fig. 1

The duality between the Hatano-Nelson (HN) model and a hyperbolic surface.

a A HN chain and its energy spectrum as a function of k. Near a vanishing (finite) K0, the effective theory in curved space is non-relativistic (relativistic). Eigenstates on the HN chain are localized at the edge, ∣ψ∣2 ∝ ∣γ∣2. b A HN chain is mapped to the shaded strip on the Poincaré half-plane, in which an eigenstate with k = 0 satisfies ∣ψ∣2 ∝ y. This shaded strip on the Poincaré half-plane with PBC in the x-direction is equivalent to a pseudosphere embedded in 3D Euclidean space. c The curvature and the inverse of the effective mass, as functions of t for a fixed t. The unites of κ and M−1 are 1/d2 and 2td2/(ℏ2), respectively. (i–v) show the dual pseudospheres of the HN model at various t > 0. A pseudosphere for t < 0 is the same as that for −t.

To derive the duality between the continuum limit of Eq. (1) at the band bottom and Eq. (2), we define with s = nd, such that the eigenstate of the HN model, ψ, defined on discrete lattice sites is extended to ψ(s) as a function of a continuous variable s. Since ψ, under OBC, includes a part that changes exponentially, i.e., , so does ψ. We thus define ϕ(s) ≡ ψ(s)e− with determining the inverse of the localization length, and ϕ(s) varies slowly with changing s. Then we have . Substituting ψ into Eq. (1) and using the Taylor expansion for ϕ(s), , we obtain . Consequently, ψ(s) satisfies

It describes a nonrelativistic particle subject to an imaginary vector potential, A ~ iq[39,43]. Unlike a real vector potential that amounts to a U(1) gauge field, here, an imaginary vector potential curves the space. Performing a coordinate transformation y/y0 = e2 and applying , , we obtain Eq. (2) up to a constant energy shift . The mapping between these two models is summarized in Table 1, which provides a dictionary translating microscopic parameters between them. For instance, ψ, the wavefunction at the n-th lattice site of the NH model is identical to ψ(y), the wavefunction on the Poincaré half-plane evaluated at . The low-energy limit of the eigenenergy of Eq. (1) is also identical to the eigenenergy of Eq. (2) as shown by Table 1.

Table 1

The mapping between the continuum limit of the HN model near the band bottom under OBC and the Poincaré half-plane.

2D Hyperbolic surface	Curvature − κ	Energy scales ℏ2/(2Md2)	Coordinate y_n=y_0e^nκd	Eigenfunctions ψ(y)∝(y/y_0)^12sinmπln(y/y_0)ln(y_N−1/y_0)	Eigenenergies E_m=κℏ^2m^2π^22Mln(y_N−1/y_0)^2
1D Non-Hermitian chain	Non-Hermiticity −4ln^2(∣γ∣)/d^2	Tunneling strength t_Rt_L	Lattice site n	ψ_n∝γ^nsinmπN−1n	E_m2t_Rt_L=−1+m^2π^22(N−1)^2+Om^4N^4

Away from the band bottom(top), similar calculations can be performed by defining using Taylor expansions of the slowly varying ϕ(s) and the same coordinate transformation y = y0e2. We obtain the effective theory near K0d ≠ 0, ±π,where , , and ± corresponds to the left and right moving waves centered near ± K0, respectively, The previously defined has been used. The eigenstate under OBC includes both the left and right moving waves and is written as with eigenenergy of , which recovers the results of the HN model near a finite K0.

Geometric interpretations of non-Hermitian phenomena

Our duality provides a natural explanation of several peculiar non-Hermitian phenomena. Firstly, the orthonormal condition of effective theories in Eq. (2) and Eq. (6) readswhere the normalization constant can be chosen freely. As a common feature of curved spaces, a finite curvature appears in the above equation. Considering a strip in the domain x0 ≤ x ≤ x0 + L, its width in the x-direction depends on y, . Thus, a wave packet traveling in the y-direction must include an extra factor to guarantee the conservation of particle numbers. In the s-coordinate, Eq. (7) is written as . Discretizing this equation with , and transforming it to the HN model, we obtain,where ∣γ∣−2 is precisely the difference between the left and right eigenvectors, or the metric operator[30]. The mapping to a curved space thus establishes an explicit physical interpretation of orthonormal conditions in non-Hermitian systems.

Secondly, the duality allows us to equate the non-Hermitian skin effect to its counterpart on the Poincaré half-plane we found recently[44]. This can be best visualized using the embedding of a hyperbolic surface in three-dimensional (3D) Euclidean space. We define , x = r0φ, where r0 is an arbitrary constant and η > 0, φ ∈ (−π, π). The embedding can then be written as

This is a parameterization of a pseudosphere with a constant negative curvature and a radius of , which satisfies . As shown in Fig. 1b, a pseudosphere features a funnel shape, since the circumference of the circle with a fixed y(η) changes with changing y(η). As previously explained, a coordinate transformation y = y0e2 maps eigenstates on the hyperbolic surface, , to , which exponentially localizes near the funneling mouth, the smaller end.

Thirdly, the collapsed energy spectrum at EP of the HN models has a natural geometric interpretation. When t = t, the pseudosphere reduces to a cylinder with a vanishing κ. For a given t( > t), κ increases with decreasing t. Increasing the non-Hermiticity thus makes the space more curved, as shown by Fig. 1 redc. Approaching EP, t → 0, κ diverges, and the localization length, , vanishes, forcing all eigenstates to coalesce. As eigenenergies read with divergent M, eigenenergies collapse to zero with a massive degeneracy. Across EP, tt < 0, and the effective mass becomes imaginary, all previous results of positive tt still apply provided that M → ±iM. Particles moving in hyperbolic spaces are thus dissipative, and stationary states no longer exist.

Lastly, similar to the HN model, changing OBC to PBC leads to drastic changes in the curved space. Eigenstates of Eq. (2) and Eq. (6) normalized to become , where ) so that ψ(y0) = ψ(y). Correspondingly, eigenenergies become complex. This can be explicitly shown from the time-dependent Schrödinger equations. For instance, at the band bottom(top), we multiply ψ*(y) to both sides of , subtract from the resultant expression its complex conjugate, and integrate over y from y0 to y. We find that the total particle number satisfies,which signifies the absence of a stationary state and explains complex eigenenergies under PBC. Using , we find . This is distinct from the result for OBC, where such that . Similar calculations can be performed for effective theories away from the band top (bottom), . Straightforward calculations show that , which explains the imaginary part of the eigenenergy, .

Despite that y0 ≠ y, these two edges of a hyperbolic surface can be identified in mathematics, since the solutions under PBC exist, as we previously discussed. In physical systems, such PBC can also be realized. In fact, the boundary condition can be continuously tuned. An onsite energy offset, V≥ 0, in one of the lattice sites of the HN model continuously changes PBC to OBC once V increase from 0 to ∞. We consider a superlattice of a lattice spacing of Nd, whose unit cell is a HN chain, as shown in Fig. 2a. Figure 2b shows eigenenergies as functions of V. Similarly, an external potential can be added to the Poincaré half-plane,where is the lattice site of the superlattice. The y-dependent amplitude of the delta-functions guarantees the scale invariance and the equivalence between each section between Y and Y. With V increasing from zero to infinity, eigenstates evolve from those under PBC to the ones under OBC. sets the energy scale of the potential, such that the larger the non-Hermiticity is, the more sensitive of the system is to the boundary condition.

Fig. 2

Changing boundary conditions.

a A constant potential, V, is added to the first site in each unit cell of the superlattice. The corresponding potential in the curved space depends on the position. b Eigenenergies for the ground state and the first two excited states as functions of V. and N = 12 are used.

Generalizations to long-range and non-uniform tunnelings

Whereas the HN model provides an illuminating example of the duality, applications of our approach to generic non-Hermitian models are straightforward. We considerwhere t and t are tunneling amplitudes from the (n ∓ m)th to nth sites. An eigenstate under OBC in the bulk is written as e, where kd ∈ [0, 2π] and q is real. Unlike the HN model, where is a constant, once beyond the nearest neighbor tunnelings exist, q becomes a function of k and defines the so-called generalized Brillouin zone (BZ) in the complex plane[10,15,27,45,46]. Near any point in the generalized BZ specified by K0d ∈ [0, 2π], we define , where ϕ(s) changes slowly as a function of s, corresponding to small deviations of the momentum in the continuum limit. Similar to discussions about the HN model, the effective theory can be formulated straightforwardly using . The Schrödinger equation satisfied by ψ(s) is written aswhere , and depend on K0, as well as t and t. When only the nearest neighbor tunnelings exist, the above equation recovers Eq. (5) at K0 = 0 and becomes real and reduces to a constant imaginary vector potential that we have discussed in the HN model. In the most generic case, provides a complex vector potential, whose real part curves the space. Using a coordinate transformation , where is the real part of , a hyperbolic surface is thus obtained in the same manner as the HN model. The only difference is that κ now is written as and depends on K0. Such K0-dependent curvature provides a geometric interpretation for the generalized BZ. Explicit calculations for a model including the next-nearest-neighbor interaction are given in Supplementary Materials (Fig. S1).

Whereas uniform chiral tunnelings lead to a hyperbolic surface with a constant curvature, we could also consider non-Hermitian models with non-uniform tunnelings,which gives rise to inhomogeneous local curvatures. For slowly varying t and t, we define and such that and . We introduce a slowly changing function ϕ(s) = eψ(s) with . This is a generalization of the uniform case, where ν(nd) reduces to a linear function of n, i.e., the previously discussed in the HN model. Using the same procedure, we obtain the effective theory of Eq. (14). For instance, the non-relativistic theory is written aswhere , g = g = 0, , and the position-dependent curvature is written as . In these expressions, s is obtained from . The constant κ of a hyperbolic surface is recovered when t and t are constants. Changing t and t then tunes local curvatures. For instance, when , , where Θ(x) is the Heaviside step function, the curvature vanishes everywhere except at a particular location, i.e., κ ~ δ(y − y*), where y* = y0 + n*d.

In addition to one dimension, many non-Hermitian models in higher dimensions can be constructed based on the HN model. For instance, 1D HN chains can be assembled to access higher dimensional curved spaces. Whereas curved spaces in higher dimensions are, in general, more complex than those in two dimensions, inter-chain couplings can be engineered to access different higher dimensional curved spaces (see Fig. S2 of Supplementary Materials).

Non-Hermitian realization of QHS in curved spaces

The duality we established has a wide range of profound applications. For instance, Fig. 3a shows a non-Hermitian realization of QHS in curved spaces. When a particle with a charge − e is subjected to a uniform magnetic field, y2∇2 in Eq. (3) is replaced by , where we have chosen the gauge with the vector potential A = (−B/(κy), 0)[47] such that k is still a good quantum number. Wavefunctions of the lowest Landau level (LLL) are written as,whose eigenenergies, , are independent of k, manifesting the degeneracy of the Landau levels. In the dual non-Hermitian systems, a finite magnetic field corresponds to an extra onsite potential in Eq. (4), V → V + V, where , as shown in Fig. 3b. The dimensionless b characterizing the strength of V relative to V is written aswhere has been used. The magnetic flux density, , is thus determined by the ratio of V to V. For a given B, a finite k shifts the position of the minimum of the total onsite potential V + V, similar to the well-known result of flat spaces where k determines the location of the minimum of the potential in the Landau gauge.

Fig. 3

Non-Hermitian realization of QHS on hyperbolic surface.

a A hyperbolic surface threaded by uniform magnetic fluxes. b An extra onsite energy in HN chain, V, encapsulates the magnetic field.

A complete description of QHS requires its gravitational responses in curved spaces. For instance, the particle density, ρ, depends on the local curvature[36],where ν is the filling factor. For integer QHS, ν = 1 and Eq. (18) for a hyperbolic surface can be straightforwardly proved using Eq. (16) (Supplementary Materials). The counterpart of Eq. (18) in the non-Hermitian lattice iswhere is the particle number at lattice site n in the non-Hermitian system. Mapping the magnetic flux, eB/(2πℏ), to the ratio between onsite potentials, , Eq. (18) and Eq. (19) are equivalent. The dependence of densities of QHS on curvatures is thus readily detectable using this non-Hermitian realization. An alternative scheme is to implement a 2D non-Hermitian lattice model as illustrated in Fig. S3 of Supplementary Materials, which serves as a non-Hermitian generalization of the Harper-Hofstadter Hamiltonian[48].

Discussion

In parallel to accessing curved spaces using non-Hermitian systems, experimentalists could also use curved spaces to study non-Hermitian physics[32,33,35]. In conventional understandings, non-Hermiticity arises when dissipations exist. While dissipations have been engineered in certain apparatuses to deliver desired non-Hermitian Hamiltonians[1,6,7], in other platforms, such engineering might be more difficult and sometimes experimentalists may have to use indirect means such as simulating non-Hermitian quantum walks[2,8,27]. Our results show that many non-Hermitian Hamiltonians are readily accessible using existing curved spaces. For instance, hyperbolic surfaces that have been created in laboratories could be used to realize the HN model directly. In particular, in contrast to current schemes used in the study of non-Hermitian physics, this method does not require engineering losses or gains. It thus provides a conceptually different protocol to access non-Hermiticity without dissipations.

The duality we have found provides insightful perspectives for both the studies of non-Hermitian physics and curved spaces. As the curvature depends on the non-Hermiticity even when the separation between any two points in the system is not physically distorted, our conventional understandings of distance may need to be reformed. We hope that our work will stimulate more interest in studying deep connections between non-Hermitian physics and curved spaces.