[论文解读] Do tensor renormalization group methods work for frustrated spin systems?

Zheng Zhu, Helmut G. Katzgraber|arXiv (Cornell University)|Mar 18, 2019

Theoretical and Computational Physics被引用 7

一句话总结

本文表明，张量重正化群（TRG）方法在如二维爱德华-艾因斯坦伊辛自旋玻璃这类几何阻挫自旋系统中失效，原因并非截断截断阈值过小，而是标准64位数据类型的数值精度有限。采用高精度算术（例如1024位）可消除非物理的负配分函数值，揭示出尽管理论上可行，TRG在自旋玻璃系统中仍存在数值不稳定性。

ABSTRACT

No. To illustrate that tensor renormalization group methods are poorly suited for frustrated magnetic systems, we study the thermodynamic properties of the two-dimensional Edwards-Anderson Ising spin-glass model on a square lattice. We show that the limited precision of standard 64-bit data types and not a small cut-off parameter is the main reason for unphysical negative partition function values in spin glasses.

研究动机与目标

• 研究张量重正化群（TRG）方法为何在阻挫自旋系统中产生非物理的负配分函数值。
• 确定根本原因是否为截断截断阈值不足，或标准数据类型中数值精度有限。
• 评估TRG在低温下研究自旋玻璃热力学性质的可行性。
• 比较TRG在铁磁伊辛模型与阻挫自旋玻璃系统中的数值行为。
• 量化TRG模拟自旋玻璃系统所需的精度水平。

提出的方法

• 使用C++语言，基于GNU MPFR实现的任意精度算术，对具有双模耦合（p = 0.5）的二维爱德华-艾因斯坦伊辛自旋玻璃在正方形晶格上进行模拟。
• 采用TRG算法，结合张量收缩与高阶奇异值分解（HOSVD）截断，以降低键维数。
• 基于玻尔兹曼权重实现局部张量分解，其中矩阵W¹与W²通过逆温度的双曲函数定义。
• 在粗粒化过程中追踪配分函数，并将负值作为失败指标。
• 在多个无序实现下，对比双精度（64位）与高精度（1024位）浮点算术的结果。
• 分析张量分量的分布，比较铁磁系统与自旋玻璃系统中的抵消效应。

实验结果

研究问题

• RQ1 为何在低温下TRG模拟二维爱德华-艾因斯坦自旋玻璃会产生非物理的负配分函数值？
• RQ2 该失败是由于截断截断阈值过小，还是标准数据类型中数值精度有限所致？
• RQ3 TRG在铁磁伊辛模型与阻挫自旋玻璃系统中的数值行为有何不同？
• RQ4 提高数据类型精度在多大程度上可降低TRG在自旋玻璃模拟中的失败率？
• RQ5 TRG模拟自旋玻璃系统时，达到物理上有意义结果所需的精度水平是多少？

主要发现

• 在使用64位双精度时，TRG模拟二维爱德华-艾因斯坦自旋玻璃的配分函数失败率在低温下急剧上升。
• 采用1024位高精度算术可完全消除失败率，表明限制精度而非截断阈值是导致负配分函数值的主要原因。
• 自旋玻璃系统中正负张量分量的近似抵消需要更高精度，而铁磁系统中所有分量均为正值，无需如此高精度。
• TRG方法在理论上仍适用于自旋玻璃系统，但因精度限制导致的数值不稳定性使其在低温模拟中不切实际。
• 失败率随数据类型精度提高而单调下降，证实精度是TRG稳定性的主导因素，而非截断截断阈值。
• 高精度需求带来显著的数值开销，使得尽管在非阻挫模型中表现成功，TRG对阻挫自旋系统仍效率低下。