[论文解读] Exact diagonalization of the $d-$dimensional confined quantum harmonic oscillator

Kunle Adegoke, Adenike Olatinwo+3 more|arXiv (Cornell University)|Apr 5, 2016

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

一句话总结

本文提出了一种d维受限量子简谐振子哈密顿量的精确矩阵表示，以动能算符的本征态为基函数。通过构建动能和势能算符的显式N^d × N^d矩阵，该方法实现了无需根查找或数值积分的直接对角化，从而获得精确的能量谱。关键贡献在于推导出一个闭式矩阵表达式，使计算更加高效，并实现了1D情形下的六阶解析近似。

ABSTRACT

In the existing literature various numerical techniques have been developed to quantize the confined harmonic oscillator in higher dimensions. In obtaining the energy eigenvalues, such methods often involve indirect approaches such as searching for the roots of hypergeometric functions or numerically solving a differential equation. In this paper, however, we derive an explicit matrix representation for the Hamiltonian of a confined quantum harmonic oscillator in higher dimensions, thus facilitating direct diagonalization.

研究动机与目标

• 开发d维受限简谐振子哈密顿量的直接、精确矩阵表示。
• 消除对超几何函数中的根查找或求解微分方程等间接数值方法的依赖。
• 通过直接对角化实现高效且精确的能量谱计算。
• 从矩阵形式推导出1D受限简谐振子能级的高阶解析近似表达式，最高至ω的六阶。
• 为具有硬壁限制的高维量子系统提供计算上可行的框架。

提出的方法

• 基态由1维动能算符本征函数构造，表达式为φ_r(x_i) = √(1/L) cos[π/2 sin²(rπ/2) - (r+1)πx_i/(2L)]，构成完备正交基。
• d维希尔伯特空间由直积态ψ_s(x) = ∏_{i=1}^d φ_{s_i}(x_i)张成，其中s索引N^d个态。
• 哈密顿量矩阵元表达为H_st = ∑_{i=1}^d c_{i_st} (T_i)_{s_i t_i} + ∑_{i=1}^d c_{i_st} (V_i)_{s_i t_i}，其中c_i为二值矩阵，用于在各维度间强制索引匹配。
• 动能矩阵T在该基下为对角矩阵，本征值为ε_s = ε ∑_{i=1}^d (s_i + 1)^2，其中ε = π²ℏ²/(8mL²)。
• 势能矩阵V通过矩阵元(V_i)_{s_i t_i} = ⟨φ_{s_i}|(1/2)mω²x_i²|φ_{t_i}⟩构造，利用三角积分计算。
• 总哈密顿量表示为H = T + V，其矩阵元以克罗内克δ结构和余弦项表示，以捕捉非对角耦合。

实验结果

研究问题

• RQ1 能否在不依赖间接数值方法的前提下，推导出d维受限简谐振子的精确矩阵表示？
• RQ2 动能和势能算符如何在受限态的有限维基中显式表示？
• RQ3 在动能算符本征态基下，d维系统的哈密顿量矩阵具有何种结构？
• RQ4 能否从矩阵形式中推导出1D受限简谐振子能量谱的高阶解析近似？
• RQ5 二值耦合矩阵c_i在将总哈密顿量分解为解耦的1D贡献中起到何种作用？

主要发现

• 哈密顿量以显式N^d × N^d矩阵形式表达，其矩阵元包含δ函数、余弦项和反平方项，支持直接对角化。
• 在1D情形下，矩阵退化为N×N形式，当s≠t时，H_st = εδ_st(s+1)^2 + (λ²ε/8)δ_st(π²/6 - 1/(s+1)^2) + (λ²ε/2)(1/(s−t)^2 + δ_st − 1/(s+t+2)^2)cos²((s−t)π/2)。
• 推导出最高至ω六阶的近似解析能量谱，表达式为E_r ≈ ∑_{m=0}^3 {λ^{2m}ε / 2^{4m−1} ∑_{n=0}^m [(-1)^n ζ(2m−2n)c_n^{(m)} / (r+1)^{2(m+n−1)}]}
• 二阶能量修正为E_r^{(2)} = (λ⁴ε/128)[ζ(4)/(r+1)^2 − 5ζ(2)/(r+1)^4 + 7/(r+1)^6]，其中ζ(2)=π²/6且ζ(4)=π⁴/90。
• 该方法避免求解微分方程或在特殊函数中进行根查找，为受限量子系统提供了数值稳定且高效的替代方案。
• 该方法可推广至任意d和N，其矩阵结构保持原始哈密顿量的克罗内克和形式。