[論文レビュー] Testing a generalized cooling procedure in the complex Langevin simulation of chiral Random Matrix Theory

Keitaro Nagata, Jun Nishimura+1 more|arXiv (Cornell University)|Nov 27, 2015

Quantum many-body systems参考文献 18被引用数 5

ひとこと要約

本稿は、フェルミオン的特異なドリフトを抑制するため、キラルランダム行列理論（chRMT）における複素ランジュバン法（CLM）に一般化されたゲージクーリング手順を提案する。軽いクォークと近いゼロ固有値を有する状況で生じる特異なドリフトを抑制するために、ディラック作用素に基づくノルムを用いる。この手法により、標準的なCLMが失敗する軽いクォーク領域でも、反復的クーリングによる原点付近の固有値密度の低減によって、正確な結果を再現する。これは、符号問題領域における収束性と安定性の向上を示している。

ABSTRACT

The complex Langevin method has been attracting much attention as a solution to the sign problem since the method was shown to work in finite density QCD in the deconfined phase by using the so-called gauge cooling procedure. Whether it works also in the confined phase with light quarks is still an open question, though. In order to shed light on this question, we apply the method to the chiral Random Matrix Theory, which describes the epsilon regime of finite density QCD. Earlier works reported that a naive implementation of the method fails to reproduce the known exact results and that the problem can be solved by choosing a suitable coordinate. In this work we stick to the naive implementation, and show that a generalized gauge cooling procedure can be used to avoid the problem.

研究の動機と目的

• 低温で軽いクォークを有する有限密度QCDにおけるフェルミオン的特異ドリフト問題に対処すること。
• 軽いクォークにおいて、D+mの固有値がゼロに近い場合にドリフト項が発散するため、通常のCLMが失敗するキラルランダム行列理論（chRMT）における問題を克服すること。
• 座標変換に依存せず、カーテシアン変数を保持する形でシミュレーションを安定化するゲージクーリング手順を開発すること。
• ディラック作用素から導かれるノルムが、力学を効果的に制御し、正確な結果への収束を保証できるかどうかを検証すること。

提案手法

• D+mの固有値分布がゼロに近い部分に敏感な2つの新しいゲージクーリングノルムを導入：$\mathcal{N}_1 = \mathrm{Tr}[(D+m)(D+m)^\dagger]$ および $\mathcal{N}_2 = \mathrm{Tr}[(D+m)^\dagger(D+m)]^{-1}$。
• 各ランジュバンステップの後に、選択したノルムを最小化する複素化された$\mathrm{GL}(N,\mathbb{C})$変換を実行することでクーリング手順を適用する。
• 勾配に基づく最適化により、変換の大きさと方向を決定し、各ステップで10回繰り返すことで、効果的な低減を確保する。
• N=30、$\sqrt{N}\mu=2$、$N_f=2$、$\tilde{m}=mN$を用い、$\Delta t = 5 \times 10^{-5}$、初期化ステップを20,000ステップ破棄するchRMTにこの手法を実装する。
• 初期化後、200ステップごとに観測量（キラル凝集、バリオン数密度）を測定し、正確な解析解と比較する。
• $\mathcal{N}_1$ では $r=0$、$\mathcal{N}_2$ では $r=0.01$ を用いてクーリング強度を制御し、$\xi=300$、$n_{\rm ev}=2$ を $\mathcal{N}_2$ に適用する。

実験結果

リサーチクエスチョン

• RQ1 ディラック作用素の固有値構造に基づく一般化されたゲージクーリング手順は、chRMTにおける軽いクォーク領域における複素ランジュバン法の安定化に有効か？
• RQ2 $\mathcal{N}_1$ や $\mathcal{N}_2$ ノルムを用いたクーリングは、D+mの固有値がゼロに近い場合に生じる危険なドリフト項を抑制できるか？
• RQ3 新規クーリング手順を用いたCLMは、$mN$が小さい領域においてキラル凝集とバリオン数密度の正確な結果を再現できるか？
• RQ4 クーリングノルムの選択が、複素ランジュバンシミュレーションにおけるD+mの固有値分布にどのように影響するか？
• RQ5 座標再パラメータ化を要件としない形で、この手法を有限密度QCDに一般化できるか？

主な発見

• 通常の複素ランジュバン法（CLM）は、$\tilde{m} = mN \lesssim 10$ の場合に、キラル凝集とバリオン数密度において正確な結果から顕著なずれを示し、失敗する。
• 新しいゲージクーリング（$\mathcal{N}_1$ および $\mathcal{N}_2$ ノルムを用いる）により、CLMは$\tilde{m}$の全範囲で正確な結果を成功裏に再現することができ、特に軽いクォーク領域でも有効である。
• 両クーリング手順において、D+mの原点付近の固有値分布が顕著に抑制される。$\mathcal{N}_1$ は実軸方向に分布を狭めるのに対し、$\mathcal{N}_2$ は固有値を原点から反発する。
• ホロモーラフィックな観測量たるキラル凝集は、クーリングノルムの種別にかかわらず正確な値に収束するため、適切に安定化された場合の方法の一貫性と普遍性が示された。
• $\tilde{m} = 3$ においても、通常のCLMが完全に失敗する状況でも、本手法は安定かつ正確に保たれるため、CLMの適用範囲が著しく拡張された。
• クーリング手順により、動的変数の虚数方向への長距離な逸脱が効果的に制御され、複素ランジュバンの正当化の破綻を防いでいる。