[论文解读] Viscous fluid dynamics in Au+Au collisions at RHIC

A. K. Chaudhuri|ArXiv.org|Jan 21, 2008

High-Energy Particle Collisions Research参考文献 1被引用 15

一句话总结

本文利用以色列-斯图尔特形式化方法，通过 η/s = 0.08 的剪切黏度模型，研究了RHIC能量下Au+Au碰撞中的黏性流体动力学。结果表明，最小黏性流体动力学在 pT ≈ 3.6 GeV 以内能更准确地再现实验的 pT 谱和椭圆流，而理想流体动力学在 pT ≈ 1.5 GeV 以上即失效，凸显了黏性在描述不同中心度下RHIC数据中的关键作用。

ABSTRACT

We have studied the space-time evolution of minimally viscous ($\fracη{s}$=0.08) QGP fluid, undergoing boost-invariant longitudinal motion and arbitrary transverse expansion. Relaxation equations for the shear stress tensor components, derived from the phenomenological Israel-Stewart's theory of dissipative relativistic fluid, are solved simultaneously with the energy-momentum conservation equations. Comparison of evolution of ideal and viscous fluid, both initialized under the similar conditions, e.g. same equilibration time, energy density and velocity profile, indicate that in viscous fluid, energy density or temperature of the fluid evolve slowly than in an ideal fluid. Transverse expansion is also more in viscous evolution. We have also studied particle production in viscous dynamics. Compared to ideal dynamics, in viscous dynamics, particle yield at high $p_T$ is increased. Elliptic flow on the other hand decreases. Minimally viscous QGP fluid, initialized at entropy density $s_{ini}$=110 $fm^{-3}$ at the initial time $τ_i$=0.6 fm, if freeze-out at temperature $T_F$=130 MeV, explains the centrality dependence of $p_T$ spectra of identified particles. Experimental $p_T$ spectra of $π^-$, $K^+$ and protons in 0-5%, 5-10%, 10-20%, 20-30%, 30-40% and 40-50% Au+Au collisions are well reproduced through out the experimental $p_T$ range. This is in contrast to ideal dynamics, where, the spectra are reproduced only up to $p_T\approx$1.5 GeV. Minimally viscous QGP fluid, also explain the elliptic flow in mid-central (10-20%, 16-23%, 20-30%) collisions. The minimum bias elliptic flow is also explained. However, the model under-predict/over-predict the elliptic flow in very central/peripheral collisions.

研究动机与目标

• 评估有限剪切黏度对RHIC中Au+Au碰撞中QGP时空演化的影响。
• 解决理想流体动力学在描述高 pT 谱和饱和椭圆流方面存在的局限性。
• 检验最小黏性流体动力学（η/s = 0.08）是否能同时解释多个中心度区间内 pT 谱和椭圆流。
• 确定黏性效应是否为解决理想流体动力学模型中差异（特别是HBT疑难和高 pT 偏差）所必需。
• 利用实验数据为QGP输运系数提供现象学约束。

提出的方法

• 采用相对论性耗散流体动力学的以色列-斯图尔特理论，求解剪切应力张量分量的弛豫方程。
• 在双胞胎不变和任意横向膨胀条件下，同时求解能量-动量守恒方程与剪切应力演化方程。
• 初始条件设定为与理想流体动力学一致：τi = 0.6 fm，s_ini = 110 fm⁻³，vx = vy = 0，且 η/s = 0.08。
• 弛豫时间 τπ 固定为玻尔兹曼气体值：τπ = 6η/sT。
• 通过调整冻结温度以拟合实验数据，发现在16–23%中心度下，T_F = 130 MeV 为椭圆流的最佳值。
• 使用相同的黏性流体动力学设置计算粒子谱和椭圆流，并与PHENIX和STAR实验数据进行比较。

实验结果

研究问题

• RQ1 最小黏性流体动力学（η/s = 0.08）是否能比理想流体动力学更准确地描述Au+Au碰撞中识别粒子的 pT 谱的中心度依赖性？
• RQ2 黏性动力学是否能改善对椭圆流的描述，特别是在理想流体动力学失效的高 pT 区域？
• RQ3 尽管中等中心度数据拟合良好，为何该模型在极中心碰撞中低估椭圆流，而在外围碰撞中高估椭圆流？
• RQ4 与理想流体动力学相比，剪切黏度在多大程度上改变了能量密度和横向膨胀的演化？
• RQ5 一组参数（η/s = 0.08，T_F = 130 MeV）是否能同时在多个中心度类别中准确描述 pT 谱和椭圆流？

主要发现

• 采用 η/s = 0.08 和 T_F = 130 MeV 的最小黏性流体动力学，在所有中心度区间（0–5% 至 40–50%）内成功再现了 π⁻、K⁺ 和质子的 pT 谱，将一致性扩展至 pT ≈ 3.6 GeV，而理想流体动力学在 pT ≈ 1.5 GeV 以上即失效。
• 该模型成功解释了中等中心度碰撞（10–20%、16–23%、20–30%）及最小偏差数据的椭圆流中心度依赖性，一致性可达 pT ≈ 3.6 GeV。
• 在极中心碰撞中，模型低估了椭圆流，表明可能需要更低的黏性；而在外围碰撞中则高估了椭圆流，提示可能需要更高的黏性。
• 黏性动力学导致能量密度和温度演化变慢，并增强了横向膨胀，与理想流体动力学相比。
• 该模型还成功再现了 π⁻、K⁺ 和质子的平均横向动量的中心度依赖性，直至 N_part ≈ 100，表明其具有稳定的流体动力学行为。
• 结果表明，即使剪切黏度极小（η/s = 0.08），也能显著改善对RHIC数据的描述，特别是在高 pT 和椭圆流观测量方面，凸显了在流体动力学建模中包含黏性效应的必要性。