来源论文: https://arxiv.org/abs/2605.05106v1 生成时间: May 07, 2026 00:08

执行摘要

在非平衡态量子统计物理领域，理解系统穿越量子临界点（Quantum Critical Point, QCP）时的动力学行为是一项核心挑战。传统的Kibble-Zurek（KZ）机制为预测缺陷密度（如磁畴壁或kink的密度）提供了普适的缩放框架。然而，KZ机制在描述高阶关联函数方面存在局限。本文基于Lakshita Jindal与Kavita Jain的最新研究，深入探讨了一维横场Ising模型（TFIM）在代数形式非线性淬火过程中的kink-kink关联特性。研究发现，除了KZ长度尺度外，在线性和亚线性淬火中，系统演化还依赖于一个额外的“去相位长度”（dephasing length）。更重要的是，去相位后的关联函数展现出一种压缩指数（compressed exponential）衰减行为，其衰减指数随淬火指数连续变化，这一发现挑战了传统的指数衰减预期，为量子动力学研究提供了全新的理论维度。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：超越KZ范式

当一个量子系统以有限速度穿越二阶相变点时，由于临界慢化（Critical Slowing Down），系统在临界区附近必然失去绝热性。Kibble-Zurek机制成功预测了缺陷密度$\rho_d$与淬火速率的缩放关系，但在描述空间关联函数（如$K(r, t)$）时，AIA（Adiabatic-Impulse-Adiabatic）近似往往过于简单，忽略了准粒子激发之间的相位干涉效应。本文的核心问题在于：在非线性淬火路径下，kink之间的空间关联是否仍然遵循简单的缩放？是否存在新的长度尺度？其衰减形式是否具有普适性？

1.2 理论基础：1D TFIM 与非线性协议

研究对象为一维横场Ising模型，其哈密顿量定义为：

$$H(t) = -J \sum_{j=1}^{N} \sigma_j^z \sigma_{j+1}^z - g(t) \sum_{j=1}^{N} \sigma_j^x$$

其中 $g(t)$ 是随时间变化的横场。该模型在 $g_c=1$ 处发生从顺磁（PM）到铁磁（FM）的相变。本工作采用代数形式的淬火协议：

$$|g(t) - g_c| \sim |\frac{t_c - t}{t_c}|^\omega$$

其中 $\omega$ 为淬火指数（$\omega=1$ 为线性，$\omega > 1$ 为超线性，$\omega < 1$ 为亚线性）。

1.3 技术难点：非线性的不可积性

在线性淬火情形下，系统可以通过Landau-Zener理论精确求解。然而，对于通用的 $\omega$，相应的动力学方程（由Jordan-Wigner变换得到的费米子系数演化方程）不再具有解析精确解。这要求研究者结合绝热摄动理论（APT）与高精度数值积分来处理复杂的动力学相位。特别是在脉冲区（Impulse Regime）累积的相位 $\Delta \phi_k$ 对最终关联函数的去相位效应至关重要。

1.4 方法细节：Jordan-Wigner 与对偶变换

研究采用了以下技术路径：

  1. 费米子化：通过Jordan-Wigner变换将自旋算符映射为无自旋费米子算符。此时，哈密顿量在动量空间解耦为一系列 $2 \times 2$ 的Bogoliubov子空间。
  2. 对偶映射：为了研究kink（即铁磁相中的磁畴壁），引入对偶算符 $\mu_i^x = \sigma_i^z \sigma_{i+1}^z$。在该变换下，原模型的kink关联函数等效于对偶模型的费米子关联函数。
  3. 系数演化：解Heisenberg方程组 $i\hbar \frac{d}{dt} \binom{u_k}{v_k} = \mathcal{M}(t) \binom{u_k}{v_k}$。初始条件设定为 $g_i > 1$ 处的瞬时基态。
  4. 关联函数分解：将总关联函数 $K(r, t)$ 分解为对角项 $K_{on}$（与激发概率 $p_k$ 相关）和非对角项 $K_{off}$（与量子干涉项 $\langle c_k^\dagger c_{-k}^\dagger \rangle$ 相关）。

2. 关键 Benchmark 体系，计算所得数据与性能分析

2.1 Benchmark 体系设置

研究采用环形链（Ring Boundary Condition），节点数 $N = 10^4$，确保其在热力学极限附近。固定 $J=1, \hbar=1$，淬火时间 $\tau$ 通常设为 $10^3$ 到 $10^4$ 量级，以保证处于慢淬火区。起始场 $g_i$ 从 2 或 10 变化到最终场 $g_f=0$。

2.2 激发概率 $p_k$ 的特征数据

通过绝热摄动理论计算发现，对于大的淬火指数 $\omega$，激发概率 $p_k$ 在动量空间展现出强烈的振荡行为（见论文Fig. 2b）。这是因为在超线性淬火中，能隙在大部分时间内保持恒定，导致基态与激发态之间发生多次转换。数据表明，当 $\omega$ 为奇整数时，$p_k$ 表现为 $X$ 的周期函数，而对于亚线性情形，则逐渐过渡到幂律衰减。

2.3 关键结论一：去相位长度 $\ell$

通过分析动力学相位 $\Delta \phi_k \sim k^2 \ell \hat{\xi}$，定义了去相位长度 $\ell$：

  • $\omega < 1$: $\ell \sim \tau^{\frac{1}{1+\omega}}$
  • $\omega = 1$: $\ell \sim \sqrt{\tau \ln \tau}$
  • $\omega > 1$: $\ell \sim \hat{\xi} \sim \tau^{\frac{\omega}{1+\omega}}$ 计算结果证明，对于 $\omega \le 1$，$\ell$ 的增长速度快于KZ长度 $\hat{\xi}$，导致系统在KZ尺度内发生明显的去相位。而对于超线性淬火，单一尺度 $\hat{\xi}$ 即可描述系统。

2.4 关键结论二：压缩指数衰减

这是本文最引人注目的发现。解析推导显示，$K_{on}(r, \tau)$ 的渐近形式为：

$$|K_{on}(r, \tau)| \sim \exp[-(r/\hat{\xi})^{\omega+1}]$$
  • 当 $\omega=1$ 时，回归到高斯衰减（$\exp(-r^2)$），这与前人对线性淬火的研究一致。
  • 当 $\omega \neq 1$ 时，衰减速度由 $\omega$ 调制。这一非普适性（Non-universal）特征在图4b中得到了数值验证，数值点完美落在理论预测的解析线上。

3. 代码实现细节，复现指南与工具包链接

3.1 核心算法：Bogoliubov-de Gennes (BdG) 演化

复现该工作的关键在于精确解 $N/2$ 个独立的 $2 \times 2$ 厄米矩阵方程。建议采用以下步骤：

  1. 动量离散：由于系统具有平移对称性，$k = \frac{(2m+1)\pi}{N}$，其中 $m=0, \dots, N/2 - 1$。
  2. 数值积分器：使用 Adaptive Runge-Kutta 4/5（如 Python 中的 scipy.integrate.solve_ivp 或 C++ 中的 Boost.Numeric.Odeint）。需要监控相位的积累，防止在高频振荡区丢失精度。
  3. 初值计算：利用瞬时对角化 $H_k(g_i)$ 得到初始的 $(u_k, v_k)$。

3.2 关联函数计算逻辑

得到演化后的 $u'_k(\tau), v'_k(\tau)$ 后，利用公式：

$$K(r, \tau) = \left| \frac{1}{\pi} \int_0^\pi dk \sin(kr) u'^*_k v'_k \right|^2 - \left| \frac{1}{\pi} \int_0^\pi dk \cos(kr) |u'_k|^2 \right|^2$$

注意：由于被积函数包含 $\sin(kr)$ 和 $\cos(kr)$，在大 $r$ 时需要极细的动量网格，建议采用 Filon 方法进行数值傅里叶变换。

3.3 推荐软件包与开源资源

  • QuSpin (Python): 虽然主要用于精确对角化，但其提供的对称性基组可以辅助构建 TFIM 模型。GitHub Link
  • Dynamical-Ising (Custom C++ Template): 建议自行编写 $2 \times 2$ 矩阵演化代码。核心循环如下（伪代码）：
for(double k : momenta) {
    SystemState state = initial_ground_state(k, g_initial);
    Integrate(state, quench_protocol, time_range);
    u_final[k] = state.u; v_final[k] = state.v;
}

4. 关键引用文献与局限性评论

4.1 关键参考文献

  1. Kibble (1976) & Zurek (1985): KZ机制的奠基之作。
  2. Dziarmaga (2005): 首次给出了1D TFIM线性淬火下缺陷密度的精确解。
  3. Cincio et al. (2007): 讨论了线性淬火下的去相位尺度，是本文对比的基准。
  4. Barankov & Polkovnikov (2008): 绝热摄动理论在量子淬火中的应用。

4.2 局限性评论

尽管本工作揭示了丰富的关联动力学，但仍存在以下局限：

  • 可积性依赖：TFIM 是典型的可积模型。在非可积系统（如增加近邻相互作用 $J_z \sigma_j^z \sigma_{j+1}^z$）中，准粒子间的散射可能会破坏这种基于单粒子相位的去相位图景。
  • 零温近似：实验中（如D-Wave退火机或里德堡原子阵列）不可避免存在热激发，热噪声如何与动力学相位竞争，本文未予讨论。
  • 一维限制：高维情况下的拓扑缺陷关联（如2D中的涡旋关联）可能具有完全不同的拓扑保护机制，其衰减指数可能受空间维度的强烈约束。

5. 补充内容：从理论到实验的跨越

5.1 实验验证路径

本研究提出的关联函数衰减行为可以通过以下实验平台验证：

  • 量子退火机 (D-Wave)：利用其可编程的横场变化，通过多次读取计算自旋关联 $\langle \sigma_i^z \sigma_{j}^z \rangle$。
  • 里德堡原子 (Rydberg Atoms)：利用光镊阵列精确控制原子间距，通过荧光成像技术直接观测 kink 的空间分布。
  • 离子阱 (Trapped Ions)：利用长程相互作用模拟 Ising 动力学，测量多体关联。

5.2 最佳淬火协议的设计

论文在第VIII节讨论了“最优淬火”。对于固定的总时间 $\tau$，存在一个最优的 $\omega$ 使得最终缺陷密度最小（见图9）。这对于绝热量子计算具有直接指导意义：并不是越平滑的淬火越好，适当的非线性加速可能在跨越临界点时更有效地抑制激发。

5.3 物理直觉：为什么是压缩指数？

传统的指数衰减对应于单一特征长度，而压缩指数（$\exp(-r^n), n>1$）通常暗示系统内部存在多尺度相位的相干叠加。在本模型中，$\omega$ 决定了脉冲区的持续时间及其对不同动量模的“激发权重”。随着 $\omega$ 增大，高动量模的激发被抑制得更快，导致关联函数的有效截断变得更加尖锐，反映在数学上即是 $n = \omega + 1$ 的增加。这一发现深化了我们对量子系统“记忆”淬火路径方式的理解。

5.4 总结与展望

这项工作标志着量子淬火研究从“计数”（计算缺陷数量）转向了“成像”（描绘缺陷间的空间排布）。未来的研究方向可能包括：多体局域化（MBL）背景下的非线性淬火、具有长程相互作用的模型、以及非厄米量子系统的临界动力学。对于量子化学工作者而言，这种对动力学相位的精细处理方法，也为研究强关联电子系统中的光致相变提供了有益的数学框架。