Extraction of the specific shear viscosity of quark-gluon plasma from two-particle transverse momentum correlations

The specific shear viscosity, $\eta/s$, of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions at RHIC and LHC is estimated based on the progressive longitudinal broadening of transverse momentum two-particle correlators, $G_2$, reported as a function of collision centrality by the STAR and ALICE experiments. Estimates are computed as a function of collision centrality using the Gavin ansatz which relates the $G_2$ longitudinal broadening to the specific shear viscosity. Freeze out times required for the use of the ansatz are computed using a linear fit of freeze out times reported as a function of the cubic root of the charged particle pseudorapidity density $({\rm d}N_{\rm ch}/{\rm d}\eta)^{1/3}$. Estimates of $\eta/s$ based on ALICE data exhibit little to no dependence on collision centrality at LHC energy, while estimates obtained from STAR data hint that $\eta/s$ might be a function of collision centrality at top RHIC energy.


Introduction
A key focus of the ultrarelativistic heavy-ion collision programs conducted at the Large Hadron Collider (LHC) and the Relativistic Heavy-Ion Collider (RHIC) involves precision measurements of the properties of the quark-gluon plasma (QGP) formed in high-energy nucleus-nucleus collisions. Of particular interest are the magnitude and temperature dependence of the specific shear viscosity of the QGP, expressed as the ratio η/s of the shear viscosity η to the entropy density s of the matter produced in the collisions. Shear viscosity characterizes the ability of a medium to transport momentum and carry deformations. Transverse particle anisotropy patterns, quantified in terms of anisotropic a e-mail: victor.gonzalez@cern.ch b e-mail: sumit.basu@cern.ch c e-mail: a.marin@gsi.de d e-mail: claude.pruneau@wayne.edu flow coefficients, measured in mid central heavy-ion collisions at RHIC and LHC are rather large and were, from the onset, relatively well reproduced by viscosity free hydrodynamical calculations thereby suggesting the QGP might be a perfect fluid, i.e., a fluid with vanishing or negligible shear viscosity [1][2][3][4][5][6]. The possibility that the high temperature, high density systems formed in the midst of heavy-ion collisions might be a perfect fluid thus generated 'quite' a bit of excitement [7]. Considerable experimental and theoretical efforts were consequently expanded to determine the specific shear viscosity of the matter produced at RHIC and more recently at the LHC [8,9] based on measurements of anisotropic flow in the collision transverse plane. Although theoretical efforts have been quite successful in reducing the range of η/s compatible with state of the art measurements of flow anisotropies, there still remains a certain degree of ambiguities owing to several technical difficulties. One of these technical difficulties involves the lack of knowledge on the initial conditions of the systems produced in A-A collisions at RHIC and LHC [5,10,11] . For instance, CGC inspired initial conditions yield larger initial spatial anisotropy than MC Glauber type initial conditions and thus require a somewhat larger level of shear viscosity to match the observed flow coefficients when used as input to viscous hydrodynamics simulations. While efforts to reduce the initial conditions ambiguity based on measurements of symmetric cumulants [12,13], in particular, have had some success, it remains of interest to identify techniques that might enable measurements of specific shear viscosity that are less susceptible to uncertainties associated with initial conditions. Such a technique exists. Proposed by Gavin et al. already more than a decade ago [14], it involves measurements of the longitudinal broadening of a transverse momentum two particle correlator, now dubbed G 2 , with increasing collision centrality. The correlator G 2 , defined in Ref. [14,15], is designed to be proportional to the covariance of momentum currents and is as such sensitive to dissipative viscous forces at play during the transverse and longitudinal expansion of the matter formed in A-A collisions. Gavin et al. showed these forces lead to a longitudinal broadening of G 2 measured as a function of the pseudorapidity difference of measured charged particles. As the matter expands, neighboring fluid cells drag one another. Fast fluid cells tend to slow down whereas slow fluid cells accelerate. This has the effect of dampening the expansion and produces a progressive broadening of the G 2 correlator with time. The longer the system lives, the longer viscous effects play a role, and the broader the G 2 correlator becomes. Gavin et al. showed the broadening, characterized in terms of the difference of the variance of the correlator observed in most central and most peripheral collisions, should be proportional to η/s and given by the following formula herein called the Gavin ansatz where σ c is the longitudinal width of the correlator measured in most central collisions whereas σ 0 is the longitudinal width of the correlator at formation time τ 0 . T c and τ c,f are respectively the critical temperature and the freeze-out time in most central collisions.
We first briefly review, in sec. 2, prior efforts to determine η/s based on the longitudinal broadening of the G 2 correlator in A-A collisions. The method and results of this work are presented in sec. 3 and discussed in sec. 4. Our conclusions are presented in sec. 5.
2 Prior estimates of η/s based on transverse momentum correlations A first estimate of the QGP viscosity based on the Gavin ansatz was reported several years ago by the STAR collaboration using a measurement of the charge independent correlator G CI 2 in Au-Au collisions at √ s NN = 0.2 TeV. STAR observed the longitudinal width of the G 2 correlator grows considerably from most peripheral to most central Au -Au collisions. Given the observed broadening might arise in part from other dynamical effects, STAR used the Gavin ansatz to estimate an upper limit and reported η/s to be in the range 0.06 -0.21 [16]. More recently, the ALICE collaboration reported precise measurements of the evolution of the longitudinal and azimuthal widths of charge independent and charge dependent two-particle transverse momentum correlators, G CI 2 and G CD 2 , respectively, as a function of the centrality of Pb -Pb collisions at √ s NN = 2.76 TeV [17]. Examining specifically the overall change of the correlator longitudinal width from most peripheral to most central collisions, the collaboration concluded that their observations favour small values of η/s, that is, values close to the KSS bound of 1/4π [18].
The η/s estimates reported by the ALICE and STAR collaborations focused on the overall change of the longitudinal width G CI 2 from peripheral to central collisions but did not utilize correlator widths observed in mid-central collisions. This omission resulted in large part from the lack of precise estimates of the system's life time in mid-central collisions. Effectively, STAR and ALICE did not consider the possibility that the viscosity might evolve with collision centrality and thus did not apply the Gavin ansatz to intermediate ranges of collision centralities. We note, however, that the viscosity might in fact become a function of the collision centrality if, in particular, the temperature or density of the produced system or other conditions affecting the viscosity evolve with centrality. It is also conceivable that other aspects of the collision dynamics, not related to viscous effects, could impact the broadening of the G 2 vs. centrality. It is thus of interest to consider what the evolution of the G 2 correlator observed by STAR and ALICE implies. Two specific questions arise. The first is concerned purely with the experimental technique used to estimate η/s while the second concerns a possible evolution of the effective shear viscosity of the system with collision centrality.
Let us first consider the experimental technique on its own merits. Is the technique sound? Are there experimental artifacts that can bias or skew the evaluation of η/s based on the Gavin ansatz? Indeed, the ansatz requires estimates of a critical temperature T c , as well as initial (formation) and freeze-out times τ 0 and τ c,f , respectively. These quantities are not evaluated in the context of the G 2 measurement and thus require external inputs. They may thus be subjected to systematic bias of their own and independent of the STAR and ALICE measurements of the G 2 correlator. Additionally, estimation of the broadening of the correlator might perhaps be biased by the finite acceptance or other artifacts of the measurement process. One might wonder, in particular, whether the width observed in most central collisions could be underestimated because of the finite rapidity width of the acceptance of the measurements. In this context, it becomes of interest to study what progressive changes of the width might imply about the strength of the specific shear viscosity, and whether, in particular, the evolution of the widths with centrality is self-consistent, that is, whether changes of the width from one fractional cross section to the next are consistent with the overall change from most peripheral to most central collisions.
The second set of concerns is of greater interest, from a physical standpoint, but perhaps more difficult to elucidate. Are viscous effects strictly proportional to the system lifetime? Can the correlator be affected by other physical effects, such as, possibly, the radial and anisotropic expansion of the collision system? Is the characteristic temperature used in the ansatz truly a constant independent of the collision centrality? And perhaps, most interestingly, could the effective shear viscosity extracted from the measurement be a function of collision centrality? Theoretical considerations suggest η/s is likely a function of the QGP temperature [13,[19][20][21][22][23][24][25]. Is it then possible that collisions at different impact parameter yield systems at different temperatures with slightly different time evolution of the shear viscosity, thereby resulting in effective or time-averaged shear viscosity that might depend on the collision centrality? Conceivably, answers to these questions may require more and better data than those available, but it is nonetheless of interest to consider what the available data can say about a possible evolution of η/s with collision centrality and system temperature. It is thus the primary objective of this work to explore how η/s values obtained with the Gavin ansatz evolve with collision centrality.

Evolution of η/s with system size
We proceed with the evaluation of η/s as a function of the cubic root of the pseudorapidity density dN ch /dη, based on the G CI 2 longitudinal widths already reported by the STAR and ALICE collaborations [16,17] in Au-Au collisions at √ s NN = 0.2 TeV, and in Pb-Pb at √ s NN = 2.76 TeV, respectively, as a function of the collision centrality using Eq. (1). However, we also need estimates of the lifetimes τ f of the system with collision centrality. Estimates of freeze-out times reported as a function of measured charged particle densities, dN ch /dη, in [26], are used. Values of τ f are obtained from two-pion Bose-Einstein measurements from AGS to LHC energies [26]. Freeze-out times relevant for each of the centrality ranges considered in this work are obtained by fitting a first degree polynomial to estimated values τ f according to Freeze-out times and associated uncertainties are listed in Tab. 1 and plotted in Fig. 1. Statistical and systematic errors reported by the E895, CERES, NA49, PHOBOS, STAR and ALICE collaborations [26][27][28][29][30][31][32][33][34][35][36] are used in the least square fit procedure. Note, however, that we could not obtain a fit with Eq. (2) that satisfactorily match all available data. We thus proceeded to use linear polynomial fits (τ = a 0 + a 1 dN ch /dη) and opted to give larger emphasis and weights to τ estimates obtained at 0.2 and 2.76 TeV given our goal is to determine η/s based on G 2 data acquired at these two energies. Several distinct fits were carried out to obtain parameterizations of the τ f dependence on the charged particle density used in our determination of η/s. Fit conditions were varied: our primary fit included all data points but we also considered fits based on data in selected energy ranges, and with or without constraining the fits to pass through the origin. All fits considered yield chi-square per degrees of freedom of the order of χ 2 /ndf = 2. For the STAR energy, 0.2 TeV, the fit that better reproduces the published results does not pass through the origin and yields A = 0.72 while for the ALICE energy, 2.76 TeV, the best fit yields a straight line that passes through the origin and a value A = 0.88. Fig. 1 shows the extrapolated decoupling times τ f corresponding to the charged particle density used in the measurements of the G 2 correlator. Error bands show the systematic uncertainties introduced by the fit procedure. Additionally, we also carried out estimates of η/s based on τ f values obtained from pion interferometry and blast wave fits to particle spectra [27,[37][38][39] and found that estimates η/s obtained with these alternative values of τ f were in agreement, within uncertainties, with the results obtained with τ f values obtained with the fit procedure described above.
Computation of the Gavin ansatz is accomplished using the canonical values T c = 160 ± 5 MeV and τ 0 = 1.0 ± 0.5 fm/c for the critical temperature and formation time, respectively [40]. The value of σ 0 is estimated by extrapolating the width of the correlator to N part = 2. Values of (dN ch /dη) 1/3 for a given centrality class are taken from [36,41]. The STAR and ALICE collaborations estimated shear viscosities, using Eq. 1, based exclusively on most central and most peripheral collisions. In this letter, the collision centrality dependencies of the longitudinal widths reported by both experiments are used to investigate whether η/s exhibits a dependence on (dN ch /dη) 1/3 . The longitudinal broadening of the G 2 correlator, defined below, is expected to be insensitive to initial state density fluctuations in the transverse plane. As such, it powerfully complements studies of η/s based on measurements of anisotropic flow that suffer in part from such a dependence [25]. Measurements of G 2 correlators additionally have a different sensitivity to non-flow effects which make them an invaluable tool in the understanding of the dynamics of A-A collisions, and as such, provide additional testing grounds of hydrodynamical and other types of theoretical models.
The ALICE measurements were reported in terms of a dimensionless variant of the G 2 correlator [14,15] defined as with where n 1 ≡ n(η 1 , ϕ 1 ) and n 2 ≡ n(η 2 , ϕ 2 ) are the number of charged particle tracks detected, in each event, within bins centered at η 1 , ϕ 1 and η 2 , ϕ 2 , respectively. Sums are carried over particle transverse momenta p T,i , i ∈ [1, n 1 ], and p T, j , j = i ∈ [1, n 2 ], respectively. The bracket notation O is used to represent event ensemble averages computed within the bins η i , ϕ i , i = 1, 2. Thus n i and p T,i represent average number of particles and average transverse momenta in bin η i , ϕ i , respectively. The reported ALICE measurement was limited to charged particles with transverse momenta in the range 0.2 ≤ p T < 2 GeV/c and pseudorapidities within |η| < 0.8. Within the context of the ALICE analysis, the G CI 2 correlator dependence on ∆ η and ∆ ϕ was parametrized with a two-component model defined as where B and a n describe the long-range mean correlation strength and azimuthal anisotropy, respectively, while the bidimensional generalized Gaussian, whose shape is determined by the parameters A, ω ∆ η , ω ∆ ϕ , γ ∆ η and γ ∆ ϕ , is used to model the correlation signal of interest. The ALICE collaboration reported longitudinal widths σ ∆ η computed as the standard deviation of the generalized Gaussian along ∆ η [17]. These values are plotted as a function of the number of collision participants estimated from Glauber models [42] in Fig. 2. Instead of using a fitting procedure, the STAR collaboration estimated the longitudinal width of measured correlators by computing the rms of one-dimensional projections of G 2 correlators onto the ∆ η axis [16]. These rms width values are plotted vs. the number of collision participants in Fig. 2.
All parameters (dN ch /dη, G CI 2 widths and τ f ) used in the computation of η/s as well as the extracted values of η/s at RHIC and LHC energies are listed in Tab. 1 as a function of collision centrality expressed in terms of the fractional cross section: the range 0-5% corresponds to most central collisions while ranges 5-10%, 10-20%, etc, represent collisions with increasingly larger impact parameters. Our analysis, based on STAR and ALICE is limited to quasiperipheral collisions up to the range 70-80%, beyond which the applicability of Gavin's model might be put into question. Estimates of η/s computed with Eq. 1 based on the above widths and freeze-out times are listed, for both RHIC and LHC energies, in the two right-most columns of Tab. 1 . The values of η/s vs ( rmdN ch /dη) 1/3 are plotted in Fig. 3. Statistical and systematic uncertainties, the last ones incorporating the uncertainties from the fit procedure for freezeout times extraction, are displayed with vertical bars and rectangular boxes, respectively.      Computed values of η/s range from 0.04 ± 0.02 sys to 0.07 ± 0.03 sys and from a value compatible with 0 to 0.16 ± 0.03 sta ± 0.07 sys for LHC energies and RHIC energies, respectively. One observes that values extracted from Pb-Pb collisions at the LHC exhibit a weak dependence on (dN ch /dη) 1/3 , while those from Au-Au collisions, measured at RHIC, show a rising trend with increasing (dN ch /dη) 1/3 , albeit with large uncertainties. Values of η/s obtained from Pb-Pb collisions are close but somewhat lower than the KSS bound of 1/4π , while those obtained from Au-Au colli-sions are compatible with vanishing viscosities in the range (dN ch /dη) 1/3 < 5 but exceed the KSS bound above (dN ch /dη) 1/3 ≈ 6. However, η/s values derived from STAR and ALICE data are compatible with one another at the one σ level at all values of (dN ch /dη) 1/3 .

Discussion
A compilation of η/s values obtained in this work as well as those reported in theoretical and phenomenological calculations is presented in Fig. 4. Shear viscosity values for ultracold Fermi gas [44], Helium [56] and water [56], evaluated at their respective critical temperatures, as well as the holographic bounds [45] and the KSS limit [18], have also been incorporated as baselines and references.
The results reported in this work are displayed in Fig. 4 with two horizontal lines (one for the LHC and one for RHIC). The ranges of the horizontal lines span the smallest to highest values obtained at each energy, including systematic uncertainties. Somewhat older compilations have also been reported [62]. A comparative analysis of the results presented in this compilation is complicated in part by the fact that viscous effects are likely to accumulate throughout a system's evolution. However, the shear and bulk viscosities may depend on the temperature, matter density, the presence of magnetic fields, and possibly other system conditions, that evolve as the QGP expands and goes through a transition into a hadron phase. The model used in this work and several of the calculations listed in the compilation neglect such a time/temperature dependence and represent the viscosity as a sin-  gle effective value, while others attempt to account for time and temperature dependencies using various prescriptions. The horizontal span of the lines displayed in Fig. 4 is thus meant to represent either the range of effective η/s values constrained by comparisons with experimental data or the ranges of η/s values considered in the models and yielding a good representation of the measured data. Although in both cases the reported uncertainties have been also incorporated into to the line length, the apparent relationship should not be considered as a statement of the precision achieved in the studies included in this compilation. Comparative studies of hydrodynamics and measured data arguably culminated with studies based on a Bayesian estimation of the properties of the QGP [25,[44][45][46] yielding most probable η/s values in rather good agreement with the results of this work. Estimates of η/s obtained in this work are also in quantitative agreement with QCD inspired calculations including, for instance, estimates based on nonperturbative gluon spectral functions at finite temperature in quenched QCD with the maximum entropy method [23] and calculations based on the Kubo formula in Yang-Mills theory [24]. Results of this work are also in qualitative agreement with estimates based on lattice QCD (LQCD) calculations of the QGP transport coefficients [47] and perturbative QCD calculations at almost NLO [48], as well as recent calculation based on the MUSIC framework that used a temperature dependent η/s computed with a QCD based approach [49]. In all cases in which the temperature evolution of η/s have been estimated the range reported in Fig. 4 goes up to T = 3T c which in most of the cases matches the published range. It is interesting to address the results reported in this work from an additional perspective. Although as was mentioned before the Gavin ansatz does not consider any temperature dependence for η/s the temperature reached by the medium produced in A-A collisions, presumably, will not be the same for the different centrality ranges. The values of η/s quoted in this work for each centrality range at both energies (i.e., RHIC and the LHC) thus correspond to an effective viscosity which condensates the whole system evolution for that centrality range. If as suggested above, the produced system reach different temperatures in each centrality range, there could be an implicit link between the results presented in this work and the evolution of η/s with system temperature.
Taking the estimates of η/s shown in Fig. 3 at face value, it is interesting to consider whether they might have any implications concerning the nature and properties of hot QCD matter produced in heavy-ion collisions at RHIC and LHC. First consider that estimates of the initial temperature reached in central Pb-Pb collisions at LHC suggest it is of the order of 300 MeV, i.e., 30% larger compared to that achieved at RHIC in central Au-Au collisions [63][64][65][66]. Also consider that the fireball formed in Pb-Pb collisions at the LHC have been estimated to live approximately 40% longer that those produced in Au-Au collisions at RHIC [26]. This implies that shear viscous forces have more time to operate in central Pb-Pb collisions at LHC than in Au-Au at RHIC. For systems of equal η/s and temperature, one would expect to observe a larger longitudinal broadening of the G 2 correlator in Pb-Pb but the observed broadening is in fact smaller than that seen in central Au-Au collisions. Taken at face value, this suggests that the effective shear viscosity per unit of entropy is smaller in Pb-Pb at 2.76 TeV. We should stress, however, that the extracted values of η/s reflect the complete evolution of systems formed in A-A collisions. It is consequently incorrect to associate and use a particular system temperature to evaluate the shear viscosity. Indeed, estimates should account for possible evolution of η/s with temperature explicitly or be based on an appropriate effective, time averaged, system temperature. The interpretation of the data is further complicated by the likely presence of kinematic narrowing associated to radial flow. The average transverse momentum, p T , is found to be approximately 10% larger at LHC energies compared to RHIC. This increase may in part result from faster radial flow at the TeV energy scale. It is well established that strong radial flow produces a sizable narrowing of two particle correlators, such as balance functions B [67][68][69][70], as well as generic number and transverse momentum correlators R 2 and P 2 [71], respectively. A similar narrowing is thus expected also for G 2 and has in fact been found to occur in Pb-Pb collisions: the G CD 2 correlator, in particular, exhibits a significant narrowing from peripheral to central Pb-Pb collisions reported by the ALICE collaboration [17]. While this narrowing is most easily and explicitly observed for unlikecharge particle pairs, it should also be occurring for likesign pairs contributing to the G CI 2 correlator. Narrowing effects associated with kinematic focusing might then partially counterbalance the broadening due to viscous forces, and thus effectively reduce values of η/s extracted from both the ALICE and STAR data. But given the radial flow is likely somewhat stronger at LHC, that could imply the difference seen between central Au-Au and Pb-Pb collisions is in part due to the presence of extra focusing at LHC energy.
Additional theoretical calculations [72] suggest that η/s should increase with decreasing collision energy within the RHIC energy domain in part as a result of an explicit dependence on the matter baryochemical potential µ B [73]. Studies of relative yields of produced hadrons indicate that the baryochemical potential is nearly vanishing at central rapidities in Pb-Pb collisions, with values of order µ B ∼0.7 reported by global thermal fits [74,75], while significantly larger values, µ B ∼20, were extracted based on Au-Au collisions at RHIC top energy [76]. Differences of η/s observed by STAR and ALICE collaborations might thus also result in part from this change of the baryochemical potential. While we note that the precision of the data is clearly insufficient to establish any firm conclusion on such a dependence, we stress that precise measurements of the G 2 correlator in the context of the second RHIC beam energy scan (BES-II) might in fact provide better grounds to seek evidence of this dependence. Studies of the G 2 correlator with RHIC beam energy scan data are thus indeed of high interest.
Other considerations are also of interest. Collisions of large nuclei at ultra-high energy, both at RHIC and LHC, are expected to produce very large magnetic fields and have been predicted to induce large vorticity and global polarization effects as well as finite out of plane charge separation associated with the chiral magnetic effect (CME). While the existence of the CME remains to be established, both STAR and ALICE collaborations have reported observations of global polarization of Λ -baryons [77][78][79] believed to re-sult from the presence of large vorticity in Au-Au and Pb-Pb collisions. It has been suggested that the presence of large magnetic fields might also have an impact on viscous effects [80] as they might strongly suppress momentum diffusion in the reaction plane or impart a "paramagnetic squeezing" effect capable of altering pressure gradients. Variations of the magnetic field strength and its time evolution as function of centrality and collision energy, may also influence the effective diffusivity and the viscosity of the QCD matter produced in these collisions. The magnitude of the effect is as of yet unknown but nonetheless worthy of additional investigations given the current uncertainties in η/s values do not allow to conclude about a possible difference at the two energies.
Additionally, in order to make progress on a full characterization of η/s as a function of temperature, collision energy, baryochemical potential, etc., additional and more precise measurements of G CI 2 at different collision energies are necessary; for example, as already mentioned, from the BES-II at RHIC. Furthermore, supporting theoretical studies in the framework of relativistic hydrodynamics will also be greatly beneficial. One needs, in particular, to establish the influence of the temperature and viscosity of the different stages of the collision (QGP, phase transition region, hadronic phase) on the longitudinal broadening of G CI 2 . The role of resonance decays and charge conservation must also be clarified in association with quantitative studies of the radial flow velocities imparted to the matter produced in A-A collisions. Ideally, these studies should be conducted for several system sizes and collision energies. In light of observations of Λ -baryon global polarization already mentioned, it shall also be of interest to examine whether the strong magnetic fields present at the onset of A-A collisions can persist long enough to have a quantitatively measurable impact on the shear viscosity in general, and on the longitudinal broadening of transverse momentum correlators in particular.

Conclusion
We presented an evaluation of the collision centrality dependence of the shear viscosity per unit of entropy, η/s of the Quark Gluon Plasma produced in A-A collisions at RHIC and LHC based on measurements of the G 2 correlator by the STAR and ALICE collaborations using the Gavin ansatz embodied in Eq. 1. Freeze-out times required to carry out the calculations were determined as a function of the cubic root of the charged particle multiplicity (or collision centrality) from two-pion Bose-Einstein measurements. Values of η/s obtained in Pb-Pb collisions, based on ALICE data, indicate the shear viscosity per unit of entropy is of the order of the KSS bound and essentially independent of collision centrality at LHC energy. By contrast, the STAR data are consis-tent with vanishing η/s values in peripheral collisions and values exceeding the KSS bound in more central collisions. However, given the large systematic uncertainties of these data, one cannot exclude η/s might be invariant with collision centrality. Likewise, one cannot readily exclude that values of η/s might also be invariant with beam energy.
The precision of our estimates of the dependence of η/s, particularly at RHIC, are limited by the accuracy of the STAR measurement of G 2 . Uncertainties are also largely determined by the various caveats associated with the Gavin ansatz discussed above, most particularly the choice of characteristic temperature used in the calculation. Clearly, a more detailed calculation along the lines of Ref. [81] are needed to improve on this work.
Finally, we stress that although measurements of G 2 are challenging, owing in particular to their sensitivity to p T dependent efficiency corrections, they are nonetheless possible as demonstrated by the recent ALICE measurement. Precise studies of the evolution of the G CI 2 correlator with collision centrality thus stand to become a discriminating gauge of not only the average magnitude of the shear viscosity per unit of entropy, η/s, but its temperature dependence also. As such, they might provide new and valuable inputs to multi-system Bayesian constraints methods. Ideally, this will require measurements of G 2 be completed based, for instance, on the beam energy scan at RHIC as well as for smaller collision systems at the LHC and RHIC.