[论文解读] Numerov numerical method applied to the Schrödinger equation

Francisco Caruso, V. Oguri|arXiv (Cornell University)|Mar 17, 2014

Quantum Mechanics and Non-Hermitian Physics参考文献 6被引用 3

一句话总结

本文将数值求解方法（Numerov方法）应用于量子力学中的本征值问题，具体针对一维谐振子和氢原子系统。该研究将标准Numerov方法扩展至处理薛定谔方程中的的一阶导数项，通过与解析解对比验证了数值结果的高精度，并采用二阶有限差分格式成功计算了能级和径向波函数。

ABSTRACT

In this paper it is shown how to solve numerically eigenvalue problems associated to second order linear ordinary differential equations, containing also terms which depend on the variable. A didactic presentation of the Numerov Method is given and, in the sequel, it is applied to two quantum non-relativistic problems with well known analytical solutions: the simple harmonic oscillator and the hydrogen atom. The numerical results are compared to those obtained analytically.

研究动机与目标

• 为在量子力学中求解二阶线性常微分方程提供一种教学性介绍，重点介绍Numerov方法。
• 将标准Numerov方法扩展，以处理波函数一阶导数相关的项。
• 对两个基准量子系统（谐振子和氢原子）的定态薛定谔方程进行数值求解。
• 通过与解析解对比，验证数值本征值和波函数的准确性。
• 展示该方法在现代量子力学计算中的鲁棒性与通用适用性。

提出的方法

• 采用Numerov方法作为二阶有限差分格式对薛定谔方程进行离散化，避免显式使用一阶导数项。
• 推导出一种改进的Numerov公式，以包含依赖于波函数一阶导数的项，从而实现对具有角动量的氢原子问题的应用。
• 利用三点模板从前向后传播波函数，从初始条件出发，确保方法的稳定性和精度。
• 在空间区域的两端（如原点和无穷远）施加边界条件，以求解本征值问题。
• 使用ROOT框架在C++中实现该算法，代码可通过公开链接获取。
• 通过将数值计算的能级和径向波函数与精确解析结果对比，对方法进行验证。

实验结果

研究问题

• RQ1 当将Numerov方法扩展以包含一阶导数项时，其在计算一维谐振子能级本征值方面的精度如何？
• RQ2 改进后的Numerov方法能否正确再现三维氢原子的径向波函数和能级？
• RQ3 与解析解相比，计算得到的本征值在数值误差方面有多大？
• RQ4 该方法在求解涉及非零角动量的薛定谔方程的奇异边值问题时表现如何？
• RQ5 Numerov方法在多大程度上可被推广并应用于超越标准教科书案例的现代量子力学问题？

主要发现

• 一维谐振子和氢原子（l=1）的数值能级与解析解高度一致，偏差小于2%。
• 对于氢原子（l=1），计算得到的能级为：E2 = -3.47 eV（解析值：-3.40 eV），E3 = -1.54 eV（解析值：-1.51 eV），E4 = -0.83 eV（解析值：-0.85 eV）。
• 数值计算得到的径向波函数在形状和节点结构上与解析解完全一致，证实了方法实现的正确性。
• 通过引入有效势能和径向导数项，改进的Numerov方法成功处理了氢原子薛定谔方程中的一阶导数项。
• 该方法在大空间区域上表现出二阶收敛性和稳定性，包括原点处的奇异性行为。
• 使用ROOT框架在C++中实现的代码运行高效，适用于量子力学领域的教学与科研应用。