REVIEW 2 major objections 4 minor 37 references
Eigenfunctions with double exponential rate of localization
T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This paper constructs nonzero eigenfunctions of uniformly elliptic divergence-form operators that decay at double-exponential speed in a cylinder, the fastest rate allowed.
desk verdict Real progress on sharp decay rates, but the construction currently rests on an unproved load-bearing lemma (Lemma 7.4) that the paper incorrectly claims is identical to Lemma 7.3. 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 central object is the building block, a transformation between separated solutions carried out in a slab of bounded duration while the coefficient matrix remains in the uniform regularity class $R(80,60)$, meaning ellipticity constant at most $80$ and $C^1$ norm at most $60$. The local heart of the construction is an explicit two-dimensional formula: for $u = \cos(kx)+s\cos(k'y)$, the vector field $V = \left(\frac{\sin(2kx)\cos(k'y)}{2k}, \frac{2\sin(k'y)\sin^2(kx)}{k'}\right)$ has divergence $\cos(k'y)$, and the matrix $A_s$ defined by $A_s\nabla u = V$ satisfies $\|A_s\| \le C(1+|s|)/k^2$ and $\|\nabla A_s\| \le C(1+|s|)/k$. Because these bounds are linear in the added amplitude $s$ and degrade only by fixed powers of $k$, infinitely many glued blocks keep the coefficient matrix uniformly $C^1$. A symmetric version of this local lemma, stated without a written proof, handles the removal step in the second half of the gluing.
What would settle it
Compute the $C^1$ norms of the explicit matrices $a_s = -\cos(k'y)\cos(kx)/k^2 + 2s\sin^2(k'y)/k^2$ and $b_s = -2\sin(k'y)\sin(kx)/(kk')$ on the torus for, say, $k'=2k$; if $\|A_s\|$ or $\|\nabla A_s\|$ can exceed $C(1+s)/k^2$ or $C(1+s)/k$ with a constant that grows with $k$, the infinite gluing collapses. A finite two-block run of the building block would expose the same failure numerically.
Extended reading notes
Core claim
The core discovery is an eigenfunction that saturates the Landis barrier: in the cylinder $\mathbb{T}^2 \times \mathbb{R}$, for every $\mu > 0$ there is a real uniformly elliptic matrix $A$ with $C^1$ entries, uniformly in all variables, and a nonzero uniformly $C^2$ solution $u$ to the eigenvalue equation, with $\sup_{\mathbb{T}^2 \times \{|t|\ge T\}} |u| \le C e^{-c e^{cT}}$. The construction builds one block at a time: inside a slice of bounded length, the solution transforms $\cos(kx)e^{-kt}$ into a constant multiple of $\cos(k'y)e^{-k't}$ with $k'$ up to $2k$, while the coefficient matrix stays in the fixed regularity class $R(80,60)$. Concatenating these blocks doubles the spatial frequency while time advances linearly, and the composition produces the compounding double-exponential decay. The eigenvalue term is added by a diagonal perturbation with entries proportional to $\mu/k_n^2$, a symmetrization step makes the solution even-reflectable at $t=0$, and the $A$-harmonic and parabolic statements are parallel versions of the same block construction.
Load-bearing premise
The load-bearing premise is that the local perturbation matrices from the explicit formula satisfy uniform bounds linear in the added amplitude, with constants independent of the frequency $k$; if either the written linear bound or its asserted symmetric analogue fails, the coefficient matrix will not stay uniformly $C^1$ across the infinite sequence of blocks.
Editorial extensions
If this is right
- If the construction is correct, quantitative unique continuation for uniformly elliptic divergence-form operators with $C^1$ coefficients is too weak to support the direct transfer of the known Anderson-localization strategy for Schr\"odinger operators with bounded random potentials.
- The cylinder examples show that a positive answer to the periodic-operator question requires periodicity in all directions: super-exponential decay is possible when only one direction is homogeneous.
- The decay rate is sharp: by the barriers quoted in the paper, no nonzero solution of these elliptic or parabolic equations can decay strictly faster than double-exponentially, so the examples hit the theoretical ceiling.
- In the full cylinder the eigenvalue case is the relevant one: an $A$-harmonic solution cannot decay in both directions without being trivial, by periodicity and Liouville's theorem, while the eigenfunction can.
- The parabolic example is complex-valued; whether a real-valued solution with the same double-exponential decay exists is left open in the paper.
Reading between the lines
- One testable extension is to push the same block lemma into higher-dimensional cylinders $\mathbb{T}^d \times \mathbb{R}$; the explicit vector field would need new harmonic bookkeeping, and the uniform-in-$k$ constants would decide whether the iteration survives.
- The paper's own remark that randomness or periodicity of coefficients is not used by current methods suggests a natural next check: construct the analogous localized eigenfunction with a periodic coefficient matrix, or prove that full periodicity rules it out.
- A direct symmetry proof of the unproved symmetric variant of the local lemma would make the gluing fully self-contained; until then the second half of the construction leans on an asserted analogue.
- For the heat-equation example, a numerical simulation of a short chain of blocks could test whether the uniform bounds on the drift $B$ hold with explicit constants, giving a concrete check of the claimed double-exponential rate.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs real-valued solutions of the elliptic eigenvalue problem -div(A grad u) = lambda u in T^2 x R with uniformly elliptic, uniformly C^1 real matrices A and double-exponential decay |u| <= C exp(-c exp(c|t|)) (Theorem 1.2), together with an A-harmonic half-cylinder analogue (Theorem 1.1) and a complex-valued parabolic analogue with uniformly bounded continuous drift (Theorem 1.3). The proof is a long reduction chain: Theorem 1.2 is reduced to a half-cylinder eigenfunction construction (Theorem 4.1), which is reduced to a block lemma (Lemma 5.3), which in turn is reduced to three technical propositions (Propositions 6.1, 6.2, 6.3). The key device is to alternate between functions cos(k_n x) exp(-k_n t) and cos(k_{n+1} y) exp(-k_{n+1} t), doubling the frequency in each block, while changing the matrix coefficients so that the gluing can be done with uniformly C^1 coefficients and no fitted parameters.
Significance. If the construction is completed, the paper gives a substantial and surprising result: eigenfunctions of uniformly elliptic divergence-form operators with C^1 coefficients can decay in the non-compact direction at the maximal double-exponential rate, sharply distinguishing this class from operators with constant principal part and bounded potential. The results also give a concrete obstruction to extending quantitative unique continuation approaches in the spirit of Bourgain-Kenig to variable-coefficient divergence-form operators, and a parabolic counterpart with continuous bounded drift. The paper's strength is that the construction is fully explicit and parameter-free: the transfer constants c2, c3, and c = exp(-k/2 + 5k'/6) are forced by the interface matching conditions, and the double-exponential rate is read off from the block iteration rather than fitted. The main weakness is that one load-bearing perturbation lemma and one half of a central proposition are asserted without proof; the omitted details appear to be repairable, so the significance is conditional on supplying them.
major comments (2)
- [Section 9.2, Remark 9.1 (Lemma 7.4)] Lemma 7.4 is stated but not proved. The assertion in Remark 9.1 that its proof is 'identical' to Lemma 7.3 is not by itself a proof: the matrix Ansatz changes from [[a,b],[b,0]] to [[0,b],[b,c]], and the target divergence changes from cos(k'y) to cos(kx), so a relabeling must be spelled out. This is a load-bearing gap because Lemma 7.4 is used in the second half of Proposition 7.2, which is the only mechanism in the proof of Proposition 7.1 for removing the initial cos(kx)exp(-k sqrt(a) t) term in Step 3. The lemma is nevertheless true, and the gap can be closed explicitly: for v = s cos(kx) + cos(k'y), take V = (2 sin(kx) sin^2(k'y)/k, sin(2k'y) cos(kx)/(2k')), which satisfies div V = cos(kx), and the matrix entries are b = -2 sin(kx) sin(k'y)/(k k') and c = -cos(kx) cos(k'y)/k'^2 - 2s sin^2(kx)/(k k'), giving the stated bounds ||A_s|| <= C(1+|s|)/k^2 and ||grad A_s|| <= C(1+|s|)/k. Since the proof as published omits this, the central reduction chain is currently unverified at a critical point.
- [Section 7.2.2, end of proof of Proposition 7.2] The final paragraph of the proof of Proposition 7.2 says that the transformation from epsilon u1 + u2 to u2 is handled 'similarly' and is not included. This is not a purely cosmetic omission: the second half has a different perturbation function s(t) = epsilon alpha(t) exp((k' sqrt(b) - k sqrt(a)) t) and a different source term beta(t) = epsilon(alpha''(t) - 2k sqrt(a) alpha'(t)) u1, so the estimates for |beta|, |s'(t)|, the endpoint regularity of the matrix, and the t-derivative estimate involving partial_s A_s must all be checked with the roles of k and k' interchanged. This second half is used exactly in Step 3 of the reduction of Proposition 7.1 to Proposition 7.2, and without it the wait-time construction collapses. The details should be written out, even if they follow the first half closely.
minor comments (4)
- [Section 5.3, Step 3] The text refers to 'By Lemma 6.3' where the statement is Proposition 6.3; please correct the cross-reference.
- [Section 9.2, Lemma 7.4] There is a typo in the lemma statement: 'sataisfies' should be 'satisfies'.
- [Section 7.3, Proposition 7.5] The reduction of Proposition 6.3 to the case t1 = 0, c1 = c2 = 1 is skipped with 'we skip it'. This reduction is standard and harmless, but a one-sentence explanation of how the duration and the regularity class are invariant under the shift and rescaling would improve readability.
- [Section 1.4] The sentence 'By claiming A ≲ B we mean that A ≤ CB for some sufficiently large constant C' appears twice, with slightly different wording; please delete the duplicate.
Circularity Check
No significant circularity: the proof is a self-contained iterative construction; the unproved Lemma 7.4 is a routine omitted verification, not a circular reduction.
full rationale
The derivation chain Theorem 1.2 -> Theorem 4.1 -> Theorem 5.1 -> Lemma 5.2 -> Lemma 5.3 -> Propositions 6.1-6.3 -> Proposition 7.2 -> Lemma 7.3 is a sequence of explicit reductions in which no parameter is fitted to the target decay and no result is imported from the authors' prior work as a load-bearing premise. The transfer constants c2, c3 and C are determined by the interface equalities (122), (140), (145) and by the explicit gluing functions; the double-exponential rate then follows from kn = 2^{n0+n-1}, tn = C(n-1), and the estimates (102)-(104). Lemma 7.3 is proved in Appendix 9.2 with an explicit vector field (343) satisfying div(V) = cos(k'y) by direct computation, giving the uniform bounds (340) that keep every glued matrix in the classes R(20,10) and R(80,60). Lemma 7.4, stated without proof in Remark 9.1, is the only missing support: the x/y-swap claimed there is not literal because k' differs from k, but the analogous vector field V = (2 sin(kx) sin^2(k'y)/k, sin(2k'y) cos(kx)/(2k')) has div(V) = cos(kx) and produces the same bounds, so the omitted proof is a routine verification rather than a circular reliance on the theorem being proved. The self-citation [14] (Filonov-Krymskii) appears only in the introduction as context on Landis-type counterexamples and is not used in any reduction. No fitted-input-as-prediction, self-definitional, or uniqueness-imported-from-authors pattern occurs.
Assumptions & free parameters
free parameters (4)
- n0 (initial frequency exponent)
- Intermediate decay coefficient b = 1/9 =
1/9
- Perturbation size epsilon = 1/k^4 =
1/k^4
- Acceleration window w = 400 =
400
assumptions (4)
- standard math Uniform ellipticity of the 3D operator u_tt + div(A grad u) together with the maximum principle controls the sup-norm of glued solutions.
- standard math Landis barrier: elliptic solutions cannot decay strictly faster than double-exponentially, and the parabolic analogue holds.
- standard math Smooth cutoff functions theta with flat endpoints exist and satisfy the derivative bounds used throughout (|theta'| <= sqrt(pi), sup |theta^{(n)}| <= C_n).
- domain assumption Unique continuation holds for C1-coefficient elliptic operators, so the constructed fast-decaying solutions must not vanish on any open set.
Cite this review
Pith. "Pith review of Eigenfunctions with double exponential rate of localization." pith.science (2026). https://pith.science/paper/LTVS4MPT
@misc{pith2026250115354,
author = {Pith},
title = {Pith review of: Eigenfunctions with double exponential rate of localization},
year = {2026},
howpublished = {\url{https://pith.science/paper/LTVS4MPT}},
note = {Machine review of arXiv:2501.15354}
}
abstract
We construct a real-valued solution to the eigenvalue problem $-\text{div}(A\nabla u)=\lambda u$, $\lambda>0,$ in the cylinder $\mathbb{T}^2\times \mathbb{R}$ with a real, uniformly elliptic, and uniformly $C^1$ matrix $A$ such that $|u(x,y,t)|\leq C e^{-c e^{c|t|}}$ for some $c,C>0$. We also construct a complex-valued solution to the heat equation $u_t=\Delta u + B \nabla u$ in a half-cylinder with continuous and uniformly bounded $B$, which also decays with double exponential speed. Related classical ideas, used in the construction of counterexamples to the unique continuation by Plis and Miller, are reviewed.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
Optimal unique continuation for periodic elliptic equations on large scales
S. Armstrong, T. Kuusi, and C. Smart,Optimal unique continuation for periodic elliptic equations on large scales(2021), https: //arxiv.org/abs/2107.14248
work page Pith review arXiv 2021
- [2]
-
[3]
J. Bourgain and C. Kenig,On localization in the continuous Anderson-Bernoulli model in higher dimension, Inventiones mathe- maticae 161 (2005), no. 2, 389-426
work page 2005
-
[4]
L. Buhovsky, A. Logunov, E. Malinnikova, and M. Sodin,A discrete harmonic function bounded on a large portion ofZ2 is constant, Duke mathematical journal171 (2022), no. 6
work page 2022
-
[5]
T. H. Colding and W. P Minicozzi II,Parabolic Frequency on Manifolds, International Mathematics Research Notices2022 (2022), no. 15, 11878-11890
work page 2022
-
[6]
B. Davey,Some Quantitative Unique Continuation Results for Eigenfunctions of the Magnetic Schr¨ odinger Operator, Communi- cations in partial differential equations39 (2014), no. 5, 876-945
work page 2014
-
[7]
, On Landis’ conjecture in the plane for potentials with growth(2023), https://arxiv.org/abs/2305.05088
work page Pith review arXiv 2023
- [8]
Show all 37 references
-
[9]
Ding and C
J. Ding and C. K. Smart,Localization near the edge for the Anderson Bernoulli model on the two dimensional lattice, Inventiones mathematicae 219 (2020), no. 2, 467-506
2020
-
[10]
Duerinckx and A
M. Duerinckx and A. Gloria,Large-scale dispersive estimates for acoustic operators: homogenization meets localization(2023), https://arxiv.org/abs/2304.14046
2023 arXiv
-
[11]
D. M. Elton,Decay rates at infinity for solutions to periodic Schr¨ odinger equations, Proceedings of the Royal Society of Edinburgh. Section A. Mathematics150 (2020), no. 3, 1113-1126
2020
-
[12]
Escauriaza, C
L. Escauriaza, C. E. Kenig, G. Ponce, and L. Vega,On Uniqueness Properties of Solutions of Schr¨ odinger Equations, Communi- cations in partial differential equations31 (2006), no. 12, 1811-1823
2006
-
[13]
Filonov,Second-Order Elliptic Equation of Divergence Form Having a Compactly Supported Solution, Journal of Mathematical Sciences 106 (2001), 3078–3086
N. Filonov,Second-Order Elliptic Equation of Divergence Form Having a Compactly Supported Solution, Journal of Mathematical Sciences 106 (2001), 3078–3086
2001
-
[14]
Filonov and S
N. Filonov and S. Krymskii,On the Landis conjecture in a cylinder(2024), https://arxiv.org/abs/2311.14491
2024 arXiv
-
[15]
Garofalo and F.-H
N. Garofalo and F.-H. Lin,Monotonicity Properties of Variational Integrals, ApWeights and Unique Continuation, Indiana Uni- versity mathematics journal35 (1986), no. 2, 245-268
1986
-
[16]
3, 347-366
, Unique continuation for elliptic operators: A geometric-variational approach, Communications on pure and applied math- ematics 40 (1987), no. 3, 347-366
1987
-
[17]
A. L. Gusarov,Liouville theorems for elliptic equations in a cylinder, Russian mathematical surveys34 (1979), no. 2, 233-234
1979
-
[18]
Kenig, Some recent quantitative unique continuation theorems, S´ eminaire´Equations aux d´ eriv´ ees partielles (Polytechnique) (2006), Exp
C. Kenig, Some recent quantitative unique continuation theorems, S´ eminaire´Equations aux d´ eriv´ ees partielles (Polytechnique) (2006), Exp. No. 20, 10 p
2006
-
[19]
Kenig, L
C. Kenig, L. Silvestre, and J.-N. Wang,On Landis’ Conjecture in the Plane, Communications in partial differential equations40 (2015), no. 4, 766-789
2015
-
[20]
Kenig, J
C. Kenig, J. Zhu, and J. Zhuge,Doubling inequalities and nodal sets in periodic elliptic homogenization, Communications in partial differential equations 47 (2022), no. 3, 549-584
2022
-
[21]
V. A. Kondratiev and E. M. Landis,Qualitative theory of second order linear partial differential equations, Partial differential equations-3, Itogi Nauki i Tekhniki. Ser. Sovrem. Probl. Mat. Fund. Napr.32 (1988), 99-215
1988
-
[22]
Kuchment,Floquet Theory for Partial Differential Equations, Operator Theory: Advances and Applications, vol
P. Kuchment,Floquet Theory for Partial Differential Equations, Operator Theory: Advances and Applications, vol. 60, Birkhauser, 2012
2012
-
[23]
3, 343-414
, An overview of periodic elliptic operators, Bulletin of the American Mathematical Society53 (2016), no. 3, 343-414
2016
-
[24]
Kukavica,Quantitative uniqueness for second-order elliptic operators, Duke mathematical journal91 (1998), no
I. Kukavica,Quantitative uniqueness for second-order elliptic operators, Duke mathematical journal91 (1998), no. 2, 225-240
1998
-
[25]
E. M. Landis,Some problems of the qualitative theory of second order elliptic equations (case of several independent variables), Russian mathematical surveys18 (1963), no. 1, 1-62
1963
-
[26]
Le Balc’h,Exponential bounds for gradient of solutions to linear elliptic and parabolic equations(2020), https://arxiv.org/ abs/2006.04582
K. Le Balc’h,Exponential bounds for gradient of solutions to linear elliptic and parabolic equations(2020), https://arxiv.org/ abs/2006.04582
2020 arXiv
-
[27]
Le Balc’h and D
K. Le Balc’h and D. Souza,Quantitative unique continuation for real-valued solutions to second order elliptic equations in the plane (2024), https://arxiv.org/abs/2401.00441
2024 arXiv
-
[28]
Li and L
L. Li and L. Zhang,Anderson–Bernoulli localization on the three-dimensional lattice and discrete unique continuation principle, Duke mathematical journal171 (2022), no. 2
2022
-
[29]
Lin and J.-N
C.-L. Lin and J.-N. Wang,Quantitative uniqueness estimates for the general second order elliptic equations, Journal of functional analysis 266 (2014), no. 8, 5108-5125
2014
-
[30]
Logunov, E
A. Logunov, E. Malinnikova, N. Nadirashvili, and F. Nazarov,The Landis conjecture on exponential decay(2020), https://arxiv. org/abs/2007.07034
2020 arXiv
-
[31]
Mandache,On a Counterexample Concerning Unique Continuation for Elliptic Equations in Divergence Form, Mathematical Physics, Analysis and Geometry1 (1998), 273–292
N. Mandache,On a Counterexample Concerning Unique Continuation for Elliptic Equations in Divergence Form, Mathematical Physics, Analysis and Geometry1 (1998), 273–292. 74
1998
-
[32]
V. Z. Meshkov,On the possible rate of decay at infinity of solutions of second order partial differential equations, Mathematics of the USSR. Sbornik72 (1992), no. 2, 343-361
1992
-
[33]
K. Miller, Nonunique continuation for uniformly parabolic and elliptic equations in self-adjoint divergence form with H¨ older continuous coefficients, Archive for rational mechanics and analysis54 (1974), no. 2, 105-117
1974
-
[34]
Plis,On non-uniqueness in Cauchy problem for an elliptic second order differential equation, Bull
A. Plis,On non-uniqueness in Cauchy problem for an elliptic second order differential equation, Bull. Acad. Polon. Sci. S´ er. Sci. Math. Astronom. Phys11 (1963), 95–100
1963
-
[35]
Rossi,The Landis conjecture with sharp rate of decay, Indiana University mathematics journal70 (2021), no
L. Rossi,The Landis conjecture with sharp rate of decay, Indiana University mathematics journal70 (2021), no. 1, 301-324
2021
-
[36]
Smart, Unique continuation for lattice Schr¨ odinger operators(2022), https://www.youtube.com/watch?v=mdSbE3fPBw0
C. Smart, Unique continuation for lattice Schr¨ odinger operators(2022), https://www.youtube.com/watch?v=mdSbE3fPBw0. ICM plenary talk
2022
-
[37]
Zhu,Quantitative uniqueness of elliptic equations, American journal of mathematics138 (2016), no
J. Zhu,Quantitative uniqueness of elliptic equations, American journal of mathematics138 (2016), no. 3, 733-762. 75
2016
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.