[论文解读] Numerical study of quantum mechanical systems using a quantum wave impedance approach

Orest Hryhorchak|arXiv (Cornell University)|Oct 16, 2020

Nonlinear Photonic Systems参考文献 15被引用 4

一句话总结

本文提出了一种基于量子波阻抗的数值方法，用于求解量子力学系统，证明了将复杂势能近似为分段常数区域可实现精确且高效的计算。该方法验证了该近似方法的有效性，量化了精度随离散化步长的变化关系，并提出了一种将零范围δ势能引入量子波阻抗框架的技术，显著提升了传统方法（如有限差分法或转移矩阵法）的计算效率。

ABSTRACT

The approximate numerical method for a calculation of a quantum wave impedance in a case of a potential energy with a complicated spatial structure is considered. It was proved that the approximation of a real potential by a piesewise constant function is also reasonable in a case of using a quantum impedance approach.The dependence of an accuracy of numerical calculations on a number of cascads by which a real potential is represented was found. The method of including into a consideration of zero-range singular potentials was developed.

研究动机与目标

• 验证在复杂势能中使用分段常数势能近似在量子波阻抗方法中的有效性。
• 量化数值计算精度随势能段数（级联数）的变化关系。
• 开发一种将零范围奇异势能（如δ函数）整合进量子波阻抗框架的方法。
• 展示迭代量子波阻抗方法相较于传统方法在计算效率和数值稳定性方面的优势。

提出的方法

• 将任意非奇异势能近似为N个分段常数势能区域的级联，各区域高度为V_i，宽度为Δx_i。
• 应用迭代量子波阻抗公式：Z(x - Δx_i) = z_i [Z(x) - z_i tanh(γ_i Δx_i)] / [z_i - Z(x) tanh(γ_i Δx_i)]，其中z_i为特征阻抗，γ_i为区域i中的波矢。
• 推导迭代公式的连续极限，以恢复量子波阻抗的一阶微分方程：Z'(x) + (im/ℏ)Z²(x) = (2i/ℏ)(E - U(x))。
• 通过在势能段上逐步求解阻抗方程，利用迭代方法计算透射系数和束缚态能级。
• 引入数值误差度量：平均精度 ε̄(N) = (1/n) Σ|T^(N)(E_j) - T^(N/2)(E_j)|，以评估随N增加的收敛性。
• 通过在段间界面处放置δ势能并相应调整匹配条件，将方法扩展至包含δ势能。

实验结果

研究问题

• RQ1 当使用量子波阻抗方法时，将复杂势能近似为分段常数函数是否能获得准确结果？
• RQ2 数值解的精度如何随近似中使用的势能段数（级联数）变化？
• RQ3 零范围奇异势能（δ函数）能否被有效整合进量子波阻抗框架？
• RQ4 与有限差分法或转移矩阵法相比，量子波阻抗方法在量子系统中的效率和精度如何？

主要发现

• 在使用量子波阻抗方法时，对真实势能进行分段常数近似是有效且准确的，且随着段数增加，精度不断提高。
• 透射系数计算的平均误差随势能段数增加而减小，收敛趋势可预测，如图2所示，ε̄(N)随log₂(N)减小。
• 通过在段边界处放置δ势能，该方法成功实现了δ势能的整合，从而可精确模拟具有点相互作用的系统。
• 量子波阻抗方法仅需N−1次迭代循环即可完成N个段的计算，显著优于转移矩阵法（需N次矩阵乘法）或有限差分法（需求解大型方程组）。
• 当段尺寸远小于电子波长时，该方法可达到与解析结果相当的高精度，验证了其在共振隧穿和束缚态计算中的适用性。
• 该方法可推广至无限和半无限周期系统，作者结合前期工作已提出此可能性。