---
title: 'Hatano-Nelson Chain: A Non-Hermitian Lattice Model'
url: https://www.emergentmind.com/topics/hatano-nelson-chain
type: topic
---

# Hatano-Nelson Chain: A Non-Hermitian Lattice Model

The Hatano-Nelson chain is a paradigmatic one-dimensional non-Hermitian lattice model with asymmetric nearest-neighbor hopping. In its standard form, nonreciprocity is introduced either through unequal left- and right-hopping amplitudes or, equivalently, through an imaginary vector potential, and this produces complex spectra, point-gap topology, exceptional points, and the non-Hermitian skin effect. Across recent work, the chain has also become a reference setting for disorder-induced spectral transitions, interaction-driven many-body phenomena, nonlinear skin modes, frequency-dependent non-Markovian transport, and experimental realizations in photonics and electric circuits [2502.10494] [2201.12653] [2603.24642].

## 1. Canonical formulation

A common interacting form of the Hatano-Nelson Hamiltonian is
\[
\hat{H}=\sum_{\ell}[(t+\gamma)\hat{c}_{\ell}^{\dagger}\hat{c}_{\ell+1}+(t-\gamma)\hat{c}_{\ell+1}^{\dagger}\hat{c}_{\ell}+U\hat{n}_{\ell}\hat{n}_{\ell+1}],
\]
where \(t\) is reciprocal hopping, \(\gamma\) is the nonreciprocal hopping parameter, and \(U>0\) is a repulsive nearest-neighbor interaction for spinless fermions [2201.12653]. An equivalent open-chain parameterization used in many-body studies is
\[
H = \sum_{n=1}^{N-1} \frac{J}{2} e^{ah} c_n^\dagger c_{n+1} + \frac{J}{2} e^{-ah} c_{n+1}^\dagger c_n + U c_n^\dagger c_n c_{n+1}^\dagger c_{n+1},
\]
with \(h\) the imaginary vector potential [2208.12017].

In the non-interacting nearest-neighbor chain, boundary conditions immediately alter the spectrum. For the formulation with \(t_r=1\) and \(t_l=\alpha^2<1\), periodic boundary conditions give
\[
\sigma(H_0)= e^{ik} + \alpha^2 e^{-ik},
\]
which is an ellipse in the complex plane, while open boundary conditions give
\[
\sigma(H_0)= 2\alpha \cos k,
\]
a real line segment between \(-2\alpha\) and \(2\alpha\) [2502.10494]. In the clean unidirectional limit, the periodic spectrum is a circle,
\[
E(k)= t_1 - t_2 e^{ik},
\]
and the eigenstates are extended [2601.07236].

These elementary spectral facts organize most subsequent generalizations. The Hatano-Nelson chain is therefore used as a prototypical one-dimensional non-Hermitian model exhibiting the non-Hermitian skin effect and complex energy-band topology, and as the fundamental non-Hermitian building block of the Su-Schrieffer-Heeger model [2603.24642] [2601.07236].

## 2. Boundary sensitivity, point-gap topology, and skin localization

The defining feature of the Hatano-Nelson chain is extreme sensitivity to boundary conditions. Under open boundary conditions, bulk eigenstates in the non-interacting limit become exponentially localized at a boundary, producing the non-Hermitian skin effect; under periodic boundary conditions, the same system instead exhibits complex spectral loops [2208.12017]. In generalized settings, this boundary sensitivity is encoded by spectral winding numbers such as
\[
\nu(E_0) = \frac{1}{2\pi i}\int_0^{2\pi} \frac{d}{d\Phi}\ln\det[H(\Phi)-E_0]\,d\Phi,
\]
or equivalent momentum-space expressions used to diagnose point-gap topology [2601.07236].

In the unidirectional chain with binary diagonal disorder, the complex spectrum forms a single loop for weak disorder, bifurcates at a critical threshold, and becomes two separated loops in the strong-disorder regime. Correspondingly, the winding number changes from \(\nu=1\) in the weak-disorder phase, through \(\nu=1/2\) at criticality, to \(\nu=0\) in the strong-disorder limit [2601.07236]. The same analysis shows that most states are exponentially localized, but at weak and critical disorder two special states remain completely delocalized, with diverging localization length. The divergence is directly correlated with non-trivial spectral winding [2601.07236].

Strictly open boundaries can eliminate the loop structure entirely in this disordered unidirectional setting. For \(\alpha=0\), the spectrum collapses into two flat bands at \(E=\pm h\), the winding number loses meaning, and the delocalized states disappear; for any \(\alpha>0\), the thermodynamic spectrum and localization length recover the periodic-boundary result [2601.07236]. This establishes that spectral topology in the Hatano-Nelson chain is not merely a bulk property in the Hermitian sense, but is inseparable from the boundary condition used to define the spectrum.

A frequent oversimplification is that the skin effect must always be unidirectional. In extensions with non-reciprocal next-nearest-neighbor hopping and structured onsite potentials, the open-boundary eigenstates can instead display a bidirectional non-Hermitian skin effect, localizing on both edges, and with a periodic onsite potential the direction of skin accumulation can reverse completely as the potential strength is varied [2403.05382].

## 3. Disorder, quasiperiodicity, impurities, and structured potentials

The Hatano-Nelson chain supports several distinct routes to spectral reconstruction. For strictly ergodic potentials, the spectrum admits an explicit Lyapunov-exponent characterization,
\[
\Sigma(g)=\mathcal{E}_0 \cup (\Sigma(0)\cap \mathcal{E}_+),
\]
with \(\mathcal{E}_0=\{E:L(E)=g\}\) and \(\mathcal{E}_+=\{E:L(E)>g\}\) [2311.09899]. This yields a sharp real-complex transition: for \(g\le g_{cr}\), the spectrum remains real; for \(g_{cr}<g<g_{er}\), it contains both real and complex parts; and for \(g>g_{er}\), it is purely complex [2311.09899]. If the Lyapunov exponent is continuous, the infinite-chain spectrum on \(\ell^2(\mathbb{Z})\) is approximated by finite-interval truncations with periodic boundary conditions, a behavior explicitly contrasted with the random-potential case [2311.09899].

Adding a periodic modulation to the onsite energies introduces a second insulating mechanism besides disorder. In the period-2 model,
\[
\mathcal{H} = \sum_{x=1}^L
\left[-\frac{t}{2}\left( e^h c^\dagger_{x+1} c_x + e^{-h} c^\dagger_{x} c_{x+1} \right) + \mu_x n_x + V \cos(\pi x) n_x \right],
\]
the periodic term opens a gap in the clean limit when \(V>V_c^{\rm o}=t\sinh(h)\), with
\[
W_{\rm G}=2\sqrt{V^2-t^2\sinh^2(h)},
\]
while disorder promotes localization and can even close a previously opened gap before full localization sets in [1009.1946]. The resulting spectral topologies include \(E\), \(LEL\), \(EGE\), \(LELGLEL\), \(LELEL\), \(L\), and \(LGL\) [1009.1946].

A single long-range impurity coupling is not a perturbatively small deformation in this model. For
\[
H = H_0 + t_{pq} c_p^\dagger c_q,
\]
the thermodynamic spectrum becomes
\[
\lim_{L\rightarrow\infty} \sigma(H) = \sigma(H_0) \cup \left\{
\lambda \in \mathbb{C} \setminus \sigma(H_0) : d_{qp}(\lambda) = -\frac{1}{t_{pq}}
\right\},
\]
and the additional solutions generate tentacle-like spectral wings emerging from the Bloch ellipse or the open-boundary non-Bloch segment [2502.10494]. The number of impurity-induced bound states is \(|p-q|\), and their localization length scales with the impurity-coupling distance rather than with system size [2502.10494]. This makes long-range impurity coupling a nonperturbative control parameter for both spectrum and wavefunction structure.

## 4. Interactions and many-body Hatano-Nelson physics

With repulsive nearest-neighbor interactions at half filling, the interacting Hatano-Nelson chain displays two \(\mathcal{PT}\) transitions as the interaction strength increases. The first transition occurs at an exceptional point between the first and second excited states in finite size and coincides with a first-order symmetry-breaking transition from a gapless phase into a gapped charge-density-wave regime; the persistent current abruptly vanishes at this point. A second transition at stronger interaction collapses all many-body eigenvalues onto the real axis, but its critical scale grows with system size and diverges in the thermodynamic limit, making it a finite-size effect [2201.12653].

Away from half filling and at strong interaction, the many-body spectrum stratifies into clusters centered at interaction energies such as \(\varepsilon_s=(N-s)U\). Under twisted boundary conditions, each cluster can exhibit a point gap with winding number
\[
\nu_s=\frac{1}{2\pi i}\int_0^{2\pi} d\phi \sum_j \partial_\phi \log\left[E_j(\phi)-\varepsilon_s\right],
\]
with \(\nu_1=\mathrm{sgn}(\gamma)N\) for the first cluster [2201.12653]. These point-gap windings indicate a many-body skin effect under open boundaries. At half filling and strong \(U\), by contrast, the spectrum collapses to open lines, the winding becomes ill-defined, and the many-body wave functions delocalize [2201.12653].

Open-chain many-body observables do not simply mirror single-particle skin localization. In the interacting model with open boundaries, the ground-state density profile becomes only slightly tilted relative to the average filling, Friedel oscillations exhibit a beating pattern, and the full counting statistics of particle number over any finite interval is Gaussian. The mean scales with the imaginary vector potential, while the variance is symmetric with respect to the chain center and independent of \(h\) [2208.12017]. This directly supports the statement that many-body effects can significantly alter and conceal the single-particle properties and the skin effect in non-Hermitian systems [2208.12017].

Nonequilibrium many-body dynamics retains this tension between non-Hermitian driving and collective low-energy structure. In the Luttinger-liquid regime, switching the imaginary vector potential on or off produces spatio-temporal Friedel oscillations and ballistic light cones propagating from the open ends, with local currents of equal magnitude for both quench protocols; the long-wavelength continuity equation remains satisfied even under non-unitary evolution [2304.09688]. For a finite-time linear ramp of the imaginary vector potential, the excess energy becomes complex-valued although the instantaneous Hamiltonian retains the same real spectrum throughout; the approach to adiabaticity is slow, with a \(\tau^{-1}\) envelope and oscillation period \(2L/v\), while ramp durations commensurate with that period realize a shortcut to adiabaticity without auxiliary controls [2408.07122].

## 5. Generalized Hatano-Nelson chains

Coupling two quasiperiodic Hatano-Nelson chains changes the character of the delocalization-localization transition. In the cross-coupled model, two critical potentials \(V_{c1}<V_{c2}\) appear:
\[
V_{c1}=2\max(|t_L-u_2|,|t_R-u_1|),\qquad
V_{c2}=2\max(|t_L+u_2|,|t_R+u_1|),
\]
with all states delocalized below \(V_{c1}\), all states localized above \(V_{c2}\), and two mobility edges symmetrically placed about \(\mathrm{Re}[E]=0\) in between [2409.04417]. In this intermediate regime, the mobility edges divide the spectrum in equal proportions, with \(50\%\) localized and \(50\%\) delocalized states [2409.04417]. The same work shows that the usual one-to-one correspondence between periodic-boundary states and open-boundary skin states breaks down in the coupled system [2409.04417].

Longer-range asymmetric hopping yields another major extension. In the generalized model
\[
H_{lr} = \sum_{n=1}^N \left( t_l \, c_n^\dagger c_{n+l} + t_{-r} \, c_n^\dagger c_{n-r} \right),
\]
with \(l\ge r\ge 1\) and \(\gcd(l,r)=1\), open boundaries can host exceptional points of arbitrary order that do not scale with system size [2403.12018]. Their order is controlled by the sublattice structure and by \(N \bmod (l+r)\); the largest possible order is \(m=l+r-1\) for \(N\equiv -1 \pmod{l+r}\) [2403.12018]. These exceptional-point eigenstates can have support on only a subset of sites, still exhibit skin localization, and can be tuned to localize at the opposite edge from all remaining states. They are robust against generic hopping perturbations and against a specific class of on-site disorder [2403.12018].

Further generalization to matrix-valued couplings leads to the non-Abelian Hatano-Nelson model with a nonreciprocal \(U(2)\) gauge field. Its Bloch Hamiltonian,
\[
H(k)= t_0\, \mathbf{d}_R\!\cdot\!\boldsymbol{\sigma}
+ t_L e^{ik}\, \mathbf{d}_L\!\cdot\!\boldsymbol{\sigma}
+ t_R e^{-ik}\, \mathbf{d}_R\!\cdot\!\boldsymbol{\sigma},
\]
produces Hopf-link-shaped complex energy braiding with braiding degree \(|v|=2\) and a bipolar skin effect in which left- and right-localized skin modes coexist for the same parameters, a phenomenon impossible in the Abelian nearest-neighbor case [2603.24642]. In a two-orbital interacting ladder with opposite non-Hermiticity on the two legs, Hermitian interchain hopping can restore a fully real non-interacting spectrum when \(V_0\ge 2\delta\), while interactions generate detached doublon branches and winding-number structures linking periodic-boundary doublons to open-boundary skin modes [2604.14533].

## 6. Nonlinear, non-Markovian, and experimental directions

Kerr nonlinearity does not destroy the skin effect but reorganizes it. In the nonlinear open-chain equation
\[
i\frac{d\psi_n}{d\tau} = C\left(\psi_{n+1} + t\psi_{n-1}\right) + \sigma|\psi_n|^2 \psi_n,
\]
families of nonlinear skin modes bifurcate from the linear skin modes at any non-reciprocal strength [2311.09139]. For focusing nonlinearity, these modes become more localized and connect continuously to discrete solitons in the anti-continuum limit; for defocusing nonlinearity, they broaden relative to the linear modes [2311.09139].

Under periodic boundary conditions, the nonlinear Hatano-Nelson model exhibits a different phenomenon. All nonlinear plane waves are modulationally unstable, and the total norm
\[
P(t)=\sum_n |\psi_n(t)|^2
\]
undergoes an initial growth followed by an abrupt halt: a dynamical growth blockade [2501.01226]. The paper interprets this as a stopping of convective motion caused by self-induced disorder, rather than by boundaries or external randomness [2501.01226].

A non-Markovian generalization replaces constant dissipation by frequency-dependent dissipation and frequency-dependent nonreciprocal hopping,
\[
t_{\pm}(z) = -g - i e^{\mp i\phi} \frac{\Gamma(z)}{2},
\]
derived microscopically by integrating out dissipative auxiliary sites with non-equilibrium Green’s functions [2511.05328]. This produces purely non-Markovian effects: unidirectional frequency blocking in the bosonic setting and a nonequilibrium dissipative quantum phase transition in the fermionic setting, both absent in Markovian theory [2511.05328].

Experimental realizations now access several of these phenomena directly. In integrated photonics, time-domain modulation in a lithium-niobate photonic molecule realizes dynamically tunable Hatano-Nelson couplings,
\[
\lambda_{12,21}=\lambda\left(1\pm\frac{\lambda \beta^2}{(4\lambda^2+\gamma_0^2)\gamma_0}\right),
\]
reaches an exceptional point, surpasses previous asymmetry records, and yields giant \(60\) dB optical contrast together with non-reciprocal \(\pi\)-phase contrast or photonic gyration [2410.10079]. In electric circuits, the non-Abelian Hatano-Nelson model has been implemented experimentally, and the Hopf-link-shaped admittance spectra and bipolar skin admittance modes have been observed [2603.24642]. Experimental relevance has also been noted for quantum dot arrays, cold atom systems, and monitored quantum circuits in interacting settings [2201.12653].

Taken together, these developments establish the Hatano-Nelson chain not as a single isolated model but as a compact framework in which nonreciprocity, spectral topology, localization, many-body correlations, nonlinear self-organization, and open-system memory effects can all be formulated with unusual precision.

Source: https://www.emergentmind.com/topics/hatano-nelson-chain