Structural constraints on mobility edges in one-dimensional quasiperiodic systems
Pith reviewed 2026-05-21 15:52 UTC · model grok-4.3
The pith
Mobility edge positions in quasiperiodic systems are fixed by an exact identity linking Lyapunov exponents across isospectral duals.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the bichromatic Aubry-André model and its isospectral dual, the Thouless formula implies an exact identity between the Lyapunov exponents at dual energies. Consequently the positions at which the Lyapunov exponent changes sign are constrained to a reduced set of energies. In the self-dual limit these positions coincide at one critical energy, and the physical Lyapunov spectrum displays linear critical scaling with the standard Aubry-André exponent ν = 1 together with a non-universal energy-dependent prefactor.
What carries the argument
the exact identity for Lyapunov exponents derived from the Thouless formula that connects energies in the bichromatic Aubry-André model to its isospectral dual
If this is right
- Mobility edge positions are restricted to a reduced set of energies rather than appearing independently.
- In the self-dual limit the mobility edges coincide at a single localization-delocalization transition.
- The physical Lyapunov spectrum obeys linear critical scaling near the self-dual point.
- The critical exponent remains the standard Aubry-André value ν = 1, while the prefactor becomes energy-dependent and non-universal.
Where Pith is reading between the lines
- The same duality-based identity may constrain mobility edges in other one-dimensional quasiperiodic models that admit isospectral duals.
- The identity offers a route to locate mobility edges by solving only one member of each dual pair.
- The structural constraint suggests that apparent fine-tuning of mobility edges in quasiperiodic systems is often an artifact of ignoring duality relations.
Load-bearing premise
The bichromatic Aubry-André model and its isospectral dual share exactly the same spectrum, allowing the Thouless formula to produce a direct identity between Lyapunov exponents at corresponding energies.
What would settle it
Numerical computation of the Lyapunov exponent at an energy in the model whose dual partner lies on the opposite side of the predicted mobility edge, showing that the exponents do not satisfy the identity.
Figures
read the original abstract
Mobility edges commonly arise in one-dimensional quasiperiodic systems once exact self-duality is broken, yet their origin is typically understood only at the level of individual Hamiltonians. Here we show that mobility edge positions are not independent spectral features of individual Hamiltonians, but are structurally constrained across quasiperiodic Hamiltonians related by an isospectral duality. Using a bichromatic Aubry--Andr\'e model as a minimal setting, we demonstrate that this constraint is encoded in an exact identity for Lyapunov exponents derived from the Thouless formula. As a consequence, the mobility edge positions are restricted to a reduced set of energies. In the self-dual limit, these mobility edge positions coincide at a single localization--delocalization transition. This structural constraint enforces a linear critical scaling of the physical Lyapunov spectrum near the self-dual point. Numerical results confirm a critical exponent consistent with the standard Aubry--Andr\'e value of $\nu = 1$, while simultaneously revealing a novel, non-universal energy-dependent prefactor.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that mobility edge positions in one-dimensional quasiperiodic systems are not independent features but are structurally constrained across Hamiltonians related by isospectral duality. In the bichromatic Aubry-André model and its dual, this constraint is encoded in an exact identity for Lyapunov exponents obtained from the Thouless formula, which restricts mobility edges to a reduced set of energies. In the self-dual limit the edges coincide at a single transition, enforcing linear critical scaling of the Lyapunov spectrum with the standard exponent ν=1 together with a non-universal energy-dependent prefactor; the latter is confirmed numerically.
Significance. If the identity is rigorously established, the work supplies a structural explanation for mobility-edge locations that goes beyond single-Hamiltonian analysis and rests on standard tools (Thouless formula plus duality). The exact character of the derivation, the recovery of the universal ν=1 exponent, and the identification of an energy-dependent prefactor constitute clear strengths. The result may usefully constrain mobility-edge studies in other quasiperiodic models.
major comments (1)
- [§3] §3 (derivation of the identity): the claim that the Thouless formula directly yields an identity linking Lyapunov exponents at corresponding energies relies on the dual pair sharing the identical integrated density of states. Please state explicitly which equation or step establishes that the integrated density of states is identical for the bichromatic model and its dual, and confirm that no model-specific cancellations are required.
minor comments (2)
- [Abstract] Abstract: the phrase 'the standard Aubry--André value of ν = 1' would benefit from a parenthetical citation to the original Aubry-André work for readers unfamiliar with the literature.
- [Numerical results] Numerical section: the energy-dependent prefactor is presented as novel; ensure that the scaling plots (presumably Fig. 4 or equivalent) explicitly label the distinct energies at which the prefactor is extracted so that the non-universality is immediately visible.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the positive assessment. We address the major comment on the derivation in §3 below.
read point-by-point responses
-
Referee: [§3] §3 (derivation of the identity): the claim that the Thouless formula directly yields an identity linking Lyapunov exponents at corresponding energies relies on the dual pair sharing the identical integrated density of states. Please state explicitly which equation or step establishes that the integrated density of states is identical for the bichromatic model and its dual, and confirm that no model-specific cancellations are required.
Authors: The integrated density of states is identical for the bichromatic Aubry-André model and its dual because the duality is constructed to be isospectral: the transformation maps the original Hamiltonian onto a partner with precisely the same spectrum. This property is established by the definition of the dual Hamiltonian (Eq. (2) in Section 2) and the explicit verification that the duality preserves the eigenvalues. The Thouless formula is then applied to this common integrated density of states, directly producing the Lyapunov-exponent identity without any additional cancellations that would be specific to the bichromatic potential. We will revise the text in §3 to include an explicit cross-reference to Eq. (2) and a short sentence confirming the isospectral character of the duality. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper's central derivation applies the standard Thouless formula to an isospectral dual pair of bichromatic Aubry-André models that share the integrated density of states, with the duality supplying an algebraic relation between transfer-matrix growth rates at corresponding energies. This produces an exact identity for Lyapunov exponents whose zero crossings constrain mobility-edge locations, without any reduction to fitted parameters, self-definitional loops, or load-bearing self-citations. The resulting linear scaling and critical exponent follow directly from the external Thouless relation and duality, rendering the argument self-contained against independent benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math The Thouless formula relates the integrated density of states to the Lyapunov exponent in one-dimensional systems.
- domain assumption The bichromatic Aubry-André model admits an isospectral duality that preserves the spectrum while relating localization properties.
Forward citations
Cited by 2 Pith papers
-
Mobility-edge-embedded Hofstadter butterfly from a tilt-induced quasiperiodic potential
Tilt-induced quasiperiodic potentials on square lattices generate a mobility-edge-embedded Hofstadter butterfly with fractal dimension 0.8-1.0, combining fractal spectra from 1D-like behavior with mobility edges from ...
-
Mobility-edge-embedded Hofstadter butterfly from a tilt-induced quasiperiodic potential
Tilt-induced quasiperiodic potential on a square lattice produces a mobility-edge-embedded Hofstadter butterfly with fractal dimension 0.8-1.0.
Reference graph
Works this paper leans on
-
[1]
E. Abrahams, P. Anderson, D. Licciardello, and T. Ra- makrishnan, Phys. Rev. Lett.42, 673 (1979)
work page 1979
- [2]
- [3]
- [4]
- [5]
- [6]
- [7]
-
[8]
M. Saha, S. K. Maiti, and A. Purkayastha, Phys. Rev. B 100, 174201 (2019)
work page 2019
-
[9]
X. Deng, S. Ray, S. Sinha, G. V. Shlyapnikov, and L. San- tos, Phys. Rev. Lett.123, 025301 (2019)
work page 2019
- [10]
- [11]
- [12]
-
[13]
S. Ganeshan, J. H. Pixley, and S. Das Sarma, Phys. Rev. Lett.114, 146601 (2015)
work page 2015
- [14]
-
[15]
M. Gon¸ calves, B. Amorim, E. V. Castro, and P. Ribeiro, Phys. Rev. B108, L100201 (2023)
work page 2023
-
[16]
C. M. Soukoulis and E. N. Economou, Phys. Rev. Lett. 48, 1043 (1982)
work page 1982
- [17]
- [18]
-
[19]
D. J. Boers, B. Goedeke, D. Hinrichs, and M. Holthaus, Phys. Rev. A75, 063404 (2007)
work page 2007
-
[20]
X. Li, X. Li, and S. Das Sarma, Phys. Rev. B96, 085119 (2017)
work page 2017
- [21]
-
[22]
M. Gon¸ calves, B. Amorim, E. V. Castro, and P. Ribeiro, Phys. Rev. Lett.131, 186303 (2023)
work page 2023
- [23]
- [24]
-
[25]
O. I. Utesov, Y. Kim, and S. Flach, Phys. Rev. B112, 174207 (2025)
work page 2025
-
[26]
S. Iyer, V. Oganesyan, G. Refael, and D. A. Huse, Phys. Rev. B87, 134202 (2013)
work page 2013
- [27]
-
[28]
M. ˇZnidariˇ c and A. Ljubotina, Proc. Natl. Acad. Sci. 115, 5345 (2018)
work page 2018
-
[29]
S. Xu, X. Li, Y.-T. Hsu, B. Swingle, and S. Das Sarma, Phys. Rev. Res.1, 032039 (2019)
work page 2019
-
[30]
E. V. H. Doggen and A. D. Mirlin, Phys. Rev. B100, 104203 (2019)
work page 2019
-
[31]
A. S. Aramthottil, T. Chanda, P. Sierant, and J. Za- krzewski, Phys. Rev. B104, 214201 (2021)
work page 2021
-
[32]
D. D. Vu, K. Huang, X. Li, and S. Das Sarma, Phys. Rev. Lett.128, 146601 (2022)
work page 2022
-
[33]
M. Gon¸ calves, B. Amorim, F. Riche, E. V. Castro, and P. Ribeiro, Nat. Phys.20, 1933 (2024)
work page 1933
-
[34]
J. D. Bodyfelt, D. Leykam, C. Danieli, X. Yu, and S. Flach, Phys. Rev. Lett.113, 236403 (2014)
work page 2014
- [35]
-
[36]
S. Lee, A. Andreanov, and S. Flach, Phys. Rev. B107, 014204 (2023)
work page 2023
-
[37]
S. Lee, S. Flach, and A. Andreanov, Chaos: An Interdis- ciplinary Journal of Nonlinear Science33, 073125 (2023)
work page 2023
- [38]
-
[39]
B. Deissler, M. Zaccanti, G. Roati, C. D’Errico, M. Fat- tori, M. Modugno, G. Modugno, and M. Inguscio, Nat. Phys.6, 354 (2010)
work page 2010
-
[40]
M. Schreiber, S. S. Hodgman, P. Bordia, H. P. L¨ uschen, M. H. Fischer, R. Vosk, E. Altman, U. Schneider, and I. Bloch, Science349, 842 (2015)
work page 2015
- [41]
-
[42]
F. A. An, E. J. Meier, and B. Gadway, Phys. Rev. X8, 031045 (2018)
work page 2018
-
[43]
F. A. An, K. Padavi´ c, S. Scherg, A. Janot, and B. Gad- way, Phys. Rev. Lett.126, 040603 (2021)
work page 2021
-
[44]
H. P. L¨ uschen, S. Scherg, T. Kohlert, M. Schreiber, P. Bordia, X. Li, S. Das Sarma, and I. Bloch, Phys. Rev. Lett.120, 160404 (2018)
work page 2018
-
[45]
T. Shimasaki, M. Prichard, L. Cao, and D. M. Weld, Nat. Phys.20, 409 (2024)
work page 2024
- [46]
-
[47]
Y. E. Kraus, Y. Lahini, Z. Ringel, M. Verbin, and O. Zil- berberg, Phys. Rev. Lett.109, 106402 (2012)
work page 2012
- [48]
-
[49]
P. Wang, Y. Zheng, X. Chen, C. Huang, Y. V. Kartashov, L. Torner, V. V. Konotop, and F. Ye, Nature577, 42 (2020)
work page 2020
-
[50]
X. Ni, K. Chen, M. Weiner, D. J. Apigo, F. Schindler, J. Elmer, A. D. Greentree, M. Hafezi, N. Malkova, and M. C. Rechtsman, Commun. Phys.2, 55 (2019)
work page 2019
-
[51]
D. J. Apigo, W. Cheng, K. F. Dobiszewski, E. Prodan, and C. Prodan, Phys. Rev. Lett.122, 095501 (2019)
work page 2019
- [52]
-
[53]
V. I. Oseledets, Trudy Moskovskogo Matematicheskogo Obshchestva19, 179 (1968)
work page 1968
- [54]
-
[55]
D. J. Thouless, J. Phys. C: Solid State Phys.5, 77 (1972)
work page 1972
-
[56]
J. Chapman and G. Stolz, inAnnales Henri Poincar´ e, Vol. 16 (Springer, 2015) pp. 405–435
work page 2015
-
[57]
F. V. Oliveira and S. L. Carvalho, Operators & Matrices 16(2022)
work page 2022
- [58]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.