跳到主要内容
Sandbox Physics
重离子物质与相对论流体

FR06 · L3 · 第 8 章,共 9 章

从粒子关联里读出共同的流动

每次碰撞偏爱的方向都不相同。把所有事件的粒子直接混在实验室坐标里平均,流动花纹可能消失。同一事件内部的粒子夹角却保留了共同方向的信息,让我们不用先知道事件朝向,也能估计这种集体起伏。

读前准备第 7 章定义了角谐波。这里需要条件独立和逐事件平均。复数记号只是把旋转单位箭头相加,取实部就回到余弦。

Sandbox Physics ·

先消去每次碰撞未知的朝向

在一次事件中,固定谐波强度和共同朝向 Ψn\Psi_n,设两粒子在这个条件下独立发射。它们各自满足 ⟨einϕ∣Ψn⟩=vneinΨn\langle e^{in\phi}\mid\Psi_n\rangle=v_ne^{in\Psi_n}。一个均值乘另一个的复共轭,共同角度就被消去了。

⟨cos⁡n(ϕ1−ϕ2)∣vn,Ψn⟩=vn2\left\langle\cos n(\phi_1-\phi_2)\mid v_n,\Psi_n\right\rangle=v_n^2
等式需要条件独立。衰变粒子对等相关发射,可能额外贡献无法拆成两个单粒子均值的项。

跨事件平均后,写成 cn{2}=⟨vn2⟩+δnc_n\{2\}=\langle v_n^2\rangle+\delta_n。额外的 δn\delta_n 叫非流,指共同单事件调制以外的角关联。流的贡献是一个小量的平方,因此非流项即使很小,也可能影响明显。

计算时必须把粒子与自己配对的项去掉

一次事件有 MM 条径迹,先加单位箭头 Qn=∑j=1MeinϕjQ_n=\sum_{j=1}^M e^{in\phi_j}。模平方把所有有序粒子对都加进去,也带上 MM 个粒子与自身的配对,每项余弦都等于一。减去它们,再除以不同粒子的有序对数。

⟨2⟩=∣Qn∣2−MM(M−1)\langle2\rangle=\frac{|Q_n|^2-M}{M(M-1)}
对保存的均匀权重样本,这个恒等式精确成立,不需要大粒子数,也不需要已知事件面。

各事件粒子数不同时,还要说明怎样给事件加权。冻结的教学基准让每个事件粒子数相同,算术平均与按粒子对数加权一致。模型也直接设定探测接受度均匀;真实分析则要校正接受度和效率。

不枚举全部四元组,也能算出同一个量

保存的 Q2nQ_{2n} 箭头和,可以帮助减掉四粒子式中重复索引的项。在 M>3M>3 时,精确分子和归一化为:

N4=∣Qn∣4+∣Q2n∣2−2Re⁡ ⁣[Q2n(Qn∗)2]−4(M−2)∣Qn∣2+2M(M−3),⟨4⟩=N4M(M−1)(M−2)(M−3).\begin{aligned}N_4={}&|Q_n|^4+|Q_{2n}|^2-2\operatorname{Re}\!\left[Q_{2n}(Q_n^*)^2\right]\\&-4(M-2)|Q_n|^2+2M(M-3),\\\langle4\rangle={}&\frac{N_4}{M(M-1)(M-2)(M-3)}.\end{aligned}
复共轭把箭头相位反向。这个代数恒等式排除了所有重复径迹,改变的是计算量,平均的不同粒子四元组仍相同。

检验时用小径迹样本直接枚举四元组,和快速公式比较。这项独立计算可以查出漏减自关联的错误。

四个粒子探测的是另一个矩

对各不相同的四个粒子平均 cos⁡n(ϕ1+ϕ2−ϕ3−ϕ4)\cos n(\phi_1+\phi_2-\phi_3-\phi_4),再跨事件平均,记作 ⟨ ⁣⟨4⟩ ⁣⟩\langle\!\langle4\rangle\!\rangle。累积量扣去两份双粒子平均的平方。

cn{4}=⟨ ⁣⟨4⟩ ⁣⟩−2⟨ ⁣⟨2⟩ ⁣⟩2,vn{4}=(−cn{4})1/4  (cn{4}<0)c_n\{4\}=\langle\!\langle4\rangle\!\rangle-2\langle\!\langle2\rangle\!\rangle^2,\qquad v_n\{4\}=(-c_n\{4\})^{1/4}\ \ (c_n\{4\}<0)
若流强度固定且粒子独立发射,四粒子矩就是流强度的四次方,累积量便等于它的负值。

有逐事件起伏时,四阶矩一般不等于二阶矩的平方。即使只有集体流,两粒子与四粒子的结果也可能不同。累积量能压低特定的低阶关联,有限粒子数、成簇发射和探测接受度仍会留下其他贡献。

若估计出的四粒子累积量为正,这套约定下就没有实数的四次根流值。应报告符号和样本误差;直接取绝对值,会悄悄换掉所估计的统计量。

完整例题

少量双粒子簇,足以移动双粒子信号

合成基准每事件有 M=128M=128 条径迹,其中四对各有两条完全相同角度的径迹。在 128⋅127128\cdot127 个有序不同粒子对中,八个来自这些双粒子簇,比例为 d=8/(128⋅127)≃0.00049213d=8/(128\cdot127)\simeq0.00049213。

E[⟨2⟩]=vn2+d(1−vn2)\mathbb E[\langle2\rangle]=v_n^2+d(1-v_n^2)
每个簇的朝向也从相同单粒子分布抽取。簇内两条径迹的夹角余弦为一;不同独立发射之间的粒子对,期望仍是流强度的平方。

取 vn=0.06v_n=0.06,纯流双粒子量为 0.00360.0036,加上这些簇后约为 0.0040900.004090。若忽略非流,表观双粒子流强度约为 0.063950.06395。非流比例相对一很小,相对平方后的信号却并不小。

先预测,再操作模型

保留同一批粒子,只把分析中假设的非流设为零,已保存的关联量会变吗?

打开输运推断,比较“生成数据时的假设”和“忽略双粒子背景”,展开“保存的事件与四粒子诊断”,检查假设改变时已保存的双关联和四阶累积量是否变化。

输运参数推断 ↗
对照解释

样本始终包含相同双粒子簇。忽略它们改的是推断假设,粒子并未改变。保存的四阶累积量仍有有限样本起伏和簇贡献。

练习与解释性答案

1. 取 M=2M=2,Q 向量的双粒子式会返回什么?

查看思路

∣Qn∣2=2+2cos⁡n(ϕ1−ϕ2)|Q_n|^2=2+2\cos n(\phi_1-\phi_2)。减去二再除以二,恰好得到唯一不同粒子对的夹角余弦。这是一个方便的小样本检验。

2. 事件面均匀随机时,直接跨事件平均单粒子余弦,为什么会丢掉花纹?

查看思路

实验室里每个偏爱方向同样可能,复相位平均为零。同一事件的粒子角度差会先消去共同相位,之后再跨事件平均就能保留信号。

来源与核读范围

  1. Gardim et al., Mapping the hydrodynamic response to the initial geometry in heavy-ion collisions, arXiv:1111.6538v2

    2012-03-05 第二版,§II–IV:偏心率、线性响应及高次谐波的非线性项。教学用指数衰减式另有定义,未由论文事件模拟拟合。

  2. Bilandzic, Snellings & Voloshin, Flow analysis with cumulants: direct calculations, arXiv:1010.0233v2

    2011-10-17 第二版,式 (3)–(4)、(11)–(18):非流、Q 向量及均匀接受度下排除自关联的精确多粒子关联。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策