Graded hopping screens nonreciprocity and reorganizes Stark asymptotics in a non-Hermitian Stark chain
Pith reviewed 2026-05-08 06:16 UTC · model grok-4.3
The pith
Graded hopping in non-Hermitian Stark chains screens nonreciprocity by turning exponential skin factors into algebraic accumulation with a threshold at |F1| = 2|F2|.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
An exact diagonal similarity transformation removes the bond asymmetry and converts the usual exponential skin factor into an algebraic boundary accumulation with exponent eta=gamma/F2. The transformed symmetric chain then reduces asymptotically to a constant-coefficient recurrence, giving the Stark threshold |F1|=2|F2|. The original right eigenstates acquire the unified envelope psi_j^R ~ j^eta phi_j, with oscillatory, double-root, and exponentially localized branches across the threshold. This form also yields two finite-size scales, one measuring the logarithmic screening of nonreciprocity and the other balancing the algebraic skin factor against the exponential Stark tail.
What carries the argument
The linearly graded hopping term that enables the diagonal similarity transformation symmetrizing the chain and reorganizing the large-position asymptotics of the eigenstates.
If this is right
- The right eigenstates of the original system have the envelope form j to the power eta times a transformed state.
- The system exhibits a Stark threshold at |F1/F2| = 2 separating algebraic skin accumulation from exponential localization.
- Edge polarization bends near the threshold and weakens beyond it on the localized side.
- The inverse participation ratio of the most localized states increases rapidly when F1/F2 exceeds 2.
- Half-chain entanglement growth after a charge-density-wave quench is enhanced above the threshold.
Where Pith is reading between the lines
- Graded hopping could serve as a general design tool to screen nonreciprocity in other non-Hermitian models beyond the Stark chain.
- Experimental realizations in optical lattices or photonic arrays with position-dependent couplings might directly observe the algebraic accumulation and the threshold.
- The two finite-size scales identified suggest that for intermediate system sizes the crossover between skin and Stark effects can be tuned independently.
Load-bearing premise
The asymptotic analysis at large positions and the reduction of the recurrence relation to constant coefficients apply across the full spectrum and for finite but large systems without new artifacts.
What would settle it
Exact diagonalization of finite chains for different values of F1 and F2 should reveal a qualitative change in the decay of eigenstate amplitudes at the ratio |F1/F2|=2, from power-law modulated to purely exponential.
Figures
read the original abstract
We study a one-dimensional non-Hermitian Stark chain in which nonreciprocal hopping, a linear potential, and linearly graded hopping act simultaneously. The central question is how boundary pumping and field-induced confinement are reorganized when the hopping amplitude itself grows with position. We show that the graded term separates the two localization channels at the level of the large-position asymptotics. An exact diagonal similarity transformation removes the bond asymmetry and converts the usual exponential skin factor into an algebraic boundary accumulation with exponent $\eta=\gamma/F_2$. The transformed symmetric chain then reduces asymptotically to a constant-coefficient recurrence, giving the Stark threshold $|F_1|=2|F_2|$. The original right eigenstates acquire the unified envelope $\psi_j^R\sim j^{\eta}\phi_j$, with oscillatory, double-root, and exponentially localized branches across the threshold. This form also yields two finite-size scales, one measuring the logarithmic screening of nonreciprocity and the other balancing the algebraic skin factor against the exponential Stark tail. A joint localization map in the $(\gamma,F_1/F_2)$ plane verifies this structure. The edge polarization bends near the Stark threshold and weakens on the localized side, while the inverse participation ratio of the most localized eigenstates rises rapidly for $F_1/F_2>2$. Using a normalized Gaussian projector appropriate for non-unitary evolution, we further show that the same threshold enhances half-chain entanglement growth after a charge-density-wave quench. These results identify graded hopping as a controlled mechanism for screening nonreciprocity, resetting Stark asymptotics, and organizing the finite-size crossover between algebraic skin accumulation and Stark localization.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper examines a one-dimensional non-Hermitian Stark chain combining nonreciprocal hopping, a linear potential, and linearly graded hopping. It claims that an exact diagonal similarity transformation eliminates bond asymmetry and converts the exponential skin factor into an algebraic boundary accumulation with exponent η=γ/F₂. The transformed symmetric chain is then asymptotically reduced to a constant-coefficient recurrence, yielding the Stark threshold |F₁|=2|F₂| that separates oscillatory, double-root, and exponentially localized branches of the unified envelope ψ_j^R ∼ j^η ϕ_j. These results are verified via a joint localization map in the (γ, F₁/F₂) plane, edge polarization behavior, inverse participation ratio of localized states, and enhanced half-chain entanglement growth after a charge-density-wave quench using a normalized Gaussian projector.
Significance. If the similarity transformation and asymptotic threshold derivation hold rigorously, the work identifies graded hopping as a tunable mechanism to screen nonreciprocity and reorganize Stark localization, providing concrete finite-size scales and a phase diagram for algebraic versus exponential accumulation. This could be significant for non-Hermitian lattice physics, particularly in controlling boundary effects and quench dynamics in open systems.
major comments (2)
- [asymptotic analysis and derivation of the Stark threshold] The asymptotic reduction to the constant-coefficient recurrence (leading to the quadratic r² + (F₁/F₂)r + 1 = 0 and the threshold |F₁|=2|F₂|) drops the E/j term and treats graded hopping as strictly linear at large j. This approximation is load-bearing for the claim that the threshold uniformly organizes the localization character of the full spectrum, yet its applicability to states with support near j=0, high-energy eigenstates, or finite-N systems is not demonstrated with bounds or exhaustive checks.
- [similarity transformation section] The exact diagonal similarity transformation is stated to map right eigenvectors to the algebraic envelope without new artifacts, but the manuscript does not explicitly verify that the localization classification (algebraic skin accumulation versus exponential Stark localization) is preserved for every eigenstate under finite-size boundary conditions.
minor comments (3)
- [notation and envelope definition] The exponent η=γ/F₂ and the form of the unified envelope should be defined with its domain of validity immediately after the transformation is introduced, rather than relying on the abstract.
- [finite-size scales discussion] The two finite-size scales (logarithmic screening of nonreciprocity and algebraic-versus-Stark balance) are mentioned but lack explicit expressions or scaling derivations that would aid reproducibility.
- [localization map figure] The joint localization map in the (γ, F₁/F₂) plane is a useful verification; ensure all axes, color scales, and contour definitions are labeled consistently with the text and that the threshold line is overlaid for direct comparison.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. We address each major point below, indicating revisions made to strengthen the manuscript.
read point-by-point responses
-
Referee: [asymptotic analysis and derivation of the Stark threshold] The asymptotic reduction to the constant-coefficient recurrence (leading to the quadratic r² + (F₁/F₂)r + 1 = 0 and the threshold |F₁|=2|F₂|) drops the E/j term and treats graded hopping as strictly linear at large j. This approximation is load-bearing for the claim that the threshold uniformly organizes the localization character of the full spectrum, yet its applicability to states with support near j=0, high-energy eigenstates, or finite-N systems is not demonstrated with bounds or exhaustive checks.
Authors: We agree that the reduction is an asymptotic approximation valid in the large-j regime, where the E/j term is subdominant and the graded hopping remains linear to leading order. This is the standard procedure for extracting the characteristic equation and the resulting threshold |F₁|=2|F₂| that governs the transition between oscillatory, double-root, and exponentially decaying branches of the envelope. The threshold organizes the leading asymptotics of the unified form ψ_j^R ∼ j^η ϕ_j. To address applicability in finite systems, we have performed exhaustive numerical diagonalizations across the full spectrum for multiple N (up to several hundred sites), including eigenstates with significant weight near j=0 and high-energy states. These confirm that the threshold continues to separate the localization regimes, with deviations controlled by the two finite-size scales already derived in the manuscript. We have added error estimates for the neglected terms and a new supplementary figure summarizing the classification fidelity versus N. revision: partial
-
Referee: [similarity transformation section] The exact diagonal similarity transformation is stated to map right eigenvectors to the algebraic envelope without new artifacts, but the manuscript does not explicitly verify that the localization classification (algebraic skin accumulation versus exponential Stark localization) is preserved for every eigenstate under finite-size boundary conditions.
Authors: The similarity transformation is constructed to be exact and diagonal, establishing a bijective mapping between the spectra and eigenvectors of the original and transformed operators with no additional artifacts introduced. Because the transformation is position-dependent but invertible for any finite N, the asymptotic classification derived from the transformed chain carries over directly to the original right eigenvectors. We have now added an explicit verification: a supplementary analysis that classifies every eigenstate in finite chains according to the algebraic versus exponential character of its envelope, confirming that the threshold |F₁|=2|F₂| organizes the localization type uniformly across the spectrum under open boundary conditions. revision: yes
Circularity Check
No significant circularity in the derivation chain
full rationale
The paper's central results are obtained via an explicit exact diagonal similarity transformation applied to the non-Hermitian Hamiltonian, which algebraically converts the skin factor, followed by direct large-position asymptotic analysis that reduces the transformed recurrence to a constant-coefficient form whose characteristic equation supplies the threshold |F1|=2|F2|. These manipulations follow from the model's defining equations without any reduction to fitted inputs, self-referential definitions, or load-bearing self-citations. The derivation remains self-contained against the stated Hamiltonian and recurrence relations.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The system is described by a one-dimensional tight-binding Hamiltonian with nonreciprocal nearest-neighbor hopping, a linear Stark potential, and linearly position-dependent hopping amplitude.
Reference graph
Works this paper leans on
-
[1]
Reduction of many-body evolution to a single-particle problem The Hamiltonian is quadratic, ˆH= X α,β c† αhαβcβ,(A1) wherehis the one-body matrix. The fermion operators satisfy {cα, c† β}=δ αβ,(A2) {cα, cβ}={c † α, c† β}= 0.(A3) A useful commutator is [c† αcβ, c† m] =c † α{cβ, c† m} − {c† α, c† m}cβ =δ βmc† α.(A4) It follows that [ ˆH, c† m] = X α,β hαβ[c...
-
[2]
We therefore de- rive the correlation matrix for a general full-rank orbital matrixU
Normalized correlation matrix for nonorthogonal occupied orbitals For a non-Hermitian Hamiltonian, the columns ofU(t) do not in general remain orthonormal. We therefore de- rive the correlation matrix for a general full-rank orbital matrixU. Define a† µ = X i Uiµc† i ,|Ψ U ⟩=a † 1a† 2 · · ·a † Np |vac⟩,(A12) The corresponding annihilation operator isa µ =...
-
[3]
QR stabilization The formula above is exact, but direct evaluation can become numerically unstable under non-unitary evolu- tion. The columns ofU(t) may develop very different norms and the overlap matrixM=U †Umay become ill-conditioned. QR stabilization removes this problem without changing the correlation matrix. The key observation is that the projecto...
-
[4]
N. Hatano and D. R. Nelson, Localization transitions in non-hermitian quantum mechanics, Phys. Rev. Lett.77, 570 (1996)
work page 1996
- [6]
-
[7]
K. Yokomizo and S. Murakami, Non-bloch band theory of non-hermitian systems, Phys. Rev. Lett.123, 066404 (2019)
work page 2019
- [8]
-
[9]
L. Xiao, T. Deng, K. Wang,et al., Non-Hermitian bulk–boundary correspondence in quantum dynamics, Nat. Phys.16, 761 (2020)
work page 2020
-
[10]
K. Kawabata, T. Numasawa, and S. Ryu, Entanglement phase transition induced by the non-hermitian skin effect, Phys. Rev. X13, 021007 (2023)
work page 2023
- [11]
- [12]
-
[13]
M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. Sa- lomon, Bloch oscillations of atoms in an optical potential, Phys. Rev. Lett.76, 4508 (1996)
work page 1996
-
[14]
S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Q. Niu, and M. G. Raizen, Observation of atomic wannier-stark ladders in an accelerating optical potential, Phys. Rev. Lett.76, 4512 (1996)
work page 1996
-
[15]
A. Chakrabarty and S. Datta, Fate of wannier-stark localization and skin effect in periodically driven non- hermitian quasiperiodic lattices, Phys. Rev. B111, 174202 (2025)
work page 2025
-
[16]
M. Gl¨ uck, M. Hankel, A. R. Kolovsky, and H. J. Korsch, Wannier-stark ladders in driven optical lattices, Phys. Rev. A61, 061402 (2000)
work page 2000
-
[17]
X. Y. Guo, Z. Y. Ge, H. Li,et al., Observation of Bloch oscillations and Wannier-Stark localization on a super- conducting quantum processor, npj Quantum Inf7, 51 (2021)
work page 2021
-
[18]
A. H. Karamlou, J. Braum¨ uller, Y. Yanay,et al., Quan- tum transport and localization in 1d and 2d tight-binding lattices, npj Quantum Inf8, 35 (2022)
work page 2022
-
[19]
Y. Peng, J. Jie, D. Yu, and Y. Wang, Manipulating the non-hermitian skin effect via electric fields, Phys. Rev. B 106, L161402 (2022)
work page 2022
- [20]
-
[21]
S. Ghosh and S. Roy, Exploring topological and localiza- tion phenomena in ssh chains under generalized aah mod- ulation: A computational approach, arXiv:2506.10195 . 11
-
[22]
H.-Z. Li, M. Wan, and J.-X. Zhong, Fate of non-hermitian free fermions with wannier-stark ladder, Phys. Rev. B 110, 094310 (2024)
work page 2024
-
[23]
M. Szyniszewski, O. Lunt, and A. Pal, Disordered moni- tored free fermions, Phys. Rev. B108, 165126 (2023)
work page 2023
- [24]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.