Quark number susceptibilities from HTL-resummed thermodynamics

We compute analytically the diagonal quark number susceptibilities for a quark-gluon plasma at finite temperature and zero chemical potential, and compare with recent lattice results. The calculation uses the approximately self-consistent resummation of hard thermal and dense loops that we have developed previously. For temperatures between 1.5 to 5 $T_c$, our results follow the same trend as the lattice data, but exceed them in magnitude by about $5-10%$. We also compute the lowest order contribution, of order $\alpha_s^3\log(1/\alpha_s)$, to the off-diagonal susceptibility. This contribution, which is not a part of our self-consistent calculation, is numerically small, but not small enough to be compatible with a recent lattice simulation.