2025-11-24T22:19:16.851165

Universal Growth of Krylov Complexity Across A Quantum Phase Transition

Grabarits, del Campo

We study the statistical properties of the spread complexity in the Krylov space of quantum systems driven across a quantum phase transition. Using the diabatic Magnus expansion, we map the evolution to an effective one-dimensional hopping model. For the transverse field Ising model, we establish an exact link between the growth of complexity and the Kibble-Zurek defect scaling: all cumulants of complexity exhibit the same power-law scaling as the defect density, with coefficients identical to the mean, and the full distribution asymptotically becomes Gaussian. These results yield general scaling arguments for the growth of complexity across arbitrary second-order quantum phase transitions.

academic

Universal Growth of Krylov Complexity Across A Quantum Phase Transition

基本信息

论文ID: 2510.13947
标题: Universal Growth of Krylov Complexity Across A Quantum Phase Transition
作者: András Grabarits, Adolfo del Campo
分类: quant-ph
发表时间: October 17, 2025
论文链接: https://arxiv.org/abs/2510.13947

摘要

本文研究了量子系统在量子相变过程中Krylov空间中扩散复杂度的统计性质。通过使用非绝热Magnus展开，将演化映射为有效的一维跳跃模型。对于横场伊辛模型，建立了复杂度增长与Kibble-Zurek缺陷标度之间的精确联系：复杂度的所有累积量都表现出与缺陷密度相同的幂律标度，系数与均值相同，完整分布渐近趋于高斯分布。这些结果为任意二阶量子相变的复杂度增长提供了一般的标度论证。

研究背景与动机

要解决的问题：
- 量子多体动力学中复杂度增长是热化、混沌和描述量子演化计算困难度的基础问题
- 时间依赖设置下算符增长和Krylov复杂度的表征仍然基本未开发
- 需要确定复杂度等物理量是否也表现出非平衡普适性
问题的重要性：
- Krylov子空间方法已成为量子化复杂度增长的通用框架
- 在量子控制和变分量子算法中具有实用价值
- 连接复杂度增长与临界量子动力学的普适特征
现有方法的局限性：
- 现有时间依赖设置的一般框架使用非局域时间的Floquet算符，在多体系统中应用困难
- 缺乏连接复杂度增长与量子相变普适性的理论框架
研究动机：
- 开发适用于时间依赖设置的精确且通用的分析框架
- 建立Krylov复杂度与Kibble-Zurek机制的联系
- 为临界点附近的复杂度增长提供普适框架

核心贡献

提出了基于非绝热Magnus算符的Krylov子空间新方法，适用于驱动量子系统的时间依赖设置
建立了复杂度增长与Kibble-Zurek缺陷标度的精确联系，证明了横场伊辛模型中所有复杂度累积量都遵循相同的幂律标度
证明了Krylov复杂度统计在KZ标度区间内收敛到普适高斯分布，完整分布表现出普适行为
提供了适用于任意二阶量子相变的一般标度论证，将结果推广到更广泛的临界系统类别

方法详解

任务定义

研究量子系统在有限速率驱动穿越量子相变时Krylov复杂度的统计性质，建立复杂度增长与相变动力学的普适联系。

模型架构

1. 非绝热Magnus算符构造

对于驱动系统 $H(t) = \sum_n E_n(t)|n(t)\rangle\langle n(t)|$ ，定义非绝热时间演化算符： $U(t) = U(t)U_{ad}(t)^\dagger$

其中 $U_{ad} = \sum_n |n(t)\rangle\langle n(0)|$ 是平行输运算符。非绝热Magnus算符定义为： $\Omega(t) = i\log(U(t))$

时间演化态可表示为： $|\psi(t)\rangle = e^{-i\Omega(t)}|GS(t)\rangle = \sum_{m=0}^{\infty} \frac{(-i)^m \Omega^m(t)}{m!}|GS(t)\rangle$

2. Krylov基构造

相应的Lanczos算法为： $|K_{n+1,t}\rangle b_{n+1,t} = \Omega(t)|K_{n,t}\rangle - a_{n,t}|K_{n,t}\rangle - b_{n,t}|K_{n-1,t}\rangle$

其中：

$a_{n,t} = \langle K_{n,t}|\Omega(t)|K_{n,t}\rangle$
$b_{n,t} = \langle K_{n-1,t}|\Omega(t)|K_{n,t}\rangle$

3. 横场伊辛模型分析

TFIM哈密顿量： $\hat{H}(t) = -J\sum_{j=1}^L [\sigma_j^z \sigma_{j+1}^z + g(t)\sigma_j^x]$

在费米子表示下，每个动量模式对应独立的二能级系统： $\hat{H} = 2\sum_k \hat{\psi}_k^\dagger H_k(t) \hat{\psi}_k$

其中 $H_k(t) = (g(t) - \cos k)\tau_z + \sin k \tau_x$ 。

技术创新点

时间局域性：不同时间的Krylov基和Lanczos系数独立生成，避免了时间非局域性问题
精确解析结果：在KZ标度区间内获得了Lanczos系数和Krylov波函数的精确表达式： $b_n \sim L^{1/2}\tau^{-1/4}\sqrt{n}, \quad a_n \sim L\tau^{-1/2}$
泊松统计：Krylov波函数在领头阶近似下表现为泊松分布： $|\phi_n(\tau)|^2 \sim \frac{(2CL\tau^{-1/2})^n e^{-2CL\tau^{-1/2}}}{n!}$

实验设置

数据集

横场伊辛模型，系统大小 $L = 200, 400$
有效系统大小 $L_{eff} = 46, 48$ （用于Lanczos算法）
驱动时间范围： $\tau = 2$ 到 $\tau = 750$

评价指标

Krylov复杂度的前三个累积量： $K_1$ （均值）、 $K_2$ （方差）、 $K_3$ （三阶中心矩）
复杂度分布的高斯拟合质量
与KZ标度 $\tau^{-1/2}$ 的符合程度

实现细节

线性驱动协议： $g(t) = t/\tau$
初始化在顺磁相 $t = -g_0\tau$
终止于 $g(0) = 0$
计算复杂度限制精确数值模拟至 $L \approx 48$

实验结果

主要结果

累积量的普适标度：所有复杂度累积量在KZ标度区间内都遵循相同的幂律： $K_q \approx 2CL\tau^{-1/2}$
这与缺陷密度的KZ标度 $n \sim \tau^{-1/2}$ 完全一致。
高斯分布收敛：在KZ标度区间 $L\tau^{-1/2} \gg 1$ 内，复杂度分布收敛到高斯分布： $P_K(n) \approx \frac{\exp[-(n-K_1)^2/(2K_2)]}{\sqrt{2\pi K_2}}$
Lanczos系数标度：数值结果验证了理论预测的标度行为，所有不同 $\tau$ 和 $L$ 的数据在适当重标度后坍塌到单一曲线。