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

技术摘要：用于量子信息诊断的 Bogoliubov 比率框架

问题陈述
描述随时间变化的二模玻色子系统（例如描述原初宇宙扰动或非退化光学参数放大器的系统）的动力学，传统上通过挤压振幅 $r_k(\eta)$ 和挤压相位 $\phi_k(\eta)$，或等效地通过 Bogoliubov 系数 $\alpha_k(\eta)$ 和 $\beta_k(\eta)$ 来表征。虽然这些参数化是完备的，但在处理一般的随时间变化的系统时，直接实现它们往往非常繁琐。变量 $r_k$ 和 $\phi_k$ 满足耦合的非线性微分方程，且相位 $\phi_k$ 可能会发生快速振荡。此外，评估信息论量（如纯度和熵）通常需要针对每个特定的模型或驱动协议来重建并对约化密度矩阵进行对角化。这种“一模型一方法”的处理方式掩盖了归一化二模对态所固有的通用信息论结构，并增加了比较不同物理实现（例如宇宙学与光学）之间的难度。

方法论
作者提出了一个基于 Bogoliubov 比率 的紧凑且统一的框架，该比率定义为 $\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.$$
   这种形式将由广义 Lanczos 系数或特定背景动力学引起的态层变形直接吸收进 $\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$ 的单个封闭 Riccati 型方程：
   $$\frac{d\lambda_k}{d\eta} = g_k(\eta)(1 - \lambda_k^2) - 2i\omega_k(\eta)\lambda_k.$$
   其中，$\omega_k$ 控制自由相位旋转，而 $g_k$ 控制对产生（pair production）。

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}$)、Rényi-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$ 的 Riccati 方程。
  2. 通用诊断： 将 $|\lambda_k|^2$ 转换为信息度量。

结果
该框架在两种不同的非平稳设置中得到了验证：

1. 原初宇宙扰动：

   • 作者将该框架应用于有效频率 $\omega_k = \sqrt{k^2 + \alpha/\eta^2}$ 和耦合 $g_k = -1/\eta$ 的模。
   • 数值解表明，背景诱导的相位旋转和频率软化（由参数 $\nu$ 控制）调节了挤压的增长。
   • 当有效频率软化时（较大的 $\nu$），相位旋转与对产生之间的竞争发生转移，导致更大的 $|\lambda_k|$、降低的单模纯度和增加的熵。
   • 结果阐明了背景动力学如何抑制相比于理想共振基准的相干挤压积累。

2. 啁啾脉冲非退化光学参数放大器：

   • 该框架被应用于具有高斯泵浦包络和随时间变化的频率啁啾（失谐）的系统。
   • 结果捕捉到了相关性产生的延迟启动（局部于脉冲中心附近）以及由于啁啾导致的挤压积累抑制。
   • 与恒定耦合极限下熵无限增长的情况不同，有限的泵浦持续时间和频率啁啾导致了纯度和熵在后期趋于有限的饱和值。
   • 对恒定共振耦合下标准 $\tanh(\tau)$ 极限的数值恢复，验证了所提议的 Riccati 方程的正负号约定和归一化。

意义与主张
论文声称，这种 Bogoliubov 比率形式化提供了一种高效、标准化的诊断工具，适用于广泛的参数驱动型二次玻色子系统。其主要意义在于：

• 统一性： 它将二模对产生中的通用信息内容与导致演化的特定物理机制（如宇宙膨胀 vs 光学泵浦）分离开来。
• 计算效率： 通过避免为每个模型重建约化密度矩阵，它简化了量子信息度量的计算。
• 物理洞察： 它提供了一个清晰的机制，用以理解背景诱导的相位旋转和频率软化如何调节态的混合度（mixedness）和挤压增长。

作者指出，该框架适用于任何演化态仍保持由厄米二次哈密顿量生成的归一化配对 Bogoliubov 表示的系统。他们明确说明，向真正的开放系统动力学（如 Lindblad 演化）或弱非高斯相互作用的扩展超出了当前工作的范围。