[论文解读] Three-Dimensional Ising Model and Transfer Matrices

Shuwen Lou, Sheng‐Jhih Wu|arXiv (Cornell University)|Mar 22, 2000

Theoretical and Computational Physics参考文献 1被引用 8

一句话总结

本文将 B. Kaufman 的转移矩阵方法应用于三维伊辛模型，利用旋量表示和直积空间中的二维旋转，推导出单位点内能的高温展开。尽管所得级数的系数与计算机图形法结果不同，但结构高度一致，表明在保持本征值简化与准确捕捉最大本征值之间存在权衡，所提出的修正因子 Φ(x) 可使结果与已知结果重新一致。

ABSTRACT

The use of a transfer matrix method to solve the 3D Ising model is straightforwardly generalized from the 2D case. We follow B.Kaufman's approach. No approximation is made, however the largest eigenvalue cannot be identified. This problem comes from the fact that we follow the choice of directions of 2-dimensional rotations in the direct product space of the 2D Ising model such that all eigenvalue equations reduce miraculously to only one equation. Other choices of directions of 2-dimensional rotations for finding the largest eigenvalue may lose this fascinating feature. Comparing the series expansion of internal energy per site at the high temperature limit with the series obtained from the computer graphic method, we find these two series have very similar structures. A possible correct via a factor Φ(x) is suggested to fit the result of the graphic method.

研究动机与目标

• 将 B. Kaufman 的方法推广至三维伊辛模型，实现二维转移矩阵方法的推广。
• 研究在直积空间中使用旋量表示与二维旋转是否能简化三维模型中的本征值方程。
• 将本方法所得的单位点内能高温级数与计算机图形法结果进行比较。
• 识别级数系数差异的根源，并探讨可能的修正因子以调和结果差异。

提出的方法

• 使用转移矩阵 V 表征三维伊辛模型的配分函数，将 V 分解为 x、y 和 z 方向最近邻相互作用的 V₁、V₂、V₃。
• 在转移矩阵形式中应用旋量表示，实现在自旋态直积空间中的二维旋转。
• 在热力学极限下推导转移矩阵的本征值方程，导出单一函数 γ(θ₁, θ₂)。
• 利用关系式 coshγ = c³/s − c(cosθ₁ + cosθ₂) + sc cosθ₁ cosθ₂ + s² sinθ₁ sinθ₂ 表达自由能与内能。
• 通过 x = tanhK 进行高温展开，使用符号计算（Maple V）将 u 表示为 x 的级数。
• 提出一个修正因子 Φ(x)，用于乘以相互作用项，以恢复与计算机图形法已知结果的一致性。

实验结果

研究问题

• RQ1 能否通过旋量表示与直积空间中的二维旋转，一致地将二维伊辛模型的转移矩阵方法推广至三维伊辛模型？
• RQ2 为何该转移矩阵方法与计算机图形法所得的内能高温级数系数存在差异，尽管结构高度相似？
• RQ3 是否存在一个修正因子 Φ(x)，可在保持本征值方程简化的同时调和级数系数的差异？
• RQ4 选择二维旋转方向在平衡本征值简化与最大本征值精度方面起到何种作用？
• RQ5 是否可通过修正的 coshγ 表达式在不损失本征值约化结构优雅性的情况下，恢复与已知结果的一致性？

主要发现

• 所得单位点内能的高温级数为 u = −3x − 8x³ − 28x⁵ − 132x⁷ − 832x⁹ + O(x¹¹)，与计算机图形法结果不同。
• 计算机图形法的结果为 u = −3x − 12x³ − 120x⁵ − 1332x⁷ − 17676x⁹ + O(x¹¹)，结构相同但系数不同。
• 首项非零项（−3x）及 x 的偶次幂的缺失在两种方法中均一致。
• 差异源于无法精确确定转移矩阵的最大本征值，可能与旋转方向选择中的相互冲突要求有关。
• 修正因子 Φ(x) = 1 − x² − (27/2)x⁴ − (249/2)x⁶ − (12325/8)x⁸ + O(x¹⁰) 可使结果与计算机图形法级数重新一致。
• 所提出的 Φ(x) 作为不同二维旋转方向的加权函数，为解决三维伊辛模型的精确解提供了可能路径。