[论文解读] Application of the variational $R$-matrix method to one-dimensional quantum tunneling

Joseph Kimeu, Roland Mai+1 more|ArXiv.org|Jul 29, 2004

Theoretical and Computational Physics参考文献 5被引用 3

一句话总结

本文将变分R-矩阵方法应用于计算五种势垒剖面——截断线性、指数、抛物线、钟形及Eckart势垒——的一维量子隧道效应的透射与反射系数。该方法被证明是一种简单、非迭代且精确的替代方案，与解析解及既有的数值方法在所有测试势垒上均表现出高度一致。

ABSTRACT

We have applied the variational $R$-matrix method to calculate the reflection and tunneling probabilities of particles tunneling through one-dimensional potential barriers for five different types of potential profiles -- truncated linear step, truncated exponential step, truncated parabolic, bell-shaped, and Eckart. Our variational results for the transmission and reflection coefficients are compared with exact analytical results and results obtained from other numerical methods. We find that our results are in good agreement with them. We conclude that the variational $R$-matrix method is a simple, non-iterative, and effective method to solve one-dimensional quantum tunneling problems.

研究动机与目标

• 评估变分R-矩阵方法在计算一维势垒隧道概率方面的有效性。
• 将该方法的结果与精确解析解及其他数值技术（如WKB、MAT和TM方法）进行比较。
• 评估该方法在多种势垒形状（包括不连续和光滑变化的势垒）下的精度与效率。
• 确定基组优化是否能改善变分R-矩阵框架中的收敛性与精度。
• 确立该方法作为现有隧道计算技术的通用、非迭代且计算高效的替代方案。

提出的方法

• 变分R-矩阵方法通过将空间划分为内区（a ≤ x ≤ b）和外区（x ≤ a，x ≥ b）来求解一维定态薛定谔方程，其中内区采用变分法处理。
• 内区波函数展开为余弦基函数 χ_i(x) = cos(κ_i x)，其系数 c_i 通过变分法确定，以满足薛定谔方程和边界条件。
• 该方法从薛定谔方程推导出广义特征值问题，形成包含动能、势能及边界导数项（A、Δ^a、Δ^b 矩阵）的矩阵方程。
• 利用对数导数 λ(a) = ψ′(a)/ψ(a) 和 λ(b) = ψ′(b)/ψ(b) 实现边界 a 和 b 处波函数及其导数的连续性匹配。
• 通过内区解的相移与归一化，从外区波函数的渐近行为中提取透射与反射系数。
• 该方法应用于五类势垒：截断线性势、指数势、抛物线势、钟形势及Eckart势，结果与精确解析表达式及数值基准进行对比。

实验结果

研究问题

• RQ1 变分R-矩阵方法在不同势垒剖面下计算一维量子隧道效应的透射与反射系数时，精度如何？
• RQ2 该方法的性能与精确解析解及其他成熟数值方法（如WKB、MAT和TM）相比如何？
• RQ3 通过优化基组参数（κ_i）在多大程度上可提升变分R-矩阵方法的精度？
• RQ4 对于WKB方法失效的具有不连续性或快速变化形状的势垒，该方法是否能有效应用？
• RQ5 变分R-矩阵方法是否可作为现有隧道计算技术的可行、非迭代且计算高效的替代方案？

主要发现

• 变分R-矩阵方法在所有五种测试势垒剖面下，所得透射与反射系数与精确解析结果高度一致。
• 对于截断线性阶跃势，当使用优化基组时，误差低于0.05%，结果与精确解及数值矩阵方法完全一致。
• 对于Eckart势，当能量E略高于势垒峰顶（2.23 hartree）时，透射系数可达约98%，与精确值极为接近。
• 在Eckart势中，该方法的能量本征求解最大误差为1.43%，预计通过基组优化可进一步改善。
• 即使在复杂势垒（如钟形势和截断抛物线势）下，该方法仍保持高精度，透射系数与精确值的偏差在百分之几以内。
• 变分R-矩阵方法为非迭代、计算高效，适用于任意势垒形状，在非缓慢变化或不连续势垒中，其精度优于半经典WKB及改进型WKB方法。