REVIEW 1 major objections 6 minor 14 references
Eigenvalues of the Magnetic Neumann Laplacian on Domains with Peaks
T0 review · 1 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read For cusp-shaped boundaries, magnetic Laplacian eigenvalues grow as λ^{2/(q+1)}, with sharpness q setting the rate.
desk verdict New asymptotic law for magnetic Neumann Laplacians on outward power-law peaks; the proof is solid, the exact-geometry assumption is the main limitation, and it deserves expert peer review. 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 carrying object is the coordinate change Φ(s,t) = (s, s^q t), which maps the rectangle I×(a_-,a_+) onto the peak section V_I = {a_- x^q < y < a_+ x^q}. It turns the magnetic form into an expression on the rectangle whose leading part, after projection onto t-independent functions, is the one-dimensional model T_λ = -d²/ds² + ((a_+-a_-)^2/12) λ² $s^{{2q}}$ + q(q-2)/(4s²). The proof partitions the cusp into three scales: (0,$λ^{{-κ}}$), ($b_1λ^{{-κ}}$, $b_2λ^{{-1/(2q)}}$), and ($b_3λ^{{-1/(2q)}}$, b). The innermost section supplies eigenvalues of order $λ^{{2/(q+1)}}$; the two outer sections are lower-bounded by operators with faster growth. Min-Max comparison (Lemma 2.1) and the IMS formula transfer these bounds to the full domain, and optimizing the cutoff parameter κ yields the stated remainder.
What would settle it
Take the innermost section I=(0,$λ^{{-κ}}$) and compute E_1(Q^I_λ) numerically for increasing λ at fixed q, a_-, a_+ (for instance q=2). If E_1(Q^I_λ)/$λ^{{2/(q+1)}}$ does not converge to (a_+-a_-)^{2/(q+1)} E_1(T) with deviations no larger than O($λ^{{2κ}}$), the theorem's leading term is false. Alternatively, perturb the upper side to y = a_+ x^q + ε $x^{{q+δ}}$ and check whether the first eigenvalue changes at leading order; if it does, the exact-form hypothesis is indispensable.
Extended reading notes
Core claim
The central claim is that the low-lying spectrum of the magnetic Neumann Laplacian on a domain with an outward peak is asymptotically governed by the local power-law cusp alone. For each fixed n, E_n(N^Ω_λ) = $λ^{{2/(q+1)}}$ (a_+-a_-)^{2/(q+1)} E_n(T) + O($λ^{{2κ}}$), where κ = (2q+3)/((q+1)(2+3q)) and T is the Friedrichs extension of -d²/ds² + $s^{{2q}}$/12 + q(q-2)/(4s²) on L²(0,∞). The exponent 2/(q+1) is less than 1 for q>1, so the eigenvalues diverge slower than the linear-in-λ growth known for smooth domains and slower than the corner-driven growth known for curvilinear polygons. The result holds for all n simultaneously, meaning the cusp controls the whole bottom of the spectrum, not just the ground state.
Load-bearing premise
The load-bearing premise is that, in a fixed neighbourhood of the origin, the domain is exactly the straight power-law wedge V_b = {a_- x^q < y < a_+ x^q}; if the boundary only approximates that shape, or is perturbed by any subleading term, the theorem's constant and even the exponent are not established by the proof.
Editorial extensions
If this is right
- Every fixed eigenvalue E_n(N^Ω_λ) has the same leading power λ^{2/(q+1)}; the peak determines the entire low-lying spectrum, not only the ground state.
- A sharper peak (larger q) gives a smaller exponent 2/(q+1), so the divergence is slower and the model operator T fixes the corresponding constant for each n.
- The section (0,λ^{-κ}) is the only region contributing at leading order; 'near the peak' and 'away from the peak' regions contribute only to the remainder.
- The remainder bound O(λ^{2κ}) with κ=(2q+3)/((q+1)(2+3q)) quantifies how quickly the asymptotic law appears and is the price paid for the stronger singularity.
Reading between the lines
- If the cusp is only asymptotically a power law, the same exponent might survive but the prefactor (a_+-a_-)^{2/(q+1)}E_n(T) should be checked against a subleading term such as ε x^{q+δ}; the paper's exact-form assumption does not cover this, and a dependence on ε at leading order would be a natural test.
- The same straightening procedure should transfer to three-dimensional peaks or conical tips, producing a model operator with an angular coordinate and a scaling exponent set by the tip's power; this is a direct extension not stated in the paper.
- In the superconductivity reading, cusped samples would be predicted to nucleate at lower field strength than smooth or cornered ones, because slower eigenvalue growth lowers the energy of the first magnetic state; this is a physical consequence the paper leaves implicit.
- The model eigenvalues E_n(T) depend on q through s^{2q}/12 + q(q-2)/(4s²); tabulating them numerically for a range of q would produce a benchmark curve that direct finite-element computations on cusped domains could confirm or refute.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the magnetic Neumann Laplacian with constant magnetic field strength lambda on bounded planar domains having an outward power-law peak V_b = {a_- x^q < y < a_+ x^q}. Theorem 1.1 claims that for the n-th eigenvalue, E_n(N_lambda^Omega) = lambda^{2/(q+1)} (a_+ - a_-)^{2/(q+1)} E_n(T) + O(lambda^{2 kappa}) with kappa = (2q+3)/((q+1)(2+3q)), where T is the one-dimensional operator -h'' + h(s^{2q}/12 + q(q-2)/(4s^2)). The proof derives a one-dimensional model operator by explicit changes of variables and gauge transformations, establishes matching upper and lower bounds (Propositions 4.2 and 4.3) after splitting the peak into three regions, and optimizes the remainder exponent. The paper is clearly written and the proof strategy is coherent; however, I find a sign error in the s^{-2} coefficient of the model operator that changes the limiting constant E_n(T) and therefore the main theorem as stated.
Significance. If corrected, this would be a genuine extension of the magnetic Neumann Laplacian asymptotics from smooth boundaries and corners to zero-angle cusps. The leading-order exponent 2/(q+1), the explicit dependence on a_+ - a_-, and the optimized remainder are concrete and falsifiable predictions. A particular strength is that the model operator is derived through unitary transformations and operator comparisons rather than fitted, and the upper and lower bounds match at the chosen kappa. The error found below affects the limiting constant, not the overall method, so the result is likely repairable by a systematic sign correction.
major comments (1)
- [Definition 2.3 and Lemma 3.4] The coefficient q(q-2)/(4s^2) in the model operator is incorrect; the calculation in Lemma 3.4 yields q(q+2)/(4s^2). For v = g(s) in the form ell_I^lambda of Lemma 3.3, the cross term coming from -q/(2s) g is -(q/(2s))(|g|^2)', and after integration by parts this contributes +q/(2s^2)|g|^2. Hence the total coefficient of s^{-2}|g|^2 is q^2/4 + q/2 = q(q+2)/4. For example, when q=2 and a_- = -a_+ = 1, the transformed form is integral |g'|^2 + (lambda^2 s^4/3 + 2/s^2)|g|^2 ds, not integral |g'|^2 + (lambda^2 s^4/3)|g|^2 ds. This error propagates into Corollary 3.5, the lower bound in Lemma 3.6 and Corollary 3.7, and the upper and lower bounds in Propositions 4.2 and 4.3, because all these arguments compare with the wrong operator T. The proof structure can likely be repaired by replacing q(q-2) with q(q+2) everywhere, but as written the constant E_n(T) in Theorem 1.1 is not the limiting constant for the actual operator.
minor comments (6)
- [Definition 4.1 / Theorem 1.1] The theorem is proved only for domains that coincide exactly with V_b in a neighbourhood of the peak. The proof uses this exact equality in the IMS partition, in the construction of the tubular neighbourhoods, and in the inclusion V_I subset Omega; merely asymptotic power-law cusps or small perturbations of the sides are not covered. This is a scope limitation rather than an inconsistency, but it should be stated explicitly in the introduction.
- [Lemmas 3.11, 3.12 and Proposition 4.3] Please state the precise hypotheses of the external results invoked: [5, Appendix B] for the unitary equivalence in curved strips, [12, Lemma 8] for the density of H^1_I(V_I), and [3, Theorem 2.1] for the lower bound on the cut-off curvilinear polygon, so that the reader can verify that they apply to the present peak geometries.
- [Lemma 2.2] In the IMS formula, the term ||grad chi_j f||^2 should be written as integral |grad chi_j|^2 |f|^2 dx, or as || |grad chi_j| f ||^2, to avoid confusion with a derivative applied to the product f.
- [Section 3 and Theorem 1.1] The symbol kappa is used both for the interval-endpoint exponent in Lemma 2.9 and Section 3 and for the optimized remainder exponent in Theorem 1.1; using distinct letters would improve readability.
- [Lemma 3.8] The proof speaks of an eigenfunction of Q_I^lambda but then computes with the unitarily equivalent form ell_I^lambda; please make the notation consistent.
- [References] Reference [12] lists an author as 'F. Sk'; this is presumably a typesetting corruption of 'F. Šk' or a similar name and should be corrected.
Circularity Check
No circularity: the model operator T is derived by unitary transformation and a matching two-sided estimate, not fitted or imported by self-citation.
full rationale
The derivation of the main term is self-contained in the sense relevant to circularity. Lemma 3.3 is an exact unitary change of variables on the peak section. Lemma 3.4 shows that the one-dimensional operator T_I^λ is obtained by restricting ℓ_I^λ to t-independent test functions, which is a legitimate variational upper bound via Lemma 2.1, not a fit: the coefficients (a_+-a_-)^2/12 and q(q-2)/4 come from explicit integration over the transverse interval, not from matching eigenvalues. The converse is proved in Lemma 3.6, where the projection P0 onto t-independent functions gives E_n(L_I^λ) ≥ E_n(T_I^λ)(1-o(1)), so the model operator is shown to control the true eigenvalues from both sides; no fitted parameter is renamed as a prediction. Proposition 2.4 derives λ^{2/(q+1)}(a_+-a_-)^{2/(q+1)} E_n(T) by unitary rescaling from T_λ, and the final remainder exponent is obtained by optimizing κ in Propositions 4.2-4.3. The paper cites [3], [5], [9], and [12] for standard IMS, tubular-neighbourhood, density, and corner estimates; none are self-citations of the author and none carry the leading-order claim. The exact equality Ω∩(-b,b)^2 = V_b in Definition 4.1 and Theorem 1.1 is a hypothesis restricting the theorem to the exact power-law peak class; this is a scope limitation, not a circular step. Therefore no circularity is present.
Assumptions & free parameters
assumptions (5)
- standard math One-dimensional Hardy inequality (2.1)
- standard math External lower bound for curvilinear polygons from [3, Theorem 2.1]
- standard math Unitary tubular-neighbourhood description from [5, Appendix B]
- standard math Density result [12, Lemma 8] for H^1_I(V_I)
- domain assumption Exact geometric form of the peak: Omega intersect (-b,b)^2 = V_b
Cite this review
Pith. "Pith review of Eigenvalues of the Magnetic Neumann Laplacian on Domains with Peaks." pith.science (2026). https://pith.science/paper/CXYPMRZB
@misc{pith2026260809469,
author = {Pith},
title = {Pith review of: Eigenvalues of the Magnetic Neumann Laplacian on Domains with Peaks},
year = {2026},
howpublished = {\url{https://pith.science/paper/CXYPMRZB}},
note = {Machine review of arXiv:2608.09469}
}
abstract
We consider the magnetic Neumann Laplacian on bounded domains in $\mathbb{R}^2$ with outward peaks. The operator is associated with a large magnetic field depending on a parameter $\lambda$, and we investigate the behaviour of the eigenvalues as $\lambda$ tends to $+\infty$. We show that their asymptotic expansion is influenced by the sharpness $q$ and geometry of the peak, and that its main term is of order $\lambda^{\frac{2}{q+1}}$. This is an extension of previous works on magnetic Neumann Laplacians in smooth domains and domains with corners.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[12]
Laplacian eigenvalues for large negative Robin parameters on domains with outward peaks
K. Pankrashkin, F. Sk, and M. Vogel. “Laplacian eigenvalues for large negative Robin parameters on domains with outward peaks”. In:Nonlinear Differential Equations and Applications NoDEA 32.129 (2025)
work page 2025
-
[1]
Onset of superconductivity in decreasing fields for general domains
A. Bernoff and P. Sternberg. “Onset of superconductivity in decreasing fields for general domains”. In:Journal of Mathematical Physics39.3 (1998), pp. 1272–1284
work page 1998
-
[2]
V. Bonnailie-No¨ el. “On the fundamental state energy for a Schr¨ odinger operator with magnetic field in domains with corners”. In:Asymptotic Analysis41 (2005), pp. 215–258
work page 2005
-
[3]
V. Bonnaillie-No¨ el and M. Dauge. “Asymptotics for the low-lying eigenstates of the Schr¨ odinger operator with magnetic field near corners”. In:Annales Henri Poincar´ e7 (2006), pp. 899–931
work page 2006
-
[4]
D.E. Edmunds and W.D. Evans.Spectral theory and differential operators. Oxford: Clarendon Press, 1987
work page 1987
-
[5]
Accurate eigenvalue asymptotics for the magnetic Neumann Lapla- cian
S. Fournais and B. Helffer. “Accurate eigenvalue asymptotics for the magnetic Neumann Lapla- cian”. In:Annales de l’Institut Fourier56.1 (2006), pp. 1–67
work page 2006
-
[6]
S. Fournais and B. Helffer.Spectral Methods in Surface Superconductivity. Birkhaeuser Boston, MA, 2010
work page 2010
-
[7]
Semiclassical analysis for the ground state energy of a Schr¨ odinger operator with magnetic wells
B. Helffer and A. Mohamed. “Semiclassical analysis for the ground state energy of a Schr¨ odinger operator with magnetic wells”. In:Journal of Functional Analysis138.1 (1996), pp. 40–81
work page 1996
Show all 14 references
-
[8]
Magnetic bottles in connection with superconductivity
B. Helffer and A. Morame. “Magnetic bottles in connection with superconductivity”. In:Journal of Functional Analysis185.2 (2001), pp. 604–680
2001
-
[9]
Bound states for the magnetic Neumann Laplacian in planar sectors
A. Kachmar and M. Sundqvist. “Bound states for the magnetic Neumann Laplacian in planar sectors”. In:arXiv2607.12600 (2026), https://arxiv.org/abs/2607.12600
2026 arXiv
-
[10]
Gauge invariant eigenvalue problems inR 2 andR 2 +
K. Lu and X. Pan. “Gauge invariant eigenvalue problems inR 2 andR 2 +”. In:Transactions of the American Mathematical Society352.3 (1999), pp. 1247–1276
1999
-
[11]
Estimates of the upper critical field for the Ginzburg–Landau equations of superconductivity
K. Lu and X.B. Pan. “Estimates of the upper critical field for the Ginzburg–Landau equations of superconductivity”. In:Physica D: Nonlinear Phenomena127.1 (1999), pp. 73–104
1999
-
[13]
Raymond.Bound States of the Magnetic Schr¨ odinger Operator
N. Raymond.Bound States of the Magnetic Schr¨ odinger Operator. European Mathematical So- ciety, 2017
2017
-
[14]
Weidmann.Spectral Theory of Ordinary Differential Operators
J. Weidmann.Spectral Theory of Ordinary Differential Operators. Springer Berlin, Heidelberg, 1987
1987
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.