[論文レビュー] Winding number expansion for the canonical approach to finite density simulations

Xiangfei Meng, Anyi Li+2 more|arXiv (Cornell University)|Nov 13, 2008

Markov Chains and Monte Carlo Methods参考文献 4被引用数 8

ひとこと要約

本稿では、有限密度 lattice QCD シミュレーションにおけるクォーク行列式の射影計算の安定化と高速化を目的として、巻き数展開法（WNEM）を導入する。時間的格子方向まわりのクォークループの巻き数を用いて行列式の対数を展開することで、k=50 に達するような大きなクォーク数に対しても、正確かつ数値的に安定した射影行列式の評価が可能となり、k≈20 を超えると不安定化する離散フーリエ変換法の欠陥を克服する。

ABSTRACT

The canonical partition function approach was designed to avoid the overlap problem that affects the lattice simulations of nuclear matter at high density. The method employs the projections of the quark determinant on a fix quark number sector. When the quark number is large, the evaluation of the projected determinant becomes numerically unstable. In this paper a different evaluation method based on expanding the determinant in terms of loops winding around the lattice is studied. We show that this method is stable and significantly faster than our original algorithm. This greatly expands the range of quark numbers that we can simulate effectively.

研究の動機と目的

• カノニカル lattice QCD シミュレーションにおける、大規模なクォーク数に対する射影行列式を計算する際の離散フーリエ変換法の数値的不安定性を解消すること。
• カノニカルアプローチにおけるクォーク行列式射影を評価するための、離散フーリエ変換法の代替として安定的かつ効率的な手法を開発すること。
• 特に大きな格子サイズや高いバリオン密度を想定した状況において、利用可能なクォーク数の範囲を拡張すること。
• 高精度な離散フーリエ変換結果と比較して、巻き数展開法（WNEM）の正確性と安定性を検証すること。

提案手法

• フェルミオン行列式の対数を、時間的格子方向まわりのクォークループの巻き数に基づいて展開し、Tr log M(U,φ) を巻き数寄与の和に分解する。
• 行列式を φ に関するフーリエ級数の指数関数としてモデル化する：exp(A₀ + Σₙ(eⁱⁿφWₙ + e⁻ⁱⁿφWₙ†)，ここで Wₙ は巻き数 n のループの重みである。
• 係数 Aₙ は、Nₔ 個の離散的 φ 点における log det M(U,φ) の離散フーリエ変換により計算され、Nₔ = 16 が高精度に十分であることが判明した。
• 射影行列式は、Aₙ 係数から解析的に6次までの級数展開を用いて計算され、高周波数のフーリエ成分の直接評価を回避する。
• 高次のフーリエ成分が指数的に減衰する対数行列式空間で作業することで、離散フーリエ変換の数値的不安定性を回避する。
• WNEM の正確性と安定性は、Nₔ = 204 の高精度な離散フーリエ変換結果と比較して検証され、k ≤ 50 の範囲で相対誤差 10⁻⁶ 以内の一致が確認された。

実験結果

リサーチクエスチョン

• RQ1 巻き数展開法は、大規模なクォーク数におけるクォーク行列式射影を計算する際、離散フーリエ変換法の代替として数値的に安定した手法を提供できるか？
• RQ2 クォーク数 k=50 までのシミュレーションで十分な正確性を達成するために、巻き数展開に必要な項数はどの程度か？
• RQ3 離散フーリエ変換法の不安定化閾値（k > 20）を超えるクォーク数に対しても、WNEM 法は正確性と安定性を維持できるか？
• RQ4 Aₙ 係数を十分な精度で計算するために必要な最小の φ 点数（Nₔ）はいくつか？
• RQ5 k が大きい場合に WNEM の計算コストは、離散フーリエ変換法と比べてどの程度か？

主な発見

• 表1の結果から、クォーク数 k > 20 の場合、離散フーリエ変換法は Nₔ の値に依存して一貫性のない結果を示し、数値的不安定性が顕著であることが示された。
• 巻き数展開法（WNEM）は、クォーク数 k = 50 まで安定的かつ正確に動作し、k が増加するに従い射影行列式が指数関数的に減少することが確認された。
• φ 点数 Nₔ = 16 のみを用いても、Aₙ 係数の精度が十分に確保され、Nₔ = 204 の結果と比較して射影行列式の相対誤差が 10⁻⁶ 未満であった。
• 6次までの WNEM 近似は、k ≤ 50 の範囲で7次近似と10⁻⁴ 未満の相対誤差で一致しており、高次項の寄与が無視できることが示された。
• WNEM は離散フーリエ変換法よりも著しく高速であり、精度を損なわず、はるかに大きなクォーク数でのシミュレーションが可能になった。
• 本手法により、特に高バリオン密度を要するシミュレーション（例：6³×4 格子における k ~ 30）において、利用可能な位相空間が著しく拡張された。