Quarkonium production in high energy pp collisions

. 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.

Hidden heavy flavour mesons are a useful tool to study the strongly interacting quark gluon plasma, which is created during high energy heavy-ion collisions.Due to the large quark mass m Q , m Q ≫ Λ QCD , with Λ QCD being the QCD cutoff, quarkonium production can be factorized into the production of a heavy quark, Q Q, pair, which can be described by perturbative QCD and a subsequent soft non-perturbative process, which describes the formation of a colorless quarkonium from the Q Q pair.For the latter part many approaches have been advanced.They include the Color-Evaporation Model , the Color-Singlet Model and the Color-Octet Model.The latter two are encompassed in the NRQCD approach.For a review we refer to [1].More recently, the rapidity and p T distributions of hidden heavy flavour mesons, produced in pp collisions, have also been well reproduced in the Wigner density matrix formalism [2][3][4].Here we will extend this formalism up to 3S and use EPOS4 to generate the initial heavy quarks.
The Wigner density formalism is based on the quantal density matrix projection in which the probability that a meson i is produced is given by P i = T r(ρ i ρ (N) ) with ρ i being the density matrix of the meson i and ρ (N) the density matrix of the N heavy quarks and antiquarks, produced in a pp collision.A partial Fourier transformation of the density matrices yields then W i is the two-body Wigner density of the bound heavy quark pair and W (N) (r 1 , p 1 , r 2 , p 2 , ..., r N , p N ) is the quantal density matrix in Wigner representation of the ensemble of N heavy quarks produced in a pp collision.r(R) and p(P) are the relative (center of mass) coordinate and momentum of the heavy quark and antiquark, which are  1.The masses, root-mean-square radius, and the Gaussian widthes σ of different charmonium and bottomonium states in vacuum.The experimental data are from Ref. [7].
bound in a quarkonium.We assume that the unknown quantal N-body Wigner density can be replaced by the average of classical phase space distributions, W (N) ≈ ⟨W (N) classical ⟩.The classical momentum space distributions of the heavy quarks is provided by EPOS4 [5,6].EPOS4, however, provides only the coordinate information of the vertex where the Q Q pair is created.For a heavy quark pair, created at the same vertex, we assume that the relative distance between Q and Q in their center-of-mass frame is given by a Gaussian distribution and the Wigner density for a Q Q pair can be expressed as where the distance is controlled by the effective width σ Q Q.Having now momenta and positions of the heavy quarks we can calculate the yield of charmonium and bottomonium via Eq.(1).We come now to the construction of the quarkonium Wigner density.The Wigner density is obtained by the Wigner-Weyl transformation of the density matrix of the quarkonia.The quarkonium wave function is the solution of the two-body Schrödinger equation, which we solve for charmonium and bottomonium with the Cornell potential, V(r) = −α/r + κr + c with α = 0.513, κ = 0.17GeV 2 , c = −0.161,and with the quark masses m c = 1.5GeV and m b = 5.2GeV.The wave functions are shown in Fig. 1, and the masses and root-meansquare radii ⟨r 2 ⟩ are shown in Table 1.We can see that the masses are very close to the experimental values.The wave function is, however, not analytical, and therefore the Wigner density can only be calculated numerically , what makes the solution of Eq. ( 1) complicated.It is therefore convenient to approximate the wave function by a 3-D isotropic harmonic oscillator wave function for the potential V(r) = 1/(2m Q σ 4 )r 2 .These wave functions are analytical and can be expressed as ψ nlm (r, θ, ϕ) = R nl (r)Y l,m (θ, ϕ), where Y l,m are the spherical harmonics.The radial part can be expressed as, where L l+1/2 n are Laguerre polynomials.The parameters of the 3-D isotropic harmonic oscillator wave functions are chosen to match the root-mean-square radius ⟨r 2 ⟩ of the real quarkonium wave function: ⟨r 2 ⟩ = 3σ 2 /2 for 1S , ⟨r 2 ⟩ = 5σ 2 /2 for 1P, ⟨r 2 ⟩ = 7σ 2 /2 for 1D and 2S , ⟨r 2 ⟩ = 9σ 2 /2 for 2P, and ⟨r 2 ⟩ = 11σ 2 /2 for 3S states.The corresponding widths are shown in Table 1 and the wave functions are shown in Fig. 1 with dashed lines.We can see that the ground states and low lying excited states can be well reproduced by the 3-D isotropic harmonic oscillator, while the difference increases for higher excited states, e.g.2S , 2P, and 3S .
and bottomonium (see Eq. ( 1)) are shown in Figs. 2 and 3. We see, as far as data are available, We can conclude that the experimentally available rapidity and transverse momentum distribution of cc and b b quarkonia can be well described in the Wigner density formalism.
only parameter which enters the calculation is the width of the distribution of the relative distance of the Q Q pair at production.The relative contribution of the different states is then exclusively given by their wave function.