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$.
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.
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 method of developing exact power-series expansions for the partition function ZN and related thermodynamic functions for the Ising model valid below the critical point is generalized to include exchange interactions between first-, second-, and third-neighbor pairs. Expansions of the spontaneous magnetization M0(T) and zero field susceptibility χ0(T) are derived through to sixth order of perturbation for the s.q. lattice, and through to fifth order of perturbation for the Δ′r, b.c.c., s.c., and f.c.c. lattices, when interactions J1SizSjz and J2SkzSlz are present between first- and second-neighbor spins, respectively (second-neighbor model). These expansions have also been obtained for the case where interactions of equal magnitude (J1 = J2 = J3) are present between first-, second-, and third-neighbor pairs (third-equivalent-neighbor model); here expansions through to fifth order of perturbation are obtained for the s.q., Δ′r, b.c.c., and s.c. lattices and through to fourth order for the f.c.c. lattice. The Padé approximant·procedure is employed to discuss the effects of an extended but finite range of interaction on the behavior of M0(T) and χ0(T) for T→Tc− as characterized by the critical exponents β and γ′, respectively. For the second-equivalent-neighbor model lattices, it is found that 0.122 ≤ β ≤ 0.134 in two dimensions, and that 0.308 ≤ β ≤ 0.328 in three dimensions; from which it is concluded that β remains unchanged from its value in the nearest-neighbor model. The corresponding limits for γ′ in three dimensions are 1.18 ≤ γ′ ≤ 1.28; from this and the results for the b.c.c. lattice in particular, it is concluded that γ′ is probably 54 and hence the transition in χ0 is symmetrical about Tc (γ′ = γ). A repetition of this analysis for the third-equivalent-neighbor model three-dimensional lattices shows a marked shift in the estimated range of β and γ′; the results are 0.345 ≤ β ≤ 0.365, and 1.01 ≤ γ′ ≤ 1.14. In each of the above cases, the corresponding high-temperature (T > Tc) expansions of χ0(T) obtained previously have been analyzed to yield estimates of the critical exponent γ. The over-all results and in particular the estimates of γ for the s.q. and b.c.c. lattices suggest that this index is unaffected by extending the range of interaction, and that if γ is a rational fraction then it is the same fraction for the n.n. model and second- and third-equivalent neighbor models. Finally the high-temperature expansions of ZN in zero field, and of χ0(T) for the second-neighbor model are used to examine the dependence of the critical temperature Tc, the critical energy (E∞ − Ec)/kTc, and the critical entropy (S∞ − Sc)/k on the relative strengths of J1 and J2 for values of J2/J1 in the range 0 to 1. It is found that the variation of the critical point is well represented by Tc(α)=Tc(0)[1+m1α],where α = J2/J1 and lies in the range 0 ≤ α ≤ 1; and Tc(0) is the critical temperature of the nearest-neighbor model. The values of m1 are 0.61, 2.47, 0.84, 1.45, and 1.35 for the f.c.c., s.c., b.c.c., s.q., and Δ′r lattices, respectively. All these calculations are compared with the corresponding results for the Heisenberg model.
For pt. I see ibid., vol. 5, no. 5, 624. 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 simple cubic lattice to v 17 , the body-centred cubic lattice to v 15 , and the face-centred cubic lattice to v 12 , inclusive. The coefficients are analysed by the ratio method and it is estimated that, subject to the implicit assumptions thereof, the critical temperatures correspond to v f =0.21813+or-0.00001 (SC), v f =0.15612+or-0.00003 (BCC) and v f =0.10174+or-0.00001 (FCC). The symptotic behaviour of the coefficients is studied in detail; for the loose-packed lattices the decay of the even-odd oscillation in the ratios is found to be consistent with the assumption that the critical index of the high temperature specific heat is very close to 1/8.