Fluctuations of $v_2$ versus fluctuations of $\varepsilon_2$: Beyond linear and cubic response

技術要約：$v_2$ のゆらぎ対 $\varepsilon_2$ のゆらぎ：線形および 3 次応答を超えて

問題提起
楕円流（$v_2$）は、超相対論的重イオン衝突における集団的挙動の主要なシグネチャーであり、クォーク・グルオン・プラズマの流体動力学的膨張に由来する。伝統的に、$v_2$ のイベントごとのゆらぎは、初期空間異方性（$\varepsilon_2$）に対する線形応答としてモデル化され、非線形効果に対する現象論的な補正は、$(\varepsilon_2)^3$ に比例する 3 次項に限定されることが多い。この線形＋3 次の枠組みは、初期状態モデル（例：TRENTo モデル）を制約する上で成功してきたが、著者らは、既存の応答理論では十分に捉えきれない偏差が存在することを特定している。具体的には、初期離心率と最終的な楕円流の確率分布の関係において、標準的な線形＋3 次近似を覆す非対称性と応答の大きさが見られる。

手法
本研究では、MUSIC コードを用い、$\sqrt{s_{NN}} = 5.02$ TeV における Pb+Pb 衝突の広範なイベントごとの流体動力学シミュレーションを採用している。決定的な手法上の特徴は、衝撃パラメータ（$\vec{b}$）の大きさと方向の両方を固定し、「固有フレーム」を定義している点である。このアプローチにより、衝撃パラメータにおける古典的なゆらぎから、量子的なゆらぎ（核子波動関数に由来するもの）を分離することが可能となり、$\vec{b}$ がゆらぐ標準的な中心度クラスのシミュレーションと比較して、初期状態と最終状態をより精密に比較できる。

解析は以下のステップで行われる：

1. シミュレーション： 特定の衝撃パラメータ（例：$b=10$ fm、これは 40% 中心度に相当）に対して、初期条件に TRENTo を、粘性流体動力学発展に特定の剪断粘性（$\eta/s = 0.08$）および体積粘性パラメータリングを用いて、約 $4 \times 10^5$ イベントをシミュレートする。
2. 累積量展開： 著者らは、固有フレーム内での累積量（$c_{mp}$）を用いて $\varepsilon_2$ と $v_2$ の確率分布を展開し、$\varepsilon_2$ と $v_2$ のゆらぎを分析する。これにより、ガウス型ゆらぎと非ガウス型の補正を体系的に分類できる。
3. 応答モデリング： 線形応答（$v_2 = \kappa \varepsilon_2$）および線形＋3 次応答（$v_2 = (\kappa_2 + \kappa'_2 |\varepsilon_2|^2)\varepsilon_2$）の妥当性を、フル流体動力学シミュレーションデータと散布図および統計モーメントを比較することで検証する。

主な結果
解析の結果、標準的な線形＋3 次応答シナリオからのいくつかの重要な逸脱が明らかになった：

• 逆の非対称性： $v_2$ のゆらぎの非対称性は、$\varepsilon_2$ のゆらぎの非対称性と逆であることが判明した。初期の離心率分布は負の歪度（スキューネス）を示すが、結果としての $v_2$ 分布は逆の非対称性を示しており、これは 3 次項のみでは再現できない特徴である。
• ゆらぎへの増強された応答： ゆらぎに対する流体動力学的応答は、平均的な幾何学的形状に対する応答よりも約 10～20% 大きい。この効果は、平均と分散を主に調整するものの、ゆらぎへの応答の全容を捉えるには至らない 3 次応答項では説明できない。
• 分布の形状： 反応平面上に投影された $\varepsilon_2$ と $v_2$ の同時確率分布は、3 次応答と一致する上向きの曲率を持つ「バナナ型」の高密度領域を示す。しかし、フル流体動力学シミュレーションは、この線に対して有意な広がりを示す一方で、線形＋3 次モデルは、（微小な $\varepsilon_{2y}$ のゆらぎを除いて）実質的に広がりを持たない、非常に狭い分布を予測する。
• 等方性： 初期離心率のゆらぎが（反応平面に対して狭い）異方的であるのに対し、$v_2$ のゆらぎはおよそ等方的であることが観察された。

意義および主張
本論文は、標準的な線形＋3 次応答モデルは、初期条件から最終的な観測量への流体動力学的写像を完全に記述するには不十分であると主張している。特定された不一致（特に、ゆらぎへの応答の増強と逆の非対称性）は、より高次またはより複雑な非線形メカニズムが作用していることを示唆している。

著者らは、これらの知見が実験データの解釈に直接的な影響を与えることを強調している。特に、ゆらぎに対する応答の未説明の増強は、2 次（$v_2\{2\}$）と 4 次（$v_2\{4\}$）のフロー・累積量の分裂を拡大させる。したがって、線形＋3 次の近似のみに依存することは、初期状態モデルの制約を不正確にする可能性がある。固有フレームにおける累積量展開を利用することで、本研究は初期条件と観測量を結びつけるためのより厳密な枠組みを提供し、既存の LHC データを用いて初期状態モデルのパラメータをより精密かつ正確に制約することを可能にする。