[论文解读] Matrix product decomposition and classical simulation of quantum dynamics in the presence of a symmetry

Surendra Singh, Hongyi Zhou+1 more|arXiv (Cornell University)|Jan 18, 2007

Quantum many-body systems被引用 23

一句话总结

本文提出了一种适用于具有 SU(2) 对称性的量子多体系统的对称性自适应矩阵乘积态（MPS）表示方法，从而显著提升了时间演化块稀释（TEBD）算法的效率。通过利用 SU(2) 对称性，基于总自旋基分解与 Clebsch-Gordan 耦合的改进型 MPS 结构，该方法在模拟临界自旋链时实现了更高的精度与效率，尤其在长程关联函数的计算中表现突出。

ABSTRACT

We propose a refined matrix product state representation for many-body quantum states that are invariant under SU(2) transformations, and indicate how to extend the time-evolving block decimation (TEBD) algorithm in order to simulate time evolution in an SU(2) invariant system. The resulting algorithm is tested in a critical quantum spin chain and shown to be significantly more efficient than the standard TEBD.

研究动机与目标

• 开发一种显式编码 SU(2) 对称性的矩阵乘积态表示方法，适用于多体量子系统。
• 将 TEBD 算法扩展以在时间演化过程中保持并利用 SU(2) 对称性，降低计算成本。
• 实现对具有长程关联函数的临界量子自旋链的精确模拟。
• 通过纯化技术将该方法推广至混合态及具有确定自旋量子数的态。

提出的方法

• 引入一种基于总自旋基（TSB）的改进型矩阵乘积态表示方法，将希尔伯特空间分解为 SU(2) 的不可约表示（irreps）。
• 利用 Clebsch-Gordan 系数对 SU(2) singlet 态进行二分分解，施加约束条件 $ j_1 = j_2 $ 且 $ m_1 = -m_2 $。
• 通过将 SU(2) 结构直接编码至 MPS 张量中，推导出保持对称性的 TEBD 算法，从而降低有效键维数。
• 采用纯化方法对临界自旋-1/2 反铁磁海森堡链进行模拟，处理混合态。
• 通过密度矩阵的 singlet 纯化，将 SU(2) 不变的混合态表示为 SU(2) MPS，从而实现对具有确定 $ j $ 的态的高效模拟。
• 采用张量网络结构，其中系数受对称性约束，确保时间演化过程中总自旋的精确守恒。

实验结果

研究问题

• RQ1 能否对矩阵乘积态表示进行优化，以在多体量子系统中显式编码 SU(2) 对称性？
• RQ2 如何修改 TEBD 算法，使其在时间演化过程中保持并利用 SU(2) 对称性？
• RQ3 与标准 MPS 相比，使用对称性自适应 MPS 在模拟临界自旋链时具有何种计算优势？
• RQ4 该方法能否以高精度准确计算长程两点关联函数？
• RQ5 该形式化方法能否推广至表示混合态及具有确定自旋量子数的态？

主要发现

• 在计算 $ r=1 $ 处的两点关联函数时，SU(2)-对称 TEBD 算法实现了 9 位有效数字的精度，而标准 MPS 在相似计算成本下仅能达到 6 位有效数字。
• 当 $ r=13,000 $ 时，该方法在渐近两点关联函数中的误差为 10%，而标准 MPS 在 $ r \approx 500 $ 时即达到 10% 误差。
• 该方法可高精度计算至 $ r=20,000 $ 的两点关联函数 $ C_2(r) $，在长距离上表现出优异的稳定性和精度。
• 在 SU(2) MPS 中，键维数 $ \chi \approx 2200 $ 的精度优于标准 MPS 中 $ \chi \approx 1000 $ 的结果，表明计算效率有显著提升。
• 该方法可实现对具有大局部希尔伯特空间的系统（如自旋-$ s=4 $ 链或多重腿自旋梯）的模拟，而这些系统在传统方法下通常难以处理。
• 该形式化方法可通过纯化与投影从纯化态的 SU(2) MPS 推广至混合态及具有确定 $ j $ 和 $ m $ 的纯态。