High-temperature series expansions are used to examine the dependence of critical-point exponents upon the presence of second-neighbor interactions. We consider the Hamiltonian ${\mathcal{H}}_{\mathrm{nnn}}=\ensuremath{-}{J}_{1}\ensuremath{\Sigma}\stackrel{\mathrm{nn}}{〈\mathrm{ij}〉}{\stackrel{\ensuremath{\rightarrow}}{\mathrm{S}}}_{i}^{(D)}\ifmmode\cdot\else\textperiodcentered\fi{}{\stackrel{\ensuremath{\rightarrow}}{\mathrm{S}}}_{j}^{(D)}\ensuremath{-}{J}_{2}\ensuremath{\Sigma}\stackrel{\mathrm{nnn}}{〈\mathrm{ij}〉}{\stackrel{\ensuremath{\rightarrow}}{\mathrm{S}}}_{i}^{(D)}\ifmmode\cdot\else\textperiodcentered\fi{}{\stackrel{\ensuremath{\rightarrow}}{\mathrm{S}}}_{j}^{(D)},$ where the first and second sums are over pairs of nearest-neighbor (nn) and next-nearest-neighbor (nnn) sites, and where the spins ${\stackrel{\ensuremath{\rightarrow}}{\mathrm{S}}}^{(D)}$ are
D-dimensional unit vectors. The two-spin correlation function, ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$, is calculated to tenth, ninth, and eighth order in
kBT1 for the Ising (
D=1), classical-planar (
D=2), and classical-Heisenberg (
D=3) models, respectively, for various values of the parameter ${R}^{\ensuremath{'}}\ensuremath{\equiv}\frac{{J}_{2}}{{J}_{1}}$ and for various cubic lattices (fcc, bcc, and sample cubic). These represent the first series expansions of the spin correlation function for nnn interactions. From ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ we obtain series for the specific heat, susceptibility, and second moment. Analysis of these series and detailed comparisons with the exactly soluble spherical model ($D=\ensuremath{\infty}$) lead us to conclude that the exponents $\ensuremath{\gamma}$ (susceptibility) and $\ensuremath{\nu}$ (correlation length) may be independent of ${R}^{\ensuremath{'}}$; this suggestion is consistent with the universality hypothesis.
Series for the reduced susceptibility $\overline{\ensuremath{\chi}}$, the reduced specific heat
CH, and second moment ${\ensuremath{\mu}}_{2}$ of the static correlation function for the three-dimensional
S=21 Ising model on both the simple cubic (sc) and face-centered cubic (fcc) lattices with different coupling strengths in different lattice directions have been analyzed to determine the crossover exponent $\ensuremath{\phi}$ describing the behavior of the critical temperature as a function of the anisotropy parameter
R in the Hamiltonian $\mathcal{H}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j}\ensuremath{\equiv}\ensuremath{-}{J}_{\mathrm{xy}}(\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}+R\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j})$. Here
si=+1, the first sum is over all nearest-neighbor pairs in the
xy plane, and the second sum is over all pairs coupled in the
z direction. The constant gap exponent we obtain for successive derivatives of $\overline{\ensuremath{\chi}}$ and
CH with respect to
R confirms the exponent predictions of scaling in the parameter
R for thermodynamic functions, while the results of the ${\ensuremath{\mu}}_{2}$ series confirm the exponent predictions of scaling with respect to
R for the two-spin correlation function. Our results agree with the predictions for $\ensuremath{\phi}$ of Abe and Suzuki, and also with rigorous relations satisfied by the exponents describing the derivatives of the various functions. Our results do not agree with previously published results, which are based on an analysis of only the susceptibility on only the sc lattice.
Of particular current interest is the critical behavior of functions on crossing over from one lattice dimensionality to another. To this end, we report high-temperature series for an Ising model with lattice anisotropy-i.e., with different exchange constants for different lattice directions. The Hamiltonian is ${\mathcal{H}}_{l,\mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j}\ensuremath{\equiv}\ensuremath{-}{J}_{\mathrm{xy}}(\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}+R\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j})$ where ${s}_{i}=\ifmmode\pm\else\textpm\fi{}1$, the first summation is over all nearest-neighbor pairs in the
xy plane, and the second sum is over pairs coupled in the
z direction. The susceptibility, second moment, and specific-heat series are explicitly presented for arbitrary
Jxy and
Jz for the simple cubic (sc) and face-centered cubic (fcc) lattices to tenth order in inverse temperature. The general-
R series are essential if one wishes to study the Riedel-Wegner crossover exponent appropriate to changing lattice dimensionality, since for $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}=0$, both the sc and fcc lattices reduce to two-dimensional square lattices, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc reduces to noninteracting linear chains.
The theory of the correlations and critical scattering of two- and three-dimensional nearest-neighbor Ising models is discussed critically. A distinction is drawn between $\ensuremath{\kappa}(T)$, the true inverse range of exponential decay of the correlations, and ${\ensuremath{\kappa}}_{1}(T)$, the effective range determined from the low-angle scattering intensity. Ten to eleven terms of appropriate high-temperature series exapnsions for $\ensuremath{\kappa}$ and ${\ensuremath{\kappa}}_{1}$ are determined for the square and simple cubic lattices, and shorter series are given for the triangular, fcc, and bcc lattices. For the former lattices, the complete correlation expansions are obtained to the same order. It is shown that $\ensuremath{\kappa}$ and ${\ensuremath{\kappa}}_{1}$ vary as ${(T\ensuremath{-}{T}_{0})}^{\ensuremath{\nu}}$ when $T\ensuremath{\rightarrow}{T}_{c}$, with $\ensuremath{\nu}=1$ for dimensionality
d=2, but $\ensuremath{\nu}=0.6430\ifmmode\pm\else\textpm\fi{}0.0025\ensuremath{\simeq}\frac{9}{14}$ for
d=3. The asymptotic decay of correlation at
T=Tc is found to be $\frac{1}{{r}^{d\ensuremath{-}2+\ensuremath{\eta}}}$, where $\ensuremath{\eta}$ is related to the exponent $\ensuremath{\gamma}$ of the divergence of the susceptibility by $(2\ensuremath{-}\ensuremath{\eta})\ensuremath{\nu}=\ensuremath{\gamma}$, Numerical values are $\ensuremath{\eta}=\frac{1}{4}$ for
d=2 and $\ensuremath{\eta}=0.056\ifmmode\pm\else\textpm\fi{}0.008\ensuremath{\simeq}\frac{1}{18} \mathrm{for} d=3$. The relative scattering intensity $\stackrel{^}{\ensuremath{\chi}}$ as a function of wave number k is given to high accuracy for all
T> Tc by $\stackrel{^}{\ensuremath{\chi}}(\mathbf{k},T)\ensuremath{\simeq}\left(\frac{a}{{r}_{1}}\right)\frac{^{2\ensuremath{-}\ensuremath{\eta}}[{({\ensuremath{\kappa}}_{1}a)}^{2}+{\ensuremath{\varphi}}^{2}{a}^{2}{K}^{2}(\mathbf{k})]^{\frac{\ensuremath{\eta}}{2}}}{[{({\ensuremath{\kappa}}_{1}a)}^{2}+\ensuremath{\psi}{a}^{2}{K}^{2}(\mathbf{k})]},$ where (i)
a is the lattice spacing, (ii) ${a}^{2}{K}^{2}=2d[1\ensuremath{-}{q}^{\ensuremath{-}1}\ensuremath{\Sigma}\mathrm{exp}(i\mathbf{k}\ifmmode\cdot\else\textperiodcentered\fi{}\mathbf{r})]\ensuremath{\simeq}{(\mathrm{ka})}^{2}$, the sum being over the
q nearest-neighbor lattice sites, (iii)
r1(T) is a slowly-varying decreasing function near
Tc, (iv) $\ensuremath{\psi}=1+\frac{1}{2}\ensuremath{\eta}{\ensuremath{\varphi}}^{2}$, and $(\mathrm{v}) \ensuremath{\varphi}(T)$ is slowly varying with a magnitude at
Tc of 0.03 for
d=2 and of 0.06 to 0.09 for
d=3. Explicit formulas are given for ${\ensuremath{\kappa}}_{1}$,
r1, and $\ensuremath{\varphi}$ as functions of
T. The correlations and the scattering are isotropic near
Tc. The critical scattering isotherm is curved for low
k according to ${\stackrel{^}{\ensuremath{\chi}}}^{\ensuremath{-}1}\ensuremath{\sim}{k}^{2\ensuremath{-}\ensuremath{\eta}}$ and it intersects the isotherms for
T>Tc. Correspondingly, $\stackrel{^}{\ensuremath{\chi}}(\mathbf{k},T)$ exhibits a maximum for fixed k, at a temperature above
Tc; for
d=2 the maxima are very well marked, but for
d=3 they are smaller and occur closer to
Tc. The theory is compared favorably with recent neutron-scattering experiments on pure beta-brass.
The spin-spin correlation functions and the critical-scattering intensity for Heisenberg models of general spin,
S=21 to $\ensuremath{\infty}$, on the sc, bcc, and fcc lattices are studied on the basis of high-temperature series expansions along the lines developed in Paper I [M. E. Fisher and R. J. Burford, Phys. Rev. 156, 583 (1967)]. Subject to increased uncertainties for low spin, it is concluded that the exponents $\ensuremath{\gamma}={1.375}_{\ensuremath{-}0.01}^{+0.02}$, $2\ensuremath{\nu}={1.405}_{\ensuremath{-}0.01}^{+0.02}$, and $\ensuremath{\eta}=0.043\ifmmode\pm\else\textpm\fi{}0.014$ describe all lattices and all spin. Explicit formulas are presented for the susceptibility/zero-angle scattering ${\ensuremath{\chi}}_{0}(T)$, for the inverse correlation length ${\ensuremath{\kappa}}_{1}(T)$, for the effective interaction range
r1(T), and using the Fisher-Burford approximant, for the total scattering $\stackrel{^}{\ensuremath{\chi}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}},T)$. The shape parameter ${\ensuremath{\varphi}}_{c}$ attains the "universal" value ${\ensuremath{\varphi}}_{c}\ensuremath{\simeq}0$. 11 for large spin but shows signs of spin dependence (and lattice dependence) for low spin. At fixed k the scattering is predicted to display a maximum above
Tc determined by $\frac{{\ensuremath{\kappa}}_{1}({T}_{max})}{k}\ensuremath{\simeq}0.10 (\mathrm{for} S\ensuremath{\gtrsim}2)$ to 0.15. A detailed study is made of the structure dependence of the critical-point correlations ${〈{S}_{0}^{z}{S}_{\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{z}〉}_{c}$ for various models. This leads to the revised, universal estimate ${\ensuremath{\varphi}}_{c}\ensuremath{\simeq}0$. 15 for all three cubic lattice, spin-\textonehalf{} Ising models. The results are compared d briefly with various experiments which support $\ensuremath{\eta}\ensuremath{\gtrsim}0.05$.
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.