Skewness, kurtosis and the 5th and 6th order cumulants of net baryon-number distributions from lattice QCD confront high-statistics STAR data

We present new results on up to $6^{th}$ order cumulants of net baryon-number fluctuations at small values of the baryon chemical potential, $mu_B$, obtained in lattice QCD calculations with physical values of light and strange quark masses. Representation of the Taylor expansions of higher order cumulants in terms of the ratio of the two lowest order cumulants, $M_B/sigma_B^2=chi_1^B(T,mu_B)/chi_2^B(T,mu_B)$, allows for a parameter free comparison with data on net proton-number cumulants obtained by the STAR Collaboration in the Beam Energy Scan at RHIC. We show that recent high statistics data on skewness and kurtosis ratios of net proton-number distributions, obtained at beam energy $sqrt{s_{_{NN}}}=54.4$ GeV, agree well with lattice QCD results on cumulants of net baryon-number fluctuations close to the pseudo-critical temperature, $T_{pc}(mu_B)$, for the chiral transition in QCD. We also present first results from a next-to-leading order expansion of $5^{th}$ and $6^{th}$ order cumulants on the line of pseudo-critical temperatures.