Shear viscosity of a gluon plasma in heavy ion collisions

1. Shear viscosity of a gluon plasma in heavy ion collisions Qun Wang Univ of Sci & Tech of China With J.W. Chen, H. Dong, K. Ohnishi, arXiv:0907.2486 Dense Matter in Heavy-Ion Collisions and Supernovae, Preprow, Oct. 10-14,2009

2. What is viscosity related to HIC viscosity = resistance of liquid to viscous forces (and hence to flow) Shear viscosity Navier 1822 Bulk viscosity

3. Shear viscosity in ideal gas and liquid ideal gas, high T liquid, low T lower bound by uncertainty principle Frenkel, 1955 Danielewicz, Gyulassy, 1985 Policastro,Son,Starinets, 2001 T. Schaefer's talk

4. η/s around phase transition Lacey et al, PRL98, 092301(2007) Csernai, et al PRL97, 152303(2006) T. Schaefer's talk

5. Previous results on shear viscosity for QGP ►PV: Perturbative and Variational approach Danielewicz, Gyulassy, Phys.Rev.D31, 53(1985) Dissipative Phenomena In Quark Gluon Plasmas Arnold, Moore and Yaffe, JHEP 0011, 001 (2000),0305, 051 (2003) Transport coefficients in high temperature gauge theories: (I) Leading-log results (II): Beyond leading log ........... ►BAMPS: Boltzmann Approach of MultiParton Scatterings Xu and Greiner, Phys. Rev. Lett. 100, 172301(2008) Shear viscosity in a gluon gas Xu, Greiner and Stoecker, Phys. Rev. Lett. 101, 082302(2008) PQCD calculations of elliptic flow and shear viscosity at RHIC ►Different results of AMY and XG for 2↔3 gluon process: ~(10-20)% (AMY) ~ (70-90)% (XG)

6. Difference: AMY vs XG Both approaches of XG and AMY are based on kinetic theory. However, the main points of differences are: 1) A parton cascade model is used by XG to solve the Boltzmann equation. Since the bosonic nature of gluons is hard to implement in real time simulations in this model, gluons are treated as a Boltzmann gas (i.e. a classical gas). For AMY, the Boltzmann equation is solved in a variation method without taking the Boltzmann gas approximation. 2) The Ng↔ (N+1)g processes, N=2,3,4,..., are approximated by the effective g↔gg splitting in AMY with 2-body-like phase space, while the Gunion-Bertsch formula for gg↔ggg process is used in XG with 3-body-like phase space.

7. Our goal and strategy Goal: to calculate the shear viscosity in a different way, to understand the nature of the difference between two results Strategy: 1) We use variational method as AMY 2) We use the Gunion-Bertsch formula for gg↔ggg process as XG 3) For evaluating collisional integrals we treat phase space for 3 gluons in two ways: (a) 3 body state as XG; (b) 2+1(soft) state, almost 2 body state, close to AMY. We call it the soft gluon approximation;

8. Boltzmann equation for gluon plasma gluon distribution function gg↔gg collision terms gg↔ggg collision terms phase-space measure matrix element delta function EM conservation [ gain - loss ]

9. q Matrix elements: gg↔gg and gg↔ggg Soft-collinear approximation gg↔ggg, factorized form, Gunion-Bertsch, PRD 25, 746(1982) k q

10. Shear viscosity: variational method perturbation in distribution function linear in χ(x,p)

11. Shear viscosity: variational method shear viscosity in terms of B(p) solve χ(x,p) by Boltzmann eq. → the constraint for B(p) S. Jeon, Phys. Rev. D 52, 3591 (1995); Jeon, Yaffe, Phys. Rev. D 53, 5799 (1996).

12. Shear viscosity: variational method Inserting eq for B(p) into shear viscosity, quadratic form in B(p) B(p) can be expanded in orthogonal polynomials orthogonal condition

13. Shear viscosity: variational method Inserting eq for B(p) into shear viscosity, quadratic form in B(p)

14. Collisional rate Boltzmann equation written in Collisional rate is defined by

15. Regulate infrared and collinear divergence for kT in gg↔ggg ■ Landau-Pomeronchuk-Migdal (LPM) effect by cutoff (used by Xu-Greiner and Biro et al) ■ Debye mass m_D as the gluon mass or regulator (used by Arnold-Moore-Yaffe)

16. Importance of phase space for gg↔ggg ■ almost 3-body (3-jet) phase space (used by Xu-Greiner) ■ almost 2-body phase space (used by Arnold-Moore-Yaffe) treated as equal footing phase space dim: ~ 3X3-4=5 soft splitting function is used phase space dim: ~ 2X3-4=2 polar and azimuthal angles, (θ,φ) colinear

17. Importance of phase space for gg↔ggg ■Soft gluon approximation in our work (as one option of our calculation) Emission of the 5th gluon does not influence the configuration of 22 very much, therefore gg↔ggg can be factorized into gg↔gg and g↔gg This is just the way Gunion-Bertsch got their formula. → Phase space dim: ~ 2X3-4=2, polar and azimuthal angles, (θ,φ) This is equivalent to exand Jacobian of δ(E1+E2-E3-E4-E5)in large √s limit and keeping the leading order. For the form of Jacobian, see Appendix D of Xu, Greiner, PRC71, 064901(2005).

18. Soft gluon approximation in cross section of gg↔ggg Eq.(D5), Xu & Greiner, PRC71, 064901 (2005) Biro, et al, PRC48, 1275 (1993) two roots: y' (forward), -y' (backward) keep only positive root for y': a factor 1/2

19. Our results

20. Leading-Log result for gg↔gg We reproduced AMY's leading-log(LL), For Boltzmann gas, LL result: Our numerical results show good agreement to LL result in weak coupling

21. η22: Bose and Boltzmann gas

22. Collisional rates

23. Shear viscosity from 22 and 23 process

24. Numerical results: shear viscosity

25. Comparison : AMY, XG, Our work

26. Concluding remarks ■ We have bridged to some extent the gap between AMY and XG. ■ To our understanding, their main difference is in the phase space for number changing processes, there are much more 3-body configurations in XG approach than in AMY, or equivalently phase space in XG for gluon emission is much larger than in AMY (about dim 5 : dim 2), causing effect of 23 for viscosity in XG is much larger than in AMY. ■Core question:Is GB formula still valid for general 3-body (3-jet) configuration? or equivalently: Does GB formula over-estimate the rate of the general 3-body (3-jet) configuration? Further study of viscosity using exact matrix element should give an answer to this question. ■Viscosity in association with 2nd order hydrodynamics can be studied in the same way. [XG, PRC 2008]

27. THANK YOU !

28. Orthogonal polynomials: B^r Dimensionless polynomial for y=1 orthogonal condition