[论文解读] Causal Viscous Hydrodynamics for Relativistic Heavy Ion Collisions

Huichao Song|ArXiv.org|Aug 25, 2009

High-Energy Particle Collisions Research参考文献 271被引用 30

一句话总结

本文采用2+1维的因果黏性流体力学模型，基于Israel-Stewart形式化方法，研究RHIC能区重离子碰撞中黏性效应对流致各向异性、粒子谱和熵产生影响。结果表明，黏性效应可将一阶相变平滑为连续转变，并为从实验数据中一致提取η/s提供框架，结果表明QGP的η/s ≈ 0.08是一个合理的值。

ABSTRACT

The viscosity of the QGP is a presently hotly debated subject. Since its computation from first principles is difficult, it is desirable to try to extract it from experimental data. Viscous hydrodynamics provides a tool that can attack this problem and which may work in regions where ideal hydrodynamics begins to fail. This thesis focuses on viscous hydrodynamics for relativistic heavy ion collisions. We first review the 2nd order viscous equations obtained from different approaches, and then report on the work of the Ohio State University group on setting up the equations for causal viscous hydrodynamics in 2+1 dimensions and solving them numerically for central and noncentral Cu+Cu and Au+Au collisions at RHIC energies and above. We discuss shear and bulk viscous effects on the hydrodynamic evolution of entropy density, temperature, collective flow, and flow anisotropies, and on the hadron multiplicity, single particle spectra and elliptic flow. Viscous entropy production and its influence on the centrality dependence of hadron multiplicities and the multiplicity scaling of eccentricity-scaled elliptic flow are studied in viscous hydrodynamics and compared with experimental data. The dynamical effects of using different versions of the Israel-Stewart second order formalism for causal viscous fluid dynamics are discussed, resolving some of the apparent discrepancies between early results reported by different groups. Finally, we assess the present status of constraining the shear viscosity to entropy ratio of the hot and dense matter created at RHIC.

研究动机与目标

• 基于Israel-Stewart形式化方法，发展适用于重离子碰撞的因果黏性流体力学框架。
• 研究剪切黏性和体积黏性对Au+Au与Cu+Cu碰撞中流体演化、流致各向异性和粒子谱的影响。
• 通过比较Israel-Stewart方程的不同版本，解决早期黏性流体力学研究中的结果不一致问题。
• 利用实验数据评估奇异夸克-胶子等离子体剪切黏性与熵密度之比η/s的约束。
• 研究黏性效应对火球中v₂/ε多重性标度及熵产生的影响。

提出的方法

• 基于热力学、输运理论和共形对称性，利用Israel-Stewart形式化方法推导2+1维二阶黏性流体力学方程。
• 采用能量密度、剪切应力π^{μν}和体积压强Π的现实初始条件，实现黏性流体力学方程的数值求解。
• 采用包含一阶相变与混合相的现实状态方程（EOS），并正确处理冻结温度与谱函数计算过程。
• 对π^{μν}和Π采用不同初始化方式，并改变弛豫时间τ_π和τ_Π，以研究初始条件与输运参数的敏感性。
• 比较“简化”与“完整”Israel-Stewart方程的结果，并与Öttinger-Grmela形式化方法进行验证。
• 执行代码验证与多组间交叉比较，确保数值结果的可靠性与一致性。

实验结果

研究问题

• RQ1 剪切黏性和体积黏性如何改变中心与非中心重离子碰撞中熵密度、温度与集体流的流体演化？
• RQ2 黏性效应在多大程度上将一阶夸克-强子相变平滑为连续转变？这对火球寿命与相变动力学有何影响？
• RQ3 Israel-Stewart方程的不同形式化方法（简化版与完整版）如何影响对流致各向异性和粒子谱的预测？
• RQ4 黏性流体力学能否一致地解释v₂/ε的多重性标度？这对QGP的η/s意味着什么？
• RQ5 从实验数据中提取η/s时，主要不确定性来源是什么？体积黏性、初始态模型（Glauber与CGC）以及数值黏性如何影响结果？

主要发现

• 剪切黏性降低纵向压强，减缓初始冷却过程，延长QGP寿命；同时增强横向压强，促进径向流的增强。
• 在边缘碰撞中，由于初始冷却减弱，黏性火球寿命长于理想情况；在中心碰撞中，后期冷却加快，反而缩短总寿命。
• 黏性效应抹平声速中的速度梯度与不连续性，有效将一阶相变转化为平滑的连续转变。
• 在黏性流体力学中，混合相仍受黏性压强梯度驱动而加速，其驻留时间短于理想流体力学中的情况。
• 黏性流体力学能良好再现v₂/ε的多重性标度，支持η/s ≈ 0.08与实验数据一致的结论。
• 不确定性分析表明，体积黏性和初始态模型（Glauber与CGC）是η/s提取中的主要误差来源，数值黏性也对结果产生影响。