2025-11-26T04:46:18.614838

Low-lying eigenvalues in the semiclassical limit of a Schr{ö}dinger operator with an inverse square potential, and non-asymptotic a-zeros of Kummer functions

Vanlaere
We provide a precise description of the bottom of the spectrum in the semiclassical limit of a harmonic-type Schr{ö}dinger operator with an inverse square potential. By exploiting the connection between the eigenfunctions of these operators and the Kummer and Whittaker functions, we derive accurate localization results for the non-asymptotic zeros of these functions with respect to their first parameter, uniformly with respect to the argument taken large and real. The spectral analysis relies on a WKB-type approach.
academic

Low-lying eigenvalues in the semiclassical limit of a Schrödinger operator with an inverse square potential, and non-asymptotic a-zeros of Kummer functions

Basic Information

  • Paper ID: 2511.20025
  • Title: Low-lying eigenvalues in the semiclassical limit of a Schrödinger operator with an inverse square potential, and non-asymptotic a-zeros of Kummer functions
  • Author: Roman Vanlaere (CEREMADE, Université Paris-Dauphine PSL, CNRS)
  • Classification: math.CA (Classical Analysis), math.SP (Spectral Theory)
  • Submission Date: November 21, 2025
  • Paper Link: https://arxiv.org/abs/2511.20025

Abstract

This paper provides a precise description of the spectrum near the bottom of a harmonic-type Schrödinger operator with an inverse square potential in the semiclassical limit. By exploiting the connection between the eigenfunctions of these operators and Kummer and Whittaker functions, the author derives precise localization results for non-asymptotic zeros of these functions with respect to their first parameter, uniformly for large real values of the argument. The spectral analysis relies on WKB-type methods.

Research Background and Motivation

Core Problem

The paper studies the spectral behavior in the semiclassical limit ξ→+∞\xi \to +\infty of a Schrödinger operator with a singular inverse square potential defined on the bounded interval (0,1)(0,1): Gξ=−∂x2+ξ2x2+ν2−1/4x2G_\xi = -\partial_x^2 + \xi^2 x^2 + \frac{\nu^2 - 1/4}{x^2} where ξ>0\xi > 0 is the semiclassical parameter and ν≥0\nu \geq 0 is a fixed parameter.

Research Significance

  1. Theoretical Importance: The inverse square potential plays an important role in quantum mechanics, but spectral analysis on bounded intervals presents technical difficulties because the singularity at the boundary makes standard semiclassical analysis methods difficult to apply directly.
  2. Special Function Connection: The eigenfunctions of this operator can be expressed explicitly using Kummer functions (confluent hypergeometric functions), establishing a deep connection between spectral theory and the theory of zeros of special functions.
  3. Literature Gap: Existing literature lacks systematic research on the zeros of Kummer functions with respect to the first parameter aa (a-zeros). The known asymptotic behavior a∼−k2π2/4za \sim -k^2\pi^2/4z does not cover all solutions, and uniform estimates with respect to the argument are lacking.

Limitations of Existing Methods

  1. Difficulties in Semiclassical Analysis: Traditional WKB methods require the potential function to have sufficient regularity at the interval boundaries, but the inverse square potential is singular at x=0x=0.
  2. Known Results on Special Function Zeros:
    • z-zeros of Kummer functions (zeros with respect to the argument) have been studied more extensively
    • Zeros with respect to parameter aa have only imprecise asymptotic formulas or numerical approximations
    • Systematic characterization of non-asymptotic zeros is lacking
  3. Difference Between Full Line and Bounded Interval: When considered on the full line R+\mathbb{R}^+, the requirement that eigenfunctions be square-integrable forces the parameter a=−ka=-k to be a negative integer, causing Kummer functions to degenerate into polynomials (Laguerre polynomials), making eigenvalues explicitly available. However, on bounded intervals, boundary conditions change the nature of the problem.

Core Contributions

The main contributions of this paper include:

  1. Establishing Precise Correspondence Between Eigenvalues and Kummer Zeros (Theorems 1.1 and 1.2):
    • Proves that λ\lambda is an eigenvalue of GξG_\xi if and only if a=−(λ−2ξ(1+ν))/(4ξ)a = -(\lambda - 2\xi(1+\nu))/(4\xi) satisfies M(a,1+ν,ξ)=0M(a, 1+\nu, \xi) = 0
  2. Two-Sided Bounds for Eigenvalues (Theorem 1.3):
    • Lower bound: λξ,k/ξ>4k+2(1+ν)\lambda_{\xi,k}/\xi > 4k + 2(1+\nu) (strict inequality)
    • Improved lower bound for large kk: λξ,k≥ck2\lambda_{\xi,k} \geq ck^2 (with c∈(0,π2)c \in (0,\pi^2))
  3. Exponentially Precise Upper Bounds for Low-Lying Eigenvalues (Theorem 1.4):
    • For k≤⌊τξ/4⌋k \leq \lfloor \tau\xi/4 \rfloor (τ∈(0,1)\tau \in (0,1)): λξ,kξ≤4k+2(1+ν)+C1e−C2ξ\frac{\lambda_{\xi,k}}{\xi} \leq 4k + 2(1+\nu) + C_1 e^{-C_2\xi} This shows that low-lying eigenvalues stabilize to 4k+2(1+ν)4k + 2(1+\nu) at exponential rate
  4. Systematic Characterization of Kummer Function a-Zeros (Theorem 1.5):
    • Proves that zeros form a strictly decreasing sequence: ⋯<aξ,k<aξ,k−1<⋯<aξ,0<0\cdots < a_{\xi,k} < a_{\xi,k-1} < \cdots < a_{\xi,0} < 0
    • Establishes precise bounds for non-asymptotic zeros: for k≤⌊τξ/4⌋k \leq \lfloor \tau\xi/4 \rfloor, −k−C1e−C2ξ≤aξ,k<−k-k - C_1e^{-C_2\xi} \leq a_{\xi,k} < -k

Detailed Methodology

Problem Setup

Operator Definition: Consider the Friedrichs extension on L2(0,1)L^2(0,1):

D(G_\xi) = \{f \in H^1_{0,\nu}(0,1), G_\xi f \in L^2(0,1)\} \\ G_\xi = -\partial_x^2 + \xi^2 x^2 + \frac{\nu^2-1/4}{x^2} \end{cases}$$ where $H^1_{0,\nu}(0,1)$ is the completion of $C^\infty_c(0,1)$ under the norm $$\|f\|_\nu := \left(\int_0^1 f'(x)^2 + \left(\xi^2 x^2 + \frac{\nu^2-1/4}{x^2}\right)f(x)^2 dx\right)^{1/2}$$ **Role of Hardy Inequality**: The key inequality $$\int_0^1 \frac{u(x)^2}{x^2} dx \leq 4\int_0^1 u'(x)^2 dx, \quad \forall u \in H^1((0,1),\mathbb{R}), u(0)=0$$ ensures the positive definiteness and self-adjointness of the operator. ### Explicit Representation of Eigenfunctions (Proof of Theorem 1.1) **Step 1: Transformation to Kummer Equation** Using the ansatz $f(x) = e^{-\xi x^2/2}g(x)$ to eliminate the harmonic potential term, then setting $g(x) = x^\alpha h(x)$ with $\alpha = 1/2 + \nu$ (the larger positive root of $-\alpha(\alpha-1) + \nu^2 - 1/4 = 0$). Finally, through the variable substitution $z = \xi x^2$, the function $\tilde{h}(z) = h(x)$ satisfies the Kummer equation: $$z\tilde{h}''(z) + (1+\nu-z)\tilde{h}'(z) + \frac{E-2\xi(1+\nu)}{4\xi}\tilde{h}(z) = 0$$ **Step 2: Boundary Condition Analysis** The two fundamental solutions of the Kummer equation are: - $\varphi_1(x) = e^{-\xi x^2/2}x^{1/2+\nu}M(a, b, \xi x^2)$ (Kummer function) - $\varphi_2(x) = e^{-\xi x^2/2}x^{1/2+\nu}\bar{M}(a, b, \xi x^2)$ (second kind solution) where $a = -(E-2\xi(1+\nu))/(4\xi)$ and $b = 1+\nu$. Through analysis: - When $\nu > 0$, $\varphi_2 \notin H^1(0,1)$, so $B=0$ is required - When $\nu = 0$, although $\varphi_1, \varphi_2 \notin H^1(0,1)$, we have $\varphi_2 \notin H^1_{0,\nu}(0,1)$ - The boundary condition $\varphi_1(1) = 0$ gives $M(a, 1+\nu, \xi) = 0$ ### Proof Method for Eigenvalue Lower Bounds **Method 1: Using z-Zero Counting Function** (Section 3.1) Citing Lemma 3.1 from Olver's handbook: For $a < 0$, the number of positive z-zeros of $M(a,b,z)$ is $p(a,b) = \lceil -a \rceil$. By the Sturm-Liouville oscillation theorem, the $k$-th eigenfunction $g_{\xi,k}$ has $k+2$ zeros on $[0,1]$. Since $g_{\xi,k}(x) = 0 \Leftrightarrow M(a_{\xi,k}, 1+\nu, \xi x^2) = 0$ (except at $x=0$), $M(a_{\xi,k}, 1+\nu, \xi x^2)$ has at least $k+1$ positive x-zeros. Therefore $p(a_{\xi,k}, b) \geq k+1$, i.e., $-a_{\xi,k} > k$, which upon expansion gives $\lambda_{\xi,k}/\xi > 4k + 2(1+\nu)$ (strict inequality). **Method 2: Min-Max Principle** (Section 3.2) Through the variable substitution $y = \sqrt{\xi}x$, $v(y) = u(x)$, the operator becomes: $$\tilde{G}_\xi = -\partial_y^2 + y^2 + \frac{\nu^2-1/4}{y^2}, \quad \text{on } (0,\sqrt{\xi})$$ satisfying $\lambda_{\xi,k} = \xi\tilde{\lambda}_{\xi,k}$. Comparing with the operator $G$ on the full line $\mathbb{R}^+$ (with domain on $\mathbb{R}^+$), whose eigenvalues are $\mu_k = 4k + 2(1+\nu)$ (Proposition 3.2). Applying the Rayleigh formula: $$\tilde{\lambda}_{\xi,k} = \max_{f \in \tilde{E}_{n,k}} \int_0^{\sqrt{\xi}} |f'|^2 + x^2|f|^2 + \frac{\nu^2-1/4}{x^2}|f|^2 dx \geq \mu_k$$ ### WKB-Type Method for Upper Bounds on Low-Lying Eigenvalues (Proof of Theorem 1.4) **Core Idea**: Construct quasi-modes $\phi_{\xi,k}$ that approximately satisfy the eigenvalue equation, then apply Lemma 4.1: $$\text{dist}(\lambda, \sigma(A)) \leq \frac{\|(A-\lambda)u\|}{\|u\|}$$ **Quasi-Mode Construction**: $$\phi_{\xi,k}(x) = \Phi_k(x) - \frac{x^{1/2+\nu}}{\xi^{1/4+\nu/2}}\Phi_k(\sqrt{\xi})$$ where $\Phi_k(x) = e^{-x^2/2}x^{1/2+\nu}\frac{k!}{(1+\nu)_k}L^{(\nu)}_k(x^2)$ is an eigenfunction of the full-line operator $G$ ($L^{(\nu)}_k$ is the Laguerre polynomial). **Key Technical Estimates**: 1. **Lemma 4.2** (Exponential Decay of Boundary Values): For $k \leq \lfloor\tau\xi/4\rfloor$, $$|\Phi_k(\sqrt{\xi})| \leq C_1\frac{\Gamma(1+\nu)}{2}e^{-C_2\xi}$$ The proof uses Plancherel-Rotach type formulas (Lemma 4.4): For $r = (4n+2\nu+2)\cosh(\theta)^2$, $$e^{-r/2}L^{(\nu)}_n(r) = (-1)^n\frac{n^{\nu/2-1/4}}{2}\frac{e^{(n+(\nu+1)/2)(2\theta-\sinh(2\theta))}}{(\pi\sinh(\theta))^{1/2}r^{\nu/2+1/4}}(1+O(1/n))$$ 2. **Lemma 4.3** ($L^2$ Norm Lower Bound): For $k \leq \lfloor\xi/4\rfloor$, $$\int_0^{\sqrt{\xi}}\Phi_k(x)^2 dx \geq (1-\delta)\frac{\Gamma(\nu+1)}{2}\frac{k!}{(1+\nu)_k}$$ The proof uses Plancherel-Rotach formula analysis in the oscillatory region. **Error Estimate**: $$\|\tilde{G}_\xi\phi_{\xi,k} - \mu_k\phi_{\xi,k}\| \leq \frac{|\Phi_k(\sqrt{\xi})|}{\xi^{1/4+\nu/2}}\|x^{5/2+\nu}\| \leq C_1\frac{\Gamma(1+\nu)}{2}\frac{e^{-C_2\xi}}{\xi^{11/4+\nu/2}(6+2\nu)}$$ Combining with the norm lower bound yields: $$\text{dist}(\mu_k, \sigma(\tilde{G}_\xi)) \leq C\xi^{11/4+3\nu/2}e^{-C_2\xi}$$ Finally, the lower bound from Theorem 1.3 is used to verify that no eigenvalues are missed. ## Technical Innovation Points ### 1. Treatment Strategy for Singular Potentials By choosing an appropriate weight $x^{1/2+\nu}$ via the Frobenius method, the main contribution of the singular term $(\nu^2-1/4)/x^2$ is absorbed into the basis functions. This is a standard but crucial technique for handling inverse square potentials. ### 2. Bidirectional Conversion Between Spectral Analysis and Special Function Zeros An exact equivalence relationship is established between operator spectral theory problems and special function zero problems (equation (13)), allowing tools from both fields to be mutually borrowed. ### 3. Multi-Scale Analysis - **Low Energy Levels** ($k \leq \lfloor\tau\xi/4\rfloor$): Eigenvalues stabilize to linear behavior $4k + 2(1+\nu)$ - **Intermediate Region** ($\lfloor\tau\xi/4\rfloor < k \leq \lfloor\xi/(4\tau)\rfloor$): Should be connected with Airy equation - **High Energy Levels** ($k > \lfloor\xi/(4\tau)\rfloor$): Should be connected with Bessel equation ### 4. Refinement of Quasi-Mode Method By subtracting boundary values to construct quasi-modes satisfying boundary conditions, direct solution of singular boundary value problems is avoided. Application of the Plancherel-Rotach formula provides exponentially precise error control. ### 5. Obtaining Strict Inequalities Through the zero counting function method, strict inequality $\lambda_{\xi,k}/\xi > 4k + 2(1+\nu)$ is obtained, which is relatively rare in spectral theory (typically only non-strict inequalities can be obtained). ## Summary of Main Results ### Theorem 1.3 (Two-Sided Bounds for Eigenvalues) For all $\xi > 0$ and $k \geq 0$: $$\frac{\lambda_{\xi,k}}{\xi} > 4k + 2(1+\nu)$$ There exist $c \in (0,\pi^2)$ and $\xi_0 > 0$ such that for $\xi \geq \xi_0$, $k \geq 0$: $$\lambda_{\xi,k} \geq ck^2$$ ### Theorem 1.4 (Exponentially Precise Upper Bounds for Low-Lying Eigenvalues) For $\tau \in (0,1)$, there exist $\xi_\tau > 0$, $C_1, C_2 > 0$ such that for $\xi \geq \xi_\tau$, $k \leq \lfloor\tau\xi/4\rfloor$: $$\frac{\lambda_{\xi,k}}{\xi} \leq 4k + 2(1+\nu) + C_1e^{-C_2\xi}$$ ### Theorem 1.5 (Characterization of Kummer Function a-Zeros) For $b \geq 1$, $\xi > 0$, the solutions to $M(a,b,\xi) = 0$ form a strictly decreasing sequence: $$\cdots < a_{\xi,k} < a_{\xi,k-1} < \cdots < a_{\xi,0} < 0$$ satisfying: 1. $a_{\xi,k} < -k$ (for all $k \in \mathbb{N}$) 2. For $k \leq \lfloor\tau\xi/4\rfloor$: $-k - C_1e^{-C_2\xi} \leq a_{\xi,k} < -k$ ## Related Work ### Semiclassical Analysis - **Classical Literature**: [voros1981spectre, simon1983semiclassical, helffer1984puits] and others established semiclassical theory for regular potentials - **Distinction of This Work**: Handles singular inverse square potentials with lack of regularity at boundaries ### Special Function Zero Theory - **z-Zeros**: [ahmed1982properties, slater1960confluent, boussaada2022some] and others studied zeros of Kummer functions with respect to the argument - **a-Zeros**: [OlverHandbook2010, buchholz2013confluent] mention asymptotic behavior $a \sim -k^2\pi^2/(4z)$ but without proof or references - **Contribution of This Work**: First systematic characterization of non-asymptotic a-zeros with uniform estimates ### Spectral Theory of Inverse Square Potentials - **Hardy Inequality**: [cannarsa2008carleman, vazquez2000hardy] studied function spaces in the critical case $\nu=0$ - **Bessel Operators**: [martinez2018cost] studied spectral properties of pure inverse square potentials (without harmonic term) ### Laguerre Polynomials and Plancherel-Rotach Formula - **Classical Results**: [szeg1939orthogonal] provides orthogonality and asymptotic formulas for Laguerre polynomials - **Application in This Work**: These formulas are used to control quasi-mode errors ## Conclusions and Discussion ### Main Conclusions 1. In the semiclassical limit $\xi \to +\infty$, the low-lying eigenvalues ($k \leq O(\xi)$) of the Schrödinger operator with inverse square potential converge to $\xi(4k + 2(1+\nu))$ at exponential rate 2. An exact correspondence is established between eigenvalues and Kummer function a-zeros, providing the first systematic characterization of non-asymptotic a-zero locations 3. For different energy level ranges, eigenvalues exhibit different asymptotic behaviors: - Low energy levels: linear behavior $\sim 4\xi k$ - High energy levels: quadratic behavior $\sim ck^2$ ### Limitations and Future Directions **Incomplete Work Explicitly Noted by the Author**: 1. **Intermediate Energy Level Analysis** ($\lfloor\tau\xi/4\rfloor < k \leq \lfloor\xi/(4\tau)\rfloor$): - Should establish connection with Airy equation - WKB analysis becomes complex due to singular terms - Uniform asymptotic expansions from [erdelyi1957asymptotic] may be useful 2. **High Energy Level Analysis** ($k > \lfloor\xi/(4\tau)\rfloor$): - Should establish connection with Bessel equation - Bessel-type expansions from [erdelyi1957asymptotic] may be applicable - Rouché's theorem needed to verify completeness of zeros 3. **Generalization to More General Potentials**: - Replace $\xi^2 x^2$ with $\xi^2 q(x)^2$ where $q(x) \sim q'(0)^2 x$ (near $x=0$) - Utilize perturbation theory for self-adjoint operators **Technical Difficulties**: - Airy-type expansions cannot be extended to complex plane, preventing use of Rouché's theorem - Remainders may be too large as $\xi \to +\infty$ - Need to ensure no zeros are missed ### Applicability of Methods This paper's methods apply to: - Operators whose eigenfunctions can be expressed explicitly using special functions - Cases where semiclassical parameters appear in the potential function - Singular Sturm-Liouville problems on bounded intervals ## In-Depth Evaluation ### Strengths 1. **Theoretical Rigor**: - All theorems have complete proofs - Two independent proof methods provided (counting function method and min-max method) - Careful treatment of differences between critical case $\nu=0$ and non-critical case $\nu>0$ 2. **Precision of Results**: - Exponentially precise upper bounds $O(e^{-C\xi})$ are very strong results in semiclassical analysis - Strict inequalities obtained (first bound in Theorem 1.3) - Uniform estimates (uniform in both $\xi$ and $k$) 3. **Cross-Disciplinary Contribution**: - Simultaneously advances spectral theory and special function theory - Fills gap in Kummer function a-zero theory - Provides new technical tools for related fields 4. **Technical Innovation**: - Clever quasi-mode construction (satisfying boundary conditions by subtracting boundary values) - Refined application of Plancherel-Rotach formula - Clear multi-scale analysis framework ### Weaknesses 1. **Completeness of Results**: - Only covers low energy levels $k \leq O(\xi)$ - Analysis of intermediate and high energy levels missing - Authors acknowledge these as future work 2. **Non-Explicit Constants**: - Constants $C_1, C_2$, etc. are not given explicit expressions - May affect operability in practical applications 3. **Lack of Numerical Verification**: - No numerical experiments verify theoretical results - Cannot intuitively observe convergence rates 4. **Limitations of Generalization**: - Methods strongly depend on explicit representation of eigenfunctions - Not easily generalizable to arbitrary potentials ### Impact Assessment **Theoretical Impact**: - Provides new paradigm for semiclassical analysis of singular potentials - First systematic study of Kummer function a-zeros may inspire follow-up work - Method connecting spectral theory and special functions is instructive **Practical Value**: - Provides precise predictions for quantum mechanical problems with inverse square potentials - Potentially applicable to models in atomic and molecular physics - Provides theoretical foundation and error estimates for numerical computation **Reproducibility**: - Proofs are complete and technical route is clear - Primarily relies on known special function properties - Theoretical results are easy to verify (though numerical experiments are lacking) ### Applicable Scenarios 1. **Quantum Mechanics**: - Quantum systems with inverse square potentials (e.g., radial equation of hydrogen-like atoms) - Refinement of semiclassical approximations 2. **Special Function Theory**: - Problems requiring Kummer/Whittaker function zero information - Asymptotic analysis of confluent hypergeometric functions 3. **Numerical Analysis**: - Verification of numerical methods for eigenvalue problems - Error estimation and convergence analysis 4. **Mathematical Physics**: - Other problems transformable to Kummer equations - Singular Sturm-Liouville problems ## Supplementary Technical Details ### Role of Hardy Inequality Hardy inequality (equation (2)): $$\int_0^1 \frac{u(x)^2}{x^2}dx \leq 4\int_0^1 u'(x)^2 dx$$ ensures: 1. Positive definiteness of operator $G_\xi$ (when $\nu^2 - 1/4 \geq -1$) 2. Equivalence between norm $\|\cdot\|_\nu$ and standard $H^1_0$ norm (when $\nu > 0$) 3. Specialness of critical case $\nu = 0$: $H^1_0(0,1) \subsetneq H^1_{0,\nu}(0,1)$ ### Application of Sturm-Liouville Theory - Simplicity of eigenvalues (spectrum is simple) - Oscillatory properties of eigenfunctions ($k$-th eigenfunction has $k+2$ zeros) - These properties are crucial for the zero counting method ### Elegance of Plancherel-Rotach Formula Lemma 4.4 provides uniform asymptotic expansions of Laguerre polynomials in different regions: - **Oscillatory Region** ($r = O(n)$): Trigonometric type - **Transition Region** ($r \approx 4n$): Requires Airy functions - **Exponential Decay Region** ($r > 4n$): Exponential type The paper cleverly exploits estimates in the exponential decay region to control $\Phi_k(\sqrt{\xi})$. ## Key References 1. **Special Functions**: - [OlverHandbook2010] NIST Handbook of Mathematical Functions - [szeg1939orthogonal] Szegő, Orthogonal Polynomials - [buchholz2013confluent] Buchholz, The Confluent Hypergeometric Function 2. **Semiclassical Analysis**: - [helffer1984puits] Helffer & Sjöstrand, Multiple wells in the semiclassical limit - [simon1983semiclassical] Simon, Semiclassical analysis of low lying eigenvalues 3. **Hardy Inequality**: - [vazquez2000hardy] Vázquez & Zuazua, The Hardy inequality and the asymptotic behaviour 4. **Related Spectral Analysis**: - [allonsius2021analysis] Allonsius et al., Analysis of the spectrum of Grushin operators --- **Overall Assessment**: This is a high-quality mathematical analysis paper that is technically rigorous with novel and precise results. Although only low-energy level analysis is completed, it establishes a solid foundation for complete theory. The paper makes substantive contributions at the intersection of spectral theory and special function theory, and is expected to inspire subsequent research. The main limitations are the incompleteness of results (intermediate and high energy levels not covered) and lack of numerical verification, but these do not diminish the value of the completed portions.