Variational Multi-Gaussian Quantum Trajectories

技術要約：変分マルチガウス量子軌跡

問題提起
開いた量子多体系、特に相互作用するボゾン系の理論的研究は、局所的なヒルベルト空間が非有界であるという極めて困難な課題に直面している。占有数やシステムサイズが増大するにつれ、厳密なシミュレーションは急速に計算不能となる。切断ウィグナー近似（TWA）やガウス近似のような半古典的アプローチは効率的ではあるものの、相互作用や非ガウス相関が顕著になる場合には機能しない。さらに、量子軌跡は連続的なモニタリングの下での確率的発展を記述するための強力な枠組みを提供するが、その厳密なシミュレーションもまた、同様の多体系の複雑性に悩まされる。したがって、半古典的な限界を超えて非ガウス物理を捉えつつ、確率的な量子軌跡を効率的にシミュレートできる手法が必要とされている。

手法：確率的時間依存変分原理（sTDVP）
著者らは、**確率的変分マルチガウス（sVMG）**フレームワークを提案している。本手法の核となるのは、時間依存変分原理（TDVP）を確率的な量子状態拡散プロセスに対応させるための一般化である。

1. 確率的動力学： システムはリンドブラッド・マスター方程式に従い、これはヘテロダイン展開（量子状態拡散）を用いて量子軌跡へと「展開」される。これにより、伊藤型確率シュレディンガー方程式（SSE）が導かれる。変分解析を適用するために、著者らはこれを、変分多様体への投影に必要な通常の連鎖律を保持するストラトノヴィッチ型SSEへと変換する。
2. 変分アンザッツ： 波動関数 $|\psi\rangle$ は、$N_c$ 個のガウス波束のコヒーレントな重ね合わせとして近似される：
   $$\langle x | \psi[\theta] \rangle = \sum_{k=1}^{N_c} G(x; \theta_k)$$
   各ガウス成分は、実パラメータ $\theta$（平均位置・運動量、共分散行列、および重みを含む）によってパラメータ化される。大規模なシステムに対しては、パラメータ数を削減するために、共分散行列を「低ランク」分解によって圧縮することができる。
3. 運動方程式： 波動関数の厳密な確率的増分を、変分多様体の接空間に投影することにより、著者らはsTDVP方程式を導出する。これにより、以下のパラメータ $\theta$ に関する確率微分方程式が得られる：
   $$d\theta = \text{Re}(S^{-1}) \text{Re} \left( V dt + \sum_m Y^{(m)} \circ dZ_m^* \right)$$
   ここで、$S$ は量子幾何テンソルであり、$V$ および $Y^{(m)}$ はドリフトおよびノイズの力ベクトルである。決定的なのは、アンザッツがガウス関数で構成されているため、$S, V, Y^{(m)}$ を組み立てるために必要な積分（多項式で重み付けされたガウス関数を含む）が、ウィックの定理を用いて完全に解析的に計算可能である点である。

主な貢献

• TDVPの一般化： 本研究は、単なるユニタリ発展やリンドブラッド発展ではなく、連続測定によって駆動される確率的動力学に対して、時間依存変分原理を拡張したものである。
• 解析的組み立て： マルチガウス・アンザッツを用いることで、変分運動方程式を解析的に導出することが可能となり、波動関数自体のモンテカルロ・サンプリングを必要としない効率的な数値積分が可能となった。
• 制御された非ガウス性： 本手法は、ガウス成分の数（$N_c$）を変化させることで、波動関数の非ガウス性を系統的に制御する。$N_c$ が増加するにつれ、近似は厳密な動力学へと収束する。

結果とベンチマーク
本手法は、駆動・散逸型ボーズ・ハバード・モデルに対してベンチマーク適用された：

1. ボーズ・ハバード・ダイマー（$N=2$）：

   • 少数の成分（$N_c=4$）を用いたsVMGアプローチは、半古典的近似（TWA、ガウス軌跡近似）が失敗する領域において、個々の軌跡および絡み合い（エンタングルメント）の動力学を忠実に再現する。
   • 観測量の相対誤差は、$N_c$ の増加に伴い系統的に減少する。
   • 本手法は、非ガウス性が減少する熱力学的極限（大きな $N$ スケーリング）における対称性の破れの相転移およびエンタングルメント転移を、平均場予測と一致して捉えることに成功している。

2. 拡張ボーズ・ハバード格子（1Dおよび2D）：

   • 対称性の破れの転移（チェッカーボード・パターン）を調査するために、本フレームワークを1次元鎖および2次元格子に適用した。
   • 次元依存性： 2次元格子において、sVMGの結果は、明確に分離したチェッカーボード・サブ格子を持つ、明瞭な対称性の破れた相の出現を示している。対照的に、1次元鎖では、ドメイン壁によって分離された局所的なチェッカーボード状の領域が見られ、平均秩序パラメータが減衰した結果となっている。これは、対称性が破れた相が2次元では存続する一方で、1次元では揺らぎの増大によって抑制されることを示唆している。
   • 本手法は、低ランクアンザッツ（$N_c=4, K=2$）を用いてこれらの異なる挙動を捉えており、ブルートフォース的な厳密対角化の手が届かない領域へのスケーラビリティを示している。

意義
本論文は、sVMGフレームワークが、相互作用するボゾン系における確率的量子動力学をシミュレートするための、汎用性が高く効率的なツールを提供すると主張している。その主要な意義は、計算コストの低い半古典的手法と、厳密だが計算不可能なシミュレーションとの間のギャップを埋める能力にある。個々の軌跡に沿った非ガウス相関を捉えることで、本手法は、標準的なガウス近似や切断ウィクター近似では不十分な領域における、エンタングルメントの動力学や対称性の破れの相転移といった現象の研究を可能にする。運動方程式を解析的に組み立てられる能力は、開いた多体系の非平衡相を探索するための実用的なアプローチとなっている。