[論文レビュー] On the Convergence of Time Splitting Methods for Quantum Dynamics in the Semiclassical Regime

François Golse, Shi Jin+1 more|arXiv (Cornell University)|Jun 8, 2019

Numerical methods for differential equations参考文献 24被引用数 21

ひとこと要約

本稿は、半古典的領域におけるフォン・ノイマン方程式に対する時間分割法（Lie-TrotterおよびStrang）について、プランク定数に依存しない誤差推定を確立する。量子密度作用素と古典的位相空間密度の間の2階 Wasserstein 距離に類似した擬距離を用いることで、数値的収束性が $\hat{\hbar}$ に依存しないことを証明し、$\hbar \to 0$ であっても $\hat{\Delta t} = O(1)$ の時間ステップが可能であることを示す。主な結果は、$\hbar$ に一様に有効な明示的な誤差境界であり、半古典的量子ダイナミクスのシミュレーションにおける長年の課題を解決する。

ABSTRACT

By using the pseudo-metric introduced in [F. Golse, T. Paul: Archive for Rational Mech. Anal. 223 (2017) 57-94], which is an analogue of the Wasserstein distance of exponent $2$ between a quantum density operator and a classical (phase-space) density, we prove that the convergence of time splitting algorithms for the von Neumann equation of quantum dynamics is uniform in the Planck constant $\\hbar$. We obtain explicit uniform in $\\hbar$ error estimates for the first order Lie-Trotter, and the second order Strang splitting methods.

研究の動機と目的

• 量子ダイナミクスのシミュレーションにおいて、$\sim \hbar$ のオーダーの小さなド・ブロイ波長を解像する計算的課題に対処すること。
• 時間分割スペクトル法における $O(1)$ 時間ステップの使用を、形式的レベルから厳密に正当化すること。
• プランク定数 $\hbar$ に一様に有効な誤差推定を確立することにより、半古典的領域における信頼性の高い数値的シミュレーションを可能にすること。
• 混合量子状態の半古典的振る舞いを捉える位相空間距離を用いて、量子ダイナミクスと古典的極限の間のギャップを埋めること。

提案手法

• 著者らは、量子密度作用素と古典的位相空間確率密度の間で定義される、2階 Wasserstein 距離に類似した擬距離を用いる。
• この距離は、時間分割法によって得られた量子状態のウィグナー変換と、古典的輸送解とを比較するために適用される。
• 解析は、量子進化を位相空間における古典的リウヴィルダイナミクスに結びつけるためにウィグナー変換に依存する。
• トレースクラスおよびヒルベルト＝シュミット作用素ノルムを用い、量子および古典的測度の間のカップリング技術を組み合わせて誤差推定が導かれる。
• 指数型分割スキーム（Lie-TrotterおよびStrang）の既知の作用素ノルムの境界を、混合状態のためのフォン・ノイマン方程式に適用する。
• 重要な技術的道具は、量子および古典的分布の間のカップリングを構築するための、 coherent states $|q,p\rangle$ を基底とする手法である。

実験結果

リサーチクエスチョン

• RQ1 フォン・ノイマン方程式に対する時間分割法が、半古典的領域において $\hbar$ に依存しない時間ステップで収束することができるか？
• RQ2 $\hbar$ に一様な誤差推定が可能となる適切な距離は存在するか？
• RQ3 Lie-TrotterおよびStrang分割スキームの誤差は、$\hbar \to 0$ のときどのように振る舞うか？
• RQ4 ウィグナー変換フレームワークを用いて、時間分割スペクトル法における $O(1)$ 時間ステップの使用を厳密に正当化できるか？
• RQ5 半古典的極限において、誤差は $\hbar$、$\Delta t$、初期データにどのように依存するか？

主な発見

• 1次精度の Lie-Trotter 分割法は、$\hbar$ に一様な誤差境界 $\Delta t + \min(\frac{\Delta t}{\hbar}, \sqrt{\hbar})$ を達成し、$\hbar = \Delta t^{4/3}$ のとき $\Delta t + \Delta t^{1/3}$ に簡略化される。
• 2次精度の Strang 分割法は、$\hbar$ に一様な誤差境界 $\Delta t + \min(\frac{\Delta t^2}{\hbar}, \sqrt{\hbar})$ を達成し、$\hbar = \Delta t^{4/3}$ のとき $\Delta t + \Delta t^{2/3}$ に簡略化される。
• 誤差推定は $\hbar$ に一様であるため、$\hbar \to 0$ であっても収束速度が劣化しない。これは、半古典的シミュレーションにおいて極めて重要である。
• 使用された擬距離は、2階 Wasserstein 距離に類似した量子-古典的アナログであり、量子密度作用素と古典的位相空間密度の間の比較を可能にする。
• 解析により、標準的な $L^2$ またはソボレフノルムでは誤差が $O((\Delta t / \hbar)^m)$ に比例するため、$O(1)$ 時間ステップが不可能であることが示されたが、本研究で選ばれた距離はこの問題を回避する。
• 本研究の結果は、半古典的領域における時間分割スペクトル法の経験的成功を裏付け、$\Delta t = O(1)$ を用いた使用の厳密な根拠を提供する。