[论文解读] Tensor renormalization group study of classical XY model on the square lattice

Ji-Feng Yu, Z. Y. Xie+7 more|arXiv (Cornell University)|Sep 19, 2013

Physics of Superconductivity and Magnetism被引用 7

一句话总结

本研究将基于高阶奇异值分解的张量重整化群（HOTRG）方法应用于正方形晶格上的经典连续XY模型，实现了在无限晶格极限下热力学性质的精确计算。该方法得到的临界相变温度为 $ T_c = 0.8921 \pm 0.0019 $，与蒙特卡罗方法和高温展开结果一致，并估算出临界指数 $ \delta = 14.5 $，与Kosterlitz的解析预测值 $ \delta = 15 $ 良好吻合。

ABSTRACT

Using the tensor renormalization group method based on the higher-order singular value decom- position, we have studied the thermodynamic properties of the continuous XY model on the square lattice. The temperature dependence of the free energy, the internal energy and the specific heat agree with the Monte Carlo calculations. From the field dependence of the magnetic susceptibility, we find the Kosterlitz-Thouless transition temperature to be 0.8921 \pm 0.0019, consistent with the Monte Carlo as well as the high temperature series expansion results. At the transition temperature, the critical exponent δis estimated as 14.5, close to the analytic value by Kosterlitz.

研究动机与目标

• 开发并应用一种新型张量重整化群方法于连续自旋模型，特别是每个格点具有无限自由度的经典XY模型。
• 通过在HOTRG中直接计算热力学极限，克服传统方法中有限尺寸外推的局限性。
• 精确确定二维XY模型中Kosterlitz-Thouless相变的临界温度和临界指数 $ \delta $。
• 将HOTRG结果与已建立的蒙特卡罗模拟和高温系列展开进行验证。
• 确立HOTRG在连续模型中的可行性与准确性，将其应用范围从离散自旋系统扩展至连续系统。

提出的方法

• 通过配分函数中余弦相互作用的傅里叶展开，利用特征展开构建张量网络表示，实现将连续自旋映射为适合HOTRG处理的离散张量网络。
• 局部张量通过修正贝塞尔函数 $ I_n(\beta) $ 定义，该函数源于配分函数中余弦相互作用的傅里叶展开。
• 采用高阶奇异值分解（HOSVD）高效收缩并粗粒化张量网络，以保留主要的纠缠和热力学信息。
• 通过粗粒化后张量网络获得的自由能对温度求导，计算热力学量——自由能、内能、比热。
• 通过分析磁化率随外场变化的依赖关系提取临界温度 $ T_c $，利用幂律拟合 $ T = T_c + a h^b $ 将峰值位置外推至零场。
• 通过在 $ T_c $ 处对磁化强度随外场变化的幂律拟合 $ m \sim h^{1/\delta} $，确定临界指数 $ \delta $，数据来自HOTRG计算。

实验结果

研究问题

• RQ1 使用HOTRG方法，二维经典XY模型中Kosterlitz-Thouless相变的临界温度 $ T_c $ 是多少？
• RQ2 HOTRG方法在计算连续XY模型的自由能、内能和比热等热力学量时，其精度如何？
• RQ3 在Kosterlitz-Thouless相变点，临界指数 $ \delta $ 的值是多少？与解析预测相比如何？
• RQ4 与蒙特卡罗模拟和高温系列展开相比，HOTRG方法在该模型上的表现如何？
• RQ5 HOTRG方法能否成功扩展至连续自旋模型？这类模型传统上需要不同的张量构造技术。

主要发现

• 二维XY模型中Kosterlitz-Thouless相变的临界温度确定为 $ T_c = 0.8921 \pm 0.0019 $，与蒙特卡罗结果（$ T_c = 0.89294(8) $）和高温系列结果（$ T_c = 0.89286(8) $）高度一致。
• 通过HOTRG计算得到的自由能、内能和比热随温度的变化，与边长为 $ 256 \times 256 $ 的蒙特卡罗模拟结果在定量上表现出极佳的一致性。
• HOTRG计算中比热的峰值出现在 $ T_{\text{peak}} = 1.04 $，略高于蒙特卡罗结果的 $ T_{\text{peak}} = 1.03 $，与以往蒙特卡罗研究结果一致。
• 磁化率表现出一个尖锐的峰值，且随着外磁场减小，峰值向低温方向移动；其在零场下的外推结果给出临界温度，不确定性仅为 $ \pm 0.0019 $。
• 通过在 $ T_c $ 处对磁化强度随外场变化的幂律拟合，估算出临界指数 $ \delta = 14.5 $，与Kosterlitz的解析预测值 $ \delta = 15 $ 非常接近。
• HOTRG方法无需有限尺寸外推即可成功计算连续XY模型的热力学性质，证明了其在连续系统中的准确性和潜力。