[論文レビュー] Finite Density Simulations with Canonical Ensemble

Anyi Li, Xiangfei Meng+2 more|ArXiv.org|Oct 14, 2008

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

ひとこと要約

この論文は、$6^3\times4$ スピン格子上でのウィルソンフェルミオンを用いた有限密度QCDのシミュレーションに正準集合法を適用し、4フレーバーの場合に明確な一次相転移を示す'S字型'の化学ポテンシャル対バリオン数の挙動を特定した。2フレーバー系では$0.83T_c$までそのような信号は観測されず、3フレーバー系の予備結果でも$T=0.92T_c$において一次相転移の兆候は認められない。計算上の課題があるものの、この手法はQCD相図の探索に実用可能であることが示された。

ABSTRACT

QCD at non-zero baryon density is expected to have a critical point where the zero-density cross-over turns into a first order phase transition. To identify this point we scan the density-temperature space using a canonical ensemble method. For a given temperature, we plot the chemical potential as a function of density looking for an "S-shape" as a signal for a first order transition. We carried out simulations using Wilson fermions with $m_π\approx 1{GeV}$ on $6^3 imes 4$ lattices. As a benchmark, we ran four flavors simulations where we observe a clear signal. In the two flavors case we do not see any signal for temperatures as low as $0.83 T_c$. Preliminary results for the three flavor case are also presented.

研究の動機と目的

• 非ゼロバリオン密度におけるQCDの一次相転移の存在を正準集合アプローチを用いて調査すること。
• ウィルソンフェルミオンと高止まりピオン質量（$m_\pi \approx 1$ GeV）を用いた有限密度格子QCDシミュレーションの実現可能性を検証すること。
• 温度とバリオン密度をスキャンすることで、QCD相図における臨界点の探索を行うこと。
• 2、3、4フレーバーの異なる数のダイナミカルフレーバー間で比較し、物理的点への傾向を評価すること。
• 正準集合法が、符号問題に悩まされるグランドカノニカルアプローチの代替手段として有効であるかを検証すること。

提案手法

• 正準集合は、グランドカノニカル分配関数のフーリエ変換により構成され、$ Z_C(V,T,k) = \frac{1}{2\pi} \int_0^{2\pi} d\phi \, e^{-ik\phi} Z(V,T,\mu)|_{\mu = i\phi T} $ という関係式を用いる。
• 固定されたネットクォーク数$k$を扱うために、$N$点での離散的近似$\phi_j = 2\pi j / N$を用いたフーリエ変換の近似により、射影された行列式$ \text{det}_k M^2(U) $を計算する。
• 計算の効率性と精度を向上させるために、直接的なフーリエ和算の代わりに巻き込み数の展開を用いる。これにより、計算コストと数値的不安定性が低減される。
• 化学ポテンシャルは$ \langle \mu \rangle_{n_B} = -\frac{1}{\beta} \ln \left( \frac{\langle \gamma(U) \rangle_o}{\langle \alpha(U) \rangle_o} \right) $ として抽出され、ここで$ \alpha(U) $ と $ \gamma(U) $ は位相重み付き行列式である。
• 共存領域の境界と臨界点の特定のため、'S字型'領域にマクスウェルの構成法を適用する。
• 3フレーバー系のシミュレーションでは、離散化誤差を低減するため、Iwasakiゲージ作用素とクローバーフェルミオン作用素を用いるが、2フレーバーおよび4フレーバー系ではウィルソンフェルミオンを用いる。

実験結果

リサーチクエスチョン

• RQ1 $ m_\pi \approx 1 $ GeV で、$6^3 \times 4$ 格子上での4つの縮重フレーバーに対して、正準集合法が有限密度QCDにおける一次相転移を検出できるか？
• RQ2 2つの縮重軽クォークを有するQCD相図に臨界点が存在するか、$0.83T_c$の温度までその存在が同定可能か？また、この手法で検出可能か？
• RQ3 3フレーバー系の相構造は、4フレーバー系および2フレーバー系と比較してどのように異なり、$ T = 0.92T_c $ で一次相転移の兆候を示すか？
• RQ4 フェルミオン離散化（ウィルソン対クローバー）およびピオン質量に依存する程度はどの程度か？
• RQ5 S字型のフィッティング手法の違いに対して、抽出された相境界および臨界点の位置はどの程度頑健か？

主な発見

• 4つの縮重フレーバーの場合、$ T_c $ より低い温度で化学ポテンシャル対バリオン数プロットに明確な'S字型'構造が観測され、一次相転移を示唆する。
• マクスウェルの構成法により共存領域の境界と臨界化学ポテンシャルが決定され、以前のステガードフェルミオン研究と一貫した結果を得たが、本研究では共存領域がより狭い。
• 2フレーバー系では、符号問題が深刻でないにもかかわらず、$ T = 0.86T_c $ および $ T = 0.83T_c $ において一次相転移の兆候は認められない。
• 3フレーバー系の予備結果（$ T = 0.92T_c $）では'S字型'構造が観測されず、この温度で一次相転移の兆候は認められない。
• 実行列式の絶対値と位相を観測量に組み込むことで、符号問題を回避し、安定なシミュレーションが可能となった。
• 共存領域幅の差異は、ピオン質量の重さやフェルミオン形式の違いに起因する可能性が示唆された。