2025-11-23T00:13:15.391562

On the eigenvalues of the spheroidal wave equation

Schmid

This paper presents some new results on the eigenvalues of the spheroidal wave equation. We study the angular and Coulomb spheroidal wave equation as a special case of a more general linear Hamiltonian system depending on three parameters. We prove that the eigenvalues of this system satisfy a first-order quasilinear partial differential equation with respect to the parameters. This relation offers a new insight on how the eigenvalues of the spheroidal wave equation depend on the spheroidal parameter. Apart from analytical considerations, the PDE we obtain can also be used for a numerical computation of spheroidal eigenvalues.

academic

On the eigenvalues of the spheroidal wave equation

基本信息

论文ID: 2402.07133
标题: On the eigenvalues of the spheroidal wave equation
作者: Harald Schmid (University of Applied Sciences Amberg-Weiden, Germany)
分类: math.AP math-ph math.CA math.MP
发表时间: 2024年2月11日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2402.07133

摘要

本文提出了关于球面波动方程特征值的新结果。作者将角球面波动方程和库仑球面波动方程作为依赖三个参数的更一般线性哈密顿系统的特例进行研究。证明了该系统的特征值满足关于参数的一阶拟线性偏微分方程。这一关系为球面波动方程的特征值如何依赖于球面参数提供了新的洞察。除了分析考虑外，所得到的偏微分方程还可用于球面特征值的数值计算。

研究背景与动机

核心问题: 球面波动方程在物理学和工程学的许多领域中都有出现，如量子力学、电磁学、信号处理等。计算其特征值是一个具有挑战性的任务，需要新的理论洞察和数值方法。
重要性:
- 角球面波动方程来源于椭球坐标系中亥姆霍兹方程的分离
- 库仑球面波动方程在天体物理学和分子物理学中提供精确的波函数，例如固定核的单电子双原子分子
- 这些方程的特征值计算对相关物理问题的求解至关重要
现有方法局限性:
- 传统的数值方法在参数值较大时会遇到舍入误差或数字消除等数值问题
- 缺乏描述特征值与参数之间解析关系的统一理论框架
- 现有计算方法的精度和效率有待提高
研究动机:
- 建立特征值与参数间的解析关系
- 提供新的数值计算途径
- 深入理解球面波动方程的数学结构

核心贡献

理论贡献: 将库仑球面波动方程转化为2×2线性哈密顿系统，建立了与原方程特征值问题的等价关系
主要定理: 证明了哈密顿系统的特征值Λ满足一阶拟线性偏微分方程： $2u_1\frac{\partial\Λ}{\partial u_1} + \left((Λ+2)(u_1+u_2^2)+Λ+u_2+u_3\right)\frac{\partial Λ}{\partial u_2} + \cdots = (1+2μ-2Λu_2)(Λ+2)-2μ$
数值方法: 提出了基于特征线方法求解偏微分方程来计算球面特征值的新途径
特殊情况简化: 对椭球面波动方程(β=0, γ²>0)的情况，进一步简化为强迫无粘Burgers方程的形式

方法详解

任务定义

研究库仑球面波动方程(CSWE)的特征值问题： $\frac{d}{dx}\left((1-x^2)\frac{d}{dx}w(x)\right) + \left(λ+βx+γ^2(1-x^2)-\frac{μ^2}{1-x^2}\right)w(x) = 0$

其中需要找到使方程在(-1,1)上有非平凡有界解的λ值（特征值）。

模型架构

1. 参数变换 引入技术参数：

$u_1 = γ^2$ , $u_2$ , $u_3$ 与 $β$ , $γ^2$ , $λ$ 的关系： $γ^2 = u_1, \quad β = -u_3-2(μ+1)u_2, \quad λ = u_3+μ(μ+1)$
$α := \frac{1}{2}(μ+1)$

2. 哈密顿系统构造 将CSWE转化为2×2微分系统： $\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}y'(z) - \begin{pmatrix} \frac{2α}{z}+\frac{α}{1-z}+2u_2 & \frac{u_3}{z}+2(u_1+u_2^2) \\ \frac{2α}{z}+\frac{α}{1-z}+2u_2 & \frac{2α}{z} \end{pmatrix}y(z) = Λ\begin{pmatrix} \frac{1}{1-z} & 0 \\ 0 & \frac{1}{z} \end{pmatrix}y(z)$

3. 自伴算子理论

证明微分表达式在端点处为极限点情况
建立 $L^2_W((0,1),\mathbb{C}^2)$ 希尔伯特空间中的自伴算子 $T(u)$
特征值为实数且孤立

技术创新点

1. 变形方法(Deformation Method)

采用文献4中的变形方法建立特征值满足的偏微分方程
构造特殊的矩阵函数 $G(z,Λ,u)$ 满足变形方程
不需要知道特征函数的显式形式

2. 连接系数方法

利用基本解在奇点处的渐近行为
通过连接系数 $Θ(Λ,u)$ 的零点确定特征值
提供递推公式(16)进行数值计算

3. 特征线方法

将偏微分方程转化为常微分方程组
沿特征曲线追踪参数变化
实现从已知特征值计算新特征值

实验设置

数值实现细节

使用引理3中的极限公式计算连接系数 $Θ(Λ,u)$
递推关系(16)用于级数系数计算，精度为 $O(k^{δ-μ-2})$
采用割线法求解 $Θ(Λ,u)=0$
使用Runge-Kutta方法求解特征线方程组

测试参数

$μ = 1$ , $a = 5$ (初始参数)
内部计算采用高精度，结果显示16位小数

实验结果

主要结果

数值算例（附录中的具体例子）：

初始条件： $μ = 1$ , $u_1(0) = 5$ , $u_2(0) = 0$
找到 $Θ̃(Λ)$ 的零点： $Λ_0 = -0.8417200168449013$
特征曲线在 $t_0 = 0.2793371978706399$ 处 $Λ(t_0) = 0$
得到球面波动方程特征值： $λ = -5.2736106330552739$
对应参数： $γ^2 = 8.7417666942941543$

理论发现

1. 等价性定理（引理1）： $Λ = 0$ 是哈密顿系统 $T(u_1,u_2,u_3)$ 的特征值当且仅当 $λ = u_3+μ(μ+1)$ 是相应CSWE的特征值。

2. 偏微分方程（定理2）：特征值 $Λ(u_1,u_2,u_3)$ 满足复杂的一阶拟线性PDE(5)。

3. 简化情况（定理5）：对椭球面情况，可简化为强迫无粘Burgers方程： $\frac{\partial ω}{\partial t} + 2ω\frac{\partial ω}{\partial u} = 2(μ+1)^2\frac{\tanh u}{\cosh^2 u} + e^{2t}\sinh(2u) + e^t\cosh u$