REVIEW 3 major objections 5 minor 34 references
Sharp Polynomial Velocity Decay Bounds for Multidimensional Periodic Schr\"odinger Operators
T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read This paper proves that in any dimension, both the asymptotic and Lieb-Robinson velocities of a periodic Schrödinger operator with large coupling decay at the sharp rate μ^{1-p0}, where p0 is the smallest period of the potential.
desk verdict Sharp velocity decay for periodic Schrödinger operators in any dimension, with a solid upper-bound proof; the sharpness claim rests on a deferred lower bound and a key lemma with a one-sentence proof that needs expansion. 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 proof uses a Rayleigh-Schrödinger perturbation expansion of the Floquet eigenvalues of εΔ+V, reorganized as sums over loops and paths on a graph. Each η^{(r)}_n is a sum over loops of length r from n to itself with rational weights in potential differences; these expansions feed a lemma asserting that the first dependence of an eigenvalue on the j-th wave-number component appears at order ε^{p_j}, with coefficient c_{j,n}(z_j+z_j^{-1}). A contour deformation in the Floquet integral then converts this polynomial vanishing into the exponential Lieb-Robinson tail, and the same expansion yields the asymptotic velocity via the group-velocity formula.
What would settle it
For a concrete non-degenerate potential with a small period, say p=(2,3) on Z^2, compute numerically the Floquet eigenvalue expansion η_n(ε,z) about ε=0 and check that the coefficient of ε^{p_j} z_j is nonzero and that all coefficients with r<p_j are independent of z_j. A single counterexample, such as the first z_j-dependence appearing at a different order, would change the decay exponent in Theorems 2.1 and 2.2.
Extended reading notes
Core claim
For a non-degenerate p-periodic potential on Z^d, the paper establishes that v_asy(H_μ) = C μ^{-p0+1} + O(μ^{-p0}) and, for every ρ0>0, a Lieb-Robinson bound with v_LR = C1 μ^{-p0+1}, both with positive constants depending only on the potential, the period, and d. The paper further shows this rate is optimal: any Lieb-Robinson bound valid for all large μ must have v_LR ≥ c μ^{-p0+1}. The same mechanism gives coordinate-wise rates: transport in the j-th direction decays like μ^{-p_j+1}, so directions with larger periods are suppressed more strongly, and if p_j=1 the corresponding velocity need not decay at all.
Load-bearing premise
The whole rate law depends on a lemma asserting that the first time an eigenvalue of the Floquet matrix feels the j-th wave-number component is at order ε^{p_j}, with a nonzero coefficient; the paper justifies this in one sentence via the loop expansion and the Feynman-Hellmann theorem, without giving the full combinatorial argument, so if that order were different for some period the exponent in Theorems 2.1 and 2.2 would change.
Editorial extensions
If this is right
- In any dimension, increasing the potential amplitude μ suppresses ballistic wave-packet spreading polynomially, with the exponent determined entirely by the shortest period of the lattice potential.
- Transport is direction-dependent: directions with larger periods see faster decay, and directions in which the potential is constant (period 1) can retain a nonzero velocity as μ grows.
- The Lieb-Robinson light-cone speed and the single-state asymptotic velocity share the same sharp scaling, making the two standard notions of quantum transport speed consistent in this regime.
- Because the argument is graph-theoretic, the same rates hold for periodic Schrödinger operators on Z^d-periodic graphs, with the period replaced by the minimal combinatorial distance between distinct vertices in the same orbit.
- The exponential tail bound holds uniformly for μ ≥ μ0, where μ0 depends on dimension, the potential separation, and the chosen decay exponent ρ0, but not on the period.
Reading between the lines
- If the same loop-length argument controls higher-order coefficients, one would predict the same μ^{-(p0-1)} scaling for related observables such as diffusion constants or current fluctuations, not just the two velocities studied here.
- The matching scaling of asymptotic and Lieb-Robinson velocities, which the paper notes was also seen in quantum walks, suggests a general principle: both speeds are governed by the same first nonconstant Floquet derivative, so any model where that derivative vanishes to a different order would show distinct velocity scalings.
- The perturbative path-sum representation is concrete enough to test numerically: for a small period such as p=(2,3), one can compute η^{(2)}_n and η^{(3)}_n explicitly and verify that the claimed z_j-dependence begins exactly at order ε^{p_j}.
- The non-degeneracy assumption (positive separation of potential values) appears to be the key boundary: degenerate potentials, where Floquet eigenvalues cross, may exhibit different or even exponential suppression, and exploring that regime would delimit the sharpness of the polynomial law.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the large-coupling dynamics of periodic Schrödinger operators H_µ = ∆ + µV on ℓ²(Z^d), with V p-periodic and non-degenerate. The main results, Theorems 2.1 and 2.2, assert that both the asymptotic velocity and the Lieb–Robinson velocity decay as µ^{-p0+1}, where p0 = min{p_1,...,p_d}, and that this rate is sharp. The upper-bound strategy is based on a Rayleigh–Schrödinger perturbation expansion for the Floquet eigenvalues (Theorem 3.1), a graph-theoretic loop expansion (Corollary 3.3), and a contour-deformation argument for the propagator blocks. The key spectral mechanism is Lemma 4.3, which asserts that the first dependence of a Floquet eigenvalue on the j-th quasimomentum appears at order ε^{p_j} with a nonzero coefficient c_{j,n}(z_j+z_j^{-1}). The lower bound for vasy follows from an explicit formula for the group velocity, while the lower bound for vLR is deferred to a generalization of a theorem from the companion paper [2].
Significance. If the stated exponents are correct, the paper proves a sharp, quantitative version of strong-coupling transport suppression in arbitrary dimensions, with the decay exponent determined by the minimal period. This is a substantial advance over one-dimensional arguments and the method is conceptually appealing: the perturbation recursion, the loop interpretation, and the contour deformation are explicit and largely self-contained. The per-direction formula (4.25) is a nice extra dividend. However, the sharp exponent rests on two components that are not fully proved in the manuscript: Lemma 4.3 (whose statement is ambiguous and whose proof is one sentence) and the higher-dimensional generalization of [2, Theorem A.1] used for the optimality of the Lieb–Robinson velocity. These gaps are load-bearing but appear fixable within the manuscript's scope, so the appropriate decision is major revision.
major comments (3)
- [Section 4.2, Lemma 4.3] Lemma 4.3(b), Eq. (4.17), is the linchpin of the exponent p0-1, but as printed it is ambiguous: it says 'For r=pj' while the right-hand side sums over j with pj=p0. If the intended statement is r=p0, it supplies the ε^{p0} term only for minimal directions; if it is meant for each j, the displayed formula does not describe the non-minimal case. The proof is a single sentence invoking Corollary 3.3 and the Feynman–Hellmann theorem. A load-bearing lemma of this kind needs a direct combinatorial proof: from Corollary 3.3, the coefficient of z_j+z_j^{-1} in η_n^{(p_j)} is a sum over loops of length p_j winding once around direction j, and one must show that no shorter loop depends on z_j and that the coefficient is nonzero. Without this, the positivity of c_{j,n} in (4.24) is not established.
- [Section 4.3, proof of Theorem 2.1] Equation (4.24) asserts ∂η_n/∂θ_i = 2 c_{i,n} sinθ_i ε^{p_i}+O(ε^{p_i+1}) for every direction i, and (4.25) asserts ∥G_i∥ = \tilde c_i ε^{p_i}+O(ε^{p_i+1}). Even if Lemma 4.3 is corrected to give the intended statement at order p0, it does not justify the per-direction asymptotics for directions with p_i>p0; those require a separate version of Lemma 4.3 at order p_i. The maximum over i in vasy only needs the minimal directions, but the paper explicitly advertises the per-direction result and uses it in the proof. Please state and prove the general per-direction lemma.
- [Section 4.3, proof of Theorem 2.2] The optimality (lower bound) of the Lieb–Robinson velocity is not proved in this manuscript; it is deferred to [2, Theorem A.1] with the remark that the proof extends to higher dimensions 'with cosmetic changes.' Since sharpness is a central claim (title, abstract, Theorem 2.2), this is load-bearing. The higher-dimensional generalization should either be proved here or stated as a theorem with full hypotheses and a proof sketch sufficient to verify the dimensional dependence of the constants. Reference [2] is a 1D result by overlapping authors and is cited as '202X', so it cannot be checked independently.
minor comments (5)
- [Section 1] Typo: 'dipersive spreading' should be 'dispersive spreading'.
- [Section 4.2, Lemma 4.3] The notation 'r=pj' conflicts with the summation index j in (4.17). Please separate the minimal-period case (r=p0) from the higher-period case (r=p_i with p_i>p0), and define c_{j,n} explicitly.
- [Theorem 2.2] In the optimality statement, clarify that vLR denotes any constant for which (2.10) holds, and that the lower bound applies to the infimum of such constants.
- [References] Reference [2] has incomplete publication data ('202X'); please provide an arXiv identifier or DOI, especially since the proof of Theorem 2.2 depends on it.
- [Section 4.3, Eq. (4.21)] The norm in ∥(G1ψ,...,Gdψ)∥ should be identified as the Euclidean norm on (ℓ²)^d to match the definition (2.5).
Circularity Check
No circular reduction: the velocity exponents are derived from the Rayleigh–Schrödinger perturbation expansion, and the only overlapping-author citation is an external lower-bound criterion, not an input of the target result.
full rationale
The derivation is self-contained for the main upper-bound argument. The exponents in Theorems 2.1 and 2.2 come from Lemma 4.3, which asserts that the first z_j-dependence of a Floquet eigenvalue appears at order epsilon^{p_j}; this is intended to follow from the paper's own Corollary 3.3 (the loop/path expansion from Theorem 3.1), not from any fitted parameter or from the velocities being bounded. No parameter is fitted to velocity data; the constants are determined by V and p, and the exponent p0 is the periodicity input that emerges from perturbation theory rather than being imposed by the conclusion. The lower-bound half of Theorem 2.2 uses (2.9) together with a 'straightforward generalization' of [2, Theorem A.1]; [2] has overlapping authors but is a prior external theorem (a general criterion relating v_LR and v_asy), so it is independent support rather than a circular re-use of the present claim. The one-sentence proof of Lemma 4.3(b) — 'The statements in the case p0 >= 2 follow from Corollary 3.3 while the statements for p0 = 1 follow from computations with the help of the Feynman–Hellmann theorem' — is an omitted combinatorial argument and therefore a proof-completeness/correctness risk, but it is not a circular reduction: no displayed equation is equivalent by construction to another, and no predicted quantity is a renamed fit. The central claim therefore does not reduce to its own assumptions.
Assumptions & free parameters
assumptions (4)
- domain assumption V is non-degenerate (injective on the fundamental cell W)
- ad hoc to paper The first dependence of a Floquet eigenvalue on the j-th quasimomentum appears at order p_j in the perturbation expansion (Lemma 4.3)
- standard math Floquet-Bloch decomposition and ballistic limit formula (2.3) / (4.22) as in [19]
- ad hoc to paper Theorem A.1 of [2] generalizes to higher dimensions
Cite this review
Pith. "Pith review of Sharp Polynomial Velocity Decay Bounds for Multidimensional Periodic Schr\"odinger Operators." pith.science (2026). https://pith.science/paper/YOCJZYOS
@misc{pith2026250904381,
author = {Pith},
title = {Pith review of: Sharp Polynomial Velocity Decay Bounds for Multidimensional Periodic Schr\"odinger Operators},
year = {2026},
howpublished = {\url{https://pith.science/paper/YOCJZYOS}},
note = {Machine review of arXiv:2509.04381}
}
read the original abstract
We investigate periodic Schr\"odinger operators in arbitrary dimensions in the large coupling regime. Our results establish that both the Lieb--Robinson velocity and the asymptotic velocity decay at an inverse polynomial rate in the coupling, with the precise exponent determined by the period of the underlying potential. In particular, we obtain sharp polynomial decay rates that capture the precise dependence on the periodic structure.
Reference graph
Works this paper leans on
-
[2]
H. Abdul-Rahman, M. Darras, C. Fischbacher, and G. Stolz. Slow propagation velocities in Schrödinger operators with large periodic potential.Ann. Henri Poincaré, 202X
-
[1]
H. Abdul-Rahman, M. Cedzich, G. Stolz, and A. H. Werner. Exponential suppression of transport in periodic electric quantum walks and skew-shift cmv matrices
-
[3]
H. Abdul-Rahman and G. Stolz. Exponentially decaying velocity bounds of quantum walks in peri- odic fields.Communications in Mathematical Physics, 403:1297–1327, 2023
work page 2023
-
[4]
A. Ahlbrecht, H. Vogts, A. H. Werner, and R. F. Werner. Asymptotic evolution of quantum walks with random coin.J. Math. Phys., 52(4):042201, 36, 2011. SHARP VELOCITIES FOR PERIODIC OPERATORS 15
work page 2011
-
[5]
M. Aizenman and S. Warzel. Absolutely continuous spectrum implies ballistic transport for quantum particles in a random potential on tree graphs.Journal of Mathematical Physics, 53(9):095205, 06 2012
work page 2012
-
[6]
J. Arbunich, J. Faupin, F. Pusateri, and I. M. Sigal. Maximal speed of quantum propagation for the hartree equation.Communications in Partial Differential Equations, 48(15):1–34, 2023
work page 2023
-
[7]
J. Asch and A. Knauf. Motion in periodic potentials.Nonlinearity, 11(1):175–200, 1998
work page 1998
-
[8]
N. W. Ashcroft and N. D. Mermin.Solid State Physics. Brooks Cole, 1976
work page 1976
Show all 34 references
-
[9]
F. Bloch. Über die Quantenmechanik der Elektronen in Kristallgittern. Zeitschrift für Physik, 52:555–600, 1929
1929
-
[10]
Boutet de Monvel and M
A. Boutet de Monvel and M. Sabri. Ballistic transport in periodic and random media. InFrom complex analysis to operator theory—a panorama, volume 291 ofOper. Theory Adv. Appl., pages 163–216. Birkhäuser/Springer, Cham, 2023
2023
-
[11]
Breteaux, J
S. Breteaux, J. Faupin, M. Lemm, D. H. Ou Yang, I. M. Sigal, and J. Zhang. Light cones for open quantum systems. arXiv:2303.08921
-
[12]
Cedzich, A
C. Cedzich, A. Joye, A. H. Werner, and R. F. Werner. Exponential tail estimates for quantum lattice dynamics, 2024. arxiv:2408.02108
2024
-
[13]
Damanik, J
D. Damanik, J. Fillman, and D. C. Ong. Spreading estimates for quantum walks on the integer lattice via power-law bounds on transfer matrices.J. Math. Pures Appl. (9), 105(3):293–341, 2016
2016
-
[14]
Damanik, J
D. Damanik, J. Fillman, and G. Young. Optimal dispersion for discrete periodic schrödinger opera- tors. arxiv:2505.14475
-
[15]
Damanik, M
D. Damanik, M. Lukic, and W. Yessen. Quantum dynamics of periodic and limit-periodic Jacobi and block Jacobi matrices with applications to some quantum many body problems.Commun. Math. Phys., 337(3):1535–1561, 2015
2015
-
[16]
Damanik, T
D. Damanik, T. Malinovitch, and G. Young. What is ballistic transport? in press
-
[17]
M. S. P. Eastham.The Spectral Theory of Periodic Differential Equations. Scottish Academic Press, 1973
1973
-
[18]
Faust and I
M. Faust and I. Kachkovskiy. Absence of flat bands for discrete periodic graph operators with generic potentials, 2025
2025
-
[19]
J. Fillman. Ballistic transport for periodic Jacobi operators onZd. In From operator theory to or- thogonal polynomials, combinatorics, and number theory—a volume in honor of Lance Littlejohn’s 70th birthday, volume 285 ofOper. Theory Adv. Appl., pages 57–68. Birkhäuser/Spring...
2021
-
[20]
Fillman, W
J. Fillman, W. Liu, and R. Matos. Irreducibility of the Bloch variety for finite-range Schrödinger operators. J. Funct. Anal., 283(10):Paper No. 109670, 22, 2022
2022
-
[21]
Fillman, W
J. Fillman, W. Liu, and R. Matos. Algebraic properties of the Fermi variety for periodic graph operators. Journal of Functional Analysis, 286(4):110286, 2024
2024
-
[22]
Fisher, W
L. Fisher, W. Li, and S. P. Shipman. Reducible Fermi surface for multi-layer quantum graphs including stacked graphene.Comm. Math. Phys., 385(3):1499–1534, 2021
2021
-
[23]
G. Floquet. Sur les équations différentielles linéaires à coefficients périodiques.Annales scientifiques de l’École Normale Supérieure, 12:47–88, 1883
-
[24]
J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade.Photonic crystals: molding the flow of light. Princeton University Press, Princeton, NJ, 2nd edition, 2008
2008
-
[25]
Kato.Perturbation Theory for Linear Operators
T. Kato.Perturbation Theory for Linear Operators. Classics in Mathematics. Springer, Berlin, Hei- delberg, 1995
1995
-
[26]
Kuchment
P. Kuchment. The mathematics of photonic crystals. InMathematical modeling in optical science, pages 207–272. SIAM, Philadelphia, PA, 2001
2001
-
[27]
Kuchment
P. Kuchment. An overview of periodic elliptic operators.Bull. Amer. Math. Soc. (N.S.), 53(3):343– 414, 2016
2016
-
[28]
Li and S
W. Li and S. P. Shipman. Irreducibility of the Fermi surface for planar periodic graph operators. Lett. Math. Phys., 110(9):2543–2572, 2020
2020
-
[29]
E. H. Lieb and D. W. Robinson. The finite group velocity of quantum spin systems.Communications in Mathematical Physics, 28(3):251–257, 1972. 16 H. ABDUL-RAHMAN, J. FILLMAN, C. FISCHBACHER, AND W. LIU
1972
-
[30]
W. Liu. Irreducibility of the Fermi variety for discrete periodic Schrödinger operators and embedded eigenvalues. Geom. Funct. Anal., 32(1):1–30, 2022
2022
-
[31]
Reed and B
M. Reed and B. Simon.Methods of Modern Mathematical Physics IV: Analysis of Operators. Aca- demic Press, 1978
1978
-
[32]
Sagiv, R
A. Sagiv, R. Kassem, and M. I. Weinstein. Dispersive decay estimates for periodic jacobi operators on the half-line. arxiv:2505.14498
-
[33]
S. P. Shipman. Reducible Fermi surfaces for non-symmetric bilayer quantum-graph operators.J. Spectr. Theory, 10(1):33–72, 2020
2020
-
[34]
M. C. Tran, A. Y. Guo, C. L. Baldwin, A. Ehrenberg, A. V. Gorshkov, and A. Lucas. Lieb-Robinson light cone for power-law interactions.Physical Review Letters, 127(16):160401, 2021. [H. Abdul-Rahman] Department of Mathematical Sciences, United Arab Emirates Uni- versity, Al Ain...
2021
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.