Elliptic flow of hadrons via quark coalescence mechanism using Boltzmann transport equation for Pb+Pb collision at =2.76 TeV
Abstract
Elliptic flow of hadrons observed at relativistic heavyion collision experiments at Relativistic HeavyIon Collider (RHIC) and Large Hadron Collider (LHC), provides us an important signature of possible deconfinement transition from hadronic phase to partonic phase. However, hadronization processes of deconfined partons back into final hadrons are found to play a vital role in the observed hadronic flow. In the present work, we use coalescence mechanism also known as Recombination (ReCo) to combine quarks into hadrons. To get there, we have used Boltzmann transport equation in relaxation time approximation to transport the quarks into equilibration and finally to freezeout surface, before coalescence takes place. A BoltzmannGibbs Blast Wave (BGBW) function is taken as an equilibrium function to get the final distribution and a powerlike function to describe the initial distributions of partons produced in heavyion collisions. In the present work, we try to estimate the elliptic flow of identified hadrons such as , , etc., produced in Pb+Pb collisions at = 2.76 TeV at the LHC for different centralities. The elliptic flow () of identified hadrons seems to be described quite well in the available range. After the evolution of quarks until freezeout time, has been calculated using BTERTA, the approach used in this paper consists of combining two or more quarks to explain the produced hadrons at intermediate momenta regions. The formalism is found to describe elliptic flow of hadrons produced in Pb+Pb collisions to a large extent.
pacs:
25.75q,12.38.Mh, 25.75.Ld, 25.75.DwI Introduction
Study of collective phenomena, among many contemporary signatures of quarkgluonplasma (QGP), continues to remain in the forefront of scientific investigations. While, on the experimental front, analysis of final hadrons’ data from RHIC and LHC experiments has enabled us to look back in time and reconstruct the flow phenomena, the phenomenological models using theoretical and numerical techniques have been able to simulate the events, starting from the point of collision of heavy ions, until the freezeout. The theoretical results have been successful in explaining experimental data to a large extent. With the advent of new techniques, the time is however ripe to be able to resolve differences in theories and experiments, and precisely determine various observables of QGP.
Earlier, attempts were made through extensive theoretical modelling and analysis of data, to reconstruct azimuthal anisotropy, (also known as elliptic flow) , of hadrons in transverse momentum plane (1). However, it is believed that azimuthal anisotropy would develop at the early phase of heavy ion collision, when bulk of the deconfined quarks and gluons from noncentral collision between two heavy ions, goes into local thermalization or quark gluon plasma (QGP) state. The geometrical asymmetry of the spatial overlap zone is transformed into momentum anisotropy of the produced particles. With the onset of the local themalization of the bulk partonic matter, the azimuthal anisotropy (mathematically, the second coefficient of the Fourier expansion of particle transverse momentum spectrum), is exhibited strongly in the collective behaviour of the quarkgluonplasma. This information on initial anisotropy is carried till freezeout, and finally reflected in the hadron spectra. However, the rapid expansion of the medium towards isotropization may smear this information to some extent. But on the other hand, hadronic medium effects may add to the partonic flow until kinetic freezeout sets in. Thus it is necessary to develop robust calculation as to able to discern factors emanating from various phases of heavy ion collision, which may affect particles’ flow. While phenomenologies of hadronic matter try to reconstruct the from the final hadronic spectra, initial anisotropy in the partons’ configuration space on the other hand affects the formation of flow and is calculated using phenomenological models such as Glauber mechanism along with perturbative QCD based calculations. However, the two phases of initial anisotropy and hadrons’ interaction remain separated by QGP phase. As mentioned earlier, QGP phase contributes to the evolution of particle flow to a great extent. Hence, It is up to the transport models which may properly bring in the QGP effects and bridge the initial anisotropy and effects of hadronic phase in the observed .
The available transport calculations are based on either hydrodynamical equations, Langevin equation, or Boltzmann transport equation. In this paper, we have used Boltzmann transport equation (BTE) in relaxation time approximation (RTA). BTERTA would transport the entire parton distribution to equilibration and then to kinetic freezeout whereafter any kinetic interaction among particles ceases completely. We will discuss our approach in detail in subsequent sections. Next, we have assumed that the final quark spectra would hadronize at the freezeout surface, using partonic coalescence mechanism. We will return to this formalism in detail in one of the formalism sections. Finally, elliptic flow, , for various hadrons is calculated and presented in the results and discussion section. We have then concluding section of our paper, which is followed by bibliography.
Let us now discuss BTERTA formalism briefly.
Ii Relaxation time approximation (RTA) of Boltzmann transport equation (BTE)
As mentioned in the introductory section, the evolution of quarks within the medium towards the freezeout surface have major effects on the observed final particle spectra. The transport calculations such as hydrodynamics, BTE etc. are commonly used as the evolution mechanisms and provide description of hadron spectra in both qualitative and quantitative manner (2); (3); (4); (5); (6); (7); (8). We know that various dynamical features ranging from multiparton interaction, inmedium energy loss, thermal, and chemical equilibrations, to dynamics at freezeout surfaces contribute extensively to the particle flow and can be studied using BTE. We also know that partons evolving through space and time undergo several collisions and thermalize. Furthermore, they continue to evolve and expand until freezeout even after hadronization. Any of these features can be studied using BTE.
The BTE in general can be written as:
(1) 
where is the distribution of particles which depends on position, momentum and time. is the velocity and is the external force. and are the partial derivatives with respect to position and momentum, respectively. is the collision term which depicts the interaction of the particles with the medium or among themselves. Earlier, BTE has also been used in RTA to study the time evolution of temperature fluctuation in a nonequilibrated system (9) and also for studying the and of various light and heavy flavours at RHIC and LHC energies (10); (11).
We have assumed and zero external forces, 0, the second and third terms of eq. (1) become zero and it reduces to,
(2) 
In RTA (12), the collision term is expressed as:
(3) 
where is Boltzmann local equilibrium distribution characterized by a freezeout temperature . is the relaxation time, the time taken by a nonequilibrium system to reach equilibrium. Using eq. (3), eq. (2) becomes
(4) 
Solving the above equation with the initial conditions i.e. at and at , in general, we get final distribution for any quark flavour as:
(5) 
where is the freezeout time. We use eq. (5) in the definition of the elliptic flow () at midrapidity, which is expressed as,
(6) 
Eq. (6) gives azimuthal anisotropy after incorporating RTA in BTE. It involves a powerlaw like function as the initial distribution of quarks and antiquarks, and BoltzmannGibbs Blast Wave (BGBW) function as their equilibrium distribution. Here, we take BGBW function, , as:
(7) 
where the particle fourmomentum is, , the fourvelocity denoting flow velocities in spacetime is given by, , while the kinetic freezeout surface is given by . Here, is the spacetime rapidity. Assuming boostinvariant scenario where we have taken Bjorken correlation in rapidity, (13) along longitudinal or beam axis. Thus, eq. (II) can be expressed as
(8) 
where . Here is the quark degeneracy factor, is the particle emission time , and is the transverse mass.
in the integrand is a transverse rapidity variable which is given by , with as a function of impact parameter, and gives the anisotropy dependence in the flow. (14); (15); (16); (17) is the radial flow, where is the maximum surface velocity and , with as the radial distance from the center of the fireball. In the blastwave model, the particles closer to the center of the fireball move slower than the ones at the edges. The average of the transverse velocity can be evaluated as (18),
(9) 
While the anisotropic parameter, is written as,
(10) 
In our calculation, we use a linear velocity profile, (), is the maximum radius of the expanding source at freezeout (), and is the radius of colliding nucleus. is the impact parameter to include the centrality dependence of anisotropy. In this paper, we have parametrized the initial distribution given by particle production using perturbative QCD leading order (pQCD LO) calculations for collision,
(11) 
Here, and are momentum fractions carried by interacting partons from their respective colliding protons and are given by,
(12) 
The parton density functions, , are taken to be CTEQ5M (27). The partonic differential scattering crosssections, is calculated from the LO processes, and . To incorporate NLO processes, we have taken a factor, , of value 2.5, and finally nuclear overlap function, and EKS98 parametrization for shadowing effects are taken into account to convert particle production crosssection from collision into , particle spectra. The obtained distribution is parametrized having a function with a powerlike structure and we fixed the parameters of the preequilibrated partons shown in Table 1. However, it is worthwhile to mention that other types of functions can be utilized to obtain the initial quark distributions.
Here too, we have assumed Bjorken correlation in rapidity, . Using Glauber model, is calculated to be for 05% centrality, and for 5060% centrality of the colliding nuclei at = 2.76 TeV. Using eqs. (8) and (II), the final distribution can be expressed as in eq. (5). This gives the final distribution for quarks. The quark masses for the initial distributions is taken to be, MeV, MeV, MeV, and GeV. Then the quarks recombine into hadrons by coalescence formalism, which we will discuss in the following section.
Iii Quark coalescence
The quark coalescence model (ReCo) is used to recombine quarks into hadrons and is found to be one of the prominent hadronization mechanisms beside parton fragmentation (19); (20). The coalescence or recombination of partons into hadrons are able to explain experimentally observed hadron spectra in the intermediate and perhaps at the low momentum regions, while parton fragmentation processes are aptly suitable in explaining hadrons with high momenta. Thus ReCo mechanism may provide another indirect evidence of partonic degrees of freedom in relativistic heavy ion collision.
The coalescence model can be applied to the quarks at the freezeout surface when two (three) quarks recombine to from mesons (baryons) (21); (23); (22). The model can be further utilized in describing observed spectra of light nuclei such as deuteron which contains a neutron and a proton (24).
The coalescence model combines two or more quark distributions using convoluting functions also known as Wigner functions. The basic equation showing the number of mesons from two combining quarks can be broadly written as,
where, and are the spatial and transverse momentum coordinates of the combining quarks and antiquarks, are the quark distribution functions. is the Wigner function convoluting two partonic distributions. in the front of eq. III, is meson degeneracy factor. We have also assumed Bjorken correlation in rapidities, and , throughout. We also assumed = 0.
We have assumed the delta functions correlation, , and . Let us define the partons in the spatial and momentum coordinates in the C.M. frame of meson, such as ,
(15) 
so that we can derive,
(16) 
Here we have also assumed that is small compared to , and thus neglected in the quark distributions, . Thus we have,
(17) 
We have now,
(18) 
where, meson transverse mass factor is given by, .
As for the Wigner function, , we can use the following relation,
(19) 
Therefore eq. 18 is transformed as,
(20) 
To simplify our equations, we convert our momentum variable into lightcone coordinates, of the interacting quarks in the momentum space of the hadron as follows,
(21) 
We also assume that the partons recombining into hadrons have their momenta almost parallel to the final hadron. So can be considered to be very small compared to and it’s dependency in the quark distribution, has also been neglected. It can be shown following eq. III, the parton momentum, , where , is the momentum fraction of the final hadron’s momentum, carried by its constituent quarks during recombination (25); (26). Putting the above conditions into the equation, and assuming the normalization,
(22) 
The final equation can be written,
(23) 
We have along the unit normal direction, at the freezeout hypersurface. The masses of the quarks at the freezeout surface are taken to be quark constituent mass with MeV, MeV, and GeV, respectively. Similarly, for the baryons, one can derive to show,
(24) 
To illustrate on our calculations, one may use eq. (8) as an example and show that
Thus one may calculate to show,
(26) 
Similarly for the baryons, one may write,
We have replaced transverse mass, , by expressions from eqs. (26) and (III), throughout our calculations. We have also assumed a general Gaussian distribution as Wigner functions, (for mesons), (for baryons), which are given by,
(28) 
Here is the width of the Gaussian function, and its small values would give us the narrow Wigner function closer to being a delta function or on other hand, its larger values would give us broad convoluting function instead. The values can be chosen according to the best fit with the particle spectra. We will resume its discussion in the results section.
(29) 
quark flavour  (GeV)  

5.615  1.127  3.73376 10  
5.999  1.099  8.73376 10  
5.579  1.434  3.6286 10  
5.953  1.401  9.1286 10  
6.523  1.892  2.6317 10  
7.250  3.287  2.32815 
Iv Results and Discussions
We would like to reiterate that in the current work, we have transported quarks of various flavours ( and ) produced promptly from initial gluon fusion, to freezeout time and then recombined them into hadrons using coalescence mechanism. We have neglected the decay contributions to final hadron spectra. We have presented the results on the elliptic flow () of various identified hadrons like pions, kaons, protons, D meson, lambda etc. for different centralities of Pb+Pb collision at = 2.76 TeV. While analysing the data, we have kept the freezeout temperature () for the hadrons at 0.095 GeV for the most central collisions (05)% and 0.11GeV for most peripheral collisions (5060)% (29). We assume that the value of is smaller for the central collisions in comparison to the peripheral collisions. The above assumption on freezeout temperature is based on the fact that the freezeout in peripheral collisions occurs quicker than in the most central collisions (29). The dependence of identified hadrons’ are studied by varying two parameters, and , using eq. 29. Based on the closest explanation of the data, we have kept the width of the Wigner function, fixed at 0.0009 for mesons and 0.04 for baryons in our calculations.
In fig. 1, we have shown spectra of various charged hadrons in the most central collisions of Pb+Pb at = 2.76 TeV. The coalescence method is employed to form hadrons from quarks at the freezeout surface. The resulting transverse momentum distributions are then drawn and compared with the experimental data from ALICE@CERN. It is found that, the discussed model in the above section explains the experimental data in the moderate region.
In fig. 2, we have shown of (). The left plot shows the variation of with for different surface velocity parameter, , while the right plot shows for different . Three different values of , keeping fixed are taken and vice versa. Generally speaking, our theoretical results match with the experimental data within errors, from the mid region to the max. shown. However, the model fails to explain the data for 1.0 GeV/c. The reason may be due to the absence of pions from decays of resonances (30). Pions also stand out as an example that shows coalescence picture should work mostly in the mid region.
In fig. 3, the elliptic flow of pions () is presented as a function of for various centralities at = 2.76 TeV for Pb+Pb collisions. Symbols are the experimental data and lines are model results. Here freezeout temperatures () are taken smaller for most central collision than to peripheral collisions. The model results are found to explain the data qualitatively beyond = 1 GeV/c for all the centralities within errorbar. However, the quarkcoalescence mechanism is not able to explain the data below = 1 GeV/c. Experimentally, it is found that for (5060)% appears to be inverse in order compared to (4050)% due to statistical fluctuations. However, the model follows the expected trend of higher for higher centralities.
In the left panel of fig. 4, we have shown elliptic flow or azimuthal anisotropy and spatial anisotropy of the pions versus . is generally defined in terms of spatial coordinates of participants nucleons in the transverse plane. It can be written as:
(30) 
In this paper, GlauberMC formalism (31) has been employed to calculate . Both and decrease with , which is expected. In the right panel of fig. 4, we show the ratio of and vs. or centrality. We find that the ratio tends to increase towards central collisions but drops suddenly for most central. This ratio approximately shows the strength of anisotropy developed as we move towards central collisions and may indicate the extent of collectivity undertaken by the bulk of the partons within quark gluon plasma. However, the sudden drop in this ratio at the most central will be investigated further in our future reports.
In fig. 5, we have presented the variations of of Kaons, () with for 5060 % centrality. The left panel is versus for different at constant , while the right panel shows versus for different at constant . Three different values of , keeping fixed are taken and vice versa. The theoretical curves tend to overestimates the data although it gives a consistent explanation as to the nature of shape of Kaons shown by the experimental data. Also, the plot on the left side shows the theoretical lines cross each other for the different values of , which shows greater sensitivity of on the surface velocity of the fireball. The theoretical line is quite close to experimental points at low which shows that large mass should have less contribution from resonance decays.
In fig. 6, we have shown of short (). The left plot shows the variation of with taking various values of . The right plot of the figure represents the variation of with taking different values of . Three different values of , keeping fixed are taken and vice versa. is a little heavier than Kaons, which is why the and values are almost similar in both the cases. Similarly, the theoretical curve tends to overestimate the data up to = 3 GeV/c. However, theoretical curve shows a gradually increasing trend and slopes down smoothly at high .
Fig. 7 represents the variations of with respect to of phi,. The left plot shows the variation of for different at constant , while the right plot shows the variation of with parameter keeping constant. Three different values of , keeping fixed are taken and vice versa. Phi meson’s results show a gradual rise in the values of with increase in as shown in the plot. Although, the data points show a very small variation after 3 GeV/c, the theoretical curves drop smoothly and continues to do so at = 6.0 GeV/c.
In fig. 8, the elliptic flow of D meson is presented as a function of for centrality 3050% at = 2.76 TeV Pb+Pb collisions. The left panel is shown for various at constant . The model shows rise in for 3 GeV/c and falls smoothly afterwards. The data points show almost a constant value and also number of data points are small to be explained satisfactorily by our model. The right plot is of D meson for different values of keeping constant.
In fig. 9, we have shown the variations of elliptic flow of with respect to for 5060% centrality of Pb+Pb collisions at = 2.76 TeV. In the left hand side of the figure, we show the for different values of keeping fixed. It is found that the model results explain the experimental data qualitatively above = 1 GeV/c for = 0.9. The right hand side of the figure is the results for various at constant . Again there is a good agreement between the model calculations and experimental data above = 1 GeV/c for = 2.2.
In fig. 10, the elliptic flow of is presented with respect to for centrality 5060% at = 2.76 TeV Pb+Pb collisions. In the left hand side of the figure, we show the for different values of keeping fixed. It is found that the model results explain the experimental data qualitatively above = 1 GeV/c for = 0.89. The right hand side of the figure is the results for various at constant . Again there is a good agreement between the model calculations and experimental data above = 1 GeV/c for = 2.4.
In fig 11, we have plotted the of hadron with its three constituent quarks, . Although the flow of the constituent quarks start much before 1.0 GeV/c unlike that of the , the magnitude is much smaller than that of the hadron. Another which is visible from the plot is that the constituent quarks follow some sort of mass ordering with up quark being the lightest has highest flow and strange quark has the lowest. In this calculation and taken from plot are kept fixed for its constituent quarks, .
In fig. 12, the correlation of with is shown for various identified hadrons observed after extracting the values from the model results on elliptic flow with the experimental data at = 2.76 TeV for peripheral Pb+Pb collisions. In this plot, we find that with the increase in , the surface velocity, of hadrons decreases. The mesons show this trend separately from the baryons as evident from the figure. Although the ranges of variations in the values of both the parameters are not large, we find a small mass dependence in the correlation as we go from lightest meson towards heavier meson. Similar trend is also being observed for baryons, and .
V Summary and Conclusions
We have used quark coalescence method for hadronization and Boltzmann transport equation in relaxation time approximation to estimate elliptic flow, for the identified hadrons in Pb+Pb collisions at = 2.76 TeV. The important findings are summarised as follows:

The quark coalescence approach is successful in explaining the elliptic flow data in the moderate transverse momentum region. However, it could not explain the data at low .

The present formalism successfully attempts to connect the particle production from prompt interaction of initially produced partons with finally produced hadrons at freezeout. The hadronic medium effects have not been taken into account. Similarly, resonance decays into observed particles particularly in pion sector have been neglected as well.

We have found a correlation between the radial part of the transverse flow and while explaining the spectra in peripheral collisions.

We have also compared elliptic flow of constituent quarks, () with the final hadron, . We find that of each quark is around of the final baryon. This actually verifies the coalescence mechanism used in the present calculations.

Higher mass quarks are found to have a lower as compared to lighter quarks. On the other hand, flow of mesons behave almost similarly in the mid region although their flow parameter, and time ratio, show correlation and a mass dependence. This is evident from the observation of monotonically decreasing flow parameter with time ratio and particle mass. This also shows that azimuthal anisotropy developed in the partonic phase plays a major role in the observed of final hadrons. Similarly, hadronization mechanism as a part of the freezeout dynamics also play a major role in this regard.
We will continue our investigation on particles’ flow with other hadronization mechanisms such as fragmentation and compare with our current coalescence/recombination model, within the framework of BTERTA mechanism.
Acknowledgements
The authors acknowledge the financial supports from ALICE Project No. SR/MF/PS01/2014IITI(G) of Department of Science & Technology, Government of India. ST acknowledges the financial support by DSTINSPIRE program of Government of India.
Footnotes
 Corresponding author:
References
 S. Tripathy, S. K. Tiwari, M. Younus and R. Sahoo, Eur. Phys. J. A 54, 38 (2018).
 R. Baier, P. Romatschke and U. A. Wiedemann, Phys. Rev. C 73, 064903 (2006).
 S. Gavin, Nucl. Phys. A 435, 826 (1985).
 K. Geiger and B. Muller, Nucl. Phys. B 369, 600 (1992).
 D. K. Srivastava and K. Geiger, Phys. Lett. B 422, 39 (1998).
 S. A. Bass, B. Muller and D. K. Srivastava, J. Phys. G 30, S1283 (2004).
 B. Zhang, M. Gyulassy and C. M. Ko, Phys. Lett. B 455, 45 (1999).
 M. Younus, C. E. ColemanSmith, S. A. Bass and D. K. Srivastava, Phys. Rev. C 91, 024912 (2015).
 T. Bhattacharyya, P. Garg, R. Sahoo and P. Samantray, Eur. Phys. J. A 52, 283 (2016).
 S. Tripathy, T. Bhattacharyya, P. Garg, P. Kumar, R. Sahoo and J. Cleymans, Eur. Phys. J. A 52, 289 (2016).
 S. Tripathy, A. Khuntia, S. K. Tiwari and R. Sahoo, Eur. Phys. J. A 53, no. 5, 99 (2017).
 W. Florkowski and R. Ryblewski, Phys. Rev. C 93, 064903 (2016).
 J. D. Bjorken, Phys. Rev. D 27, 140 (1983).
 P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, Phys. Lett. B 503, 58 (2001).
 E. Schnedermann, J. Sollfrank and U. W. Heinz, Phys. Rev. C 48, 2462 (1993).
 P. BraunMunzinger, J. Stachel, J. P. Wessels and N. Xu, Phys. Lett. B 344, 43 (1995).
 Z. Tang et al., Chin. Phys. Lett. 30, 031201 (2013).
 K. Adcox et al. [PHENIX Collaboration], Phys. Rev. C 69, 024904 (2004).
 D. Molnar and S. A. Voloshin, Phys. Rev. Lett. 91, 092301 (2003).
 Z. w. Lin and D. Molnar, Phys. Rev. C 68, 044901 (2003).
 V. Greco, C. M. Ko and P. Levai, Phys. Rev. C 68, 034904 (2003).
 K. J. Sun and L. W. Chen, Phys. Rev. C 95, 044905 (2017).
 Y. He and Z. W. Lin, Phys. Rev. C 96, no. 1, 014910 (2017).
 X. Yin, C. M. Ko, Y. Sun and L. Zhu, Phys. Rev. C 95, no. 5, 054913 (2017).
 T. Peitzmann, Acta Phys. Hung. A 27, 363 (2006).
 G. Kuipers, Advisors: E. Laenen and T. Peitzmann, (2005). https://www.nikhef.nl/pub/theory/masterstheses/gerben_kuipers.pdf
 H. L. Lai, J. Huston, S. Kuhlmann, F. I. Olness, J. F. Owens, D. E. Soper, W. K. Tung and H. Weerts, Phys. Rev. D 55, 1280 (1997).
 B. Abelev et al. [ALICE Collaboration], Phys. Rev. C 88, 044910 (2013).
 B. Abelev et al. [ALICE Collaboration], Phys. Rev. C 88, 044910 (2013).
 V. Greco and C. M. Ko, Phys. Rev. C 70, 024901 (2004).
 M. L. Miller, K. Reygers, S. J. Sanders and P. Steinberg, Ann. Rev. Nucl. Part. Sci. 57, 205 (2007).
 B. B. Abelev et al. [ALICE Collaboration], Phys. Rev. Lett. 111, 222301 (2013).
 B. B. Abelev et al. [ALICE Collaboration], Phys. Lett. B 736, 196 (2014).
 J. Adam et al. [ALICE Collaboration], Phys. Rev. C 95, 064606 (2017).
 B. B. Abelev et al. [ALICE Collaboration], JHEP 1506, 190 (2015).
 J. Adam et al. [ALICE Collaboration], JHEP 1603, 081 (2016).
 B. Abelev et al. [ALICE Collaboration], Phys. Rev. Lett. 111, 102301 (2013).