An analytical solution of a quantum system with non-Markovian behavior: The Bixon-Jortner system in time domain

技术摘要：Bixon-Jortner 系统在时域内的解析解

问题陈述
本文旨在解决获取量子系统非马尔可夫（non-Markovian）动力学精确解析描述的挑战。虽然非马尔可夫行为通常在二分系统（系统加环境）中进行研究，其中对环境求迹会导致记忆效应，但标准方法往往依赖于近似，如 Lindblad 主方程（马尔可夫极限）或微扰修正（例如 Nakajima-Zwanzig 形式体系）。作者专注于 Bixon-Jortner (BJ) 模型，该模型由一个离散态 $|\phi\rangle$ 与一个无限“阶梯”态 $|n\rangle$ 耦合组成。与标准的系统-环境模型不同，BJ 的希尔伯特空间是直和而非张量积。尽管 BJ 模型已知会表现出指数衰减（Wigner-Weisskopf 极限）及随后的衰减-复现（decay-revival）动力学，但以往的解析处理依赖于拉普拉斯变换，这掩盖了在整数时间间隔观察到的“踢”（kicks，即非平滑过渡）的物理起源。本文的目标是在不借助近似的情况下，推导出精确的时域动力学，并明确展示该系统的非马尔可夫本质。

方法论
作者通过以下步骤直接在时域内处理时间演化问题：

1. 公式化： 从相互作用绘景下的含时薛定谔方程出发，推导出离散态系数 $b_\phi(t)$ 的积分-微分方程。该方程通过对 $b_\phi$ 的过去历史进行积分，消除了阶梯态的系数。
2. 转化为时滞微分方程 (DDE)： 通过针对 BJ 模型参数（均匀阶梯间距 $\Delta$ 且 $E_\phi=0$）进行特化，并应用泊松求和公式，将积分-微分方程转化为时滞微分方程。在这种形式下，$b_\phi$ 在时间 $\tau$ 的时间导数取决于 $b_\phi$ 在先前时刻 $\tau - m$（其中 $m$ 为整数）的值，使得非马尔可夫记忆效应变得显而易见。
3. 通过假设进行解析求解： 为了求解所得的 DDE，作者采用了加法分解假设。他们定义了一个变换函数 $f(\tau) = b_\phi(\tau)e^{\beta\tau}$，并假设解可以写成在每个整数时间 $\tau = n$ 时“介入”的加法项 $f_n(\tau-n)$ 之和。
4. 递推关系与生成函数： 通过分析特定时间区间 $(k, k+1)$ 内的方程，他们推导出了加法分量的递推关系。该关系通过生成函数方法进行求解，将解与拉盖尔多项式 $L_n(x)$ 联系起来。

主要结果

• 精确的非马尔可夫方程： 本文推导出了离散态振幅的精确时滞微分方程（方程 19）。该方程明确显示，系统的演化受马尔可夫指数衰减项和一系列时滞项的控制。这些时滞项在每个整数时间单位 $\tau = n$ 引入“踢”，导致动力学偏离简单的指数衰减并引起部分复现。
• 闭式解： 作者得到了概率振幅 $b_\phi(\tau)$ 的精确解析解（方程 56）。该解表示为涉及拉盖尔多项式的求和形式：
  $$b_\phi(\tau) = e^{-\beta\tau} \left( 1 + \sum_{n=1}^{\infty} H(\tau-n) e^{\beta n} [L_n(2\beta(\tau-n)) - L_{n-1}(2\beta(\tau-n))] \right)$$
  其中 $H$ 是 Heaviside 阶跃函数，$\beta$ 是无量纲耦合参数。
• 衰减与复现动力学： 解证实了对于 $0 \le \tau < 1$，系统经历纯指数衰减，其速率为 $\Gamma = 2\pi v^2 / \Delta$（与费米黄金定则一致）。在 $\tau = 1, 2, 3, \dots$ 时，时滞项被激活，修改了衰减速率并导致概率表现出非平滑过渡和部分复现，这在针对不同 $\beta$ 值的绘图结果中得到了可视化验证。

意义与主张
作者声称，其主要贡献在于提供了一个显式、易懂且精确的量子系统非马尔可夫性的案例研究。

• 透明度： 不同于以往产生解析公式但掩盖了“踢”之物理机制的拉普拉斯变换法，这种时域时滞微分方程方法使系统的非马尔可夫性质变得完全透明。记忆效应直接体现为对整数间隔处过去值的依赖。
• 易读性： 所使用的求解技术（从积分-微分方程到 DDE 的转化、简单的假设以及生成函数）呈现得直观且易于理解，为非马尔可夫动力学如何产生并能被解析求解提供了一个清晰的教学范例。
• 模型效用： 作者将 BJ 系统定位为一个理解衰减与复现现象的价值极高的“纯净”模型，它剥离了更复杂的轻-物质相互作用中的数学复杂性，同时保留了核心物理特性。

文章结论较为谦逊，指出虽然针对这一特定的直和系统推导了精确解，但该方法有可能推广到其他系统，例如具有非恒定耦合的阶梯系统或传统的张量积二分系统，以便与标准主方程形式体系进行比较。