REVIEW 4 major objections 5 minor 35 references
Complex Band Structure and localisation transition for tridiagonal non-Hermitian k-Toeplitz operators with defects
T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read The paper shows that the complex band structure of tridiagonal k-Toeplitz operators yields explicit exponential decay rates for all localised eigenvectors, with a single arccosh formula separating skin localisation, bulk localisation, and…
desk verdict Useful defect-mode decay formulas, but the Toeplitz-operator eigenvector results are broken by the missing first-row boundary condition. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The key object is the complex band structure of the symbol $f(z)=A_{-1}z^{-1}+A_0+A_1 z$ of the $k$-Toeplitz operator, organised by the determinant identity $\det(f(e^{-i(\alpha+i\beta)})-\lambda I)=A e^{-i(\alpha+i\beta)}+B e^{i(\alpha+i\beta)}+g(\lambda)$. After writing $\beta=r+\tilde{\beta}$, the reality condition becomes $2A e^r \cos(\alpha)\cosh(\tilde{\beta})+g(\lambda)=0$, so band functions have $\tilde{\beta}=0$ and gap functions have $\alpha\in\{0,\pi\}$ with $\tilde{\beta}(\lambda)=\operatorname{arccosh}\left(-g(\lambda)/(2A e^r)\right)$. The argument works by gluing together two quasiperiodic eigenvector extensions, one for each branch; the imaginary part of the quasimomentum then appears directly as the exponential decay rate of the eigenvector entries.
What would settle it
For a finite monomer resonator chain with a single defect and known parameters, measure the defect eigenmode on both sides of the defect and fit exponentials with rates $r-\tilde{\beta}$ and $r+\tilde{\beta}$, where $\tilde{\beta}=\operatorname{arccosh}\left((a-\omega^2)/(2\sqrt{bc})\right)$; if the fitted rates do not sum to $2r$ or differ from the predicted values beyond numerical error, the central claim is false. A second check is to compute the eigenvalues of a large finite $k$-Toeplitz matrix inside the winding region and test whether every one lies in the union of symbol eigenvalues at fixed $\beta=r$, as Theorem 2.1 requires.
Extended reading notes
Core claim
The central claim is that, for non-Hermitian tridiagonal $k$-Toeplitz operators and their defected counterparts, the exponential decay of every localised eigenvector is determined by the complex band structure. Inside the open-limit spectrum the eigenvector entries satisfy $|u(i+k)|/|u(i)| = e^{-r}$; outside the open spectrum but inside the winding region the rate is $O(e^{-(r-\tilde{\beta}(\lambda))})$; and for a defected Laurent operator with $\lambda$ outside the winding region the rates are $e^{-(r-\tilde{\beta}(\lambda))}$ leading up to the defect and $e^{-(r+\tilde{\beta}(\lambda))}$ past it, with $\tilde{\beta}(\lambda)=\operatorname{arccosh}\left(-g(\lambda)/(2A e^r)\right)$. The proof constructs the eigenvectors as linear combinations of quasiperiodic extensions of the symbol's eigenvectors, and the same function $\tilde{\beta}$ controls convergence of finite-matrix defect eigenvalues and the accuracy of truncated eigenvectors as pseudoeigenvectors. This yields a quantitative account of skin localisation, bulk localisation, and the transition between them in non-Hermitian resonator chains and tight-binding models.
Load-bearing premise
The framework stands on two imported ingredients: the generalised Brillouin zone theorem that identifies the open-limit spectrum of the finite matrices, and a uniform bound on the condition number of the eigenbasis; if either is false for non-Hermitian $k$-Toeplitz operators, the predicted decay rates and finite-size convergence estimates do not follow.
Editorial extensions
If this is right
- Inside the open-limit spectrum, all eigenmodes decay at the uniform rate $e^{-r}$, fixing the skin-effect decay length from the coupling coefficients alone.
- Defect modes at frequencies outside the winding region are exponentially localised around the defect with asymmetric rates $r-\tilde{\beta}$ and $r+\tilde{\beta}$ on the two sides, giving a quantitative bulk localisation length.
- For finite truncated chains, truncated exact eigenvectors are exponentially good pseudoeigenvectors with error $O(e^{-B\lfloor N/2k\rfloor})$ for $B=r-\tilde{\beta}$ or $B=\tilde{\beta}-r$, so the infinite-system predictions transfer to realistic finite systems.
- Defect eigenfrequencies of finite matrices converge to their infinite-limit values with error $O(e^{-\tilde{\beta} N})$.
- In the Hatano-Nelson tight-binding Hamiltonian, the skin-to-bulk localisation transition at defect strength $d=\pm 2\sinh\gamma$ is explained by the gap branch of the complex band structure.
Reading between the lines
- A directly testable signature of the asymmetric decay law is the amplitude ratio $|u(m+q)|/|u(m-q)|$ around a single-site defect, which the paper's rates predict grows like $e^{2\tilde{\beta} q}$; an experiment on a defected resonator chain could measure this ratio.
- Because $\tilde{\beta}(\lambda)$ depends on the eigenvalue only through the scalar polynomial $g(\lambda)$, the same formula should apply to any compact defect entering the diagonal of the $k$-Toeplitz matrix, not just the single-site defects treated in the paper.
- The rate formula suggests an inverse-design route: choosing a defect parameter that puts a target frequency at a prescribed gap-branch value $\tilde{\beta}$ fixes the localisation length without solving the finite eigenproblem.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a complex band structure method for tridiagonal k-Toeplitz operators, constructing eigenvectors as linear combinations of quasiperiodic Bloch modes with complex quasimomentum. It claims sharp exponential decay rates for eigenvectors in the winding region outside the open limit spectrum, for defect modes of Laurent operators, and for finite truncations via pseudoeigenvectors. The results are applied to non-Hermitian resonator chains and Hatano–Nelson tight-binding models, with numerical confirmation and openly available Matlab code.
Significance. If established rigorously, the paper would give an explicit, symbol-only formula for localisation lengths of eigenvectors in a broad class of non-Hermitian tridiagonal problems, going beyond qualitative winding-number arguments. The combination of analytic decay rates, explicit pseudoeigenvector bounds, and numerical validation in physical models is appealing, and the open code is a strength. However, the central derivations in Section 2 currently leave key boundary-condition and convergence arguments unchecked, so the main quantitative claims are not yet fully supported.
major comments (4)
- [§2.1, Corollary 2.8] The proof of Corollary 2.8 does not verify that the constructed linear combination satisfies the first-row boundary condition of the semi-infinite Toeplitz operator, nor does it establish existence of the coefficients a,b. I checked the scalar case k=1: for λ∈σwind\σopen the two characteristic roots satisfy ρ_+ρ_- = c/b < 1, and the first-row equation reduces to Aρ_- + Bρ_+ = 0, which has nonzero real solutions, so an ℓ² eigenvector does exist. Thus the specific objection that no such eigenvector exists is not valid. Nevertheless, the manuscript should state this verification; for general k the boundary condition reduces to a scalar equation only because the off-diagonal blocks A_±1 in a tridiagonal k-Toeplitz matrix have rank one, and this structural fact is neither mentioned nor used.
- [§2.2, Corollary 2.12 and Theorem 2.13] The residual estimate in Theorem 2.13 is not derived rigorously: the expression in (2.42) contains undefined or spurious terms, and the proof does not specify precisely where the truncated infinite-operator eigenvector fails to solve the finite matrix equation. The bound on the tail sum in (2.43) also depends on applying (2.29) to a vector whose existence is not fully established. Corollary 2.12 similarly jumps from a pseudospectral inclusion to a convergence rate without a detailed argument that the truncated defect eigenvector has exponentially small residual; it also relies on an unproved uniform condition-number bound κ(V)≤C from [3, Theorem 3.5], which is load-bearing for the exponential rate.
- [§2.2, Corollary 2.12] The statement of Corollary 2.12 refers to 'the defect eigenfrequency of the Toeplitz operator T(f)', but T(f) was defined in Section 2 as an undefected Toeplitz operator and the paper does not define a defected Toeplitz operator. The manuscript later defines a defected Laurent operator L(f) and uses it in Corollary 2.9, so the convergence result should be stated for the defected operator actually under study, with matching notation.
- [§2.1, Corollary 2.9] The proof of Corollary 2.9 constructs a candidate eigenvector by splicing the left half of u1 and the right half of u2 but does not check the defect-row equation. For k=1 this can be verified directly, but for general k the defect row is a k-dimensional condition with only two free coefficients; the argument must show that the rank-one structure of the defect makes the system consistent. As written, the proof does not establish that the proposed u is an exact eigenvector of the defected Laurent operator.
minor comments (5)
- [§2.1, Theorem 2.3 proof] In the proof of Theorem 2.3, the statement 'it must hold that α={0,π}' should read 'α∈{0,π}'; the notation is otherwise misleading.
- [§2.1, Corollary 2.8 proof] The phrase 'in the poof of Theorem 2.6' is a typo; it should refer to the proof of Corollary 2.6.
- [§2.2, Theorem 2.11] The reference to Bauer–Fike is written as '[31? , Theorem 2.3]' with a stray question mark; the citation should be cleaned up.
- [Figure 2.1 caption] The caption of Figure 2.1 has garbled text ('!:=L 0 :=L' and '·;-.') that should be corrected.
- [§3.4.1, Theorem 3.11] The name 'Chebychev' is spelled inconsistently; standard spelling is 'Chebyshev'.
Circularity Check
No significant circularity: decay rates are computed from the symbol and cross-checked by explicit inverses; the self-cited prior theorems are external evidence, not definitional inputs.
full rationale
The claimed derivation is not circular. The complex band functions are introduced in Definition 2.2, and the constraint alpha in {0,pi} or beta = r is proved in Theorem 2.3 from the characteristic expansion (2.9); the gap decay parameter beta~(lambda) = arccosh(-g(lambda)/(2 A e^r)) is obtained by solving (2.16), not by fitting eigenvector data. Corollary 2.8 and Corollary 2.9 construct eigenvectors as linear combinations of the two quasiperiodic modes; the decay rates are immediate from the mode exponents, but the first-row/defect-row matching is a genuine additional condition. For scalar k = 1, the matching is possible: with z_+ = sqrt(c/b) e^{beta~} and z_- = sqrt(c/b) e^{-beta~}, choosing coefficients proportional to (z_+, -z_-) satisfies the first-row equation since z_+ + z_- = (lambda - a)/b, and the vector is square-summable with asymptotic ratio z_+ = e^{-(r - beta~)}. Thus the decay bound is not defined into existence. The same rates are independently recovered from the explicit Chebyshev inverse formula in Corollaries 3.12-3.13, confirming that the band-structure prediction has content beyond its ansatz. The paper does rely on prior theorems, notably [4, Theorem 4.1] for the generalized Brillouin zone and [3, Lemma 2.6 and Theorem 3.5] for the determinant expansion and uniform condition number bound; [3] has an overlapping author, but these are published external theorems with stated assumptions, not merely a self-citation chain, so under the rubric they are real evidence rather than circularity. No parameter is fitted and renamed as a prediction, and the central finite-chain defect-mode results are corroborated numerically against independently computed Green's-function decays.
Assumptions & free parameters
assumptions (5)
- domain assumption b_i c_i > 0 for all 1≤i≤k
- standard math Generalized Brillouin zone theorem (open-limit spectrum equals union over α of σ(f(e^{-i(α+ir)})))
- standard math Uniform condition number bound κ(V)≤C for the eigenbasis of T_N(f(e^r·))
- standard math Confluent eigenvectors can be recovered by the method of [3, Theorem 2.10]
- domain assumption Subwavelength resonances converge to eigenvalues of the quasiperiodic capacitance matrix as δ→0
Cite this review
Pith. "Pith review of Complex Band Structure and localisation transition for tridiagonal non-Hermitian k-Toeplitz operators with defects." pith.science (2026). https://pith.science/paper/YL3FBLME
@misc{pith2026250523610,
author = {Pith},
title = {Pith review of: Complex Band Structure and localisation transition for tridiagonal non-Hermitian k-Toeplitz operators with defects},
year = {2026},
howpublished = {\url{https://pith.science/paper/YL3FBLME}},
note = {Machine review of arXiv:2505.23610}
}
read the original abstract
Using the Bloch-Floquet theory, we propose an innovative technique to obtain the eigenvectors of tridiagonal k-Toeplitz operators. This method offers a more extensive and quantitative basis for describing localised eigenvectors beyond the non-trivial winding zone, yielding sharp decay bounds. The validity of our results is confirmed numerically in one-dimensional resonator chains, showcasing non-Hermitian skin localisation, bulk localisation, and tunnelling effects. We conclude the paper by analysing non-Hermitian tight binding Hamiltonians, illustrating the broad applicability of the complex band structure.
Figures
Figures from the paper (11 more)
Reference graph
Works this paper leans on
-
[1]
H. B. Al Ba’ba’a. Brillouin-zone definition in non-reciprocal willis monatomic lattices.JASA Express Letters, 3(12):120001, 12 2023
work page 2023
- [2]
- [3]
- [4]
- [5]
- [6]
- [7]
- [8]
Show all 35 references
-
[9]
C. M. Bender and S. Boettcher. Real spectra in non-hermitian hamiltonians having pt symmetry. Phys. Rev. Lett., 80:5243–5246, 6 1998
1998
-
[10]
D. S. Borgnia, A. J. Kruchkov, and R.-J. Slager. Non-hermitian boundary modes and topology.Phys. Rev. Lett., 124:056802, 2 2020
2020
-
[11]
Böttcher and B
A. Böttcher and B. Silbermann.Introduction to Large Truncated Toeplitz Matrices. Universitext. Springer-Verlag, New York
-
[12]
P. W. Brouwer, P. G. Silvestrov, and C. W. J. Beenakker. Theory of directed localization in one dimension. Phys. Rev. B, 56:R4333–R4335, 8 1997
1997
-
[13]
Y. D. Bruijn and E. O. Hiltunen. Complex brillouin zone for localised modes in hermitian and non-hermitian problems, 2025
2025
-
[14]
da Fonseca and J
C. da Fonseca and J. Petronilho. Explicit inverse of a tridiagonal k-toeplitz matrix.Numerische Mathematik, 100(3):457–482, 2005
2005
-
[15]
Davies, S
B. Davies, S. Barandun, E. O. Hiltunen, R. V. Craster, and H. Ammari. Two-scale effective model for defect-induced localization transitions in non-hermitian systems.Phys. Rev. B, 111:035109, 2025
2025
-
[16]
De Bruijn and E
Y. De Bruijn and E. O. Hiltunen. Complex band structure for subwavelength evanescent waves.Studies in Applied Mathematics, 154(2):e70022, 2025
2025
-
[17]
Feppon, Z
F. Feppon, Z. Cheng, and H. Ammari. Subwavelength resonances in one-dimensional high-contrast acoustic media. 83(2):625–665
-
[18]
F. M. Fernández. On a class of non-hermitian hamiltonians with tridiagonal matrix representation. Annals of Physics, 443:169008, Aug. 2022
2022
-
[19]
Ghatak, M
A. Ghatak, M. Brandenbourger, J. Van Wezel, and C. Coulais. Observation of non-Hermitian topology and its bulk–edge correspondence in an active mechanical metamaterial. 117(47):29561–29568
-
[20]
Gohberg, M
I. Gohberg, M. A. Kaashoek, and S. Goldberg.Classes of Linear Operators. Birkhäuser Basel, 1993
1993
-
[21]
I. S. Gradshteyn and I. M. Ryzhik.Table of Integrals, Series, and Products. Academic Press, San Diego, CA, 7th edition, 2014
2014
-
[22]
Hatano and D
N. Hatano and D. R. Nelson. Localization transitions in non-hermitian quantum mechanics.Phys. Rev. Lett., 77:570–573, 7 1996
1996
-
[23]
Hatano and D
N. Hatano and D. R. Nelson. Non-hermitian delocalization and eigenfunctions. Phys. Rev. B, 58:8384–8390, 10 1998
1998
-
[24]
Jana and L
S. Jana and L. Sirota. Invisible tunneling through non-hermitian barriers in nonreciprocal lattices. Phys. Rev. B, 111:L100301, 2025
2025
-
[25]
S. Longhi. Probing non-hermitian skin effect and non-bloch phase transitions. Phys. Rev. Res., 1:023013, 9 2019
2019
-
[26]
S. Longhi. Non-hermitian skin effect beyond the tight-binding models.Phys. Rev. B, 104:125109, 2021
2021
-
[27]
R. Nabben. Decay rates of the inverse of nonsymmetric tridiagonal and band matrices.SIAM Journal on Matrix Analysis and Applications, 20(3):820–837, 1999
1999
-
[28]
Okuma, K
N. Okuma, K. Kawabata, K. Shiozaki, and M. Sato. Topological origin of non-hermitian skin effects. Phys. Rev. Lett., 124:086801, 3 2020
2020
-
[29]
Okuma and M
N. Okuma and M. Sato. Non-hermitian topological phenomena: A review.Annual Review of Condensed Matter Physics, 14:83–107, 2023
2023
-
[30]
Reichel and L
L. Reichel and L. N. Trefethen. Eigenvalues and pseudo-eigenvalues of toeplitz matrices.Linear Algebra and its Applications, 162-164:153–185, 1992
1992
-
[31]
L. N. Trefethen and M. Embree.Spectra and Pseudospectra. Princeton University Press, Princeton, NJ
-
[32]
Z. Yang, K. Zhang, C. Fang, and J. Hu. Non-hermitian bulk-boundary correspondence and auxiliary generalized brillouin zone theory.Phys. Rev. Lett., 125:226402, 11 2020
2020
-
[33]
Yao and Z
S. Yao and Z. Wang. Edge states and topological invariants of non-hermitian systems.Phys. Rev. Lett., 121:086803, 8 2018
2018
-
[34]
Zhang, Z
K. Zhang, Z. Yang, and C. Fang. Correspondence between winding numbers and skin modes in non-hermitian systems. Phys. Rev. Lett., 125:126402, 9 2020. 20 Y. De Bruijn and E. O. Hiltunen
2020
-
[35]
Zhang, X
X. Zhang, X. Song, S. Zhang, T. Zhang, Y. Liao, X. Cai, and J. Li. Solvable non-hermitian skin effects and real-space exceptional points: non-hermitian generalized bloch theorem.Journal of Physics A: Mathematical and Theoretical, 57(12):125001, 3 2024. Yannick De Bruijn Depart...
2024
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.