$D^0$-$D_s^+$ Elliptic-Flow Splitting under Event-Shape Engineering: A Probe of Sequential Charm Hadronization

要約

巨大で超高温の、小さな粒子からなる「スープ」を想像してみてください。これは、2つの重い鉛原子が光速に近い速度で衝突したときに生成されます。この「スープ」は「クォーク・グルーオン・プラズマ（QGP）」と呼ばれます。このスープの中では、「チャームクォーク」と呼ばれる重い粒子が泳ぎ回っています。スープが冷却されるにつれて、これらのクォークはより軽い粒子と結びつき、新しい安定した粒子である「ハドロン」（具体的には、2種類のD中間子：D⁰とD⁺ₛ）を形成します。

長い間、科学者たちはこれらすべての新しい粒子が、まるで人々が一斉に建物から出てくる時のように、全く同じ瞬間に形成されると考えてきました。しかし、この論文は異なる物語、すなわち**「逐次的なハドロン化（sequential hadronization）」**を提示しています。

以下に、著者らが発見した内容を、日常的な比喩を用いて簡単に解説します。

1. 2つの物語：一斉の退場 vs. 時差のある退場

• 古い物語（同時）： コンサート会場から人々が退場する場面を想像してください。全員が全く同時に出口へ向かって歩いていきます。もし、2つの異なるグループ（例えば、赤い帽子を被った人と青い帽子を被った人）を見たとしても、両グループとも群衆によって同じように押し流されます。
• 新しい物語（逐次的）： コンサートが終わり、出口が混雑している場面を想像してください。
  • D⁺ₛ粒子は、「VIPパス」を持っている人々のようなものです（彼らは強く結合しています）。彼らは、部屋がまだ非常に熱く混沌としている時期（臨界温度の約1.2倍）に、群衆をすり抜けて早く脱出します。
  • D⁰粒子は、一般の参加者のようです。彼らは最後までスープの中に留まり、スープが最後（臨界温度 $T_c$）に達するまで泳ぎ続けます。
  • 結果： D⁰はスープの中に長く留まったため、群衆の渦巻く流れによってより強く押し流されました。彼らは、早く脱出したD⁺ₛよりも多くの「スピン（流れ）」を拾い上げました。

2. 問題点：どのようにして違いを見分けるのか？

科学者は、これらの粒子がどれくらい「回転（スピン）」しているか（これは楕円流と呼ばれます）を測定することができます。しかし、落とし穴があります。スピンの量は2つの要素に依存します。

1. 衝突がどのように始まったか： 衝突は完璧な正面衝突だったのか、それともかすめるような衝突だったのか？（これは「イベントの形状」です）。
2. いつ脱出したか： 彼らは早く脱出したのか、それとも遅く脱出したのか？

もし、すべての衝突を混ぜて見てしまうと、ある粒子がスピンを多く持っている理由が、「遅く脱出したから」なのか、それとも単に「衝突の形状がスピンを生み出しやすい形だったから」なのかを判別することが困難になります。それは、あるランナーが足が速い理由が、「生まれつきの優れたアスリートだから」なのか、それとも単に「追い風があったから」なのかを推測するようなものです。

3. 解決策：「イベント形状エンジニアリング（ESE）」（風洞実験）

著者らは、**「イベント形状エンジニアリング（ESE）」**という巧妙なトリックを用いました。これは、風洞実験のようなものです。

• 彼らは数千回の衝突を取り上げ、2つのグループに分類しました。
  • 大きな $q^2$（強風）： 非常に強く歪んだ形状から始まった衝突。
  • 小さな $q^2$（弱風）： より丸く穏やかな形状から始まった衝突。
• これら2つのグループを比較することで、粒子が衝突の幾何学的形状による「風」に対してどのように反応するかを見ることができました。

4. 発見：「傾き」が真実を語る

データを確認したとき、著者らは「逐次的（Sequential）」な物語が正しいことを証明する決定的な証拠を見つけました。

• 「傾き」 ($\chi$)： 粒子が得るスピンの量が、風が強くなるにつれてどのように変化するかをグラフに描く場面を想像してください。
  • 逐次的な物語（D⁰が長く留まる場合）では、D⁰粒子は風に対して非常に敏感です。風が強くなると、彼らのスピンは劇的に増加します。一方で、早く脱出したD⁺ₛは、それほど反応しません。
  • ルール： D⁰の「感度の傾き」は、D⁺ₛよりも急になります。
  • 同時の物語（彼らが一緒に脱出する場合）では、両方の粒子は同じように反応します。つまり、両者の傾きは同一になります。

論文によれば、セミセントラル衝突（スープが十分に長く存在しつつ、かつ歪みも残っている「スイートスポット」）において、実際にD⁰の傾きはD⁺ₛよりもはるかに急であることが示されています。これは、D⁰がより多くの流れを捉えるために、より長くスープの中に留まっていることを証明しています。

5. なぜ単なる数字の問題ではないのか

著者らは、これが単なる数字のトリック（例えば、特定の衝突においてD⁰がD⁺ₛよりも多いといったこと）ではないことも確認しました。彼らは、D⁺ₛとD⁰の比率を調べました。

• 発見： その比率は、風が強いか弱いかにかかわらず一定でした。
• 意味： これにより、スピンの違いは、どちらかの粒子が多いことによるものではなく、純粋に彼らがいつスープを脱出したかという動力学的な効果によるものであることが確認されました。

まとめ

この論文は、重い粒子はすべて熱いスープから同時に脱出するわけではないという説を提唱しています。「VIP」粒子（D⁺ₛ）は早く去り、一方で「一般」の粒子（D⁰）はより長く留まり、より強く押し流されます。

衝突の形状によって衝突を分類する手法（イベント形状エンジニアリング）を用いることで、著者らはユニークな指紋を見つけました。それは、「一般」の粒子は「VIP」粒子よりも衝突の形状に対してはるかに強く反応するというものです。この反応の差こそが、彼らが異なるタイミングでスープを脱出したことの証明であり、初期宇宙における物質の形成という隠されたタイムラインを明らかにしているのです。

技術概要

技術要約：イベントシェイプ・エンジニアリング下における $D^0$–$D_s^+$ 楕円流分裂

問題提起
クォーク・グルーオン・プラズマ（QGP）におけるオープンチャーム・ハドロンのハドロン化メカニズムは、相対論的重イオン物理学における極めて重要な未解決問題である。従来のモデルは通常、すべてのチャーム・ハドロン種が臨界温度（$T_c$）において共通の超面上に形成されるという同時ハドロン化を仮定している。しかし、近年の理論的提案では、より強く結合した種（例：$D_s^+$）はより高い温度（例：$\sim 1.2T_c$）で早期に形成され、結合が弱い種（例：$D^0$）は後から$T_c$で形成されるという逐次ハドロン化が示唆されている。

これまで、楕円流の差 $\Delta v_2 \equiv v_2(D^0) - v_2(D_s^+)$ はこれらのシナリオを区別するためのプローブとして提案されてきたが、包括的な測定には根本的な限界がある。重いフレーバー・ハドロンの楕円流は、初期衝突幾何学と部分的な進化の時間の両方に依存する。この縮退により、ハドロン化時間の順序という特定の信号を、グローバルな幾何学的効果から分離することが困難となっている。

手法
著者らは、以下の要素を組み合わせたイベントごとのシミュレーションであるSHELLフレームワークを用いている：

1. (3+1)D 粘性流体力学： CLViscコードを用いてQGP媒体を進化させ、局所的な温度 $T(x)$ と流速 $u^\mu(x)$ を提供する。
2. ランジュバン輸送： 定数空間拡散係数 $2\pi T D_s = 2.5$ を用いてチャームクォークの拡散をモデリングし、ドラッグ、拡散、および媒体誘起グルーオン放射（高次ツイスト・アプローチ）を含む。
3. ハイブリッド・ハドロン化： コアレセンス（合体）とフラグメンテーション（断片化）を組み合わせたモデル。
   • 逐次シナリオ： $D_s^+$ の形成は、$T \approx 1.2T_c$ におけるより早い超面上で評価される。もしチャームクォークが $D_s^+$ へとコアレセンスできなかった場合、そのクォークは $D^0$（または他の種）を形成するために $T_c$ まで伝播を続け、チャーム数保存則を強制する。
   • 同時ベースライン： すべての種が同一の $T_c$ 超面上に形成される。
4. イベントシェイプ・エンジニアリング (ESE)： 固定された中心度クラス（0–10% および 30–50%）内で、低次の第2次フローベクトル $q_2$ に基づいてイベントを選択する。イベントは「大-$q_2$」（強い初期偏心率）と「小-$q_2$」（弱い初期偏心率）のクラスに分類される。これにより、全体の中心度に依存することなく、ハドロン化時間の応答を微分的に測定することが可能となる。

主な貢献と結果

• 幾何学と時間のデカップリング： 本研究は、ESEが包括的な測定よりも逐次ハドロン化と同時ハドロン化を鋭敏に識別できることを示している。固定された中心度内で $q_2$ を変化させることにより、バルクの幾何学的応答から種固有のハドロン化時間の応答を分離できる。
• フロー分裂 ($\Delta v_2$)：
  • 逐次シナリオでは、中間横運動量領域（$p_T \approx 3\text{--}4$ GeV/c）において、正の分裂 $\Delta v_2(D^0 - D_s^+) > 0$ が現れる。これは、$D_s^+$ が早期に凍結する一方で、$D^0$ の親クォークが $1.2T_c \to T_c$ の間に追加の楕円流を蓄積するためである。
  • 同時ベースラインでは、分裂はゼロに近いか負となり、正の増強は見られない。
  • 決定的なことに、逐次シナリオにおける正の分裂は、$q_2$ とともに系統的に増大する。
• 応答スロープ ($\chi$)： 著者らは、$v_2$ を $q_2$ の線形関数（$v_2 = \chi q_2 + v_2^{(0)}$）としてフィッティングすることで、応答スロープ $\chi$ を定義している。
  • 逐次予測： 明確な階層関係 $\chi(D^0) > \chi(D_s^+)$ が観察される。$D^0$ は後期に形成されるため、後期の幾何学的進化に対してより敏感である。
  • 同時予測： $\chi(D^0) - \chi(D_s^+)$ はほぼゼロまたは負となる。
  • このスロープは、全体のフローの正規化に影響されない、ロバストでモデルに依存しない識別子として機能する。
• 最適な中心度窓： 30–50% の半中心度クラスが、観測のための最適な窓として特定された。中央衝突（0–10%）はQGPの寿命は長いが、初期偏心率が小さいため、後期の相互作用を楕円流へと変換する能力が限定される。周辺衝突は偏心率は大きいが寿命が短い。30–50% クラスは、後期のフロー変換を最大化するために必要な非単調な相互作用を提供している。
• 化学的収量の安定性： $D_s^+/D^0$ の収量比およびその $q_2$ ダブル比（$R_{q_2}$）は、すべての選択において 1 に近い値を保っている。これは、観察されたフロー分裂が、イベントシェイプによる化学的収量やコアレセンス分率の変化によるものではなく、動力学的なフロー効果であることを裏付けている。

意義と主張
本論文は、イベントシェイプ・エンジニアリング下での $\Delta v_2(D^0 - D_s^+)$ および 応答スロープ $\chi$ が、QCD転移温度付近におけるチャーム・ハドロン化の時空構造に対する相補的な微分プローブを構成すると主張している。

主な意義は、以下の点を確立したことにある：

1. 逐次ハドロン化は、従来の同時ベースラインと比較して、$D^0$–$D_s^+$ のフロー分裂の符号を自然に反転させること。
2. ESEは、幾何学駆動の増強からハドロン化時間の応答を分離する。これは、幾何学（$q_2$ による）と時間（形成温度による）という「2つのノブ」を持つフレームワークを提供し、これはオープンチャームに限らず他の種にも一般化可能であること。
3. 階層関係 $\chi(D^0) > \chi(D_s^+)$ は、検証可能な定量的予測を提供する。もし将来のALICEデータの解析において $\chi(D^0) \approx \chi(D_s^+)$ が見出された場合、逐次的な描像は否定されることになる。
4. 結果は、半中心度領域がこれらメカニズムを区別するための最もクリーンな環境であることを検証している。これは、QGPの寿命と初期偏心率の相互作用が、観測可能な信号を最大化するためである。

著者らは、これらの知見が彼らの先行研究 [19] で導入された特定の逐次ハドロン化の理論的枠組みに基づいていることを強調しており、この効果を実験的に確認したと主張しているのではなく、そのようなテストを行うために必要な理論的ツールと具体的な観測量を提供している。