Low-momentum pion enhancement from schematic hadronization of a gluon-saturated initial state

We study the particle production in the early stage of the ultrarelativistic heavy-ion collisions. To this end the Boltzmann kinetic equations for gluons and pions with elastic rescattering are considered together with a simple model for the parton-hadron conversion process (hadronisation). It is shown that the overpopulation of the gluon phase space in the initial state leads to an intermediate stage of Bose enhancement in the low-momentum gluon sector which due to the gluon-pion conversion process is then reflected in the final distribution function of pions. This pattern is very similar to the experimental finding of a low-momentum pion enhancement in the ALICE experiment at CERN LHC. Relations to the thermal statistical model of hadron production and the phenomenon of thermal and chemical freeze-out are discussed in this context.


Introduction
One of the issues which can be addressed by the kinetic approach is the question of a low-momentum pion enhancement in heavy ion collisions [1]. There are several solutions proposed to explain this effect as, e.g., the hadronization and freeze-out in a chemical non-equilibrium [2][3][4], the separate freeze-out for strange particles [5], Bose-Einstein condensate of pions [6][7][8][9][10], established by elastic rescattering in the final stage [10,11]. However, none of them is commonly accepted yet [8]. We believe, an explanation linked to the presence of non-equilibrium physics and a precursor of pion condensation in heavy ion collisions should be the favorable one, especially after the recent analysis of particle correlations performed by the ALICE collaboration is showing a coherent fraction of charged π−meson emission that is reaching 23% [1,9]. Such formation of Bose condensate is usually described by the introduction of additional non-equilibrium parameters to the statistical approach [10,12], see also [2,8,13].
An alternative scheme may rely on the Boltzmann kinetic equation for gluons and pions with elastic rescattering and a simple model for the parton-hadron conversion process (hadronisation). There are deep physical reasons for the non-equilibrium and pion condensation at the LHC. It can be due to fast expansion and overcooling of the QGP, or due to gluon condensation in the color glass condensate (CGC) initial state preceeding subsequent hadronization of the low-momentum gluons into low-momentum pions [8]. A scenario with an initial state dominated by gluons which subsequently hadronize, eventually via a quarkless evolution through a first order phase transition, has recently been considered in Ref. [14].
In this short communication we investigate the idea that a certain oversaturation of the purely gluonic initial state could lead by elastic rescattering to a precursor of Bose condensation in the gluon sector in the form of a low-momentum gluon enhancement which, however, should be depopulated by the gluon-pion conversion process and thus appear as low-p pion enhancement in the pion sector. The gluon-pion conversion process is assumed with a constant matrix element which may be pictured as the local limit of a quark one-loop diagram for the case of large quark mass (quark confinement). We demonstrate the evolution of the coupled gluon and pion distribution functions in this case within a schematic model of coupled kinetic equations.

Kinetic equation approach to thermalization and hadronization
We start with the kinetic equation in the form of a Boltzmann-Nordheim equation, which for a single particle distribution function f = f ( x, p, t) can be written as and C[ f ] represents the collision integral. In this study we restrict ourselves to the case of a uniform (∂ f /∂ x = 0) system in a non-expanding box ( F = d p/dt = 0), therefore only the explicit time-dependence remains: On the other hand, the collision integral for the 1 + 2 → 3 + 4 process is defined as: so that the Eq. (1) will take the following form: describing elastic scattering of the system of particles of one type, e.g. gluons. Here for the process 1 + 2 → 3 + 4 we define as f i the distribution function of particle i with 4-momentum P i = (E i , p i ), |M| as the transition amplitude of the process, and represents the gain and loss terms in the collision integral. In the current study we consider the distribution function to be isotropic through the whole evolution. Moreover, the matrix elements of all the processes involved are taken to be constant: following [11], where the case of a system of pions was considered. Albeit this work describes the academic study with only constant matrix elements, the ongoing project involving momentum-and angle-dependent transition amplitudes is discussed in the section 4. As we consider an isotropic, uniform, non-expanding system and constant matrix elements for the processes, and taking into account the 4-momentum conservation (P 1 + P 2 = P 3 + P 4 ), the equation (4) takes the form where D = min{p 1 , p 2 , p 3 , p 4 } and p i are now the radial components of the three-momenta. Details of the derivation are shown in the Appendix A. For future investigations it is helpful to rewrite the equation (6) in terms of a momentum integration, as we would like to extend the approach to the angle-dependent collision integral, as well as to non-uniform systems. Therefore, in the current work we use the following formula: Obviously, elastic scattering is necessary but not sufficient to achieve low-p pion enhancement. The second required process which needs to be accounted for is hadronization. In this exploratory work we connect the gluon sector directly with the pion one. For such system there are three contributing channels: ππ → ππ, gg → gg, and gg ↔ ππ. Therefore at the end we have a coupled system of equations: ε 3 ε 4 dp 3 dp 4 DF[ f g ] where M gg→ππ and M ππ→gg are matrices for hadronization channels. Note, that due to the momentum conservation p 2 = p 3 + p 4 − p 1 in Eq. (8). In this study we set M ππ→gg = 0, which is motivated by the threshold for this process due to the large value of the gluon mass: m g = 0.7 GeV. The value of M gg→ππ is set to be constant and should be seen as an academic example.
As the initial condition of the system we take an oversaturated gluon distribution given by a step-like function [15,16] inspired by the CGC picture of the initial state which is assumed to have no pions By Q s we denote the saturation scale. However, in order to avoid numerical problems that would occur with the step-function distribution, we use instead the following smooth function [15] f g (t, p) to define the initial conditions. We keep our model simple and therefore do not introduce an extra timescale for the start of hadronization. However, we keep in mind that the underlying microphysical process is, e.g., a quark-box diagram, which consists of the Breit-Wheeler type process of 2g → qq and subsequent hadronization cross section qq → ππ . In the future we plan to investigate the problem of the gluon-to-pion conversion in detail, for instance within a Nambu-Jona-Lasinio model [17][18][19][20] and/or by exploiting dynamical schemes of hadronization that would address the confinement aspect as well [21][22][23][24].

Results
In Fig. 1 we show the evolution of the gluon distribution function from a CGC motivated initial (over-)saturated gluon state to a thermal distribution due to elastic scattering according to the gg → gg process. The timescale to reach a thermalized final state is of the order of t final ∼ 250 fm/c and thus exceeds the typical duration evolution towards freeze-out of the fireball created in a heavy-ion collision. This is mainly due to the fact that the value of the matrix element taken in this example calculation as |M| = 4.5 is unrealistically small. In Fig. 2 we show the same evolution of the gluon distribution function for three different values of the matrix element. The value M = 140 leads to a thermalization time scale which nicely corresponds to the result of a calculation by Shuryak [25].
When the coupling to the pion sector is switched on, the gluon conversion proceeds and the initially empty pion phase space gets populated at the expense of the gluon one. Due to the relation of the gluon and pion masses the reverse process (the pion annihilation to two gluons) does practically not take place. In Fig. 3 the evolution from the initially pure gluon saturated state to the thermal pion state without gluons is shown. The pion distribution shows clearly the low-momentum enhancement typical for a precursor of Bose condensation. This is the fact observed in the ALICE experiment at CERN for which we wanted to give a qualitative explanation with the simple kinetic model presented here. It should be noted that here we used as a test the equal values for the three transition amplitudes: |M gg→gg | = |M ππ→ππ | = |M gg→ππ | = 4.5.  Our simplified model shows, under the assumption of gluon dominance in the initial state, the quarkless evolution of the system towards a pion gas with low-momentum pion enhancement as a precursor of Bose condensation. According to (8) both particle species (gluons and pions) undergo two main processes: conversion and elastic scattering. Both of them are responsible for low-momentum (low-p) pion enhancement.
The fist process turns π-mesons to gluons and vice versa and its rate is defined by two matrix elements M ππ→gg , M gg→ππ , which in the simplest case considered here are constant numbers. Particle conversion can take place only when energy of incoming particles is at least equal to mass of outgoing ones. Consequently, in case of massless gluons kinematics restricts gg → ππ reaction to higher energy region of a spectrum, making the whole process slow.
The impact of the second process (elastic scattering ) is more subtle. It lowers momentum of particles through subsequent collisions, leading them to "pile-up" near zero momentum mode. The effect is especially strong for bosons due to the statistical factor (1 + f i ) in (8) and allows pre-condensate formation even before thermalization. In normal circumstances for long enough times the distribution should become an equilibrium Bose function. However, in our model massive gluons undergo a complete conversion to pions before thermalization because m g > m π (m g = 0.7 GeV, while m π = 0.14 GeV).

Discussion
The present model, albeit quite simple, shows the formation of the pion condensation precursor emerging from an oversaturated purely gluonic state. The process takes place before the system reaches equilibrium. The model can be improved by the use of non-constant matrix elements and thus taking into account scattering angle in collision kinematics. Such an improved model would allow us to discuss the different scales and their evolution, e.g., the Debye scale, the UV and IR scale, see Refs. [7,11,15,26].
These improved matrix elements should also bear the confining aspects of gluon-gluon interactions which ultimately should be responsible for the absence of gluons from the final state.
The assumption of a constant gluon mass, exceeding the value of the pion mass is a rather schematic realization of this concept which provides ample room for improvement. Here it would be beneficial to make a comparison with the study in Ref. [15], where a system of massless gluons undergoes the evolution due to elastic scattering with similar restrictions as used in the current paper. However, the equation (7) will no longer be valid in the case of non-constant matrix elements and angle-dependence, and thus will need to be rederived.
Another room for advancement lies in direct handling of the kinetics of Bose condensation (see, e.g., Ref. [10]). One way to do that is the separation of the distribution function into two parts: where the first term represents the "gas" and the second describes BEC. This ansatz has been discussed for the oversaturated pion gas in Ref. [11] and recently also for the gluon plasma in Ref. [27]. We hope to achieve manifest energy and particle number conservation with such an improved formulation of the particle kinetics in the presence or precursory development of a Bose condensate in the system. The model can be extended towards a more realistic description of a hadronizing gluon-dominated initial state in high-energy heavy-ion collisions by including more hadronic species as they are observed in those experiments in good agreement with the thermal statistical model [28]. This calls then for an extension of the collision integrals in our kinetic model to other classes of processes than just 2 → 2 processes as, e.g., the three-meson conversion to a baryon-antibaryon pair and its reverse [29].
Last, but not least we want to mention that the assumed absence of dynamical quarks is only a simplifying assumption. In an improved model, their kinetics shall be coupled to that of the gluons and all considered hadron species. Their absence in the final state shall be realised due to a confining mechanism. The one already tested in the framework of a kinetic theory is the Gribov-Zwanziger confinement realized by an infrared-divergent selfenergy [22][23][24]. We shall come back to these issues in a subsequent, more elaborate work on the subject.