  In this work, we present a semi-analytical 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.
