Recursive perturbation approach to time-convolutionless master equations: Explicit construction of generalized Lindblad generators for arbitrary open systems

技术摘要：时间无卷积（TCL）主方程的递归摄动方法

问题陈述
开放量子系统的理论需要稳健的框架，以描述在现实环境中由于系统-环境相互作用导致退相干和不可逆行为的动力学过程。虽然 Lindblad 主方程被广泛使用，但它依赖于马尔可夫动力学、弱耦合以及清晰的时间尺度分离假设，这些假设在具有强耦合或结构化谱密度的机制中往往失效。在这些非马尔可夫机制中，时间无卷积（Time-Convolutionless, TCL）主方程通过利用考虑环境反馈作用的时间局部生成元，提供了更优越的描述。然而，TCL 生成元的解析推导通常局限于可解模型（例如高斯环境或线性相互作用）。对于涉及非高斯库或非线性耦合的复杂场景，由于嵌套时间积分和环境相关函数在更高阶时会出现复杂度指数级增长，传统的摄动方法变得难以处理。此外，虽然 TCL 生成元可以写成类 Lindblad 形式，但将其分解为相干（哈密顿量）和耗散部分在本质上是非唯一的，这导致在定义有效系统哈密顿量和耗散率时存在歧义。

方法论
作者为以左作用和右作用系统算符显式表达的 TCL 生成元开发了一种递归摄动展开法。该方法假设一个初始不相关状态（$\rho_{SE}(0) = \rho_S \otimes \rho_E$）的一般二分系统，其总哈密顿量为 $H_t = H_{S,t} + H_{E,t} + \lambda A_t \otimes B_t$。

1. 递归展开： 该方法建立在文献 [16] 中先前开发的递归摄动展开基础之上。动力学映射按耦合强度 $\lambda$ 进行展开。生成元 $L_t$ 通过将动力学映射的时间导数与其逆映射进行复合来构建（$L_t = \dot{\Phi}_t \circ \Phi_t^{-1}$）。
2. 算符表示： 展开式使用超算符 $A^L[\rho] = A\rho$ 和 $A^R[\rho] = \rho A$ 进行表达。这种表示法允许系统算符与环境算符的系统性分离。展开式的系数被定义为“广义累积量”（$D$），它们整合了多时环境相关函数和时间排序。
3. 递归累积量： 一个关键的技术创新是这些广义累积量的递归定义（公式 27）。这使得高阶项可以由低阶项构建，从而显著减少了计算冗余，并实现了摄动项的图示化表示。
4. 规范分解： 作者将得到的生成元映射到广义 Lindblad 形式。为了解决分离相干部分和耗散部分固有的歧义，他们引入了最小耗散原则。这是通过要求 Lindblad 算符是无迹的（具体为 $V_i \to V_i - \langle V_i \rangle_{1/d} \mathbb{1}$）来实现的。这一约束唯一地指定了分解方式，确保了耗散部分在希尔伯特-施密特平均范数下是最小的，并产生了一个定义良好的有效哈密顿量。

主要贡献与结果

• 广义 Lindblad 形式： 本文提供了一个在规范 Lindblad 类结构下，在任意阶次下的 TCL 生成元的系统性表达式：
  $$L[X] = -i[K, X] + \sum_{ij} \gamma_{ij} \left( L_i X L_j^\dagger - \frac{1}{2} \{L_j^\dagger L_i, X\} \right)$$
  此处，速率 $\gamma_{ij}$ 和算符 $L_i$ 均随时间变化。负速率的出现被识别为非 CP-可分性的特征，表明了非马尔可夫动力学。
• 有效哈密顿量构建： 作者推导出了有效系统哈密顿量 $K(t)$ 的递归公式（公式 45）。该项捕捉了由于环境相互作用导致的系统能级重整化。该公式区分了偶数阶和奇数阶摄动，表明哈密顿量的偶数阶贡献源自累积量的虚部，而奇数阶贡献则源自实部。
• 显式计算： 为了验证该方法，作者显式计算了直到四阶的生成元：
  • 一阶： 产生一个正比于 $\langle B_t \rangle$ 的有效驱动项。
  • 二阶： 恢复了涉及两点浴相关函数的标准二阶 TCL 结果。
  • 三阶和四阶： 作者提供了系数及由此产生的哈密顿量贡献的显式表达式。他们证明了在假设一阶项消失（例如 $\langle B_t \rangle = 0$）的情况下，递归公式（公式 27）相比于直接积分，能显著简化四阶系数的计算。

意义与主张
本文声称，这种递归摄动方法为推导任意开放量子系统的 TCL 生成元提供了一种系统且可处理的方法，将其扩展到了高斯或线性相互作用模型的限制之外。通过将生成元表述为具有无迹算符的规范 Lindblad 形式，该方法解决了定义有效哈密顿量和耗散率时的歧义，即使在强耦合和非马尔可夫机制中也能提供具有物理可解释性的分解。

作者指出，该框架能够研究记忆效应显著的复杂环境。直到四阶的显式推导展示了该方法在处理非高斯效应和非线性耦合方面的有效性。这项工作被呈现为研究自旋系统相干动力学的工具，并用于佐证数值证据，提供了一条在保持生成元物理可解释结构的同时，系统研究开放量子系统的路径。本文并非提出新的实验装置，而是提供了一种分析现有模型并解释非马尔可夫现象（如能级分裂重整化和非 CP-可分性）的理论工具。