Effects of the centrality determination method for the equation of state and nucleonic observables from Au+Au collisions at $\sqrt{s_{NN}}$ = 2.4 GeV

本研究は、$\sqrt{s_{NN}}$ = 2.4 GeVにおけるAu+Au衝突のUrQMDシミュレーションを用いて、中心度決定手法の選択が核物質状態方程式の抽出に著しい影響を与えることを示し、Glauberモンテカルロに基づくアプローチは状態方程式そのものよりも大きな不確かさを導入する一方で、幾何学的解釈は動力学的な多重度選択と一致していることを明らかにしている。

要約

あなたは、非常に高密度で高温のスープ（核物質）のレシピを理解しようとしているところだと想像してください。その方法は、2つの巨大なボウルが互いに激突する様子を観察することです。物理学者はこれを「重イオン衝突」と呼びます。これらのボウル（金原子）が衝突すると、中の物質が押しつぶされ、中性子星の核や初期宇宙に見られるような条件が作り出されます。

この論文の目的は、この「核のスープ」がどれくらい「硬い」か、あるいは「柔らかい」かを解明することです。この性質は**状態方程式（EoS）**と呼ばれます。これは、スープがゼリー（柔らかい）なのか、それともコンクリート（硬い）なのかを問うようなものです。

問題点：どのように「衝撃」を測定するか？

実際の実験では、衝突の中心を直接見ることはできません。それはあまりにも速く、あまりにも小さすぎるからです。そのため、科学者たちは、飛び出してきた「破片（粒子）」の数を数えることで、衝突がどれほど激しかったのかを推測しなければなりません。

この論文の著者たちは、次のような単純な問いを投げかけました。「衝撃の激しさを推測する方法は重要だろうか？ その推測によって、スープがゼリーかコンクリートかという結論が変わってしまうのだろうか？」

彼らは、衝突の「激しさ（中心度）」を推測するために、3つの異なる方法をテストしました。

1. 破片カウント (Mch): 「より多くの粒子が見えるほど、より激しい衝突だったと仮定する。」
2. 幾何学的推測 (bf): 「ボールがどのように重なっているかという単純な図に基づき、衝突の激しさを判断する。」
3. コンピュータモデルによる推測 (br): 「有名なコンピュータ・シミュレーション（グローバー・モデル）を用いて、関与した人数（粒子）に基づいて衝突の激しさを推測する。」

実験：仮想衝突テスト

研究者たちは、実際の原子を衝突させたわけではありません。彼らはUrQMDと呼ばれるスーパーコンピュータのシミュレーションを使用して、特定のエネルギー（2.4 GeV）で金原子を衝突させました。彼らは衝突を2回実行しました。

• 一度は「ゼリー・スープ（軟らかいEoS）」を用いて。
• もう一度は「コンクリート・スープ（硬いEoS）」を用いて。

その後、上述の3つの異なる推測方法を用いて、結果を検証しました。

調査結果：「推測ゲーム」が重要である

彼らが発見したことを、簡単な比喩を用いて説明します。

1. 「破片カウント」対「幾何学的推測」 (Mch vs. bf)
破片カウントを使用し、それを単純な幾何学的推測と比較したとき、結果は驚くほど似ていました。

• 比喩: 車の衝突がどれほど激しかったかを、壊れたヘッドライトの数を数えて推測することを想像してください。ヘッドライトを数えれば、衝突の深刻さについてかなり良い推測ができ、それは車の重なり具合を示す単純な図ともよく一致します。
• 結果: 「推測」における不確かさは、スープの硬さの測定を台無しにしませんでした。「ゼリー」と「コンクリート」の結果は明確に区別されたままでした。

2. 「破片カウント」対「コンピュータモデル」 (Mch vs. br)
ここから事態は混乱してきました。破破片カウントを使用し、それを**コンピュータモデル（グローバー）**と比較したとき、結果は大きく異なりました。

• 比喩: 衝突の深刻さを推測するために、ハイテクな天気アプリを使うことを想像してください。そのアプリは、「部品の破片が多い＝激しい衝突」という前提を持っていますが、それは大きな高速道路での衝突には機能しても、小さな低速の接触事故には機能しません。
• 結果: この特定のエネルギーレベル（2.4 GeV）において、コンピュータモデルは間違っていました。それは衝突の激しさを誤認していました。この誤った推測のために、「ゼリー・スープ」と「コンクリート・スープ」の違いが曖昧になってしまいました。実際、間違った推測方法によって生じた誤差は、ゼリーとコンクリートの実際の差よりも大きかったのです！

主な教訓

この論文は、核物質の「硬さ」を正確に測定したいのであれば、次のように結論付けています。

• 古いコンピュータモデル（グローバー）を信じてはいけません。 それは、都市の地図を使って森の中をナビゲートするようなものです。ルールが適用されません。
• 幾何学的な論理、または直接的な破片カウントに従ってください。 これらの手法は、このエネルギーレベルにおいては現実とより良く一致しています。

要約すると： もし衝突を測るための間違った定規を使えば、実際には「コンクリート」であるのに「ゼリー」だと判断したり、あるいはその逆になったりする可能性があります。衝突を定義する方法は、あなたが測定しようとしている物理現象と同じくらい重要です。著者たちは、レシピを正しく導き出すためには、「破片カウント」を「実際の衝突の幾何学」と一致させる一貫した方法が必要であると警告しています。

技術概要

技術要約：$\sqrt{s_{NN}} = 2.4$ GeV における Au+Au 衝突における核子観測量への中心度決定の影響

問題提起
中間エネルギー重イオン衝突（HIC）の解析において、衝突の中心度（centrality）の決定は主要な系統誤差の要因であり続けている。これは、超飽和密度における核物質状態方程式（EoS）を調査する際に特に重要である。理論的なシミュレーションでは、衝突パラメータ（$b$）が直接的な入力として利用される一方で、実験的な中心度は、理論モデルを用いて最終状態の観測量（例：荷電粒子多重度、$M_{ch}$）から推測される。衝突パラメータとこれらの観測量との間には、厳密な一対一の対応関係は存在しない。先行研究では、選択されたイベントの真の衝突パラメータ分布が持つ幅が、集団流（collective flow）に大きな影響を与えることが示されている。さらに、低エネルギー（例：2.4 GeV）におけるグローバー・モンテカルロ（MC）モデルのような標準的な中心度決定手法の適用可能性は疑わしい。なぜなら、このモデルは粒子多重度が参加核子数に比例することを前提としているが、この仮定はより高エネルギーにおいてより妥当なものであるからである。本研究は、異なる中心度決定手法に関連する不確実性を定量的に評価し、それらが最終状態の EoS 感応的観測量に与える影響を評価することを目的とする。

手法
本研究では、$\sqrt{s_{NN}} = 2.4$ GeV（ビームエネルギー 1.23A GeV に相当）における Au+Au 衝突をシミュレートするために、Ultra-relativistic Quantum Molecular Dynamics (UrQMD) モデルを使用する。

• 状態方程式 (EoS): 2 つの運動量依存型 EoS シナリオを採用する。ソフトな EoS（$K_0 = 200$ MeV、SM と表記）およびハードな EoS（$K_0 = 380$ MeV、HM と表記）。
• 中心度決定法: 特定の中心度クラス（0–10%、10–20%、20–30%、30–40%）内のイベントサンプルを選択するために、3 つの異なる手法を比較する：
  1. $M_{ch}$: 全ての荷電粒子の多重度に基づく選択。
  2. $b_f$: 入力された衝突パラメータの範囲に基づく幾何学的な解釈。
  3. $b_r$: 多重度を衝突パラメータの範囲にマッピングする、グローバー MC モデル（文献 [7] に基づく）から導出された手法。
• 観測量: 本研究では、自由陽子の収量およびラピディティ分布、直接流（$v_1$）および楕円流（$v_2$）係数（具体的には傾き $v_{11}$ および中間ラピディティ値 $v_{20}$）、横運動量（$p_T$）分布、および平均横運動量（$\langle p_T \rangle$）を解析する。

主な結果

• 衝突パラメータ分布: $M_{ch}$、$b_f$、および $b_r$ によって選択されたイベントサンプルの真の衝突パラメータ分布の間には、顕著な相違が存在する。$M_{ch}$ と $b_f$ は比較的整合性のある分布を示すが、$M_{ch}$ と $b_r$ の間の偏差は大きい。特に、周辺衝突（例：30–40%）において、$b_r$ サンプルに含まれる期待範囲外の衝突パラメータを持つイベントの割合は、$M_{ch}$ サンプルと比較して最大 60% 異なる。
• 陽子収量: 自由陽子の収量は、EoS の硬さと中心度決定法の両方に敏感である。
  • 幾何学的解釈（$b_f$）を用いる場合、中心度選択によって導入される不確実性は、EoS によって誘起される効果よりも弱い。
  • グローバー MC ベースの手法（$b_r$）を用いる場合、中心度に関連する不確実性の影響は、特に周辺衝突において、EoS の効果よりも顕著になる。収量の $M_{ch}$ と $b_r$ 選択間の比は、中心度が増すにつれて大幅に増加する。
• 集団流: 直接流の傾き（$v_{11}$）および楕円流（$v_{20}$）は、主に EoS の硬さに支配される。$v_{11}$ は中心度決定法に対して概ね鈍感であるが、$v_{20}$ はわずかな感度を示し、$b_r$ で選択されたイベントは $M_{ch}$ で選択されたイベントよりもわずかに大きな値を示す。
• 横運動量: $p_T$ 分布および平均横運動量は、収量と同様の傾向に従う。$p_T$ 分布は、$b_r$ を用いた場合の中心度定義に強く依存するが、平均横運動量 $\langle p_T \rangle$ は EoS および中心度決定法の両方に対して比較的鈍感であることが判明した。

意義および結論
本論文は、高密度核 EoS に対する定量的な制約を課すためには、荷電粒子多重度（$M_{ch}$）と衝突パラメータとの間の厳密かつ一貫したマッピングが不可欠であると結論付けている。

• 手法の妥当性: SIS18 エネルギー（$\sqrt{s_{NN}} = 2.4$ GeV）において、中心度の幾何学的解釈（$b_f$）は妥当であり、動力学的な多重度選択と一致している。対照的に、グローバー MC ベースの中心度決定（$b_r$）は、このエネルギー領域においては信頼性が低くなる。これは、粒子多重度が参加核子数に厳密に比例するという仮定の崩壊によるものと考えられる。
• EoS 制約への影響: 本研究は、不適切な中心度決定法（特に低エネルギーにおけるグローバー MC）を使用すると、物理的な EoS 効果に匹敵するか、あるいはそれを上回る系統誤差が導入される可能性があることを示している。したがって、将来の核 EoS 制約においては、使用される特定の中心度選択法を考慮し、実験的な観測量と幾何学的な中心度との間のマッピングが、調査対象のエネルギー領域において物理的に正当化されていることを確認する必要がある。

今後の展望
著者らは、将来の研究において、変動（fluctuations）や異なる中心度クラス間の相関に焦点を当て、シミュレーションと実験データとの間の系統的な偏差を定量的に評価するために、ベイズ推論手法を用いるべきであると示唆している。また、低・中間エネルギーにおける EoS 研究の信頼性を向上させるために、グローバー MC ベースの手法の下にあるエネルギー依存性とモデルの仮定に関する、より詳細な調査が必要である。