Closed-Form Expectation Values of the Damped Kerr Oscillator

技术摘要：阻尼克尔振子的闭式期望值

问题陈述
克尔非线性振子由哈密顿量 $H = \omega a^\dagger a + \chi (a^\dagger a)^2$ 描述，是典型的非线性量子光学和超导电路量子电动力学模型。虽然无阻尼系统由于特定的代数性质对于相干初始态是精确可解的，但通过 Lindblad 主方程引入的耗散（光子损耗）通常会导致无法获得期望值的闭式解。标准方法依赖于在截断的 Fock 空间内进行数值积分，这种方法随着希尔伯特空间维度的增加，计算成本会迅速上升，并引入截断误差。此外，尽管阻尼克尔振子在 Liouville 空间中存在精确传播子，但它们并不能直接给出实验相关的正规序期望值 $\langle (a^\dagger)^\ell a^k \rangle$ 的闭式表达式。

方法论
作者通过利用系统 Liouvillian 的底层李代数结构，推导出了一个解析解。该方法通过以下步骤进行：

1. 向量化表示： 将密度矩阵 $\rho$ 映射为双希尔伯特空间（热场表示）中的向量 $|\rho\rangle$。这使得 Lindblad 主方程可以被视为一个线性演化方程 $d|\rho\rangle/dt = \mathcal{L}|\rho\rangle$，其中 $\mathcal{L}$ 作为一个算符作用。
2. 李代数解离： 将 Liouvillian $\mathcal{L}$ 表示为 $su(1,1)$代数生成元$K_+$, $K_-$, $K_3$ 以及一个中心算符 $K_0$ 的组合。作者利用 Wei–Norman 分解法，将时间演化算符 $e^{\mathcal{L}t}$ 解离为一个指数乘积形式：$e^{\Gamma_-(t)K_-}e^{\Gamma_3(t)K_3}$。
3. 精确传播子： 通过求解关于时间相关系数 $\Gamma_-(t)$ 和 $\Gamma_3(t)$ 的微分方程，推导出了作用于 Liouville Fock 基底上的精确传播子。该传播子具有普适性，适用于任意初始密度算符。
4. 向相干态特化： 推导过程将针对初始相干态 $|\alpha\rangle$ 进行特化。相干态在 Fock 基底中的泊松分布结构使得在计算期望值时，剩余的无穷级数可以通过组合恒等式进行解析求和。

主要贡献与结果
本文的主要结果是得到了初始化为相干态的阻尼克尔振子任意正规序矩 $\langle (a^\dagger)^\ell a^k \rangle$ 的闭式解析表达式。该期望值由下式给出：

$$\langle (a^\dagger)^\ell a^k \rangle = \alpha^\ell (\alpha^*)^k e^{-|\alpha|^2} e^{\left(\frac{\kappa}{2} - i(\omega-\chi)(k-\ell)\right)t} e^{\frac{\lambda t}{2}(k+\ell+1)} \exp\left[ |\alpha|^2 e^{\lambda t} \left( 1 - \frac{\kappa}{\lambda}(e^{-\lambda t} - 1) \right) \right]$$

其中 $\lambda = -2i\chi(k-\ell) - \kappa$。

该结果的关键特性包括：

• 普适性： 该公式适用于任何非负整数 $k$ 和 $\ell$，涵盖了场正交分量（例如 $\langle X(t) \rangle, \langle P(t) \rangle$）及高阶矩。
• 高效性： 该解仅依赖于系统参数 ($\omega, \chi, \kappa$) 和初始振幅 $\alpha$。它不需要数值积分，也避免了与截断 Fock 空间相关的误差。
• 验证： 利用 QuTiP 对 Lindblad 主方程进行的直接数值模拟验证了这些解析表达式，结果在各种非线性相互作用强度下均显示出一致性。

意义与主张
本文声称，阻尼克尔振子对于相干态的可解性并非源于矩闭合（这在非线性系统中通常会失效），而是源于 Liouvillian 的 $su(1,1)$ 结构与相干态的泊松结构之间的兼容性。这使得无穷阶耦合矩可以被精确重求和。

作者将这一闭式解定位为计算高阶矩时的一种高效替代方案，因为在截断 Fock 空间中，计算高阶矩的成本会变得极高。此外，论文指出，底层的李代数结构可能为获取更广泛类别的开放量子光学系统的精确解提供出发点。该工作范围较为适中，严格专注于这些特定期望值的推导与验证，并未提出超出这些解析公式直接效用之外的新实验装置或应用。