## Exceptional points in SSH-like models with hopping amplitude gradient

David S. Simon, 1, 2, ∗ Christopher R. Schwarze, 2, † Abdoulaye Ndao, 2, 3, ‡ and Alexander V. Sergienko 2, 4, §

1 Dept. of Physics and Astronomy, Stonehill College, 320 Washington Street, Easton, MA 02357

2 Dept. of Electrical and Computer Engineering &amp; Photonics Center, Boston University, 8 Saint Mary's St., Boston, MA 02215, USA 3

Department of Electrical and Computer Engineering,

University of California San Diego, La Jolla, CA, USA 590 Commonwealth Ave., Boston, MA 02215, USA

4 Dept. of Physics, Boston University,

The Su-Schrieffer-Heeger (SSH) system is a popular model for exploring topological insulators and topological phases in one dimension. Recent interest in exceptional points has led to re-examination of non-Hermitian generalizations of many physical models, including the SSH model. In such nonHermitian systems, singular points called exceptional points (EPs) appear that are of interest for applications in super-resolution sensing systems and topological lasers. Here, a non-Hermitian and nonPT -symmetric variation of the SSH model is introduced, in which the hopping amplitudes are non-reciprocal and vary monotonically along the chain. It is found that, while the existence of the EPs is due to the nonreciprocal couplings, the number, position, and order of the EPs can all be altered by the addition of the hopping amplitude gradient, adding a new tool for tailoring the spectrum of a non-Hermitian system.

## I. INTRODUCTION

Traditionally, research in quantum mechanics has focused on Hermitian systems, since they come with real eigenvalues that can serve as physical observables. But the discovery [1-3] that the weaker condition of parity-time ( PT ) invariance is sufficient to ensure real eigenvalues has caused researchers to re-examine the physics of non-Hermitian Hamiltonians [4, 5]. Eigenvalue degeneracies in Hermitian systems typically correspond to singularities called diabolic points (DPs). At the diabolic point, two or more of the eigenvalues become equal, while the dimension of the space of eigenvectors remains unchanged. As the system is perturbed away from the DP, the splitting between the energy levels is proportional to the perturbation strength /epsilon1 , for small values of /epsilon1 . A simple example of a DP occurs in the Bohr model of hydrogen. In the absence of any external fields, the 2s and 2p orbitals of the atom are degenerate. However when an electric field is turned on, the new eigenstates become 2s-2p hybrid states, and the energy levels split. The electric field magnitude plays the role of the perturbation parameter, and for small fields the energy splittings increase proportional to /epsilon1 . The /epsilon1 = 0 point is the diabolical point: at this point the energies are degenerate, but the eigenstates remain distinct.

Non-Hermitian systems, having inherently complex spectra, allow for the possibility of an additional type of singularity, called an exceptional point (EP) or non-Hermitian degeneracy [6-11], that does not appear in Hermitian quantum systems. As the EP is approached in parameter space, not only do two or more eigenvalues degenerate, but the eigenvectors coalesce as well, causing a drop in the dimension of the eigenspace. Generally, this is accompanied by a transition from real to complex eigenvalues. Initially, these EPs were viewed as a nuisance, since they form obstacles that limit the range of validity of perturbative expansions. But, more recently, they have come to be of great interest due to the rapid variation in eigenvalue splitting when perturbed away from the EP by a small amount [12-16]. This strong response to very small perturbations makes them useful for high-sensitivity sensors [12, 17-21]. Exceptional point sensors have now been developed in a wide variety of physical platforms. Examples include microcavity particle detectors [22], magnetic field sensors [23], and optical gyroscopes [24, 25], and other photonic sensing systems [4, 16, 17, 26]. In addition, PT symmetry and its breaking have been used to develop single-mode microring lasers [27-29] and to generate optical orbital angular momentum (OAM) states [30, 31].

The number of eigenvectors that coalesce is called the order of the exceptional point [11]. The simplest case is a second-order EP, which has an energy level splitting that scales ∼ /epsilon1 1 / 2 as a function of the perturbation strength /epsilon1 ; the

∗ e-mail: simond@bu.edu

† e-mail: crs2@bu.edu

‡ e-mail: a1ndao@ucsd.edu

§ e-mail: alexserg@bu.edu

low exponent implies that the levels vary faster than any polynomial-order perturbation for small /epsilon1 . More generally, an n th-order EP has scaling that is enhanced scale to give an overall sensitivity of /epsilon1 1 /n .

Although there has been some debate over whether the improved sensitivity leads to real sensing improvement once the similarly-enhanced noise is taken into account [32-39], evidence (based on increased Fisher information arguments) seems to be that the improvement remains significant [40] under some circumstances. In any case, improved understanding of EPs and the ability to produce them in different physical platforms continue to be topics of great interest.

Periodic modulations and their effects on topological properties in both Hermitian [41-45] and non-Hermitian [46] Floquet systems have been well-studied, but the effect of monotonic parameter changes has been less examined. Here it is shown that even such monotonic variations can lead to useful effects related to the topological and singularity structure of the system. This will be demonstrated using a non-Hermitian variation on the one-dimensional SuSchreiffer-Heeger (SSH) system, in which it will be shown that allowing a nonzero gradient can alter the singularity structure, and in particular, it can change number and position of EPs and can cause two lower-order EPs to merge into a higher-order EP.

The linear parameter gradient is a fairly general situation in the sense that if the hoppings vary according to any analytic function of position, they can always be approximated for slow spatial variation by the first few terms of a Taylor expansion. The linear term will be the leading non-trivial contribution, and perturbations to this term will control the overall behavior of the system near singular points. So, the behavior seen here should be generic for slow parameter variations.

In the next section, we give a brief review of the SSH model and some of its non-Hermitian generalizations. Then in Section 3 we introduce the specific model studied in this paper, in which both non-reciprocal couplings and a variation of the couplings are introduced along the chain. Our main result is that the introduction of the coupling gradient can alter the EP structure, causing EPs to collide or bifurcate, as well as changing their locations in parameter space. In particular, the collision of two second order EPs can result in the creation of a fourth order EP. We briefly potential experimental implementations in optical and electronic circuit platforms in Section 4, before discussing some conclusions in Section 5.

## II. THE SSH MODEL AND ITS EXTENSIONS

The standard Hermitian Su-Schrieffer-Heeger (SSH) model [47-49] is widely used as a simple example of a onedimensional topological insulator. It was originally proposed as a model of the polyacetylene molecule [47], but has since been realized in physical platforms ranging from cold atoms and photonic systems to electric circuits [50-55]. SSH models have gained increased importance since lasing has been demonstrated in a one-dimensional lattice of polariton micropillars that realizes the SSH Hamiltonian [56]. Lasing has since been shown to also exist in other non-Hermitian systems [57].

The model discussed here consists of a chain of unit cells, each consisting of two distinct lattice sites. Labelling the unit cells by integer m , the two lattice sites within are labelled A and B . Nearest-neighbor hopping amplitudes between the sites alternate, with v being the amplitude to hop between B and A sites within the same cell, and w the amplitude to hop between A of cell m and the B site of the adjacent cell, m -1 (Figure 1). Hermiticity requires that the hopping amplitudes be reciprocal: the amplitude to hop to the left between a pair of cells should be the complex conjugate of the amplitude to hop to the right between the same pair, so that the SSH Hamiltonian takes the form

$$\begin{array} {$$

In the Hermitian case, it is possible, by redefining the phase of the wavefunctions if necessary, to take the amplitudes to be real, v = v ∗ and w = w ∗ . Assuming real amplitudes, the SSH model Hamiltonian can be put into matrix form,

$$H = \left ( \begin{array} { c c c c c c } 0 & v &$$

where the odd-numbered rows and columns represent the A sites at consecutive m values, while the even-numbered rows and columns represent the corresponding B sites. (Note that the signs, or more generally the phases, of the

FIG. 1: The SSH model: each unit cell consists of two distinct subsites or substates, A and B . v is the amplitude for transition leftward from B to A within the same unit cell, while w describes leftward transitions of A in cell m to B of cell m -1. Hermiticity requires that the transitions in the opposite directions along the same bonds to be v ∗ and w ∗ . In the non-Hermitian case this constraint will no longer be hold.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-e2024550-a090-4726-aa65-e1937af7302a.jpg)

couplings can easily be controlled in optical implementations by changing the distance between lattice sites, by inserting phase shifters, or by altering the dispersion of the material between the sites. For simplicity, we keep the couplings real but not necessarily positive.)

The SSH model has a chiral sublattice symmetry, Γ H Γ † = -H , given by Γ = ˆ P A -ˆ P B , where P A,B are projectors onto the A and B sublattices. Because of the chiral symmetry, all nonzero eigenvectors come in opposite sign pairs. More importantly, it has implications for the topological properties of the system. Up to an overall additive constant, the momentum space Hamiltonian can be written as H k ( ) = d ( k ) · σ , where k is wavenumber and σ j are the Pauli matrices. The chiral symmetry forces the complex vector d ( k ) to lie in the x y -plane. In order for the material to be an insulator, with a nonzero energy gap, the remaining two-dimensional vector must avoid the origin. So as k varies across the Brillouin zone, d k ( ) traces out a closed curve, which circles the origin either once, or not at all. In other words, the system has a winding number of either V = 1 or V = 0 about the origin. The system's properties are unusually stable with respect to perturbations, because transitions from one topological phase to the other cannot occur unless the energy gap closes, which can only happen if the chiral symmetry is broken. It can easily be shown that the V = 0 phase occurs when | v | &lt; w | | , while V = 1 occurs when | v | &gt; w | | .

In the model studied here, we will for simplicity assume that all the parameters in the Hamiltonian are real. The fact that the Hamiltonian will then be real will ensure the existence of time reversal invariance. The chiral and time reversal invariances will together guarantee that the energy spectrum is symmetric under reflections across both the real and imaginary axes. When the parameters are allowed to be complex, new effects could potentially appear, but this is a topic for later investigation.

For periodic boundary conditions or for infinite-length chains, the choice of which parameter to call w and which to call v is arbitrary. But for finite chains with open boundary conditions (the case considered here), v is always the amplitude on the last link at either end of the chain. The topologically-nontrivial case is then the one in which the hoppings across those end links are weak, encouraging isolation of the lattice sites at the ends of the chain.

Several non-Hermitian variations on the SSH model have been studied [58-62]. Ref. [60] looks at non-Hermitian models with longer range hopping (longer than nearest neighbor) and Ref. [61] looks at non-Hermitian, PT -symmetric models with complex boundary potentials. Ref. [62] examined models which are non-reciprocal (left-moving and right-moving hoppings are not equal) but not PT -symmetric; these were shown to share several properties with PT-symmetric models. In the non-Hermitian case, the amplitudes for hopping to the left and to the right are not necessarily complex conjugates of each other,

$$H \ = \ \sum _ { m } \left ( v _ { L } | A, m$$

Here, we look at finite-length SSH-like models with nearest-neighbor hoppings and open boundary conditions at the ends. Two parameters are added to the model: one to control the amount of non-reciprocity in hopping amplitudes, and one to control the rate of variation of hopping amplitudes with position along the chain. We assume the simplest case, with constant hopping gradient, or in other words, an increase or decrease of hopping amplitude that varies linearly along the length of the chain. Exceptional points (EPs) will form continuous families in the parameter space, and for long chains the structure of the spectrum and the distribution of the EPs can become highly complex, with the number and order of the EPs grows as the SSH chain becomes longer. The nonzero gradient breaks the discrete translational symmetry of the system; as has been seen in other systems [14, 63-66], the breaking of spatial symmetries often leads to the appearance of EPs.

with w L = w ∗ R and v L = v ∗ R . Non-Hermiticity commonly arises in open systems, where gain and loss can occur.

/negationslash

/negationslash

There is a rich landscape of variations of related models. Here we look at a simple model that seems to be fairly representative, focusing on short chain lengths for which analytic expressions can be found for the eigenvectors and spectrum. Because of the lack of translational invariance, it is more convenient to work in position space, rather than in momentum space.

FIG. 2: Examples of the spectrum for | v | ≤ | w | . Here the values v = 0 1, . w = 1 were used. Energies are plotted versus non-reciprocity parameter a for several values of gradient parameter g . The solid blue curve is the real part of the energy and the dashed yellow curve is the imaginary part.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-88bc96cf-2cdf-4b0b-8513-c3a849c961c3.jpg)

## III. THE NON-RECIPROCAL LINEAR GRADIENT MODEL

Take a one-dimensional lattice with hoppings between nearest-neighbor sites, and assume a position-space Hamiltonian of the form

$$H = \begin{pmatrix} 0 & v + g & 0 & 0 & 0 & 0$$

/negationslash

For g = 0 the model is non-reciprocal, with different probabilities for hopping to the left and to the right. This model is neither Hermitian nor PT -symmetric for a = 0, but reduces back to the SSH model for a = g = 0. Because of the asymmetry of the hopping coefficients, the amplitude will tend (for open boundary conditions) to collect at one end of the lattice, in accord with the non-Hermitian skin effect [58, 67]. Such boundary states exhibit new effects not seen in Hermitian systems [68-71]. The additional parameter a controls how fast the hopping rates vary along the chain. We will restrict ourselves here to the case where all parameters ( v , w a , , and g ) are real. For a = 0 and v = w , the model examined here reduces to the Hatano-Nelson [72, 73]. Recently, Ref. [74] demonstrated a method of generating non-Hermitian Hamiltonians with real spectra from Hermitian Hamiltonians; this approach can generate both spatially varying hopping coefficients and nonreciprocality. The current paper can be viewed as a case of this more general setup with vanishing on-site potential and a particular choice of hopping coefficient variations. Other non-reciprocal models have been studied on occasion in platforms ranging from photonic and acoustic systems to mechanical systems [75-78].

/negationslash

For now, focus on the simplest nontrivial case, of four lattice sites, spanning two unit cells:

$$H = \begin{pmatrix} 0 & v + g & 0 & 0 \\ v - g$$

Representative examples of the spectrum versus g are shown in Figs. 2 and 3 for several values of a . It can be seen that the spectrum is very different for the cases | v | ≤ | w | and | v | ≥ | w | , which would correspond to the topologically trivial and nontrivial cases for the standard SSH model. The structure varies considerably as the gradient parameter a is varied, with exceptional points appearing, disappearing, colliding, and merging. In addition, there are bifurcations of the real and imaginary parts at points well removed from the E = 0 axis. The goal below is to classify this complicated behavior in a relatively simple manner.

/negationslash

First, consider the nonreciprocal case ( g = 0) with no gradient ( a = 0). Setting a = 0 in Eq. 5 the eigenvalues are readily found to be of the form

$$E = \pm \sqrt { A \pm B },$$

where

2

2

2

$$A \ & = \ - 3 / 2 g ^ { 2 } + v ^ { 2 } + 1 / 2 w ^ { 2 } \\ B \ & = \ \sqrt { A ^ { 2 } - ( g ^ { 2 } - v ^ { 2 } ) }.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-890868e1-bd7f-438f-aa46-7b1a189edb6f.jpg)

1.5

1.5

FIG. 3: Examples of the spectrum versus g for | v | ≥ | w | . The values v = 1 0, . w = 0 1 were used. . The solid blue curve is the real part of the energy and the dashed yellow curve is the imaginary part.

There are two potential sources of EPs, leading to three opposite-sign pairs of EPs. One possibility occurs when the expression inside the square root of Eq. 6 passes through zero. This happens when A 2 = B 2 , or equivalently when

$$g = \pm v.$$

These EPs always occur at zero energy, and will be referred to as type I points.

The second possibility (type II points) is for the expression inside the square root of Eq. 8 to pass through zero; this happens when

$$g = \pm v \sqrt { \frac { 2 } { 5 } ( 2 v ^ { 2 } + \frac { 1 } { 2 } w ^ { 2 } ) },$$

or

$$g = \pm w.$$

The energies at these points are E = ± A which can be either zero or nonzero.

Now allow a = 0. Conceptually, it is expected that the addition of hopping gradients should lead to interesting behavior. This is because, as one moves along the chain, there will be points where the magnitudes of the effective local hopping coefficients will cross each other, with transitions between regions with | v -ma &gt; w | | -ma | and regions with | v -ma &lt; w | | -ma | as the site index m increases. Thus, the gradient would be expected to cause transitions between topological phases in the chain's bulk.

/negationslash

With the gradient turned on, the energies are still of the form of Eq. 6, but with A and B now being given by

$$A \ & = \ \left ( \ - \frac { 3 } { 2 } g ^ { 2 } + v ^ { 2 } + \frac { 1 } { 2 } w ^ { 2 } \right ) + \left ( \frac { 5 } { 2 } a ^ { 2 } - 2 a v - a w \right ) \\ D \ & = \ \overline { ( \lambda 2 \cdots ( \lambda 2 \cdots 2 \vee 2 + \lambda \cdots ( \lambda \cdots ( \lambda 2 \cdots 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \veed 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \veew 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vew 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vees 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \ve. 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \ve� 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \veo 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vea 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \veea 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vei 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \veex 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \ve e 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vee 2 \vema \$$

It is interesting to note that the energies only depend on the geometric means of the Hamiltonian entries on opposite sides of the diagonals, i.e. opposite-direction hopping amplitudes on the same link. To be more specific, define the geometric means

$$B \ = \ \sqrt { A ^ { 2 } - ( g ^ { 2 } - v ^ { 2 } ) ^ { 2 } + 4 a ( a - v ) ( g ^ { 2 } - v ^ { 2 } ) }, \\ \intertext { a n o t h o n o i r a n d o n t h o m o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t h o m o t i o n }$$

$$m _ { 1 } ^ { 2 } \ & = \ ( v + g ) ( v - g ) \\ m _ { 2 } ^ { 2 } \ & = \ ( w + a - a \vee w - a - a \vee$$

$$m _ { 3 } ^ { 2 } \, = \, ( v + g - 2 a ) ( v - g - 2 a ) ;$$

$$m _ { 2 } ^ { \bar { 2 } } \, & = \, ( w + g - a ) ( w - g - a ) \\ m ^ { 2 } \, & = \, ( w + a - 2 a \nmid n - a - 2 a \nmid \cdot }$$

then the constants A and B can be written as

$$A \, & = \, \frac { 1 } { 2 } \left ( m _ { 1 } ^ { 2 } + m _ { 2 } ^ { 2 } + m _ { 3 } ^ { 2 } \right ) \\ D \, & = \, \sqrt { \frac { \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 - \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 + \lambda 2 } }$$

$$B \ = \ \bar { \sqrt { A ^ { 2 } + ( m _ { 3 } ^ { 2 } - 2 m _ { 1 } ^ { 2 } ) m _ { 1 } ^ { 2 } } }$$

With a = 0, the type I EPs split into two pairs of points separated by a distance proportional to a :

/negationslash

$$g = \pm v, \, \pm ( 2 a - v ).$$

,

4

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-fd237ce7-be3a-49f5-ae8a-49faa19c1d83.jpg)

4

4

4

FIG. 4: A and B versus non-reciprocity parameter g for an example with v = 1 0, . w = 0 1. . The solid blue and red curves and the dashed black curve are, respectively, Re B ( ), Im B ( ), and A . As the gradient parameter a varies, EPs are created or removed when the A and B lines touch each other or the horizontal axis. A more detailed description of the behavior can be found in the text.

The type II points become more complicated, but can still be written in analytic form:

$$g \, & = \, \pm \sqrt { Q \pm 2 \sqrt { P } }, \\ P \, & - \epsilon \Omega _ { n } ^ { 2 } \, \mp \, ( n, 2 \, - \, 7 n ^ { 2 } ) ^ { 2 } \, \mp \, ( n, 2 \, \mp \, n ^ { 2 } ) ^ { 2 } }$$

2

$$P \ = \ - 6 9 a ^ { 2 } + ( v ^ { 2 } - 7 a ^ { 2 } ) ^ { 2 } + ( w ^ { 2 } + a ^ { 2 } ) ^ { 2 } \\ + \Delta a ( \Omega a ^ { 2 } \Omega + a ^ { 2 } \Omega - w ^ { 3 } - 2 \Omega \Omega \Omega + w ^ { 2 } \Omega + w ^ { 2 } - w ^ { 3 } )$$

$$\begin{array} {$$

a w

2

v

3

2

avw

+

v w

2

+

vw

2

w

3

)

,

Note the double square root structure of Eq. 20 compared to the single square root of Eq. 10; this is what opens the possibility of fourth order exceptional points occurring when a is nonzero.

As can be seen from Fig. 4, as a is varied the curves of A and B versus g move up and down, while the real part of B can also develop a double-well shape. As a result, the number of intersections of the A and Re B ( ) curves changes. At these transition points, the imaginary part of B can become nonzero, creating an exceptional point.

To study the behavior of the system near the exceptional points more quantitatively, it is necessary to introduce perturbations. So consider perturbations in the gradient, a → a + , with real /epsilon1 /epsilon1 &lt;&lt; a . Under such a perturbation, it is found that

$$E \ & = \ \pm \left \{ A _ { 0 } + ( 5 a - 2$$

$$\begin{array} { r c l } A & = & A _ { 0 } + ( 5$$

$$\sqrt { B _ { 0 } ^ { 2 } + 4 ( 2 a - v ) ($$

$$B \ = \ \sqrt { B _ { 0 } ^ { 2 } + 4 ( 2 a$$

where A and B have been expanded to first-order in epsilon, and zero subscripts indicate the unperturbed expressions. The Type II case above can now be split into two distinct cases, so that there are now four cases total (including the trivial case):

## Case 0: B 0 = 0 , A 2 0 = B 2 0 :

In this case, expanding the square roots in Eqs. 23 and 25, the result will be linear in /epsilon1 , so that no EP occurs.

/negationslash

/negationslash

## Case I: B 0 = 0 , A 2 0 = B 2 0 :

If A 0 vanishes, the constant ( -independent) term inside the square root of Eq. 23 vanishes, leaving a leading term /epsilon1 of order √ /epsilon1 . There is therefore a second-order EP at the parameter values given by Eq. 19.

/negationslash

## Case IIA: B 0 = 0 , A 2 0 = B 2 0 :

Now, using the binomial approximation, the energy is

/negationslash

$$E \ = \ \pm \sqrt { E _ { 0 } ^ { 2 } \pm 4 ( 2 a - v ) ( g ^ { 2 } - v ^ { 2 } ) \sqrt { \epsilon } } + \mathcal { O } ( \epsilon ) \\ r \quad 4 ( 2 a - n ) ( a ^ { 2 } - n ^ { 2 } ) \quad _ { - 1 }$$

$$\approx \stackrel { \cdot } { \pm E _ { 0 } } \left [ 1 \pm \frac { 4 ( 2 a - v ) ( g ^ { 2 } - v ^ { 2 } ) } { 2 E _ { 0 } ^ { 2 } } \sqrt { \epsilon } \right ],$$

giving a second-order EP at the solutions to B 0 = 0.

$$\text{Case IIB} \colon A _ { 0 } ^ { 2 } = B _ { 0 } ^ { 2 } = 0 \colon$$

or

$$a \ & = \ \frac { 1 } { 5 } ( 2 v + w ) \pm \sqrt { ( 2 v + w ) ^ { 2 } - 5 ( w ^ { 2 } - v ^ { 2 } ) }, & ( 3 1 ) \\ g \ & = \ \pm v. & \quad \text{unused in Table 1 Note that when $a-0$ case IIR sequence at the locational transition not}$$

These cases are summarized in Table I. Note that when a = 0, case IIB occurs at the topological transition point of the standard SSH model, | v | = | w | = | g | .

Now consider Case II ( B 0 = 0). One possibility is that g = ± w , in which case the middle link becomes one-way, as in the bottom row of Figure 5. In this latter case, the amplitude eventually becomes localized in the last two links of the chain at one end or the other. Here, both the real and imaginary parts of the energy bifurcate, breaking degeneracies. One part (real or imaginary) is already zero at the bifurcation point, while the other is nonzero. At these points, two of the eigenvectors vanish.

A coherent picture of these cases can now be constructed. Again consider a = 0 first. When B 2 0 = A 2 0 (Case I), we have g = ± v , in which case either the v + g or v -g terms in H vanish. Then the links at the end of the chain become one-way, as in the top row of Figure 5. As a result, the amplitude tends to accumulate at one site at the left or right of the chain. The parameter-space points where these 'traps' or attractors occur turn out to be especially interesting when a becomes nonzero. For a = 0, the real parts of the energies cross the E = 0 axis, but no EP occurs at the crossing points, since the imaginary parts between the two points remain zero. But as a starts to move away from zero, these two points split into four: two points remain fixed at g = ± v , while the other pair, at g = ± ( v -2 a ), moves away from the first pair. As these two pairs separate, an arc of nonzero imaginary part connects v to v -2 a and another such curve connects -v to -( v -2 a ), so that all four points become EPs. Turning on the gradient term therefore creates four new exceptional points at g = ± v, ± ( v -2 a . )

The other possibility for case II at a = 0 is when the nonreciprocity parameter becomes g = √ 2 5 (2 v 2 + 1 2 w 2 ), which

/negationslash then implies that E = ± √ A 0 = √ 1 5 ( w 2 -v 2 ). Assuming that A 0 = 0, these are points at which bifurcations occur at nonzero energies; the real and imaginary parts separately bifurcate, one starting from zero and one starting from a nonzero value. These points are second-order EPs, at which two eigenvectors merge.

/negationslash

/negationslash

TABLE I: Summary of the exceptional points for the length-4 chain, in terms of the non-reciprocity and gradient parameters, g and a .

| Case                         | EPs                                                                            |
|------------------------------|--------------------------------------------------------------------------------|
| (0) B 0 = 0, A 2 0 = B 2 0   | No EP                                                                          |
| (I) B 0 = 0, A 2 0 = B 2 0   | 2nd-order EPs at g = ± v, ± (2 a - v ) ( a arbitrary)                          |
| (IIA) B 0 = 0, A 2 0 = B 2 0 | 2nd-order EPs at solutions to B 0 = 0                                          |
| (IIB) B 0 = 0, A 2 0 = B 2 0 | 4th-order EPs at g = ± (2 a - v ), a = 1 7 [ 4 v - w ± √ 9 v 2 - 8 vw +8 w 2 ] |

Now,

$$E \ = \ \pm \sqrt { ( 5 a - 2 v - w ) \epsilon \pm 4 ( 2 a - v ) ( g ^ { 2 } - v ^ { 2 } ) \sqrt { \epsilon } } \\ \approx \ \pm \sqrt { \pm 1 } \sqrt { 4 ( 2 a - v ) ( g ^ { 2 } - v ^ { 2 } ) } \cdot \sqrt { \epsilon } + \mathcal { O } ( \epsilon ^ { 3 / 4 } ), \\ \text{EPs at the solutions to } B _ { 0 } = 0. \text{ Each fourth-order EP is formed in the collision of two second-order}$$

giving fourth-order EPs at the solutions to B 0 = 0. Each fourth-order EP is formed in the collision of two second-order EPs.

Case IIB is the overlap of cases I and II. So the situation can be readily summarized as follows: if a point in parameter space falls into either category I or II there is a second-order EP at that point, while if it falls into both classes simultaneously (IIB), it becomes fourth-order.

Case (IIB) occurs when the parameter values are given by either

$$a & \, = \, \frac { 1 } { 7 } \left [ 4 v - w \pm \sqrt { 9 v ^ { 2 } - 8 v w + 8 w ^ { 2 } } \right ], \\ g & \, = \, \pm ( v - 2 a ),$$

/negationslash

/negationslash

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-b78f85c1-0352-4b3e-aa17-609004d58f0a.jpg)

| g = + v                                | g = - v                                |
|----------------------------------------|----------------------------------------|
| (v-g = + 0, v+g = 2v) v+2g v+g w+g w-g | (v-g = + 2v, v+g = 0) v-g v-2g w+g w-g |

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-98753f4d-4411-405a-a9a4-bfa178085199.jpg)

/negationslash

FIG. 5: Cases for B 2 0 = A 2 0 = 0 . In the top row, g = ± v , causing one of the end links to become one-way and leaving all of the amplitude funneled to the site at the left or right end of the chain. Similarly, in the bottom row g = ± w , causing the middle link to become one-way. In the latter case, the amplitude still becomes localized at one end or the other, but oscillating between the last two sites.

FIG. 6: Real (a) and imaginary (b) parts of the energy as function of nonreciprocity g and gradient a , for a chain of length N = 4. The hopping parameters were taken to be v = 1 0, . w = 0 1 . . In (a), the curves marked C are examples of lines of second order EPs. The intersections of two of these at points such as the one marked A lead to fourth order EPs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-6bcaa297-1a5c-401f-91e2-38590ec326d1.jpg)

/negationslash

If a = 0, the additional constraint A 0 = 0 cannot be satisfied concurrently with B 0 = 0 (except in the trivial case v = w = 0), but allowing nonzero a opens access to points at which Case IIB occurs. At such points, the E = 0 bifurcations mentioned in the previous paragraph will now occur at E = 0, and two pairs of second-order EPs collide, forming fourth-order exceptional points. These collisions and mergers can be seen clearly, for example, between the second and third images of Figure 2, where pairs of pitchfork-shaped imaginary parts move onto the real axis and coincide with each other. As a is increased further, the merged points separate again; after the separation, they remain on the real axis, forming two second-order EPs.

The full evolution of the eigenvalues can be seen when plotted as a function of both g and a (Figure 6). It can be observed that for N = 4, the real and imaginary parts both consist primarily of three sheets, which split for some parameter ranges into opposite sign subsheets. The bifurcation points form curves of exceptional points in the parameter space; two examples of such curves are marked C in Figure 6(a). The fourth-order EP emerges when two of these curves intersect, for example at the point marked A in the Figure.

The pattern persists for longer length chains, where the number of principle sheets is N -1, with each sheet splitting at exceptional point curves. The N = 6 case is shown as an example in Figure 7. While, for the N = 4 case there are four second-order EPs for a = 0, which then split or merge as a is varied, in the N = 6 case there are eight EPs at a = 0, which then undergo similar evolution as a changes. As the chain length increases, the number of spectral sheets and the number of EPs both increase, with the spectrum rapidly becoming extremely complicated.

## IV. POSSIBLE EXPERIMENTAL IMPLEMENTATIONS

Optical and photonic platforms (see Refs. [79-81], among others) and electric circuits [82] have been implemented that provide physical embodiments of the standard SSH system. Many of these can be altered to implement the

FIG. 7: (a) Real part of the energy as function of nonreciprocity g and gradient a , for a chain of length N = 6. The hopping parameters were taken to be v = 1 0, . w = 0 1 . . (b) An example of a cross-section of this plot for a = 0 35. . The solid blue curves are the real parts and the dashed yellow curves are imaginary parts.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-8914bb87-8118-40c2-aa47-d744a426b730.jpg)

FIG. 8: Optical implementation of non-reciprocal hopping. The circles on the two sides are optical circulators. The upper and lower lines in the left-hand figure contain loss represented by α j and a phase shift φ j . α j is the fraction of amplitude that remains after the loss. The left-hand figure will be used as shorthand for the arrangement on the right.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-f891704b-cae4-490e-8afe-76548d268756.jpg)

model of the previous section, providing a means of easily producing a highly-controllable set of EPs in systems whose behavior can be altered in real time.

For example, this model can be realized with linear optics, augmented by optical circulators. Non-reciprocal hopping between adjacent lattice sites can be arranged using a pair of circulators, a pair of phase shifters, and some loss. The circulators make use of magnetic fields, and therefore break time reversal symmetry. The hopping amplitudes can be implemented in an arrangement of the form shown in Fig. 8, where the circles on the two sides are optical circulators. The circulators ensure that amplitude can only flow to the right in the upper branch and only toward the left in the lower branch. Denoting the upper and lower branches by labels 1 and 2, respectively, the amplitude in each line encounters a phase shifter with phase φ j for j = 1 2. , The losses are represented by real, positive parameter α j . The amplitude for transmission to the right is then α e 1 iφ 1 , with transmission to the left being given by α e 2 iφ 2 . The losses can be implemented by inserting Mach-Zehnder interferometers with one output port used to eject light from the system; by altering the phase in the interferometer the amount of loss through this port can be controlled. The magnitudes of hopping amplitudes can also be easily controlled, again by use of Mach-Zehnder interferometers (or by temperature control if the system is implemented on an integrated chip). So this optical implementation allows all parameters to be tuned in real time via experimenter-controlled phase shifts.

Electric circuits provide another viable platform. A circuit implementation of the standard SSH model [82] is shown in Fig. 9(a). Each unit cell consists of pairs of inductors and capacitors, with the capacitance values alternating from one cell to the next. A hopping gradient could easily be introduced into this system by varying the capacitances from one cell to the next. Loss and nonreciprocality are readily introduced by the addition of resistors and diodes (Fig. 9(b)). The use of variable resistors and capacitors allows the loss and phase shifts to be controlled in real time.

Both optical and electrical implementations can easily be constructed (either on a table top or in integrated form) with low resource cost and allow real time control of all system parameters. This greatly simplifies experimental study of the topological and singularity structures of the system, and allows band gaps to be opened or closed at will. Different experiments may be run on the same apparatus, with only a change of a few resistance and capacitance values, or the change of a few phase shifts in the optical case. In either case, this is achieved with the turn of a knob, rather than a physical alteration of the apparatus.

It is possible that a nanopillar array (Fig. 10) such as that used to implement SSH-based topological lasers in Ref. [56], can also be arranged to implement the desired gradient, simply by gradually increasing the spacing between

FIG. 9: (a) An electric circuit implementation of the standard Hermitian SSH model. (b) Electrical analog of the module in Fig. 8. The phase shifts can be implemented by appropriate inductor-capacitor combinations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-d150dfc4-84f5-4178-952a-e4015f073ad0.jpg)

FIG. 10: Lasing occurs in staggered nanopillar SSH arrays [56]. By gradually increasing the spacing between pillars from left to right, a gradient in the hopping amplitudes can be introduced. Heating of the substrate can lead to a hopping amplitude perturbation. The shaded lobes represent the amplitudes of p -state polariton excitations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00879v1-figure-5d742e53-bcd7-4f43-99f8-7a51b1d9663a.jpg)

adjacent staggered pillars. Coupled microring cavities could be used in a similar manner. Introducing non-reciprocal transitions is more difficult, but may be possible using recent advances in integrated photonics, such as the integrated isolator developed in Ref. [83].

## V. CONCLUSIONS

In the previous sections, a non-Hermitian model has been presented and analyzed which generalizes the SSH model by introducing both nonreciprocal couplings and spatial variation of the coupling constants. The nonreciprocity creates exceptional points, while the introduction of the hopping gradients affects the number and nature of the EPs. By implementing this system in a physical platform that allows the hopping gradients to be controlled, this introduces a means of tailoring the EPs to have desired properties or be located at desired points in parameter space. In particular, if the gradients can be varied in a user-controllable manner, it allows experimental study of new situations such as the braiding of exceptional points around each other in parameter space.

The model studied here is the simplest in a class of models that has a rich and unexplored expanse of variations. By allowing the hopping coefficients to vary in a nonlinear manner, or by allowing the nonreciprocity parameter g to also be spatially-varying, additional control over the spectrum can be introduced, with the potential of uncovering further interesting and useful effects.

## Funding

This research was supported by the Air Force Office of Scientific Research MURI Award No. FA9550-22-1-0312, DOD/ARL Grant No. W911NF2020127, and the Beckman Young Investigator Award, from the Arnold and Mabel Beckman Foundation.

Disclosures

The authors declare no conflicts of interest.

## Data Availability Statement

All data generated for this paper is available upon reasonable request.

## References

- [1] C. M. Bender, S. Boettcher, 'Real Spectra in Non-Hermitian Hamiltonians Having PT Symmetry', Phys. Rev. Lett. 80 , 5243 (1998)
- [2] C. M. Bender, D. C. Brody, and H. F. Jones, 'Complex Extension of Quantum Mechanics,' Phys. Rev. Lett. 89 , 270401 (2002)
- [3] C. M. Bender, PT Symmetry: In Quantum And Classical Physics (World Scientific 2018)
- [4] S. K. Ozdemir, S. Rotter, F. Nori, L. Yang, 'Parity-time symmetry and exceptional points in photonics,' Nature Mat. 18 , 783 (2019)
- [5] L. Feng, R. El-Ganainy, L. Ge, 'Non-Hermitian photonics based on parity-time symmetry,' Nature Phot. 11 , 752 (2017)
- [6] T. Kato, Perturbation Theory of Linear Operators (Springer, Berlin 1966)
- [7] M. V. Berry, 'Pancharatnam, virtuoso of the Poincar´ sphere,' Curr. Sci. e 67 , 220-231 (1994)
- [8] M. V. Berry, D. H. J. O'Dell, 'Diffraction by volume gratings with imaginary potentials,' J. Phys. A: Math. Gen. 31 , 2093-2101 (1998)
- [9] W. D. Heiss, 'Phases of wave functions and level repulsion,' Eur. Phys. J. D 7 , 1-4 (1999)
- [10] W. D. Heiss, 'Repulsion of resonance states and exceptional points,' Phys. Rev. E 61 , 929 (2000)
- [11] W. D. Heiss, 'The physics of exceptional points,' J. Phys. A: Math. Theor. 4, 444016 (2012)
- [12] H. Hodaei, A. Hassan, S. Wittek, et al., 'Enhanced sensitivity at higher-order exceptional points,' Nature 548 , 187 (2017)
- [13] M. P. Hokmabadi, A. Schumer, D. N. Christodoulides, and M. Khajavikhan, 'Non-Hermitian ring laser gyroscopes with enhanced Sagnac sensitivity,' Nature 576 , 70 (2019)
- [14] J.-H. Park, A. Ndao, L. Hsu, et al., 'Symmetry-breaking-induced plasmonic exceptional points and nanoscale sensing,' Nat. Phys. 16 , 462 (2020)
- [15] H. Zhao, Z. Chen, R. Zhao, L. Feng, 'Exceptional points engineered glass slide for microscopic thermal mapping,' Nat. Commun. 9 , 1764 (2018)
- [16] M.-A. Miri and A. Al`, 'Exceptional points in optics and photonics,' Science 363, eaar7709 (2019) u
- [17] J. Wiersig, 'Review of exceptional point-based sensors,' Photonics Res. 8 , 1457 (2020)
- [18] J. Wiersig, 'Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: application to microcavity sensors for single-particle detection,' Phys. Rev. Lett. 112 , 203901 (2014)
- [19] J. Wiersig, 'Sensors operating at exceptional points: general theory,' Phys. Rev. A 93 , 033809 (2016)
- [20] Z.-P. Liu, et al., 'Metrology with PT-symmetric cavities: enhanced sensitivity near the PT-phase transition,' Phys. Rev. Lett. 117 , 110802 (2016)
- [21] W. Chen, S. K. Ozdemir, G. Zhao, J. Wiersig, L. Yang, 'Exceptional points enhance sensing in an optical microcavity,' Nature 548 , 192 (2017)
- [22] F. Vollmer and L. Yang, 'Label-free detection with high-Q microcavities: a review of biosensing mechanisms for integrated devices,' Nanophotonics 1 , 267 (2012)
- [23] L. Rondin, J.-P. Tetienne, T. Hingant, et al., 'Magnetometry with nitrogen-vacancy defects in diamond,' Rep. Prog. Phys. 77 , 056503 (2014)
- [24] W. W. Chow, J. Gea-Banacloche, L. M. Pedrotti, et al., 'The ring laser gyro,' Rev. Mod. Phys. 57 , 61 (1985)
- [25] S. Sunada and T. Harayama, 'Design of resonant microcavities: application to optical gyroscopes,' Opt. Express 15 , 16245 (2007)
- [26] R. El-Ganainy, K. G. Makris, M. Khajavikhan, , et al., 'NonHermitian physics and PT symmetry,' Nat. Phys. 14 , 11 (2018)
- [27] L. Feng, Z. L. Wong, R. M. Ma, Y. Wang, X. Zhang, 'Single-mode laser by parity-time symmetry breaking,' Science 346 , 972-975 (2014)
- [28] H. Hodaei, M. A. Miri, M. Heinrich, D. N. Christodoulides, M. Khajavikhan, 'Parity-time-symmetric microring lasers,' Science 346 , 975-978 (2014)
- [29] P. Miao, Z. Zhang, J. Sun, et al., 'Orbital angular momentum microlaser,' Science 353 , 464-467 (2016)
- [30] T. Lin, A. Liu, H. Li, et al., 'Exclusive generation of orbital angular momentum modes in parity-time symmetry fiber gratings,' Appl. Phys. Lett. 114, 141103 (2019)
- [31] C. N. Alexeyev, E. V. Barshak, B. P. Lapin, and M. A. Yavorsky, 'Generation of optical vortices in non-parity-timesymmetric chiral-core optical fibers,' Opt. Lett. 46 , 4474-4477 (2021)
- [32] W. Langbein, 'No exceptional precision of exceptional-point sensors,' Phys. Rev. A 98 , 023805 (2018)

- [33] H.-K. Lau and A. A. Clerk, 'Fundamental limits and nonreciprocal approaches in non-Hermitian quantum sensing,' Nat. Commun. 9 , 1 (2018)
- [34] C. Chen, L. Jin, and R.-B. Liu, 'Sensitivity of parameter estimation near the exceptional point of a non-Hermitian system,' New J. Phys. 21 , 083002 (2019)
- [35] M. Zhang, W. Sweeney, C. W. Hsu, et al., 'Quantum Noise Theory of Exceptional Point Amplifying Sensors,' Phys. Rev. Lett. 123 , 180501 (2019)
- [36] H. Wang, Y.-H. Lai, Z. Yuan, M.-G. Suh, and K. Vahala, 'Petermann-factor sensitivity limit near an exceptional point in a Brillouin ring laser gyroscope,' Nat. Commun. 11, 1 (2020)
- [37] D. Smith, H. Chang, E. Mikhailov, and S. Shahriar, 'Beyond the Petermann limit: Prospect of increasing sensor precision near exceptional points,' Phys. Rev. A 106, 013520 (2022)
- [38] J. Wiersig, 'Robustness of exceptional-point-based sensors against parametric noise: The role of Hamiltonian and Liouvillian degeneracies,' Phys. Rev. A 101, 053846 (2020)
- [39] R. Kononchuk, J. Cai, F. Ellis, R. Thevamaran, and T. Kottos, 'Enhanced signal-to-noise performance of EP-based electromechanical accelerometers,' Preprint ArXiv:2201. 13328 (2022)
- [40] D. Anderson , M. Shah, and L. Fan, 'Clarification of the Exceptional-Point Contribution to Photonic Sensing,' Phys. Rev. Appl. 19 , 034059 (2023)
- [41] M. Bukov, L. D'Alessio, A. Polkovnikov, 'Universal high-frequency behavior of periodically driven systems: From dynamical stabilization to Floquet engineering,' Adv. Phys., 64 , 139-226 (2015).
- [42] M. Holthaus, 'Floquet engineering with quasienergy bands of periodically driven optical lattices,' J. Phys. B At. Mol. Opt. Phys. 49 , 013001 (2015)
- [43] A. Eckardt, 'Atomic quantum gases in periodically driven optical lattices,' Rev. Mod. Phys. 89 , 011004 (2017)
- [44] P. Weinberg, M. Bukov, L. D'alessio, et al., 'Adiabatic perturbation theory and geometry of periodically-driven systems,', Phys. Rep. 688 , 3 (2017)
- [45] T. Oka, S. Kitamura, 'Floquet Engineering of Quantum Materials,' Annual Rev. Condens. Matter Phys. 10 , 387-408 (2019)
- [46] L. Zhou, D. J. Zhang, 'Non-Hermitian Floquet Topological Matter-A Review,' Entropy 25 , 1401 (2023)
- [47] W. P. Su, J. R. Schrieffer, A. J. Heeger, 'Solitons in Polyacetylene,' Phys. Rev. Lett. 42 , 1698 (1979)
- [48] J. K. Asb´th, L. Oroszl´ny ,A. P´lyi, o a a A Short Course on Topological Insulators: Band Structure Topology and Edge States in One- and Two-Dimensional Systems (Springer, Berlin, 2016)

[49]

N. Batra,

G. Sheet,

'Physics with

Coffee and Doughnuts:

Understanding the Physics Behind Topological

Insulators

Through Su-Schrieffer-Heeger Model,' Resonance

25

,

765 (2020)

- [50] M. Atala, M. Aidelsburger, J. T. Barreiro, et al., 'Direct measurement of the Zak phase in topological Bloch bands,' Nat. Phys. 9 , 795-800 (2013)
- [51] M. Karski, . F¨rster, J.-M. Choi, et al., 'Quantum Walk in Position Space with Single Optically Trapped Atoms,' Science o 325 , 174-177 (2009)
- [52] M. A. Broome, A. Fedrizzi, B. P. Lanyon, et al., 'Discrete Single-Photon Quantum Walks with Tunable Decoherence,' Phys. Rev. Lett. 104 , 153602 (2010)
- [53] A. Schreiber, K. N. Cassemiro, V. Potoˇek, et al., 'Photons Walking the Line: c A Quantum Walk with Adjustable Coin Operations,' Phys. Rev. Lett 104 , 1050502 (2010)
- [54] T. Kitagawa, M. A. Broome, A. Fedrizzi, et al., 'Observation of topologically protected bound states in photonic quantum walks,' Nat. Commun. 3:882 (2012)
- [55] S. Liu, W. Gao, Q. Zhang, et al., 'Topologically Protected Edge State in Two-Dimensional Su-Schrieffer-Heeger Circuit,' Research, 2019 , 8609875 (2019)

[56]

P. St-Jean, V. Goblot, E. Galopin, et al., 'Lasing in topological edge states of a one-dimensional lattice,' Nature Photonics

11

,

651-656 (2017)

- [57] L. Ge, Z. Gao, and L. Feng, 'Non-Hermitian gauged laser arrays with localized excitations: Anomalous threshold and generalized principle of selective pumping,' Phys. Rev. B 108 , 104111 (2023)
- [58] S. Yao, Z. Wang, 'Edge States and Topological Invariants of Non-Hermitian Systems,' Phys. Rev. Lett. 121 , 086803 (2018)
- [59] F. Song, S. Yao, Z. Wang, 'Non-Hermitian Skin Effect and Chiral Damping in Open Quantum Systems,' Phys. Rev. Lett. 123 , 170401 (2019)
- [60] Y. He, C. C. Chien, 'Non-Hermitian generalizations of extended Su-Schrieffer-Heeger models,' J. Phys.: Condens. Matter 33 085501 (2021)
- [61] B. Zhu, R. L¨, S. Chen, 'PT symmetry in the non-Hermitian Su-Schrieffer-Heeger model with complex boundary potenu tials,' Phys. Rev. A 89 , 062102 (2014)
- [62] S. Lieu, 'Topological phases in the non-Hermitian Su-Schrieffer-Heeger model,' Phys. Rev. B 97, 045106 (2018)
- [63] H. X. Cui, Xue-Wei Cao, M. Kang, et al., 'Exceptional points in extraordinary optical transmission through dual subwavelength metallic gratings,' Opt. Express 21 , 13368-13379 (2013)
- [64] A. Kodigala, T. Lepetit, B. Kant´, 'Exceptional points in three-dimensional plasmonic nanostructures,' Phys. Rev. B e 94 , 201103(R) (2016)
- [65] H. Lourenco-Martins, P. Das, L. H. G. Tizei, R. Weil, M. Kociak, 'A Self-hybridization within non-Hermitian localized plasmonic system,' Nat. Phys. 14 , 360-364 (2018)
- [66] A. Ndao, J. Park, L. Hsu, and B. Kant´, 'Observation of plasmonic exceptional point and nano-scale sensing,' in Optical e Sensors and Sensing Congress, OSA Technical Digest (Optica Publishing Group, 2020), paper SM3B.2.

- [67] V. M. Martinez Alvarez, J. E. Barrios Vargas, L. F. Torres, 'Non-Hermitian robust edge states in one dimension: Anomalous localization and eigenspace condensation at exceptional points,' Phys. Rev. B 97 , 121401 (2016)
- [68] B. Qi, L. Zhang, L. Ge, 'Defect States Emerging from a Non-Hermitian Flatband of Photonic Zero Modes,' Phys. Rev. Lett. 120 , 093901 (2018)
- [69] L. Ge, 'Symmetry-protected zero-mode laser with a tunable spatial profile,' Phys. Rev. A 95 , 023812 (2017)
- [70] J. D. H. Rivero, C. Fleming, B. Qi, L. Feng, L. Ge, 'Robust Zero Modes in Non-Hermitian Systems without Global Symmetries,' Phys. Rev. Lett. 131 , 223801 (2023)
- [71] F. Hentinger, et al., 'Direct observation of zero modes in a non-Hermitian optical nanocavity array,' Phot. Res. 10 , 574 (2022)
- [72] N. Hatano and D. R. Nelson, 'Vortex pinning and non-Hermitian quantum mechanics,' Phys. Rev. B 56 , 8651 (1997)
- [73] N. Hatano and D. R. Nelson, 'Non-Hermitian delocalization and eigenfunctions,' Phys. Rev. B 58 , 8384 (1998)
- [74] L. Ge, 'Non-Hermitian Systems with a Real Spectrum and Selective Skin Effect,' Innov. Discov. 1 , 4 (2024)
- [75] Z. Gao, et al., 'Two-Dimensional Reconfigurable Non-Hermitian Gauged Laser Array,' Phys. Rev. Lett. 130 , 263801 (2023)
- [76] S. Longhi, D. Gatti, G. Della Valle, 'Robust light transport in non-Hermitian photonic lattices,' Sci. Rep. 5 , 13376 (2015)
- [77] L. Zhang, et al., 'Acoustic non-Hermitian skin effect from twisted winding topology,' Nature Comm. 12 , 6297 (2021)
- [78] W. Wang, X. Wang, G. Ma, 'Non-Hermitian morphing of topological modes,' Nature 608 , 50 (2022)
- [79] R. S. Savelev, M. A. Gorlach, 'Topological states in arrays of optical waveguides engineered via mode interference,' Phys. Rev. B 102 , 161112(R)(2020)
- [80] T. G. Rappoport, Y. V. Bludov, F. H. L. Koppens, N. M. R. Peres, 'Topological graphene plasmons in a plasmonic realization of the Su-Schrieffer-Heeger model,' ACS Photonics 8 , 1817 (2021)
- [81] M. Pieczarka, 'Topological phase transition in an all-optical exciton-polariton lattice,' Optica 8 , 1084-1091 (2021)
- [82] C. H. Lee, S. Imhof, C. Berger et al, 'Topolectrical Circuits,' Comm. Phys. 1 , article no. 39 (2018)
- [83] H. Tian, J. Liu, A. Siddarth, et al., 'Magnetic-free silicon nitride integrated optical isolator,' Nat. Photon. 15 , 828 (2021)
David S. Simon Christopher R. Schwarze Abdoulaye Ndao Alexander V. Sergienko
Exceptional points in SSH-like models with hopping amplitude gradient
The Su-Schrieffer-Heeger (SSH) system is a popular model for exploring
topological insulators and topological phases in one dimension. Recent interest
in exceptional points has led to re-examination of non-Hermitian
generalizations of many physical models, including the SSH model. In such
non-Hermitian systems, singular points called exceptional points (EPs) appear
that are of interest for applications in super-resolution sensing systems and
topological lasers. Here, a non-Hermitian and non-PT-symmetric variation of the
SSH model is introduced, in which the hopping amplitudes are non-reciprocal and
vary monotonically along the chain. It is found that, while the existence of
the EPs is due to the nonreciprocal couplings, the number, position, and order
of the EPs can all be altered by the addition of the hopping amplitude
gradient, adding a new tool for tailoring the spectrum of a non-Hermitian
system.