The CMS Collaboration / Physics Letters B 789 (2019) 643–665 647
Fig. 3. Ratios of higher order cumulant elliptic-flow harmonics with values obtained from the moments of the unfolded p(v
2
) distributions. Both statistical (lines) and
systematic (gray bands) uncertainties are shown. Hydrodynamic predictions for 2.76 TeV collisions from Ref. [33]are presented as a dark color band and are compared to
the measured v
2
{6}/v
2
{4} ratio. In addition, higher order cumulant ratios reported by the ATLAS Collaboration for 2.76TeV collisions [37]with 0.5 < p
T
< 20.0GeV/c and
|η| < 2.5are compared to the 5.02 TeV measurement. The error bars on the ATLAS measurement represent the quadratic sum of statistical and systematic uncertainties and
points are offset horizontally for clarity.
in Ref. [52]. The cumulant results exhibit the previously observed
v
2
{2} > v
2
{4} ≈ v
2
{6} ≈ v
2
{8} behavior. The centrality-dependent
ratios for the elliptic-flow coefficients obtained for different cumu-
lant
orders are shown in Fig. 3. For most centrality ranges, the
ratios indicate a rank ordering of the cumulants, with differences
on the order of a few percent and with v
2
{4} > v
2
{6} > v
2
{8},
that is qualitatively inconsistent with a pure Gaussian fluctuation
model of flow harmonics. The differences increase as the collisions
become more peripheral. The calculated v
2
{6}/v
2
{4} ratio based
on an event-by-event hydrodynamic calculation using Monte Carlo
Glauber initial conditions [53] and an η/s value of 0.08 is shown
by the shaded band. This simulation is for pions with 0.2 < p
T
<
3.0GeV/c in PbPb collisions at
√
s
NN
= 2.76 TeV [33]. Also shown
are results from the ATLAS Collaboration [37]for PbPb collisions
at 2.76 TeV and for charged particles with 0.5 < p
T
< 20.0GeV/c
and
|η| < 2.5. The calculation is consistent with the experimental
results found at both beam energies. The similarity between exper-
imental
results with 2.76 and 5.02 TeV is consistent with the small
changes in the initial-state eccentricities expected between these
energies [54] and the expectation that the cumulant flow harmonic
ratios follow those of the corresponding eccentricity ratios [33].
Fig. 4 shows the centrality dependence of the standardized
skewness γ
exp
1
. Finite values are found for the standardized skew-
ness
for collisions with centralities greater than ≈15%. The hy-
drodynamic
predictions for the γ
exp
1
values for PbPb collisions at
2.76 TeV from Ref. [33]are also shown and found to be consis-
tent
with the current measurements. Within the hydrodynamic
model and allowing for a finite skewness of the event-by-event v
2
distribution, the small splitting between the cumulant orders is ex-
pected
to follow the relationship (v
2
{6} −v
2
{8})/(v
2
{4} −v
2
{6}) =
0.091 [33]. Experimentally, we find a value for this splitting ra-
tio
of 0.143 ± 0.008 (stat) ± 0.014 (syst) for 20–25% central events,
with the ratio increasing to 0.185 ±0.005(stat) ±0.012(syst) as the
centrality increases to 55–60%. The observed values might suggest
higher order terms in a cumulant expansion of the v
2
distribu-
tion
are required to account for the skewness. This relationship
was recently examined by the ALICE collaboration in Ref. [55]us-
ing
a q-cumulant analysis, with results comparable to the findings
in this paper when considering systematic uncertainties and a dif-
ferent
kinematic range for the ALICE measurement.
Both
elliptic power and Bessel–Gaussian parametrizations used
for fits such as shown in Fig. 1 assume a linear response between
eccentricity and flow, but only the elliptic power law allows for a
finite skewness. For a Bessel–Gaussian distribution, the skewness
is equal to zero. This feature results in the elliptic power function
being in better agreement with the observed fluctuation behavior
Fig. 4. The skewness estimate with respect to the reaction plane determined using
the elliptic-flow harmonic based on different cumulant orders. Both statistical and
systematic uncertainties are shown, where statistical uncertainties are smaller than
the data points. Hydrodynamic model predictions for 2.76 TeV PbPb collisions from
Ref. [33]are shown as a colored band.
than the Bessel–Gaussian parametrization, yielding χ
2
/dof values
on the order of unity. To avoid bin-to-bin correlations introduced
by the unfolding procedure, goodness of fit values are obtained
by refolding the fitted distributions with the response matrix and
comparing to the measured distribution. The elliptic power χ
2
/dof
values
vary between 0.8 and 1.5 from central to peripheral col-
lisions,
while the Bessel–Gaussian χ
2
/dof values vary between 3
and 9. Point-by-point systematic uncertainties on the unfolded dis-
tributions
are correlated and are thus not considered in the fits.
The
fit parameters for the elliptic power function are shown in
Fig. 5 for the different centrality bins. As also found in Ref. [15],
the fits do not converge for central collisions where the distribu-
tions
become very close to a Bessel–Gaussian form. Consequently,
the parameters are shown for centralities >15%. The experimen-
tal
k
2
values show only a weak centrality dependence. Viscous
hydrodynamic calculations indicate that deviations from thermal
equilibrium should lead to a reduced correspondence between the
initial-state geometry and the flow signal in peripheral collisions
[27,28]. This effect is suggested in Fig. 5 by the decrease in the
k
2
value with increasing centrality, although the systematic uncer-