The specific heat can be calculated from the mean expansion order, and the mean of the squared expansion order. It is proportional to <n^2>-<n>^2-<n>, when n is the order. You will need to look up the prefactors in the SSE references. I don't them with me right now. Matthias On 16 Jan 2014, at 05:10, Tadeusz Wasiutyński <tadeusz.wasiutynski@ifj.edu.pl> wrote:

Is it possible to calculate specific heat within spin model by dirloop_sse? From observables list it is not. It looks that <E> is observable but <E2> is not. Am I right? Tadeusz