Extended series expansions for the high temperature zero-field susceptibility of the Ising model are given in powers of the usual high temperature counting variable v=tanh K; for the triangular lattice to v 16 , for the square lattice to v 21 and for the honeycomb lattice to v 32 , inclusive. The asymptotic behaviour of the ferromagnetic and antiferromagnetic susceptibility is studied. It is concluded that the ferromagnetic singularity is not exactly factorizible. The antiferromagnetic susceptibility of the square and honeycomb lattices has a singularity of the same type as the energy at the antiferromagnetic critical temperature.
For Pt.1 see ibid., vol.8, no.8, p.1236 of 1975. Amplitudes for a number of three-dimensional Ising model series have been estimated from recently extended series expansions. The amplitudes have then been used to re-estimate the field and temperature scaling parameters that are fundamental to the Generalized Law of Corresponding States. These new scaling parameters have been used (i) to estimate the amplitudes of the next most singular term in the high-temperature susceptibility for the four common three-dimensional lattices, and (ii) to show that an earlier conjecture, that the field scaling parameter n x is the same for the Ising and spherical models, is unlikely to hold.
We have extended the series for the zero-field susceptibility of the spin-
s Ising model through tenth order in the reduced inverse temperature
K on the square, triangular, simple-cubic, body-centered-cubic, and face-centered-cubic lattices. The series coefficients
hn(s) are expressed as simple polynomials in
X=s(s+1). Using extended methods of analysis we have estimated the nature of the leading singularities on the face-centered lattice and conclude with good confidence that the susceptibility exponent ${\ensuremath{\gamma}}_{1}$ equals 1.25, independent of
s. The exponent of the leading correction term is estimated to be ${\ensuremath{\gamma}}_{2}\ensuremath{\simeq}0.75\ifmmode\pm\else\textpm\fi{}0.08$ in good agreement with renormalization-group theory. For
s=21 only, the amplitude of the confluent correction apparently vanishes. We have also studied the leading singularities on the triangular lattice and conclude that ${\ensuremath{\gamma}}_{1}=\frac{7}{4}$, independent of
s, with, however, much stronger corrections than in three dimensions. These results provide a very strong corroboration of the universality hypothesis.
We investigate the complex-temperature singularities of the susceptibility of the 2D Ising model on a square lattice. From an analysis of low-temperature series expansions, we find evidence that, as one approaches the point u=u s =-1 (where u=e '4K ) from within the complex extensions of the FM or AFM phases, the susceptibility has a divergent singularity of the form X approximately A s '(1+u)(- gamma s ') with exponent gamma s '=3/2. The critical amplitude A s ' is calculated. Other critical exponents are found to be alpha s '= sigma s =0 and beta s = 1/4 , so that the scaling relation alpha s '+2 beta s + gamma s ' is satisfied. However, using exact results for beta s on the square, triangular, and honeycomb lattices, we show that universality is violated at this singularity: beta s is lattice-dependent. Finally, from an analysis of spin-spin correlation functions, we demonstrate that the correlation length and hence susceptibility are finite as one approaches the point u=-1 from within the symmetric phase. This is confirmed by an explicit study of high-temperature series expansions.
The high-temperature susceptibility of the spin-
S Ising model on the triangular net is analyzed as a function of
S using tenth-order high-temperature series expansions. In particular the spin dependence of the critical temperature and the amplitude of the dominant singularity are obtained by assuming a universal value $\ensuremath{\gamma}=\frac{7}{4}$ for the critical exponent of the susceptibility. The results, for the critical temperature particularly, provide a valuable guidepost in judging the convergence of recent nonlinear renormalization-group calculations for these models.