Abstract:We investigate the charmonium and bottomonium production in $pp$ collisions using the Wigner densities formalism. The Wigner density of the quarkonia is approximated by analytical 3-D isotropic harmonic oscillator Wigner densities with the same root-mean-square radius given by the solution of the Schrödinger equation. This approach reproduces quite well the available experimental transverse momentum and rapidity distributions.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the production mechanisms of charmonium and bottomonium in high - energy proton - proton (pp) collisions. Specifically, the authors use the Wigner density formalism to study the production of these heavy - flavor mesons and attempt to reproduce the experimentally observed transverse momentum and rapidity distributions through this method.
### Main research contents:
1. **Theoretical background**: Heavy - flavor mesons (such as charmonium and bottomonium) are important tools for studying the strongly - interacting quark - gluon plasma (QGP). Since the mass of heavy quarks is much larger than the QCD scale parameter \(\Lambda_{\text{QCD}}\), the production process of heavy - flavor mesons can be decomposed into two parts:
- **Hard process**: Produce a heavy - quark pair \(Q\bar{Q}\), and this part can be described by perturbative QCD.
- **Soft process**: Form colorless heavy - flavor mesons from the heavy - quark pair, and this part is non - perturbative and requires other models to describe, such as the Color - Evaporation Model, the Color - Singlet Model, and the Color - Octet Model.
2. **Method**: The authors use the Wigner density formalism to describe the production of heavy - flavor mesons. The Wigner density is obtained by the projection of the quantum density matrix, where the probability \(P_i\) of producing a certain meson can be expressed as:
\[
P_i=\text{Tr}(\rho_i\rho^{(N)})
\]
where \(\rho_i\) is the density matrix of meson \(i\), and \(\rho^{(N)}\) is the density matrix of \(N\) heavy quarks and antiquarks produced in pp collisions. Through partial Fourier transformation, we can get:
\[
\frac{dP_i}{d^3R d^3P}=\sum_j\int\frac{d^3r d^3p}{(2\pi)^6}W_i(r, p)\prod_{j > 2}\int\frac{d^3r_j d^3p_j}{(2\pi)^3(N - 2)}W^{(N)}(r_1, p_1, r_2, p_2,\ldots, r_N, p_N)
\]
3. **Model assumptions**: The authors assume that the unknown multi - body Wigner density can be replaced by the average value of the classical phase - space distribution:
\[
W^{(N)}\approx\langle W^{(N)}_{\text{classical}}\rangle
\]
The classical momentum - space distribution of heavy quarks is provided by EPOS4, but EPOS4 only provides vertex coordinate information. For heavy - quark pairs produced at the same vertex, assuming that their relative distances in the center - of - mass system follow a Gaussian distribution, the Wigner density can be expressed as:
\[
W^{(2)}(r, p)\sim r^2\exp\left(-\frac{r^2}{2\sigma^2_{Q\bar{Q}}}\right)f_{\text{EPOS4}}^{Q\bar{Q}}(p)
\]
where \(\sigma_{Q\bar{Q}}\) is the effective width.
4. **Results**: Through numerical calculations, the authors obtained the Wigner densities in different states and compared them with experimental data. The results show that the Wigner density formalism can well describe the experimentally observed transverse momentum and rapidity distributions, especially when choosing appropriate width parameters \(\sigma_{c\bar{c}} = 0.4\,\text{fm}\) and \(\sigma_{b\bar{b}} = 0.2\,\text{fm}\).
### Conclusion:
This paper successfully describes the production mechanisms of heavy - flavor mesons (charmonium and bottomonium) in high - energy pp collisions through the Wigner density formalism and is in good agreement with experimental data.