Charmonia at Finite Momentum and Spatial Correlators in Quark-Gluon Plasma

技術的要約：有限運動量におけるチャームオニウムとクォーク・グルーオン・プラズマ中の空間相関関数

問題提起
クォーク・グルーオン・プラズマ（QGP）におけるクォークオニアの空間相関関数は、有限温度格子QCD（lQCD）を用いて高精度に計算することができる。しかし、これらの結果を微視的なスペクトル関数として解釈することは困難であることが証明されている。解釈を妨げる主な課題は以下の2点である：

1. 運動量依存性： 束縛状態の全3運動量（$P$）に対する依存性を信頼性高く計算することは、ポテンシャルモデルにおいては非自明な問題である。なぜなら、ポテンシャルモデルは真空において本質的にローレンツ共変ではないためである。これにより、真の媒質効果を不明瞭にする偽の運動量依存性が生じることがしばしばある。
2. フーリエ変換における振動： 運動量空間から空間相関関数へ移行するために必要なフーリエ変換は、大きな$P$に対して非常に振動的な積分を伴うため、収束が遅く、数値的な不安定性を引き起こす。

近年の光フレーバー部門における進展では、図式的なアンザッツを用いた成果が得られているが、両方の課題を同時に解決するクォークオニウムの空間相関関数に関する微視的モデル計算は欠けていた。

手法
著者らは、ローレンツ共変な散乱フレームワークを構築するために、熱力学的T行列アプローチを採用している。手法の詳細は以下の通りである：

• 真空定式化： 研究は、二体散乱振幅に関する共変な四次元ベテ・サルト（Bethe-Salpeter, BS）方程式から始まる。ローレンツ共変性を維持しつつ相互作用を簡略化するため、分離可能ポテンシャル近似（$V(p', p) = C v(p) v(p')$）を用いている。このポテンシャルは、運動量カットオフ $\Lambda$ によって正則化されている。
• 次元削減： BS方程式は、ブランケンベクレ・スラー（Blankenbecler-Sugar）スキームを用いて三次元へと削減される。これにより、二体プロパゲーターの正しいカット構造が保持され、T行列が不変質量平方（$s = E^2 - \vec{P}^2$）のみに依存することが保証され、偽の運動量依存性が回避される。
• 媒質による修正： システムは、以下の2つの主要な媒質効果を組み込むことでQGP内に埋め込まれる：
  1. 重いクォークの自己エネルギー： 重いクォールおよび反クォークのプロパゲーターは、lQCDの自由エネルギーおよび相関関数によって制約された「強結合シナリオ（SCS）」から導出された自己エネルギー（$\Sigma$）によってドレッシングされる。これらの自己エネルギーには、sQGPの特徴である大きな虚部（衝突幅 $\sim 0.5$ GeV）が含まれている。
  2. ポテンシャルの遮蔽： 相互作用ポテンシャルは、頂点関数へのデバイ型の修正と、温度依存の結合定数 $C(T)$ を通じて遮蔽される。これらのパラメータは、ゼロ運動量におけるモデルのユークリッド相関関数比（ECR）をlQCDデータに一致させることで制約されている。
• スペクトル関数の構成： 媒質中のスペクトル関数はT行列から計算される。極めて重要な点として、著者らは2つの異なる寄与を含めている：
  1. ユニタリティ・カット： 標準的な二粒子（$Q\bar{Q}$）連続体および束縛状態の寄与。
  2. ゼロモード（粒子・ホール励起）： 重いクォークの自己エネルギーにエンコードされた、重いクォークと軽いクォークの散乱（粒子・ホール励起）から生じる寄与。この項は、有限の全運動量 $\vec{P}$ に対してのみ非ゼロとなる。
• 相関関数の評価：
  • ユークリッド相関関数： スペクトル関数のラプラス変換を介して計算される。
  • 空間相関関数： フーリエ変換を介して計算される。高度に振動する被積分体の数値的な課題を克服するため、著者らはスペクトル関数が運動量に依存しないと近似される運動量カットオフ（$P_{cut} = 250$ GeV）を導入しており、この戦略は安定した結果が得られることが検証されている。

主な結果

• T行列の挙動： 全運動量 $P$ が増加するにつれ、T行列における束縛状態のピークはより高い不変質量へとシフトし（これは、遮蔽の強化による結合エネルギーの減少を示唆している）、また、$P \gtrsim 10$ GeV の高運動量では低温において、あるいは低運動量では高温において、ピークの先鋭化（narrowing）が見られる。
• ユークリッド相関関数比（ECR）： モデルは、運動量が増加するにつれて大きなユークリッド時間（$\tau$）でユークリッド相関関数が著しく増大するというlQCDの観測結果を再現することに成功した。著者らは、この増大の主要な要因としてゼロモードの寄与（10〜20%の増加を占める）を特定しており、ユニタリティ・カットの役割は二次的であるとしている。
• 空間相関関数（SCR）： 計算された空間相関関数は、温度の上昇とともに大幅な抑制を示す。
  • ユニタリティ・カットは主要な抑制を駆動しており、大きな空間間隔（$z$）においてより顕著になるものの、緩やかな温度依存性を示す。
  • ゼロモードは、強い逆の温度依存性を示し（ユニタリティ・カットの抑制に対して、低 $z$ で相関関数を増大させる負の寄与を生む）、絶対的な大きさにおいてはユニタリティ・カットよりも劣る（サブリーディングである）。
  • モデルは、アンクエンチドlQCDデータに見られる抑制傾向を定性的に再現しているが、最低温度（$T=194$ MeV）における抑制を過大評価している。

意義と主張
本論文は、ローレンツ共変な運動量依存性とフーリエ変換の数値的困難の両方に同時に対処する、クォークオニアの空間相関関数に関する初の微視的モデル計算を提示していると主張している。

著者らは、彼らのアプローチがlQCDの知見を定性的（場合によっては半定量的）に再現していることを強調している：

1. ゼロモードは、以前のポテンシャルモデルでは捉えることが困難であった、ユークリッド相関関数の運動量依存的な増大をもたらす不可欠なメカニズムとして特定された。
2. 空間相関関数の抑制は、主に媒質中の束縛状態の性質（質量シフトと広がり）およびユニタリティ・カットによって駆動されており、ゼロモードは負の符号と逆の温度依存性を特徴とする二次的な役割を果たしている。
3. 本研究は、明示的な運動量依存性を持つスペクトル関数を与える枠組みを提供することで、lQCDデータと輸送計算の間の溝を埋めるものである。これは、重イオン衝突におけるキネティック理論への必要な入力となる。

著者らは、定量的精度に関しては謙虚な姿勢を保っており、空間相関関数の抑制の正確な起源には、より現実的な相互作用ポテンシャルや更新された重いクォークの自己エネルギーを伴うさらなる調査が必要である可能性を示唆している。彼らは、束縛状態の観測された一様な上方への質量シフトが、観察された挙動の妥当な候補であると考えている。