Charm quark production in heavy-ion collisions as a signature of pre-equilibrium

本論文は、現在の理論的な不確実性にもかかわらず、重イオン衝突における全チャーム生成の精密な測定が、初期のハード散乱の計算の改善と組み合わせることで、非平衡段階の性質を推論するためのシグネチャーとして機能し得ることを提案している。

要約

重イオン衝突（2つの重い原子核を光速に近い速度で衝突させること）を、激しい爆発とともに始まり、やがて穏やかな群衆へと落ち着いていく、大規模で混沌としたパーティーとして想像してみてください。

設定：「パーティー前」の混沌
これらの原子核が衝突すると、瞬時に滑らかで熱い粒子のスープ（クォーク・グルーオン・プラズマ、またはQGPと呼ばれるもの）になるわけではありません。落ち着く前に、短く混沌とした「パーティー前」のフェーズが存在します。この間、圧力は偏っており（前方よりも横方向への押し出しが強い）、成分（グルーオンやクォーク）はまだ均一に混ざり合っていません。科学者たちはこれを**非平衡フェーズ（pre-equilibrium phase）**と呼んでいます。

通常、科学者たちは、チャームクォークと呼ばれる重い粒子は、衝突の最初の瞬間における非常に短い「激しい衝突（ハードクラッシュ）」の際にのみ生成されると考えています。これは、金槌が金敷を叩いたときに飛び散る火花のようなものです。一度この初期の衝突が終わると、チャームクォークの数は一定に保たれると考えられています。

新しいアイデア：「パーティー前」の火花
この論文は、シンプルな問いを投げかけています。この混沌とした「パーティー前」のフェーズもまた、これらの重いチャームクォークを生成しているのではないか？ ということです。

著者らは、この非平衡フェーズが極めて高密度かつ高エネルギーであるため（後の穏やかなフェーズよりもさらに激しい）、実際にチャームクォークの工場になり得るのではないかと示唆しています。彼らは、軽い粒子（ジレプトン）がこのフェーズ中に生成されることが知られていることと比較しています。もし軽い粒子がここで作られるのであれば、重い粒子も作れるのではないか、というわけです。

実験：シミュレーションの実行
これをテストするために、著者らは複雑なコンピュータ・シミュレーション（高精度の気象モデルのようなものですが、対象は亜原子粒子です）を使用しました。彼らは、混沌としたパーティー前のフェーズを2つの異なるアプローチでモデル化しました。

1. 「現実的な」モデル： 粒子がどのように跳ね返り、相互作用するかを詳細にシミュレートしたもの（QCDキネティック理論）。
2. 「簡略化された」モデル： 混沌が特定のパターンに従うと仮定した、より滑らかで計算しやすいバージョン（ロマツクェ・ストリックランド・モデル）。

彼らは、システムが冷却される前の、この短い混沌とした窓口の間に、どれだけのチャーム・反チャーム対が誕生するかを計算しました。

結果：驚くべき寄与
結果は興味深いものでした。

• 実際に起こっている： 非平衡フェーズはチャームクォークを生成しています。それは単なる微々たるものではなく、「無視できない」量です。
• タイミング： 軽い粒子がイベント全体を通じて生成される可能性があるのに対し、重いチャームクォークは主に、混沌がピークに達する非常に早い段階で生成されます。
• 規模： 衝突の具体的な条件によりますが、この「パーティー前」の生成は、最終的な残骸に見られる全チャームクォークの**10%から50%**を占める可能性があります。これはかなりの割合です！

問題点：霧がかかった測定
ここには落とし穴があります。数学的にはこの追加の生成が存在すると示されていますが、現在の私たちはこれを現実世界のデータで証明することができません。

なぜでしょうか？それは、現在の重イオン衝突における全チャームクォーク生成量の測定値には、巨大な「不確実性の霧」が存在するからです。それは、メインスピーカー（初期のハードクラッシュ）が叫んでいる部屋の中で、ささやき声（非平衡チャーム）を聞き取ろうとしているようなものです。しかも、私たちはそのメインスピーカーが本来どの程度の大きさであるべきかさえ、正確には分かっていません。メインスピーカーに関する理論的な計算には大きな誤差が含まれており、その「ささやき」が実際に存在しているのか、それとも単なるノイズの一部なのかを判別することを不可能にしています。

解決策：より優れたマイクロフォン
この論文は、この隠れた「パーティー前」のチャームを見つけるためには、より精密な測定が必要であると結論付けています。

• 重イオン衝突における全チャーム生成量を、陽子衝突におけるものと同じ精度で測定する必要があります。
• 「核環境」が生成率をどのように変化させるかを、より深く理解する必要があります。

結論
この論文は、重イオン衝突の混沌とした初期の瞬間が、重いチャームクォークの隠れた工場であることを提案しています。現在の測定の不確実性のためにまだ明確に見ることはできませんが、将来の実験（次世代のALICE 3やLHCbのアップグレードなど）が十分に精密になれば、全チャームクォークの数を探偵の道具として使い、まさに「パーティー前」の混沌がどのように振る舞い、巨大な衝突の後に宇宙がいかに熱平衡化していくのかを解明できるはずです。

技術概要

技術要約：重イオン衝突におけるチャームクォーク生成と前平衡状態のシグネチャー

問題提起
超相対論的重イオン衝突において、クォーク・グルオン・プラズマ（QGP）の形成には、縦方向の圧力異方性と化学的非平衡（グルオン優勢）を特徴とする前平衡相が先行する。電磁放射（特に2–5 GeV/$c^2$の質量範囲におけるジレプトン）は、この相のプローブとして確立されているが、前平衡段階における重いクォーク（チャーム）の生成については、これまでほとんど探索されてこなかった。従来、チャーム生成は初期のハード散乱に起因すると考えられており、全チャーム収量はその後、保存されると仮定されている。しかし、原子核・原子核衝突における全チャーム生成断面積は、クォークオニウム抑制を解釈する上での主要な不確念の源であり続けている。本論文では、前平衡相におけるチャーム生成が全収量に対して無視できない割合を占めるのか、また、初期段階のグルオン豊富さに敏感な、ジレプトンに対する補完的なプローブとなり得るかを調査する。

手法
著者らは、摂動論的QCDと有効運動論を組み合わせた理論的枠組みを用いて、前平衡のダイナミクスをモデル化している。

• 生成率の計算： チャーム・反チャーム（$c\bar{c}$）対生成率は、摂動論の一次まで計算される。この生成率には、不変質量 $Q$ および入射粒子の相対速度に依存する行列要素を用いた、クォーク・反クォークの対消滅とグルオン融合からの寄与が含まれる。
• 前平衡ダイナミクス： グルーオンとクォークの位相空間分布の進化は、QCDの有効運動論（EKT）を用いてモデル化される。系は、ビョルケン時間 $\tau$ と剪断粘性係数 $\eta/s$ を用いたユニバーサルなスケーリング変数 $\tilde{w} = \tau T_{\text{eff}} / (4\pi\eta/s)$ によって記述される共形ビョルケン流として記述される。
• スケーリング解析： 著者らは、微分収量を無次元比（$m_c/T_{\text{eq}}$, $Q/T_{\text{eq}}$）の関数として表すことで、生成率の普遍性をテストしている。異なる結合強度（$\lambda = 5, 10, 20$）におけるEKTシミュレーションの結果を、運動量異方性とクォーク抑制をパラメータ化した簡略化されたRomatschke-Strickland（RS）モデルと比較している。
• 積分と比較： 瞬時生成率を時間および横断面積にわたって積分することで、全前平衡収量を得る。これらの結果を、$\sqrt{s_{NN}} = 5.02$ TeVにおけるPb+Pb衝突における初期ハード散乱からのチャーム生成の推定値と比較する。ハード散乱のベースラインは、nPDF（EPPS21）を用いたNLO pQCD計算（MadGraph）を用い、ALICEの陽子・陽子データに正規化して導出される。データ駆動型の代替案も提示されており、これは陽子・原子核修飾因子（$R_{pA}$）を組み合わせることで、Pb+Pbにおける初期ハード散乱収量を推定するものである。

主な結果

• 無視できない寄与： 本研究では、前平衡チャーム生成が全チャーム収量のかなりの割合を占めることが示された。ラピディティや特定の理論的仮定（結合強度、$\eta/s$）に応じて、この寄与はおよそ10%から、初期ハード散乱による収量に匹例する値までの範囲をとる。
• 初期条件への感度： 後期の段階でユニバーサルなスケーリングを示すジレプトン生成とは異なり、チャーム生成は初期条件への感度を保持している。積分された収量のスケーリング関数は、異なるEKT結合強度およびRSモデルの間で、特に質量がチャーム質量付近またはそれ以上の領域において、顕著な不一致を示している。これは、重いフレーバーの生成が、系が初期状態の記憶を失う前の、より早い時刻まで及んでいることに起因する。
• モデル依存性： RSモデルは一般にEKTシミュレーションよりも大きな収量を予測する。これは、EKTシミュレーションが初期時刻 $\tau = 1/Q_s$ から開始するのに対し、RSのパラメータ化が実質的に $\tau=0$ から開始しているためである。
• ラピディティ分布： 前平衡収量は、ジレプトン生成の知見と一致して、より低い $\eta/s$（例：0.16 vs 0.32）において大きくなることが判明した。
• 不確実性の支配： 前平衡シグナルを孤立させる能力は、現在、初期ハード散リングのベースライン計算における大きな理論的不確実性によって制限されている。nPDFおよびスケール変動に起因する核修飾因子（$R_{AA}$）の不確実性は極めて大きく（境界間で最大5倍）、その結果、中間ラピディティにおける前平衡の寄与は、ハード散乱の推定値の中央値に対して10〜50%を占める可能性があり、理論的な不確実性の高い境界内では、全収量の最大100%に達する可能性もある。

意義と主張
本論文は、初期ハード散乱の理論計算の精度向上と組み合わせることで、重イオン衝突における全チャーム生成の精密な測定を通じて、前平衡相に関する情報を推論できると結論付けている。著者らは、チャーム生成が前平衡相におけるグルオンの豊富さにも敏感であることから、ジレプトン生成に対して補完的なものであると主張している。

しかし、著者らは現在の抽出能力については慎重な姿勢を見せている。彼らは、原子核・原子核衝突における全チャーム生成断面積の大きな不確実性のために、前平衡チャーム生成の直接的な抽出はまだ不可能であると述べている。本論文は、今後の進展が以下の2点にかかっていると位置づけている：

1. 初期ハード散乱計算（特にnPDFとスケール変動）における理論的不確実性の低減。
2. 陽子・陽子および陽子・原子核系と同等の精度を持つ、重イオン衝突における全チャーム収量の精密測定の達成。

著者らは、ALICE 3、ALICE 2、LHCb Upgrade 1、およびLHCb Upgrade 2 (U1/U2) プログラムを、前平衡シグナルを孤立させるために必要な精度を達成するための重要なプログラムとして特定している。もしこれらの実験的目標が達成されれば、測定された原子核・原子核の収量と、陽子・陽子および陽子・原子核のデータから推定された収量との差が、重イオン衝突における熱化プロセスへの新たな直接的なアクセスポイントを提供することになる。