High order quark number susceptibilities in hot QCD from lattice electrostatic QCD

Building on the experience of [A. Hietanen et al., Phys. Rev. D 79, 045018 (2009)], we develop a formalism to construct operators for higher derivatives of the pressure in hot QCD with respect to the quark chemical potential mu. We provide formulas for the operators up to the sixth derivative, and obtain continuum-extrapolated results from lattice electrostatic QCD (EQCD) at zero and finite mu and at six different pairs of temperature T and number of massless quark flavors nf. Our data is benchmarked against full-QCD lattice and perturbative results, allowing us to judge the quality of the perturbative series expansion in EQCD and the dimensional reduction procedure as a whole.Peer reviewe