Pith. sign in

REVIEW 5 minor 2 cited by

Eigenvalues of the Neumann magnetic Laplacian in the unit disk

T0 review · 0 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read The lowest eigenvalue of the magnetic Laplacian in a disk is governed by a strictly ordered sequence of crossings with an exact formula and explicit asymptotics.

desk verdict Solid, careful proof of the Saint-James formula for the Neumann magnetic Laplacian on the disk, with new interlacing and asymptotic results; deserves serious refereeing. read the letter →

arxiv 2411.11721 v4 pith:X5425KJW submitted 2024-11-18 math.SP

classification math.SP MSC 35P1581Q1034L2047A75
keywords magneticLaplacianNeumannboundaryconditionsunitdiskeigenvaluecrossingsSaint-JamesformulaDeGennesconstantstrongdiamagnetismsurfacesuperconductivity
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper studies the lowest eigenvalue $\lambda(\beta)$ of the Neumann magnetic Laplacian in the unit disk as the magnetic field strength $\beta$ varies. It proves that the angular-momentum branches $\lambda(n,\beta)$ meet pairwise at unique crossing points $\beta_n$, which form a strictly increasing sequence, and that between consecutive crossings the lowest eigenvalue is exactly $\lambda(n,\beta)$. At each crossing, an exact identity due to the physicist Saint-James connects $\beta_n$, $n$, and the reduced eigenvalue $\eta = \lambda/\beta$. Combining this identity with high-order expansions of a boundary model operator yields the complete large-field asymptotics $\beta_n = 2n + \xi_1 n^{1/2} + \kappa_0 + O(n^{-1/2})$, so the low-lying spectrum is fully organized by these ordered crossings.

What carries the argument

The load-bearing object is the sequence of crossing points $\beta_n$ satisfying $\lambda(n,\beta_n)=\lambda(n+1,\beta_n)$, together with the Saint-James identity, an exact algebraic relation that holds at each such point. The proof of the identity uses the Kummer-function representation of the eigenfunctions and their recursion relations, expressing the intersection condition as the vanishing of a $2\times2$ determinant. For the asymptotics, the central mechanism is the uniform high-order expansion of the boundary model operator $h(\delta,\beta)$ as a sum $\Theta_0 + \lambda_1 \beta^{-1/2} + \lambda_2(\delta)\beta^{-1} + \cdots$ with polynomial coefficients $\lambda_j(\delta)$; inserting this expansion into the exact crossing identity produces the full series for $\beta_n$ and the derivative formulas.

What would settle it

Compute the crossing points $\beta_n$ numerically to high precision for $n$ up to, say, $10^4$ by solving the implicit equation $\Phi(\nu,n)=0$ given in Section 6, and compare $\beta_n - 2n - \xi_1 n^{1/2} - \kappa_0$ with the predicted $O(n^{-1/2})$ decay; if the leading coefficient of $n^{1/2}$ differs from $-2^{3/2}\xi_0$, or if any $\eta^*_n$ exceeds the De Gennes constant $\Theta_0 \approx 0.590106$, the paper's main asymptotic and conjectural claims would be falsified.

Watch

Extended reading notes

Core claim

The central claim is that the infimum $\lambda(\beta) = \inf_{n\in\mathbb{N}}\lambda(n,\beta)$ is realized by exactly one branch $n$ on each interval between successive crossings, with unique, strictly increasing crossing points $\beta_n$ (Theorem 1.2). At every crossing the exact Saint-James identity $\beta = 2\eta + 2n + 1 + \sqrt{(2\eta+1)^2 + 8n\eta}$ holds (Theorem 1.1). In the large-field limit the crossings have the explicit expansion $\beta_n = 2n - 2^{3/2}\xi_0\, n^{1/2} + (1 - 2\delta_0 + 2\xi_0^2) + O(n^{-1/2})$, where $\xi_0$ is the minimizer of the De Gennes model and $\delta_0 = \frac12\Theta_0^{-1/2} C_1$ is a computable constant (Theorem 4.9). Consequently the paper obtains the asymptotics of $\lambda(\beta_n)$ and of the one-sided derivatives of $\lambda$ at $\beta_n$, which converge to distinct positive limits $\Theta_0 \mp \frac{3}{2}C_1|\xi_0|$.

Load-bearing premise

The arbitrary-order asymptotic expansion of the model operator $h(\delta,\beta)$, stated as Proposition 4.3, is assumed to hold uniformly for $\delta$ in bounded sets, so that the crossing-point expansions of Theorems 4.7 and 4.9 rest on it.

Editorial extensions

If this is right

  • On each interval $[\beta_{n-1},\beta_n]$, the lowest eigenvalue is exactly $\lambda(n,\beta)$, so the crossings give a complete staircase description of the low-lying spectrum.
  • The asymptotic expansion $\beta_n = 2n + \xi_1 n^{1/2} + \kappa_0 + O(n^{-1/2})$ implies the crossings are asymptotically spaced by 2, with $\beta_{n+1}-\beta_n$ eventually decreasing.
  • The one-sided derivatives at $\beta_n$ converge to $\Theta_0 \pm \frac32 C_1|\xi_0|$, so the right derivative stays positive in the large-$n$ limit, supporting strong diamagnetism.
  • If the sequence $\eta^*_n = \eta(n,\beta_n)$ is increasing (Conjecture 1.4), then $\lambda(\beta) < \Theta_0\beta$ for all $\beta>0$ (Conjecture 1.3).

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same crossing analysis may extend to other radially symmetric planar domains or to balls in higher dimensions, where analogous exact identities could be derived from the special-function representation of eigenfunctions.
  • The newly computed constant $\delta_0$ makes it possible in principle to generate the higher coefficients $\hat\kappa_j$ in the $\beta_n$ expansion explicitly, a step the paper does not carry out.
  • A direct high-precision computation of $\eta^*_n$ for much larger $n$ (e.g., up to $10^5$) would test Conjecture 1.4 more stringently than the tabulated range $n\le 400$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

0 major / 5 minor

Summary. The paper studies the lowest eigenvalue λ(β) of the Neumann magnetic Laplacian in the unit disk. Using the Fourier decomposition and Sturm–Liouville theory, the authors prove a formula due to Saint-James expressing the intersection point β of the eigenvalue curves for angular momenta n and n+1 in terms of the common eigenvalue. They combine this formula with the asymptotic analysis of Fournais–Helffer to prove uniqueness and monotonicity of the crossing points, to describe λ(β) through an ordered sequence of crossings, and to derive complete asymptotic expansions for the crossing points β_n and for the left and right derivatives of λ at those points. A new identity for the constant δ0 is obtained, and extensive numerical computations are presented as evidence for the conjectures concerning the De Gennes constant and strong diamagnetism.

Significance. The paper gives a complete qualitative description of the low-lying spectrum of the Neumann magnetic Laplacian in the disk, with explicit constants in the leading asymptotic terms. The derivation of the Saint-James formula via a two-by-two determinant is clean, and its combination with the independent asymptotics of Fournais–Helffer yields new results such as δ0 = (1/2)Θ0^{-1/2}C1 and the expansion β_n = 2n + ξ1 n^{1/2} + κ0 + O(n^{-1/2}). The numerical results are consistent with the proven asymptotics and provide falsifiable predictions. The proofs are detailed, and the paper is a valuable contribution to the spectral theory of magnetic Laplacians.

minor comments (5)
  1. [Table 3] Table 3 appears to list values of β_n rather than η*_n: the entries coincide with the numbers β_n^{(1)} in Table 1, while the definition of ε_n refers to η*_n; the column header or the definition of ε_n should be corrected.
  2. [Proof of Theorem 4.7, after Eq. (4.30)] The sentence following Eq. (4.30) cites Eq. (4.27) as an input for Eq. (4.31); Eq. (4.31) follows directly from (4.1) and (4.30), so the reference to (4.27) appears to be a typo.
  3. [Eq. (4.24)] The remainder term in Eq. (4.24) is written as O(β^{-(N+1)/2}) after multiplication by β; to be consistent with (4.20), the remainder inside the parentheses should be O(β^{-(N+3)/2}) or the index N should be shifted.
  4. [Proposition 4.3] Proposition 4.3 states an arbitrary-order expansion with only a proof sketch; since the main asymptotic results rely on this expansion, a precise citation of the exact statement in Fournais–Persson [17] would be helpful to the reader.
  5. [Section 4.3, proof of Theorem 4.9] The induction step for the higher-order coefficients κ_j is stated as 'by recursion' without details; making the recursion explicit would improve readability, although the procedure is standard.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the Saint-James formula is proven from the exact eigenvalue equations in Section 3.4, and Theorem 4.9 follows from consistency between that formula and published Fournais-Helffer expansions; no fitted parameter is renamed as a prediction (minor arithmetic slip in Remark 4.8 noted, not circular).

full rationale

The derivation chain is self-contained at every link. Theorem 1.1 (the Saint-James formula) is proven in Section 3.4 by rewriting the two exact eigenvalue equations (3.5) and (3.6), which express the Neumann condition (g_n,β)'(β^{1/2}) = 0 for the two angular modes, as a singular 2x2 linear system in X = M(ν,n+1,x) and Y = M(ν,n+2,x), and then setting the determinant to zero; equation (3.14) is thus a consequence of the paper's own eigenvalue equations, not an input. Theorem 1.2 then follows from Theorem 1.1 together with the published Dauge-Helffer monotonicity results [8,9]; the sign computations in Proposition 3.9 use only (3.12) and the bound β* > 2(n+1). Section 4 rests on the model-operator expansion (4.20), taken from the published, peer-reviewed works [13,17]; those results are externally falsifiable, state assumptions that do not include the crossing theorems proved here, and are further supported in Proposition 4.3 by an in-paper quasimode construction (ψ_N built from rapidly decaying u_j) giving the uniform O(β^{-(N+1)/2}) remainder for δ in bounded sets - this is the flagged weakest assumption, and it is adequately supported. The new content of Theorem 4.7, δ(n,β_n) = δ0 - 1/2 + O(β_n^{-1/2}), is derived by comparing the two equal expansions (4.24) and (4.26) at the crossing and using Δδ = 1 from (4.1); no fitting is involved. Theorem 4.9 and the identity δ0 = (1/2)Θ0^{-1/2}C1 (Remark 4.8) are obtained by inserting the expansions into the exact Saint-James formula (3.15) and matching coefficients: the β^{1/2} terms cancel since Θ0^{1/2} = -ξ0, and matching the constant terms yields δ0 = C1/(2Θ0^{1/2}) as a genuine consistency condition. No fitted parameter is renamed as a prediction; C1, ξ0, Θ0, δ0 enter as inputs from prior work. Two flags that are not circularity: first, Remark 4.8 quotes '~0.0975', which is inconsistent with its own formula and C1 ~ 0.254 (which give ≈ 0.1653), and Table 1 data (δ(400,β_400) ≈ -0.34 versus δ0 - 1/2) supports δ0 ≈ 0.165; this is an arithmetic slip. Second, the Section 6 agreement between two numerical methods is an internal consistency check because both derive from the same exact equations, so it does not independently confirm the (already proven) formula. Self-citation is heavy (Helffer in [8,9,12,13,14,17,20,21]) but load-bearing citations are real, published evidence and do not presuppose the target conclusions.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters are fitted by the paper. The central results rest on standard Sturm-Liouville theory, known asymptotic expansions from Fournais-Helffer [13] and Fournais-Persson [17], properties of Kummer functions, and numerical values of universal constants from prior work. No new entities are postulated.

assumptions (4)
  • standard math Sturm-Liouville eigenvalue variation theorems of Dauge-Helffer ([8,9]) for boundary value problems with variable interval length.
    Used in Proposition 3.3 to locate minima and monotonicity of η(n,β).
  • domain assumption Uniform high-order asymptotic expansion of the model operator eigenvalues e_{δ,β} (Proposition 4.3) from Fournais-Helffer [13] and Fournais-Persson [17].
    Needed for Corollary 4.5 and all subsequent asymptotic derivations; not fully reproved in this paper.
  • standard math Kummer function identities and representation formulas from DLMF Chapter 13, including the derivative formula ∂_z M(a,b,z)=a/b M(a+1,b+1,z).
    Used in derivation of the Saint-James formula in Section 3.4.
  • domain assumption Numerical values of universal constants Θ0≈0.590106, C1≈0.254, ξ0≈-0.768, taken from prior literature ([5],[13],[27]).
    Used for numerical illustrations and evaluation of limits; not derived in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Eigenvalues of the Neumann magnetic Laplacian in the unit disk." pith.science (2026). https://pith.science/paper/X5425KJW

@misc{pith2026241111721,
  author       = {Pith},
  title        = {Pith review of: Eigenvalues of the Neumann magnetic Laplacian in the unit disk},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/X5425KJW}},
  note         = {Machine review of arXiv:2411.11721}
}
abstract

In this paper, we study the first eigenvalue of the magnetic Laplacian with Neumann boundary conditions in the unit disk $\mathbb D$ in $\mathbb R^2$. There is a rather complete asymptotic analysis when the constant magnetic field tends to $+\infty$ and some inequalities seem to hold for any value of this magnetic field, leading to rather simple conjectures. Our goal is to explore these questions by revisiting a classical picture of the physicist D. Saint-James theoretically and numerically. On the way, we revisit the asymptotic analysis in light of the asymptotics obtained by Fournais-Helffer, that we can improve by combining them with a formula stated by Saint-James.

Figures

Figures reproduced from arXiv: 2411.11721 by the authors.

Figure 1
Figure 1. Plot of β 7→ η(n, β) = λ(n,β) β for n from 0 to 20. The black dotted line and the black dashed line are the constants 1 and Θ0 respectively. The constants C β n and Dβ n , depending on n and β, are chosen so that Z 1 0 fn,β(r) 2 r dr = 1. The eigenvalues of (2.1) (and in particular the first eigenvalue, λ(n, β)) are the roots of the equation (fn,β) ′ (1) = 0 (2.5) (where the left-hand side is understood as the funct… view at source ↗
Figure 2
Figure 2. Saint-James picture recovered by M. Persson-Sundqvist [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Zoom on the plot of β 7→ η(n, β) for n from 0 to 20. The black dashed line is the constant Θ0. The black dots mark the points (β (1) n , η∗ n (1)). Let us recall that the sequence γn := βn+1−βn converges to 2 and, more precisely, admits 28 [PITH_FULL_IMAGE:figures/full_fig_p028_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. On the magnetic Dirichlet to Neumann operator on the disk -- strong diamagnetism and strong magnetic field limit--

    math.AP 2024-11 conditional novelty 7.0 of 10

    For the unit disk with a constant magnetic field of strength 2b, the ground state energy of the magnetic Dirichlet-to-Neumann operator satisfies λ_DN(b)=α√b-(α²+2)/6+O(b^{-1/2}) with α=0.76495..., the unique negative ...

  2. Three-dimensional magnetic Schr\"odinger operator with the potential supported in a tube

    math.SP 2025-05 reject novelty 4.0 of 10

    For a 3D magnetic Schrödinger operator with potential on a locally deformed tube, the paper proves essential spectrum stability and gives an insufficiently supported sufficient condition for the absence of discrete spectrum.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages · cited by 2 Pith papers

  1. [13]

    Fournais and B

    S. Fournais and B. Helffer. On the third critical field in Ginzburg-Landau theory. Comm. Math. Phys. 266, no. 1, 153–196 (2006)

  2. [17]

    Fournais and M

    S. Fournais and M. Persson. Strong diamagnetism for the ball in three dimensions. Asymptotic analysis, Vol. 72, Issue 1–2, 77–123 (2011)

  3. [1]

    Bauman, D

    P. Bauman, D. Phillips, and Q. Tang. Stable Nucleation for the Ginzburg-Landau system with an applied magnetic field. Arch. Rational Mech. Anal. 142, 1–43 (1998)

  4. [2]

    Baur and T

    M. Baur and T. Weidl. Eigenvalues of the magnetic Dirichlet Laplacian with constant magnetic field on disks in the strong field limit. Anal. Math. Phys. 15, art. 9 (2025)

  5. [3]

    Bernoff and P

    A. Bernoff and P. Sternberg. Onset of superconductivity in decreasing fields for general domains. J. Math. Phys. 39, 1272–1284 (1998)

  6. [4]

    Bolley and B

    C. Bolley and B. Helffer. An application of semi-classical analysis to the asymptotic study of the supercooling field of a superconducting material. Annales de l’I.H.P. Physique th´ eorique 58 (2), 189–233 (1993)

  7. [5]

    Bonnaillie

    V. Bonnaillie. Analyse math´ ematique de la supraconductivit´ e dans un domaine ` a coins : m´ ethodes semi-classiques et num´ eriques. Th` ese de Doctorat, Universit´ e Paris 11 (2003). https://theses.hal.science/tel-00005430

  8. [6]

    Colbois, C

    B. Colbois, C. L´ ena, L. Provenzano and A. Savo. Geometric bounds for the magnetic Neumann eigenvalues in the plane. Journal de Math´ ematiques Pures et Appliqu´ ees 179, 454–497 (2023)

Show all 30 references
  1. [7]

    Colbois, C

    B. Colbois, C. L´ ena, L. Provenzano and A. Savo. A reverse Faber-Krahn inequality for the magnetic Laplacian. Journal de Math´ ematiques Pures et Appliqu´ ees 192, 103632 (2024). 31

  2. [8]

    Dauge and B

    M. Dauge and B. Helffer. Eigenvalues Variation. I. Neumann Problem for Sturm- Liouville Operators. Journal of Differential Equations 104 (2), 243–262 (1993)

  3. [9]

    Dauge and B

    M. Dauge and B. Helffer. Eigenvalues Variation. II. Multidimensional Problems. Journal of Differential Equations 104 (2), 263–297 (1993)

  4. [10]

    de Gennes

    P.-G. de Gennes. Boundary effects in superconductors. Rev. Mod. Phys. 36, 225 (1964)

  5. [11]

    NIST Digital Library of Mathematical Functions , Release 1.2.0 of 2024-03-15. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, eds. NIST (2024). https://dlmf.nist.gov/

  6. [12]

    Fournais and B

    S. Fournais and B. Helffer. Accurate eigenvalue asymptotics for the magnetic Neumann Laplacian. Ann. Inst. Fourier 56 (1), 1–67 (2006)

  7. [14]

    Fournais and B

    S. Fournais and B. Helffer. Strong diamagnetism for general domains and applications. Annales de l’Institut Fourier, Tome 57, no. 7, 2389–2400 (2007)

  8. [15]

    Fournais and B

    S. Fournais and B. Helffer. Spectral Methods in Surface Superconductivity. Progress in Nonlinear Differential Equations and Their Applications, Vol. 77, Birkh¨ auser (2010)

  9. [16]

    Fournais and B

    S. Fournais and B. Helffer. Inequalities for the lowest magnetic Neumann eigenvalue. Letters in Mathematical Physics, Volume 109, pages 1683–1700 (2019)

  10. [18]

    Fournais and M

    S. Fournais and M. Persson Sundqvist. Lack of diamagnetism and the Little–Parks effect. Commun. Math. Phys. 337, 191–224 (2015)

  11. [19]

    Fournais and M

    S. Fournais and M. Persson Sundqvist. Calculations for the disc. Unpublished (2017). Personal communication

  12. [20]

    Helffer and A

    B. Helffer and A. Morame. Magnetic bottles in connection with superconductivity. J. Funct. Anal. 185 (2), 604–680 (2001)

  13. [21]

    Helffer and A

    B. Helffer and A. Morame. Magnetic bottles for the Neumann problem: curvature effects in the case of dimension 3 (general case). Ann. Sci. ´Ec. Norm. Sup´ er. (4) 37 (1), 105–170 (2004)

  14. [22]

    Helffer and M

    B. Helffer and M. Persson Sundqvist. On nodal domains in Euclidean balls. Proc. Amer. Math. Soc. 144, no. 11, 4777–4791 (2016). 32

  15. [23]

    Kachmar, V

    A. Kachmar, V. Lotoreichik and M. Sundqvist. On the Laplace operator with a weak magnetic field in exterior domains. Anal. Math. Phys. 15, art. 5 (2025)

  16. [24]

    Kachmar and G

    A. Kachmar and G. Miranda. The magnetic Laplacian on the disc for strong magnetic fields. J. Math. Anal. Appl. 546 (2), 129261 (2025)

  17. [25]

    Lu and X.B

    K. Lu and X.B. Pan. Eigenvalue problems of Ginzburg-Landau operator in bounded domains. J. Math. Phys. 40 (6), 2647–2670 (1999)

  18. [26]

    Magnus, F

    W. Magnus, F. Oberhettinger and R.P. Soni. Formulas and theorems for the special functions of mathematical physics, 3rd enlarged ed., Grundlehren der Mathematischen Wissenschaften, Volume 52, Springer (1966)

  19. [27]

    Persson Sundqvist

    M. Persson Sundqvist. Magnetic model operators. A short review and something new. Lecture at conference in honour of the 70th birthday of Bernard Helffer. April 2019

  20. [28]

    Saint-James

    D. Saint-James. Etude du champ critique HC3 dans une g´ eom´ etrie cylindrique. Physics Letters 15 (1), 13–15 (1965)

  21. [29]

    Saint-James and P-G

    D. Saint-James and P-G. de Gennes. Onset of superconductivity in decreasing fields. Physics Letters 7 (5), 306–308 (1963)

  22. [30]

    Soojin Son

    S. Soojin Son. Spectral problems on triangles and discs: extremizers and ground states. PhD thesis, University of Illinois Urbana-Champaign (2014). https://www.ideals. illinois.edu/items/49400. 33

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.