[论文解读] Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions

Frank Verstraete, J. I. Cirac|arXiv (Cornell University)|Jul 2, 2004

Quantum many-body systems被引用 503

一句话总结

本文将投影纠缠对态（PEPS）作为矩阵乘积态（MPS）在二维及以上维度的自然推广，实现了对量子多体系统高效变分模拟。作者开发了一种可扩展的算法以计算关联函数，并将其应用于基态和虚时演化，即使在较小的键维数D下也实现了高精度，如在二维海森堡模型和受挫自旋系统中的展示所示。

ABSTRACT

We describe quantum many--body systems in terms of projected entangled--pair states, which naturally extend matrix product states to two and more dimensions. We present an algorithm to determine correlation functions in an efficient way. We use this result to build powerful numerical simulation techniques to describe the ground state, finite temperature, and evolution of spin systems in two and higher dimensions.

研究动机与目标

• 开发一种可扩展的数值方法，用于模拟二维及以上维度的量子多体系统，其中传统方法如DMRG和蒙特卡洛方法面临局限。
• 通过辅助系统的纠缠对，将矩阵乘积态（MPS）推广至高维，提出投影纠缠对态（PEPS）的概念。
• 在PEPS中实现关联函数的高效计算，这对变分方法和时间演化模拟至关重要。
• 展示基于PEPS的算法在寻找基态及模拟自旋系统（包括受挫模型）的真实或虚时演化方面的有效性。
• 建立一种适用于各种几何结构、有限温度及具有退相干的开放量子系统的框架。

提出的方法

• 提出PEPS作为MPS的推广，其中每个物理自旋由键维数为D的辅助系统之间的纠缠对投影而成。
• 将量子态表示为张量网络，其中局部张量（投影算符）将辅助自由度映射到物理自旋，相邻位点之间存在纠缠。
• 通过张量网络收缩开发高效算法以计算PEPS中的关联函数，从而支持变分优化。
• 使用Trotter分解实现时间演化算法，将时间步长划分为子步，使得每次仅一个相互作用项激活。
• 在每个子步后通过最优逼近减小键维数，逐行迭代优化投影算符，以保持计算成本较低。
• 采用多项式时间收缩技术处理计算复杂度，确保随系统尺寸N和键维数D的可扩展性。

实验结果

研究问题

• RQ1 能否将矩阵乘积态推广至二维及以上维度，以实现对量子多体系统的可扩展模拟？
• RQ2 能否在PEPS框架内开发出高效计算关联函数的算法，以支持变分方法？
• RQ3 基于PEPS的算法能否准确模拟二维自旋系统的基态和虚时演化，包括受挫模型？
• RQ4 对于如二维海森堡模型等系统，PEPS近似的精度如何随键维数D的增加而变化？
• RQ5 PEPS形式化能否扩展至有限温度态、真实时间演化及具有退相干的系统？

主要发现

• PEPS形式化成功地将类似DMRG的变分方法推广至二维及以上维度，实现了对强关联系统的精确模拟。
• 对于4×4海森堡反铁磁体，当D=3时，变分能量与精确基态能量的偏差仅为0.004%，显示出随D快速收敛。
• 即使D=2时，PEPS算法也与精确结果高度一致，能量偏差从D=1时的35%降至D=2时的2%。
• 该算法能快速收敛至基态，且随着D增大，受挫海森堡模型的能量间隙被更准确地捕获。
• 在10×10晶格上的模拟表明，该方法具有良好的可扩展性，D=3在基态能量估计上显著优于D=2。
• 该方法对反铁磁和受挫相互作用均表现出鲁棒性，证明了其通用性及在更广泛应用中的潜力。