300 likes | 322 Views
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, Phys.Lett.B685 , 277 ( 2010 ) The 8- th conference of chinese high energy physics society, Nanchang, April 16-21, 2010.
E N D
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, Phys.Lett.B685, 277(2010) The 8-th conferenceof chinese high energy physics society, Nanchang, April 16-21,2010
What is viscosity related to HIC viscosity = resistance of liquid to viscous forces (and hence to flow) Shear viscosity Navier 1822 Bulk viscosity
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
η/s around phase transition Lacey et al, PRL98, 092301(2007) Csernai, et al PRL97, 152303(2006)
ζ/s around phase transition Bernard et al, (MILC) PRD 2007, Cheng et al, (RBC-Bielefeld) PRD 2008, Bazavov et al, (HotQCD),arXiv:0903.4379 Karsch, Kharzeev, Tuchin, PLB 2008 Noronha *2, Greiner, 2008, Chen, Wang, PRC 2009, B.C.Li, M. Huang, PRD2008, ......
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)
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.
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;
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 ]
q Matrix elements: gg↔gg and gg↔ggg Soft-collinear approximation gg↔ggg, factorized form, Gunion-Bertsch, PRD 25, 746(1982) k q
Shear viscosity: variational method perturbation in distribution function linear in χ(x,p)
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).
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
Shear viscosity: variational method Inserting eq for B(p) into shear viscosity, quadratic form in B(p)
Collisional rate Boltzmann equation written in Collisional rate is defined by
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)
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
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).
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
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
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.
Orthogonal polynomials: B^r Dimensionless polynomial for y=1 orthogonal condition