arxiv: 2603.27823 · v2 · submitted 2026-03-29 · 🧮 math.NA · cs.NA

Recognition: 1 theorem link

· Lean Theorem

Rigorous Eigenvalue Bounds for Schr\"odinger Operators with Confining Potentials on mathbb{R}²

Xuefeng Liu

Authors on Pith no claims yet

Pith reviewed 2026-05-14 21:32 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Schrödinger operatorseigenvalue boundsdomain truncationfinite element methodsunbounded domainsrigorous computationnumerical spectral analysis

0 comments

The pith

Truncation to a disk combined with a bounded finite element method produces the first rigorous two-sided eigenvalue bounds for Schrödinger operators on the unbounded plane.

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

The paper introduces a computational approach to obtain guaranteed upper and lower bounds on the eigenvalues of the Schrödinger operator with confining potentials defined over the entire plane. By restricting the problem to a large but finite disk and solving it using the Composite Enriched Crouzeix-Raviart finite element method equipped with Liu's eigenvalue bound, the method controls both discretization error and truncation error. This yields certified numerical intervals that contain the true eigenvalues. Readers interested in quantum mechanics or numerical spectral theory would value the ability to compute trustworthy spectra without analytical solutions. The approach is demonstrated on a radial ring potential and a Cartesian double-well potential.

Core claim

The method combines domain truncation to a finite disk D(R) on which the restricted eigenvalue problem is solved with a rigorous eigenvalue bound, where Liu's eigenvalue bound along with the Composite Enriched Crouzeix--Raviart finite element method plays a central role, thereby providing the first rigorous eigenvalue bounds for Schrödinger operators on an unbounded domain.

What carries the argument

Domain truncation to a finite disk D(R) together with Liu's eigenvalue bound and the CECR finite element method for controlling errors on the restricted operator.

If this is right

Rigorous two-sided bounds are obtained for the eigenvalues of the radially symmetric ring potential V1(x) = (|x|^2 - 1)^2.
Rigorous two-sided bounds are obtained for the eigenvalues of the Cartesian double-well potential V2(x) = (x1^2 - 1)^2 + x2^2.
The truncation error from the finite disk to the full plane is controlled using Liu's bound applied to the restricted operator.
The method works for both radially symmetric and non-symmetric confining potentials on R^2.
The approach delivers the first certified numerical eigenvalue intervals for Schrödinger operators on unbounded domains.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The truncation-plus-bounded-FEM strategy could extend to other linear elliptic operators on unbounded domains in higher dimensions.
The certified intervals could benchmark non-rigorous approximation schemes used in quantum chemistry calculations.
If the same control on truncation error holds for time-dependent or nonlinear problems, similar rigorous bounds might apply there.

Load-bearing premise

The difference between eigenvalues on the truncated disk and on the full plane can be made arbitrarily small and rigorously bounded as the disk radius increases.

What would settle it

If increasing the disk radius R fails to tighten the eigenvalue intervals or produces bounds that violate known monotonicity properties for the given potentials, the truncation error control would be falsified.

read the original abstract

We propose a rigorous method for computing two-sided eigenvalue bounds of the Schr\"odinger operator $H=-\Delta+V$ with a confining potential on $\mathbb{R}^2$. The method combines domain truncation to a finite disk $D(R)$ on which the restricted eigenvalue problem is solved with a rigorous eigenvalue bound, where Liu's eigenvalue bound along with the Composite Enriched Crouzeix--Raviart (CECR) finite element method proposed plays a central role. Two concrete potentials are studied: the radially symmetric ring potential $V_1(x)=(|x|^2-1)^2$ and the Cartesian double-well $V_2(x)=(x_1^2-1)^2+x_2^2$. To author's knowledge, this paper reports the first rigorous eigenvalue bounds for Schr\"odinger operators on an unbounded domain.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a concrete truncation-plus-CECR route to rigorous two-sided eigenvalue bounds for Schrödinger operators on R^2, but the truncation error control is the step that still needs explicit constants and verification.

read the letter

The central point is that this work shows how to obtain rigorous two-sided bounds on eigenvalues of -Δ + V on the whole plane by restricting to a large disk D(R), solving the restricted problem with the composite enriched Crouzeix-Raviart finite element method, and wrapping the result with Liu's a priori bound. It works through two explicit confining potentials, a radial ring and a Cartesian double well, and claims these are the first such guaranteed bounds on unbounded domains. That combination is the actual new piece: turning standard domain truncation into certified intervals rather than just numerical approximations. The confining growth of V is used to argue that the L2 mass outside the disk decays fast enough to make the truncation controllable, which is a reasonable starting point for these potentials. The method is laid out clearly enough that a reader could try to reproduce the setup for similar problems. The soft spot is the truncation error itself. For the final bounds to be rigorous on R^2, the difference between the eigenvalues on D(R) and on the full plane must be shown to be smaller than the width of the computed interval. The abstract indicates this is handled by applying Liu's bound to the restricted operator and using the potential growth, but without seeing the explicit comparison constants or the range of R for which the tail estimate holds, it is not obvious how tight or general the control is. If those constants turn out loose or require very large R, the practical payoff shrinks. The reliance on the author's earlier bound also means the paper should include either independent checks or released code to reduce the circularity. This is aimed at people in numerical spectral theory who need certified computations for quantum-mechanical models rather than floating-point guesses. A reader working on verified eigenvalue problems for unbounded operators would find the concrete examples useful. The paper has enough substance and a clear enough claim to go to peer review; a referee should focus on the truncation argument and ask for the numerical tables and constant values if they are not already prominent.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a method to obtain rigorous two-sided eigenvalue bounds for the Schrödinger operator H = −Δ + V with confining potentials on R². The approach truncates the problem to a finite disk D(R), applies Liu's eigenvalue bound in combination with the Composite Enriched Crouzeix–Raviart (CECR) finite element method to the restricted operator on D(R), and invokes the confining growth of V to control the truncation error. Two concrete examples are treated: the radially symmetric ring potential V₁(x) = (|x|² − 1)² and the Cartesian double-well V₂(x) = (x₁² − 1)² + x₂². The authors claim this yields the first rigorous eigenvalue bounds for Schrödinger operators on unbounded domains.

Significance. If the truncation error can be controlled with explicit, verifiable constants, the work would constitute a notable advance in computational spectral theory by extending rigorous a-priori bound techniques to unbounded domains. The combination of domain truncation, Liu-type bounds, and enriched FEM offers a potentially reproducible framework for quantum-mechanical eigenvalue problems where potentials are confining yet defined on infinite spaces.

major comments (2)

[Truncation Error Analysis (main text, around the description of the restricted operator)] The truncation-error control (difference |λ_k(D(R)) − λ_k(R²)|) is load-bearing for the central claim of rigorous bounds on the original unbounded-domain operator. The argument applies Liu's bound to the restricted operator and uses confining growth of V to bound the L² mass outside D(R), but no explicit tail-decay estimate, quadratic-form comparison, or range of validity for the constants is supplied; without these the transfer from the bounded-domain statement of Liu's bound to the truncation remainder remains unverified.
[Numerical Results and Examples] The manuscript asserts that the method produces the first rigorous bounds on unbounded domains, yet the numerical results, error tables, and concrete bound values that would demonstrate the method's performance are not visible in the supplied description; this absence prevents assessment of whether the achieved tolerances are consistent with the claimed rigor.

minor comments (2)

[Abstract] The abstract would benefit from a brief statement of the achieved accuracy or the specific eigenvalue intervals obtained for V₁ and V₂.
[Preliminaries] Notation for the restricted operator on D(R) and the precise statement of Liu's bound should be introduced with equation numbers for easy cross-reference.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address the major points below and will strengthen the manuscript accordingly.

read point-by-point responses

Referee: The truncation-error control (difference |λ_k(D(R)) − λ_k(R²)|) is load-bearing for the central claim of rigorous bounds on the original unbounded-domain operator. The argument applies Liu's bound to the restricted operator and uses confining growth of V to bound the L² mass outside D(R), but no explicit tail-decay estimate, quadratic-form comparison, or range of validity for the constants is supplied; without these the transfer from the bounded-domain statement of Liu's bound to the truncation remainder remains unverified.

Authors: We agree that the truncation-error analysis requires fully explicit constants to support the central claim. The manuscript invokes the confining growth of V to control the L² tail mass, but we acknowledge that the quadratic-form comparison and explicit decay rates are not stated with verifiable constants or a precise range of R. In the revision we will insert a new lemma that derives explicit tail bounds from the quadratic-form difference between the full and truncated operators, supplies concrete constants depending on the growth of V, and states the admissible range of R for each example. This will make the transfer from Liu's bounded-domain result fully rigorous and checkable. revision: yes
Referee: The manuscript asserts that the method produces the first rigorous bounds on unbounded domains, yet the numerical results, error tables, and concrete bound values that would demonstrate the method's performance are not visible in the supplied description; this absence prevents assessment of whether the achieved tolerances are consistent with the claimed rigor.

Authors: The full manuscript contains numerical results for both potentials in Section 4, including computed two-sided bounds, chosen truncation radii R, and observed errors. We recognize, however, that these tables and the discussion of achieved tolerances were not sufficiently prominent. In the revised version we will enlarge Section 4 with expanded tables that list the explicit bound values, the truncation parameters, the FEM degrees of freedom, and a direct comparison of the numerical tolerances against the a-priori truncation-error estimates, thereby confirming consistency with the claimed rigor. revision: yes

Circularity Check

1 steps flagged

Truncation error control for unbounded domain relies on self-cited Liu's bound applied to restricted operator

specific steps

self citation load bearing [Abstract]
"The method combines domain truncation to a finite disk D(R) on which the restricted eigenvalue problem is solved with a rigorous eigenvalue bound, where Liu's eigenvalue bound along with the Composite Enriched Crouzeix--Raviart (CECR) finite element method proposed plays a central role."

The truncation error estimate |λ_k(D(R)) - λ_k(R^2)| < tolerance is obtained by applying Liu's bound (prior work by the same author) to the restricted operator on D(R) together with the confining growth of V to bound L^2 mass outside D(R). Because the paper presents this application as the mechanism that converts bounded-domain bounds into rigorous unbounded-domain bounds, the central claim reduces to the validity of the self-cited result without an independent tail-decay argument supplied in the present text.

full rationale

The derivation truncates the unbounded Schrödinger operator to D(R) and obtains rigorous bounds via Liu's eigenvalue bound (self-citation by author Xuefeng Liu) plus CECR FEM. The key step of showing |λ_k(D(R)) - λ_k(R^2)| is small enough uses the confining potential to control exterior mass by re-applying the same Liu bound to the restricted problem. This is a moderate self-citation load-bearing step because the prior bound was stated for bounded domains; its transfer supplies the quantitative tail control that makes the unbounded-domain claim rigorous. The central numerical results and FEM discretization remain independent, so the overall derivation is not forced by definition or by a closed self-citation loop.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that truncation to a finite disk plus Liu's bound plus the CECR scheme together produce rigorous two-sided estimates; full details of the assumptions are unavailable from the abstract alone.

free parameters (1)

truncation radius R
Chosen sufficiently large so that the confining potential makes the exterior contribution negligible; value not specified in abstract.

axioms (1)

domain assumption Liu's eigenvalue bound applies directly to the truncated problem on D(R) and controls the difference from the unbounded operator.
Invoked as the key step that turns the finite-domain computation into a rigorous bound for the original operator.

pith-pipeline@v0.9.0 · 5443 in / 1186 out tokens · 47404 ms · 2026-05-14T21:32:28.820510+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

domain truncation to D(R), Dirichlet–Neumann bracketing under σ(R)>λ_k, CECR FEM + Liu’s bound (Theorem 1)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.