[論文レビュー] Gauge cooling for the singular-drift problem in the complex Langevin method --- an application to finite density QCD

Keitaro Nagata, Hideo Matsufuru+2 more|arXiv (Cornell University)|Nov 24, 2016

High-Energy Particle Collisions Research参考文献 8被引用数 8

ひとこと要約

本研究では、$4^3\times8$ スピン格子上での有限密度 QCD に対して、ゲージクーリングを用いた複素ランジュバン法を適用し、シミュレーションの安定化とユニタリティノルムの制御を図った。バリオン数密度の onset は $\mu a \approx 0.55$ で発生し、平均場予測と整合的であり、位相クエンチド QCD よりも高い値である。これは、$\mu a \lesssim 0.65$ までにドリフト項の確率分布の指数的減衰基準により妥当性が確認された。特異ドリフト問題の影響が小さいにもかかわらず、ゲージクーリングがシミュレーションを安定化させることを示唆しており、閉じた相における QCD のシミュレーションへの応用の可能性を示している。

ABSTRACT

We study full QCD at finite density and low temperature with light quark mass using the complex Langevin method. Since the singular drift problem turns out to be mild on a $4^3 imes 8$ lattice we use, the gauge cooling is performed only to control the unitarity norm in this exploratory study. We report on our preliminary data obtained from the complex Langevin simulation up to certain Langevin time. While the data are still noisy due to lack of statistics, the onset of the baryon number density seems to occur at larger $μ$ than half the pion mass, which is the value for the phase quenched QCD. The validity of our simulation is tested by the recently proposed criterion based on the probability distribution of the drift term.

研究の動機と目的

• 有限バリオン密度および低温における完全 QCD における複素ランジュバン法の実用可能性を調査すること。特に、閉じた相における応用を目的とする。
• 有限密度 QCD におけるフェルミオン行列の固有値がゼロに近づくことによる特異ドリフト問題に起因する課題を解決すること。
• 最近提案されたドリフト項の確率分布の指数的減衰に基づく基準を用いて、複素ランジュバンシミュレーションの妥当性を検証すること。
• ユニタリティノルムの制御を焦点に、ゲージクーリングが特異ドリフト問題が弱い場合でもシミュレーションを安定化させることを検討すること。

提案手法

• SU(3) ゲージリンク変数の複素化により SL(3,\mathbb{C}) に拡張し、複素ランジュバン方程式による確率的量論化を可能にする。
• ユニタリティノルムの制御のためゲージクーリングを適用するが、本研究では特異ドリフト問題の解消に直接利用しない。
• プラケット作用素とステガードフェルミオン形式を用い、化学ポテンシャル $\mu$ を導入して複素作用を生成する。
• 離散化ステップサイズ $\epsilon$ を用いたランジュバンダイナミクスの数値的シミュレーションを実施し、ガウスノイズと複素化作用から導かれるドリフト項を含む。
• 分割関数のホロモーラル拡張を用いて、バリオン数密度やクリア・コンデンサートなどの物理量を測定する。
• ドリフト項の確率分布が収束のための正しくは指数的以上に減衰する必要があるという基準を用いて、シミュレーション結果の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ゲージクーリングを用いた複素ランジュバン法は、特異ドリフト項の可能性があるにもかかわらず、閉じた相における有限密度 QCD に対して信頼できる結果をもたらすか？
• RQ2 クォーク化学ポテンシャル $\mu$ がどの値でバリオン数密度の onset が発生するか？ また、位相クエンチド QCD と比較するとどうなるか？
• RQ3 有限 $N_s$ の格子 QCD において、最近提案されたドリフト項の確率分布の指数的減衰に基づく基準は有効か？
• RQ4 フェルミオン行列の固有値分布は $\mu$ の増加に伴いどのように変化するか？ また、特異ドリフト問題にどのような含意をもたらすか？
• RQ5 本研究の $4^3\times8$ 格子設定のように特異ドリフト問題が弱い状況において、ゲージクーリングはランジュバンダイナミクスをどの程度安定化させるか？

主な発見

• バリオン数密度は $\mu a \approx 0.4$ から $0.5$ の間で増加を始め、$\mu a = 0.5$ で $\sim 0.01$ に達し、これは格子体積あたり約1つのバリオンに相当する。
• 完全 QCD におけるバリオン数密度の onset は $\mu a \approx 0.55$ で発生し、位相クエンチド QCD の $\mu a \approx 0.265$ よりも高い値であり、平均場予測と整合的である。
• フェルミオン行列の固有値分布は $\mu a = 0.4$ で明確なギャップを示すが、$\mu a = 0.65$ ではそのギャップが薄れ、固有値が原点に近づき、特異ドリフトのリスクが高まっていることが示唆される。
• ドリフト項の確率分布は $\mu a \lesssim 0.65$ で指数的に減衰しており、最近提案された複素ランジュバン法の収束正しさを示す基準を満たしている。
• $\mu a \gtrsim 0.4$ で自己相関時間が著しく増加しており、特異ドリフト問題の発現に起因する動的問題の兆候が示唆される。
• 指数的減衰基準を満たす限り、$\mu a \lesssim 0.65$ までシミュレーションは安定かつ妥当であり、初期ノイズや有限統計の影響があるものの、このパrameter 範囲における手法の実用性が示された。