Laguerre-Sobolev orthogonal Polynomials and Boundary Value Problems on a semi-infinite domain

技术摘要：Laguerre–Sobolev 正交多项式与半无限域上的边值问题

问题陈述
本文研究了半无限区间 $(0, +\infty)$ 上二阶边值问题（BVP）的数值解，特别关注具有奇异势能（如反距离势 $V(x) = 1/x$）的定常薛定谔型方程。虽然基于 Laguerre 多项式的经典谱方法在处理半无限域上的光滑解时已非常成熟，但在处理无穷远处的边界条件或微分算子中的奇异系数时，效率可能会降低。作者试图开发一种能够将导数信息和边界行为直接纳入近似空间的框架，以更有效地处理这些奇异性。

方法论
作者引入了一族 Laguerre–Sobolev 正交多项式，记作 $\{S_n(x)\}_{n \geq 0}$，该多项式与源自目标 BVP 变分形式的特定 Sobolev 内积相关联。内积定义为：
$$\langle u, v \rangle_\lambda = \lambda \int_0^{+\infty} \frac{u(x)v(x)}{x} dx + \int_0^{+\infty} u'(x)v'(x) dx$$
测试函数选自空间 $\mathcal{P} \times x e^{-x/2}$，确保其在两个端点（$0$和$+\infty$）处均消失。

该方法的核心包括：

1. 构造：通过对序列 $\{x^k x e^{-x/2}\}$ 关于上述 Sobolev 内积应用 Gram-Schmidt 正交化过程，生成正交基 $\{S_n(x)\}$。
2. 解析推导：建立新的 Sobolev 多项式 $S_n(x)$ 与经典 Laguerre 多项式 $L_n^{(1)}(x)$ 之间的显式联系公式。这包括推导连接系数的递推关系，并利用精细化的 Laguerre 多项式比值渐近性来分析其渐近行为。
3. 生成函数：推导出涉及 Bessel 函数的 Sobolev 多项式的闭式生成函数。
4. 谱算法：开发一种完全对角化的谱方法。通过利用连接公式，可以递归地计算解的 Fourier–Sobolev 系数。这种方法避免了求解线性方程组，因为所得算法仅涉及对经典 Laguerre 多项式积分的求值。

主要贡献

• 新型正交族：定义并分析了一类专门用于处理 $(0, +\infty)$ 上具有奇异势能的 BVP 的 Laguerre–Sobolev 正交多项式。
• 连接公式：推导了将 $S_n(x)$ 与 $L_n^{(1)}(x)$ 联系起来的显式关系，包括连接系数 $a_n$ 的递推关系及其渐近展开式（$a_n \sim 1 - \sqrt{4\lambda/n}$）。
• 生成函数：提供了用第一类 Bessel 函数表示的 Sobolev 多项式的生成函数。
• 对角化算法：制定了一种递归谱方法，无需求解线性方程组即可计算解的展开系数，从而在降低计算复杂度的同时，保持了针对奇异问题的稳定性。

结果
论文通过数值实验求解了一个具有反距离势（$V(x)=1/x$）和 Dirichlet 边界条件的非齐次薛定谔方程。

• 指数收敛性：对于表现出指数衰减特性的解（$u(x) = x \cos x e^{-x}$），该方法展示了谱精度和指数收敛性。随着项数的增加，Sobolev 范数下的误差迅速下降。
• 局限性：作者指出，对于不具备指数衰减特性（例如有理衰减）的解，近似函数的表现较差，这符合谱方法在半无限域上的预期表现。
• 效率：算法的递归性质使得计算成本与标准 Laguerre 方法相当，同时为处理此类特定奇异问题提供了更高的稳定性。

意义
本文声称，所提出的 Laguerre–Sobolev 谱方法为处理涉及奇异系数和无穷远边界条件的类问题提供了一种极具价值的替代方案。通过将导数信息和边界行为直接嵌入近似空间，该方法实现了完全对角化的系统。这消除了求解线性方程组的需求，显著降低了计算复杂度。该工作提供了严谨的解析基础（连接公式、渐近性、生成函数），为这些方法在处理具有奇异势能的薛定谔型问题中的实际应用提供了支持。