We review the applications of the integral over anticommuting Grassmann variables (nonquantum fermionic fields) to the analytic solutions and the field-theoretical formulations for the 2D Ising models. The 2D Ising model partition function
Q is presentable as the fermionic Gaussian integral. The use of the spin-polynomial interpretation of the 2D Ising problem is stressed, in particular. Starting with the spin-polynomial interpretation of the local Boltzmann weights, the Gaussian integral for
Q appears in the universal form for a variety of lattices, including the standard rectangular, triangular, and hexagonal lattices, and with the minimal number of fermionic variables (two per site). The analytic solutions for the correspondent 2D Ising models then follow by passing to the momentum space on a lattice. The symmetries and the question on the location of critical point have an interesting interpretation within this spin-polynomial formulation of the problem. From the exact lattice theory we then pass to the continuum-limit field-theoretical interpretation of the 2D Ising models. The continuum theory captures all relevant features of the original models near
Tc. The continuum limit corresponds to the low-momentum sector of the exact theory responsible for the critical-point singularities and the large-distance behaviour of correlations. The resulting field theory is the massive two-component Majorana theory, with mass vanishing at
Tc. By doubling of fermions in the Majorana representation, we obtain as well the 2D Dirac field theory of charged fermions for 2D Ising models. The differences between particular 2D Ising lattices are merely adsorbed, in the field-theoretical formulation, in the definition of the effective mass.
We present and solve the replica symmetric equations in the context of the replica cluster variational method for the 2D random bond Ising model (including the 2D Edwards–Anderson spin-glass model). First, we solve a linearized version of these equations to obtain the phase diagrams of the model on the square and triangular lattices. In both cases, the spin-glass transition temperatures and the multicritical point estimations improve largely over the Bethe predictions. Moreover, we show that this phase diagram is consistent with the behavior of inference algorithms on single instances of the problem. Finally, we present a method to consistently find approximate solutions to the equations in the glassy phase. The method is applied to the triangular lattice down to T = 0, also in the presence of an external field.