Lindblad-driven quarkonium production in heavy-ion collisions

技術的概要：重イオン衝突における Lindblad 駆動型クォークニウム生成

問題提起
重いクォークニウム状態は、超相対論的重イオン衝突で形成されるクォーク・グルーオンプラズマ（QGP）の精密プローブとして機能する。その生成は、重いクォークポテンシャルのディバイ遮蔽と媒質内でのコヒーレンス喪失の相互作用によって支配され、結合エネルギーに基づいた階層的な抑制をもたらす。媒質内ポテンシャルが複素数値（遮蔽を記述する実部とランダウ減衰を符号化する虚部）を持つことは知られているが、第一原理から既成状態の抑制と熱平衡化された重いクォークの再結合（凝集）の両方を同時に記述する統一的な理論枠組みは欠如していた。従来のアプローチでは、抑制と再生成が独立した現象論的メカニズムとして扱われることが多かった。

手法
著者らは、QGP バスに結合した重いクォーク・反クォーク（$Q\bar{Q}$）密度行列の進化を記述するために、Lindblad 方程式に基づく開放量子系枠組みを採用する。手法は以下の手順で進行する。

1. 媒質内ポテンシャル: 本研究は、ガウス法則モデル（Lafferty および Rothkopf）を用いて、複素媒質内ポテンシャル $V(r, T) = \text{Re } V + i \text{Im } V$ をパラメータ化する。このモデルは、真空のコーネルポテンシャルとハード・サーマル・ループ（HTL）摂動論を組み合わせ、温度依存性パラメータである単一のディバイ質量 $m_D(T)$ によって支配される。虚部は、弦項における非物理的发散を除去するために厳密に正則化される。
2. スペクトル関数と解離特性: 位置依存の虚数ポテンシャルを用いた時間依存シュレーディンガー方程式を解く（数値的に困難である）代わりに、著者らは周波数空間におけるシュレーディンガー方程式を解く。追跡された遅延グリーン関数の虚部を計算することで、媒質内スペクトル関数を抽出する。
   • 結合エネルギー（$E_n$）: 構成クォークの質量を差し引いた後のスペクトル関数のピーク位置から決定される。
   • 熱的減衰幅（$\Gamma_n$）: ピークの半値全幅（FWHM）から抽出され、連続閾値近傍の非対称性を考慮するために歪んだ Breit-Wigner プロファイルが用いられる。
   • 解離温度（$T_d$）: スペクトルピークが連続体へと融合し（結合エネルギーが消失し）、特定される温度として定義される。
3. 生存確率: ビョルケン膨張を受ける系に対して、既成状態の生存確率は、温度依存減衰幅 $\Gamma(T(t))$ を時間に対して積分することで計算される。
4. 再結合モデル: 枠組みは、Lindblad 方程式の確率的ジャンプ演算子を束縛状態部分空間に射影することで再結合を含むように拡張される。断熱近似（有効ハミルトニアンが量子ジャンプの時間スケールと比較して緩やかに進化する場合）の下で、凝集率が導出される。この率は、自由な重いクォークの数（$N^2$）の二乗と媒質体積（$V$）に比例し、重いクォークが熱平衡に達するのにかかる有限時間を考慮する緩和因子によって変調される。

主な貢献

• 統一的枠組み: 本論文は、$Q\bar{Q}$ 密度行列に対する単一の運動方程式（Lindblad 方程式）から、抑制と再結合の両方のメカニズムを導出する。これにより、減衰を駆動するポテンシャルの虚部と、再結合を駆動するジャンプ演算子の間の内部的整合性が保証される。
• 第一原理からの導出: このアプローチは、独立したモデルのアドホックな組み合わせを回避する。複素ポテンシャルは熱的減衰幅を決定し、それは揺動散逸定理を通じて生存確率と再結合率の正規化の両方を直接決定する。
• 断熱凝集: 著者らは、Lindblad 方程式から特定の凝集モデルを導出し、媒質内での重いクォークの非瞬時的熱平衡化を補正する緩和因子を組み込んでいる。

結果
本研究は、この枠組みを $\sqrt{s_{NN}} = 5.02$ TeV における Pb-Pb 衝突におけるチャームニウム（$J/\psi, \psi(2S)$）およびボトムニウム（$\Upsilon(nS)$）状態に適用する。

• 解離温度:
  • チャームニウム: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV（臨界温度付近または以下で解離）。
  • ボトムニウム: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV。
• 核変換因子（$R_{AA}$）:
  • $J/\psi$: このモデルは、LHC エネルギーにおける高いチャームクォーク多重度によって駆動される中心衝突（$N_{part} \gtrsim 200$）において、有意な再結合寄与を予測する。この再生成は抑制を部分的に補償し、ALICE および CMS データと一致する総 $R_{AA}$ をほぼ 1 に近づける。
  • $\Upsilon(1S)$: ボトムクォークの多重度がチャームのそれよりも約 2 桁小さいため、再結合項は無視できる。$\Upsilon(1S)$ の $R_{AA}$ はほぼ完全に抑制によって支配され、モデル結果は追加の fitting パラメータなしで CMS データと合理的に一致する。

意義と主張
著者らは、この仕事がクォークニウム生成の「第一原理に基づく統一記述」を提供すると主張する。同じ微視的ダイナミクス（複素ポテンシャルおよび Lindblad ジャンプ演算子）から抑制と再結合の両方を導出することにより、この枠組みはこれらの競合する効果の独立した現象論的調整を不要にする。論文は、複素ポテンシャルが（ガウス法則モデルを通じて）固定されれば、クォークニウム収量に対するすべての媒質効果は追加の自由パラメータなしに追随すると主張している。これは、抑制と再生成がしばしば手動で結合される従来の輸送モデルに対する、より理論的に裏付けられた代替案を提供し、重イオン衝突におけるクォークニウムの真の第一原理記述に向けた重要な一歩を表している。