[論文レビュー] Finite density simulations using a determinant estimator

Andrei Alexandru, Anyi Li+1 more|ArXiv.org|Nov 16, 2007

Markov Chains and Monte Carlo Methods参考文献 3被引用数 5

ひとこと要約

本稿では、正確な行列式計算の$V^3$スケーリングに起因する計算コストの問題を克服するため、行列式推定器を用いたハイブリッドノイジィモンテカルロアルゴリズムを導入し、より大きなスケールの格子QCDシミュレーションを可能にする。この手法は、Pade近似と確率的ノイズ場を用いて複素数のフェルミオン行列式比を推定し、$4^4$格子上での計算で正確な計算結果と一致する高精度な結果を得た。計算コストは$V^2$スケーリングを示し、非ゼロバリオン密度における相転移の研究をより大きなスケールで行う道筋を示した。

ABSTRACT

Previous investigations have shown that the canonical approach to simulating QCD at finite density is promising. The algorithm we used in our earlier work employs an exact calculation of the fermionic determinant which limits the size of the lattices we can simulate. Interesting questions can only be answered if we simulate at larger volume. In this paper we explore an algorithm, Hybrid Noisy Monte Carlo, that employs a determinant estimator rather than an exact calculation. We first present the technical aspects of the estimator, check that the algorithm is correct by comparing it with our previous study, and then discuss its merits. We will also discuss the challenges faced when simulating larger lattice volumes.

研究の動機と目的

• 正確な行列式計算が計算的に不可能になる、より大きな格子体積における有限密度QCDシミュレーションを可能にすること。
• ギャンブル・カノニカルな手法に内在する重なり問題を回避するため、カノニカル分配関数に基づく行列式推定器手法を実装・検証すること。
• 推定器手法の性能とスケーリング特性を正確な行列式計算と比較すること。
• 非ゼロバリオン密度における相転移を、より大きな格子上でシミュレートする基盤を築くこと。

提案手法

• フェルミオン行列式を$\det M(U)^2$と比推定子$f_k(U,\xi)$に分離するため、補助場$\xi$を用いてカノニカル分配関数を再定式化する。
• 行列式の対数$\ln M$のPade近似（$K=30$）を用い、$Z(4)$ノイズベクトルを用いて$\ln M_{\phi} - \ln M$のトレースの不偏推定器を構築する。
• 対角成分を除いた行列を差し引くことで、対数差行列の非対角構造を模倣する分散低減技術を適用する。
• 2段階更新を用いたハイブリッドノイジィモンテカルロアルゴリズムを実装：ゲージ場の更新は$e^{-S_G(U)}\det M(U)^2$で重み付けされ、補助場の再設定は比$|f_k(U,\xi')/f_k(U,\xi)|$を用いる。
• $\phi$位相にわたる$\mathrm{Tr} \ln M_\phi$の観察されたフラクチュエーションに基づき、$n=20$のブレイクアップレベルを設定し、大多数の推定値が$n$未満に収まるようにする。
• 50%の受容率を達成するHMC更新を用い、ポリコフ・ループや化学ポテンシャルなどの観測量を測定することで、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 行列式推定器手法は、有限密度QCDシミュレーションにおいて正確な行列式計算の結果を正確に再現できるか？
• RQ2 推定器手法の計算コストは、正確な行列式計算と比較して、格子体積にどのように依存するか？
• RQ3 推定器の数値的安定性と効率性を維持するための最適なブレイクアップレベルは何か？
• RQ4 この推定器手法により、$6^3 \times 4$のようなより大きな格子でも実行可能なシミュレーションが可能になるか？

主な発見

• $4^4$格子上では、推定器手法がポリコフ・ループの絶対値とバリオン化学ポテンシャルの両方を、正確な行列式計算の結果と非常に良好に再現した。
• 推定器手法の計算コストは$V^2$スケーリングを示し、初期の予想である$V$スケーリングよりも著しく優れていたが、正確な行列式計算の$V^3$スケーリングよりも顕著に良好であった。
• $4^4$格子上では、推定器手法は100秒で2回の評価を完了したが、正確な手法は140秒を要した。
• 推定器のフラクチュエーションの平均と分散は、格子体積に比例して増加し、$V^2$スケーリングが、増大するフラクチュエーションに応じてブレイクアップレベルを増加させる必要があることによって駆動されていることを確認した。
• 過去の正確なシミュレーションとの比較に基づき、手法は正しいかつ安定しており、系統的バイアスの兆候は一切認められなかった。
• 著者らは、より大きな格子上での有限密度QCDシミュレーションのための最も現実的な道筋として、この推定器手法を結論づけ、$6^3 \times 4$格子上での相転移の研究を計画している。