Exact WKB method for radial Schrödinger equation

Basic Information

• Paper ID: 2510.11766
• Title: Exact WKB method for radial Schrödinger equation
• Authors: Okuto Morikawa (RIKEN iTHEMS), Shoya Ogawa (Kyushu University)
• Classification: quant-ph nucl-th
• Publication Date: October 13, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.11766

Abstract

This paper revisits the exact WKB quantization of the radial Schrödinger problem from the perspective of modern resurgence theory, with emphasis on how to select and interpret "physically meaningful" quantization paths. By employing connection formulas at simple turning points and regular singularities, the authors demonstrate that non-trivial monodromy data yields the energy spectrum. Specifically, for the three-dimensional harmonic oscillator and three-dimensional Coulomb potential, the authors explicitly compute closed contours starting from +∞, indenting into the r<0 region around the origin, and returning to +∞. They propose that appropriate slicing of the closed path at r=0 provides a physical local basis for open paths from the origin to infinity. This analysis clarifies how mathematical monodromy data and physical boundary conditions fit together in the radial setting, thereby resolving disputes over path selection in resurgence-theory-based quantization.

Research Background and Motivation

1. Core Problem

The core problem addressed in this paper is how to select "physically meaningful" quantization paths in the exact WKB method, particularly for radial Schrödinger equations with regular singularities.

2. Problem Significance

• Theoretical Unification: Resurgence theory provides a unified framework for handling Stokes analysis, monodromy, and Voros data, yet conceptual tension exists between bound state problems (traditionally quantized via closed loops) and scattering/quasinormal mode problems (via boundary-to-boundary conditions)
• Contemporary Applications: Recent advances in quasinormal mode (QNM) analysis suggest that spectral conditions can be formulated as open connection problems along physically meaningful real axes, contradicting traditional closed-loop quantization
• Existing Method Limitations:
  • Lack of clear physical interpretation of the regular singularity at r=0 in radial problems
  • Absence of explicit proof of equivalence between different path choices (closed vs. open paths)
  • Systematic analysis lacking for how angular momentum phase (Maslov contribution l+1/2) enters the quantization condition

3. Research Motivation

By reinterpreting the exact WKB method through modern resurgence theory, establish equivalence between seemingly different path choices and provide a unified theoretical framework for radial problems.

Core Contributions

1. Explicit Proof of Path Equivalence: Demonstrates that given non-trivial monodromy data, the choice of basepoint and detailed path shape are gauge choices, showing equivalence between matrix (open path) conditions and closed-loop quantization conditions
2. Construction of Local Basis at Origin: Proposes a simplified method for constructing a local basis at r=0 through path equivalence, where the basis at |r|=ϵ connects via another transverse Borel direction from subleading wavefunctions at infinity
3. Explicit Calculations for Concrete Systems: Performs explicit computations for the three-dimensional harmonic oscillator and three-dimensional Coulomb potential, proving that period integrals of higher-order WKB corrections vanish, reducing quantization to leading action plus discrete Maslov phase
4. Unification of Physical Boundary Conditions and Mathematical Monodromy: Clarifies how mathematical monodromy data and physical boundary conditions fit together in the radial setting

Detailed Methodology

Problem Setup

Study the radial part of the three-dimensional Schrödinger equation:

[-ℏ²/2 d²/dr² + V(r) + ℏ²l(l+1)/(2r²) - E]φₗ(r) = 0

where V(r) is the potential, l is the angular quantum number, and the goal is to determine energy eigenvalues E.

WKB Ansatz and Riccati Equation

Introduce WKB ansatz:

φₗ(r,ℏ) = exp(∫ʳ dr' S(r',ℏ))
S(r,ℏ) = Σ(i=-1 to ∞) ℏⁱSᵢ(r)

Substituting into the radial Schrödinger equation yields the Riccati equation:

S² + dS/dr = ℏ⁻²Q(r)

where Q(r) = Q₀(r)/r + ℏ²Q₂(r)/r², Q₀(r) = 2rV(r)-E, Q₂(r) = l(l+1).

Stokes Analysis and Connection Formulas

Turning Point Connection Matrices:

M•₊ = (1  i)    M•₋ = (1  0)
      (0  1)          (i  1)

Regular Singularity Connection Matrices:

M◦₊ = (1  2i cos(π√(1+4Q₂)))
      (0  1              )

M◦₋ = (1                        0)
      (2i cos(π√(1+4Q₂))  1)

where √(1+4Q₂) = √(1+4l(l+1)) = 2l+1.

Quantization Condition

Through path analysis, obtain the unified quantization condition:

cosh(1/ℏ ∮ Sₒdd) + cos(π(2l+1)) = 0

where the second term encodes the angular momentum phase provided by the regular singularity at r=0.

Experimental Setup

Objects of Study

Two classical integrable systems are selected for detailed analysis:

1. Three-dimensional Harmonic Oscillator: V(r) = r²/2
2. Three-dimensional Coulomb Potential: V(r) = -e²/r

Analysis Methods

• Construct Stokes diagrams and determine turning points
• Calculate closed contour integrals
• Compare quantization conditions for open and closed paths
• Verify vanishing of higher-order WKB corrections

Technical Tools

• AKT (Aoki-Kawai-Takei) connection formulas
• Borel summation and transverse Borel directions
• Monodromy analysis and Voros period calculations

Experimental Results

Three-dimensional Harmonic Oscillator

Stokes Diagram Structure:

• Turning points: a₁ = -√(2E), a₂ = √(2E)
• Quantization condition obtained via closed path around origin

Energy Spectrum Result:

E = ℏ(2nᵣ + l + 3/2)

where nᵣ is the radial quantum number.

Key Findings:

• Contour integrals of higher-order WKB corrections Sᵢ(i≥0) vanish
• Langer correction l(l+1) ≃ (l+1/2)² emerges naturally
• Angular momentum phase l+1/2 originates from monodromy at the origin

Three-dimensional Coulomb Potential

Stokes Diagram Structure:

• Turning point: a = -e²/E > 0 (for ReE < 0)
• Origin r=0 is a regular singularity

Energy Spectrum Result:

E = -e⁴/(2ℏ²n²), n = nᵣ + l + 1

Maslov Index Analysis:

• Turning point contribution: 1/2
• Angular momentum phase contribution from origin: l+1/2
• Total phase: nᵣ + 1/2 + (l+1/2) = n

Path Equivalence Verification

For both systems, verified:

1. Closed path quantization conditions
2. Open path connection matrix conditions
3. Identical energy spectra from both methods

Related Work

Development of Exact WKB Theory

• Voros's complex WKB method as foundation
• Rigorous mathematical framework by AKT school
• Applications of modern resurgence theory

Quasinormal Mode Research

• Applications in black hole perturbation theory
• Success of open path quantization
• Use of Siegert boundary conditions

Path Integrals in Quantum Mechanics

• Modern understanding of Bohr-Sommerfeld quantization
• Maslov index theory
• Shape invariance and integrable systems

Conclusions and Discussion

Main Conclusions

1. Path Choice is a Gauge Choice: Once non-trivial monodromy data is specified, the energy spectrum is independent of basepoint choice and detailed path shape
2. Unified Understanding of Origin Phase: The angular momentum phase l+1/2 can be encoded either in the momentum (Langer correction) or in the boundary phase on the right-hand side of the quantization condition
3. Equivalence of Open and Closed Paths: Under correct Stokes data, the boundary condition pair (regular at origin, subleading at infinity) is equivalent to the closed-loop condition

Limitations

1. Restriction to Integrable Systems: Analysis focuses on integrable systems where higher-order WKB corrections vanish
2. Simple Turning Point Cases: Primarily addresses systems with simple single or double turning points
3. Absence of Numerical Verification: Lacks numerical verification for non-integrable systems

Future Directions

1. Multi-Turning-Point Systems: Extension to more complex potential functions
2. Confluent Limits: Study cases where turning points and poles merge
3. Scattering Problem Applications: Apply methods to scattering and QNM settings
4. Higher-Order Corrections: Quantify higher-order Voros periods in non-integrable cases

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Based on mature AKT theoretical framework with rigorous mathematical derivations
2. Conceptual Clarity: Successfully resolves conceptual confusion over path selection, providing unified understanding
3. Computational Completeness: Complete explicit calculations for classical systems
4. Physical Insight: Clarifies connections between mathematical formalism and physical intuition

Weaknesses

1. Limited Scope: Primarily restricted to integrable systems; applicability to general systems requires further verification
2. Numerical Aspects: Lacks numerical calculations to verify theoretical predictions
3. Experimental Connection: Insufficient direct connection to practical physical experiments

Impact

1. Theoretical Contribution: Provides new theoretical perspective for exact WKB methods
2. Methodological Value: Proof of path equivalence offers guidance for analyzing other quantum systems
3. Cross-Disciplinary Applications: Methods potentially applicable to similar problems in condensed matter physics and gravitational theory

Applicable Scenarios

1. Quantum Mechanical Bound States: Quantum systems with radial symmetry
2. Quasinormal Mode Analysis: Black hole physics and open quantum systems
3. Mathematical Physics: Differential equation problems requiring exact WKB analysis

References

The paper cites 15 important references, covering:

• Foundational resurgence theory literature 1-3
• Recent exact WKB advances 4-8
• Quasinormal mode related work 9-15
• Classical Siegert boundary condition literature 16

Summary: This paper successfully unifies different quantization methods for the radial Schrödinger equation through modern resurgence theory, providing important conceptual clarification and technical advances for exact WKB theory. While primarily focused on integrable systems, its theoretical framework and methodology hold significant value for broader quantum system analysis.