2025-11-17T14:04:12.377243

Clifford circuits Augmented Matrix Product States for fermion systems

Huang, Qian, Qin
Clifford circuits Augmented Matrix Product States (CAMPS) was recently proposed to leverage the advantages of both Clifford circuits and Matrix Product States (MPS). Clifford circuits can support large entanglement and can be efficiently simulated classically according to the Gottesman-Knill theorem. So in CAMPS, MPS needs only to handle the so-called Non-stabilizerness Entanglement Entropy which significantly improves the simulation accuracy for a given bond dimension. In this work, we generalize CAMPS to study the Fermion system by taking advantage of the Jordan-Wigner transformation which can map the studied Fermion system to a spin system. We benchmark the method on both the spinless $t-V$ model and the spinful Hubbard model. Our test results show significant improvement of the accuracy of CAMPS over MPS, especially when the interactions are strong. Fermionic CAMPS provides a useful tool for the accurate study of many-body fermion systems in the future and has the potential to help resolve long-standing issues.
academic

Clifford circuits Augmented Matrix Product States for fermion systems

基本信息

  • 论文ID: 2501.00413
  • 标题: Clifford circuits Augmented Matrix Product States for fermion systems
  • 作者: Jiale Huang, Xiangjian Qian, Mingpu Qin
  • 分类: cond-mat.str-el quant-ph
  • 发表时间: January 3, 2025
  • 论文链接: https://arxiv.org/abs/2501.00413

摘要

Clifford circuits Augmented Matrix Product States (CAMPS) was recently proposed to leverage the advantages of both Clifford circuits and Matrix Product States (MPS). Clifford circuits can support large entanglement and can be efficiently simulated classically according to the Gottesman-Knill theorem. So in CAMPS, MPS needs only to handle the so-called Non-stabilizerness Entanglement Entropy which significantly improves the simulation accuracy for a given bond dimension. In this work, we generalize CAMPS to study the Fermion system by taking advantage of the Jordan-Wigner transformation which can map the studied Fermion system to a spin system. We benchmark the method on both the spinless t−Vt-V model and the spinful Hubbard model. Our test results show significant improvement of the accuracy of CAMPS over MPS, especially when the interactions are strong. Fermionic CAMPS provides a useful tool for the accurate study of many-body fermion systems in the future and has the potential to help resolve long-standing issues.

研究背景与动机

研究问题

本文要解决的核心问题是如何高效、准确地模拟强关联多体费米子系统。传统的数值方法在处理这类系统时面临着显著的计算挑战，特别是在强相互作用区域。

问题重要性

强关联多体费米子系统的研究是现代凝聚态物理学的基石，对理解超导、磁性和量子相变等现象至关重要。这些系统的准确模拟对于：

  1. 理解高温超导体的物理机制
  2. 探索新型量子材料
  3. 解决长期存在的理论问题具有重要意义

现有方法局限性

  1. 量子蒙特卡罗方法：存在符号问题，在某些参数区域难以应用
  2. 密度矩阵重整化群(DMRG)/矩阵乘积态(MPS)：在高维系统中受到纠缠结构限制，计算成本过高
  3. 二维张量网络方法：如PEPS、MERA等虽能编码二维纠缠结构，但计算成本通常很高，难以达到大键维

研究动机

CAMPS方法已在自旋系统中显示出显著优势，能够有效减少纠缠熵并提高模拟精度。本文的动机是将这一成功的方法扩展到费米子系统，利用Jordan-Wigner变换建立费米子-自旋系统的映射关系。

核心贡献

  1. 方法扩展：首次将CAMPS方法从自旋系统推广到费米子系统，通过Jordan-Wigner变换实现费米子到自旋算符的映射
  2. 算法实现：提出了适用于费米子系统的CAMPS算法框架，保持了Clifford电路的高效经典模拟特性
  3. 基准测试：在无自旋t-V模型和有自旋Hubbard模型上进行了系统的基准测试，验证了方法的有效性
  4. 精度提升：实验结果显示，特别是在强相互作用区域，CAMPS相比传统MPS方法具有显著的精度提升
  5. 理论洞察：为费米子系统的Non-stabilizerness纠缠熵概念提供了实际应用，展示了其在量子多体系统经典模拟中的重要作用

方法详解

任务定义

本文的任务是开发一种高效的数值方法来计算强关联费米子系统的基态性质，包括基态能量和纠缠结构。输入为费米子哈密顿量参数，输出为基态波函数的高精度近似和相关物理量。

模型架构

Jordan-Wigner变换

费米子产生算符ci†c_i^†和湮灭算符cic_i满足反对易关系：

{c_i, c_j} = {c_i^†, c_j^†} = 0, {c_i, c_j^†} = δ_{ij}

Jordan-Wigner变换将这些算符映射为字符串算符：

c_i^† = σ_1^z σ_2^z ⋯ σ_{i-1}^z σ_i^+
c_i = σ_1^z σ_2^z ⋯ σ_{i-1}^z σ_i^-

其中σi±=(σix±iσiy)/2σ_i^± = (σ_i^x ± iσ_i^y)/2，密度算符ni=ci†cin_i = c_i^†c_i映射为(1+σiz)/2(1 + σ_i^z)/2。

CAMPS架构

CAMPS波函数定义为：

|CAMPS⟩ = C|MPS⟩

其中C表示Clifford电路，MPS ansatz为：

|MPS⟩ = ∑_{σ_i} Tr(A_1^{σ_1} A_2^{σ_2} ⋯ A_N^{σ_N})|σ_1σ_2⋯σ_N⟩

观测量计算

任何观测量O=∑iPiO = ∑_i P_i的期望值计算为：

⟨O⟩ = ⟨CAMPS|O|CAMPS⟩ = ⟨MPS|O'|MPS⟩

其中O′=C†OCO' = C^†OC。由于Clifford电路将Pauli字符串映射为另一个Pauli字符串，O′O'保持相同的计算复杂度。

技术创新点

  1. 费米子-自旋映射：巧妙利用Jordan-Wigner变换，使得原本只适用于自旋系统的CAMPS方法能够处理费米子系统
  2. 纠缠分离：通过Clifford电路处理可稳定化的纠缠部分，MPS只需处理Non-stabilizerness纠缠熵(NsEE)，显著提高了给定键维下的模拟精度
  3. 计算效率：根据Gottesman-Knill定理，Clifford电路可以在经典计算机上高效模拟，避免了指数级的计算复杂度
  4. 二维扩展：通过蛇形映射将二维晶格转换为一维链，使得方法能够处理高维费米子系统

实验设置

测试模型

2D t-V模型

哈密顿量为：

H = -t ∑_{⟨i,j⟩} c_i^† c_j + h.c. + V ∑_{⟨i,j⟩} (n_i - 1/2)(n_j - 1/2)

测试参数：

  • 系统尺寸：6×6, 8×8
  • 相互作用强度：V = 0.5, 1, 3
  • 填充：半填充
  • 边界条件：开边界

Hubbard模型

哈密顿量为：

H = -t ∑_σ ∑_{⟨i,j⟩} c_{i,σ}^† c_{j,σ} + h.c. + U ∑_i n_{i,↑}n_{i,↓} - μ ∑_i n_i

测试参数：

  • 系统尺寸：1×32, 2×16, 4×8
  • 相互作用强度：U = 8
  • 掺杂：1/8空穴掺杂
  • 边界条件：开边界

评价指标

  1. 相对误差：基态能量的相对误差 (E−Eref)/∣Eref∣(E - E_{ref})/|E_{ref}|
  2. 纠缠熵：MPS部分在中心键的纠缠熵
  3. 改进比率：定义为 (EMPS−Eref)/(ECAMPS−Eref)(E_{MPS} - E_{ref})/(E_{CAMPS} - E_{ref})

参考基准

使用大键维DMRG计算得到的高精度基态能量作为参考值，通过外推确保参考值的可靠性。

实验结果

主要结果

t-V模型结果

对于6×6系统：

  • V = 0.5时，CAMPS在D = 2000时改进比率达到约4倍
  • V = 1时，改进比率约为3.5倍
  • V = 3时，改进比率达到约4倍，且随键维增加而提升

对于8×8系统：

  • V = 0.5时，改进比率约为2.5倍
  • V = 1时，改进比率约为3倍
  • V = 3时，改进比率达到4倍，显示出更强的相互作用下CAMPS优势更明显

Hubbard模型结果

  • 1×32系统：改进比率达到约7.5倍
  • 2×16系统：改进比率约为5倍
  • 4×8系统：改进比率约为2.75倍

关键发现

  1. 强相互作用优势：在所有测试案例中，相互作用越强，CAMPS相对于MPS的优势越明显
  2. 纠缠熵减少：CAMPS在MPS部分的纠缠熵显著低于纯MPS方法
  3. 可扩展性：改进比率随着键维D的增加而提升，显示出良好的可扩展性
  4. 系统尺寸依赖性：对于不同系统尺寸，CAMPS都表现出一致的改进效果

精度提升量化

以8×8 t-V模型V=3为例：

  • 在相同键维D=1000下，CAMPS的相对误差比MPS低约一个数量级
  • 纠缠熵从MPS的约0.58减少到CAMPS的约0.48

相关工作

张量网络方法发展

  1. MPS/DMRG：在一维和准一维系统中表现优异，但受限于面积定律
  2. 二维张量网络：PEPS、MERA、PESS等方法能处理二维系统，但计算成本高
  3. 增强型MPS：FAMPS等方法通过解纠缠器增强MPS的纠缠编码能力

CAMPS方法演进

  1. 原始CAMPS：最初针对自旋系统的基态计算
  2. 时间演化扩展：结合TDVP实现时间演化模拟
  3. 有限温度推广：扩展到有限温度计算
  4. 临界系统应用：用于解纠缠临界量子链

费米子-自旋映射

  1. Jordan-Wigner变换：经典的费米子-自旋映射方法
  2. Bravyi-Kitaev变换：在某些情况下更高效的替代方案
  3. 其他映射技术：近年来发展的各种改进映射方法

结论与讨论

主要结论

  1. 成功扩展：成功将CAMPS方法从自旋系统扩展到费米子系统，验证了方法的普适性
  2. 显著改进：在强相互作用区域，CAMPS相比传统MPS方法实现了2-7倍的精度提升
  3. 纠缠优化：通过Clifford电路有效处理可稳定化纠缠，使MPS专注于处理NsEE
  4. 实用价值：为研究强关联费米子系统提供了新的高效数值工具

局限性

  1. 映射依赖：方法的效果依赖于Jordan-Wigner变换，在二维系统中会产生长程字符串算符
  2. 对称性限制：当前实现中未强加U(1)对称性，需要通过化学势调节填充因子
  3. 边界条件：目前测试仅限于开边界条件，周期边界条件的效果有待验证
  4. 系统尺寸：受限于计算资源，测试的系统尺寸相对较小

未来方向

  1. 对称性集成：在Clifford电路中引入U(1)对称性以提高计算效率
  2. 映射方案比较：测试Bravyi-Kitaev等其他费米子-自旋映射方案的效果
  3. 方法扩展：推广到实时演化和有限温度计算
  4. 应用拓展：应用于解决高温超导等长期存在的物理问题

深度评价

优点

  1. 方法创新性：首次将CAMPS成功扩展到费米子系统，开辟了新的研究方向
  2. 理论基础扎实：基于Gottesman-Knill定理的理论基础确保了方法的可行性
  3. 实验充分性：在多个模型和参数设置下进行了系统的基准测试
  4. 结果说服力强：数值结果清晰展示了方法的优势，特别是在强相互作用区域
  5. 写作清晰：论文结构清晰，技术细节描述准确

不足

  1. 系统尺寸限制：测试的系统尺寸相对较小，大系统的效果有待验证
  2. 边界条件单一：仅测试了开边界条件，缺乏周期边界条件的结果
  3. 对称性处理：未充分利用系统的对称性，可能影响计算效率
  4. 理论分析不足：缺乏对方法为何在强相互作用下表现更好的深入理论分析

影响力

  1. 学术贡献：为强关联费米子系统的数值研究提供了新工具
  2. 实用价值：有望应用于高温超导等重要物理问题的研究
  3. 方法论意义：展示了量子信息方法在凝聚态物理中的应用潜力
  4. 可复现性：方法描述清晰，易于其他研究者复现和扩展

适用场景

  1. 强关联费米子系统：特别适用于强相互作用区域的计算
  2. 二维量子材料：可用于研究二维电子系统的量子相变
  3. 高温超导研究：为理解铜氧化物超导体提供新的计算工具
  4. 量子模拟：为量子模拟器的基准测试提供经典算法支持

参考文献

本文引用了49篇重要文献，涵盖了强关联电子系统、张量网络方法、量子信息理论等多个领域的关键工作，为研究提供了坚实的理论基础和方法论支撑。

总体评价：这是一篇高质量的方法论论文，成功地将CAMPS方法扩展到费米子系统，为强关联多体系统的数值研究提供了重要的新工具。虽然在系统尺寸和对称性处理方面还有改进空间，但其创新性和实用价值使其成为该领域的重要贡献。