Questions concerning the universality between different two-dimensional Ising-like systems have recently been raised, with evidence for the existence of several different subclasses of Ising-like systems being presented. Some new facts concerning similarity between border model and low-temperature spin-1 Ising model series analyses are noted. For the spin-1 case an alternative determination of the critical temperature determines unequivocally that this model is in the Ising universality class and it is suggested that a similar independent result could lead to a similar conclusion in the border model case. The possibility of an analytic correction to scaling in the spin-1 model is proposed.
We have extended through ${\ensuremath{\beta}}^{23}$ the high-temperature expansion of the second field derivative of the susceptibility for Ising models of general spin, with nearest-neighbor interactions, on the simple cubic and the body-centered cubic lattices. Moreover, the expansions for the nearest-neighbor correlation function, the susceptibility, and the second correlation moment have been extended up to ${\ensuremath{\beta}}^{25}.$ Taking advantage of these new data, we can improve the accuracy of direct estimates of critical exponents and of hyperuniversal combinations of critical amplitudes such as the renormalized four-point coupling
gr or the quantity usually denoted by ${R}_{\ensuremath{\xi}}^{+}.$ In particular, we obtain $\ensuremath{\gamma}=1.2371(1),$ $\ensuremath{\nu}=0.6299(2),$ ${\ensuremath{\gamma}}_{4}=4.3647(20),$
gr=1.404(3), and ${R}_{\ensuremath{\xi}}^{+}=0.2668(5).$ We have used a variety of series extrapolation procedures and, in some of the analyses, we have assumed that the leading correction-to-scaling exponent $\ensuremath{\theta}$ is universal and roughly known. We have also verified, to high precision, the validity of the hyperscaling relation and of the universality property both with regard to the lattice structure and to the value of the spin.
High-temperature series expansions are used to examine the dependence of critical-point exponents upon lattice anisotropy (different interaction strengths in different directions in the lattice). The two-spin correlation function ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ is calculated to tenth order in (
kBT1) for the Ising Hamiltonian ${\mathcal{H}}_{l \mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}\stackrel{\mathrm{xy}}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}\stackrel{z}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}$ for a wide range of anisotropy parameters $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}$ and for both the sc and fcc lattices; here the first summation is over all pairs of nearest-neighbor sites whose relative displacement vector ${\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{\mathrm{ij}}$ has no
z component, while the second summation is over all other pairs of nearest-neighbor sites. Hence for
R=0, both the sc and fcc lattices reduce to the two-dimensional square lattice, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc becomes a one-dimensional linear chain and the fcc becomes two noninteracting three-dimensional bcc lattices. The series for ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ are then used to obtain series of corresponding lengths for the specific heat, susceptibility, and second moment. Analysis of these series yields results consistent with the universality hypothesis of critical-point exponents. Specifically, it is found that when lattice anisotrophy is introduced, the critical-point exponents studied (the susceptibility exponent $\ensuremath{\gamma}$ and the correlation length exponent $\ensuremath{\nu}$) do not appear to change from their values for an isotropic lattice. The problem of next-nearest-neighbor interactions is treated using similar methods in Paper II of this series (and briefly discussed in this paper).
Using a renormalized linked-cluster-expansion method, we have extended to order ${\ensuremath{\beta}}^{23}$ the high-temperature series for the susceptibility $\ensuremath{\chi}$ and the second-moment correlation length $\ensuremath{\xi}$ of the spin-1/2 Ising models on the sc and the bcc lattices. A study of these expansions yields updated direct estimates of universal parameters, such as exponents and amplitude ratios, which characterize the critical behavior of $\ensuremath{\chi}$ and $\ensuremath{\xi}.$ Our best estimates for the inverse critical temperatures are ${\ensuremath{\beta}}_{c}^{\mathrm{sc}}=0.221654(1)$ and ${\ensuremath{\beta}}_{c}^{\mathrm{bcc}}=0.1573725(6).$ For the susceptibility exponent we get $\ensuremath{\gamma}=1.2375(6)$ and for the correlation length exponent $\ensuremath{\nu}=0.6302(4).$ The ratio of the critical amplitudes of $\ensuremath{\chi}$ above and below the critical temperature is estimated to be ${C}_{+}{/C}_{\ensuremath{-}}=4.762(8).$ The analogous ratio for $\ensuremath{\xi}$ is estimated to be ${f}_{+}{/f}_{\ensuremath{-}}=1.963(8).$ For the correction-to-scaling amplitude ratio we obtain ${a}_{\ensuremath{\xi}}^{+}{/a}_{\ensuremath{\chi}}^{+}=0.87(6).$
The high-temperature expansion coefficients of the ordinary and the higher susceptibilities of the spin-1/2 nearest-neighbor Ising model are calculated exactly up to the 20th order for the general d-dimensional (hyper)simple-cubical lattices. These series are analyzed to study the dependence of critical parameters on the lattice dimensionality. Using the general d expression of the ordinary susceptibility, we have more than doubled the length of the existing series expansion of the critical temperature in powers of 1/d.
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.