q=1/(1-exp(-E/kT))