▎ 摘　　要

In this work, we present a semianalytical expression for the temperature dependence of a spin-resolved dynamical density-density response function of massless Dirac fermions within the random phase approximation. This result is crucial in order to describe thermodynamic properties of the interacting systems. In particular, we use it to make quantitative predictions for the paramagnetic spin susceptibility of doped graphene sheets. We find that, at low temperatures, the spin susceptibility behaves like T-2, which is completely different from the temperature dependence of the magnetic susceptibility in undoped graphene sheets.