Power-series expansions of the spin-pair correlation functions of the square, sc, and bcc Ising lattices have been obtained using the semi-invariant approach for general field
H and temperature
T. In three dimensions P\'ade-approximant analyses indicate $2{v}^{\ensuremath{'}}=1.285\ifmmode\pm\else\textpm\fi{}0.020$, in agreement with scaling, but $2{v}^{c}\ensuremath{\simeq}0.835\ifmmode\pm\else\textpm\fi{}0.02$ which is some (1-3)% above the scaling prediction ${v}^{c}=\frac{v}{\ensuremath{\beta}\ensuremath{\delta}}\ensuremath{\simeq}0.823$. However, ratio techniques reveal that this discrepancy can be attributed to significant confluent critical singularities. Cubic and quintic parametric representations of the critical equation of state, and corresponding expressions for the correlation length, ${\ensuremath{\xi}}_{1}(H, T)$, are developed, which are considerable improvements over the linear model. The universality and spherical symmetry of the critical scattering intensity $\stackrel{^}{\ensuremath{\chi}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}, H, T)$ is confirmed to within (1-2)% by estimating suitable invariant combinations of amplitudes. The deviations from Ornstein-Zernike theory for general
H and
T are found to be considerably greater than for
H=0,
T>Tc. Complete parametric scaling representations of $\stackrel{^}{\ensuremath{\chi}}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}, H, T)$ are developed for three dimensions; corresponding scaling approximants are constructed in two dimensions but only for
T=Tc and for $T\ensuremath{\gtrless}{T}_{c}$ with
H=0.
The Ising-model correlation function $C({R}_{12})=〈{\ensuremath{\mu}}_{1}{\ensuremath{\mu}}_{2}〉$ is studied in terms of a novel
N-fold integral representation. This formula stems from a procedure proposed by Montroll and Berlin. The integral is estimated by maximizing the integrand, an approximation related to the spherical-model assumptions. The correlation function is not of the Ornstein-Zernike type, just above the critical point, but rather $C(R)\ensuremath{\propto}{R}^{\ensuremath{-}1\ensuremath{-}\ensuremath{\eta}}$ for $R\ensuremath{\ll}\frac{1}{\ensuremath{\kappa}}$, and $C(R)\ensuremath{\propto}{\ensuremath{\kappa}}^{\ensuremath{\eta}}{R}^{\ensuremath{-}1}\mathrm{exp}(\ensuremath{-}\ensuremath{\kappa}R)$ for $R\ensuremath{\gg}\frac{1}{\ensuremath{\kappa}}$. The correlation length $\frac{1}{\ensuremath{\kappa}}$ becomes infinite at the critical point. The calculated value $\ensuremath{\eta}=0.646$ is too large, reflecting the omission of important terms in the evaluation of the integral. The unusual mechanism inducing the nonclassical behavior is carefully examined.
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.
A version of the Ising model is developed in which the spin variables can be treated accurately in the continuum approximation. The perturbation series, both above and below the critical temperature
Tc, is examined, and it is shown that there is a shift of
Tc from its mean-field value proportional to ${q}^{\ensuremath{-}2}\mathrm{ln}q$, as well as the well-known shift proportional to ${q}^{\ensuremath{-}1}$; here
q is the number of mutually interacting particles. It is shown, using renormalization theory, that there is a perturbation series in ${q}^{\ensuremath{-}1}{|T\ensuremath{-}{T}_{c}|}^{\ensuremath{-}\frac{1}{2}}$ for which all terms are finite in the limit $q\ensuremath{\rightarrow}\ensuremath{\infty}$, if the shift of
Tc is put in correctly. For the two-dimensional model, the shift is shown to be proportional to ${q}^{\ensuremath{-}1}\mathrm{ln}q$. Conditions are derived for a finite system to display critical behavior characteristic of three, two, one, or zero dimensions. It is shown how similar results can be obtained for a model similar to the Heisenberg model and for the standard Ising and Heisenberg models with interactions extending over many neighbors. A comparison is made between previously calculated numerical results for
Tc and the asymptotic forms derived here.
We consider a one-dimensional model with infinite-range interaction, a two-dimensional model, and a three-dimensional model, whose free energies can be expressed in terms of the largest eigenvalue of an integral equation. High- and low-temperature expansions in powers of the reciprocal of the range of the exponential part of the interaction, with the classical Curie-Weiss theory as leading term, are developed and studied in the critical region. We find that to leading order in the critical region the resummed high- and low-temperature expansions are analytic at the classical critical point but are nonanalytic at a displaced critical point. The modified singularities, which are no longer of Curie-Weiss type, give critical exponents which are identical with those obtained by Brout and others, and are almost surely not the true exponents. The technique, however, suggests a possible general method of successive approximation to true critical behavior.
The behavior of the spherical model of a ferromagnet with an interaction energy between the magnetic spins which varies with distance as $\frac{1}{{r}^{d+\ensuremath{\sigma}}}$ (where
d is the dimensionality of the lattice and $\ensuremath{\sigma}>0$) is analyzed. It is shown that the model exhibits a ferromagnetic transition in one and two dimensions, providing $0<\ensuremath{\sigma}<d$. (The usual spherical model with nearest-neighbor interactions does not have a transition in one and two dimensions.) The critical-point behavior is investigated. It is found that the singularities in the specific heat and susceptibility are dependent on $\ensuremath{\sigma}$ and
d, but the behavior of the magnetization is independent of $\ensuremath{\sigma}$ and
d. In three dimensions the susceptibility diverges as ${(T\ensuremath{-}{T}_{c})}^{\ensuremath{-}\ensuremath{\gamma}}$, where $\ensuremath{\gamma}=1$ for $0<\ensuremath{\sigma}<\frac{3}{2}$, $\ensuremath{\gamma}=\frac{\ensuremath{\sigma}}{(3\ensuremath{-}\ensuremath{\sigma})}$ for $\frac{3}{2}<\ensuremath{\sigma}<2$ and $\ensuremath{\gamma}=2$ for $\ensuremath{\sigma}>2$. The asymptotic form of the spin-spin correlation function $\ensuremath{\Gamma}(\mathrm{r})$ is studied in the neighborhood of the critical temperature
Tc. At
T=Tc, $\ensuremath{\Gamma}(\mathrm{r})$ decays for large
r as $\frac{1}{{r}^{d\ensuremath{-}\ensuremath{\sigma}}}$. Several two-dimensional models with long-range interactions falling off as
r21 in certain directions only are also investigated.