Event-by-event fluctuations of elliptic flow in ultrarelativistic O+O collisions

技術要約：超相対論的O+O衝突における楕円流のイベントごとのゆらぎ

問題提起
本論文は、$\sqrt{s_{NN}} = 5.36$ TeVにおける酸素・酸素（O+O）衝突における異方的流、特に楕円流（$v_2$）の性質を調査している。大きな原子核（例：Pb+Pb）における支配的な見解では、楕円流は衝撃パラメータによって決定される重なり領域の幾何学的非対称性に由来するが、著者らはO+O衝突がこれとは異なる領域を代表している可能性を提示している。小規模な系においては、核子の数が限られているため、エネルギー堆積のイベントごとのゆらぎが初期幾何学を支配し、その結果、O+O衝突は、フローの起源に関して重イオン衝突よりも小規模系（p+p）に類似した性質を持つ可能性がある。本研究は、これらのゆらぎを定量化し、削減されたゆらぎモードのセットが楕円流の統計的性質を正確に再現できるかどうかを判断することを目的としている。

手法
著者らは、McDipper（初期状態）、MUSIC（流体力学発展）、およびiSS（粒子化）を組み合わせた、完全な3次元数値フレームワークを用いている。

• 初期状態： $^{16}$O核の核子構成は、核格子有効場理論（NLEFT）を用いて生成される。これらの構成はランダムに回転され、ランダムな衝撃パラメータと共に衝突させられる。Mcディッパー（McDipper）コードは、カラーグラス凝縮有効場理論（GBWモード）に基づいており、サブ核子レベルのホットスポットと厚さのゆらぎを含む3次元エネルギー密度プロファイル $\Phi(x, y, \eta_s)$ を生成する。
• 流体力学： これらのプロファイルは、一定の比剪断粘性（$\eta/s = 0.16$）と温度依存のバルク粘性を持ち、凍結温度155 MeVまで進化するMUSICの初期条件として用いられる。
• 分解フレームワーク： 本手法の主要な革新は、3次元初期状態のエネルギー密度を、平均プロファイル $\bar{\Psi}$ と、係数 $c_l$ を持つ一連の無相関なゆらぎモード $\{\Psi_l\}$ へと分解することである。これは、エネルギー密度のゆらぎの自己相関行列の固有値分解を通じて行われる。このアプローチは、非ブースト不変な3次元初期状態に対して初めて適用された。
• 線形応答： 著者らは、観測量（非心度 $\varepsilon_2$ およびフロー $v_2$）をゆらぎモードの線形関数として近似する線形応答フレームワークを利用している。彼らは、「イベントごと」のシミュレーション（中心度クラスごとに5000イベント）と、減少したモードのガウス分布に従う係数をサンプリングして生成された「線形化された」観測量とを比較している。

主な貢献と結果

1. ゆらぎの支配性： 本研究は、O+O衝突において、初期の楕円非心度（$\varepsilon_2$）が衝撃パラメータによる有限の幾何学的形状ではなく、主にエネルギー堆積のイベントごとのゆらぎによって駆動されていることを示している。その結果、非中心的な衝突（20–30% 中心度）においても、平均的な初期状態 $\bar{\Psi}$ はほぼ円形であり、平均的な非心度は非常に小さい（$\varepsilon_2 \approx 0.087$）。これは、平均的な状態が顕著にアーモンド型を示す重イオン衝突とは対照的である。
2. ゆらぎの重み： O+Oにおける平均状態の相対的な重み（$\bar{w}$）はわずか $\sim 2\%$ であり、Pb+Pbで観察される $\sim 13\%$ よりも大幅に低いことが判明した。これは、O+Oにおける初期幾何学へのゆらぎの寄与が支配的であることを定量的に裏付けている。
3. モード分解の精度： 著者らは、少数のゆらぎモード（具体的には10個のモード）を用いることで、$\varepsilon_2$ と $v_2$ の同時確率分布を妥当な精度で再現できることを示している。「線形モデル」および「解析的ガウスゆらぎ近似」（G.f.a.）は、$\varepsilon_2 \lesssim 0.6$ において、フルイベントごとのシミュレーションと極めて近い一致を示す。
4. 非線形性の限界： 大きな非心度（$\varepsilon_2 > 0.6$）において、線形近似は非線形流体応答（具体的には $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$ という関係における3次の項）により、イベントごとのデータから逸脱する。線形モデルは、分布の「バナナ型」の形状や、これらの非線形性によって引き起こされる $v_2$ のより広い裾野を捉えることができない。
5. 縦方向のデコレーション： 本論文は、積分されたフローは少数のモードで捉えられる一方で、非心度の縦方向のデコレーション（$r_2$）のような縦方向の構造に敏感な観測量は、イベントごとの結果を再現するためにはるかに多くのモード（数百個）を必要とすることを指摘している。これは、モードスペクトルの高次に出現するラピディティ奇モードの重要性に起因している。

意義と主張
本論文は、O+O衝突が単なる重イオン衝突の「スケールダウン版」ではないと主張している。むしろ、O+Oにおける楕円流の起源は根本的に異なり、グローバルな衝撃パラメータの幾何学ではなく、初期状態のゆらぎによって支配されている。このことは、O+O衝突が、ゆらぎが主要な駆動源となる集団現象の発現を研究するためのユニークな実験室を提供することを示唆している。

著者らは、自らのモードごとの分解フレームワークが、楕円流のバルク統計的挙動を記述するための、フルイベントごとの3次元流体シミュレーションに代わる、コストの低い実行可能な代替案であると控えめに主張している。彼らは、限定されたセットのゆらぎモードが、平均的な幾何学に対するゆらぎの支配性を考慮すれば、$v_2$ と $\varepsilon_2$ の平均値および確率分布を再現するのに十分であることを実証した。しかし、線形近似には高非心度における非線形応答に関する限界があること、また、より複雑な観測量（縦方向のデコレーションなど）には、より大規模なモード基底が必要であることを認めている。