The transport properties of the strongly coupled quark-gluon plasma created in ultrarelativistic heavy-ion collisions are extracted by Bayesian parameter estimate methods with the latest collision beam energy data from the CERN Large Hadron Collider. This Bayesian analysis includes sophisticated flow harmonic observables for the first time. We found that the temperature dependence of specific shear viscosity appears weaker than in the previous studies. The results prefer a lower value of specific bulk viscosity and a higher switching temperature to reproduce additional observables. However, the improved statistical uncertainties both on the experimental data and hydrodynamic calculations with additional observables do not help to reduce the final credibility ranges much, indicating a need for improving the dynamical collision model before the hydrodynamic takes place. In addition, the sensitivities of experimental observables to the parameters in hydrodynamic model calculations are quantified. It is found that the analysis benefits most from the symmetric cumulants and nonlinear flow modes, which mostly reflect nonlinear hydrodynamic responses, in constraining the temperature dependence of the specific shear and bulk viscosities in addition to the previously used flow coefficients.Peer reviewe
