[論文レビュー] QCD equation of state via the complex Langevin method

Felipe Attanasio, Benjamin Jäger+1 more|arXiv (Cornell University)|Mar 24, 2022

Quantum Chromodynamics and Particle Interactions被引用数 5

ひとこと要約

本稿では、符号問題を克服するために複素ランジュバン法を用いて、有限のバリオン化学ポテンシャルにおける2フラバーQCDの状態方程式のab initio格子QCDシミュレーションを提示する。圧力、エネルギー、エントロピーの状態方程式をT–μ平面で報告し、低温・低密度領域におけるシルバー・ブレーズ現象の数値的証拠を示す。

ABSTRACT

We present lattice simulations on the phase diagram of Quantum Chromodynamics (QCD) with two light quark flavours at finite chemical potential $μ$. To circumvent the sign problem we use the complex Langevin method. In this study, we have carried out finite density lattice computations for a pion mass of $\sim 480$ MeV. We report on the pressure, energy and entropy equations of state in ab-initio lattice QCD calculations, as well as the observation of the Silver Blaze phenomenon.

研究の動機と目的

• 有限のバリオン化学ポテンシャルμおよび温度TにおけるQCD状態方程式を計算すること、特に符号問題が顕著になる低温・高密度領域を対象とする。
• 動的なクォークを伴う格子QCDシミュレーションにおける符号問題に対処するため、ゲージ場をSL(3,C)に拡張する複素ランジュバン法を採用すること。
• パイオン質量が ≲ 480 MeVである範囲において、T–μ平面全域における熱力学的観測量（圧力、エネルギー密度、エントロピー密度）の振る舞いを調査すること。
• 低温・低密度領域において、T=0でμ < m_N/3のときクォーク数密度がゼロのままとなるシルバー・ブレーズ現象の数値的証拠を提供すること。
• 中性子星の安定性に関連する低温における状態方程式の硬さ（stiffness）を調査すること。

提案手法

• 複素ランジュバン法を用いて、有限の化学ポテンシャルにおける2つの縮退した軽いクォークを有するQCDをシミュレートし、ゲージ場を非コンパactsな群SL(3,C)に拡張することで符号問題を回避する。
• シミュレーションにはウィルソンプラケットゲージ作用素と、時間反転対称性を破る化学ポテンシャル項を含む非改善ウィルソンフェルミオン作用素を用いる。
• 適応的ステップサイズを用いた確率的時間発展とゲージ冷却を組み合わせ、複素化された多様体上でのダイナミクスを安定化させる。
• SU(3)部分群に近い位置を維持するための動的安定化項を適用し、手法の収束性と正しさを保証する。
• 圧力のμ関数としての多項式フィットに数値微分を施し、熱力学的観測量を計算する。
• フィットされた圧力から数値的にトレース異常を計算し、熱力学的関係式を用いてエネルギー密度およびエントロピー密度を導出する。

実験結果

リサーチクエスチョン

• RQ1 低温・有限密度領域におけるQCD状態方程式（特に圧力、エネルギー密度、エントロピー密度）は、温度とバリオン化学ポテンシャルにどのように依存するか？
• RQ2 符号問題が顕著な状況下でも、複素ランジュバン法がQCD状態方程式に対して信頼できる結果をもたらすか？
• RQ3 低温および小化学ポテンシャル領域における格子QCDシミュレーションで、シルバー・ブレーズ現象の数値的証拠を観測できるか？
• RQ4 状態方程式の硬さは温度にどのように依存し、中性子星の安定性にどのような示唆を与えるか？
• RQ5 固定温度Tにおけるボトン数密度および圧力はμの関数としてどのように振る舞うか、特に中性子質量閾値近傍で注目すべき点は何か？

主な発見

• 圧力はμに伴い増加し、低温ではμが中性子質量m_Nを超えた以降に顕著な増加を示す。これはバリオン物質の始発に一致する。
• エネルギー密度およびエントロピー密度はμ ≳ m_Nで明確な始発を示し、エネルギー密度のゼロからの乖離がμとともに徐々に増大するため、熱力学的活動性の増大が示唆される。
• 低温領域では状態方程式が硬くなる。これはμに伴うボトン数密度の上昇が遅やかであることに起因し、中性子星の安定性要件と整合的である。
• シルバー・ブレーズ現象の数値的証拠が観測された：T = 50 MeVでμ < m_N/3の範囲ではクォーク数密度はほぼゼロのまま維持され、低温における期待される振る舞いが確認された。
• シミュレーションではT ≈ 50 MeVでバリオン数密度が約15n₀、T ≈ 200 MeVで約20n₀に達しており、パイオン質量が約480 MeVであるため、高密度領域へのアクセスが可能である。
• 各固定温度における圧力状態方程式は、μに関して立方多項式フィットでよく記述され、導出観測量のための高精度な数値微分が可能である。