[論文レビュー] Numerical path integral approach to quantum dynamics and stationary quantum states

Tapio T. Rantala|arXiv (Cornell University)|Aug 17, 2016

Quantum, superfluid, helium dynamics参考文献 22被引用数 5

ひとこと要約

本稿では、コherentな時間発展演算子（Trotter核）と、新しい非coherentな時間発展演算法を用いて、実時間における量子ダイナミクスのシミュレーションおよび定常量子状態の計算を目的とした数値的経路積分法を提示する。モンテカルログリッドが規則的グリッドよりも数値的アーティファクトを回避することを示し、モンテカルロサンプリングを用いた非coherentな実時間発展演算により、高精度な基底状態エネルギーの計算が可能であることが示された。系統的誤差は統計的ノイズよりもTrotter近似に起因していることが明らかになった。

ABSTRACT

Applicability of Feynman path integral approach to numerical simulations of quantum dynamics in real time domain is examined. Coherent quantum dynamics is demonstrated with one dimensional test cases (quantum dot models) and performance of the Trotter kernel as compared with the exact kernels is tested. A novel approach for finding the ground state and other stationary sates is presented. This is based on the incoherent propagation in real time. For both approaches the Monte Carlo grid and sampling are tested and compared with regular grids and sampling. We asses the numerical prerequisites for all of the above. ar

研究の動機と目的

• 1次元系におけるフェインマン経路積分の実時間量子ダイナミクスシミュレーションへの適用可能性を評価すること。
• 経路積分シミュレーションにおける規則的グリッドとモンテカルログリッドの性能を比較すること。
• coherent時間発展演算におけるTrotter核近似の正確さを検証すること。
• 定常状態および基底状態を求めるために、新しい非coherentな実時間発展演算法を開発・検証すること。
• 安定的かつ正確なシミュレーションに必要な数値的条件（グリッドサイズ、時間刻み）を評価すること。

提案手法

• ハミルトニアンの時序指数関数を用いた数値的経路積分定式化により、時間発展演算子を計算する。
• coherentダイナミクスの近似時間発展演算子としてTrotter核を用いる。
• 空間グリッドのモンテカルロサンプリングを用いて、波動関数分布の非coherentな実時間発展演算を実装する。
• 配置空間をモンテカルログリッドでサンプリングし、規則的グリッドに見られる一般的な回折アーティファクトを回避する。
• 初期状態および発展演算後の波動関数分布を用いて、期待値（例えばポテンシャルエネルギー）を評価する。
• 標準偏差（σ）および平均の標準誤差（SEM）を用いた統計的分析により、結果の精度と正確さを推定する。

実験結果

リサーチクエスチョン

• RQ1 1次元量子ドットモデルにおけるcoherent量子ダイナミクスのシミュレーションにおいて、Trotter核近似はどの程度正確か？
• RQ2 経路積分シミュレーションにおいて、モンテカルログリッドが規則的グリッドを上回る数値的利点は何か？
• RQ3 モンテカルロサンプリングを用いた非coherentな実時間発展演算は、特に基底状態に安定して収束するか？
• RQ4 グリッドサイズおよび時間刻みは、基底状態エネルギー計算の正確さと安定性にどのように影響するか？
• RQ5 シミュレーションにおいて、Trotter近似に起因する系統的誤差が統計的誤差よりも支配的である程度はどの程度か？

主な発見

• モンテカルログリッド法により、coherent発展演算における規則的グリッドに見られる数値的回折アーティファクトが効果的に回避された。
• 調和ポテンシャルを有する1次元調和振動子（ODHO）において、Trotter核は基底状態エネルギーを高精度に再現し、3×10⁴グリッドと∆t=1を用いることで、偏差∆Vが30×10⁻⁶ auまで低下した。
• 統計的不確実性（σ）はサンプリング数の増加に伴い減少するが、実際の正確さは統計的ノイズではなくTrotter近似に起因する系統的誤差によって制限されている。
• 長時間シミュレーションでは、実際の正確さが統計的誤差推定値（SEM）よりも悪くなることが示され、時間離散化に起因する系統的誤差が支配的であることが明らかになった。
• 非coherent発展演算法は、波動関数分布および期待値の安定的収束を示し、定常状態の探索に有効であることが裏付けられた。
• より大きなグリッドサイズは正確さの向上に寄与し、時間刻み∆t=3がより短い時間刻みよりも良好な結果をもたらすという驚きの傾向が観察された。これは、Trotter核が時間刻みが小さい極限において示す挙動と整合的であった。