A Bogoliubov-ratio framework for quantum-information diagnostics of time-dependent two-mode Boson Hamiltonian

技術要約：量子情報診断のためのボゴリューボフ比フレームワーク

問題提起
原始宇宙論的摂動や非退化光パラメトリック増幅器を記述するような、時間依存する二モード・ボゾン系のダイナミクスは、従来、スクイージング振幅 $r_k(\eta)$ とスクイージング位相 $\phi_k(\eta)$、あるいは同等にボゴリューボフ係数 $\alpha_k(\eta)$ および $\beta_k(\eta)$ によって特徴付けられる。これらのパラメータ化は完全ではあるものの、一般的な時間依存系における直接的な実装はしばしば煩雑である。変数 $r_k$ と $\phi_k$ は結合した非線形微分方程式に従い、また位相 $\phi_k$ は急速な振動を起こすことがある。さらに、情報理論的な量（純粋度やエントロピーなど）を評価するには、特定のモデルや駆動プロトコルごとに簡約密度行列を再構成し、対角化する必要がある。このようなモデルごとのアプローチは、正規化された二モード対状態に内在する普遍的な情報理論的構造を不明瞭にし、異なる物理的実現（例：宇宙論的 vs 光学的）間の比較を困難にする。

手法
著者らは、$\lambda_k(\eta) \equiv \beta_k(\eta)/\alpha_k(\eta)$ と定義されるボゴリューボフ比に基づく、コンパクトで統一されたフレームワークを提案している。これは規格化条件 $|\alpha_k|^2 - |\beta_k|^2 = 1$ および収束領域 $|\lambda_k| < 1$ に従うものである。

1. 状態表現: 著者らは、正規化された一般化二モードスクイーズド状態が、単一の複素変数 $\lambda_k$ を用いて、対となる数基底 $|n, n\rangle$ において正確に表現できることを確立した：
   $$|\psi_k(\eta)\rangle = \sqrt{1 - |\lambda_k(\eta)|^2} \sum_{n=0}^{\infty} \lambda_k^n(\eta) |n_k, n_{-k}\rangle.$$
   この定式化は、状態レベルの変形（一般化ランチョス係数や特定の背景ダイナミクスから生じるもの）を $\lambda_k$ の中に直接吸収している。

2. ダイナミクスの簡約化: 次のような一般的な時間依存二次二モードハミルトニアン：
   $$\hat{H}_k(\eta) = \omega_k(\eta) (\hat{a}^\dagger_k \hat{a}_k + \hat{a}^\dagger_{-k} \hat{a}_{-k} + 1) + i g_k(\eta) (\hat{a}^\dagger_k \hat{a}^\dagger_{-k} - \hat{a}_k \hat{a}_{-k}),$$
   に対して、シュレディンガー方程式は $\lambda_k$ に関する単一の閉じたリッカチ型方程式へと簡約される：
   $$\frac{d\lambda_k}{d\eta} = g_k(\eta)(1 - \lambda_k^2) - 2i\omega_k(\eta)\lambda_k.$$
   ここで、$\omega_k$ は自由位相回転を制御し、$g_k$ は対生成を制御する。

3. 普遍的な診断: パートナーモード（$-k$）をトレースアウトすると、$\lambda_k$ の位相は消失する。一モード簡約密度行列のスペクトルは、完全に $\lambda_k$ の絶対値の二乗 $q_k(\eta) = |\lambda_k(\eta)|^2$ によって決定される：
   $$p_n(\eta) = [1 - q_k(\eta)] q_k^n(\eta).$$
   その結果、標準的な情報理論的指標（純粋度 $\mu_k$、線形エントロピー $S_{L,k}$、レニー2エントロピー $S_{2,k}$、およびフォン・ノイマン・エントロピー $S_{vN,k}$）は、$q_k$ の明示的かつモデルに依存しない関数として導出される。

主な貢献

• 変数の削減: 本フレームワークは、二モードスクイージングのダイナミクスの記述を、二つの結合した実変数（$r_k, \phi_k$）から、単一の複素変数（$\lambda_k$）を制御する一次微分方程式へと削減する。
• モデルに依存しない公式: 著者らは、簡約状態の純粋度およびエントロピーの明示的な公式を、$|\lambda_k|^2$ のみに依存する形で提供している。これにより、新しい時間依存ハミルトニアンごとに簡約密度行列を再構成したり対角化したりする必要がなくなる。
• ダイナミクスと診断の分離: 本フレームワークは、問題を以下の二つのステップに明確に分解する：
  1. モデル依存のダイナミクス: 特定の入力 $\omega_k(\eta)$ および $g_k(\eta)$ を用いて $\lambda_k$ のリッカチ方程式を解く。
  2. 普遍的な診断: $|\lambda_k|^2$ を用いて、固定された普遍的な公式により情報量へと変換する。

結果
本フレームワークは、二つの異なる非定常セットアップにおいて実証されている：

1. 原始宇宙論的摂動:

   • 著者らは、有効周波数 $\omega_k = \sqrt{k^2 + \alpha/\eta^2}$ および結合 $g_k = -1/\eta$ を持つモードに本フレームワークを適用した。
   • 数値解は、背景による位相回転と周波数の軟化（パラメータ $\nu$ によって制御される）が、スクイージングの成長を制御することを示している。
   • 有効周波数が軟化するにつれ（$\nu$ が大きいほど）、位相回転と対生成の間の競争関係が変化し、より大きな $|\lambda_k|$、減少した一モード純粋度、および増加したエントロピーをもたらす。
   • これらの結果は、背景のダイナミクスがいかにして、理想的な共鳴ベンチマークと比較して、コヒーレントなスクイージングの蓄積を抑制するかを明らかにしている。

2. チャープパルス非退化光パラメトリック増幅器:

   • ガウス型ポンプエンベロープと時間依存の周波数チャープ（デチューニング）を持つシステムに本フレームワークを適用した。
   • 結果は、相関生成の遅延した立ち上がり（パルス中心付近に局在）と、チャープによるスクイージング蓄積の抑制を捉えている。
   • エントロピーが無限に増大する定数結合の極限とは異なり、有限のポンプ持続時間と周波数チャープは、有限の遅延時刻における純粋度およびエントロピーの飽和値をもたらす。
   • 定数共鳴結合における標準的な $\tanh(\tau)$ 極限の数値的な回収により、提案されたリッカチ方程式の符号の慣習と規格化が検証された。

意義と主張
本論文は、このボゴリューボフ比の定式化が、広範なクラスのパラメトリックに駆動される二次ボゾン系に対する効率的で標準化された診断ツールを提供することを主張している。その主な意義は以下の通りである：

• 統一: 二モードの対生成に伴う普遍的な情報内容を、その進化を引き起こす具体的な物理的メカニズム（宇宙膨張か光ポンプ駆動か）から分離することにある。
• 計算効率: 各モデルに対して簡約密度行列を再構成することを回避することで、量子情報指標の計算を合理化する。
• 物理的洞察: 背景による位相回転と周波数の軟化が、いかにして状態の混合度（mixedness）とスクイージングの成長を制御するかを理解するための明確なメカニズムを提供する。

著者らは、本フレームワークが、エルミートな二次ハミルトニアンによって生成される正規化された対ボゴリューボフ表現を保持するあらゆる系に適用可能であることを述べている。また、真の開放系ダイナミクス（例：リンドブラッド発展）や弱い非ガウス相互作用への拡張は、現在の研究の範囲外であることを明示している。