# Elliptic Flow of Hadrons via Quark Coalescence Mechanism Using the Boltzmann Transport Equation for Pb+Pb Collision atsNN=2.76TeV

Authors: Mohammed Younus; Sushanta Tripathy; Swatantra Kumar Tiwari; Raghunath Sahoo
Journal: Advances in High Energy Physics
 (2020)

---

## Abstract

Elliptic flow of hadrons observed at relativistic heavy ion collision experiments at relativistic heavy ion collider (RHIC) and large hadron collider (LHC) provides us an important signature of possible deconfinement transition from the 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 a coalescence mechanism also known as recombination (ReCo) to combine quarks into hadrons. To get there, we have used the Boltzmann transport equation in relaxation time approximation to transport the quarks into equilibration and finally to freeze-out the surface, before coalescence takes place. A Boltzmann-Gibbs blast wave (BGBW) function is taken as an equilibrium function to get the final distribution and a power-like function to describe the initial distributions of partons produced in heavy ion collisions. In the present work, we try to estimate the elliptic flow of identified hadrons such asπ, K, and p, produced in Pb+Pb collisions at sNN=2.76TeV at the LHC for different centralities. The elliptic flow (v2) of identified hadrons seems to be described quite well in the available pT range. After the evolution of quarks until freeze-out time has been calculated using BTE-RTA, 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 the elliptic flow of hadrons produced in Pb+Pb collisions to a large extent.


---

## Body

## 1. Introduction

The study of collective phenomena, among many contemporary signatures of quark-gluon plasma (QGP), continues to remain in the forefront of scientific investigations [1, 2]. While on the experimental front [3–9] 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 to the freeze-out. 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)v2, of hadrons in a transverse momentum plane [10–13]. 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 the 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 quark-gluon plasma. This information on initial anisotropy is carried till freeze-out 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 freeze-out sets in. Thus, it is necessary to develop robust calculation to be able to discern factors emanating from various phases of heavy ion collision, which may affect the particles’ flow. While phenomenologies of hadronic matter try to reconstruct the v2 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 the Glauber mechanism along with perturbative QCD-based calculations. However, the two phases of initial anisotropy and hadrons’ interaction remain separated by the QGP phase. As mentioned earlier, the QGP phase contributes to the evolution of the 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 the hadronic phase in the observed v2 [14]. The transport models help us in studying collision centrality dependency of QGP properties. They not only shed light on properties of hot and dense matter viz. average temperature and momentum reached by equilibrated system and their dependency on collision centrality but also provide us with vast information on transport properties such as radial flow coefficients, momentum broadening, drag and diffusion coefficients, and electrical and thermal conductivity of QGP matter [15–19].The available transport calculations are based on either the hydrodynamical equation, Langevin equation, or Boltzmann transport equation. In this paper, we have used the Boltzmann transport equation (BTE) in relaxation time approximation (RTA). BTE-RTA would transport the entire parton distribution to equilibration and then to freeze-out surface whereafter kinetic or chemical interaction among particles ceases completely. Using the BTE-RTA approach in the present work, we have attempted to study parameters associated with particle production such as the radial flow and relaxation time of the final-state particles. Neglecting the effects of the hadron medium on the particles’ flow as approximation, the current work focuses on the interplay of the various parameters and mechanisms on the partonic states. We will discuss our approach in detail in subsequent sections. We have also assumed that the final quarks would hadronize into mesons and baryons using the partonic coalescence mechanism at the hadronization hypersurface. We will discuss this formalism in one of the following sections. Finally, the elliptic flow,v2, for various hadrons is calculated and presented in the results and discussion section. We have also presented figures on our study of the parameters and their interdependencies. We have then a concluding section in our paper, which is followed by bibliography.Let us now discuss BTE-RTA formalism briefly.


## 2. Relaxation Time Approximation (RTA) of Boltzmann Transport Equation (BTE)

As mentioned in the introductory section, the evolution of quarks within the medium towards the freeze-out surface has major effects on the observed final particle spectra. The transport calculations such as hydrodynamics and BTE are commonly used as the evolution mechanisms and provide description of hadron spectra in both qualitative manner and quantitative manner [20–26]. We know that various dynamical features ranging from multiparton interaction, in-medium energy loss, thermal and chemical equilibrations, to dynamics at freeze-out 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 freeze-out even after hadronization. Any of these features can be studied using BTE. The BTE in general can be written as
(1)dfx,p,tdt=∂f∂t+v→·∇xf+F→·∇pf=Cf,where fx,p,t is the distribution of particles which depends on position, momentum, and time. v→ is the velocity, and F→ is the external force. ∇x and ∇p are the partial derivatives with respect to position and momentum, respectively. Cf 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 [27] and also for studying the RAA and v2 of various light and heavy flavours at RHIC and LHC energies [28, 29].We have considered the evolution of particle momentum distribution with time. We have taken∇xf=0 assuming particle distribution to be homogeneous in space and the configuration space distribution or spatial distribution has been parametrized accordingly. There are no external forces acting on the system F→=0. Hence, the second and third terms of equation (1) become zero and it reduces to
(2)dfx,p,tdt=∂f∂t=Cf.The full kernel of the collision termCf contains microscopic interaction cross-sections of particles. For any transport models such as AMPT and UrQMD containing microscopic Boltzmann equation, the full interaction kernel along with space and time evolution of the system becomes important. In our calculation owing to assumed homogeneous spatial distribution, the spatial variables have been parametrized.In BTE-RTA [30], which is an effective model, the collision term is however expressed as
(3)Cf=−f−feqτ,where feq is Boltzmann local equilibrium distribution characterized by a freeze-out temperature T. τ is the relaxation time, the time taken by a nonequilibrium system to reach equilibrium. Using equation (3), equation (2) becomes
(4)∂fq∂t=−fq−feqqτ.Solving the above equation with the initial conditions,i.e., at t=0,f=fi and at t=tf,f=ff, in general, we get final distribution for any quark flavour as
(5)ffq=feqq+fiq−feqqe−tf/τ,where tf is the freeze-out time parameter. The initial distribution fi at t=0 is taken as power-like distribution. We will come back to this later. BTE-RTA computes to give the final distribution ff as a function of parameter tf/τ. If the system is given enough time or tf is large compared to τ, ff might converge to feq.We use equation (5) in the definition of the elliptic flow (v2) at midrapidity, which is expressed as
(6)v2qpT=∫ffq×cos2ϕdϕ∫ffqdϕ.Equation (6) gives azimuthal anisotropy after incorporating RTA in BTE. The Boltzmann-Gibbs blast wave (BGBW) function has been taken as the equilibrium distribution function, feq, as
(7)feqqpT=C·exp−pμuμT,dNq/q¯eqpTdpTdy=∫d3σμpμfeqq/q¯pT,where the particle four-momentum is pμ=mTcoshy,pTcosϕ,pTsinϕ,mTsinhy and the four-velocity denoting flow velocities in space time is given by uμ=coshρcoshη,tanhρcosϕr,tanhρsinϕr,sinhη, while the kinetic freeze-out surface is given by d3σμ=coshη,0,0,−sinhητrdrdηdϕr. Here, η is the space-time rapidity. Assuming boost-invariant scenario where we have taken Bjorken correlation in rapidity, i.e., y=η [31] along the longitudinal or beam axis. Thus, equation (7) can be expressed as
(8)dNq/q¯eqpTdpTdy=12π·D∫0R0rdr∫0∞coshyexp−mTcoshycoshρTdy∫02πexppTsinhρcosϕTdϕ,where D=gqtfmT/2π2. Here, gq is the quark degeneracy factor, tf is the particle emission time, and mT=pT2+mq2 is the transverse mass.ρ in the integrand is a transverse rapidity variable which is given by ρ=tanh−1βr+ρabcos2ϕ, with ρa as a function of impact parameter, b, and gives the anisotropy dependence in the flow. βr=βsξn[32–35] is the radial flow, where βs is the maximum surface velocity and ξ=r/R0, with r as the radial distance from the center of the fireball. In the blast wave 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 [36]
(9)<βr>=∫βsξn+1dξ∫ξdξ=22+nβs.While the anisotropic parameter,ρa is written as
(10)ρab=1−ζ2−1−ζ1−ζ2+1−ζ,ζ=b2Ra.In our calculation, we use a linear velocity profile, (n=1), R0 is the maximum radius of the expanding source at freeze-out 0<ξ<1, and RA is the radius of colliding nucleus. b 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 p+p collision,
(11)dσpp→qq¯d2pTdy1dy2=2x1x2∑1,2fpx1,Q2·fpx2,Q2·dσ^12→qq¯dt^+1↔2×11+δ12.Here,x1 and x2 are momentum fractions carried by interacting partons from their respective colliding protons and are given by
(12)x1=2mTsexp−y1+exp−y2,x2=2mTsexp−y1−exp−y2.ApT cut of 2 GeV/c is taken for the jet production following other event generators like PYTHIA and HIJING [37, 38]. The parton density functions, fix,Q2, are taken to be CTEQ5M [39]. The partonic differential scattering cross-sections, dσ^/dt^, is calculated from the LO processes, gg⟶qq¯ and qq¯⟶qq¯. To incorporate NLO processes, we have taken a factor, ′K′, of value 2.5, and finally, nuclear overlap function, TAAb and EKS98 parametrization for shadowing effects are taken into account to convert particle production cross-section from p+p collision into A+A, particle spectra. Equation (11) is parametrized using a function with a power-like structure (Juttner distr.), 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.
(13)dNpp→qq¯d2pTdy1dy2=TAAb·dσpp→qq¯d2pTdy1dy2=TAAb·K·C·1+mTB−αfiq/q¯pT=1tf·π·RA2·mT·coshy−ηdNq/q¯id2pTdyTable 1
Extracted parameters of equation (11) at LO pQCD calculations of p+p collision at sNN=2.76TeV.
Quark flavour
α
B (GeV)
Cfm4
u
5.615
1.127
3.73376×103
u¯
5.999
1.099
8.73376×102
d
5.579
1.434
3.6286×103
d¯
5.953
1.401
9.1286×102
s=s¯
6.523
1.892
2.6317×102
c=c¯
7.250
3.287
2.32815Here too, we have assumed Bjorken correlation in rapidity, i.e.,y=η. Using the Glauber model, TAAb is calculated to be 260.50fm−2 for 0-5% centrality and 13.1fm−2 for 50-60% centrality of the colliding nuclei at sNN=2.76TeV. Using equations (8) and (13), the final distribution can be expressed as in equation (5). This gives the final pT distribution for quarks. The quark masses for the initial distributions are taken to be mu=2.3MeV, md=4.5MeV, ms=95MeV, and mc=1.25GeV. After the transport quark coalescence, formalism has been used to combine the quarks into hadrons. This will be discussed next.


## 3. 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 [40, 41]. In Refs. [42, 43], the authors have used a two-component behaviour of hadronic spectra. For low-pT (pT<5GeV/c), they have used the recombination mechanism for thermalized partons, while for pT>5GeV/c, a power law-like distribution with fragmentation has been used. In our work, we have concentrated our transport approach to low and intermediate pT (<5 GeV/c), where instead of adopting thermalized distribution directly, we have allowed jet distribution of partons to relax or thermalize and then proceed to hadronization. The idea was to use the BTE-RTA equation to study the applicability of interpolation of jet distribution with the blast wave equation at the intermediate pT region. The coalescence or recombination of partons into hadrons has been 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. And thus, the ReCo mechanism has been used for final transported quark distributions in the present work. The ReCo mechanism also highlights the major contribution of partonic degrees of freedom in the observed hadron flow. The process such as gg⟶gg has been neglected as gluon contributions are mostly at low pT<1GeV/c, while high pT gluons contribute to hadrons via a fragmentation mechanism. In the intermediate-momentum region, constituent quark counting becomes important for recombination process. At the hadronization surface, only constituent quarks behave as effective degrees of freedom with mass. However, it must be noted that for net entropy and energy density calculation, gluon contribution is most vital [42–46]. In the present work, we have not included gluon contribution to the hadron production which is one of the main differences from the earlier works and it is most visible at the low momenta where gluon contribution is important. Unlike earlier works, we have not included the fragmentation mechanism for high pT particles and only focused our observations to the intermediate pT. Another difference in the current work is the absence of the hadron decay mechanism which is important at a low-momentum region. However, these are out of the scope of the current work and could constitute a better prospective for future research.The coalescence model can be applied to the quarks at the hadronization surface when two (three) quarks recombine to form mesons (baryons) [47–49]. The model can be further utilized in describing observed spectra of light nuclei such as deuteron which contains a neutron and a proton [50].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(14)NM=gM∫mT1coshy1−η1d3r1×mT2coshy2−η2d3r2dp→T1dy1dp→T2dy2×fqr1;p1fq¯r2;p2WMr1,r2;p1,p2,where r→1,r→2 and p→T1,p→T2 are the spatial and transverse momentum coordinates of the combining quarks and antiquarks and fq/q¯ are the quark distribution functions. WMr1,r2;p1,p2 is the Wigner function convoluting two partonic distributions. gM in the front of equation (14) is meson degeneracy factor. We have also assumed Bjorken correlation in rapidities, y1=η1 and y2=η2, throughout. We also assumed ∣y1∣=∣y2∣≤0.5. This ensures a close phase space for quarks in both momentum and configuration spaces.We have assumed the delta functions correlation,δ3p→−p→1−p→2, and δ32→R−r→1−r→2. We have defined the partons in the spatial and momentum coordinates in the C.M. frame of meson, such as
(15)R→=r→1+r→22,r→=r→1−r→2,p→=p→1+p→2,q→=p→2−p→12,so that we can derive
(16)fqr1,p1⟶fqR→+r→2,p→2+q→,fq¯r2,p2⟶fq¯R→−r→2,p→2+q→,WMr→1−r→2;p→1−p→2⟶WMr,q.Here, we have also assumed that∣r→∣ is small compared to ∣R→∣ and thus neglected ∣r→∣ in the quark distributions, fq/q¯. Thus, we have
(17)NM=gM∫d3rd3R2π6∫d2qd2pT2π6mT1·mT2·fqR→,p→T2−q→fq¯R→,p→T2+q→WMr,q.We have now,(18)dNMd2pT=gM∫d3R2π3∫d2qd3r2π6mT1.mT2.fqfq¯R→,p→T2+q→WMr,q,where meson transverse mass factor is given by, MT=pT2+M2.As for the Wigner function,WM, we can use the following relation,
(19)WMq=∫d3rWMr,q.Therefore, equation (18) is transformed as
(20)dNMd2pT=gM∫d3R2π3∫d2q2π3mT1·mT2·fqR→,p→T2−q→fq¯R→,p→T2+q→WMq.To simplify our equations, we convert our momentum variable into light-cone coordinates,kμ of the interacting quarks in the momentum space of the hadron as follows:
(21)q→=p→2−k→,so that, d3q=d3k and d3k=dk+d2k⊥,
(22)k±=k0±k32,k⊥2=2k+k−−k2,andk+=x·p+.We also assume that the partons recombining into hadrons have their momenta almost parallel to the final hadron. Sok⊥ can be considered to be very small compared to k+ and its dependency in the quark distribution, fq/q¯, has also been neglected. It can be shown following equation (22) that the parton momentum, k≈x.p, where x is the momentum fraction of the final hadron’s momentum, carried by its constituent quarks during recombination [51, 52]. Putting the above conditions into the equation, and assuming the normalization,
(23)∫dxd2k⊥p+2π3WMx,k⊥2=1.Finally, we can write(24)dNMd2pT=gM∫d3R2π3∫01dxfqR→,xpTfq¯R→,1−xpTWMx.We haved4R=pμ.dσμ along the unit normal direction, uR=1,0,0,0 at the freeze-out hypersurface. Similarly, for the baryons, one can derive to show
(25)dNBd2pT=gB∫d3R2π3∫01dx1∫01dx2fqR→,x1pTfqR→,x2pTfqR→,1−x1−x2pTWBx1,x2.To illustrate on our calculations, one may use equation (8) as an example and show that
(26)fqeqR,xp·fq¯eq1−xp=e−p1+p2·uR/T,p1+p2·uR=μTMx,pTcoshρ−pTsinhρcosϕr−ϕp.Thus, one may calculate to show(27)μTMx,pT=mq2+x2pT2+mq¯2+1−x2pT2.Similarly for the baryons, one may write(28)μTBx1,x2,pT=mq2+x12pT2+mq2+x22pT2+mq2+1−x1−x22pT2.We have replaced transverse mass,mT, by expressions from equations (26) and (28), throughout our calculations. We have also assumed a general Gaussian distribution as Wigner functions, WM (for mesons), WB (for baryons), which are given by
(29)WMx=e−x−0.52/2σM2,WBx=e−x1−x22+x1+x2−0.662/2σB2.Here,2σ2 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 a 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.Thus, using equation (4), in equations (24) and (25), we calculate v2 of the final hadrons at midrapidity as
(30)v2pTy=0=∫dNM/B/d2pT×cos2ϕdϕ∫dNM/B/d2pTdϕ.


## 4. Results and Discussions

We would like to reiterate that in the current work, using BTE in RTA, we have transported quarks of various flavors (u,d,s, and c) produced promptly (assumed to be out of equilibrium) from initial gluon fusion to thermalization and freeze-out time and recombined them into hadrons using the coalescence mechanism. We have neglected the decay contributions to final hadron spectra as well as effects due to hadronic interactions. We have tried to extract the correlation between parameters such as the radial flow, βs, and the ratio of freeze-out time and relaxation/thermalization time, tf/τ, We have presented the results on the elliptic flow (v2) of various identified hadrons like pions, kaons, protons, D meson, and lambda for different centralities of Pb+Pb collision at sNN=2.76TeV. While analysing the data, we have kept the freeze-out temperature (Tf) for the hadrons at 0.095 GeV for the most central collisions (0-5)% and 0.11 GeV for most peripheral collisions (50-60)% [53]. We assume that the value of Tf is smaller for the central collisions in comparison to the peripheral collisions. The above assumption on the freeze-out temperature is based on the fact that the freeze-out in peripheral collisions occurs quicker than in the most central collisions [53]. The dependence of identified hadrons’ v2 is studied by varying two parameters, βs and tf/τ, using equation (29). Based on the closest explanation of the data, we have kept the width of the Wigner function 2σ2 fixed at 0.0009 for mesons and 0.04 for baryons in our calculations. As discussed earlier that the constituent quarks for recombination process occupy a close phase space, we needed a narrow Gaussian function and not a delta function so as to avoid the collinear divergences as well as to satisfy the above condition. We have not included the flow from hadronic medium as they are most visible for particles below pT<2GeV [54, 55] where our results focus on pT≥2GeV.In Figure1, we have shown pT spectra of various charged hadrons in the most central collisions of Pb+Pb at sNN=2.76TeV. The coalescence method is employed to form hadrons from quarks at the freeze-out surface. The resulting transverse momentum distributions are then drawn and compared with the experimental data from ALICE@CERN [57–59]. It is found that the discussed model in the above section explains the experimental data in the moderate pT region.Figure 1
The transverse momentum spectra of (π++π−), (K++K−), (p+p¯), and D0 mesons versus pT for most central Pb+Pb collisions at sNN=2.76TeV. Symbols are experimental data points [56, 57], and lines are the model results from equations (24) and (25).In Figure2, we have shown v2 of (π++π−). The left plot shows the variation of v2 with pT for different surface velocity parameters, βs, while the right plot shows for different tf/τ. Three different values of βs keeping tf/τ fixed are taken and vice versa. Generally speaking, our theoretical results match with the experimental data within errors for the discussed particles, from the mid-pT region to the max pT shown. However, the model fails to explain the data for pT<1.0GeV/c. The reason may be due to the absence of pions from decays of resonances [60]. Pions also stand out as an example that shows coalescence picture should work mostly in the mid-pT region.Figure 2
The elliptic flow (v2) of (π++π−) versus pT at constant tf/τ and βs for peripheral collisions (50-60)% at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.In Figure3, the elliptic flow of pions (π++π−) is presented as a function of pT for various centralities at sNN=2.76TeV for Pb+Pb collisions. Symbols are the experimental data, and lines are model results. Here, freeze-out temperatures (Tf) are taken smaller for most central collision than to peripheral collisions. The model results are found to explain the data qualitatively beyond pT=1GeV/c for all the centralities within the error bar. However, the quark coalescence mechanism is not able to explain the data below pT=1GeV/c. Experimentally, it is found that v2 for (50-60)% appears to be inverse in order compared to (40-50)% due to statistical fluctuations. However, the model follows the expected trend of higher v2 for higher centralities.Figure 3
The elliptic flow (v2) of (π++π−) versus pT for various centralities for Pb+Pb collisions at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.In the left panel of Figure4, we have shown the elliptic flow or azimuthal anisotropy v2 and spatial anisotropy ε2 of the pions versus Npart ε2 is generally defined in terms of spatial coordinates x,y of participants’ nucleons in the transverse plane. It can be written as
(31)ε2=x2−y2x2+y2.Figure 4
(a) The elliptic flow (v2) and spatial anisotropy (ε2) of (π++π−) versus Npart for Pb+Pb collisions at sNN=2.76TeV. (b) The ratio of v2 and ε2 of pions versus Npart.
(a)
(b)In this paper, Glauber-MC formalism [61] has been employed to calculate ε2. Both v2 and ε2 decrease with Npart, which is expected. In the right panel of Figure 4, we show the ratio of v2 and ε2 vs. Npart 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 Figure5, we have presented the variations of v2 of Kaons, (K++K−) with pT for 50-60% centrality. The left panel is v2 versus pT for different βs at constant tf/τ, while the right panel shows v2 versus pT for different tf/τ at constant βs. Three different values of βs keeping tf/τ fixed are taken and vice versa. The theoretical curves tend to overestimate the data although it gives a consistent explanation as to the nature of the shape of Kaons v2 shown by the experimental data. Also, the plot on the left side shows that the theoretical lines cross each other for the different values of βs, which shows greater sensitivity of v2 on the surface velocity of the fireball. The theoretical line is quite close to experimental points at low pT which shows that large mass should have less contribution from resonance decays.Figure 5
The elliptic flow (v2) of (K++K−) versus pT at constant tf/τ and βs for peripheral collisions (50-60)% at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.In Figure6, we have shown v2 of K-short (KS0). The left plot shows the variation of v2 with pT taking various values of βs. The right plot of the figure represents the variation of v2 with pT taking different values of tf/τ. Three different values of βs keeping tf/τ fixed are taken and vice versa. KS0 is a little heavier than Kaons, which is why the tf/τ and βs values are almost similar in both the cases. Similarly, the theoretical curve tends to overestimate the data up to pT=3GeV/c. However, the theoretical curve shows a gradually increasing trend and slopes down smoothly at high pT.Figure 6
The elliptic flow (v2) of KS0 versus pT at constant tf/τ and βs for peripheral collisions (50-60)% at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.Figure7 represents the variations of v2 with respect to pT of phi, ϕ. The left plot shows the variation of v2 for different βs at constant tf/τ, while the right plot shows the variation of v2 with parameter tf/τ keeping βs constant. Three different values of βs keeping tf/τ fixed are taken and vice versa. Phi meson’s results show a gradual rise in the values of v2 with an increase in pT as shown in the plot. Although, the data points show a very small variation after pT>3GeV/c, the theoretical curves drop smoothly and continues to do so at pT=6.0GeV/c.Figure 7
The elliptic flow (v2) of phi-meson (ϕ) versus pT at constant tf/τ and βs for peripheral collisions (50-60)% at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.In Figure8, the elliptic flow of D meson is presented as a function of pT for centrality 30-50% at sNN=2.76TeV Pb+Pb collisions. The left panel shows elliptic flow for various βs at constant tf/τ. The model shows a rise in v2 for pT<3GeV/c and falls smoothly afterwards. The data points show almost a constant v2 value and also, a number of data points are small to be explained satisfactorily by our model. The right plot is v2 of D meson for different values of tf/τ keeping βs constant.Figure 8
The elliptic flow (v2) of D meson versus v2 at constant tf/τ and βs for peripheral collisions (30-50)% at sNN=2.76TeV. Symbols are experimental data points [59], and lines are the model results.In Figure9, we have shown the variations of the elliptic flow of p+p¯ with respect to pT for 50-60% centrality of Pb+Pb collisions at sNN=2.76TeV. In the left hand side of the figure, we show the v2 for different values of βs keeping tf/τ fixed. It is found that the model results explain the experimental data qualitatively above pT=1GeV/c for βs=0.9. The right hand side of the figure is the results for various tf/τat constant βs. Again, there is a good agreement between the model calculations and experimental data above pT=1GeV/c for tf/τ=2.2.Figure 9
The elliptic flow (v2) of (p+p¯) versus pT at constant tf/τ and βs for peripheral collisions (50-60)% at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.In Figure10, the elliptic flow of Λ+Λ¯ is presented with respect to pT for centrality 50-60% at Pb+Pb collisions. In the left hand side of the figure, we show the v2 for different values of βs keeping tf/τ fixed. It is found that the model results explain the experimental data qualitatively above pT=1GeV/c for βs=0.89. The right hand side of the figure is the results for various tf/τ at constant βs. Again, there is a good agreement between the model calculations and experimental data above pT=1GeV/c for tf/τ=2.4.Figure 10
The elliptic flow (v2) of (Λ+Λ¯) versus pT at constant tf/τ and βs for peripheral collisions (50-60)% at sNN=2.76TeV. Symbols are experimental data points [58], and lines are the model results.In Figure11, we have plotted the v2 of Λ hadron with its three constituent quarks, u,d, and s. Although the flow of the constituent quarks starts long before pT<1.0GeV/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 βs and tf/τ taken from Λ v2 plot are kept fixed for its constituent quarks, u,d,ands.Figure 11
Comparison ofv2 of constituent quarks with the final Λ hadron v2 for peripheral Pb+Pb collisions at sNN=2.76TeV.In Figure12, the correlation of tf/τ with βs is shown for various identified hadrons observed after extracting the values from the model results on the elliptic flow with the experimental data at sNN=2.76TeV for peripheral Pb+Pb collisions. In this plot, we find that with the increase in tf/τ, the surface velocity βs 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 the lightest π-meson towards heavier D0 meson. A similar trend is also being observed for baryons, p and Λ.Figure 12
The plot ofβs versus tf/τ of various particles for peripheral Pb+Pb collisions at sNN=2.76TeV.


## 5. Summary and Conclusions

We have used the quark coalescence method for hadronization and the Boltzmann transport equation in relaxation time approximation to estimate the elliptic flow,v2, for the identified hadrons in Pb+Pb collisions at sNN=2.76TeV. The important findings are summarised as follows:
(1)
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 lowpT
(2)
The present formalism successfully attempts to connect the particle production from prompt interaction of initially produced partons with finally produced hadrons at hadronization hypersurface. For intermediatepT ranges, the present formalism may successfully interpolate nonequilibrium or jet-like quarks into blast wave distribution. The hadronic medium effects have not been taken into account. Similarly, resonance decays into observed particles particularly in the pion sector have been neglected as well
(3)
We have found a correlation between the radial part of the transverse flow andtf/τ while explaining the v2 spectra in peripheral collisions
(4)
We have also compared elliptic flow of constituent quarks,u,d,s with the final hadron, Λ. We find that v2 of each quark is around 1/3 of the final Λ baryon. This actually verifies the coalescence mechanism used in the present calculations
(5)
Higher mass quarks are found to have a lowerv2 as compared to lighter quarks. On the other hand, the flow of mesons behaves almost similarly in the mid-pT region although their flow parameter, βs, and time ratio, tf/τ, 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 v2 of final hadrons. Similarly, the hadronization mechanism as a part of the freeze-out dynamics also plays a major role in this regardWe 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 BTE-RTA mechanism.



---

*Source: 4728649-2020-03-23.xml*