[論文レビュー] 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

ひとこと要約

この論文は、2次元エドワーズ＝アンドリュース・イジンスピンガラスのような frustrated スピン系において、標準的な64ビットデータ型の数値精度の制限が原因で、テンソル縮約群（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）手法が、frustrated スピン系において物理的に不適切な負の分配関数値を生じる理由を調査すること。
• その根本的原因が、不十分な切断カットオフパラメータにあるのか、それとも標準的なデータ型における数値精度の制限にあるのかを特定すること。
• 低温におけるスピンガラスの熱力学的性質を研究するためのTRGの可能性を評価すること。
• 強磁性イジンスモデルとfrustrated スピンガラスにおけるTRGの数値的挙動を比較すること。
• スピンガラスの信頼できるTRGシミュレーションに必要な精度要件を定量化すること。

提案手法

• C++を用いて、二項分布結合（p = 0.5）を有する正方形格子上の2次元エドワーズ＝アンドリュース・イジンスピンガラスをシミュレートし、任意精度算術にGNU MPFRを使用した。
• テンソルの縮約と高次特異値分解（HOSVD）による切断を用いて、ボンド次元を低減するTRGアルゴリズムを実装した。
• 局所的テンソル分解を、逆温度の双曲線関数を用いて定義された行列W¹とW²に基づいて行った。
• 粗化過程中に分配関数を追跡し、負の値が生じるかを失敗の指標として監視した。
• 複数の不規則性実現について、倍精度（64ビット）と高精度（1024ビット）浮動小数点算術の結果を比較した。
• 強磁性系とスピンガラス系におけるテンソル成分のキャンセル効果を比較するために、テンソル成分の分布を分析した。

実験結果

リサーチクエスチョン

• RQ1 なぜ2次元エドワーズ＝アンドリュース・スピンガラスのTRGシミュレーションは、低温で物理的に不適切な負の分配関数値を生じるのか？
• RQ2 この失敗は、小さな切断カットオフパラメータのためなのか、それとも標準データ型における数値精度の制限のためなのか？
• RQ3 強磁性イジンスモデルとfrustrated スピンガラスにおけるTRGの数値的挙動にはどのような違いがあるか？
• RQ4 データ型の精度を向上させることで、スピンガラスシミュレーションにおけるTRGの失敗率はどの程度低下するか？
• RQ5 スピンガラスのTRGにおいて物理的に意味のある結果を得るためには、どの程度の精度レベルが必要か？

主な発見

• 64ビット倍精度を用いた2次元エドワーズ＝アンドリュース・スピンガラスのTRGシミュレーションにおいて、低温域で分配関数の失敗率が急激に上昇する。
• 1024ビットの高精度算術を用いることで、失敗率は完全に解消され、負の分配関数値の原因は精度の制限にあることが示唆される。
• スピンガラスでは正と負のテンソル成分の近接キャンセルが生じ、これには高い精度が要求されるが、強磁性系ではすべての成分が正であるため、そのような要求はない。
• TRG手法は理論的にはスピンガラスに適用可能ではあるが、精度の制限による数値的不安定性のため、低温シミュレーションには実用的ではない。
• 失敗率は、データ型の精度が向上するにつれて単調に減少し、TRGの安定性において精度が切断カットオフよりも主要因であることが確認された。
• 高精度要件に伴い、大きな数値的オーバーヘッドが生じ、非frustrated モデルで成功を収めた一方で、frustrated スピン系ではTRGは非効率である。