源于斯塔克局域化的量子姆潘巴效应

论文信息

标题: Quantum Mpemba effect from Stark localization

作者: Nico Albert, Masudul Haque, Shovan Dutta

发布日期: 2026-09-21

arXiv ID: 2609.24993v1

PDF 链接: 下载 PDF

3 分钟速览

• 研究问题：如何在量子开放系统中设计一种可调、时间尺度高度分离的 Mpemba 效应，使能量更高的初态反而更快冷却到稳态。
• 核心方法：利用斯塔克局域化（Stark localization）构造能量随位置单调变化的 “量子能量景观”，并通过靠近基态的局域非相干跳跃作为冷却机制。
• 关键结果：单粒子局域初态的弛豫时间呈超指数分离，近似为 $\tau_j \approx \frac{1}{\gamma}[(j-2)!]^2(\Delta/J)^{2(j-2)}$；离耗散点越远的低能态反而越慢冷却（见论文公式 (2) 和图 1(b)）。
• 主要局限：机制依赖较强的斯塔克局域化，通常要求 $J/\Delta \lesssim 1$；论文主要为理论推导与数值模拟，未报告新实验数据。
• 适合读者：研究开放量子系统、非平衡统计力学、冷原子量子模拟以及 Mpemba 效应的研究生和科研人员。

背景与研究动机

Mpemba 效应最初指某些经典系统中较热物体反而更快冷却的现象。近年来，这一概念被推广到封闭和开放量子系统。对于经典系统，崎岖势能景观给出了一个直观机制：系统若陷入局部势能最小值，需要先 “爬坡” 越过势垒才能向全局最小值弛豫，因此靠近局部极小点的低能态可能比某些高能态更慢冷却。

但这种景观图像通常难以直接搬到量子系统。原因是量子能量本征态未必局域在势能景观的确定位置上，局域耗散也未必产生沿能谱的简单定向流动。本文的目标是构造一个量子对应版本：让能量本征态在空间上局域，从而恢复 “位置—能量—耗散路径” 的对应关系。

核心方法和技术细节

论文考虑一条长度为 $L$ 的玻色子紧束缚链。第 $1$ 个格点孤立，势能 $V_1=0$；其余格点施加线性倾斜势

$$V_{j\ge 2}=\varepsilon+\Delta(L-j)$$

其中 $\Delta$ 为倾斜斜率，$\varepsilon$ 是最后一个格点与第 $1$ 个格点之间的势能差。哈密顿量为

$$H=\sum_{j=2}^{L-1}J(b_{j+1}^{\dagger}b_j+\text{h.c.})+\sum_{j=1}^{L}V_j b_j^{\dagger}b_j.$$

第 $2$ 个格点到第 $1$ 个格点的相干跳跃不写入 $H$，而由 Lindblad 跳算子 $\Gamma=\sqrt{\gamma}b_1^{\dagger}b_2$ 描述。整个动力学满足 Lindblad 主方程

$$\dot\rho=-i[H,\rho]+\Gamma\rho\Gamma^{\dagger}-\frac12\{\Gamma^{\dagger}\Gamma,\rho\}.$$

该非相干跳跃只发生在靠近基态的位置，因此构成 “局域冷却”：粒子一旦在第 $2$ 个格点被 “检测” 到，就会跳到第 $1$ 个格点，最终进入局域在格点 $1$ 的稳态。

在倾斜势下，哈密顿量本征态经历瓦尼尔—斯塔克局域化。对于 $J/\Delta$ 不太大的情况，本征态局域在某个格点附近，本征能量近似跟随 $V_j$。这意味着势能景观直接映射到能量本征态位置：格点 $2$ 附近能量最高，格点 $L$ 附近能量较低，而格点 $1$ 通过选择足够大的 $\varepsilon$ 成为真正基态。

单粒子弛豫时间由局域本征态在耗散格点上的权重决定。对于 $\gamma\ll\Delta$，论文在 End Matter 中给出

$$\tau_j \approx \frac{1}{\gamma}\frac{1}{|\psi_j(x=2)|^2} \approx \frac{1}{\gamma}[(j-2)!]^2(\Delta/J)^{2(j-2)}.$$

因此，初始在格点 $j$ 的态越靠右，能量越低，但在格点 $2$ 的波函数尾巴越以超指数方式变小，冷却越慢。这个时间尺度分离是参数性的，并且来自量子局域化尾巴的指数小量。对于 $\gamma\gg\Delta$，强耗散把格点 $2$ 积掉，产生等效跳

$$\gamma_{\text{eff}}=\frac{4J^2}{\gamma},$$

弛豫时间近似变为

$$\tau_j \approx \frac{\gamma}{4J^2}\frac{1}{|\psi_j(x=3)|^2} \approx \frac{\gamma}{4J^2}[(j-3)!]^2(\Delta/J)^{2(j-3)}.$$

论文图 1(d) 显示，这两个极限估计与 Lindbladian 谱数值结果吻合；中间 $\gamma$ 区用 Padé 插值也能很好地描述。

Lindbladian 谱的结构也直接反映了这种局域化。论文图 2 显示，除了最快衰减的模式，右本征模在能量本征基中非常局域。实本征值对应布居衰减，复本征值对应相干衰减。慢模衰减率与 $\tau_j$ 对应，形成分层的亚稳态结构。

创新点与贡献

本文的主要贡献是把 “崎岖能量景观” 这一经典 Mpemba 机制量子化。经典景观图像通常依赖势垒控制的热激活逃逸，而本文的机制利用瓦尼尔—斯塔克局域化产生的指数小波函数尾巴，得到比经典扩散强得多的弛豫时间分离。

另一个重要贡献是明确比较了经典对应过程。论文图 4 构造了一个经典 Markov 跳跃模型：粒子在右侧格点间以速率 $C$ 扩散，从格点 $2$ 到格点 $1$ 以速率 $R\gg C$ 冷却。结果显示，经典演化主要是扩散传播，寿命按 $(j-2)^2$ 增长，相邻初态的寿命比趋于 $1$，因此没有明显的 Mpemba 效应。这反衬出量子版本中的超指数时间尺度分离并非简单的经典随机行走效应。

论文还进一步分析了两个相互作用玻色子。引入在位相互作用 $U$ 后，系统出现三种动力学区域：在吸引相互作用区和强排斥相互作用区，两个粒子以 doublon 形式集体跳跃；中间区域则表现为两个粒子分别弛豫。数值结果表明，吸引相互作用可以强烈增强热初态的 Mpemba 信号，而强排斥相互作用则抑制该效应（见论文图 5(c)–(d)）。

数值结果分析

论文采用单粒子局域初态、能量本征态和热 Gibbs 态进行数值检验。图 1(b) 显示，不同初始格点的能量曲线出现明显交叉：初始在格点 $2$ 附近的高能态更快进入基态，而初始在右端低能态的曲线衰减极慢。这种曲线交叉正是 Mpemba 效应的直接表现。

对于热初态

$$\rho^{(T)}\propto e^{-H/T},$$

论文图 3(a) 用迹距离刻画弛豫。迹距离在 Lindblad 动力学下单调下降，因此适合作为 Mpemba 检测量。在该图中，不同温度初态的迹距离曲线出现交叉，说明较热初态在特定参数范围内更快接近稳态。

论文还使用最慢衰减左本征模与热态的重叠

$$c_{\text{SM}}(T)=\mathrm{tr}[l_{\text{SM}}^{\dagger}\rho^{(T)}]$$

作为判据。图 3(b) 显示该重叠随 $T$ 先增后减，最大值约在 $T\sim \mathcal{O}(\varepsilon)$ 附近。这个非单调性来自两个竞争因素：温度升高最初让系统布居更均匀，增加慢模权重；但温度继续升高后，该低能慢模的布居又被稀释。

在两玻色子情形中，图 5(b) 显示慢模本征值在中间区域不随 $U$ 变化，而在集体跳跃区域随 $|U|$ 增大而减小。图 5(d) 显示慢模重叠峰值在吸引相互作用区域最高，在强排斥区域没有明显峰值。这表明相互作用可以调节 Mpemba 效应的显著程度。

局限与待解决问题

该机制的有效性依赖于较强的瓦尼尔—斯塔克局域化。论文提到分析在 $J/\Delta$ 达到 $\mathcal{O}(1)$ 时仍然稳健，但如果局域长度过大，单粒子能量与位置的对应关系会减弱，超指数时间尺度分离可能不再成立。论文主要给出单粒子和两玻色子结果，对更多粒子、更高维或更一般相互作用下的行为没有系统展开。

此外，本文是理论数值工作，没有提供新的实验观测。作者推测在当前冷原子平台中可能实现，尤其是 $J/\Delta=0.5$ 时交叉发生在约 $10^5$ 次 Bloch 振荡内；但这仍属于实验可行性判断，需要通过实际测量检验。退相干、外场不均匀性、高激发态布居等因素在论文中部分提及，但未做完整处理。

还有一个方法上的限制：中间耗散强度区的寿命估计依赖 Padé 插值，而不是严格的微扰展开。虽然图 1(d) 显示该插值与数值吻合良好，但它属于近似描述，并非控制性推导。最后，经典类比说明该效应高度依赖量子相干局域化；如果实验中存在强退相干，可能在某种程度上退化为经典扩散行为，削弱 Mpemba 效应。这些都需要后续工作进一步澄清。