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

Basic Information

• Paper ID: 2310.12627
• Title: Energy dependence of $J/\psi$ production in pp collisions with the PACIAE model
• Authors: 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
• Categories: hep-ph (High Energy Physics - Phenomenology), nucl-th (Nuclear Theory)
• Publication Date: February 9, 2024
• Paper Link: https://arxiv.org/abs/2310.12627v2

Abstract

This study investigates $J/\psi$ meson production in proton-proton collisions at center-of-mass energies $\sqrt{s}$ of 2.76, 5.02, 7, 8, and 13 TeV using the Parton and Hadron Cascade model PACIAE 2.2a. The model is based on PYTHIA but extends the framework by incorporating parton scattering before hadronization and hadronic scattering after hadronization. Quantum chromodynamic processes for $J/\psi$ production are specifically selected within the PYTHIA module, and corresponding bias factors are introduced. The calculated total cross sections, transverse momentum and rapidity differential cross sections for $J/\psi$ in the forward rapidity region (2.5 < y < 4) show good agreement with corresponding experimental measurements. In the mid-rapidity region, the double differential cross sections at $\sqrt{s}=5.02, 7,$ and 13 TeV also agree well with experimental data. Additionally, the study provides interpolated double differential cross sections and total cross sections for the mid-rapidity region at $\sqrt{s}=8$ TeV for future experimental verification.

Research Background and Motivation

1. Core Problem: The $J/\psi$ meson, as the lightest vector charmonium state, is crucial for understanding the production mechanism and studying the quark-gluon plasma (QGP). In relativistic heavy-ion collisions, $J/\psi$ suppression is considered an important probe for QGP formation.
2. Problem Significance: To distinguish between hot medium effects and cold nuclear matter effects, it is essential to first understand the $J/\psi$ production mechanism in proton-proton collisions, which serves as an important baseline for studying nucleus-nucleus collisions.
3. Limitations of Existing Methods:
   • Theoretical approaches such as color singlet model, non-relativistic QCD, and color evaporation model cannot simultaneously describe $J/\psi$ polarization, transverse momentum spectra, and energy dependence of cross sections
   • PYTHIA 6.4 is insufficient in describing the correlation between $J/\psi$ yield and charged particle multiplicity
   • Lack of comprehensive models considering both parton and hadronic scattering effects
4. Research Motivation: Systematically investigate the energy dependence of $J/\psi$ production in pp collisions at different energies using the PACIAE 2.2a model, which incorporates parton and hadronic scattering processes beyond the PYTHIA framework.

Core Contributions

1. Model Application: First systematic study of $J/\psi$ production in pp collisions at multiple LHC energies using the PACIAE 2.2a model
2. Bias Factor Methodology: Introduction of a bias factor method for handling specially selected $J/\psi$ production processes
3. Comprehensive Comparison: Systematic comparison with ALICE experimental data in both forward rapidity region (2.5 < y < 4) and mid-rapidity region (-0.9 < y < 0.9)
4. Theoretical Predictions: Provision of theoretical predictions for $J/\psi$ production cross sections in the mid-rapidity region at $\sqrt{s}=8$ TeV

Methodology Details

Task Definition

Investigation of the total cross section, transverse momentum differential cross section, and rapidity differential cross section for $J/\psi$ meson production in proton-proton collisions as functions of collision energy.

Model Architecture

PACIAE 2.2a Model Structure:

1. Parton Initialization Stage: Generation of hard scattering, initial-state and final-state radiation based on PYTHIA 6.4
2. Parton Scattering Stage: Application of leading-order perturbative QCD parton-parton interaction cross sections with K-factor to account for higher-order effects
3. Hadronization Stage: Conversion of partons to hadrons through string fragmentation or coalescence models
4. Hadronic Scattering Stage: Consideration of elastic and inelastic scattering processes between hadrons

$J/\psi$ Production Process Selection:

• Settings: MSEL=0 and MSUB(i)=1 (i=86,87,88,89,104,105,106)
• Corresponding color singlet processes: $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$

Technical Innovations

1. Bias Factor Method: $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}$
2. K-factor Optimization: Determination of optimal K-factor through $\chi^2$ minimization: $\chi^2 = \sum_{i=1}^{N} \frac{(y_i - f(p_T^i; K))^2}{(\Delta_i)^2}$
3. Scattering Processes: Detailed consideration of $J/\psi$-hadron inelastic scattering processes, such as $J/\psi + \pi \to D + D^*$, etc.

Experimental Setup

Dataset

• Experimental Data: $J/\psi$ production data in pp collisions published by the ALICE Collaboration
• Energy Range: $\sqrt{s} = 2.76, 5.02, 7, 8, 13$ TeV
• Rapidity Intervals: Forward rapidity (2.5 < y < 4) and mid-rapidity (-0.9 < y < 0.9)
• Simulation Samples: $10^6$ events generated for each energy point

Evaluation Metrics

• Total cross section $\sigma_{J/\psi}$
• Transverse momentum differential cross section $d^2\sigma_{J/\psi}/dp_T dy$
• Rapidity differential cross section $d\sigma_{J/\psi}/dy$
• Goodness-of-fit $\chi^2/\text{ndf}$

Comparison Method

Direct comparison with ALICE experimental measurements, with reference to results from PYTHIA 6.4, PYTHIA 8.2, EPOS3, and other models

Implementation Details

• K-factor range: 1.1-1.9 with step size 0.1
• Optimal K-factor determined through $\chi^2$ minimization
• Inclusion of radiative decay contributions from $\chi_c$ states to $J/\psi$

Experimental Results

Main Results

1. Forward Rapidity Region Performance:
   • Transverse momentum differential cross sections show good agreement with experimental data (Figure 2)
   • Rapidity differential cross sections accurately reproduce experimental measurements (Figure 3)
   • Total cross section energy dependence trend is consistent with experiments (Figure 5)
2. Mid-Rapidity Region Results:
   • Double differential cross sections at $\sqrt{s} = 5.02, 7, 13$ TeV agree well with experimental data (Figure 7)
   • Total cross sections agree with experiments within error margins (Figure 8)
3. K-factor Optimization Results:
   • 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$)

Bias Factor Analysis

• Bias factors show slight decrease with increasing energy
• Bias factors for rapidity distributions are approximately 1.5 times higher than those for transverse momentum distributions
• Energy dependence can be described by $B = a + b\ln(\sqrt{s}/\text{TeV})$

Prediction Results

Theoretical predictions are provided for $J/\psi$ double differential cross sections and total cross sections in the mid-rapidity region at $\sqrt{s} = 8$ TeV, with total cross section approximately 15 μb.

Related Work

Theoretical Model Development

1. Color Singlet Model: The earliest theoretical framework for $J/\psi$ production
2. Non-relativistic QCD: Incorporation of color octet contributions
3. Color Evaporation Model: Assumption that all $c\bar{c}$ pairs have fixed probability of forming $J/\psi$

Monte Carlo Simulation Studies

• Application of PYTHIA series in multiplicity correlation studies
• Use of EPOS3 and UrQMD in $J/\psi$ suppression research
• Present work represents the first systematic application of PACIAE model to $J/\psi$ production studies

Conclusions and Discussion

Main Conclusions

1. The PACIAE 2.2a model successfully describes $J/\psi$ production in pp collisions at multiple LHC energies
2. The bias factor method effectively handles sample bias from special process selection
3. The model shows good agreement with experimental data in both forward and mid-rapidity regions
4. Provides verifiable theoretical predictions for $\sqrt{s} = 8$ TeV

Limitations

1. Polarization effects of $J/\psi$ are not considered
2. Only color singlet production processes are considered, neglecting color octet contributions in high $p_T$ regions
3. Contributions from $\psi(2S) \to J/\psi$ decay are not included
4. Physical interpretation of K-factor requires further investigation

Future Directions

1. Simultaneous description of $J/\psi$ polarization and transverse momentum spectra
2. Separate investigation of contributions from parton and hadronic scattering effects
3. Extension to higher transverse momentum regions and color octet processes
4. Application to $J/\psi$ suppression studies in heavy-ion collisions

In-Depth Evaluation

Strengths

1. Methodological Innovation: First application of PACIAE model to $J/\psi$ production research, incorporating complete scattering processes
2. Strong Systematicity: Coverage of major LHC operating energies with comprehensive energy dependence study
3. Technical Rigor: Reasonable bias factor methodology and standardized K-factor optimization procedure
4. Sufficient Experimental Comparison: Comprehensive comparison and validation with ALICE experimental data

Weaknesses

1. Theoretical Depth: Insufficient analysis of physical mechanisms underlying scattering effects
2. Applicable Range: Limited to low and intermediate $p_T$ regions; high $p_T$ behavior not addressed
3. Model Limitations: Polarization and other important physical quantities not considered
4. Prediction Verification: Predictions at $\sqrt{s} = 8$ TeV await experimental verification

Impact and Significance

1. Academic Value: Provides new theoretical tools for understanding $J/\psi$ production mechanisms in pp collisions
2. Practical Significance: Provides important baseline calculations for heavy-ion collision research
3. Method Generalization: Bias factor methodology can be extended to other charmonium studies

Applicable Scenarios

• $J/\psi$ production research in pp collisions within LHC energy range
• Baseline calculations for $J/\psi$ suppression in heavy-ion collisions
• Reference for investigating production mechanisms of other heavy quark bound states

References

This paper cites 50 important references covering theoretical models, experimental measurements, and Monte Carlo simulations of $J/\psi$ production, providing a comprehensive literature foundation for the field.