Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

P\'olya's Conjecture for the Neumann Eigenvalues on Euclidean Balls

T0 review · 3 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

Pith's one-line read This paper claims a proof of Pólya's 1954 Neumann counting conjecture for every Euclidean ball: for all d≥2, R>0, E≥0, the number of Neumann eigenvalues below E is at least the Weyl term (R√E)^d/(2^d Γ(d/2+1)^2).

desk verdict A serious, well-architected proof of Pólya's Neumann-ball conjecture in all dimensions whose only load-bearing soft spot is an unverified exact-rational appendix — worth refereeing, conditional on CAS verification. read the letter →

arxiv 2607.25958 v3 pith:XXLNRZBF submitted 2026-07-28 math.SP math.APmath.CA

classification math.SPmath.APmath.CA MSC 35P1533C1011P21
keywords NeumannLaplacianPólyainequalityEuclideanballBesselzerosDiniboundaryconditionvariationalprinciplespectralcountingfunctionsWeylterm
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

The paper claims to settle Pólya's 1954 conjecture for the Neumann Laplacian on Euclidean balls: in every dimension d≥2, at every radius, and at every energy, the number of Neumann eigenvalues below that energy is at least the Weyl term with the sharp constant. This is a one-sided spectral bound that holds at all energies, not just asymptotically, and it has the correct leading-order coefficient. The proof separates variables into angular sectors, uses a strict comparison between the physical Neumann boundary condition and a Robin problem whose frequencies are Bessel derivative zeros, and then reduces the problem to comparing a weighted staircase sum with the Weyl integral. Low frequencies are handled by variational trial spaces, intermediate frequencies by finite radial-level estimates and beta-moment bounds, and high frequencies by an explicit scalar tail criterion; the disk in dimension two is treated separately. If correct, the ball case of Pólya's Neumann conjecture is closed in every dimension.

What carries the argument

The load-bearing bridge is the strict Robin transfer: for d≥3 it turns physical Neumann eigenvalues into strictly smaller Robin eigenvalues with frequencies j'_{ν,k} (zeros of the derivative of a Bessel function), so a phase estimate for those zeros yields a lower bound for the strict physical count. In frequency variables the bound becomes a weighted quarter-shifted floor sum over the inverse of the action function G_x(z)=(√(x²-z²)-z arccos(z/x))/π, whose exact integral is the Weyl term; comparing that sum with the integral is the task carried out by variational trial spaces, a four-level radial estimate, beta-moment convex quadrature, and the scalar tail criterion.

What would settle it

Run an independent exact-arithmetic verification of the appendix's finite inequalities, starting with Proposition C.1: check whether P_{3,2}(t)-1, P_{3,3}(t)-1, and P_{3,4}(t)-1 are positive on the stated rational t-intervals, and check the table residuals in (C.4) and (H.11)–(H.14). A single negative value would refute the intermediate comparison and with it Theorem 1.1. A complementary check is to evaluate N_d^<(x)-W_d(x) on a fine grid of x for d=2,3,4,5 using the phase-sum representation and the strict Robin lower bound; any negative difference would disprove the paper's central claim. The

Watch

Extended reading notes

Core claim

The central claim is Theorem 1.1. For the unit ball, separation of variables shows that the physical Neumann radial boundary condition is not J'_ν(k)=0 but the Dini equation k J'_ν(k)-δ_d J_ν(k)=0 with δ_d=(d-2)/2. The paper proves a strict Robin comparison: each physical Neumann eigenvalue lies strictly below the corresponding Robin eigenvalue whose frequencies are the positive zeros of J'_ν. Feeding a published phase estimate for those zeros into a sum over angular multiplicities gives a lower bound P_d(x) for the strict Neumann count N_d^<(x). The rest of the proof is a comparison of that weighted staircase with the exact inverse-action integral W_d(x), carried out through variational sta

Load-bearing premise

The theorem stands or falls on the correctness of the appendix's long chains of exact rational inequalities that close the intermediate frequency ranges—for example, the claims that the polynomials P_{3,2}, P_{3,3}, P_{3,4} exceed 1 on their respective bands—together with the validity of the single imported phase estimate for zeros of J'_ν.

Editorial extensions

If this is right

  • For every n≥1, the (n+1)-st Neumann eigenvalue of a ball of radius R satisfies μ^N_{n+1}(B_R^d) ≤ 4π² n^{2/d}/(ω_d |B_R^d|)^{2/d}.
  • Any bounded Lipschitz domain that tiles a ball by finitely many congruent copies inherits the Pólya Neumann lower bound; this includes half-balls and orthant sectors.
  • The bound is valid with the sharp Weyl constant at all frequencies, including the low- and middle-frequency ranges that are invisible to two-term Weyl asymptotics.
  • The strict counting convention at equality faces means the inequality remains true when an eigenvalue sits exactly on the bound.
  • The proof supplies explicit rational thresholds where each of the four methods takes over, making the all-energy statement fully explicit in every dimension.

Reading between the lines

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

  • If the proof is verified, the Neumann-ball version of Pólya's conjecture is closed in all dimensions; the disk case was already known, and dimensions d≥3 were the open part.
  • The same strict Robin transfer could be attempted on other rotationally symmetric domains, such as spherical shells or sectors, where angular multiplicities are explicit, though the weighted-staircase-versus-integral comparison would need a new argument.
  • The appendix's exact rational inequalities are independently checkable by exact-arithmetic computation; turning them into machine-verified certificates would remove the main residual doubt quickly.
  • A natural testable extension is whether the method degrades gracefully for nearly spherical domains, where the Bessel phase structure is lost but a perturbative version of the staircase might still yield a Pólya-type bound with an explicit error.
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

3 major / 4 minor

Summary. The paper claims to prove Pólya's conjectured Neumann lower bound for every Euclidean ball in every dimension d ≥ 2. The proof combines a strict Robin comparison (Prop. 2.1) with the published FLPS derivative-zero phase estimate to reduce the Neumann counting function to a weighted sum of shifted Bessel-phase floors; it then compares this discrete sum to the Weyl inverse-action integral. The comparison is carried out range-by-range: a uniform scalar high-frequency criterion (§2.4), a correlated dimension-three tail (Thm. 4.1), a four-dimensional beta-moment/phase estimate (§5), and for d ≥ 5 a low-frequency variational staircase, a four-level estimate, a beta-moment estimate, and a scalar tail at 5d^{3/2}. The compact-range inequalities are stated as exact rational assertions in a long proof appendix (Props. C.1, D.1, F.1, G.1–G.3, H.1, I.1 and associated tables). If correct, the result proves the Neumann-ball case of Pólya's conjecture and gives the eigenvalue corollary (1.2).

Significance. The main theorem is a major result if the proof is correct: it settles a prominent open case of Pólya's 1954 conjecture, for which only the Dirichlet-ball and planar Neumann cases were previously known by the FLPS work. The global architecture is strong: the strict Robin transfer, the exact discrete phase-sum representation (2.11)–(2.14), the inverse-action identity (2.18), and the scalar tail criterion are mathematically coherent and, as far as I could check, internally consistent. The paper also contains several genuinely useful auxiliary estimates, including the quarter-shift quadrature lemma, exact inverse moments, and the correlated dimension-three inequality. However, the final acceptance hinges on an extensive layer of exact-rational inequalities that is not independently verifiable from the submitted text alone; that gap is load-bearing, because every finite frequency window depends on it.

major comments (3)
  1. [§H.2 and proof appendix (Tables C.4, D.3, F.12, H.11–H.14)] The finite exact-rational layer is the principal load-bearing component of the proof, yet it is not independently checkable as submitted. The text in §H.2 references the file verification/generated/dge7_aggregate_finite.txt as recording exact curvature residuals, but that file is not part of the preprint. The displayed tables are too large for hand audit. A spot check near (6.48)–(6.49) did not reproduce the claimed rational residual; the signs survived, but the exact numbers did not. Since every intermediate frequency band is closed by one of Props. C.1, D.1, E.1, F.1, G.1–G.3, H.1, I.1, a single arithmetic slip in any of these tables can invalidate the corresponding band and hence Theorem 1.1. The author should supply a machine-checkable certificate or a CAS-verified transcript for the entire exact-rational appendix, or independently prove the unreferenced tables within the text.
  2. [Abstract vs. appendix, 'No numerical approximation or executable certificate is used'] The abstract says no executable certificate is used as a premise, but the appendix itself refers to a generated verification file for exact residuals. The manuscript therefore appears to rely on an external finite computation that is neither included nor described. This is not a numerical-approximation issue, but it is a verification gap: the displayed tables may be the output of a program, and without the program or its certified output the reader cannot confirm them. At minimum, the statement in the abstract should be qualified, and the missing file should be distributed or replaced by independently checkable certificates.
  3. [Proposition 2.1/§2.1.2, use of the FLPS phase estimate] The only imported nonstandard inequality is [1, Prop. 3.1], used at (2.10) and Lemma 3.1. This is a legitimate published result and I do not see circularity. The strict Robin transfer is carefully argued and I verified the key monotonicity direction. My concern is narrower: the phase estimate [1, Prop. 3.1] is stated as an external theorem, and the present proof inherits all of its endpoint and convention subtleties. Given that the entire high-frequency argument rests on this estimate, the author should spell out precisely how the endpoint convention in [1] aligns with the strict/non-strict counting used here, especially in the Appendix J disk argument where a non-strict count is converted to a strict count by a limiting procedure.
minor comments (4)
  1. [§J.3] Typo: 'convenient convenient value' should read 'convenient value'.
  2. [§1.4/§2.3] The notation α_n is introduced in §2.3 and again used in §3 and §6 with the same meaning; this is acceptable but a one-line reminder in the later sections would improve readability.
  3. [Table 1] The table lists 'x=20/3' for d=2, but the finite staircase extends to 2√12; the text later explains the overlap. A short remark in the caption would prevent confusion.
  4. [Appendix H, §H.2] The deterministic rule defining k_d is terse; please include one worked example (e.g., d=7) showing how the rational endpoints e_d^± are obtained from k_d, so the reader can verify the rule without reverse-engineering it.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the Neumann-ball Pólya inequality is derived from an external phase estimate and self-contained reductions; the target is never assumed.

full rationale

The derivation chain is self-contained apart from the published derivative-zero phase estimate [1, Prop. 3.1] (Filonov–Levitin–Polterovich–Sher, external to the author). That estimate is used at (2.10) and Lemma 3.1 to lower-bound the number of derivative zeros; the paper's own Proposition 2.1 (strict Robin transfer) converts this to a strict physical Neumann count via (2.8); (2.11)–(2.14) is an exact level-set identity; (2.18) is the exact inverse-action representation of the Weyl term. The target inequality N^< ≥ W is never assumed as an input. The remaining finite-window verifications (Props. C.1–I.1) are asserted exact rational inequalities; the missing certificate and a failed spot-check are a verification/correctness gap of the proof appendix, not a circular dependency. The paper does not fit parameters to the predicted quantity, does not rename a known result, and contains no load-bearing self-citation: all cited work is by other authors. Hence no circular step can be exhibited, and the appropriate score is 0.

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

The theorem statement contains no free parameters, and the proof introduces no invented objects (no new forces, particles, or conserved quantities). All constants (c₀, α_n, β_d, D_d, h_d, γ) are derived from π and Γ via exact identities. The listed free parameters are proof-internal frequency-range partition points whose endpoint inequalities are (purportedly) verified exactly; they do not enter the result and none is fitted to numerical data. The single significant external input is the published FLPS phase estimate, acknowledged explicitly in §1.3 and the appendix.

free parameters (5)
  • d=2 range cutoffs x=20/3 and 2√12
    Boundary between the finite Ritz staircase and the analytic tail in the disk proof (§3.2–§3.4); chosen so that (20/3)²=400/9<48 and the tail margin closes at 61/3000. Proof-internal; verified exactly.
  • d=3 cutoffs x=5 and x=14
    Junctions of the variational staircase, the three truncated-sum bands (5≤x≤14, Prop. 4.4), and the correlated tail (x≥14) (§4).
  • d=4 cutoffs x=38/5, 26, 42
    Junctions of the enlarged Ritz range, the phase bands, the beta-moment band, and the scalar tail (§5); the tail entry point 42 is where (2.67) closes with margin ≈43/300−0.1427.
  • d≥5 cutoffs a_d^- (41/50 for d=5, 4/5 for d≥6), 5/2, 3, 5, and tail 5d^{3/2}
    Uniform dimension-scaled partition of the frequency axis (§6). a_5^- is chosen because the four-level estimate provably fails at a=4/5 (Remark after Thm. 6.2 and (F.1): L_{5,4}(4/5)<999/1000). All proof-internal; not in the theorem statement.
  • Explicit scalar tail threshold X_p=(24p)^{3/2} (p=d−1), sharpened to 5d^{3/2}
    Closed-form cutoff for the scalar high-frequency criterion (2.32) and Prop. 2.4; the sharper 5d^{3/2} entry is proved by the two-parity induction (I.1).
assumptions (5)
  • domain assumption FLPS derivative-zero phase estimate: for ν>0, #{k: j'_{ν,k} ≤ x} ≥ ⌊G_x(ν)+3/4⌋ with G_x as in (2.9)
    Imported from Filonov–Levitin–Polterovich–Sher [1, Prop. 3.1, eq. (3.3)]; used at (2.10) and Lemma 3.1 for orders ν = ℓ+(d−2)/2 ≥ 1/2. Published (Invent. Math. 234, 2023) and peer-reviewed; external to the present author, so treated as independent support.
  • standard math Separation of variables gives the full Neumann spectrum of the ball from angular sectors with the Dini condition kJ'_ν(k) − δ_d J_ν(k) = 0
    §2.1.1; regular solutions r^{-δ_d}J_{ν_{d,ℓ}}(kr); standard Bessel identities and zero conventions per DLMF [15] / Watson [16].
  • standard math Min–max principle; Friedrichs form-domain theory at the singular endpoint; analytic type-(B) family of Robin forms with positive Hellmann–Feynman derivative
    Proposition 2.1 and its proof; Kato [12, Ch. VII §4]. Used to obtain the strict bound λ^N_{ℓ,k} < (j'_{ν,k})², the load-bearing spectral transfer.
  • standard math Log-convexity of the Gamma function; beta–Gamma duplication identities; monotonicity of the ratios h_d and D_d along parity classes
    Used throughout: exact inverse moments (2.52)–(2.54), dimension factors (2.77), Gamma estimates (6.19)–(6.21), (H.38), and the tail induction (I.3).
  • standard math Zero-counting and endpoint conventions for Bessel derivatives in the disk case (j'_{0,1}=0, j'_{0,k+1}=j_{1,k}; formal zeros J'_m(0)=0 for m≥2 do not create eigenmodes)
    §3.1, (3.3)–(3.4). Needed for the planar phase transfer and for the strict/non-strict counting distinction in d=2.

how reviews work

0 comments
Cite this review

Pith. "Pith review of P\'olya's Conjecture for the Neumann Eigenvalues on Euclidean Balls." pith.science (2026). https://pith.science/paper/XXLNRZBF

@misc{pith2026260725958,
  author       = {Pith},
  title        = {Pith review of: P\'olya's Conjecture for the Neumann Eigenvalues on Euclidean Balls},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XXLNRZBF}},
  note         = {Machine review of arXiv:2607.25958}
}
abstract

We prove P\'olya's conjectured lower bound for the Neumann eigenvalue counting function of Euclidean balls. If $B_R^d\subset\mathbb R^d$ is the ball of radius $R$, then, for every $d\ge2$, $R>0$, and $E\ge0$, $$ N_{B_R^d}^{<}(E) \ge \frac{\omega_d}{(2\pi)^d}|B_R^d|E^{d/2} = \frac{(R\sqrt E)^d}{2^d\Gamma(\frac d2+1)^2}, $$ where $\omega_d$ is the volume of the unit $d$-ball and $N_{B_R^d}^{<}(E)$ counts Neumann eigenvalues strictly below $E$. Combined with the Dirichlet theorem for balls, this settles both P\'olya inequalities for Euclidean balls in every dimension $d\ge2$. In the disk case, the proof replaces a computer-assisted finite-frequency step by explicit Rayleigh--Ritz estimates. In dimensions $d\ge3$, the radial Neumann condition is a Dini condition rather than a derivative-zero Bessel condition. A strict comparison with an auxiliary Robin problem transfers a derivative-zero Bessel phase estimate to the physical Neumann spectrum. The problem then becomes a comparison between a multiplicity-weighted phase staircase and an integral equal to the Weyl term. Variational trial spaces control low frequencies; finitely many radial levels and beta-integral estimates cover the intermediate range; and a uniform phase estimate treats high frequencies. All finite computations for $2\le d\le6$ are printed in the paper. For $d\ge7$, one compact two-parameter estimate is verified in exact rational arithmetic by the ancillary program.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. P\'{o}lya's conjecture for higher-dimensional Neumann balls

    math.SP 2026-07 accept novelty 7.0 of 10

    For every d≥3 and every λ≥0, the Neumann counting function of the unit d-ball satisfies N(λ) ≥ w_d λ^d, so Pólya's conjecture holds for all Euclidean balls.

Pith tools

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