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

TL;DR

This paper models early-stage particle production in heavy-ion collisions, showing that gluon overpopulation causes low-momentum pion enhancement, aligning with CERN LHC experimental results.

Contribution

It introduces a schematic hadronization model linking gluon phase space overpopulation to pion enhancement, providing insights into initial state effects in heavy-ion collisions.

Findings

Gluon overpopulation leads to Bose enhancement in low-momentum gluons.

Gluon-pion conversion reflects in low-momentum pion enhancement.

Model aligns with experimental data from CERN LHC ALICE experiment.

Abstract

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.

Figures

Equations

$$\frac{\mathrm{d}f}{\mathrm{d}t}=C[f],$$

$$\frac{\mathrm{d}f}{\mathrm{d}t}=\frac{\partial f}{\partial \vec{x}}\frac{\mathrm{d}\vec{x}}{\mathrm{d}t}+\frac{\partial f}{\partial \vec{p}}\frac{\mathrm{d}\vec{p}}{\mathrm{d}t}+\frac{\partial f}{\partial t}$$

$$C[f(t,{\vec{p}}_1)]=\frac{(2\pi {)}^4}{2{\mathrm{E}}_1}\int {\delta }^4(\sum _i{P}_i)|\mathrm{M}{|}^{2}F[f]\prod _{k=2}^4\frac{{\mathrm{d}}^3{\vec{p}}_k}{(2\pi {)}^{3}2E_k},$$

$${\partial }_{t}f(t,{\vec{p}}_1)=\frac{(2\pi {)}^4}{2{\mathrm{E}}_1}\int {\delta }^4(\sum _i{P}_i)|\mathrm{M}{|}^{2}F[f]\prod _{k=2}^4\frac{{\mathrm{d}}^3{\vec{p}}_k}{(2\pi {)}^{3}2E_k},$$

$$|M{|}_{12\to 34}=\mathrm{const},$$

$${\partial }_{t}f({\epsilon }_1)$$

$${\partial }_{t}f(p_1)$$

$$\begin{array}{rl}\frac{\partial f_{\pi }}{\partial t}(t,{\vec{p}}_1)& =\int \int \frac{{\left|M_{\pi \pi \to \pi \pi }\right|}^2}{64{\pi }^3{\epsilon }_1}\frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}F[f_{\pi }]\\ & +\left(1+f_{\pi }(t,p_1)\right)\int \int \frac{{\left|M_{gg\to \pi \pi }\right|}^2}{64{\pi }^3{\epsilon }_1}\frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}\left(1+f_{\pi }(t,p_2)\right)f_g(t,p_3)f_g(t,p_4)\\ & -f_{\pi }(t,p_1)\int \int \frac{{\left|M_{\pi \pi \to gg}\right|}^2}{64{\pi }^3{\epsilon }_1}\frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}{f}_{\pi }(t,p_2)\left(1+f_g(t,p_3)\right)\left(1+f_g(t,p_4)\right)\end{array}$$

$$\begin{array}{rl}\frac{\partial f_g}{\partial t}(t,{\vec{p}}_1)& =\int \int \frac{{\left|M_{gg\to gg}\right|}^2}{64{\pi }^3{\epsilon }_1}\frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}F[f_g]\\ & +\left(1+f_g(t,p_1)\right)\int \int \frac{{\left|M_{\pi \pi \to gg}\right|}^2}{64{\pi }^3{\epsilon }_1}\frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}\left(1+f_g(t,p_2)\right)f_{\pi }(t,p_3)f_{\pi }(t,p_4)\\ & -f_g(t,p_1)\int \int \frac{{\left|M_{gg\to \pi \pi }\right|}^2}{64{\pi }^3{\epsilon }_1}\frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}{f}_g(t,p_2)\left(1+f_{\pi }(t,p_3)\right)\left(1+f_{\pi }(t,p_4)\right)\end{array}$$

$$\displaystyle f_{\pi}(t,p)\bigg{|}_{t=0}$$

$$\displaystyle f_{g}(t,p)\bigg{|}_{t=0}=f_{0}\,\left[\theta(1-p/Q_{s})+\theta(p/Q_{s}-1)e^{-a\,(p/Q_{s}-1)^{2}}\right],\quad a=10$$

$${\tilde{f}}_{\pi }(p)=f_{\pi }(p)+(2\pi {)}^3{n}_c^{\pi }\delta (p)$$

$${\tilde{f}}_g(p)=f_g(p)+(2\pi {)}^3{n}_c^g\delta (p)$$

$$C_1[f]=\frac{(2\pi {)}^4}{2{\mathrm{E}}_1}\int {\delta }^4(\sum _i{P}_i)|\mathrm{M}{|}^{2}F[f]\prod _{k=2}^4\frac{{\mathrm{d}}^3{\vec{p}}_k}{(2\pi {)}^{3}2E_k}$$

$${\delta }^3(\sum {\vec{p}}_i)=\int \exp (i(\vec{\lambda },{\vec{p}}_1+{\vec{p}}_2-{\vec{p}}_3-{\vec{p}}_4))\cdot \frac{{\mathrm{d}}^3\vec{\lambda }}{(2\pi {)}^3},$$

$$\mathrm{d}{\vec{p}}_i=\mathrm{d}{\phi }_i\mathrm{d}\cos {\theta }_i{p}_i^2\mathrm{d}{p}_i={\epsilon }_i{p}_i\mathrm{d}{\Omega }_i\mathrm{d}{\epsilon }_i$$

$$C_1[f]=\frac{{\left|M\right|}^2}{64{\pi }^3{\epsilon }_1}\int \delta ({\epsilon }_1+{\epsilon }_2-{\epsilon }_3-{\epsilon }_4)\mathrm{D}F[f]\mathrm{d}{\epsilon }_3\mathrm{d}{\epsilon }_4\mathrm{d}{\epsilon }_2,$$

$$\mathrm{D}=\frac{p_2{p}_3{p}_4}{64{\pi }^5}\int {\lambda }^2\mathrm{d}\lambda \int e^{i({\vec{p}}_1,\vec{\lambda })}\mathrm{d}{\Omega }_{\lambda }\int e^{i({\vec{p}}_2,\vec{\lambda })}\mathrm{d}{\Omega }_2\int e^{i({\vec{p}}_3,\vec{\lambda })}\mathrm{d}{\Omega }_3\int e^{i({\vec{p}}_4,\vec{\lambda })}\mathrm{d}{\Omega }_4.$$

$$\displaystyle\begin{split}\int e^{i(\vec{p}_{1},\vec{\lambda})}\mathrm{d}\Omega_{\lambda}&=\int e^{i(p_{1}\lambda\cos\theta_{\lambda})}\mathrm{d}\Omega_{\lambda}=\int_{0}^{2\pi}\mathrm{d}\varphi\int_{-1}^{1}\mathrm{d}\cos\theta e^{i(p_{1}\lambda\cos\theta_{\lambda})}=\\ &=\frac{2\pi}{ip_{1}\lambda}e^{ip_{1}\lambda x}\big{|}^{x=1}_{x=-1}=\frac{2\pi}{p_{1}\lambda}\frac{e^{ip_{1}\lambda}-e^{-ip_{1}\lambda}}{2i}\cdot 2=\frac{4\pi}{p_{1}\lambda}\sin(p_{1}\lambda)\end{split}$$

$$\mathrm{D}=\frac{4}{\pi p_1}\int \frac{\mathrm{d}\lambda }{{\lambda }^2}\sin (p_1\lambda )\sin (p_2\lambda )\sin (p_3\lambda )\sin (p_4\lambda )$$

$$\begin{array}{rl}& -\sqrt{\frac{\pi }{2}}w\mathrm{Sign}(w)=\frac{1}{\sqrt{2\pi }}{\int }_{-\infty }^{\infty }\frac{1}{x^2}{e}^{iwx}\mathrm{d}x=\\ & \frac{1}{\sqrt{2\pi }}\left({\int }_0^{\infty }\frac{1}{x^2}{e}^{iwx}\mathrm{d}x+{\int }_0^{\infty }\frac{1}{(-x{)}^2}{e}^{-iwx}\mathrm{d}x\right)=\frac{1}{\sqrt{2\pi }}{\int }_0^{\infty }\frac{1}{x^2}\left(e^{iwx}+e^{-iwx}\right)\mathrm{d}x,\end{array}$$

$$\displaystyle\begin{split}\mathrm{D}&=\frac{4}{\pi p_{1}}\int^{\infty}_{0}\frac{\mathrm{d}\lambda}{\lambda^{2}}\frac{e^{ip_{1}\lambda}-e^{-ip_{1}\lambda}}{2i}\frac{e^{ip_{2}\lambda}-e^{-ip_{2}\lambda}}{2i}\frac{e^{ip_{3}\lambda}-e^{-ip_{3}\lambda}}{2i}\frac{e^{ip_{4}\lambda}-e^{-ip_{4}\lambda}}{2i}=\\ &=\frac{1}{4\pi p_{1}}\int^{\infty}_{0}\frac{\mathrm{d}\lambda}{\lambda^{2}}\Big{(}(e^{i\lambda(p_{1}+p_{2})}-e^{i\lambda(p_{2}-p_{1})}--e^{i\lambda(p_{1}-p_{2})}+e^{i\lambda(-p_{1}-p_{2})})(e^{i\lambda(p_{3}+p_{4})}-\\ &-e^{i\lambda(p_{4}-p_{3})}-e^{i\lambda(p_{3}-p_{4})}+e^{i\lambda(-p_{3}-p_{4})})\Big{)}=\\ &=\frac{1}{4\pi p_{1}}\int^{\infty}_{0}\frac{\mathrm{d}\lambda}{\lambda^{2}}\Big{(}e^{i\lambda(p_{1}+p_{2}+p_{3}+p_{4})}+e^{-i\lambda(p_{1}+p_{2}+p_{3}+p_{4})}-e^{i\lambda(p_{3}+p_{4}+p_{2}-p_{1})}\\ &-e^{-i\lambda(p_{3}+p_{4}+p_{2}-p_{1})}-e^{i\lambda(p_{3}+p_{4}+p_{1}-p_{2})}-e^{-i\lambda(p_{3}+p_{4}+p_{1}-p_{2})}+e^{i\lambda(p_{3}+p_{4}-p_{1}-p_{2})}+\\ &+e^{-i\lambda(p_{3}+p_{4}-p_{1}-p_{2})}-e^{i\lambda(p_{1}+p_{2}+p_{4}-p_{3})}-e^{-i\lambda(p_{1}+p_{2}+p_{4}-p_{3})}+e^{i\lambda(p_{4}-p_{3}+p_{2}-p_{1})}+\\ &+e^{-i\lambda(p_{4}-p_{3}+p_{2}-p_{1})}+e^{i\lambda(p_{4}-p_{3}+p_{1}-p_{2})}+e^{-i\lambda(p_{4}-p_{3}+p_{1}-p_{2})}-e^{i\lambda(p_{4}-p_{3}-p_{1}-p_{2})}\\ &-e^{-i\lambda(p_{4}-p_{3}-p_{1}-p_{2})}\Big{)}=\\ &=\frac{1}{4\pi p_{1}}(-\pi)(|p_{1}+p_{2}+p_{3}+p_{4}|-|p_{3}+p_{4}+p_{2}-p_{1}|-|p_{3}+p_{4}+p_{1}-p_{2}|\\ &+|p_{3}+p_{4}-p_{1}-p_{2}|-|p_{1}+p_{2}+p_{4}-p_{3}|+|p_{4}-p_{3}+p_{2}-p_{1}|\\ &+|p_{4}-p_{3}+p_{1}-p_{2}|-|p_{4}-p_{3}-p_{1}-p_{2}|)=\\ &=\frac{1}{4\pi p_{1}}(-\pi)(-4\min\{p_{1},p_{2},p_{3},p_{4}\})=\frac{\min\{p_{1},p_{2},p_{3},p_{4}\}}{p_{1}}\end{split}$$

$$C[f({\epsilon }_1)]=\frac{{\left|M\right|}^2}{64{\pi }^3{\epsilon }_1}\int \int \mathrm{d}{\epsilon }_3\mathrm{d}{\epsilon }_4\mathrm{D}F[f]$$

$$C[f(p_1)]=\frac{{\left|M\right|}^2}{64{\pi }^3{\epsilon }_1}\int \int \frac{p_3{p}_4}{{\epsilon }_3{\epsilon }_4}\mathrm{d}{p}_3\mathrm{d}{p}_4\mathrm{D}F[f].$$

Full text

\pubvolume

xx \issuenum1 \articlenumber5

\historyReceived: date; Accepted: date; Published: date \TitleLow-momentum pion enhancement from schematic hadronization of a gluon-saturated initial state \AuthorElizaveta Nazarova1,2\orcidA, Łukasz Juchnowski2, David Blaschke2,3,4,∗ and Tobias Fischer2 \AuthorNamesElizaveta Nazarova, Lukasz Juchnowski, David B. Blaschke, Tobias Fischer

\corresCorrespondence: david.blaschke@gmail.com

1 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 Abelev et al. (2014). There are several solutions proposed to explain this effect as, e.g., the hadronization and freeze-out in a chemical non-equilibrium Kataja and Ruuskanen (1990); Gavin and Ruuskanen (1991); Begun (2015), the separate freeze-out for strange particles Begun et al. (2014a), Bose-Einstein condensate of pions Begun et al. (2014b); Blaizot et al. (2012); Begun (2016); Begun and Florkowski (2015); Voskresensky (1996), established by elastic rescattering in the final stage Voskresensky (1996); Semikoz and Tkachev (1997). However, none of them is commonly accepted yet Begun (2016). 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 $\pi-$meson emission that is reaching 23% Abelev et al. (2014); Begun and Florkowski (2015). Such formation of Bose condensate is usually described by the introduction of additional non-equilibrium parameters to the statistical approach Voskresensky et al. (1995); Voskresensky (1996), see also Kataja and Ruuskanen (1990); Begun (2016); Koch and Rafelski (1986).

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 Begun (2016). 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. Stoecker et al. (2016).

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.

2 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(\vec{x},\vec{p},t)$ can be written as

[TABLE]

where

[TABLE]

and $C[f]$ represents the collision integral. In this study we restrict ourselves to the case of a uniform ($\partial f/\partial\vec{x}=0$) system in a non-expanding box ($\vec{F}=\mathrm{d}\vec{p}/\mathrm{d}t=0$), therefore only the explicit time-dependence remains: $\mathrm{d}f/\mathrm{d}t=\partial f/\partial t\equiv\partial_{t}f$.

On the other hand, the collision integral for the $1+2\rightarrow 3+4$ process is defined as:

[TABLE]

so that the Eq. (1) will take the following form:

[TABLE]

describing elastic scattering of the system of particles of one type, e.g. gluons. Here for the process $1+2\rightarrow 3+4$ we define as $f_{i}$ the distribution function of particle $i$ with 4-momentum $P_{i}=(E_{i},\vec{p}_{i})$, $|M|$ as the transition amplitude of the process, and $F[f]=(1+f_{1})(1+f_{2})f_{3}f_{4}-f_{1}f_{2}(1+f_{3})(1+f_{4})$ 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:

[TABLE]

following Semikoz and Tkachev (1997), 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

[TABLE]

where $\mathrm{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:

[TABLE]

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: $\pi\pi\rightarrow\pi\pi$, $gg\rightarrow gg$, and $gg\leftrightarrow\pi\pi$. Therefore at the end we have a coupled system of equations:

[TABLE]

[TABLE]

where $\mathrm{M}_{gg\rightarrow\pi\pi}$ and $\mathrm{M}_{\pi\pi\rightarrow 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 $\mathrm{M}_{\pi\pi\rightarrow 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 $\mathrm{M}_{gg\rightarrow\pi\pi}$ 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 Blaizot et al. (2013, 2016) inspired by the CGC picture of the initial state which is assumed to have no pions

[TABLE]

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 Blaizot et al. (2013)

[TABLE]

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\to q\bar{q}$ and subsequent hadronization cross section $q\bar{q}\to\pi\pi$ . 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 Rehberg et al. (1999, 1996); Friesen et al. (2014); Marty et al. (2016) and/or by exploiting dynamical schemes of hadronization that would address the confinement aspect as well Feng et al. (2017); Florkowski et al. (2015, 2016a, 2016b).

3 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\to gg$ process. The timescale to reach a thermalized final state is of the order of $t_{\rm final}\sim 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 Shuryak (1992).

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\rightarrow gg}|=|M_{\pi\pi\rightarrow\pi\pi}|=|M_{gg\rightarrow\pi\pi}|=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 $\pi$-mesons to gluons and vice versa and its rate is defined by two matrix elements $M_{\pi\pi\rightarrow gg},~{}M_{gg\rightarrow\pi\pi}$, 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\rightarrow\pi\pi$ 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_{\pi}$ ($m_{g}=0.7$ GeV, while $m_{\pi}=0.14$ GeV).

4 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. Semikoz and Tkachev (1997); Blaizot et al. (2012, 2013); Zhou et al. (2017).

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. Blaizot et al. (2013), 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. Voskresensky (1996)). One way to do that is the separation of the distribution function into two parts:

[TABLE]

[TABLE]

where the first term represents the "gas" and the second describes BEC. This ansatz has been discussed for the oversaturated pion gas in Ref. Semikoz and Tkachev (1997) and recently also for the gluon plasma in Ref. Xu et al. (2015). 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 Andronic et al. (2018). This calls then for an extension of the collision integrals in our kinetic model to other classes of processes than just $2\to 2$ processes as, e.g., the three-meson conversion to a baryon-antibaryon pair and its reverse Seifert and Cassing (2018).

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 Florkowski et al. (2015, 2016a, 2016b). We shall come back to these issues in a subsequent, more elaborate work on the subject.

\funding

This research was funded by the Polish Narodowe Centrum Nauki under grant number UMO-2014/15/B/ST2/03752 (E.N., L.J., D.B.) and grant number UMO-2016/23/B/ST2/00720 (E.N.,T.F.).

Acknowledgements.

We acknowledge discussions with Viktor Begun, Marcus Bleicher, Wojciech Florkowski, Carsten Greiner, Gerd Röpke, Ludwik Turko and Dmitry Voskresensky on the topic of the nonequilibrium nature of the pion distribution emerging from an oversaturated initial state. Further we thank Niels-Uwe Friedrich Bastian for his support in numerical and algebraic questions. Brent Harrison gave a nice presentation of his ongoing work with Andre Peshier on “Bose-Einstein Condensation from a Gluon Transport Equation” at the Symposium on “Nonequilibrium Phenomena in Strongly Correlated Systems” in Dubna, 16. April 2018, which inspired us to prepare this communication. \appendixtitlesyes \appendixsectionsmultiple

Appendix A Collision integral derivation

Here we will simplify the collision integral (3):

[TABLE]

Using the identity:

[TABLE]

and separating the angle integrations:

[TABLE]

the integral takes the following form:

[TABLE]

where D is defined as follows:

[TABLE]

Taking into account that

[TABLE]

we can rewrite D as:

[TABLE]

Using the Fourier transformation:

[TABLE]

we can simplify the integral in the formula for D:

[TABLE]

The last step in this equation (21) (changing to the minimum function between the 4 momenta) can be easily done by checking one of the possibilities, - for example, the case when $p_{1}=\min\{p_{1},p_{2},p_{3},p_{4}\}$ (or $p_{1}<p_{2}<p_{3}<p_{4}$). Taking into account the 4-momentum conservation: $P_{1}+P_{2}=P_{3}+P_{4}$, we get the final result:

[TABLE]

where $D=\frac{1}{p_{1}}\min\{p_{1},p_{2},p_{3},p_{4}\}$. In order to change the formula (22) to the integration over momentum, we can use the connection between energy and momentum in the relativistic case: $\varepsilon^{2}=p^{2}+m^{2}\rightarrow\varepsilon\mathrm{d}\varepsilon=p\mathrm{d}p$, so that the equation (22) takes form:

[TABLE]

\reftitle

References