[論文レビュー] 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

ひとこと要約

本稿では、5つのポテンシャルプロファイル（切り捨て線形、指数関数的、放物線型、ベル型、Eckart障壁）を対象に、変分R行列法を用いて一次元量子トンネルにおける透過係数と反射係数を計算する。この手法は、解析的および数値的手法と比較して、単純で反復処理を要せず、高い正確性を示し、すべてのテストされたポテンシャルで正確な結果および既存の数値的手法と強い一致を示した。

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 および b における波動関数と微分の連続性を満たすために、対数微分 λ(a) = ψ′(a)/ψ(a) および λ(b) = ψ′(b)/ψ(b) が用いられる。
• 透過係数および反射係数は、外側領域における波動関数の漸近的挙動から、内側領域の解の位相シフトおよび正規化を用いて抽出される。
• この手法は、切り捨て線形、指数関数的、放物線型、ベル型、Eckart という5種類のポテンシャルに適用され、正確な解析的式および数値的ベンチマークと比較された。

実験結果

リサーチクエスチョン

• RQ1 変分R行列法は、多様なポテンシャルプロファイルにおける一次元量子トンネルの透過係数と反射係数をどの程度正確に計算できるか？
• RQ2 この手法の性能は、正確な解析解およびWKB、MAT、TM法などの既存の数値的手法と比較してどうなるか？
• RQ3 基底セットパラメータ（κ_i）の最適化によって、変分R行列法の正確性をどの程度向上させられるか？
• RQ4 WKB法が失敗する不連続性や急激に変化する形状を有するポテンシャルに対しても、この手法は効果的に適用可能か？
• RQ5 変分R行列法は、実用的で反復処理を要せず、計算的に効率的な既存のトンネル計算手法の代替として有効か？

主な発見

• 変分R行列法は、5つのテストされたポテンシャルプロファイルすべてにおいて、正確な解析的結果と良好な一致を示す透過係数および反射係数を生成した。
• 切り捨て線形ステップポテンシャルに対しては、最適化された基底セットを用いることで、誤差が0.05%未満に抑えられ、正確な解および数値的行列法の結果と一致した。
• Eckartポテンシャルに対しては、エネルギーEが障壁ピークよりわずかに高い場合（2.23ハートリー）、透過係数が約98%に達し、正確な値とよく一致した。
• Eckartポテンシャルにおいて、エネルギー固有値の最大誤差は1.43%にとどまり、基底セットの最適化によりさらなる改善が期待された。
• ベル型や切り捨て放物線型障壁のような複雑なポテンシャルに対しても、この手法は高い正確性を維持し、透過係数が正確な値と数パーセント以内で一致した。
• 変分R行列法は反復処理を要せず、計算的に効率的であり、任意のポテンシャル形状に適用可能で、非緩徐変化または不連続なポテンシャルに対しては、半古典的WKB法や修正WKB法を上回る正確性を示した。