[論文レビュー] 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行列を構築することにより、根の探索や数値積分を必要とせず、直接対角化が可能となり、高精度なエネルギー準位が得られる。主な貢献は、計算を効率化する閉形式の行列式であり、1次元の場合の6次までの解析的近似が可能となる。

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次元の閉じ込められた調和振動子ハミルトニアンの直接的かつ正確な行列表現の開発。
• 超幾何関数の根の探索や微分方程式の数値的解法といった間接的数値手法の排除。
• 直接対角化による、効率的かつ高精度なエネルギー準位の計算の実現。
• 1次元の閉じ込められた調和振動子のエネルギー準位に対する6次までの解析的近似式の導出。
• 硬い壁による閉じ込めを伴う高次元量子系に対して、計算的に取り扱いやすいフレームワークの提供。

提案手法

• 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は次元間のインデックス一致を制御する2値行列である。
• 運動エネルギー行列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として構成され、Kroneckerデルタ構造と余弦項を用いて非対角結合を捕捉する形で行列要素が表現される。

実験結果

リサーチクエスチョン

• RQ1 間接的数値手法に依存せずに、d次元の閉じ込められた調和振動子の正確な行列表現を導出可能か？
• RQ2 閉じ込められた状態の有限次元基底において、運動エネルギーとポテンシャルエネルギー演算子を明示的に表現できるか？
• RQ3 d次元系において、運動エネルギー固有状態の基底におけるハミルトニアン行列の構造はいかなるものか？
• RQ4 行列形式から、1次元の閉じ込められた調和振動子のエネルギー準位に対する高次の解析的近似を導出可能か？
• RQ5 2値結合行列c_iは、全ハミルトニアンを1次元の独立した寄与に分解可能にする役割を果たすか？

主な発見

• ハミルトニアンは、δ関数、余弦項、逆2乗項を含む閉形式の要素を持つ明示的なN^d × N^d行列として表現され、直接対角化が可能となる。
• 1次元の場合、行列は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) が成り立つ。
• ωについて6次までの解析的エネルギー準位の近似式が導出され、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)}] } となる。
• 2次エネルギー補正はE_r^{(2)} = (λ⁴ε/128)[ζ(4)/(r+1)^2 − 5ζ(2)/(r+1)^4 + 7/(r+1)^6] であり、ζ(2)=π²/6 および ζ(4)=π⁴/90 である。
• 微分方程式の解法や特殊関数における根の探索を回避し、閉じ込められた量子系に対して数値的に安定かつ効率的な代替手法を提供する。
• 本手法は任意のdおよびNに一般化可能であり、元のハミルトニアンのKronecker和形式を保持する行列構造を維持する。