2025-11-28T18:52:19.570674

Opening Krylov space to access all-time dynamics via dynamical symmetries

Loizeau, Buča, Sels
Solving short and long time dynamics of closed quantum many-body systems is one of the main challenges of both atomic and condensed matter physics. For locally interacting closed systems, the dynamics of local observables can always be expanded into (pseudolocal) eigenmodes of the Liouvillian, so called dynamical symmetries. They come in two classes - transient operators, which decay in time and perpetual operators, which either oscillate forever or stay the same (conservation laws). These operators provide a full characterization of the dynamics of the system. Deriving these operators, apart from a very limited class of models, has not been possible. Here, we present a method to numerically and analytically derive some of these dynamical symmetries in infinite closed systems by introducing a naturally emergent open boundary condition on the Krylov chain. This boundary condition defines a partitioning of the Krylov space into system and environment degrees of freedom, where non-local operators make up an effective bath for the local operators. We demonstrate the practicality of the method on some numerical examples and derive analytical results in two idealized cases. Our approach lets us directly relate the operator growth hypothesis to thermalization and exponential decay of observables in chaotic systems.
academic

Opening Krylov space to access all-time dynamics via dynamical symmetries

基本信息

  • 论文ID: 2503.07403
  • 标题: Opening Krylov space to access all-time dynamics via dynamical symmetries
  • 作者: Nicolas Loizeau, Berislav Buča, Dries Sels
  • 机构: Niels Bohr Institute (Copenhagen), Université Paris-Saclay, University of Oxford, New York University, Flatiron Institute
  • 分类: quant-ph (量子物理)
  • 发表时间: 2025年8月22日 (v2版本)
  • 论文链接: https://arxiv.org/abs/2503.07403v2

摘要

本文解决了闭合量子多体系统短时和长时动力学的核心难题。对于局域相互作用的闭合系统,局域可观测量的动力学可以展开为Liouvillian的(伪局域)本征模,即所谓的动力学对称性。这些算符分为两类:随时间衰减的瞬态算符和永久振荡或保持不变的永恒算符(守恒律)。本文提出了一种通过在Krylov链上引入自然涌现的开边界条件来数值和解析地推导这些动力学对称性的方法,将Krylov空间划分为系统和环境自由度,其中非局域算符构成局域算符的有效浴。

研究背景与动机

核心问题

  1. 量子多体动力学的计算挑战: 求解闭合量子多体系统的短时和长时动力学是原子物理和凝聚态物理的主要挑战之一
  2. 动力学对称性的提取困难: 虽然理论上局域可观测量的动力学可以完全用Liouvillian的本征模(动力学对称性)来表征,但除了极少数模型外,推导这些算符一直不可能

问题重要性

  • 完整动力学刻画: 动力学对称性提供了系统动力学的完整表征,对于任意算符O和初态ρ,其时间演化可表示为: Tr(O(t)ρ) = ∑_ω e^(iωt)μ_ω
  • 非平衡现象理解: 对理解多体局域化、时间晶体、量子疤痕、碎片化等非平凡动力学行为至关重要
  • 热化机制: 在混沌系统中关联算符增长假设与热化及可观测量指数衰减

现有方法局限

  1. 有限系统限制: 有限闭合系统的Liouvillian谱是纯实的,无法探测热力学极限性质
  2. 伪模式展开方法: 之前的方法(Refs 96,97)局限于接近平衡态的动力学
  3. 伪局域性未研究: 以往研究未考察动力学对称性的(伪)局域性,而这正是影响非平衡物理的关键

研究动机

提出一种更严格的方法,通过在Krylov链上引入开边界条件实现系统-环境分解,使得:

  • 局域量在系统中支撑
  • 非局域量构成环境
  • 能够识别哪些动力学对称性是伪局域的并因此相关

核心贡献

  1. 新颖的开边界条件方法: 提出在Krylov链特定位置L引入耗散边界条件,避免波函数的非物理反弹,形式为: ∂t φ_L = (b_L + b{L+1})φ_ - 2b_{L+1}φ_L
  2. 系统-环境划分框架: 建立Krylov空间的自然分解,其中:
    • 左侧(系统):简单的局域算符,跳跃系数精确已知
    • 右侧(环境):复杂的非局域算符,构成有效浴
  3. 理想情况的解析解: 对两种理想化情况给出精确解:
    • 线性增长 b_n = n: 本征向量用Meixner多项式表示
    • 平方根增长 b_n = √n: 本征向量用Hermite多项式表示
  4. 算符增长与热化的直接联系: 证明在线性增长情况下(混沌系统),所有可观测量以速率-2α指数衰减,其中α是Lanczos系数增长率
  5. 数值验证与应用: 在XXZ链和横场Ising链等实际自旋系统上验证方法有效性,成功提取动力学对称性

方法详解

任务定义

输入:

  • Hamiltonian H
  • 初始算符 O_0 (或初态密度矩阵ρ)

输出:

  • 动力学对称性 A_ω 满足 H, A_ω = -ωA_ω
  • 频率ω(可为复数):实部对应振荡,虚部对应衰减

约束:

  • 只关注伪局域动力学对称性(影响局域物理)
  • 适用于局域相互作用系统和短程关联初态

模型架构

1. Krylov空间构造 (Lanczos算法)

从种子算符O_0开始,递归应用Liouvillian L = H, ·构造正交归一基{O_n}:

O_1 = LO_0/b_1 = [H,O_0]/b_1,  b_1 = ||LO_0||

对于 n > 2:
O'_n = LO_{n-1} - b_{n-1}O_{n-2}
O_n = O'_n/b_n
b_n = ||O'_n||

其中范数定义为 ||O||² = 1/(2N) TrO²

2. 一维单粒子映射

时间演化在Krylov基中映射为1D单粒子跳跃问题:

O(t) = 1/(2N) ∑_n i^n φ_n(t)O_n

系数演化方程: ∂t φ_n = b_n φ - b_{n+1} φ_{n+1}, φ_n(0) = δ_

3. 开边界条件

在位置L截断链,假设φ_n足够光滑可局域近似为线性函数: φ_{L+1} ≈ φ_L + (φ_L - φ_)

得到耗散边界条件: ∂t φ_L = (b_L + b{L+1})φ_ - 2b_{L+1}φ_L

4. 修正的Liouvillian矩阵

截断的Liouvillian具有三对角+边界修正形式:

L̃ = i × [三对角矩阵,最后一行为特殊形式]
最后一行: [0, 0, ..., b_L + b_{L+1}, -2b_{L+1}]

该矩阵非厄米,允许复本征值

技术创新点

1. 边界条件的物理意义

  • 避免参数引入: 不同于Refs 96,97的外加耗散+外推,本方法无需额外参数
  • 光滑性假设: 基于Krylov波函数导数连续性的自然假设
  • 局域-非局域分离: Krylov链右侧算符主要由k>L的k-局域Pauli串主导

2. 与传统方法的区别

  • Dirichlet边界: φ_{L+1} = 0导致非物理反弹
  • 本方法: 允许流入环境并返回,更真实地反映算符传播

3. 设计合理性

在理想情况(b_n = n或√n)下:

  • 开边界条件下动力学与无限链几乎无法区分(图2)
  • 实部谱几乎不变,仅在虚部引入物理的衰减(图3)
  • 线性增长情况:所有非平凡根满足Im(ω) = -2i (数值验证至n=2^12)

实验设置

模型系统

1. XXZ链 (可积模型)

H = ∑_i (s^x_i s^x_{i+1} + s^y_i s^y_{i+1} + Δs^z_i s^z_{i+1} + hs^z_i)
参数: Δ = -1/2, h = 2
初始算符: O_0 = Q_3 = ∑_i (s^+_i s^+_{i+1}s^+_{i+2} + s^-_i s^-_{i+1}s^-_{i+2})

已知存在动力学对称性,频率ω = 12

2. 混沌链

H = ∑_i (s^x_i s^x_{i+1} - 1.05s^z_i + 0.5s^x_i)
初始算符: O_0 = ∑_i (1.05s^x_i s^x_{i+1} + s^z_i)

选择与H无重叠的算符以排除平凡守恒量

数值实现

  • 工具: Julia包 PauliStrings.jl (利用Pauli串表示高效计算嵌套对易子)
  • 平移不变性: 只存储单胞而非广延算符
  • Pauli截断增强: 在每步Lanczos算法中保留M个最大权重的Pauli串(算法1)

对比方法

  1. 矩阵乘积算符(MPO)模拟: 作为基准
  2. 标准Dirichlet边界: 展示非物理反弹
  3. Teretenkov等方法: 线性外推+人工耗散(补充材料比较)

评价指标

  • 动力学对称性频率精度: 与已知结果(XXZ)或解析预测(混沌)比较
  • 时间演化保真度: TrO_0(0)O_0(t) vs 精确/MPO结果
  • 局域性度量: 本征向量在Krylov基中的分布 |ψ_0|²(越局域越暗)

实验结果

主要结果

1. 理想情况验证 (图2-3)

线性增长 b_n = n:

  • 开边界条件下动力学与无限链无法区分
  • 所有本征值满足Im(ω) = -2i (除两个平凡根ω = -i, -3i)
  • Dirichlet边界导致t≈5时显著非物理反弹

平方根增长 b_n = √n:

  • 开边界条件完美重现无限链动力学
  • 谱的实部几乎不变,虚部引入小的衰减

2. XXZ链结果 (图4a,c)

  • 动力学对称性恢复: 在ω = 12附近成功提取已知动力学对称性
    • L=26: Re(ω) ≈ 12.0, Im(ω) ≈ -0.5
    • Pauli截断增强(L=50): 收敛性显著改善
  • 时间演化: TrO_0(0)O_0(t)展示持续振荡(与MPO吻合至t≈3)
  • 局域性: 仅实值动力学对称性是局域的(|ψ_0|²集中在小n)

3. 混沌链结果 (图4b,d)

  • 无永恒对称性: 所有模式均衰减,无纯实频率
  • 衰减率聚集: 本征值集中在Im(ω) ≈ -0.72i附近
  • 与理论预测一致: -0.72 ≈ -2λ,其中λ是Lanczos系数增长率
  • 指数衰减: TrO_0(0)O_0(t)快速衰减至0(与MPO吻合)
  • 非局域性: 所有模式极度非局域(|ψ_0|²均匀分布)

4. 稳定态淬火 (图5)

从|ψ_0⟩ = 2^{-N/2}(|0⟩+|1⟩)^⊗N淬火:

  • Q_1 = ∑_i (s^+_i + s^-_i): 快速衰减
  • Q_5 (3阶5-局域动力学对称性展开): 长时振荡
  • 不同L(24-30)结果收敛良好
  • 标准截断(灰线)出现非物理行为

消融实验

Pauli截断的影响 (补充图4)

  • 无Pauli截断: L_max = 26 (精确计算极限)
  • 有Pauli截断 (保留2^22个串): L_max = 50
    • 动力学对称性频率: Re(ω)从≈12.1改善到≈12.05
    • Im(ω)从≈-0.8改善到≈-0.2

收敛性分析 (补充图2-3)

  • XXZ模型: Re(ω)在L=10后稳定,Im(ω)持续改善
  • 混沌链: 谱在L=30后基本收敛,深色(大L)本征值密集在Im(ω)=-2λ线附近

与Teretenkov方法比较 (补充图6-7)

混沌链:

  • Teretenkov方法(γ=0.2-0.4)的本征值不收敛到解析结果
  • 本方法准确再现Im(ω)=-2λ的聚集

XXZ链:

  • 两种方法对局域动力学对称性(ω=12)精度相似
  • 非局域模式的预测显著不同

实验发现

  1. 算符增长与热化的定量关系: 在线性增长(混沌)情况,衰减率精确为-2α,直接验证了理论预测
  2. 局域性判据: Krylov表示中|ψ_0|²自动提供动力学对称性局域性的量化指标
  3. 方法适用性:
    • 混沌系统(线性增长):效果优异,光滑性假设良好
    • 可积系统(平方根增长):可提取动力学对称性,但收敛较慢
  4. 计算瓶颈: 主要限制是计算高阶Lanczos系数(算符变得非局域),Pauli截断提供有效缓解

相关工作

Krylov空间方法

  • 算符复杂度: Parker等105用于探测量子混沌中的算符增长
  • 流体力学: Refs 97, 117-119研究输运性质
  • Floquet系统: Refs 120-124应用于周期驱动系统

伪模式展开

  • Mori投影: 早期计算时间相关关联函数的方法92
  • 近期发展: Refs 96, 97建议添加耗散+外推或研究Ruelle-Pollicott共振98-104
  • 局限: 限于接近平衡态,未考察伪局域性

动力学对称性理论

  • 永恒平衡: Buča 82的理论框架,状态形式ρ(t) = Z^{-1}exp(∑_u μ_u e^{iλ_u t}A_u)
  • 谱生成代数: Refs 33, 34, 90, 91相关但不必伪局域
  • 时间晶体与疤痕: Refs 12, 29-44中动力学对称性的应用

本文优势

  1. 无参数: 不引入人工耗散率等额外参数
  2. 不需外推: 直接使用已知Lanczos系数
  3. 伪局域性: 明确识别相关的伪局域对称性
  4. 远离平衡: 适用于强淬火等非平衡情况

结论与讨论

主要结论

  1. 方法有效性: 开边界条件成功实现无限闭合系统Liouvillian谱的数值计算
  2. 理论洞察:
    • 混沌系统:线性Lanczos增长→所有可观测量以-2α速率衰减
    • 可积系统:平方根增长→存在永恒振荡模式
  3. 实用价值: 结合Pauli串方法,提供计算多体动力学的强大工具

局限性

  1. 光滑性假设: 方法依赖于Krylov表示中φ_n的光滑性
    • 混沌系统:线性增长保证光滑性
    • 其他情况:需要验证
  2. 计算复杂度:
    • 高阶算符O_n变得极度非局域
    • XXZ模型仅能精确计算至L=26
    • Pauli截断引入近似但有效
  3. 长时精度:
    • 动力学对称性虚部未完全收敛导致长时误差
    • 需要更大L或改进截断策略
  4. 热态淬火困难: 从热态ρ = e^{-βO_0}淬火需要计算非对易算符和的指数,数值上困难(仅能处理稳定态)

未来方向

  1. 量子电路: 推广到Liouvillian为上Hessenberg矩阵的情况107
  2. 其他伪局域对称性:
    • 半局域动力学对称性137-139
    • 量子疤痕的伪局域对称性30, 31, 36, 82
  3. NMR实验模拟: 直接计算Tr(Z_tot(t)Z_tot(0))的Fourier变换140
  4. 改进收敛:
    • 更精细的截断策略
    • 自适应Pauli截断
    • 递归算法(补充材料算法)优化

深度评价

优点

1. 方法创新性 ⭐⭐⭐⭐⭐

  • 原创边界条件: 基于波函数光滑性的自然物理假设,不同于人工添加耗散
  • 系统-环境分解: 在算符空间层面实现局域-非局域分离的巧妙框架
  • 理论-数值结合: 理想情况解析解(Meixner/Hermite多项式)与实际系统数值验证相辅相成

2. 理论深度 ⭐⭐⭐⭐⭐

  • 算符增长-热化联系: 首次定量证明线性增长→-2α衰减的普适关系
  • 伪局域性量化: 通过|ψ_0|²自动识别相关动力学对称性
  • 永恒平衡框架: 与Buča的理论82深度整合

3. 实验充分性 ⭐⭐⭐⭐

  • 多模型验证: 可积(XXZ)和混沌链两个典型案例
  • 多重对比: MPO、Dirichlet边界、Teretenkov方法
  • 收敛性分析: 详细的L依赖性和Pauli截断研究
  • 不足: 缺少更多实际物理系统(如Hubbard模型)的测试

4. 技术实现 ⭐⭐⭐⭐⭐

  • 开源代码: 提供Julia实现(GitHub + PauliStrings.jl)
  • 算法优化: 平移不变性、Pauli截断、递归算法等实用技巧
  • 可复现性: 详细的补充材料和参数设置

不足

1. 光滑性假设的严格性

  • 缺乏理论证明: 何时φ_n足够光滑可线性近似?
  • 失效情况: 强碎片化或特殊边界条件下可能失效
  • 建议: 需要更严格的数学论证或适用性判据

2. 计算可扩展性

  • L上限: XXZ仅达26(50需截断),限制高精度应用
  • 系统尺寸: 未讨论N(系统大小)的依赖性
  • 2D推广: 仅1D链,高维系统困难

3. 实验验证有限

  • 淬火类型: 仅稳定态淬火,未实现热态淬火
  • 时间尺度: 长时精度受限(t<5对XXZ)
  • 物理量: 主要是算符自关联,缺少其他可观测量

4. 与Teretenkov方法比较不充分

  • 仅补充材料: 主文未详细讨论
  • 参数敏感性: 未系统研究γ和外推方案的影响
  • 优劣分析: 何时本方法更优?何时等价?

影响力

对领域的贡献 ⭐⭐⭐⭐⭐

  • 范式转变: 将开放系统思想引入闭合系统动力学研究
  • 普适工具: 适用于广泛的量子多体问题
  • 理论基础: 为理解热化、时间晶体、疤痕等提供统一框架

实用价值 ⭐⭐⭐⭐

  • 数值方法: 可直接用于计算实际系统动力学
  • 实验指导: NMR等实验的理论支持
  • 局限: 计算成本仍限制超大系统应用

可复现性 ⭐⭐⭐⭐⭐

  • 代码开源: GitHub完整实现
  • 详细文档: 补充材料包含所有算法细节
  • 参数透明: 所有实验设置明确

适用场景

理想场景

  1. 1D局域相互作用系统: 自旋链、晶格玻色子等
  2. 混沌系统: Lanczos线性增长保证方法有效性
  3. 短中期动力学: t<10τ(特征时间尺度)
  4. 局域可观测量: 测量k-局域算符

需谨慎场景

  1. 强碎片化系统: 光滑性假设可能失效
  2. 长程相互作用: Krylov空间结构可能不同
  3. 极长时动力学: 虚部收敛不足导致误差累积
  4. 高维系统: 计算成本急剧增加

不适用场景

  1. 开放系统: 已有Lindbladian框架更合适
  2. 强测量系统: 需要不同的理论工具
  3. 经典极限: 过度复杂,经典方法更简单

技术细节补充

关键公式推导

边界条件的主方程形式 (补充材料S20)

精确时间演化Ȯ = iH,O,硬边界条件对应:

Ȯ = iH,O - 1/(2N) Tr(O_L† iH,O)O_L

代入开边界条件:

Ȯ = iH,O - 1/(2N) Tr(O_L† iH,O)O_L + ib_{L+1}Tr(O_† O)O_L - 2b_{L+1}Tr(O_L† O)O_L

展示了边界条件如何平衡流入环境的算符流

Meixner多项式解 (补充材料S3-S7)

线性增长b_n = n时,本征向量:

φ_n(ω) = i^n M_n(ω)/n!

其中M_n满足递推:M_{n+1}(x) = xM_n(x) - n²M_(x)

边界条件给出频率约束: (ω + 2i(L+1))M_L(ω) = (2L+1)LM_(ω)

平凡根ω = -i, -3i对应φ_n ∝ 1和φ_n ∝ 2n-1

数值技巧

递归算法 (补充材料算法)

改进动力学对称性提取:

  1. 随机初始化k-局域O_0
  2. 对角化L̃,选择频率≈ω的本征模A_ω
  3. 截断A_ω至k-局域,更新O_0 ← truncate(A_ω)
  4. 重复至收敛

XXZ模型中k=3,4步收敛Re(ω)从12.14→12.02

参考文献 (精选)

方法论基础

  • 82 B. Buča, PRX 13, 031013 (2023) - 永恒平衡理论
  • 105 D. Parker et al., PRX 9, 041017 (2019) - Krylov复杂度
  • 96,97 A. Teretenkov, O. Lychkovskiy et al., PRB (2024) - 伪模式展开

应用领域

  • 12 B. Buča et al., Nat. Commun. 10, 1730 (2019) - 时间晶体
  • 29-44 量子疤痕系列文献
  • 119 N. Loizeau et al., SciPost (2025) - PauliStrings.jl

理论背景

  • 1 L. D'Alessio et al., Adv. Phys. 65, 239 (2016) - 热化综述
  • 86-89 伪局域性理论(Prosen, Doyon等)

总体评价: 这是一篇在量子多体动力学领域具有重要创新的优秀论文(评分4.5/5)。核心方法(开边界Krylov链)简洁优雅且物理图像清晰,理论推导严谨(理想情况解析解),数值验证充分。主要贡献在于建立了算符增长与热化的定量联系,并提供了实用的数值工具。局限性主要在于光滑性假设的适用范围需进一步明确,以及计算可扩展性。对于研究量子多体系统非平衡动力学、时间晶体、量子疤痕等前沿问题的研究者,这是一篇必读文献。