1 / 53

Yu. Sinyukov, BITP, Kiev

Spectra and space-time scales in relativistic heavy ion collisions : the results based on HydroKinetic Model (HKM). Yu. Sinyukov, BITP, Kiev. Based on: Yu.S., I. Karpenko, A. Nazarenko J. Phys. G: Nucl. Part. Phys. 35 104071 (2008);

loyal
Download Presentation

Yu. Sinyukov, BITP, Kiev

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Spectra and space-time scales in relativistic heavy ion collisions: the results based on HydroKinetic Model (HKM) Yu. Sinyukov, BITP, Kiev Based on:Yu.S., I. Karpenko, A. Nazarenko J. Phys. G: Nucl. Part. Phys. 35 104071 (2008); S.V. Akkelin, Y. Hama, Iu. Karpenko, Yu.S., PR C 78, 034906 (2008) Yu.S. Two lectures in Acta Phys Polon. 40 (2009) GDRE - 2009

  2. Heavy Ion Experiments LHC FAIR E_lab/A (GeV)

  3. Thermodynamic QCD diagram of the matter states • The thermodynamic arias occupied by different forms of the matter Theoretical expectations vs the experimental estimates

  4. UrQMD Simulation of a U+U collision at 23 AGeV

  5. Soft Physics” measurements A x t ΔωK Landau/Cooper-Frye prescription A p=(p1+ p2)/2 q= p1- p2 Tch and μchsoon after hadronization (chemical f.o.) (QS) Correlation function Space-time structure of the matter evolution, e.g., Radial flow

  6. Collective transverse flows P T Initial spatial anisotropydifferent pressure gradients momentum anisotropy v2 NPQCD-2007

  7. Energy dependence of the interferometry radii Energy- and kt-dependence of the radii Rlong, Rside, and Rout forcentral Pb+Pb(Au+Au) collisions from AGS to RHICexperiments measurednear midrapidity. S. Kniege et al. (The NA49 Collaboration), J. Phys. G30, S1073 (2004).

  8. Expecting Stages of Evolution in Ultrarelativistic A+A collisions t Relatively small space-time scales (HBT puzzle) 10-15 fm/c Early thermal freeze-out: T_th Tch 150 MeV 7-8 fm/c Elliptic flows 1-3 fm/c Early thermalization at 0.5 fm/c 0.2?(LHC) or strings 8

  9. HBT PUZZLE & FLOWS • Possible increase of the interferometry volume with due to geometrical volume grows is mitigated by more intensive transverse flows at higher energies: ,  is inverse of temperature • Why does the intensity of flow grow? More more initial energy density  more (max) pressure pmax BUT the initial acceleration is ≈ the same HBT puzzle Intensity of collective flows grow Time of system expansion grows Initial flows (< 1-2 fm/c) develop

  10. Interferometry radii Borysova, Yu.S., Akkelin,Erazmus, Karpenko: PRC 73, 024903 (2006)

  11. Duration of particle emission is taken into account by means ofenclosedfreeze-out hypersurface: volume emission surface emission vi=0.35 fm/c NPQCD-2009

  12. Ro/Rs ratio and initial flows NPQCD-2009

  13. Basic ideas for the early stage: developing of pre-thermal flows Yu.S. Acta Phys.Polon. B37 (2006) 3343; Gyulassy, Yu.S., Karpenko, Nazarenko Braz.J.Phys. 37 (2007) 1031. For nonrelativistic gas For thermal and non-thermal expansion at : Hydrodynamic expansion: gradient pressure acts So, even if : and Free streaming: Gradient of density leads to non-zero collective velocities In the case of thermalization at later stage it leads to spectra anisotropy NPQCD -2009

  14. Comparision of flows at free streaming and hydro evolution NPQCD-2009

  15. Boost-invariant distribution function at initial hypersurface McLerran-Venugopalan CGC effective FT (for transversally homogeneous system) is the variance of a Gaussian weight over the color charges of partons A.Krasnitz, R.Venugopalan PRL84 (2000) 4309; A. Krasnitz, Y. Nara, R. Venugopalan: Nucl. Phys.A 717 (2003) 268, A727 (2003) 427;T. Lappi: PRC 67 (2003) 054903, QM 2008 (J.Phys. G, 2008) Transversally inhomogeneous system: <transverse profile> of the gluon distribution proportional to the ellipsoidal Gaussian defined from the best fit to the density of number of participants in the collisions with the impact parameter b. If one uses the prescription of smearing of the -function as , then . As the result the initial local boost-invariant phase-space density takes the form 15

  16. Developing of collective velocities in partonic matter at pre-thermalstage (Yu.S. Acta Phys. Polon. B37, 2006) • Equation for partonic free streaming in hyperbolic coordinates between • Solution where 16

  17. Collective velocity developed at pre-thermal stage from proper time tau_0 =0.3 fm/c by supposed thermalization time tau_i = 1 fm/c for scenarios of partonic freestreaming and free expansion of classical field. The results are compared with thehydrodynamic evolution of perfect fluid with hard equation of state p = 1/3 epsilon startedat tau_0 .Impact parameter b=0. Yu.S., Nazarenko, Karpenko: Acta Phys.Polon. B40 1109 (2009)

  18. Collective velocity developed at pre-thermal stage from proper time tau_0 =0.3 fm/c by supposed thermalization time tau_i = 1 fm/c for scenarios of partonic freestreaming. The results are compared with thehydrodynamic evolution of perfect fluid with hard equation of state p = 1/3 epsilon startedat tau_0 .Impact parameter b=6.3 fm. Yu.S., Nazarenko, Karpenko: Acta Phys.Polon. B40 1109 (2009)

  19. Initial parameters becomes to be isotropic at =1 fm/c: The "effective" initial distribution - the one which being used in the capacity of initial condition bring the average hydrodynamical results for fluctuating initial conditions: is a fitting parameter where For central Au+Au (Pb+Pb) collisions

  20. Hydro-Evolution: Equation of State EoS from LattQCD (in form proposed by Laine & Schroder, Phys. Rev. D73, 2006) MeV The EoS accounts for gradual decays of the resonances into expanding hadronic gasconsistiong of 359 particle species with masses below 2.6 GeV.The EoS in this non chemically equilibrated system depends now on particle number densities of all the 359 particle species. Since the energy densities in expanding system do not directly correlate with resonance decays, all the variables in the EoS depends on space-time points and so an evaluation of the EoS is incorporated in the hydrody-namic code. We calculate the EoS below in the approximation of ideal multi-component hadronic gas. MeV

  21. Hydro-evolution Chemically equilibrated evolution of the quark-gluon and hadron phases. For : MeV The system evolves as non chemically equilibra- ted hadronic gas. The conception of the chemical freeze-out imply that mainly only elastic collisi- ons and the resonance decays among non-elastic reactions takes place because of relatively small densities allied with a fast rate of expansion at the last stage. Thus, in addition to (1) the equa-tions accounting for the particle number conser- vation and resonance decays are added: MeV where denote the average number of i-th particles coming from arbitrary decay of j-th resonance, is branching ratio, is a number of i-th particles produced in decaychannel.

  22. Cooper-Frye prescription (CFp) t t z r • CFp gets serious problems: • Freeze-out hypersurface contains non-space-like • sectors • artificial discontinuities appears across • Sinyukov (1989), Bugaev (1996), Andrelik et al (1999); • cascade models show that particles escape from the system about whole time of its evolution. Hybrid models (hydro+cascade) and the hydro method of continuous emission starts to develop.

  23. Hybrid models: HYDRO + UrQMD (Bass, Dumitru (2000)) t t t UrQMD HYDRO z r The initial conditions for hadronic cascade models should be based on non-local equilibrium distributions • The problems: • the system just after hadronization is not so dilute to apply hadronic cascade models; • hadronization hypersurface contains non-space-like sectors (causality problem: Bugaev, PRL 90, 252301, 2003); • hadronization happens in fairly wide 4D-region, not just at hypersurface , especially in crossover scenario. 20 May 2009 23

  24. Basic ideas for the late stage Yu.S., Akkelin, Hama: PRL. 89, 052301 (2002); + Karpenko: PRC 78 034906 (2008). Hydro-kinetic approach Continuous emission t • is based on combination of Boltsmann equation and hydro for finite expanding system; • provides evaluation of escape probabili- ties and deviations (even strong) of distri-bution functions from local equilibrium; • accounts for conservation laws at the particle emission; PROVIDE earlier (as compare to CF-prescription) emission of hadrons, because escape probability accounts for whole particle trajectory in rapidly expanding surrounding (no mean-free pass criterion for freeze-out) x F. Grassi,Y. Hama, T. Kodama

  25. Hydro-kinetic approach • MODEL • provides evaluation of escape probabilities and deviations (even strong) • of distribution functions [DF] from local equilibrium; • is based on relaxation time approximation for relativistic finite expanding system; • (Shear viscosity ). • accounts for conservation laws at the particle emission; • Complete algorithm includes: • solution of equations of ideal hydro [THANKS to T. Hirano for possibility to use • code in 2006] ; • calculation of non-equilibrium DF and emission function in first approximation; • [Corresponding hydro-kinetic code: Tytarenko, Karpenko,Yu.S.(to be publ.)] • Solution of equations for ideal hydro with non-zero left-hand-side that accounts • for conservation laws for non-equlibrated process of the system which radiated • free particles during expansion; • Calculation of “exact” DF and emission function; • Evaluation of spectra and correlations. Is related to local *

  26. Boltzmann equations and Escape probabilities Boltzmann eqs (differential form) and are G(ain), L(oss) terms for p. species Escape probability (for each component ) 26

  27. Spectra and Emission function Boltzmann eqs (integral form) Index is omitted everywhere Spectrum Method of solution (Yu.S. et al, PRL, 2002) 27

  28. Saddle point approximation Spectrum Emission density where Normalization condition Eqs for saddle point : Physical conditions at

  29. Cooper-Frye prescription Spectrum in new variables Emission density in saddle point representation Temporal width of emission Generalized Cooper-Frye f-la

  30. Generalized Cooper-Frye prescription: t 0 Escape probability r Yu.S. (1987)-particle flow conservation; K.A. Bugaev (1996) (current form) 30

  31. Momentum dependence of freeze-out Here and further for Pb+Pb collisions we use: initial energy density Pt-integrated EoS from Lattice QCD when T< 160 MeV, and EoS of chemically frozen hadron gas with 359 particle species at T< 160 MeV. RANP08

  32. The pion emission function for different pT in hydro-kinetic model (HKM)The isotherms of 80 MeV is superimposed.

  33. The pion emission function for different pT in hydro-kinetic model (HKM). The isotherms of135 MeV (bottom) is superimposed.

  34. Transverse momentum spectrum of pi− in HKM, compared with the suddenfreeze-out ones at temperatures of 80 and 160 MeV with arbitrary normalizations.

  35. Conditions for the utilization of the generalized Cooper-Frye prescription • For each momentum p, there is a region of r where the emission function has a • sharp maximum with temporal width . ii) The width of the maximum, which is just the relaxationtime ( inverse of collision rate), should be smaller than the corresponding temporal homogeneitylength of the distribution function: 1% accuracy!!! iii) The contribution to the spectra from the residual region of r where the saddlepoint methodis violated does not affect essentially theparticle momentum spectrum. iiii) The escape probabilities for particles to be liberated just from the initial hyper-surface t0 are small almost in the whole spacial region (except peripheral points) Then the momentum spectra can be presented in Cooper-Frye form despite it is, in fact, not sadden freeze-out and the decaying region has a finite temporal width . Also, what is very important, such a generalized Cooper-Frye representation is related to freeze-out hypersurface that depends on momentum p and does not necessarilyencloses the initially dense matter.

  36. Inititial conditions in Pb+Pb and Au+Au colliisions for SPS, RHIC and LHC energies • We use 0.2 (the average transverse flow = 0.17) for SPS energies and 0.35 ( = 0.26) for all the RHIC and LHC energies. Note that this parameter include also a correction of underestimated transverse flow since we did not account for the viscosity effects. • For the top SPS energy GeV/fm3 = 4.7 GeV/fm3), for the top RHIC energy GeV/fm3 = 8.0 GeV/fm3 ) We also demonstrate results at = 40 GeV/fm3 and = 50 GeV/fm3 that probably can correspond to the LHC energies 4 and 5.6 A TeV. As for latter, according to Lappy, at the top LHC energy recalculated to the time 1 fm/c is 0.07*700=49 GeV/fm3 . The chosen values of approximately correspond to the logarithmic dependence of the energy density on the collision energy.

  37. Energy dependence of pion spectra

  38. Energy dependence of the interferometry scales, R_long

  39. Energy dependence of the interferometry scales, R_side

  40. Energy dependence of the interferometry scales, R_out

  41. Energy dependence of ratio

  42. Pion emission density at different energies in HKM

  43. The ratio as function on in-flow and energy Suppose, that the fluid elements, that finally form almost homogeneous in time emission with duration time are initially situated within radial interval shifted to the periphery of the system. If the extreme internal fluid element has initialtransverse velocity and undergo mean acceleration, then .

  44. Pion and kaon interferometry Long- radii (preliminary) NPQCD-2009

  45. Pion and kaon interferometry Side- radii (preliminary) NPQCD-2009

  46. Pion and kaon interferometry Out- radii (preliminary) NPQCD-2009

  47. Hydrodynamics with initial tube-like fluctuations (formation of ridges, initially 10 small tubes)

  48. Hydrodynamics with initial tube-like fluctuations (formation of ridges, initially 1 small tube)

  49. Conclusions-1 49 • The following factors reduces space-time scales of the emission and Rout/Rside ratio. • developing of initial flows at early pre-thermal stage; • more hard transition EoS, corresponding to cross-over; • non-flat initial (energy) density distributions, similar to Gaussian; • early (as compare to standard CF-prescription) emission of hadrons, because escape probability account for whole particle trajectory in rapidly expanding surrounding (no mean-free pass criterion for freeze-out) • Viscosity [Heinz, Pratt] • The hydrokinetic approach to A+A collisions is proposed. It allows one to describe the continuous particle emission from a hot and dense finite system, expanding hydrodynamically into vacuum, in the way which is consistent with Boltzmann equations and conservation laws, and accounts also for the opacity effects.

  50. Conclusions-2 The CFp might be applied only in a generalized form, accounting for the direct momentum dependence of the freeze-out hypersurface corresponding to the maximum of the emission functionat fixed momentum p in an appropriate regionof r.

More Related