[论文解读] Estimating Spectral Density Functions for Sturm-Liouville problems with two singular endpoints

Charles T. Fulton, D. B. Pearson+1 more|arXiv (Cornell University)|Mar 12, 2013

Spectral Theory in Mathematical Physics参考文献 18被引用 5

一句话总结

本文提出了一种新颖的数值方法，用于计算具有两个奇异端点的Sturm-Liouville问题的谱密度函数，扩展了以往针对正则左端问题的算法。通过结合在 x=0 处的正则奇异点附近的Frobenius幂级数、用于数值解的分段三角双曲样条，以及通过在无穷远处求解 (P,Q,R) 的初值问题来建立Titchmarsh-Weyl m-函数的两点边值公式，该方法在径向氢原子和贝塞尔方程中实现了接近机器精度的准确性，其精度和速度显著优于SLEDGE软件。

ABSTRACT

In this paper we consider the Sturm-Liouville equation -y"+qy = lambda*y on the half line (0,infinity) under the assumptions that x=0 is a regular singular point and nonoscillatory for all real lambda, and that either (i) q is L_1 near x=infinity, or (ii) q' is L_1 near infinity with q(x) --> 0 as x --> infinity, so that there is absolutely continuous spectrum in (0,infinity). Characterizations of the spectral density function for this doubly singular problem, similar to those obtained in [12] and [13] (when the left endpoint is regular) are established; corresponding approximants from the two algorithms in [12] and [13] are then utilized, along with the Frobenius recurrence relations and piecewise trigonometric - hyperbolic splines, to generate numerical approximations to the spectral density function associated with the doubly singular problem on (0,infinity). In the case of the radial part of the separated hydrogen atom problem, the new algorithms are capable of achieving near machine precision accuracy over the range of lambda from 0.1 to 10000, accuracies which could not be achieved using the SLEDGE software package.

研究动机与目标

• 将谱密度函数 f(λ) = ρ′(λ) 的表征扩展至具有两个奇异端点的Sturm-Liouville问题，其中 x=0 为正则奇异点，x=∞ 为非振荡极限点。
• 建立 f(λ) 的表示形式，即作为无穷远处一阶ODE系统 (P,Q,R)ᵀ 的解的有理函数，且满足特定渐近条件。
• 通过在有限匹配点 x₀(λ) 处进行高精度初值求解，将 [12,13] 中针对正则左端问题的数值算法适配至双重奇异情形。
• 开发并验证一种新型数值代码（AutoB），可在大 λ 范围内实现高精度与高效率的谱密度计算，涵盖径向氢原子和贝塞尔方程。
• 证明新方法克服了SLEDGE软件的局限性，特别是在大 λ 值下的精度和计算时间方面。

提出的方法

• 谱密度 f(λ) 表示为 f(λ) = 1 / [π(Pϕ² + Qϕϕ′ + Rϕ′²)]，其中 (P,Q,R)ᵀ 求解一阶线性ODE系统 dU/dx = A(x,λ)U，矩阵 A 依赖于 λ 和 q(x)，并在 x=∞ 处设定初值条件。
• 利用渐近条件 limₓ→∞ (P,Q,R)ᵀ = (√λ, 0, 1/√λ)ᵀ 作为边界条件，为每个 λ 定义唯一解。
• 在 x=0 附近，通过精确的Frobenius幂级数展开计算解 ϕ(x,λ) 和 ϕ′(x,λ)，确保在有限点 x₀(λ) > 0 处初始条件的高精度。
• 采用基于分段三角双曲样条的射击法，从 x₀(λ) 积分至匹配点 x，同时求解Sturm-Liouville方程和 (P,Q,R)ᵀ 系统，实现 f(λ) 的高精度计算。
• 该方法利用Titchmarsh-Kodaira公式将谱密度与 m-函数关联，使已知的解析性质可被纳入数值方案中。
• 所提出的算法通过 AutoB 代码实现，并在径向氢原子和贝塞尔方程上进行测试，与 SLEDGE 进行对比以验证其准确性和性能。

实验结果

研究问题

• RQ1 能否通过类似正则左端情形的Titchmarsh-Weyl m-函数两点边值公式，表征具有两个奇异端点的Sturm-Liouville问题的谱密度函数？
• RQ2 如何将专为正则左端问题设计的数值算法，适配至在 x=0 处具有正则奇异点、在 x=∞ 处为极限点的双重奇异情形？
• RQ3 结合Frobenius级数、样条积分与渐近边界条件的新型数值方法，在 f(λ) 计算中可达到的精度与计算效率如何？
• RQ4 该新方法是否能在大 λ 值下，特别是在径向氢原子等量子力学问题中，超越 SLEDGE 软件在精度和速度方面的表现？
• RQ5 该方法在长 λ 范围内，对谱密度函数的计算能达到多大程度的接近机器精度？

主要发现

• AutoB 代码在径向氢原子的谱密度函数计算中实现了接近机器精度的准确性，λ ∈ [0.1, 10000] 范围内相对误差低于 10⁻⁸，显著优于 SLEDGE。
• 对于 ν=1 的贝塞尔方程，AutoB 在 τ=10⁻⁴ 时达到 1.30×10⁻⁷ 的相对误差，在 τ=10⁻¹⁰ 时达到 1.73×10⁻¹³ 的相对误差，τ=10⁻¹⁰ 时计算时间低于 3 秒。
• SLEDGE 在氢原子问题中 λ > 1000 时无法达到 10⁻³ 的精度，且在 λ=1000 时耗时超过 177 秒，而 AutoB 在不到 4 秒内以 10⁻⁸ 精度完成相同计算。
• 该方法在多种势能场中表现出鲁棒性，包括 ℓ=1 的氢原子和贝塞尔方程，在所有测试的 λ 范围和容差水平下均保持一致的高精度。
• 利用 Frobenius 级数处理 x₀(λ) 处的初始条件，结合分段三角双曲样条，实现了高精度积分，使该方法适用于两端均具有强奇异性的复杂问题。
• 径向氢原子在 λ=10000 时的谱密度计算值达到 144274264.3，AutoB 与精确值在小数点后 8 位内一致，而 SLEDGE 的相对误差超过 10⁻³。