Statistics of non-conserved observables in Lindblad master equations

技術的要約：リンドブラッド・マスター方程式における非保存観測量の統計量

問題提起
閉じた量子系において、保存則は観測量 $O$ とハミルトニアン $H$ の交換関係（$[H, O] = 0$）によって厳密に制御されており、これにより初期状態によらず期待値 $\langle O \rangle$ が一定に保たれる。しかし、マルコフ的な環境を介したリンドブラッド・マスター方程式で記述される開いた量子系においては、この条件だけでは不十分である。散逸項はハミルトニアンによる保存則を破り、かつて保存されていた物理量を変化させる可能性がある。期待値の保存の喪失については十分に理解されているが、全統計分布、特に二次のモーメント $\langle O^2 \rangle$ および分散 $\text{Var}(O)$ への影響については、より微細な分析が必要となる。本論文は、ハミルトニアンによって保存されるものの、全開放系ダイナミクスによっては保存されない観測量の統計に対して、崩壊演算子（collapse operators）がどのように影響を与えるかを扱う。

手法
著者は、観測量 $O$ の一次および二次のモーメントの時間発展に関する基本式を導出するために、随伴リンドブラッド形式（adjoint Lindblad formalism）を用いる。

1. 随伴生成子（Adjoint Generator）: 密度行列 $\rho$ に対する標準的なリンドブラッド・マスター方程式から、期待値の時間微分は随伴生成子 $\mathcal{L}^\dagger$ を通じて次のように表される：
   $$\frac{d}{dt}\langle O \rangle = \langle \mathcal{L}^\dagger(O) \rangle$$
   ここで、$\mathcal{L}^\dagger(O) = i[H, O] + \mathcal{D}^\dagger(O)$ であり、$\mathcal{D}^\dagger$ は散逸による寄与を表す。
2. モーメント方程式: この形式を二次のモーメント $\langle O^2 \rangle$ および分散 $\text{Var}(O) = \langle O^2 \rangle - \langle O \rangle^2$ へと拡張する。論文では、分散生成率 $\Pi_O(t)$ を以下のように定義して導出する：
   $$\Pi_O(t) = \langle \mathcal{L}^\dagger(O^2) \rangle - 2\langle O \rangle \langle \mathcal{L}^\dagger(O) \rangle$$
   ハミルトニアンによって保存される観測量（$[H, O]=0$）の場合、この生成率は散逸部分 $\mathcal{D}^\dagger$ のみによって生成される。
3. 解析モデル: 導出された公式を4つの特定の解析モデルに適用し、異なるダイナカルな挙動を示す。
   • 一量子ビット振幅減衰モデル
   • 局所的な崩壊を伴う二量子ビット交換モデル
   • 運動量拡散モデル
   • 異なる結合形式を持つジェーンズ・カミングス・モデル

主要な貢献と結果

• 平均値の保存と分散の保存の区別: 本論文は、期待値の保存が観測量の全統計の保存を意味しないことを確立している。以下の3つの異なるケースを特定している：

  1. 完全な保存: $\mathcal{L}^\dagger(O) = 0$ かつ $\mathcal{L}^\dagger(O^2) = 0$。平均値と分散の両方が一定である。
  2. 平均のドリフト: $\mathcal{L}^\dagger(O) \neq 0$。平均値が変化し、通常は分散の変化を伴う。
  3. 平均保存を伴う分散生成: $\mathcal{L}^\dagger(O) = 0$ だが $\mathcal{L}^\dagger(O^2) \neq 0$。平均は一定に保たれるが、ゆらぎが増大する。

• 解析例:

  • 一量子ビット振幅減衰: $O=\sigma_z$ かつ崩壊演算子が $L=\sigma_-$ の場合、平均 $\langle \sigma_z \rangle$ は基底状態へと減衰する。$\sigma_z^2 = I$ であるため、二次のモーメントは自明であり、分散は平均の減衰のみによって進化する。
  • 二量子ビット交換モデル: 全励起数 $N$ を保存するハミルトニアンに対し、局所的な崩壊演算子が $L_j = \sigma_-^{(j)}$ である場合、平均 $\langle N \rangle$ は減衰する。極めて重要な点として、システムが $N$ の固有状態から始まったとしても（分散がゼロ）、散逸ダイナミクスによって、システムが励起セクター間（$|11\rangle \to |10\rangle/|01\rangle \to |00\rangle$）を遷移することに伴い、動的に非ゼロの分散が生成される。
  • 運動量拡散: ハミルトニアンが $H=p^2/2m$ で、位置に結合した散逸子 $L \propto x$ を持つ粒子のモデルでは、平均運動量は保存される（$\mathcal{L}^\dagger(p)=0$）が、分散は線形に増大する（$\mathcal{L}^\dagger(p^2) \neq 0$）。これは、一次モーメントが保持される一方で分布が広がるケースを明示的に示している。
  • ジェーンズ・カミングス・モデル: 論文では、全励起数 $N$ に関する2つのシナリオを対比させている。
    • デフェージング結合: レザーバーが $N$ 自体を通じて結合する場合、$N$ と $N^2$ はともに保存される（$\mathcal{L}^\dagger(N)=\mathcal{L}^\dagger(N^2)=0$）。この場合、レザーバーは $N$ セクター間のコヒーレンスを破壊するが、$N$ の統計量自体は変化させない。
    • 散逸結合: 標準的な光子損失（$L_a=a$）および原子崩壊（$L_\sigma=\sigma_-$）は、$N$ の平均および二次のモーメントの両方の保存を破る。

• ゆらぎ散逸関係との関連: 分散生成項 $\Pi_O(t)$ は、散逸によるドリフトとゆらぎ生成を結びつける自然な対象として特定される。熱モデルにおいて、この形式は標準的なゆらぎ散逸バランスを回収する。すなわち、減衰を引き起こすのと同じ環境結合が、平衡状態における分散をも決定するのである。

意義と主張
本論文は、開放系におけるハミルトニアン保存観測量の過渡的な統計に対する「単純かつ操作的な」診断法を提供すると主張している。その主な意義は、平均のドリフトとは独立して、ゆらぎを生成する役割としての随伴生成子の散逸部分を孤立させている点にある。

著者は、この研究が不変演算子や定常部分空間の新しい分類を行うものではなく、むしろリンドブラッド対称性の厳密な理論を補完するものであり、解析的に透明性の高いモデルに焦点を当てたものであることを強調している。さらに、リンドブラッド源項と理想化された波動関数の収縮モデルとの間の類似性を指摘している。これによれば、収縮モデル（あるいはリンドブラッド演算子によるもの）において、局所的な保存則はイベントごとに破られる可能性があるが、平均化すれば統計的には保持される可能性がある。この文脈において、二次のモーメントと分散は、平均値が保存されているように見える場合であっても、潜在的な確率論的な非保存イベントに関する情報を運ぶ担い手となるのである。