High order quark number susceptibilities in hot QCD from lattice EQCD

Building on the experience of [1], 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 formulae for the operators up to the sixth derivative, and obtain continuum-extrapolated results from lattice EQCD at zero and finite $\mu$ and at six different pairs of temperature T and number of massless quark flavors $n_f$. Our data is benchmarked against full-QCD lattice and perturbative results, allowing to judge the quality of the perturbative series expansion in EQCD and the dimensional reduction procedure as a whole.