Efficient time-evolution of matrix product states using average Hamiltonians

技術要約：平均ハミルトニアンを用いた行列積状態の効率的な時間発展

問題提起
時間依存ハミルトニアンによって支配される量子多体系（QMBS）のシミュレーションは、ヒルベルト空間の指数関数的な増大により、依然として大きな課題となっている。テンソルネットワーク、特に行列積状態（MPS）は一次元系に対して効率的な枠組みを提供するが、標準的な時間発展アルゴリズム（TEBDやMPO $W_{II}$など）は、主に時間独立なハミルトニアンに対して最適化されている。時間依存のシナリオに適用する場合、標準的な実装はしばしば「リーマン・ステッパー（Riemann stepper）」に依存する。これは、時間間隔の開始時における瞬時ハミルトニアン（$H(0)$）を時間発展演算子に代入する手法である。このアプローチは高次インテグレータの精度を低下させ、時間依存の場合、発展演算子の収束率を（単一ステップにおける）3次から2次へと減少させてしまうため、与えられたステップサイズに対するシミュレーションの精度を制限する要因となる。

手法
著者らは、テンソルネットワーク・ライブラリの基本アーキテクチャを変更することなく、時間依存ハミルトニアンに対する高次の精度を回復させるための、既存のMPS時間ステッパーの修正案を提案している。核心となる手法は、単一ステップのプロパゲーター内で使用される瞬時ハミルトニアンを、数値積分理論から導出された、精緻に構成された「平均ハミルトニアン」に置き換えることである。

具体的には、著者らは標準的なリーマン近似の代わりに、シンプソンの公式を用いて計算されたハミルトニアンを代入する：
$$H_{\text{Simp}}(\Delta t) \equiv \frac{1}{6} \left( H(0) + 4H(\Delta t/2) + H(\Delta t) \right)$$
この平均ハミルトニアンは、標準的な2次TEBDまたはMPO $W_{II}$アルゴリズムに供給される。理論的解析によれば、リーマン・ステッパー（$H(0)$）は誤差界 $O(\Delta t^2)$ の発展演算子をもたらすのに対し、シンプソン・ステッパーは基礎となるインテグレータが本来持つ3次の局所精度（$O(\Delta t^3)$）を維持する。この証明は、Hölder-$\infty$ ノルムを用いて、厳密な時間発展演算子と近似演算子の差を境界付けることに基づいており、シンプソン近似がリーマン近似と比較して積分誤差項を大幅に低減することを示している。

主要な貢献と結果
提案された手法は、ダイヤモンド中の窒素空孔（NV）センターを表すスピン-1サイトの鎖の数値シミュレーションを通じて検証された。この系は、実験データから導出されたパラメータを用いて、ドリフト項、相互作用項、および制御項を含むハミルトニアンでモデル化されている。

1. 精度の向上： シミュレーションでは、様々なシステムサイズ（3から8個のNVセンター）および積分ステップ数にわたって、「シンプソン・ステッパー」を標準的な「リーマン・ステッパー」と比較した。

   • 小規模な系（例：3個のNVセンター）では、中程度のステップサイズにおいて、シンプソン・ステッパーはリーマン・ステッパーと比較して平均誤差を約1000分の1に低減した。
   • より大きな系（最大8個のNVセンター）においても、本手法は一貫して低い伝播誤差を示し、ステップ数やシステムサイズに応じて、平均誤差の改善率は50倍から4000倍以上に達した。
   • 誤差比（$E_{\text{Riemann}} / E_{\text{Simpson}}$）はテストされたすべての構成において高い値を維持しており、本手法の堅牢性を裏付けている。

2. 計算コスト： 精度の向上に伴い、実行時間はわずかに増加する。シンプソン・ステッパーは、リーマン・ステッパーの1回に対し、1ステップあたり3回のハミルトニアン評価（$t=0$, $t=\Delta t/2$, $t=\Delta t$）を必要とする。著者らは、この実行時間の増加は軽微であり、純粋に追加の関数評価に起因するものであると述べており、精度に対するトレードオフは非常に有利であるとしている。

3. スケーラビリティ： 本手法は、システムサイズが増加しても安定した性能向上を示した。結合次元が切断誤差を無視できるほど十分に大きい限り、新たな切断誤差を導入することなく、標準的なアプローチに対する優位性を維持した。

意義と限界
本論文は、このアプローチが既存のMPSソルバーに対する、単純かつ非常に効率的な拡張であることを主張している。古典的な数値積分のアイデア（シンプソンの公式）を活用することで、時間依存のダイナミクスをシミュレートする際に失われがちな、インテグレータの高次収束特性を回復させている。これは、量子制御、分光学、あるいはフロケ・エンジニアリングのように、時間依存の駆動場が遍在するアプリケーションにおいて特に重要である。

著者らは、いくつかの制限事項と実用上の留意点を明示している：

• 滑らかさの要件： 本手法は、各区間内におけるハミルトニアン $H(t)$ の第4次までの微分が存在し、かつ有界であることを前提としている。
• ハミルトニアン構成のコスト： 時間依存ハミルトニアンの構築に計算コストが高い場合、1ステップあたり3回の評価が必要となることが計算コストを大幅に増大させる可能性がある。
• もつれの限界： 本手法は積分誤差を改善するが、大規模な多体系におけるエンタングルメント（量子もつれ）の増大による制限を排除するものではない。これらはMPS表現における切断誤差によって支配される。

著者らは、シンプソンの平均化の概念は適応可能であり、将来の研究において、さらに高次の有効ハミルトニアンを構築するために、より高次の求積則（例：ガウス＝ルジャンドル公式）や多段組成法と組み合わせることが可能であると示唆している。