By combinatorial reasons, for 2D toric Ising model of fixed size I need lowtemperature expansion for complex temperature, magnetic field, etc. Is there are some references? More precisely, energy of spin configuration $x_i=+1,1$ given by $H=\sum_{edges}(N+a)(x_i+x_j)Nx_ix_j$ for some $a\in i\mathbb{R}$, where sum going over all edges in latticegraph on torus surface, and all what i need is asymptotics in first term of partition function $\sum_{states} e^{H}$ when $N\to +\infty$
I don't think Onsager's solution will help you in the presence of a magnetic field. In the low temperature regime there are rigorous expansions for pretty much all the quantities of interest. These go under the names of low temperature cluster expansions, contour expansions, PirogovSinai theory (the most general theory for this kind of things). You can find presentations of these tools especially in the archetypal example of the Ising model in just about any book on Gibbs measures: e.g. the books by Malyshev and Minlos, Preston, Georgii, Sinai, etc. A good place to start is the lecture notes by Velenik. If you do not read French then look up the references in these lecture notes. 


This problem is known in literature due to the Onsager's solution of the 2d Ising model. You can already find this on an old paper by Fisher and Ferdinand. Some more recent papers, available also from arxiv are the following: http://arxiv.org/abs/condmat/0009054 http://arxiv.org/abs/condmat/0110287 A paper just available with subscription is Partition function of a finite Ising model on a torus T Morita 1986 J. Phys. A: Math. Gen. 19 L1191. 

