Hydrodynamic models of particle production

Hydrodynamic models of particle production


Viscous hydrodynamics gives a satisfactory description of the transverse momentum spectra, of the elliptic and triangular flow, and of the femtoscopic correlations for particles produced in relativistic heavy-ion collisions. On general grounds, a similar collective behavior has been predicted for proton-lead (p-Pb) collisions at the LHC. We present results of the hydrodynamics calculation of the elliptic and triangular flow in p-Pb. We discuss the mass dependence of flow coefficients and of the average transverse momentum for identified particles.

1 Hydrodynamic model

The observation of the jet quenching and of the strong azimuthally asymmetric flow shows that in a high energy heavy-ion collision at ultrarelativistic energies a dense and hot droplet of matter is formed. Event-by-event hydrodynamic calculations [1, 2, 3, 4, 5, 6] are used to provide the response of the collective flow to the deformations of the fluctuating initial state. The eccentricity and triangularity of the initial density result in the formation of the elliptic and triangular flow in the spectra [3, 7]. Not only the average, but the whole distribution of the harmonic flow coefficients [8] can be reproduced [9] for the IP-Glasma initial conditions. The event-by event distribution of the hydrodynamic response shows that and follow the initial distribution of and , with nonlinearities in the response appearing for the harmonics and higher. Non-flow correlations between particles with soft momenta can be included in the simulation giving correlations between unlike-sign particles mainly at small relative angle and relative pseudorapidity [10]. The odd and even components of the directed flow as a function of pseudorapidity, appear from the expansion of a tilted, fluctuating fireball [11, 12, 13, 14].

The harmonic event-plane for each order can be determined from the initial density [15, 16], while the hydrodynamic evolution that involves nonlinear coupling between harmonics yields additional correlations between the event planes, in good agreement with experiment [16]. Fluctuations may lead to a decorrelation of the event planes at different rapidities [17, 18, 19], a prediction that could be tested experimentally. Fluctuations of the volume of the initial fireball lead to event-by-event fluctuations of the average transverses momentum [20], in an analogous way as for the response of the flow harmonics to the fireball anisotropies. The hydrodynamic model gives predictions for the spectra of identified particles that are in satisfactory agreement with experiment at LHC energies [21, 22, 23, 24, 25]. Advanced hydrodynamic calculations are used to extract the properties of the quark-gluon plasma at different temperatures [26, 27, 28, 9, 29] and barionic densities [30, 31]. The hard equation of state of the quark-gluon plasma leads to a rapid build up of the transverse flow, which implies a strong dependence of the femtoscopic radii on the pion pair momentum [32, 33].

The p-Pb collisions have been proposed as a laboratory to study the initial state effects and to obtain reference data for the Pb-Pb collisions [34]. On the other hand, it has been predicted that collective expansion of the fireball formed in central p-Pb collisions is significant and can lead to observable elliptic and triangular flow [35]. Experiments with p-Pb collisions at the LHC measure two-particle correlation function in relative pseudorapidity and relative azimuthal angle [36, 37, 38]. The two-dimensional correlation function presents two ridge-like structures, elongated in pseudorapidity, with pairs collimated in the same direction (same-side ridge, ) and in the opposite direction (away-side ridge, ). The azimuthal structure in the two-dimensional correlation function can be understood as an effect of the azimuthally asymmetric collective flow and of the transverse momentum conservation [39], in a similar way as in heavy-ion collisions [40, 41]. The formation of the two ridges is predicted in the color glass condensate approach as well [42, 43, 44]. In Ref. [35] it has been proposed to look for collectivity in small systems using d-A collisions. In central d-A collisions the density profile in the transverse plane is determined by the deuteron wave-function, hence the eccentricity is large , which gives a large elliptic flow [35]. The analysis of correlation functions in d-Au collisions at RHIC shows a strong component [45], but quantitative conclusions are difficult due to large non-flow contributions in the low multiplicity environment at the RHIC energies.

2 Transverse flow in p-Pb interactions

The multiplicity in central p-Pb collisions is comparable to peripheral Pb-Pb collisions. The source size can be estimated using the Glauber Monte Carlo model [46, 47]. The value of the root mean square radius, 1-2 fm, depends on the details of the energy deposition in the model, and grows for more central collisions.

Figure 1: Left: The freeze-out isotherms in the plane (solid) and the plane (dashed) in a sample event. Right: Mean transverse momentum of identified particles as a function of charged particle density from hydrodynamics (panel (a)) and HIJING  (panel (b)), ALICE Collaboration data [48].

The expansion of the large energy density deposited in the small volume generates large transverse collective flow [46, 47, 49, 50]. In the left panel of Fig. 1 is shown the freeze-out hypersurface for an event with 19 participant nucleons. The transverse size increases during the freeze-out and can be measured using interferometry radii. The values of the interferometry radii grow as a power of the multiplicity [51]. The slope of this dependence is different for A-A and for p-p collisions. It may suggest that the dynamics of a typical p-p interaction is different than in peripheral A-A collisions. The hydrodynamical model predicts for the p-Pb system interferometry radii that are close to those corresponding to A-A interactions [49].

The average transverse momentum of particles produced in the p-p, p-Pb and A-A interactions increases with the event multiplicity [48, 52]. For p-p interactions this increase can be explained as a color reconnection effect [53]. Modeling the particle production in p-Pb and A-A collisions as a superposition of nucleon-nucleon interactions one finds the average transverse momentum below the measured value [54, 48]. This effect is visible in the right panel of Fig. 1, where the HIJING model, based on the superposition of nucleon-nucleon collisions, underpredicts the average transverse momentum. This difference leaves room for additional collective transverse velocity to be generated during the expansion phase in p-Pb collisions. An important characteristic of the collective transverse flow is the mass hierarchy. The average transverse momentum from collective flow is larger for heavier particles. The data for , , and p transverse momenta can be reproduced naturally in the hydrodynamic framework [55] (Fig. 1 right panel). The consistency of the hydrodynamic calculations with the experimental data validates the collective flow interpretation, while we note that the mass hierarchy of the average transverse momentum can also be understood as coming from geometrical scaling [56]. The transverse momentum spectra in p-Pb collisions contain a soft component coming from the collective expansion with statistical emission at freeze-out [35, 50, 57]. The spectra can be reproduced in the region of intermediate in the EPOS LHC model [50, 57, 52].

3 Elliptic and triangular flow in p-Pb interactions

The most important evidence for the collective expansion in the A-A collisions is the observation of elliptic and triangular flow. The elliptic flow coefficient has been measured using 2- and 4-particle cumulants in p-Pb collisions [58, 59]. The data is consistent with the predictions of the hydrodynamic model [35, 46, 60, 61] (Fig. 2 right panel).

Figure 2: Left: and for charged particles from the hydrodynamic calculation, CMS Collaboration data  [59]. Right: for pions, kaons and protons from the hydrodynamic model, ALICE Collaboration data [62]).

In p-Pb interaction the model uncertainty on the value of the initial eccentricity is a more severe problem than for the A-A system [46, 54]. The small size and short lifetime of the system make the final results depend significantly on other parameters of the model as well: the initial time of the expansion, the freeze-out temperature, the shear viscosity. An additional difficulty in comparing to the experimental data is related to the large contribution of non-flow correlations that are not implemented in the models. The relativistic viscous hydrodynamics can be reliably applied if the velocity gradients in the system are moderate and if the mean free path is smaller than the size of the system. In small systems the contribution of non-thermal corona is relatively more important, which can reduce the value of the elliptic flow. These assumptions become more and more questionable when applying the model to more peripheral p-Pb collisions. All of the above mentioned factors mean that the agreement of the calculation with the data can be at best semi-quantitative.

An important experimental result is the observation of mass splitting of the elliptic flow in p-Pb interactions [62]. The smaller value of the elliptic flow coefficient for heavy particles is a characteristic of the collective elliptic flow. This feature is reproduced by the calculation [55, 61] (Fig. 2 right panel). The observation of triangular flow in p-Pb interactions is a strong argument for the collective expansion scenario. The predictions of the hydrodynamic model, with fluctuating initial conditions are consistent with measurements [35, 46, 60] (Fig. 2 right panel). This means that the model captures correctly the fluctuations of the initial state and describes realistically its collective expansion. The mass splitting for the triangular flow in p-Pb collisions is smaller than for , moreover, it is distorted by resonance decays.

4 Collectivity in p-Pb

The formation of a fireball of strongly interacting fluid in ultrarelativistic heavy-ion collisions is well established. The observation of elliptic and triangular flow in agreement with predictions of nearly prefect fluid hydrodynamics, the quantitative description of transverse momentum spectra, and jet quenching are strong evidence of the creation of the quark-gluon plasma. Quantitatively the hydrodynamic model can be reliably applied in central A-A collisions. Many experimental efforts have been devoted to the study the properties of the matter in the fireball when changing the system size and the energy. The data indicate that a collectively expanding system is formed in peripheral A-A collisions and at different energies of the RHIC beam energy scan.

It has been generally expected that in p-Pb collisions at the LHC, the final state interactions would be negligible [34], but the hydrodynamic model applied to initial conditions extrapolated to the p-Pb system at TeV energies predicted a significant collective expansion [35] that could be observed as the elliptic and triangular flow. Many experimental data from p-Pb interactions are consistent with the collective expansion scenario, the double-ridge [39], the elliptic and triangular flow [35, 46, 60, 61], the mass hierarchy of average transverse momentum [55, 50], and of the elliptic flow [55, 61]. The observed appearance of the flow in low multiplicity p-Pb events [59] and general arguments suggest that the p-Pb system can be used as a testing ground for the onset of collectivity in small systems.


PB and WB acknowledge the support of the Polish National Science Centre, grant DEC-2012/06/A/ST2/00390 and PL-Grid infrastructure. GT acknowledges the financial support received from the Helmholtz International Centre for FAIR within the framework of the LOEWE program (Landesoffensive zur Entwicklung Wissenschaftlich-Ökonomischer Exzellenz) launched by the State of Hesse, and from DOE under Grant No. DE-FG02-93ER40764.



  1. Schenke B, Jeon S and Gale C 2011 Phys. Rev. Lett. 106 042301
  2. Petersen H, Qin G Y, Bass S A and Muller B 2010 Phys.Rev. C82 041901
  3. Gardim F G, Grassi F, Luzum M and Ollitrault J Y 2012 Phys. Rev. C85 024908
  4. Bożek P and Broniowski W 2012 Phys. Rev. C85 044910
  5. Qiu Z, Shen C and Heinz U 2012 Phys. Lett. B707 151
  6. Pang L, Wang Q and Wang X N 2012 Phys.Rev. C86 024911
  7. Niemi H, Denicol G, Holopainen H and Huovinen P 2012 (Preprint 1212.1008)
  8. Aad G et al. (ATLAS Collaboration) 2013 (Preprint 1305.2942)
  9. Gale C, Jeon S, Schenke B, Tribedy P and Venugopalan R 2013 Phys.Rev.Lett. 110 012302
  10. Bożek P and Broniowski W 2012 Phys.Rev.Lett. 109 062301
  11. Bożek P and Wyskiel I 2010 Phys. Rev. C81 054902
  12. Csernai L, Magas V, Stocker H and Strottman D 2011 Phys.Rev. C84 024914
  13. Teaney D and Yan L 2011 Phys. Rev. C83 064904 1010.1876)
  14. Retinskaya E, Luzum M and Ollitrault J Y 2012 Phys.Rev.Lett. 108 252302
  15. Teaney D and Yan L 2012 Phys.Rev. C86 044908
  16. Heinz U, Qiu Z and Shen C 2013 Phys.Rev. C87 034913
  17. Bożek P, Broniowski W and Moreira J 2011 Phys. Rev. C83 034911
  18. Petersen H, Bhattacharya V, Bass S A and Greiner C 2011 Phys.Rev. C84 054908
  19. Xiao K, Liu F and Wang F 2013 Phys.Rev. C87 011901
  20. Broniowski W, Rybczynski M, Obara L and Chojnacki M 2010 Acta Phys.Polon.Supp. 3 513
  21. Shen C, Heinz U, Huovinen P and Song H 2011 Phys. Rev. C84 044903
  22. Bożek P and Wyskiel-Piekarska I 2012 Phys. Rev. C85 064915
  23. Karpenko I, Sinyukov Y and Werner K 2013 Phys.Rev. C87 024914
  24. Werner K, Karpenko I, Bleicher M, Pierog T and Porteboeuf-Houssais S 2012 Phys.Rev. C85 064907
  25. Abelev B et al. (ALICE Collaboration) 2013 (Preprint 1303.0737)
  26. Shen C, Bass S A, Hirano T, Huovinen P, Qiu Z et al. 2011 J.Phys. G38 124045
  27. Niemi H, Denicol G S, Huovinen P, Molnar E and Rischke D H 2011 Phys. Rev. Lett. 106 212302
  28. Bożek P 2010 Phys. Rev. C81 034909
  29. Luzum M and Ollitrault J Y 2013 Nucl.Phys.A904-905 2013 377c
  30. Steinheimer J, Koch V and Bleicher M 2012 Phys.Rev. C86 044903
  31. Karpenko I et al. 2013 these proceedings
  32. Broniowski W, Chojnacki M, Florkowski W and Kisiel A 2008 Phys. Rev. Lett. 101 022301
  33. Pratt S 2009 Phys. Rev. Lett. 102 232301
  34. Salgado C A 2011 J.Phys.G G38 124036
  35. Bożek P 2012 Phys. Rev. C85 014911
  36. Chatrchyan S et al. (CMS Collaboration) 2013 Phys. Lett. B718 795
  37. Abelev B et al. (ALICE Collaboration) 2013 Phys.Lett. B719 29
  38. Aad G et al. (ATLAS Collaboration) 2013 Phys.Rev.Lett. 110 182302
  39. Bożek P and Broniowski W 2013 Phys. Lett. B718 1557
  40. Takahashi J, Tavares B, Qian W, Andrade R, Grassi F et al. 2009 Phys. Rev. Lett. 103 242301
  41. Luzum M 2011 Phys.Lett. B696 499
  42. Dusling K and Venugopalan R 2013 Phys. Rev. D 87(5) 051502
  43. Dusling K and Venugopalan R 2013 Phys. Rev. D87 094034
  44. Dusling K and Venugopalan R 2013 Phys. Rev. D 87(5) 054014
  45. Adare A et al. (PHENIX Collaboration) 2013 (Preprint 1303.1794)
  46. Bożek P and Broniowski W 2013 Phys. Rev. C88 014903
  47. Bzdak A, Schenke B, Tribedy P and Venugopalan R 2013 Phys. Rev. C87 064906
  48. Abelev B B et al. (ALICE Collaboration) 2013 (Preprint 1307.1094)
  49. Bożek P and Broniowski W 2013 Phys. Lett. B720 250
  50. Pierog T, Karpenko I, Katzy J, Yatsenko E and Werner K 2013 (Preprint 1306.0121)
  51. Kisiel A (ALICE Collaboration) 2011 J.Phys. G38 124008
  52. Chatrchyan S et al. (CMS Collaboration) 2013 (Preprint 1307.3442)
  53. Ortiz A, Christiansen P, Cuautle E, Maldonado I and Paic G 2013 Phys.Rev.Lett. 111 042001
  54. Bzdak A and Skokov V 2013 (Preprint 1306.5442)
  55. Bożek P, Broniowski W and Torrieri G 2013 (Preprint 1307.5060)
  56. McLerran L, Praszalowicz M and Schenke B 2013 (Preprint 1306.2350)
  57. Abelev B B et al. (ALICE Collaboration) 2013 (Preprint 1307.6796)
  58. Aad G et al. (ATLAS Collaboration) 2013 (Preprint 1303.2084)
  59. Chatrchyan S et al. (CMS Collaboration) 2013 Phys.Lett. B724 213
  60. Qin G Y and Müller B 2013 (Preprint 1306.3439)
  61. Werner K, Bleicher M, Guiot B, Karpenko I and Pierog T 2013 (Preprint 1307.4379)
  62. Abelev B et al. (ALICE Collaboration) 2013 (Preprint 1307.3237)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description