[论文解读] Power series Schrodinger eigenfunctions for a particle on the torus

M. Encinosa, Babak Etemadi|ArXiv.org|Apr 30, 2002

Crystallography and Radiation Phenomena参考文献 3被引用 3

一句话总结

本文提出一种基于角坐标 θ 的傅里叶展开的幂级数方法，用于计算在无闭式解的环面（torus）上运动的量子粒子的本征值与本征态。该方法通过傅里叶系数的递推关系求解薛定谔方程，实现了对 m=0 与 m≠0 态能量本征值与波函数的高精度数值计算，并通过与龙格-库塔法解的对比验证了收敛性。

ABSTRACT

The eigenvalues and a series representation of the eigenfunctions of the Schrodinger equation for a particle on the surface of a torus are derived.

研究动机与目标

• 推导在无闭式解的环面表面上运动的量子粒子的本征值与本征态。
• 解决尽管在受限量子力学与纳米管物理中具有重要性，但环面量子系统在教科书中仍缺乏标准解的问题。
• 基于傅里叶级数展开，发展一种系统性的数值方法，以求解环面上任意 m 量子数的薛定谔方程。
• 通过与微分方程的数值解对比，验证幂级数方法的有效性。

提出的方法

• 使用标准几何坐标 (R, a, θ, φ) 参数化环面，得到 θ 与 φ 坐标下的度量与拉普拉斯算符。
• 将哈密顿量表示为 θ 与 φ 坐标的函数，得到关于 θ 的二阶微分方程，并施加周期性边界条件。
• 将波函数 ψ(θ) 展开为 z = e^{iθ} 的傅里叶级数，导出系数 c_n 的递推关系。
• 对于 m=0，推导出 c_n 的三阶递推关系；对于 m≠0，推导出五阶递推关系。
• 通过施加周期性与宇称性（偶/奇）条件，约束系数关系，确保级数收敛。
• 通过在大 n 处对递推系数行列式条件求解，确定本征值，以保证非平凡解的存在。

实验结果

研究问题

• RQ1 能否构建一种系统性的幂级数方法，以在无闭式解的情况下求解环面上粒子的薛定谔方程？
• RQ2 m≠0 态的本征值与本征态在结构与收敛性方面，与 m=0 基态相比有何异同？
• RQ3 与微分方程的数值积分相比，傅里叶级数方法在多大程度上能实现高精度收敛？
• RQ4 宇称性与系数递推在确保级数解收敛过程中起到何种作用？
• RQ5 该方法能否推广至其他曲面几何或具有周期性边界条件的量子系统？

主要发现

• 对于 m=0，幂级数方法与龙格-库塔法数值解相比，收敛精度达到至少四位小数，其中 α=0.5。
• 计算并列出了 m=0、m=1 与 m=5 态的本征值，结果与数值解一致，显示良好收敛性。
• 通过施加 c_n = (-1)^n c_{-n} 与 c_n = (-1)^{n+1} c_{-n} 的关系，该方法成功处理了偶宇称与奇宇称解。
• 对于 m≠0，通过初始系数选择 (c₀=1, c₁=1) 与 (c₀=1, c₁=-1)，生成两个线性无关解，本征值通过求解行列式条件获得。
• m=2 的最低激发态 ψ₂₂ 的波函数形状与 m=0 基态 ψ₁₀ 定性相似，表明不同量子数下波函数结构具有相似性。
• 当级数截断至 n=6 时，精度达到至少四位有效数字，且在实际应用中通常使用更少项数即可满足需求。