2025-11-16T19:13:13.606657

Krylov Complexity Meets Confinement

Jiang, Halimeh, Srivatsa

In high-energy physics, confinement denotes the tendency of fundamental particles to remain bound together, preventing their observation as free, isolated entities. Interestingly, analogous confinement behavior emerges in certain condensed matter systems, for instance, in the Ising model with both transverse and longitudinal fields, where domain walls become confined into meson-like bound states as a result of a longitudinal field-induced linear potential. In this work, we employ the Ising model to demonstrate that Krylov state complexity--a measure quantifying the spread of quantum information under the repeated action of the Hamiltonian on a quantum state--serves as a sensitive and quantitative probe of confinement. We show that confinement manifests as a pronounced suppression of Krylov complexity growth following quenches within the ferromagnetic phase in the presence of a longitudinal field, reflecting slow correlation dynamics. In contrast, while quenches within the paramagnetic phase exhibit enhanced complexity with increasing longitudinal field, reflecting the absence of confinement, those crossing the critical point to the ferromagnetic phase reveal a distinct regime characterized by orders-of-magnitude larger complexity and display trends of weak confinement. Notably, in the confining regime, the complexity oscillates at frequencies corresponding to the meson masses, with its power-spectrum peaks closely matching the semiclassical predictions.

academic

Krylov Complexity Meets Confinement

基本信息

论文ID: 2511.03783
标题: Krylov Complexity Meets Confinement
作者: Xuhao Jiang, Jad C. Halimeh, N. S. Srivatsa
机构: Ludwig Maximilian University of Munich, Technical University of Munich, Max Planck Institute of Quantum Optics, Munich Center for Quantum Science and Technology, Kyung Hee University
分类: cond-mat.stat-mech, cond-mat.str-el, hep-lat, hep-th, quant-ph
发表时间: 2025年11月5日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2511.03783

摘要

本文研究了量子多体系统中禁闭（confinement）现象与Krylov复杂度之间的关系。作者使用具有横向场和纵向场的Ising模型，证明了Krylov态复杂度可以作为禁闭现象的敏感且定量的探针。研究发现：在铁磁相内进行量子淬火时，纵向场的引入会显著抑制Krylov复杂度的增长，反映了禁闭导致的关联动力学变慢；在顺磁相内淬火时，复杂度随纵向场增加而增强，表明不存在禁闭；跨越临界点的淬火表现出数量级更大的复杂度和弱禁闭趋势。重要的是，在禁闭区域，复杂度的振荡频率对应于介子质量，其功率谱峰值与半经典预测精确匹配。

研究背景与动机

要解决的核心问题

禁闭现象的量子信息刻画：禁闭是高能物理中的基本现象（如夸克禁闭），也出现在某些凝聚态系统中。如何用量子信息理论的工具定量刻画禁闔动力学是一个开放问题。
Krylov复杂度作为探针：虽然Krylov算符复杂度已被提出用于研究大N规范理论中的禁闭-去禁闭相变，但Krylov态复杂度在量子淬火后如何反映实时禁闭特征仍未被系统研究。

问题的重要性

跨学科意义：连接了高能物理、凝聚态物理和量子信息理论三个领域
实验相关性：Ising模型中的禁闭现象已在准一维化合物和冷原子系统中被观测到，甚至在IBM量子计算机上得到验证
理论价值：为理解量子信息传播如何受到禁闭约束提供新视角

现有方法的局限性

传统探针的限制：纠缠熵、关联函数等传统量度虽能反映禁闭，但不够直接和定量
算符依赖性：传统观测量的功率谱可能因矩阵元消失而遗漏某些能量差
缺乏统一框架：需要一个既能定量刻画信息传播又能编码束缚态谱的统一工具

研究动机

作者受到以下观察启发：禁闭对量子多体系统中信息和关联的传播施加严格约束。因此，量子复杂度这一量化量子信息扩散的量度可能对禁闭现象特别敏感。

核心贡献

建立Krylov复杂度作为禁闭探针：首次系统展示Krylov态复杂度可以作为禁闭现象的敏感且定量的探测工具
揭示不同淬火区域的复杂度特征：
- 铁磁相内淬火：纵向场导致复杂度显著抑制
- 顺磁相内淬火：复杂度随纵向场增加而增强
- 跨临界点淬火：复杂度数量级增大，显示弱禁闭趋势
Krylov谱学方法：证明复杂度的功率谱峰值精确对应介子束缚态质量，与半经典预测高度一致
理论优势分析：阐明Krylov复杂度相比传统算符期望值的优势——它是态特定的、算符独立的，能扫描所有可及的频率成分

方法详解

任务定义

输入：

初态 $|\Psi_0\rangle$ （全极化铁磁态或顺磁态）
淬火后的哈密顿量参数 $(h_x, h_z)$

输出：

Krylov态复杂度 $C_k(t)$ 的时间演化
复杂度的功率谱 $S_k(\omega)$

目标：通过分析 $C_k(t)$ 的增长行为和振荡特征识别禁闭信号

模型架构

1. 物理模型：横纵场Ising模型

哈密顿量为： $\hat{H} = -J\sum_{j=1}^{N}\left[\hat{\sigma}^z_j\hat{\sigma}^z_{j+1} + h_x\hat{\sigma}^x_j + h_z\hat{\sigma}^z_j\right]$

其中 $J=1$ （能量单位）， $h_x$ 是横向场， $h_z$ 是纵向场。

物理图像：

$h_z=0$ $h_{z} = 0$ ：模型可积，存在量子相变于 $h_x=1$ $h_{x} = 1$
- $h_x<1$ ：铁磁相，畴壁自由传播
- $h_x>1$ ：顺磁相
$h_z\neq 0$ ：破坏可积性和 $Z_2$ 对称性，产生线性禁闭势 $\chi = 2Jh_z\bar{\sigma}$ ，其中 $\bar{\sigma}=(1-h_x^2)^{1/8}$ 是自发磁化

2. Krylov复杂度构造

Krylov子空间生成：通过Lanczos算法构造正交归一基 $\{|K_n\rangle\}$ ： $\mathcal{K}_n(\hat{H}, |\Psi_0\rangle) = \text{span}\{|\Psi_0\rangle, \hat{H}|\Psi_0\rangle, \ldots, \hat{H}^{n-1}|\Psi_0\rangle\}$

三对角化： $\hat{H}|K_n\rangle = \alpha_n|K_n\rangle + \beta_{n+1}|K_{n+1}\rangle + \beta_n|K_{n-1}\rangle$

其中 $\alpha_n$ （对角）和 $\beta_n$ （非对角）是Lanczos系数。

时间演化：将演化态展开为 $|\Psi(t)\rangle = \sum_n \psi_n(t)|K_n\rangle$ ，满足： $i\partial_t\psi_n(t) = \alpha_n\psi_n(t) + \beta_n\psi_{n-1}(t) + \beta_{n+1}\psi_{n+1}(t)$

这等价于半无限链上的紧束缚模型。

复杂度定义： $C_k(t) = \sum_n n|\psi_n(t)|^2$

物理意义：量化波函数在Krylov空间中的扩散程度，类似于粒子在链上的平均位置。

3. 功率谱分析

$S_k(\omega) = \left|\int_{-\infty}^{\infty} e^{i\omega t}C_k(t)dt\right|^2$

通过离散傅里叶变换（DFT）计算，时间步长 $\Delta t=0.1$ ，总时间 $t=100$ 。

技术创新点

1. 状态特定的探针

与传统算符期望值 $\langle\hat{O}(t)\rangle$ 不同，Krylov复杂度仅依赖于初态和哈密顿量，避免了矩阵元 $\langle E_b|\hat{O}|E_a\rangle$ 消失导致的信息丢失。

2. 几何解释

将复杂度动力学映射到半无限链上的量子粒子运动，提供了信息扩散的几何图像。

3. 全频谱扫描

复杂度动力学扫描从初态可达的所有频率成分，提供完整的能级结构信息。

4. Full Orthogonalization Lanczos (FOL)

采用FOL算法确保数值稳定性，每个新向量对所有已生成基向量正交化，避免精度损失。终止条件： $\beta_n < 10^{-12}$ 。

实验设置

系统参数

系统尺寸：主要结果使用 $L=14$ 个格点
边界条件：周期边界条件
能量单位： $J=1$

淬火协议

协议1：铁磁相内淬火

初态：全极化铁磁态（所有自旋向上）
淬火目标： $h_x=0.25$ （铁磁相内）
纵向场范围： $h_z \in [0, 0.4]$

协议2：顺磁相内淬火

初态：顺磁态（ $h_x=2$ ）
淬火目标： $h_x=1.75$ （顺磁相内）
纵向场范围： $h_z \in [0, 0.4]$

协议3：跨临界点淬火

初态：顺磁态（ $h_x=2$ ）
淬火目标： $h_x=0.25$ （铁磁相内，跨越 $h_x=1$ 临界点）
纵向场范围： $h_z \in [0, 0.4]$

介子质量计算

半经典方法：将两个畴壁视为处于线性禁闭势中的费米子，有效两体哈密顿量： $H_{\text{two-body}} = \epsilon(p_1) + \epsilon(p_2) + \chi|x_2-x_1|$

其中单粒子色散关系： $\epsilon(p) = 2J\sqrt{1+h_x^2-2h_x\cos p}$

通过Bohr-Sommerfeld量子化条件求解束缚态能量： $2E_n(P)p_a - \int_{-p_a}^{p_a}dp\,\omega(p;P) = 2\pi\chi\left(n-\frac{1}{4}\right)$

其中 $\omega(p;P) = \epsilon(p+P/2) + \epsilon(p-P/2)$ 。

具体参数：

$h_x=0.25, h_z=0.2$ ：两个介子 $m_1=4.025J, m_2=4.702J$
$h_x=0.25, h_z=0.1$ ：四个介子 $m_1=3.662, m_2=4.127, m_3=4.480, m_4=4.769$

数值实现细节

Lanczos算法：FOL变体，重复Gram-Schmidt正交化确保数值稳定性
时间演化：求解三对角哈密顿量的薛定谔方程
功率谱：使用DFT，时间窗口 $t\in[0,100]$ ，步长 $\Delta t=0.1$

$h_z=0$ ：复杂度显示大振幅振荡，振幅随系统尺寸增加（见插图a），符合关联的光锥传播
$h_z>0$ ：
- 复杂度增长被显著抑制
- 抑制程度随 $h_z$ 增强
- 振荡频率增加但振幅减小
- 系统尺寸效应消失（插图b），表明关联传播被阻止

物理解释：纵向场使介子获得大有效质量，淬火能量不足以使其运动，只能在静止状态下产生介子。

2. 顺磁相内淬火（图2）

关键发现：

$h_z=0$ ：复杂度被抑制，仅有小振幅涨落（自由费米子行为）
$h_z>0$ ：复杂度振幅逐渐增长

对比结论：与铁磁相完全相反的趋势，表明顺磁相不支持禁闭，纵向场引入的相互作用增强了量子混沌动力学。

3. 跨临界点淬火（图3）

关键发现：

复杂度达到数量级更大的值（比相内淬火大几个数量级）
初始随 $h_z$ 增加，但在更大场强下开始下降
叠加随机涨落

物理解释：

大复杂度：跨临界点激发宽连续谱，在Krylov空间强烈退局域化
转折行为：可能标志弱禁闭的出现，但信号较弱（可能因介子有限传播速度）

Krylov谱学：介子质量提取（图4）

实验参数： $h_x=0.25, h_z=0.2$

功率谱峰位：

高频峰精确对应半经典预测的介子质量
- $m_1 = 4.025J$
- $m_2 = 4.702J$
谱还捕获介子质量间隔 $m_{12} = m_2 - m_1$

更多介子情况（补充材料图S2）： $h_x=0.25, h_z=0.1$ 时，谱成功分辨四个介子质量：

$m_1=3.662, m_2=4.127, m_3=4.480, m_4=4.769$
质量间隔 $m_{12}, m_{13}$ 也被捕获

验证精度：功率谱峰位与半经典计算的偏差在数值精度范围内。

有限尺寸效应分析

铁磁相（ $h_z=0$ ）：

显著的有限尺寸效应（图1插图a）
反映自由畴壁的弹道传播

禁闭区域（ $h_z=0.4$ ）：

有限尺寸效应被抑制（图1插图b）
确认关联传播的缺失

Krylov复杂度是禁闭的敏感探针：
- 铁磁相内淬火：纵向场导致显著抑制
- 顺磁相内淬火：复杂度增强，无禁闭
- 跨临界点淬火：复杂度极大，显示弱禁闭趋势
Krylov谱学方法：
- 功率谱峰精确对应介子质量
- 与半经典预测高度一致
- 提供完整的能级结构信息
理论优势：
- 态特定、算符独立
- 扫描所有可及频率成分
- 无矩阵元消失问题

局限性

系统尺寸限制：
- 主要结果使用 $L=14$
- 更大系统的计算成本高（Krylov空间维度指数增长）
弱禁闭信号：
- 跨临界点淬火中的禁闭信号较弱
- 可能需要更精细的分析
模型特定性：
- 仅研究了Ising模型
- 其他禁闭模型（如规范理论）有待探索
温度效应：
- 仅考虑零温度淬火
- 有限温度算法55有待应用

未来方向

作者明确提出以下开放问题：

格点规范理论：
- Krylov复杂度如何反映规范约束导致的禁闭39-48？
- 能否捕获弦断裂51,52等非微扰现象？
其他非微扰现象：
- 假真空衰变49,50
- 量子多体疤痕53,54
算符增长：
- 研究禁闭如何影响Krylov空间中的算符增长
- 使用有限温度算法探测接近基态的现象
实验实现：
- 冷原子系统中的直接测量
- 量子模拟器上的验证
理论扩展：
- 高维系统
- 长程相互作用
- 拓扑效应

跨学科桥梁：巧妙连接量子信息（Krylov复杂度）与高能/凝聚态物理（禁闭）
新探针开发：首次将Krylov态复杂度系统应用于禁闭研究
理论深度：提供了复杂度-禁闭关系的几何和物理直觉

2. 方法严谨性（★★★★★）

多协议验证：三种不同淬火协议全面测试复杂度行为
半经典对照：介子质量计算提供独立验证
数值稳定性：FOL算法确保高精度
系统尺寸检验：通过有限尺寸效应分析确认物理图像

3. 结果说服力（★★★★★）

定量一致性：功率谱峰与半经典预测的精确匹配（图4）
定性对比鲜明：铁磁相与顺磁相的相反趋势清晰展示禁闭/非禁闭
数量级差异：跨临界点淬火的极大复杂度提供强信号

4. 物理洞察（★★★★★）

统一理解：将复杂度抑制与关联传播受限直接关联
谱学优势：阐明Krylov方法相比传统算符的本质优势
机制解释：介子有效质量增大导致静态产生的图像清晰

5. 写作质量（★★★★☆）

逻辑清晰：从动机到方法到结果层次分明
图表优质：图1-4直观展示核心发现
补充材料：详细算法和额外数据增强可复现性

不足

1. 系统尺寸限制（★★★☆☆）

计算瓶颈： $L=14$ 相对较小，可能存在有限尺寸效应
标度分析缺失：未系统研究复杂度如何随 $L$ 标度
热力学极限：无法确认 $L\to\infty$ 时的行为

2. 跨临界点分析不充分（★★★☆☆）

弱禁闭信号：仅观察到"趋势"，缺乏定量刻画
随机涨落：未深入分析涨落的物理起源
机制不明：为何复杂度先增后减的详细机制未充分讨论

3. 模型覆盖有限（★★★☆☆）

单一模型：仅研究Ising模型
规范理论缺失：虽提出作为未来方向，但本文未涉及
普适性未知：方法是否适用于其他禁闭系统需验证

4. 实验连接薄弱（★★☆☆☆）

测量方案缺失：未讨论如何在实验中测量Krylov复杂度
噪声效应：未考虑实验噪声对复杂度的影响
可行性分析：冷原子或量子模拟器实现的具体方案缺乏

5. 理论深度可提升（★★★☆☆）

解析结果有限：主要依赖数值，缺乏解析理论
因果关系：复杂度抑制与禁闭的因果链条可更严格建立
临界行为：未研究接近临界点 $h_x\to 1$ 时的标度行为

影响力评估

对领域的贡献（★★★★★）

方法学贡献：
- 为禁闭研究提供新工具
- 扩展Krylov复杂度的应用范围
- 建立量子信息-高能物理新联系
理论意义：
- 深化对量子信息传播受限的理解
- 为非微扰现象提供量子信息视角
启发性：
- 开启多个后续研究方向（规范理论、疤痕态等）
- 可能激发实验组开发新测量方案

实用价值（★★★★☆）

量子模拟：
- 为冷原子实验设计提供理论指导
- 量子计算机上的禁闭验证
诊断工具：
- 可用于识别未知系统中的禁闭
- 介子谱提取的新方法
局限：
- 当前仅限于小系统数值模拟
- 实验实现路径不明确

可复现性（★★★★☆）

优点：

模型简单明确
Lanczos算法标准且详细描述（补充材料）
参数完整给出

改进空间：

代码未公开（虽然算法清晰）
某些数值细节（如时间演化方法）可更详细

适用场景

理想应用场景

低维量子自旋链：
- 禁闭势明确的系统
- 可精确数值模拟的尺寸
冷原子量子模拟：
- 可控淬火协议
- 长时间相干演化
量子计算机验证：
- 小规模禁闭模型
- 算法基准测试

限制场景

大系统/高维：计算成本指数增长
强耗散系统：需要开放系统推广24
有限温度：需要新算法55
复杂规范理论：方法有效性待验证

技术细节亮点

1. Full Orthogonalization Lanczos

创新点：采用FOL而非标准Lanczos
优势：避免数值正交性丢失
实现：重复Gram-Schmidt，终止条件 $\beta_n<10^{-12}$

2. 功率谱计算

参数选择： $t\in[0,100], \Delta t=0.1$ 平衡分辨率与计算成本
峰识别：清晰分辨介子质量和质量间隔

3. 半经典理论

Bohr-Sommerfeld量子化：提供独立理论预测
数值求解：转化为寻找转折点的非线性方程

参考文献（精选）

核心相关文献

16 Kormos et al., Nat. Phys. 13, 246 (2017)：Ising模型中禁闭的开创性研究
18 Rabinovici et al., arXiv:2507.06286：Krylov复杂度综述
33 Anegawa et al., JHEP 2024, 119：大N规范理论中的算符复杂度
38 Rutkevich, J. Stat. Phys. 131, 917 (2008)：介子谱的半经典理论

实验相关

17 Vovrosh & Knolle, Sci. Rep. 11, 11577 (2021)：IBM量子计算机上的禁闭
10-12：准一维化合物中的实验观测
47 Zhang et al., Nat. Phys. (2024)：冷原子中的拓扑θ角调控

方法学

19,20：Krylov复杂度的几何解释
35,36：Lanczos算法经典文献

总体评分

维度	评分	说明
创新性	9/10	首次系统建立Krylov复杂度-禁闭联系
严谨性	8/10	方法严谨，但系统尺寸和模型覆盖有限
影响力	9/10	跨学科贡献，开启新研究方向
实用性	7/10	理论价值高，实验路径待明确
写作质量	8/10	清晰流畅，图表优质
综合评价	8.2/10	优秀工作，推荐发表于高水平期刊

对研究者的建议

后续工作：尽快扩展到规范理论验证普适性
实验合作：与冷原子组合作开发测量方案
理论深化：发展解析理论理解复杂度-禁闭的深层联系
代码开源：提高可复现性和影响力