2025-11-23T21:07:17.387352

Energy dependence of $J/Ï$ production in pp collisions with the PACIAE model

Ye, Wang, Shi et al.

In this work we investigate the $J/Ï$ production in proton-proton collisions at the center-of-mass energy ($\sqrt{s}$) equal to 2.76, 5.02, 7, 8 and 13 TeV with a parton and hadron cascade model PACIAE 2.2a. It is based on PYTHIA but extended considering the partonic and hadronic rescatterings before and after hadronization, respectively. In the PYTHIA sector the $J/Ï$ production quantum chromodynamics processes are selected specially and a bias factor is proposed correspondingly. The calculated total cross sections, the differential cross sections as a function of the transverse momentum and the rapidity of $J/Ï$ in the forward rapidity region reproduce the corresponding experimental measurements reasonably well. In the mid-rapidity region, the double differential cross sections at $\sqrt{s}=$ 5.02, 7 and 13 TeV are also in a good agreement with the experimental data. Moreover, we interpolate the double differential cross section as well as the total cross section of $J/Ï$ in the mid-rapidity region at $\sqrt{s}=$ 8 TeV, which could be validated if the experimental data is available.

academic

pp衝突におけるJ/ψ生成のエネルギー依存性：PACIAEモデルによる研究

基本情報

論文ID: 2310.12627
タイトル: Energy dependence of $J/\psi$ production in pp collisions with the PACIAE model
著者: Kai-Fan Ye, Qiang Wang, Jia-Hao Shi, Zhi-Ying Qin, Wen-Chao Zhang, An-Ke Lei, Zhi-Lei She, Yu-Liang Yan, Ben-Hao Sa
分類: hep-ph（高エネルギー物理-現象論）、nucl-th（核理論）
発表日: 2024年2月9日
論文リンク: https://arxiv.org/abs/2310.12627v2

要約

本研究では、部分子および強子カスケードモデルPACIAE 2.2aを用いて、重心エネルギー $\sqrt{s}$ が2.76、5.02、7、8、13 TeVの陽子-陽子衝突におけるJ/ψ中間子の生成を調査した。本モデルはPYTHIAに基づいているが、強子化前の部分子散乱および強子化後の強子散乱を考慮するように拡張されている。PYTHIA部分ではJ/ψ生成の量子色力学過程を特別に選択し、それに応じてバイアス因子を導入した。計算された全断面積、前方ラピディティ領域 $J/\psi$ の横運動量および速度微分断面積は、対応する実験測定結果とよく一致している。中央ラピディティ領域では、 $\sqrt{s}=5.02、7、13$ TeVの二重微分断面積も実験データとよく一致している。さらに、本研究は $\sqrt{s}=8$ TeV中央ラピディティ領域のJ/ψ二重微分断面積および全断面積を内挿し、実験データが利用可能になった時点で検証することができる。

研究背景と動機

中心課題: J/ψ中間子は最も軽いベクトルチャーム偶素中間子として、その生成機構の理解はクォーク-グルーオンプラズマ(QGP)の研究に不可欠である。相対論的重イオン衝突において、J/ψの抑制はQGP形成の重要なプローブと考えられている。
問題の重要性: 熱媒質効果と冷核物質効果を区別するためには、まずpp衝突におけるJ/ψの生成機構を理解する必要があり、これは核-核衝突研究の重要なベンチマークである。
既存手法の限界:
- 色一重項モデル、非相対論的QCD、色蒸発モデルなどの理論的手法は、J/ψの偏極、横運動量スペクトル、および断面積のエネルギー依存性を同時に記述することができない
- PYTHIA 6.4はJ/ψ収率と荷電粒子多重度の相関を記述する際に不十分である
- 部分子および強子散乱効果を考慮した完全なモデルが欠けている
研究動機: PACIAE 2.2aモデルを用いて、異なるエネルギーでのpp衝突におけるJ/ψ生成のエネルギー依存性を系統的に研究する。本モデルはPYTHIAの基礎上に部分子および強子散乱過程を追加している。

核心的貢献

モデル応用: PACIAE 2.2aモデルを複数のLHCエネルギーでのpp衝突におけるJ/ψ生成研究に初めて系統的に適用
バイアス因子法: 特別に選択されたJ/ψ生成過程を処理するためのバイアス因子法を提案
包括的比較: 前方ラピディティ領域(2.5 < y < 4)および中央ラピディティ領域(-0.9 < y < 0.9)でALICE実験データとの系統的比較を実施
予測結果: $\sqrt{s}=8$ TeV中央ラピディティ領域のJ/ψ生成断面積の理論予測を提供

部分子初期化段階: PYTHIA 6.4に基づいて硬散乱、初期状態および最終状態放射を生成
部分子散乱段階: 主要次微扰QCD部分子-部分子相互作用断面積を使用し、K因子を導入して高次効果を考慮
強子化段階: 弦断裂または凝集モデルを通じて部分子を強子に変換
強子散乱段階: 強子間の弾性および非弾性散乱過程を考慮

J/ψ生成過程の選択:

MSEL=0およびMSUB(i)=1 (i=86,87,88,89,104,105,106)を設定
対応する色一重項過程: $gg \to J/\psi g$ 、 $gg \to \chi_{0c}g$ 、 $gg \to \chi_{1c}g$ 、 $gg \to \chi_{2c}g$ 、 $gg \to \chi_{1c}$ 、 $gg \to \chi_{2c}$ 、 $gg \to J/\psi\gamma$

技術的革新点

バイアス因子法: $B = \frac{\int_{p_T} \frac{d^2\sigma_{J/\psi}}{dp_T dy}\Big|_{\text{exp}} dp_T}{\int_{p_T} \frac{d^2\sigma_{J/\psi}}{dp_T dy}\Big|_{\text{sim}} dp_T}$
K因子最適化: $\chi^2$ 最小化法により最適K因子を決定: $\chi^2 = \sum_{i=1}^{N} \frac{(y_i - f(p_T^i; K))^2}{(\Delta_i)^2}$
散乱過程: $J/\psi + \pi \to D + D^*$ などの詳細なJ/ψ-強子非弾性散乱過程を考慮

実験設定

データセット

実験データ: ALICE協力団体が発表したpp衝突J/ψ生成データ
エネルギー範囲: $\sqrt{s} = 2.76、5.02、7、8、13$ TeV
ラピディティ区間: 前方ラピディティ(2.5 < y < 4)および中央ラピディティ(-0.9 < y < 0.9)
シミュレーション試料: 各エネルギーポイントで $10^6$ イベントを生成

評価指標

全断面積 $\sigma_{J/\psi}$
横運動量微分断面積 $d^2\sigma_{J/\psi}/dp_T dy$
ラピディティ微分断面積 $d\sigma_{J/\psi}/dy$
$\chi^2/\text{ndf}$ 適合度

比較方法

ALICE実験測定結果との直接比較、PYTHIA 6.4、PYTHIA 8.2、EPOS3などのモデル研究結果を参照

実装の詳細

K因子範囲: 1.1-1.9、ステップ0.1
最適K因子は $\chi^2$ 最小化により決定
$\chi_c$ 状態からJ/ψへの放射減衰寄与を考慮

実験結果

主要結果

前方ラピディティ領域の性能:
- 横運動量微分断面積は実験データとよく一致(図2)
- ラピディティ微分断面積は実験測定を正確に再現(図3)
- 全断面積のエネルギー増加傾向は実験と一致(図5)
中央ラピディティ領域の結果:
- $\sqrt{s} = 5.02、7、13$ TeVの二重微分断面積は実験データと一致(図7)
- 全断面積は誤差範囲内で実験と符合(図8)
K因子最適化結果:
- 2.76 TeV: K = 1.6 ( $\chi^2/\text{ndf} = 3.46/6$ )
- 5.02 TeV: K = 1.7 ( $\chi^2/\text{ndf} = 10.92/6$ )
- 7 TeV: K = 1.3 ( $\chi^2/\text{ndf} = 16.62/6$ )
- 8 TeV: K = 1.5 ( $\chi^2/\text{ndf} = 13.40/6$ )
- 13 TeV: K = 1.7 ( $\chi^2/\text{ndf} = 54.86/7$ )