[论文解读] Accurate calculation of eigenvalues and eigenfunctions. I: Symmetric potentials

Francisco M. Fernández|ArXiv.org|Jul 3, 2008

Quantum, superfluid, helium dynamics参考文献 18被引用 6

一句话总结

本文提出Riccati-Padé方法（RPM），用于精确计算具有对称势的单维量子系统（如非谐振子和双阱振子）的本征值与本征函数。通过将波函数的对数导数展开为泰勒级数，并求解相应的递推关系，该方法可快速收敛得到能量本征值与期望值的上下界，对非谐振子和双阱势系统的结果精度可达20位以上。

ABSTRACT

We describe a method for the calculation of accurate energy eigenvalues and expectation values of observables of separable quantum-mechanical models. We discuss the application of the approach to one-dimensional anharmonic oscillators with symmetric potential-energy functions.

研究动机与目标

• 开发一种系统化方法，用于计算可分离的量子力学系统中对称势的精确能量本征值与期望值。
• 提供一种替代标准微扰法或矩阵法的方法，避免计算矩阵元的需要。
• 通过符号计算展示该方法在非谐振子与双阱振子中的高收敛速度与可靠性。
• 强调该方法在使用计算机代数系统（CAS）时对研究生或高年级本科生的教育价值。

提出的方法

• 该方法基于从薛定谔方程导出的Riccati方程，以对数导数 $ f(x) = \frac{s}{x} - \frac{\psi'(x)}{\psi(x)} $ 作为核心变量。
• 假设 $ f(x) $ 在 $ x=0 $ 处进行泰勒级数展开，令 $ z = x^2 $，从而导出系数 $ f_j(E) $ 的递推关系。
• 递推系数 $ f_j(E) $ 是能量 $ E $ 的多项式，可由此构造Padé逼近以定位对应于本征值的根。
• 本征值通过由 $ f_j $ 系数构成的Hankel行列式求根得到，从而获得上下界 $ E_n^{[D,0]} < E_n < E_n^{[D,1]} $。
• 该方法应用于对称势如 $ V(x) = x^4 + \beta x^2 $，并研究了不同 $ \beta $ 值下的收敛性。
• 期望值如 $ \langle x^2 \rangle $ 可从级数系数中提取，例如，基态时 $ \langle x^2 \rangle = e^{(1)} $。

实验结果

研究问题

• RQ1 Riccati-Padé方法是否能在不计算矩阵元的情况下，为对称量子势提供高度精确的本征值？
• RQ2 与传统的矩量法相比，RPM在非谐振子中的收敛速度如何？
• RQ3 当双阱势的势垒深度变化时，该方法的行为表现如何？
• RQ4 为何在RPM中会出现虚假根与多重根？如何识别并过滤这些根？
• RQ5 该方法在多大程度上可直接从级数展开中计算期望值？

主要发现

• 对于势 $ V(x) = x^4 + \beta x^2 $ 的非谐振子，当 $ \beta = 0 $ 时，RPM在本征值上实现了20位精度，得到 $ \langle x^2 \rangle = 0.3620226487886768452 $。
• 当 $ \beta = -15 $ 时，基态能量被严格界定为 $ -50.8413872843819547 < E_0 < -50.8413872843819543 $，第一激发态也具有类似精度。
• 该方法可提供能量本征值的严格上下界，满足 $ E_n^{[D,0]} < E_n < E_n^{[D,1]} $，且随着阶数 $ D $ 增大，收敛性显著提升。
• RPM的收敛速度显著快于矩量法，尤其在深双阱势情况下更为明显。
• 在根查找过程中会出现虚假根与多重根，但正确的本征值始终可作为物理解释的根被稳定识别。
• 期望值如 $ \langle x^2 \rangle $ 可直接从级数系数中获得，基态时 $ \langle x^2 \rangle = e^{(1)} $。