Exact bounds for dynamical critical exponents of transverse-field Ising chains with a correlated disorder

基本信息

• 论文ID: 2007.07439
• 标题: Exact bounds for dynamical critical exponents of transverse-field Ising chains with a correlated disorder
• 作者: Tatsuhiko Shirai (早稻田大学), Shu Tanaka (庆应义塾大学)
• 分类: quant-ph cond-mat.dis-nn
• 发表时间: 2021年4月7日 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2007.07439

摘要

本研究调查了横向场下无序伊辛链的动力学临界指数，以检验相关无序对量子相变的影响。相关无序是指位点上的横向场强度依赖于连接该位点的最近邻耦合强度，这与非相关无序给出了定性不同的结果。在非相关无序情况下（横向场在位点上均匀或与最近邻耦合强度无关地随机选择），动力学临界指数是无穷大的。相比之下，在存在相关无序的情况下，作者解析地证明了动力学临界指数是有限的，并且依赖于横向场强度的调节过程。

研究背景与动机

问题定义

1. 核心问题：研究横向场伊辛链中相关无序对量子相变(QPT)性质的影响，特别是动力学临界指数z的行为。
2. 问题重要性：
   • 量子相变是凝聚态物理和统计物理中的重要现象，发生在绝对零温度下
   • 动力学临界指数z决定了临界点附近特征能量尺度与长度尺度的关系：Δ ∼ ξ^(-z)
   • 理解无序对量子相变的影响对量子绝热计算等应用具有重要意义
3. 现有方法局限性：
   • 在非相关无序情况下，由于Griffiths-McCoy奇异性的存在，动力学临界指数z = ∞
   • 这种发散的z值限制了量子绝热算法的性能
   • 缺乏对相关无序情况的精确解析结果
4. 研究动机：探索通过引入相关无序来抑制z值，使其从无穷大降低到有限值，从而改善量子计算性能。

核心贡献

1. 首次解析证明：在相关无序情况下，动力学临界指数z是有限的，这与非相关无序的z = ∞形成鲜明对比。
2. 精确边界结果：
   • 弱无序情况：z = 1（与无无序情况相同）
   • 强无序情况：max(D(1/2 + |s-1/2|) + 1/2, 1) ≤ z ≤ D + 1
3. 参数依赖性分析：证明了z依赖于调节参数s和无序强度D，提供了调控z值的理论基础。
4. 数值验证：通过精确对角化方法验证了理论预测，特别是在强无序情况下的边界行为。

方法详解

任务定义

研究一维N自旋伊辛链在横向场下的量子相变性质，哈密顿量为：

H = -∑[i=1 to N-1] J_i σ_i^z σ_{i+1}^z - ∑[i=1 to N] Γ_i σ_i^x

其中相关无序通过以下关系定义：

ln(Γ_1/Γ) = (1-s)ln(J_1)
ln(Γ_i/Γ) = s ln(J_{i-1}) + (1-s)ln(J_i)  (2 ≤ i ≤ N-1)
ln(Γ_N/Γ) = s ln(J_{N-1})

模型架构

1. 无序分布模型

• 弱无序：均匀分布π_w(J_i) = (1-J^(0))^(-1) for J^(0) < J_i ≤ 1
• 强无序：幂律分布π_s(J_i) = (1/D)J_i^(-1+1/D) for 0 < J_i ≤ 1

2. 费米化映射

通过Jordan-Wigner变换将自旋系统映射到自由费米子系统：

H = ∑∑ [c_i† [A]_{ij} c_i + [B]_{ij}/2 (c_i† c_j† + c_j c_i)] + const.

3. 本征值问题

能隙问题归结为矩阵本征值问题：

M φ_n = Λ_n² φ_n

其中M = (A-B)(A+B)，能隙Δ = Λ_1。

技术创新点

1. 相关无序设计

与传统非相关无序不同，横向场强度Γ_i依赖于相邻的耦合强度J_和J_i，参数s控制这种依赖关系。

2. 精确边界方法

• 下界：通过矩阵不等式M ≥ M_0和Weyl单调性定理获得能隙下界
• 上界：利用Rayleigh-Ritz变分原理和精心构造的试探波函数获得能隙上界

3. 矩估计技术

通过分析能隙m阶矩Δ^m_ave的系统尺寸依赖性来确定z：

[Δ^m]_ave ∼ N^(-mz)

实验设置

数据集

• 系统尺寸：N = 10² 到 10⁴
• 无序实现：每个数据点平均1000次无序实现
• 参数范围：s ∈ 0,1, D = 1（主要结果）

评价指标

• 动力学临界指数z通过能隙的幂律标度行为确定
• 相关长度ξ和典型相关长度ξ̄的临界指数ν和ν̄

对比方法

• 与非相关无序情况(z = ∞)对比
• 与无无序情况(z = 1)对比
• 理论边界与数值结果对比

实现细节

• 使用精确对角化方法计算矩阵M的本征值
• 临界点确定：Γ_c = 1（通过条件ΔJ = ΔΓ）
• 拟合范围：N ≥ 10³和10² ≤ N ≤ 10³

实验结果

主要结果

1. 弱无序情况

• 理论结果：z = 1（精确）
• 物理意义：与无无序情况相同，相关无序不改变临界行为

2. 强无序情况

• 理论边界：max(D(1/2 + |s-1/2|) + 1/2, 1) ≤ z ≤ D + 1
• 数值结果：在D = 1时，z在s = 0.5处最小，在s ≤ 0.3和s ≥ 0.7时z ≈ 1 + |s-1/2|

3. 参数依赖性

• z在s和1-s处对称
• s = 0.5时z达到最小值
• 强无序边界在极端s值处得到很好验证

消融实验

• s参数效应：s = 0.5时有限尺寸效应最大，需要更大系统研究
• 边界紧致性：在s接近0或1时，数值结果接近理论下界
• 相关长度：ν = ν̄ = 1，与非相关无序情况不同

实验发现

1. 临界行为质变

相关无序导致z从∞降至有限值，这是定性的改变。

2. 调节机制

通过参数s可以有效调节z值，为实际应用提供了控制手段。

3. 相关长度统一

平均相关长度和典型相关长度在相关无序情况下重合，消除了样本间涨落。

相关工作

主要研究方向

1. 无无序横向场伊辛链：z = 1，完全可解
2. 非相关无序情况：z = ∞，由Griffiths-McCoy奇异性导致
3. 强无序重整化群：用于研究无序量子自旋链的理论框架

本文与相关工作关系

• 区别于Fisher等人的非相关无序研究
• 扩展了Hoyos等人的数值工作，提供精确解析结果
• 与Knysh等人关于量子退火加速的研究相呼应

优势对比

• 首次提供相关无序情况的精确解析边界
• 系统研究了参数依赖性
• 连接了理论分析与数值验证

结论与讨论

主要结论

1. 相关无序的保护效应：相关无序有效抑制了Griffiths-McCoy奇异性，使z从∞降至有限值
2. 参数可调性：通过调节s和D可以控制z的具体值
3. 应用前景：有限的z值有助于提升量子绝热计算性能

局限性

1. 一维限制：结果仅适用于一维系统
2. 特定无序形式：相关无序的具体形式可能影响普适性
3. 强无序精确值：在某些参数区域未能确定z的精确值

未来方向

1. 高维推广：研究d ≥ 2维情况下的相关无序效应
2. 其他无序类型：探索不同相关无序形式的影响
3. 实验实现：寻找实现相关无序的实验方案

深度评价

优点

1. 理论严谨性：提供了精确的数学证明和严格的边界分析
2. 方法创新性：巧妙结合矩阵不等式和变分方法获得紧致边界
3. 实用价值：为量子计算中的参数优化提供了理论指导
4. 完整性：理论分析与数值验证相互印证

不足

1. 边界间隙：强无序情况下上下界之间仍有间隙
2. 有限尺寸效应：在某些参数区域有限尺寸效应较大
3. 物理机制：对相关无序抑制奇异性的微观机制解释不够深入

影响力

1. 理论贡献：为无序量子系统理论提供了新的视角
2. 应用潜力：对量子退火和绝热计算具有重要指导意义
3. 方法价值：边界分析方法可推广到其他无序系统

适用场景

1. 量子退火优化：设计具有有限z值的量子退火协议
2. 无序系统设计：工程化设计具有特定临界性质的量子系统
3. 理论研究：作为研究其他无序量子系统的方法论基础

参考文献

论文引用了该领域的重要工作，包括：

• Sachdev的量子相变理论基础
• Fisher关于随机横向场伊辛链的开创性工作
• Hoyos等人关于相关无序的数值研究
• 量子绝热计算相关的理论和应用研究

这篇论文在理论物理和量子计算交叉领域做出了重要贡献，通过精确的数学分析揭示了相关无序对量子相变的保护效应，为实际应用提供了有价值的理论指导。