REVIEW 3 major objections 5 minor 57 references
Mixed spectra and partially extended states in a two-dimensional quasiperiodic model
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper introduces a two-dimensional Aubry-André model in which localised and partially extended single-particle states coexist at any modulation strength, with no mobility edge separating them.
desk verdict A genuinely new 2D quasiperiodic model with a plausible mixed-spectrum mechanism; the core claim is real but the asymptotic extensivity of the partially extended states is under-supported and needs a finite-size scaling test. 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 load-bearing object is the factorisation identity $V_{nm} = 2\lambda J \cos(2\pi\beta n)\cos(2\pi\beta m) = \tilde{\lambda}_n J \cos(2\pi\beta m)$. It turns the two-dimensional model, when horizontal hopping is switched off, into a stack of independent one-dimensional Aubry-André chains, one per vertical line, with line-dependent disorder amplitudes $\tilde{\lambda}_n = 2\lambda \cos(2\pi\beta n)$. The one-dimensional critical value $|\tilde{\lambda}|=2$ marks which lines are extended and which are localised. Restoring horizontal hopping hybridises only the extended lines into a network, while the exponential localisation of the other lines acts as a shield that prevents hybridisation between the two species. The Aubry duality of the model, a Fourier transform combined with a 45-degree rotation, maps the strong-modulation line network onto a diagonal network at weak modulation.
What would settle it
Compute the eigenstates of the 2DAA model at $\lambda=40$ on lattices up to $L\approx 1000$ and check whether all states with large participation ratio are concentrated on the lines with $|2\lambda\cos(2\pi\beta n)|<2$; finding an extended state with weight away from those lines, or finding that localised states hybridise into extended ones as horizontal hopping is switched on continuously, would refute the claim. In an experiment, a site placed just off a low-disorder line should stay exponentially localised; seeing non-exponentially small probability reach the network on any timescale would also falsify the shielding picture.
Extended reading notes
Core claim
The central discovery is that the two-dimensional Aubry-André Hamiltonian is self-dual like its one-dimensional parent, yet its spectrum is much richer. The potential factorises as $V_{nm} = 2\lambda J \cos(2\pi\beta n) \cos(2\pi\beta m) = \tilde{\lambda}_n J \cos(2\pi\beta m)$, so with horizontal hopping removed each vertical lattice line is an independent one-dimensional Aubry-André chain with effective disorder amplitude $\tilde{\lambda}_n = 2\lambda \cos(2\pi\beta n)$. Because $\tilde{\lambda}_n$ is itself quasiperiodic, lines with $|\tilde{\lambda}_n| < 2$ exist for any $\lambda$, and these lines carry extended eigenstates that hybridise into a mesh of partially extended states. Lines with $|\tilde{\lambda}_n| > 2$ stay exponentially localised, and that exponential localisation shields them from the delocalising mesh, so localised and partially extended states coexist interspersed in energy with no mobility edge. The same structure yields ballistic expansion along the low-disorder network while generic initial sites remain localised.
Load-bearing premise
The central claim rests on the assumption that restoring the horizontal hopping couples only the already extended lines and cannot turn exponentially localised segments on the other lines into extended channels; if that shielding fails, the mixed-spectrum effect collapses.
Editorial extensions
If this is right
- At any $\lambda>2$, a finite fraction roughly $2/(\pi\lambda)$ of lattice lines are in the extended one-dimensional phase, so partially extended states exist for arbitrarily strong quasiperiodic modulation.
- Localised and partially extended states occur at the same energies, interspersed in the spectrum, so no mobility edge forms.
- A particle released on a low-disorder line expands ballistically along the network of such lines, while a particle released elsewhere remains exponentially localised under the same parameters.
- The effect is robust to weak random disorder and to small tilts of the modulation away from 45 degrees, and the ground-state localisation transition belongs to the same universality class as a continuum eightfold-symmetric model.
Reading between the lines
- If this line-shielding mechanism is generic, then any two-dimensional quasiperiodic potential that factorises into a product of one-dimensional incommensurate functions should show a similar mixed spectrum; the exact 45-degree geometry is sufficient but probably not necessary.
- The roughly $2/(\pi\lambda)$ fraction of extended lines implies that the number of partially extended states grows linearly with system size in the thermodynamic limit, so their prevalence is extensive even though individual states look quasi-one-dimensional on intermediate scales.
- A clean experimental signature would be higher-order moments of the expansion: Manhattan-distance moments should stay linear longer than Euclidean moments, because reaching a line intersection changes Euclidean distance but not Manhattan distance along the network.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces a two-dimensional, non-separable generalisation of the Aubry-André model (2DAA) and studies its single-particle localisation properties by exact diagonalisation. The central claims are that (i) the model is self-dual under a generalised Aubry transformation, (ii) at strong quasiperiodic modulation λ ≫ 2 a finite fraction of eigenstates remain partially extended on a sparse network of low-disorder lattice lines, while other states are exponentially localised, (iii) these two types of states are interspersed in energy so that no mobility edge forms, and (iv) dynamics initiated on the low-disorder network are ballistic along that network. The mechanism is explained by rewriting the potential as independent 1D Aubry-André chains on each lattice line with effective disorder amplitude λ̃_n = 2λ cos(2πβ n); lines with |λ̃_n| < 2 are extended and are argued to hybridise into a 2D network, while lines with |λ̃_n| > 2 remain localised and shield the extended states. The paper also studies tilted and non-self-dual variants, the effect of random disorder, and the universality of the ground-state transition with a continuum quasiperiodic model.
Significance. If the central claim holds, the paper establishes a new and rather general mechanism for coexistence of localised and partially extended states in a two-dimensional quasiperiodic system, with direct relevance to cold-atom experiments and to discussions of rare regions in many-body localisation. The paper has several concrete strengths: the duality transformation in Appendix A is exact; the line-decomposition argument after Eq. (4) is a clean, falsifiable physical picture; Fig. 3 provides a sharp, checkable prediction that partially extended states are confined to the critical 1DAA bandwidth; and the dynamics in Sec. III and Appendix B give an experimentally accessible signature. The main weakness is that the asymptotic extensivity of the partially extended states (PR ∝ L^2) is an extrapolation from data showing PR exponents only up to about 1.7, and the shielding of localised states from hybridisation is asserted rather than quantitatively tested in the coupled system. These are load-bearing for the paper's strongest claims, so the manuscript needs additional finite-size analysis before the central conclusion is fully established.
major comments (3)
- [Sec. II, paragraph after Eq. (4)] The central claim that partially extended states have extensive participation ratio is not yet established by the data. The text states that for sufficiently large systems PR scales as L^2 because a finite fraction 2/(πλ) of lines is extended, but Fig. 3 (L = 239) shows log PR/log L reaching only about 1.7, and Fig. 1 (L = 99) shows the funnel at values much closer to 1. No finite-size scaling is presented showing that the exponent increases toward 2 as L grows. Please add a finite-size analysis at fixed λ (for example λ = 5, 10, 40 across L = 99, 239, and larger accessible sizes) with confidence intervals; without it, the data are also compatible with a quasi-1D interpretation (PR ∝ L), which would change the physical conclusion.
- [Sec. II, paragraph after Eq. (4)] The shielding assumption that protects localised states from hybridisation is the load-bearing step for the mixed-spectrum and no-mobility-edge conclusions, but it is only argued qualitatively. The paper asserts that exponential localisation of lines with |λ̃| > 2 suppresses their statistical weight and prevents hybridisation with the extended network, yet no calculation in the coupled 2D system tests this. Please provide a direct numerical test, for example the inverse participation ratio of eigenstates whose weight is concentrated on high-disorder lines as a function of L at fixed λ, or the overlap of such states with the low-disorder network; this would either substantiate or refute the assumed protection.
- [Sec. III and Fig. 3] The no-mobility-edge claim is inferred from a single-size spectrum (L = 239 in Fig. 3). At finite size, an energy-dependent PR distribution could also arise from a finite-size mobility edge whose edges sharpen or shift with L. Since the abstract claims that no mobility edge emerges, please show the system-size dependence of the PR distribution within the critical 1DAA bandwidth for several λ values; this would also help separate the finite-size crossover from the proposed asymptotic mixed-spectrum behaviour.
minor comments (5)
- [Sec. II, text around Fig. 1] The text says the funnel of partially extended states has PR ∼ N_sites^{0.5}, while Fig. 3 and its caption quote values L^1 to L^1.7; please make the scaling statement consistent across the two figures.
- [Appendix B] The appendix states that the wave functions are shown 'for λ = 0' in Fig. 8, but the main text and the figure caption describe λ = 40; please correct this typo.
- [Sec. III, Eq. (6)] For degenerate eigenstates, the diagonal ensemble in Eq. (5) and the participation ratio in Eq. (6) need a precise definition of the eigenbasis; please specify how degeneracies are handled.
- [Sec. IV C and Fig. 7] The finite-size scaling claim for the universality of the ground-state transition would be strengthened by reporting the fitting procedure, the number of data points, and goodness-of-fit measures such as χ² or residuals.
- [Sec. IV A] The statement that the non-self-dual 'wide channel' construction is 'manifestly more robust' than the self-dual one is qualitative; a quantitative comparison of PR versus λ at matched system sizes would be more informative.
Circularity Check
No circularity: the 2DAA mechanism is derived from an independent 1DAA benchmark and checked against exact diagonalisation; self-citations are motivational, not load-bearing.
full rationale
The central derivation is self-contained. The partially extended states are identified directly from exact-diagonalisation participation ratios (Eq. 3, Figs. 1-3), and the explanatory mechanism in Sec. II (Eq. 4) is an independent rewriting of the potential plus the textbook 1DAA localisation transition [21], not an input fitted to the 2D spectra. The energy-window test (Fig. 3) is not circular: the critical 1DAA bandwidth [38] is computed independently, while the states plotted are all eigenstates near E=0 and are labelled by their PR, not by their energy being inside the window. The L^2-extensivity statement for asymptotically large systems is an extrapolation from the observed log PR/log L values around 1.7 at L=239; that is a correctness/robustness concern, not a definitional reduction. The universality section invokes the authors' earlier 1D result [44] and their experiment [12] as motivation, but the same-universality-class conclusion rests on the finite-size scaling collapse of the ground-state curvature computed for both models in this paper, with exponents fitted to those data; the self-citations are not load-bearing. No uniqueness theorem, fitted parameter renamed as prediction, or ansatz smuggled via citation was found. The low score reflects only the presence of minor non-load-bearing self-citations and an asymptotic extrapolation, neither of which constitutes circularity.
Assumptions & free parameters
assumptions (4)
- domain assumption The one-dimensional Aubry-André model has a localisation transition at lambda = 2: extended for lambda < 2, localised for lambda > 2.
- domain assumption In two dimensions, any amount of random disorder leads to exponential localisation with exponentially large localisation lengths.
- domain assumption The lattice model (2) and the continuum Hamiltonian (8) are continuously connected so that their ground-state localisation transitions share the same universality class.
- standard math Aubry duality for the 2D model maps lambda to 4/lambda via a Fourier transform on the square lattice.
Cite this review
Pith. "Pith review of Mixed spectra and partially extended states in a two-dimensional quasiperiodic model." pith.science (2026). https://pith.science/paper/IGZZESCJ
@misc{pith2026190902048,
author = {Pith},
title = {Pith review of: Mixed spectra and partially extended states in a two-dimensional quasiperiodic model},
year = {2026},
howpublished = {\url{https://pith.science/paper/IGZZESCJ}},
note = {Machine review of arXiv:1909.02048}
}
read the original abstract
We introduce a two-dimensional generalisation of the quasiperiodic Aubry-Andr\'e model. Even though this model exhibits the same duality relation as the one-dimensional version, its localisation properties are found to be substantially more complex. In particular, partially extended single-particle states appear for arbitrarily strong quasiperiodic modulation. They are concentrated on a network of low-disorder lattice lines, while the rest of the lattice hosts localised states. This spatial separation protects the localised states from delocalisation, so no mobility edge emerges in the spectrum. Instead, localised and partially extended states are interspersed, giving rise to an unusual type of mixed spectrum and enabling complex dynamics even in the absence of interactions. A striking example is ballistic transport across the low-disorder lines while the rest of the system remains localised. This behaviour is robust against disorder and other weak perturbations. Our model is thus directly amenable to experimental studies and promises fascinating many-body localisation properties.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
P. G. Harper, Proc. Phys. Soc. A68, 874 (1955)
1955
-
[2]
M. Ya. Azbel’, Sov. Phys. JETP19, 634 (1964)
work page 1964
-
[3]
D. R. Hofstadter, Phys. Rev. B14, 2239 (1976)
1976
- [4]
-
[5]
Roati, C
G. Roati, C. D’Errico, L. Fallani, M. Fattori, C. Fort, M. Za- ccanti, G. Modugno, M. Modugno, and M. Inguscio, Nature 453, 895 (2008)
2008
-
[6]
B. Deissler, M. Zaccanti, G. Roati, C. D’Errico, M. Fattori, M.Modugno,G.Modugno, andM.Inguscio,Nat.Phys. 6,354 (2010)
work page 2010
- [7]
-
[8]
M. Schreiber, S. S. Hodgman, P. Bordia, H. P. Lüschen, M. H. Fischer,R.Vosk,E.Altman,U.Schneider, andI.Bloch,Science 349, 842 (2015)
work page 2015
Show all 57 references
-
[9]
Schneider, Phys
P.Bordia,H.P.Lüschen,S.S.Hodgman,M.Schreiber,I.Bloch, and U. Schneider, Phys. Rev. Lett.116, 140401 (2016)
2016
-
[10]
Bordia, H
P. Bordia, H. P. Lüschen, S. Scherg, S. Gopalakrishnan, M. Knap, U. Schneider, and I. Bloch, Phys. Rev. X7, 41047 (2017)
2017
-
[11]
H. P. Lüschen, P. Bordia, S. Scherg, F. Alet, E. Altman, U.Schneider, andI.Bloch,Phys.Rev.Lett. 119,260401(2017)
2017
-
[12]
K.Viebahn,M.Sbroscia,E.Carter,J.-C.Yu, andU.Schneider, Phys. Rev. Lett.122, 110404 (2019)
2019
-
[13]
T. A. Corcovilos and J. Mittal, Appl. Opt.58, 2256 (2019)
2019
-
[14]
Kohlert, S
T. Kohlert, S. Scherg, X. Li, H. P. Lüschen, S. Das Sarma, I. Bloch, and M. Aidelsburger, Phys. Rev. Lett.122, 170403 (2019)
2019
-
[15]
Kohmoto, L
M. Kohmoto, L. P. Kadanoff, and C. Tang, Phys. Rev. Lett.50, 1870 (1983)
1983
-
[16]
Ostlund, R
S. Ostlund, R. Pandit, D. Rand, H. J. Schellnhuber, and E. D. Siggia, Phys. Rev. Lett.50, 1873 (1983)
1983
-
[17]
Kohmoto and J
M. Kohmoto and J. R. Banavar, Phys. Rev. B34, 563 (1986)
1986
-
[18]
Matter3, 7255 (1991)
J.Q.You,J.R.Yan,T.Xie,X.Zeng, andJ.X.Zhong,J.Phys.: Condens. Matter3, 7255 (1991)
1991
-
[19]
J.H.Han,D.J.Thouless,H.Hiramoto, andM.Kohmoto,Phys. Rev. B50, 11365 (1994)
1994
-
[20]
F. Liu, S. Ghosh, and Y. D. Chong, Phys. Rev. B91, 14108 (2015)
2015
-
[21]
Aubry and G
S. Aubry and G. André, Ann. Israel Phys. Soc.3, 133 (1980)
1980
-
[22]
Devakul and D
T. Devakul and D. A. Huse, Phys. Rev. B96, 214201 (2017)
2017
-
[23]
Sutradhar, S
J. Sutradhar, S. Mukerjee, R. Pandit, and S. Banerjee, Phys. Rev. B99, 224204 (2019)
2019
-
[24]
Vojta, J
T. Vojta, J. Low Temp. Phys.161, 299 (2010)
2010
-
[25]
Knap, Phys
S.Gopalakrishnan,K.Agarwal,E.A.Demler,D.A.Huse, and M. Knap, Phys. Rev. B93, 134206 (2016)
2016
-
[26]
119, 075702 (2017)
V.Khemani,D.N.Sheng, andD.A.Huse,Phys.Rev.Lett. 119, 075702 (2017)
2017
-
[27]
W.DeRoeckandF.Huveneers,Phys.Rev.B 95,155129(2017)
2017
-
[28]
Y.E.Kraus,Y.Lahini,Z.Ringel,M.Verbin, andO.Zilberberg, Phys. Rev. Lett.109, 106402 (2012)
2012
-
[29]
Senechal,Quasicrystals and geometry (Cambridge Univer- sity Press, Cambridge, 1995)
M. Senechal,Quasicrystals and geometry (Cambridge Univer- sity Press, Cambridge, 1995)
1995
-
[30]
111, 226401 (2013)
Y.E.Kraus,Z.Ringel, andO.Zilberberg,Phys.Rev.Lett. 111, 226401 (2013)
2013
-
[31]
Zilberberg, S
O. Zilberberg, S. Huang, J. Guglielmon, M. Wang, K. P. Chen, Y. E. Kraus, and M. C. Rechtsman, Nature553, 59 (2018)
2018
-
[32]
Rossignolo and L
M. Rossignolo and L. Dell’Anna, Phys. Rev. B99, 054211 (2019)
2019
-
[33]
I. M. Suslov, Sov. Phys. JETP56, 612 (1982)
1982
-
[34]
48,1043 (1982)
C.M.SoukoulisandE.N.Economou,Phys.Rev.Lett. 48,1043 (1982)
1982
-
[35]
S. Ya. Jitomirskaya, Ann. Math.150, 1159 (1999)
1999
-
[36]
N. Mott, J. Phys. C20, 3075 (1987)
1987
-
[37]
S.Lang, Introduction to Diophantine Approximations(Springer- Verlag, 1995). 10
1995
-
[38]
Forβ = 1/ √ 2,wefindtheseboundsbyexactdiagonalisationto be≈± 2.703J; see black lines in Fig. 3
-
[39]
A.-M. Guo, X. C. Xie, and Q.-f. Sun, Phys. Rev. B89, 075434 (2014)
2014
-
[40]
Li, J.-p
H.-j. Li, J.-p. Dou, and G. Huang, Phys. Rev. A89, 033843 (2014)
2014
-
[41]
A.ChandranandC.R.Laumann,Phys.Rev.X 7,031061(2017)
2017
-
[42]
Huang and W
B. Huang and W. V. Liu, Phys. Rev. B100, 144202 (2019)
2019
-
[43]
Kariyado and A
T. Kariyado and A. Vishwanath, Phys. Rev. Res.1, 033076 (2019)
2019
-
[44]
Szabó and U
A. Szabó and U. Schneider, Phys. Rev. B98, 134201 (2018)
2018
-
[45]
Abrahams, P
E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. Lett.42, 673 (1979)
1979
-
[46]
H. P. Lüschen, S. Scherg, T. Kohlert, M. Schreiber, P. Bordia, X.Li,S.DasSarma, andI.Bloch,Phys.Rev.Lett. 120,160404 (2018)
2018
-
[47]
E. H. Lieb, R. Seiringer, and J. Yngvason, Phys. Rev. B66, 134529 (2002)
2002
-
[48]
Roth and K
R. Roth and K. Burnett, Phys. Rev. A68, 23604 (2003)
2003
-
[49]
For continuum sys- tems, p0 is set tokΘ/Lπ
For lattice systems, the twist can be implemented by adding a complex phase Θ/L to all hopping terms. For continuum sys- tems, p0 is set tokΘ/Lπ
-
[50]
M. N. Barber, inPhase Transitions and Critical Phenomena , Vol.8,editedbyC.DombandJ.L.Lebowitz(AcademicPress, London, 1983) pp. 145–266
1983
- [51]
- [52]
-
[53]
Bellissard, A
J. Bellissard, A. Bovier, and J.-M. Ghez, Rev. Math. Phys.04, 1 (1992)
1992
-
[54]
Weinberg, C
M. Weinberg, C. Ölschläger, C. Sträter, S. Prelle, A. Eckardt, K.Sengstock, andJ.Simonet,Phys.Rev.A 92,043621(2015)
2015
-
[55]
Reitter, J
M. Reitter, J. Näger, K. Wintersperger, C. Sträter, I. Bloch, A. Eckardt, and U. Schneider, Phys. Rev. Lett.119, 200402 (2017)
2017
-
[56]
Bordia, H
P. Bordia, H. Lüschen, U. Schneider, M. Knap, and I. Bloch, Nat. Phys.13, 460 (2017)
2017
-
[100]
rarelines
and atλ = 0.1 (linear scale to the right). The former is strongly localised with a very small localisation length (ξ≈ 0.1); the latter is fully extended. (c) Median energy state (E≈ 0) at λ = 40. The wave function is concentrated on a few lattice lines, along which the quasipe...
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.