Generalized Brillouin Zone (GBZ): Non-Hermitian Analysis

• Generalized Brillouin Zone (GBZ) is a complex-momentum framework for non-Hermitian systems that replaces the traditional unit circle with curves defined in the complex plane.
• It is constructed by solving the characteristic polynomial of non-Bloch Hamiltonians and equating moduli of ordered roots to faithfully reproduce open-boundary spectra and skin effects.
• The GBZ framework extends to disordered, higher-dimensional, and many-body models, providing critical insights into spectral, topological, and boundary phenomena in non-Hermitian lattices.

The generalized Brillouin zone (GBZ) is the central momentum-space object of non-Bloch band theory for non-Hermitian lattices. It replaces the conventional Brillouin zone by a closed curve, or more generally a set of loci, in the complex plane of a deformation parameter $\beta$, where $\beta$ generalizes $e^{ik}$ by allowing $|\beta|\neq 1$. In one-dimensional settings, open-boundary-condition (OBC) bulk solutions take the form $\psi_n\propto \beta^n$, and the modulus of $\beta$ controls non-Hermitian skin accumulation through the decay length $\xi(E)=1/|\ln|\beta||$: $|\beta|=1$ gives Bloch-like extended states, while $|\beta|\neq 1$ gives boundary-localized states. The GBZ is therefore the non-Hermitian replacement of the unit circle for reconstructing OBC spectra, bulk states, and bulk-boundary correspondence (Liu et al., 2022, Yang et al., 2019).

1. Foundational definition and non-Bloch formulation

In a generic translationally invariant one-dimensional non-Hermitian tight-binding model with finite-range hopping, the generalized Bloch Hamiltonian is written as

$H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$

with $\beta$0. The OBC spectral problem is encoded by

$\beta$1

After clearing poles at $\beta$2, this becomes

$\beta$3

where $\beta$4 is polynomial in $\beta$5 and $\beta$6, and $\beta$7 is the pole order. The roots $\beta$8 of $\beta$9 are ordered by modulus, and the GBZ is the locus on which the middle roots satisfy the equal-modulus condition

$e^{ik}$0

This is the basic non-Bloch replacement of the Hermitian condition $e^{ik}$1 (Yang et al., 2019).

For nearest-neighbor non-reciprocal chains, the construction is especially transparent. If

$e^{ik}$2

then the characteristic equation is quadratic in $e^{ik}$3, and the GBZ is circular: $e^{ik}$4 The OBC spectrum is then reproduced by evaluating the non-Bloch Hamiltonian along that circle. This simple case already exhibits the central physical content of the GBZ: the modulus $e^{ik}$5 encodes skin accumulation, while the phase $e^{ik}$6 still parameterizes a band-like dispersion (Lu et al., 2 Jun 2026).

A rigorous mathematical formulation for non-reciprocal tridiagonal $e^{ik}$7-Toeplitz systems replaces real quasiperiodicities by complex quasiperiodicities $e^{ik}$8. In that framework, the OBC GBZ is a shifted contour, the periodic-boundary-condition (PBC) zone remains the classical real Brillouin zone, and the full semi-infinite Toeplitz spectrum is decomposed over a two-dimensional complex quasiperiodicity region. This establishes, in operator-theoretic terms, that complex deformation of the Brillouin zone is required to capture unidirectional decay and correct spectral convergence in non-reciprocal systems (Ammari et al., 2024).

2. Standard construction in translationally invariant one-dimensional systems

The standard GBZ construction proceeds from the characteristic polynomial $e^{ik}$9. For a model with left/right hopping ranges $|\beta|\neq 1$0 and $|\beta|\neq 1$1, one solves for all roots $|\beta|\neq 1$2, orders them by modulus, and imposes

$|\beta|\neq 1$3

As $|\beta|\neq 1$4 varies, the corresponding $|\beta|\neq 1$5 values trace a closed loop in the complex $|\beta|\neq 1$6-plane. Evaluating $|\beta|\neq 1$7 on that loop reproduces the OBC bulk spectrum, thereby restoring non-Bloch bulk-boundary correspondence (Liu et al., 2022).

For multiband systems, the structure is more intricate. An $|\beta|\neq 1$8-band non-Hermitian Hamiltonian is constituted by $|\beta|\neq 1$9 distinct sub-GBZs, each of which is a piecewise analytic closed loop. The full GBZ is the union of these sub-GBZs, and the OBC spectrum is reconstructed as

$\psi_n\propto \beta^n$0

The piecewise nature reflects changes in the ordering of $\psi_n\propto \beta^n$1 at branch points, where analytic arcs are stitched together into closed contours (Yang et al., 2019).

A systematic analytic construction is provided by the auxiliary generalized Brillouin zone (aGBZ). One introduces

$\psi_n\propto \beta^n$2

and eliminates $\psi_n\propto \beta^n$3 to obtain a real algebraic curve

$\psi_n\propto \beta^n$4

The aGBZ encodes all equal-modulus pairings of roots; the physical GBZ is then selected by keeping the arcs corresponding to the appropriate root-ordering label $\psi_n\propto \beta^n$5. This resultant-based method turns GBZ construction into an algebraic-geometry problem and avoids direct large-$\psi_n\propto \beta^n$6 OBC diagonalization (Yang et al., 2019).

Topological invariants are defined on the GBZ rather than on the unit circle. For chiral-symmetric models with off-diagonal blocks $\psi_n\propto \beta^n$7, a GBZ winding can be written in terms of zeros and poles inside sub-GBZs, and more general point-gap windings can be expressed as contour integrals over the GBZ. This is the non-Bloch analogue of conventional Bloch topological classification (Yang et al., 2019).

3. Analytic structure, singularities, and disconnected topology

The GBZ is often depicted as a single smooth loop, but the analytic structure is more varied. Even in translationally invariant models, the GBZ is only piecewise analytic in general, because the equal-modulus condition segments the algebraic curve into arcs across which the ordering of $\psi_n\propto \beta^n$8 is fixed. Cusps, self-intersections, and branch points occur where that ordering changes or where the map $\psi_n\propto \beta^n$9 becomes non-invertible (Yang et al., 2019).

A first major refinement is that the GBZ can become singular without destroying energetic bulk-boundary correspondence. In non-Hermitian SSH-type chains, the GBZ can collapse to a single point such as $\beta$0, yet the non-Bloch continuum bands $\beta$1 computed on that degenerate GBZ still reproduce the OBC spectrum in the thermodynamic limit. By contrast, the non-Bloch winding number can become ill-defined in such regimes. The distinction is explicit: energetic recovery of OBC spectra can survive GBZ singularity, whereas topological invariants can fail when the GBZ collapses, when one of $\beta$2 becomes constant, or when analyticity is lost (Guo et al., 2021).

A second refinement is topological disconnectedness. Sufficient conditions for connectivity have been proved for single-band models: if the OBC spectrum is simply connected, then the GBZ is connected; and if the coupling range satisfies $\beta$3 and $\beta$4, the GBZ is connected. These are sufficient, not necessary, conditions. Beyond them, disconnected GBZs occur, and the number of connected components can exceed the number of bands. This challenges the common assumption that the GBZ can always be indexed monotonically by $\beta$5 like the Hermitian Brillouin zone (Wang et al., 8 Oct 2025).

The topology of disconnected GBZs has concrete spectral consequences. In single-band models, one connected component can cover some OBC spectral arcs twice while another disconnected component is required to cover complementary arcs with opposite orientation. In sublattice-symmetric two-band constructions, this can close an OBC line gap without changing the point-gap topology. The corresponding GBZ contour must then be understood as a union of components with specific orientations rather than as a single loop (Wang et al., 8 Oct 2025).

Generalized boundary conditions introduce another kind of GBZ topology. In non-reciprocal chains with boundary hoppings and onsite boundary terms, the boundary matrix selects generalized momenta $\beta$6, and under intermediate boundary conditions between OBC and PBC multiple disjoint GBZs can appear. Topological phase transitions are then characterized by generalized-momentum touching of GBZs, which manifests as exceptional points. In that setting, the relevant winding numbers are constructed from meromorphic boundary data $\beta$7, rather than solely from bulk Bloch data (Verma et al., 2023).

A common misconception is therefore that the GBZ is necessarily a unique, smooth, bulk-defined loop. The accumulated results show instead that it can be piecewise analytic, singular, disconnected, or boundary-sensitive, while still retaining predictive power for OBC spectra in appropriately defined regimes (Guo et al., 2021, Wang et al., 8 Oct 2025).

4. Disorder, broken translational symmetry, and modified GBZ theory

Standard GBZ theory presumes spatial homogeneity, because the equal-modulus rule is derived from a translationally invariant characteristic equation. Once disorder breaks translational symmetry, $\beta$8 in the conventional sense is no longer directly available, and the clean equal-modulus construction ceases to be adequate. This difficulty is already visible at the level of bulk-bulk correspondence (BBC): PBC and OBC spectra no longer match in non-Hermitian systems with disorder, even though OBC eigenvalues remain insensitive to certain similarity transforms (Zhang et al., 2022).

A determinant-based optimization framework reconstructs BBC in disordered non-Hermitian chains. For a nearest-neighbor disordered chain with asymmetric hoppings $\beta$9, a diagonal similarity transform rescales the hoppings as

$\xi(E)=1/|\ln|\beta||$0

while preserving the OBC determinant. The exact sitewise optimum is

$\xi(E)=1/|\ln|\beta||$1

which makes $\xi(E)=1/|\ln|\beta||$2. A global effective parameter $\xi(E)=1/|\ln|\beta||$3 is obtained from

$\xi(E)=1/|\ln|\beta||$4

and the NHSE strength is quantified by $\xi(E)=1/|\ln|\beta||$5. This replaces the clean GBZ contour by an optimization-selected effective deformation (Zhang et al., 2022).

For on-site disorder, a modified GBZ theory constructs a $\xi(E)=1/|\ln|\beta||$6-deformed periodic Hamiltonian $\xi(E)=1/|\ln|\beta||$7 by the mapping $\xi(E)=1/|\ln|\beta||$8 while leaving on-site disorder untouched. The deformation parameter is selected by minimizing

$\xi(E)=1/|\ln|\beta||$9

For the disordered Hatano–Nelson chain, this reduces to

$|\beta|=1$0

with

$|\beta|=1$1

In the clean limit, minimizing $|\beta|=1$2 reproduces the conventional result

$|\beta|=1$3

but with disorder the energy-dependent term $|\beta|=1$4 becomes comparable and produces a plateau $|\beta|=1$5 of nearly minimal $|\beta|=1$6 values rather than a unique minimum (Liu et al., 2022).

The resulting prescription is state-resolved. For each OBC eigenvalue $|\beta|=1$7, one first finds the plateau $|\beta|=1$8 from the near-minima of $|\beta|=1$9, and then chooses

$|\beta|\neq 1$0

to minimize $|\beta|\neq 1$1. This selects the deformation closest to the no-skin point $|\beta|\neq 1$2, yields a faithful reconstruction of the OBC eigenstate from the $|\beta|\neq 1$3-deformed periodic eigenstate via a diagonal similarity transform, and defines the skin depth through

$|\beta|\neq 1$4

If the plateau crosses $|\beta|\neq 1$5, NHSE is suppressed for that state (Liu et al., 2022).

This modified construction also resolves two issues that are inaccessible to the conventional clean GBZ. First, in long-range hopping models, clean systems typically require an energy-dependent $|\beta|\neq 1$6, but disorder broadens the minima into plateaus that can intersect at a single unified $|\beta|\neq 1$7, thereby restoring a global BBC. Second, in quasi-one-dimensional ladders with Peierls phases $|\beta|\neq 1$8, magnetic flux can create plateaus containing $|\beta|\neq 1$9 for a subset of states, which suppresses NHSE even though the global minimum of $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$0 can remain near the clean $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$1. This explains why the conventional GBZ can appear insensitive to flux while the disorder-aware modified GBZ correctly captures flux-induced suppression of skin modes (Liu et al., 2022).

5. Boundary sensitivity, finite-size effects, and GBZ without conventional skin effect

The GBZ is often identified with the non-Hermitian skin effect itself, but that identification is not universal. A parity-induced even-odd phenomenon in a reciprocal bipartite chain can generate a nontrivial GBZ without conventional NHSE. In the reciprocal SSH$H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$2 chain with alternating complex-conjugate bonds, odd system size $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$3 yields real $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$4 and a conventional Brillouin zone, whereas even $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$5 forces complex solutions of

$H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$6

These complex generalized momenta reproduce the even-$H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$7 OBC spectrum, yet the eigenstates remain delocalized: the inverse participation ratio scales as $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$8 and the participation ratio grows as $H(\beta)=\sum_{m=-M}^{N} T_m \beta^m,$9. The resulting GBZ is therefore a finite-size, parity-driven boundary-sensitive structure rather than a manifestation of boundary accumulation (Felski, 29 May 2026).

Boundary Hamiltonians can also be used to manipulate the GBZ directly. In non-Hermitian electric circuits modeled by a Hatano–Nelson-type asymmetric chain with boundary hoppings $\beta$00 and boundary onsite terms $\beta$01, the admissible generalized momenta are selected by a $\beta$02 boundary matrix

$\beta$03

with

$\beta$04

Together with the bulk constraint $\beta$05, the condition $\beta$06 produces boundary-driven GBZ deformation. Multiple separated GBZ manifolds can then appear, containing both decaying and growing wave functions, and topological transitions occur when zeros and poles of the boundary functions cross the contour. The corresponding non-Bloch winding is formulated from

$\beta$07

This is an explicitly boundary-governed GBZ topology rather than a purely bulk-governed one (Kim et al., 21 Feb 2025).

Conditional boundary conditions provide another boundary-sensitive variant. By allowing only one direction of inter-edge hopping, one can selectively collapse spectral sectors of a PBC spectrum onto OBC sectors according to their spectral winding number, thereby producing a composite contour built from parts of the conventional Brillouin zone and parts of the GBZ. In that framework, Bloch points are the intersections between the PBC zone and the OBC GBZ, and the resulting “composite GBZ” depends explicitly on boundary control (Fu et al., 20 Jun 2025).

These results support a narrower interpretation of the GBZ: it is the correct complex-momentum construct for boundary-sensitive spectral reconstruction, but it is not by itself a sufficient diagnostic of conventional NHSE. Depending on symmetry, parity, and explicit boundary engineering, GBZ structures can arise with delocalized states, with separated decaying and growing manifolds, or with selective collapse of only part of the PBC spectrum (Felski, 29 May 2026, Kim et al., 21 Feb 2025).

6. Extensions: inhomogeneous, higher-dimensional, many-body, and entanglement-based GBZ frameworks

A major limitation of conventional GBZ theory is its reliance on a site-independent complex Bloch factor. For spatially inhomogeneous non-Hermitian systems, that assumption fails because non-Hermitian pumping competes with position-dependent energy scales and interference. A phase-space GBZ formalism replaces the global $\beta$08 by a local $\beta$09, with the local relation

$\beta$10

in a single-component inhomogeneous chain. The corresponding local branches can bifurcate, and eigenstates may jump between different GBZ branches at real-space positions $\beta$11. The jump locations act as an emergent degree of freedom, stabilize real spectral tails, and protect a class of zero modes in two-component chains when suitable domain-wall criteria are satisfied (Li et al., 16 Jan 2025).

In two dimensions, the GBZ becomes energy-resolved rather than a single universal contour. For a fixed OBC eigenvalue $\beta$12, the two-dimensional asymptotic GBZ is assembled from geometry-independent Bloch or non-Bloch Fermi points and geometry-dependent non-Bloch equal-frequency contours that connect them. A key result is that the region covered by the OBC spectrum on the complex-energy plane is independent of the open-boundary geometry, even though the equal-frequency contours and localization directions depend on boundary orientation and shape. This extends non-Bloch bulk-boundary correspondence beyond one dimension while keeping the basic idea that complexified momenta rather than real Bloch momenta control OBC bulk spectra (Xu et al., 2023).

The interacting many-body problem can also be incorporated, but at present only in restricted regimes. For circular GBZs, a diagonal many-body similarity transform

$\beta$13

maps the original OBC problem to a quasi-reciprocal many-body Hamiltonian under PBC,

$\beta$14

In the non-Hermitian SSH chain with density-density interactions, this construction preserves locality, enables exact diagonalization on the quasi-reciprocal model, and supports GBZ-based many-body Zak phases, biorthogonal structure factors, and entanglement-spectrum diagnostics. The method, however, relies on the GBZ being circular; non-circular GBZs generally generate long-range interactions under the inverse transform (Lu et al., 2 Jun 2026).

Entanglement provides a further GBZ-based diagnostic. In the non-Hermitian SSH model, a quasi-reciprocal lattice constructed by Fourier transformation on the GBZ allows one to define non-Bloch entanglement entropy. When the GBZ is circular, the resulting von Neumann entropy is real and positive-definite over wide parameter regions except near exceptional points, and each Fermi point contributes precisely $\beta$15 to the central charge of the logarithmic scaling. At the exceptional point, the central charge becomes negative due to the exceptional bound state. For non-circular GBZs, long-range hopping emerges in the quasi-reciprocal lattice and the bulk von Neumann entropy is no longer real, but the non-Bloch edge entanglement entropy remains real and continues to respect bulk-boundary correspondence as a topological indicator (Yang et al., 2024).

A rigorous operator-theoretic perspective complements these physical constructions. For non-reciprocal tridiagonal $\beta$16-Toeplitz systems, the generalized Brillouin zone is a complex quasiperiodicity region

$\beta$17

with the OBC limit selecting the shifted contour $\beta$18, the PBC limit retaining $\beta$19, and the full semi-infinite Toeplitz spectrum decomposed over the full two-dimensional GBZ. This yields a mathematically explicit distinction between OBC, PBC, and semi-infinite spectral limits (Ammari et al., 2024).

Across these extensions, the core idea remains unchanged: the conventional Brillouin zone is replaced by a complex deformation dictated by boundary conditions and non-Hermitian transport. What changes is the object that carries that deformation—an algebraic loop, a plateau of $\beta$20 values, a disconnected set of components, a position-dependent branch structure, or an energy-dependent higher-dimensional manifold.