[論文レビュー] A proposal for simulating QCD at finite chemical potential on the lattice

B. Allés, E. M. Moroni|arXiv (Cornell University)|Jun 24, 2002

Quantum Chromodynamics and Particle Interactions被引用数 5

ひとこと要約

本稿では、フェルミオン行列式の実部 Re(det D) をボルツマン重みとして用いることで、有限の化学ポテンシャルにおける格子QCDシミュレーション手法を提案する。この手法により、特定の離散的変換のもとで虚数部が対になって相殺されることで符号問題を回避し、CPT不変な物理量のモンテカルロシミュレーションが可能になる。複素数のフェルミオン行列式であっても、期待値が実数となることが保証される。

ABSTRACT

An algorithm to simulate full QCD with 3 colours at nonzero chemical potential on the lattice is proposed. The algorithm works for small values of the chemical potential and can be used to extract expectation values of CPT invariant operators.

研究の動機と目的

• 非ゼロの化学ポテンシャルを伴う有限密度格子QCDシミュレーションにおける符号問題に対処すること。
• 複素数のフェルミオン行列式であっても、CPT不変な物理量のモンテカルロシミュレーションを可能にする手法の開発。
• CPT不変なオペレーターの期待値に対して、Re(det D) が有効な実数値ボルツマン重みとして成立することの証明。
• 小スケールのμにおける物理的観測量（プラケット、き裂凝集体、ポリakov線）の数値的解析を可能にすること。
• RHIC や LHC などの重イオン実験に関連するシナリオ（Tが大きく、μが小さい）におけるシミュレーションフレームワークの提供。

提案手法

• 各ゲージ配置 C を別の配置 C^S に写像する離散的変換 S を導入し、経路積分重みが w(C^S) = w*(C) と変換されることを保証する。
• O(C^S) = O(C)、w(C^S) = w*(C)、dU(C^S) = dU(C) という条件を課すことにより、期待値が実数となることを保証する。
• CPT不変なオペレーターに対する正しいボルツマン重みが exp(−S_g) × Re(det D) であることを導出する。ここで det D は有限μにおけるフェルミオン行列式である。
• 変換 S が時間反転と空間反転の特定の組み合わせに対応することを示し、ゲージ場の測度を保存することを確認する。
• 小スケールのμでは Re(det D) が主に正であることが示され、モンテカルロシミュレーションにおける重要度サンプリングが可能になる。
• 今後の課題として、Re(det D) = det Δ を満たす行列Δの特定を提案し、これによりフェルミオンアルゴリズムの高速化が可能になる可能性を示唆する。

実験結果

リサーチクエスチョン

• RQ1 CPT不変な物理量の符号問題は、重みが複素共役となる対の配置を組み合わせることで回避可能か？
• RQ2 小化学ポテンシャルにおける物理的観測量のモンテカルロシミュレーションにおいて、Re(det D) は有効かつ十分なボルツマン重みか？
• RQ3 複素数のフェルミオン行列式であっても、提案手法により期待値の実数性が保たれるか？
• RQ4 本稿で特定された変換以外に、重みの虚数部を対で相殺する他の離散的変換は存在するか？
• RQ5 既存のアルゴリズムを用いて、本手法は実用的に実装可能か？特に det D の評価にかかる計算コストを考慮して。

主な発見

• 有限化学ポテンシャルにおけるCPT不変なオペレーターの正しいボルツマン重みは、複素数の det D ではなく、exp(−S_g) × Re(det D) である。
• 配置 C と C^S を対にして、それらの重みが複素共役となることで、期待値が実数となることが保証される。
• 4⁴のステアガードフェルミオンを用いた数値的実験から、小スケールのμでは Re(det D) が多数の配置で正であることが示された。
• 中程度のμに対しても、正の重みを持つ配置の割合が高く維持されるため、RHIC や LHC の関連領域への応用可能性が裏付けられる。
• det D の明示的計算が必要なこと、および Re(det D) が一部の配置で負値を取るという点で、計算コストと数値的制限が存在する。
• 本手法の対象となる配置クラスにおいて、本稿で特定した変換以外に、重みを複素共役に写像する簡単な離散的変換は見つからず、本手法の独自性が示された。