[論文レビュー] Exploring Complex-Langevin Methods for Finite-Density QCD

D. K. Sinclair, J. B. Kogut|arXiv (Cornell University)|Oct 21, 2015

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

ひとこと要約

本稿は、2つの軽いクォーカーを含む12⁴格子上で、有限密度QCDをゼロ温度でシミュレートするために、ゲージクーリングを用いた複素Langevinダイナミクスの応用を調査する。ランアウェイ行動の課題にもかかわらず、適応的更新とゲージクーリングによりアルゴリズムは安定化し、定性的に正しい結果が得られた：μ < m_N/3ではクォーク数密度はゼロから上昇し、飽和時にクリティカルコンデンゲートはゼロへ減少し、プラケットはブリッジド値へ上昇するが、平衡化の問題により核物質転移付近では完全な収束に至らない。

ABSTRACT

QCD at non-zero chemical potential ($μ$) for quark number has a complex fermion determinant and thus standard simulation methods for lattice QCD cannot be applied. We therefore simulate this theory using the Complex-Langevin algorithm with Gauge Cooling in addition to adaptive methods, to prevent runaway behaviour. Simulations are performed at zero temperature on a $12^4$ lattice with 2 quarks which are light enough that $m_N/3$ is significantly larger than $m_π/2$. Preliminary results are qualitatively as expected. The quark-number density is close to zero for $μ&lt; m_N/3$, beyond which it increases, eventually reaching its saturation value of $3$ for $μ$ sufficiently large. The chiral condensate decreases as $μ$ is increased approaching zero at saturation, while the plaquette increases towards its quenched value. We have yet to observe the transition to nuclear matter at $μ\approx m_N/3$, presumably because the runs for $μ$ between $m_N/3$ and saturation have yet to equilibrate.

研究の動機と目的

• 複素Langevinダイナミクスにゲージクーリングを適用した手法が、ゼロ温度における有限密度QCDに対して信頼できる結果をもたらすかどうかを検証すること。
• 化学ポテンシャルμの関数としてのクォーク数密度、クリティカルコンデンゲート、プラケットの挙動を調査すること。
• μ ≈ m_N/3で核物質への転移を捉えられるかどうかを特定すること。
• 複素作用と非コンパクトゲージ場を伴う状況下で、複素Langevinアプローチの収束性と正しさを評価すること。
• より軽いクォーク質量、より大きな格子、有限温度領域へのこの手法の拡張可能性を検討すること。

提案手法

• QCDに有限μを適用する際、ゲージ場をSU(3)からSL(3,C)に複素化することで、複素フェルミオン行列式のシミュレーションが可能になる。
• 確率的時間発展は、実ガウスノイズを用いたLangevin方程式に従い、数値的安定性を高めるために部分的2次精度積分器が用いられる。
• ゲージクーリングは、ユニタリティノルム F(U) = (1/4V)ΣTr[U†U + (U†U)⁻¹ - 2] を繰り返し最小化することで実装され、発散するゲージ場の動きが抑制される。
• 適応的時間ステップ法によりドリフト力の制御が行われ、力が増大するとdtを小さくすることで安定性が向上する。
• 各更新後にクォーク数密度、クリティカルコンデンゲート、プラケットといった観測量が測定され、90,000～500,000スイープにわたりデータ収集が行われる。
• シミュレーションは、μ ≈ m_N/3近辺での準安定状態と平衡化の挙動を調査するため、SU(3)多様体上に初期化されたゲージリンクを用いる。

実験結果

リサーチクエスチョン

• RQ1 ゲージクーリングを用いた複素Langevinは、ゼロ温度における有限密度QCDに対して安定的かつ収束的結果をもたらすか？
• RQ2 μ ≈ m_N/3を超えると、クォーク数密度が期待通りに増加するか、核物質の始まりを示唆するか？
• RQ3 クリティカルコンデンゲートはμの増加に伴い単調に減少し、飽和時にゼロに近づくか？これはクォーク- gluonプラズマ相でのスピン対称性の回復と整合的である。
• RQ4 μ ≈ m_N/3で一次相転移の兆候が認められ、遅い平衡化により準安定状態に閉じ込められる状態が続くか？
• RQ5 大μ領域で2フラーバーの色超伝導体の形成が検出可能か？また、リウェイト法や重クォーク法との予測と比較してどうか？

主な発見

• μ < m_N/3ではクォーク数密度はゼロ付近を維持し、その後μが大きくなると飽和値3へ上昇し、核物質の始まりが示唆される。
• クリティカルコンデンゲートはμの増加に伴い単調に減少し、飽和時にゼロに近づく。これはクリティカル対称性の回復と整合的である。
• プラケットはμの増加に伴い上昇し、ブリッジド値に近づく。これは高密度媒体におけるゲージ場強度の上昇を示唆する。
• μ > m_N/3では平衡化が遅く、一次相転移に起因する準安定性が原因で、現在の計算では完全な収束に至らない。
• 適応的時間ステップ法とゲージクーリングによりアルゴリズムは安定化したが、冷却を施してもユニタリティノルムがSU(3)からわずかにずれる傾向があり、残留不安定性を示唆する。
• 予備的な結果は定性的に正しいが、特にμ ≈ m_N/3近辺で完全な平衡化が得られていないため、誤差バーはまだ入手できない。