Action on the Sphere: An Interfering Mean-Field Propagator for the Bose-Hubbard Dimer

技術要約：球面上での作用：ボース・ハバード・ダイマーのための干渉型平均場プロパゲーター

問題提起
ボース・ハバード模型、特に二モード（ダイマー）の場合、自己トラッピング、崩壊、リバイバル、および多粒子トンネリングを含む複雑な多体ダイナミクスを示す。SU(2)コヒーレント状態から導かれる平均場近似は古典的極限を捉えることができるが、リバイバルやトンネリングといった量子干渉効果を再現するには不十分である。既存の半古典的アプローチには限界がある。例えば、切断されたウィグナー近似（Truncated Wigner Approximation）は崩壊を捉えることができるが位相コヒーレンス（リバイバル/トンネリング）を逃し、Glauberコヒーレント状態に基づくヘルマン・クルック（Herman-Kluk）プロパゲーターは正確ではあるものの計算負荷が高く、単一モード系を超えた拡張が困難である。また、結合コヒーレント状態法は原理的には厳密であるが、計算コストが急速に増大するという課題がある。著者らは、結合方程式の重いオーバーヘッドなしに、崩壊、リバイバル、およびトンネリング現象を捉えるための位相情報を保持した、計算効率の高い半古典的プロパゲーターの開発を目指している。

手法
著者らは、球面上におけるSU(2)コヒーレント状態の過完全性に基づいた干渉型平均場（Interfering Mean-Field: IMF）プロパゲーターを構築している。手法の手順は以下の通りである：

1. 初期値表現（Initial Value Representation）： 正確な時間発展演算子を、初期のSU(2)コヒーレント状態 $|z_0, \phi_0\rangle$ の積分として表現する。
2. 位相を伴う平均場近似： 全量子状態を時間発展させる代わりに、各初期コヒーレント状態を、時間依存変分原理（SCS-TDVP）から導かれる古典的平均場運動方程式を用いて発展させる。極めて重要な点は、各軌道に対応する平均場作用量から、各軌道の位相を追跡することである。
3. ダイナミクスにおける古典的極限： 不確定性の過剰カウントを避けるため、著者らはフル量子期待値ではなく、ハミルトニアンの厳密な古典的極限（$H = \lim_{j\to\infty} \langle \hat{H} \rangle / j$）を用いて軌道を生成する。
4. 数値実装： 球面上の積分は、ガウス・ルジャンドル求積法を用いて離散化される。
5. タイムスライシング（Time-Slicing）への拡張： 単一ステップのIMFでは（長時間の位相コヒーレンスを必要とする現象である）多粒子トンネリングを捉えられないという課題に対処するため、著者らはタイムスライス型（TS）IMFプロパゲーターを導入している。これは、短い時間間隔 $\delta t$ ごとにIMFプロパゲーターを適用し、得られた状態をコヒーレント状態基底で再展開し、このプロセスを繰り返すものである。
6. スケーリング因子の補正： 半古典的極限における純粋な相互作用項の解析的な検討により、タイムスライシングの手順が、相互作用強度 $\kappa$ に対するスケーリング因子として現れる累積誤差を導入することが明らかになった。著者らは、粒子数 $N$ が大きい場合（スピン $j=N/2$）の漸近的なスケーリング因子 $\kappa_{eff}/\kappa = (1 + 1/2j)^{-1}$ を導出し、より小さな $N$ に対する最適な数値的係数を決定した。

主な貢献と結果

• 崩壊とリバイバルの再現： IMFプロパゲーターは、ボース・ハバード・ダイマーの崩壊およびリバイバル・ダイナミクスを正常に再現する。これには、厳密な結果とほぼ区別がつかないほど微細な量子構造を示す、時間発展した全ウィグナー関数も含まれる。これは、比較的少ない粒子数（例：$j=10$）であっても達成される。
• タイムスライシングによるトンネリングの捕捉： 標準的なIMFは、自己トラッピング状態間の多粒子トンネリングを再現できない。タイムスライス型（TS-IMF）は、このトンネリング挙動を高精度に回復することに成功している。
• スケーリング因子の発見： 著者らは、TS-IMFプロパゲラーが、厳密なダイナミクスと一致させるために相互作用パラメータ $\kappa$ の再スケーリングを必要とすることを実証した。
  • $N$ が大きい場合、最適なスケーリングは $\kappa_{eff} = \kappa (1 - \frac{1}{2j+1})$ として解析的に導かれる。
  • $N$ が小さい場合、数値的な最適化により、半古典的極限から遠く離れた状態でも厳密なダイナミクスを再現できる特定のスケーリング因子（例：$j=50$ で 0.991、$j=2.5$ で 1.018）が得られる。
• 計算効率： 結合コヒーレント状態法とは異なり、IMFおよびTS-IMF法は成分状態間の結合運動方程式を含まないため、高度に並列化が可能である。また、タイムスライシング法を用いることで、単一ステップの長時間伝播と比較して、より小さなグリッドサイズの使用が可能となる。

意義と主張
本論文は、「単純な初期値コヒーレント状態プロパゲーター」を提供し、単純な平均場近似と、計算負荷の高いフル量子または結合コヒーレント状態法の間のギャップを埋めるものであると主張している。著者らは以下のように述べている：

• 本手法は非相互作用系に対して厳密であり、相互作用系において崩壊とリバイバルを捉える上で非常に効果的である。
• タイムスライシング拡張と相互作用スケーリング因子を用いることで、半古典的近似が通常失敗するような小さな粒子数においても、多粒子トンネリング・ダイナミクスを正確に複製できる。
• 本論文の焦点はボース・ハバード・ダイマーであるが、本手法は一般的なSU(2)ハミルトニアンへ容易に拡張可能であり、理論的にはSU(M)系へも拡張可能である（ただし、高次元ではグリッド実装が困難になるため、モンテカルロ・サンプリングが代替案として示唆されている）。

著者らは、プロパゲーターの非ユニタリな性質（ノルムの減少を招く）については、半古典的プロパゲーターに共通する特徴であるとして控えめに言及しているが、期待値およびウィグナー関数は定量的に正確であることを強調している。