Proof of Factorization of $\chi_{cJ}$ Production in Non-Equilibrium QCD at RHIC and LHC in Color Singlet Mechanism

Recently we have proved the factorization of NRQCD S-wave heavy quarkonium production at all orders in coupling constant. In this paper we extend this to prove the factorization of infrared divergences in $\chi_{cJ}$ production from color singlet $c{\bar c}$ pair in non-equilibrium QCD at RHIC and LHC at all orders in coupling constant. This can be relevant to study the quark-gluon plasma at RHIC and LHC.


I. INTRODUCTION
The factorization of infrared divergences in nonrelativistic QCD (NRQCD) color octet Swave heavy quarkonium production at high energy colliders at all orders in coupling constant is recently proved in [1]. In this paper we extend this formalism to non-equilibrium QCD by using the closed-time path integral formulation to prove the factorization of infrared divergences in χ cJ production from the color singlet cc pair in non-equilibrium QCD at all orders in coupling constant at RHIC and LHC. We also predict the correct definition of the non-perturbative matrix element of the χ cJ production from color singlet cc pair in nonequilibrium QCD at RHIC and LHC. This can be relevant to study the quark-gluon plasma (QGP) at RHIC and LHC.
At very high temperature (≥ 200 MeV) the normal hadronic matter becomes a new state of matter known as the QGP. About 10 −12 seconds after the big bang our universe was filled with the QGP which makes it important to produce it in the laboratory at RHIC and LHC by colliding two heavy nuclei at very high energy [2]. Since the confinement in QCD prevents us to detect the QGP directly at RHIC and LHC, various indirect signatures (such as the heavy quarkonium production/suppression [3]) are proposed for its detection.
Since the center of mass energy √ s = 200 GeV (5.5 TeV) of Au-Au (Pb-Pb) collisions at RHIC (LHC) is very high, the two nuclei at RHIC (LHC) travel almost at the speed of light creating the non-equilibrium quark-gluon plasma just after the heavy-ion collisions.
Because of the very small hadronization time scale in QCD (∼ 10 −24 seconds) there may not be enough secondary partonic collisions to bring this non-equilibrium QGP to equilibrium.
Hence the QGP at RHIC (LHC) may be in non-equilibrium where one can not define a temperature.
The hard (high p T ) parton production at RHIC and LHC can be calculated by using pQCD but the soft parton production calculation needs non-perturbative QCD which is not solved yet. This implies that there remains uncertainty in determining the soft partons production at RHIC and LHC. Note that the soft partons play an important role in determining the bulk properties of the QGP at RHIC and LHC.
It should be mentioned here that the study of hadronization from non-equilibrium QGP at RHIC and LHC is one of the most difficult and important problem because the confinement problem in QCD is not solved yet due to the lack of our understanding of non-perturbative QCD. This implies that the first principle calculation of hadron production from non-equilibrium partons at RHIC and LHC is not known.
Because of these reasons one finds that in order to detect the QGP at RHIC and LHC by using the first principle calculation one needs to study the nonequilibrium-nonperturbative QCD by using the closed-time path integral formalism which is not easy [4][5][6][7]. If one does not perform the exact first principle nonequilibrium-nonperturbative QCD calculation then the comparison of the theoretical calculation with the experimental data at RHIC and LHC becomes questionable. For example, some of the limitations of the present theoretical approaches are listed below.
The lattice QCD at finite temperature [8] is a common tool to study the properties of the QGP. However, for the reasons explained above, the actual QGP at RHIC and LHC may be in non-equilibrium where one can not define a temperature. Hence the lattice QCD at finite temperature has no application in non-equilibrium QGP at RHIC and LHC.
Similarly the hydrodynamics [9] is not applicable in non-equilibrium QGP at RHIC and LHC. Another limitation of the hydrodynamics [9] is that it does not answer the question how the partons become hadrons from first principle. As shown in [10] the parton to hadron fragmentation function in QCD in vacuum can not be used to study the hadrons production from partons from the quark-gluon plasma at RHIC and LHC. It is important to observe that even if the experimental data at RHIC and LHC is explained by using the hydrodynamics [9] it does not prove that the QGP is in equilibrium. In order to make sure that the QGP is in equilibrium at RHIC and LHC one has to prove that the same experimental data can not be explained by using the non-equilibrium QGP for which one has to study the nonequilibrium-nonperturbative QCD by using the closed-time path integral formalism.
As far as the actual physics at RHIC and LHC heavy-ion collisions is concerned the AdS/CFT based studies [11] and the supersymmetric Yang-Mills plasma based studies [12] have nothing to do it because of the lack of experimental verification of the string theory and the supersymmetry.
Regarding the initial condition for the QGP formation and the color glass condensate (CGC) [13], as discussed above, the hard (high p T ) parton production at RHIC and LHC can be calculated by using the pQCD but the soft parton production can only be correctly calculated from the first principle by using the non-perturbative QCD which is yet to be solved.
The jet quenching study, see for example [14,15], directly/indirectly uses the parton to hadron fragmentation function in QCD in vacuum. This is not possible because unlike the leading order perturbative gluon propagator in non-equilibrium QCD the non-perturbative fragmentation function in non-equilibrium QCD can not be decomposed into the vacuum part and the medium part [10].
Hence from the above discussions one finds that, although a lot of experimental data is available at RHIC and LHC heavy-ion colliders, but there exists no exact first principle theoretical calculation to explain these experimental data. It is almost impossible to make an exact first principle theoretical calculation at RHIC and LHC without studying the nonequilibrium-noperturbative QCD by using the closed-time path integral formalism.
The first principle way to study non-equilibrium quantum field theory is the Schwinger-Keldysh closed-time path (CTP) formalism [4,5]. Although the non-equilibrium QED is usually studied by using the canonical quantization formalism, the closed-time path integral formalism is useful to study the nonequilibrium-nonperturbative QCD due to the self gluon interactions and the hadronization.
As mentioned earlier, the heavy quarkonium is one of the indirect signature for the detection of QGP [3]. Both j/ψ and χ cJ are measured by various collaborations at the RHIC and LHC heavy-ion collider experiments. In order to study heavy quarkonium production from the QGP at RHIC and LHC one needs to prove factorization of infrared divergences, otherwise one will predict infinite cross section for the heavy quarkonium production.
The infrared divergences issue in the case of P-wave heavy quarkonium production is more complicated than that of the j/ψ production. This is because there are no uncanceled infrared divergences due to eikonal gluons exchange in the case of S-wave heavy quarkonium (j/ψ) production in the color singlet mechanism whereas there are uncanceled infrared divergences due to eikonal gluons exchange in case of P-wave heavy quarkonium (χ cJ ) production in the color singlet mechanism [16].
Recently we have shown that these uncanceled infrared divergences can be factored into the correct definition of the color singlet P-wave heavy quarkonium non-perturbative matrix element by supplying the eikonal lines or the gauge links [17]. In this paper we will extend this to the non-equilibrium QCD by using the closed-time path integral formalism. We will prove the factorization of infrared divergences in the χ cJ production from the color singlet cc pair in non-equilibrium QCD at RHIC and LHC at all orders in coupling constant. We will predict the correct definition of the non-perturbative matrix element of the χ cJ production from the color singlet cc pair in non-equilibrium QCD at RHIC and LHC. This can be relevant to detect the QGP at RHIC and LHC.
The paper is organized as follows. In section II a brief discussion on the generating functional in non-equilibrium QCD is presented. In section III we discuss the non-canceling infrared divergences in color singlet χ cJ production. In section IV we show that the infrared divergences due to eikonal gluons exchange can be studied by using the SU(3) pure gauge.
In section V we prove the factorization of infrared divergences in the χ cJ production from color singlet cc pair in non-equilibrium QCD at RHIC and LHC at all orders in coupling constant. In section VI we predict the correct definition of the non-perturbative matrix element of the χ cJ production from color singlet cc pair in non-equilibrium QCD at RHIC and LHC. We conclude in section VII.

II. CLOSED-TIME PATH INTEGRAL FORMALISM AND THE GENERATING FUNCTIONAL IN NON-EQUILIBRIUM QCD
Since we will use the background field method of QCD in this paper we denote the gluon field by Q λd (x) and the background field by A λd (x) where λ = 0, 1, 2, 3 and d = 1, ..., 8.
The generating functional in non-equilibrium QCD (without the background field) in the closed-time path integral formalism is given by [6,7] where δ = 0, 1, 2, 3 and we have included the heavy quark. In eq. (1) the symbol k = 1, 2, 3 = u, d, s stands for up, down and strange quark with mass m k and field ψ k . The heavy quark field is Ψ and the heavy quark mass is M. The initial density of states is denoted by ρ, the arbitrary gauge fixing parameter is α, the |0, corresponds to the state at the initial time and and similarly for the − index where +, − stand for the closed-time path indices. Note that we do not introduce ghost fields as we directly work with the ghost determinant det( in eq. (1).
The corresponding non-equilibrium QCD generating functional in the closed-time path integral formalism of the background field method of QCD is given by [6,7,[18][19][20]] where the background gauge fixing depends on the background field A λd (x). In eq. (3) and we do not have any ghost fields because we directly work with the ghost determinant For the type I gauge transformation we have [19,20] T d A ′λd where the light-like gauge link or the light-like eikonal line in the fundamental representation of SU (3) is given by [1,21,22] Φ where l λ is the light-like four-velocity.
In this paper we will use the generating functionals from eqs. (1) and (3)  The non-canceling infrared divergences were found in the higher order pQCD calculation of the annihilation of heavy quark-antiquark pair to light partons in the hadronic decay of the color singlet P-wave heavy quarkonium [16]. For example, in the partonic processes [16] of the hadronic decay of χ cJ and h c respectively, one finds the non-canceling infrared divergences due to real soft gluons (eikonal gluons) emission/absorption [16,23,24]. Now let us discuss the hadroproduction of χ cJ from color singlet cc pair at high energy colliders. If the factorization theorem is valid [1,21,22,[25][26][27] then the χ cJ production from the color singlet cc pair at high energy colliders is given by where dσ kj→CC[ 3 P J ]+X(P T ) is the partonic level cross section for the cc production in 3 P J state. This partonic level cross section can be calculated by using pQCD where k, j = q,q, g.
The parton distribution function f k/p (x, Q) of the parton k inside the proton p is a nonperturbative quantity in QCD. The non-perturbative matrix element of χ cJ production from the color singlet cc pair is denoted by < 0|O χ cJ |0 >.
As mentioned above the non-canceling infrared divergences were found in the hadronic decay of the color singlet P-wave heavy quarkonium [16,23,24]. Similarly, the non-canceling infrared divergences were also found in the hadroproduction of the color singlet P-wave heavy quarkonium [24].
Note that for S-wave and P-wave color singlet heavy quarkonium the infrared divergences occur due to coulomb gluon and eikonal gluon exchanges. The infrared divergence due to Coulomb gluon exchange is analogous to the infrared divergence due to the Coulomb photon exchange in QED, see [28]. This Coulomb gluon infrared divergence is also known as the 1 v → ∞ divergence where v is the relative velocity of the heavy quark-antiquark which is normally absorbed into the normalization of the bound state wave function [16] similar to that in QED [28].
In case of j/ψ production the infrared divergences due to the eikonal gluons interacting with charm quark exactly cancel with the corresponding infrared divergences associated with the charm antiquark [16]. Hence there is no uncanceled infrared divergences due to eikonal gluons exchange in case of j/ψ production. That is why there are no gauge links in the definition of the j/ψ wave function [17].
However, in case of χ cJ production the non-canceling infrared divergences occur due to the eikonal gluons [16]. At NLO in coupling constant the non-canceling infrared divergence due to the eikonal gluons exchange is found in the quark-antiquark fusion process [24] qq → χ cJ g.
Because of the existence of these non-canceling infrared divergences, we have shown in [17] that the gauge links are necessary in the definition of the color singlet P-wave nonperturbative matrix element of the heavy quarkonium production. These gauge links make the non-perturbative matrix element gauge invariant and cancel these non-canceling infrared divergences.
Hence the correct definition of the non-perturbative matrix element of the χ c0 production from color singlet cc pair at high energy colliders which is consistent with the factorization of infrared divergences at all orders in coupling constant in QCD is given by [17] < 0|O χ c0 |0 >=< 0|ζ † Φ∇Φ † ξa † χ c0 · a χ c0 ξ † Φ∇Φ † ζ|0 > (11) where ζ (ξ) is the two component Dirac spinor field that creates (annihilates) a heavy quark and In eq. (11) the a † χ c0 is the creation operator of the χ c0 , the < 0|O χ c0 |0 > is evaluated at the origin and is the light-like gauge link or the light-like eikonal line in the fundamental representation of SU(3).

IV. INFRARED DIVERGENCE DUE TO EIKONAL GLUON AND THE SU(3) PURE GAUGE BACKGROUND FIELD
As mentioned earlier the real gluon emission/absorption is the source of the non-canceling infrared divergences in case of P-wave heavy quarkonium production/deacy [16,23,24].
In this section we will briefly discuss the infrared divergence due to real gluon emission/absorption which can be described by eikonal Feynman rules in QCD. Let us first discuss the eikonal Feynman rules in QED before proceeding to QCD as the eikonal Feynman rules in QCD is similar to that in QED.
In QED the Feynman diagram contribution for an electron emitting a real photon is given where r λ (k λ ) is the momentum of electron (photon). Eq. (14) has both eikonal part r · ǫ(k) r · k u(r) → ∞ when k λ → 0 (15) and the non-eikonal part The eikonal part is the source of the infrared divergence as eq. (15) diverges in the infrared limit k λ → 0. The non-eikonal part in eq. (16) does not diverge in the infrared limit k λ → 0. This implies that the infrared divergence due to the emission of real photon from the electron can be studied by using only the eikonal term r·ǫ(k) r·k u(r) without taking into account the non-eikonal term k ǫ(k) 2r·k u(r) in the Feynman diagram contribution in eq. (14). Now we will show that the study of the infrared divergences due to the eikonal photons at all order in coupling constant in QED can be enormously simplified when the electron is light-like (r 2 = 0). The effective lagrangian density of the photon in the presence of current density K λ (x) in quantum field theory is given by [1] where the generating functional Z[K] in the path integral formulation involving the photon field Q λ (x) is given by From eq. (15) the eikonal contribution gives the eikonal current density where l λ is the light-like four-velocity (l 2 = 0) of the electron.
Using eq. (20) in (17) we find that From eq. (21) we find that the light-like eikonal current produces pure gauge field in quantum field theory at all space-time points except at the positions perpendicular to the direction of motion of the charge at the time of closest approach, a result which agrees with the classical mechanics [26,30,31].
Hence we find from eq. (21) that the calculation of infrared divergences due to the real photons emission from the light-like electron can be simplified by using the pure gauge field in QED. This can also be seen from Grammer-Yennie approximation [29] as follows. We write the photon polarization as the sum of the transverse (physical) polarization plus the longitudinal (pure gauge) polarization to find [29] ǫ λ (k) = ǫ λ physical (k) + ǫ λ pure gauge (k) (22) where which contributes to the physical (finite) cross section and which does not contribute to the physical (finite) cross section but contributes to the the infrared divergence. This can be explicitly seen by using eq. (22) in the eikonal part in eq.
From eq. (16) the non-eikonal contribution gives the non-eikonal current density where r λ is light-like (r 2 = 0) or non-light-like (r 2 = 0) momentum of the electron. Using eqs. (20) and (30) in eq. (17) we find that the interaction between the (light-like or non-lightlike) non-eikonal line with four-momentum r λ and the gauge field generated by the light-like eikonal line with four-velocity l λ (l 2 = 0) gives the interaction (effective) lagrangian density From eq. (31) we find that, in quantum field theory, the interaction between the noneikonal line and the gauge field generated by the light-like eikonal line does not contribute to the interaction (effective) lagrangian density. Since the light-like eikonal line produces pure gauge field in quantum field theory (see eq. (21)) we find from eqs. (31) and (28) that the light-like eikonal line does not modify the finite physical cross section.
Hence we find from eqs. (21), (31), (25), (26), (27) and (28) that the study of infrared divergences in QED due to real photons emission from the light-like electron can be enormously simplified by using the pure gauge field without modifying the finite value of the cross section.
We have shown in eqs. (21) and (31) that the light-like electron produces pure gauge field in QED. This result in QED agrees with classical mechanics [26,30,31]. Hence we find that the infrared divergences at all orders in coupling constant due to the real photons emission from the light-like electron in quantum field theory can be studied by using the path integral formulation of the background field method of quantum field theory in the presence of pure gauge background field [1,21,22,32,33].
In QED the U(1) pure gauge field A λ (x) is given by A λ (x) = ∂ λ ω(x) and in QCD the SU(3) pure gauge field A λd (x) is given by [1,21,22] as given by eq. (32) [1,21,22]. Note that the path integral technique is suitable to study the properties of the non-perturbative quantities in QCD. It should be mentioned here that the properties of a non-perturbative function may not always be correctly studied by using the perturbative method no matter how many orders of perturbation theory is used. Take, for example, a non-perturbative function The Taylor series at g = 0 gives f (g) = 0 to all all orders in perturbation theory but f (g) = 0 for g = 0.
Having considered the points mentioned above, one should note that perturbative QCD entered a new phase when the cancelation of the leading-order (LO) renormalons between the QCD potential and the pole masses of quark and antiquark was discovered (see for example [34]). Convergence of the perturbative series improved dramatically and much more accurate perturbative predictions became available. Hence, in some later works (see, for example, [35]) it was shown that perturbative predictions in QCD agree well with phenomenological QCD results (determined from heavy quarkonium spectroscopy) and lattice QCD calculations.
For recent developments on color potential produced by the color charge of the quark, see [30,31].
In this paper we will use the path integral formulation of the background field method of QCD to predict the correct definition of the non-perturbative matrix element of the χ cJ production from color singlet cc pair in non-equilibrium QCD which is gauge invariant and is consistent with the factorization of infrared divergences at all orders in coupling constant.
In the closed-time path integral formulation the generating functional in non-equilibrium QCD is given by eq. (1). Hence from eq. (1) we find that the heavy quark-antiquark non-perturbative correlation function of the type < in|Ψ non-equilibrium QCD is given by [6,7,20,36] where r, s = +, − are the closed-time path indices in non-equilibrium QCD (the repeated closed-time path indices r, s in eq. (34) are not summed) and |in > is the ground state in non-equilibrium QCD.
In the closed-time path integral formulation in non-equilibrium the generating functional in the background field method of QCD is given by eq. (3). Hence from eq. (3) we find that the heavy quark-antiquark nonequilibrium-nonperturbative correlation function of the type |in > A in the background field method of QCD is given by [6,7,[18][19][20] From eq. (35) we find where Φ(x) is the light-like gauge link or the light-like eikonal line in the fundamental representation of SU(3) given by eq. (13).

VII. CONCLUSIONS
Recently we have proved the factorization of NRQCD S-wave heavy quarkonium production at all orders in coupling constant. In this paper we have extended this to prove the factorization of infrared divergences in χ cJ production from color singlet cc pair in nonequilibrium QCD at RHIC and LHC at all orders in coupling constant. This can be relevant to study the quark-gluon plasma at RHIC and LHC.