Quantum Brownian Motion as a Classical Stochastic Process in Phase Space

技術要約：位相空間における古典的確率過程としての量子ブラウン運動

問題提起
量子ブラウン運動（QBM）のシミュレーションは、統計物理学における極めて困難な課題であり続けている。カルデイラ・レゲット・モデルは、調和振動子のバスに結合した粒子の微視的に一貫した記述を提供するが、ヒルベルト空間の次元の指数関数的な増大が、厳密な数値シミュレーションを阻害している。既存の手法には、以下のような重大な制限がある：

• 高温近似： マスター方程式（例：リンドブラッド、レドフィールド）は、多くの場合、高温極限（$T \gg \hbar\Lambda$）に限定されており、零点揺らぎのような本質的な量子特性を捉えることができなかったり、低温において完全正値性を損なったりすることがある。
• 厳密な手法： パス積分解は、任意の状態準備や中間介入を伴う場合、解析的に煩雑である。階層的方程式（HEOM）や確率的シュレディンガー方程式は、数値的に厳密ではあるものの、バスの相関の代数的な減衰が遅い場合や、状態空間の次元が大きい場合、特に低温において計算コストが増大する。
• 核心的な対立： 二次ポテンシャルを持つカルデイラ・レゲット・モデルは可積分であるが、任意の温度におけるその量子シミュレーションは依然として非自明である。あらゆる温度領域において、厳密かつ計算効率の高い簡約記述が必要とされている。

手法
著者らは、カルデイラ・レゲット・モデル（高々二次の外部ポテンシャルを持つもの）の厳密な量子力学と、位相空間における古典的な非マルコフ的確率過程との間の厳密なマッピングを確立した。

1. 理論的枠組み：

   • 粒子とバスとの間の相関のある熱平衡状態から出発し、著者らはウィグナー関数表現を利用する。
   • 二次ハミルトニアンの場合、ウィグナー関数のための量子リウヴィル方程式は、形式的に古典リウヴィル方程式と同一である。これにより、全系の発展は、全位相空間における決定論的な軌跡のアンサンブルによって表現することが可能になる。
   • バス自由度を解析的に排除することで、粒子の動力学が以下の一般化ランジュバン方程式（GLE）に従うことが示される：
     $$m\ddot{r}(t) = -V'(r(t)) - \int_{-\infty}^{t} M(t-\tau)\dot{r}(\tau)d\tau + \xi_W(t)$$
   • 「量子性」は、環境の統計量の中に完全にエンコードされている。このノイズ $\xi_W(t)$ は平均ゼロのガウス型であり、その自己相関関数は、量子揺らぎ散逸定理（FDT）から導出され、フル量子スペクトル密度と温度依存性（零点揺らぎを含む）を組み込んでいる。

2. 状態準備と介入の取り扱い：

   • 本形式は、高度に非古典的な重ね合わせを含む任意の初期量子状態を、そのウィグナー関数を通じて受け入れることができる。
   • 任意の時刻 $t_k$ における介入（測定、ユニタリ・キック）は、準備演算子のウィグナー表現に基づき、軌跡に複素重みを割り当てることで、準備関数 $\lambda$ を用いてモデル化される。
   • 期待値は、軌跡のアンサンブルに対する重み付き平均として計算され、これにより複雑な駆動・散逸プロトコルのシミュレーションが可能となる。

3. 数値実装：

   • 著者らは、以下の要素を含むモンテカルロ法を開発した：
     • ノイズ生成： 量子FDTを満たすように、周波数領域でカラーノイズの軌跡を合成する。
     • 平衡初期化： 任意の初期条件からの長時間平衡化、または応答関数を用いた明示的な生成を通じて、相関のある熱平衡状態に対応する軌跡を生成する。
     • 重み付きサンプリング： 状態準備および介入を考慮するために、再重み付けおよびリサンプリングの手順を適用する。

主要な結果
本論文は、既知の解析解に対して本手法の妥当性を検証し、特定の領域において高温近似に対する優位性を実証している：

• 検証： 自由粒子の系におけるガウス型状態準備（$T=0$）において、この確率的手法は座標および運動量の分散に関する厳密な解析式を再現する。同様に、調和振動子とオーミック・バスの平衡状態において、本手法は量子的な振る舞いの特徴である相関関数の代数的減衰を正しく捉える。これは、古典的モデルが予測する指数関数的減衰とは対照的である。
• デコヒーレンス動力学： シュレディンガーの猫状態（空間的に分離した波束の重ね合わせ）のシミュレーションにおいて、厳密な量子力学（フル量子ノイズスペクトルを使用）は、高温マスター方程式や古典的なホワイトノイズモデルの予測と比較して、著しく速い初期デコヒーレンスを示す。これは、古典的な近似で無視される真空揺らぎや低温度のバスモードの決定的な役割を浮き彫りにしている。
• 熱化： 本手法は、高周波のバスモード（真空揺らぎ）によって駆動される、短時間スケール（$t \sim \varepsilon$）での運動量分散の急速な初期増加を正しく再現する。これはホワイトノイズ近似には存在しない特徴である。
• 温度独立性： 本アプローチの数値的な複雑さは温度に本質的に依存せず、古典的な高温領域から深い量子的なゼロ温度限界まで、シームレスな運用を可能にする。

非二次ポテンシャルへの拡張
滑らかな非二次ポテンシャルに対して、著者らは自然な微小パラメータとして、平衡コヒーレンス長 $\lambda = \hbar / \sqrt{\langle p^2 \rangle_{eq}}$ を特定している。

• バスのスペクトル遮断周波数が増加するにつれ、$\lambda$ は縮小し、簡約密度行列は座標表示において強く準対角的になる。
• 潜在的なポテンシャルが $\lambda$ のスケールに対して緩やかに変化する場合（すなわち、$\epsilon = \lambda/L \ll 1$、ここで $L$ はポテンシャルの特性長さスケール）、確率的アプローチは誤差 $O(\epsilon)$ の制御された近似を提供する。これは、本手法が厳密に解ける二次モデルを超えて拡張可能であることを示唆している。

意義と主張
本論文は、QBM問題に対する「深遠な簡略化」を提供すると主張している。カルデイラ・レゲット・モデルの厳密な量子力学を古典的な確率過程へとマッピングすることで、無限の自由度を持つ開いた量子系のシミュレーションを、扱えるレベルの古典的軌跡の生成と平均化というタスクへと還元している。

主な貢献は以下の通りである：

• 厳密性： 相関のある物理的に正しい熱平衡状態から出発する場合、二次ポテンシャルに対して、この対応関係は任意の温度において厳密であることが証明されている。
• 汎用性： 本手法は、任意の初期量子状態（非古典的な重ね合わせを含む）および一連の介入を、ウィグナー重み付き軌跡アンサンブルを通じて扱うことができる。
• 効率性： 本手法は、低温における既存の手法の数値的ボトルネックを回避し、複雑な現代的手法に対する概念的に明快な代替手段を提供する、統一された計算ツールを提供する。
• 物理的洞察： 本フレームワークは、量子FDTから導出された量子ノイズ統計が、古典的な近似が逃してしまう零点揺らぎや非マルコフ効果を捉えるためにいかに十分であるかを明示的に示している。

著者らは、このフレームワークが、厳密な量子力学と古典的確率過程の間の溝を埋める、散逸量子動力学を研究するための多用途かつ全温度対応のツールを提供すると結論付けている。