Grassmann Algebra and Fermions at Finite Temperature
  For any d-dimensional self-interacting fermionic model, all coefficients in the high-temperature expansion of its grand canonical partition function can be put in terms of multivariable Grassmann integrals. A new approach to calculate such coefficients, based on direct exploitation of the grassmannian nature of fermionic operators, is presented. We apply the method to the soluble Hatsugai-Kohmoto model, reobtaining well-known results. 
I.C. Charret, E.V. Corr\^ea Silva, S.M. de Souza, O. Rojas Santos and   M.T. Thomaz