The low-density linked-cluster expansion of the ferromagnetic Ising model $F(H,T)=\ensuremath{-}mH\ensuremath{-}\frac{1}{2}qJ\ensuremath{-}{k}_{B}T\ensuremath{\Sigma}\stackrel{\ensuremath{\infty}}{l=1}{L}_{l}(u){\ensuremath{\mu}}^{l}$ ($u={e}^{\frac{\ensuremath{-}4J}{{k}_{B}T}},\ensuremath{\mu}={e}^{\frac{\ensuremath{-}2\mathrm{mH}}{{k}_{B}}T}$) is numerically analyzed. We find that all the polynomials
Ll are positive between 0 and ${u}_{c}={e}^{\frac{\ensuremath{-}4J}{{k}_{B}{T}_{c}}}$,
Tc being the critical temperature, and that
uc can be calculated with good accuracy as the limit point of the zeros of these polynomials
Ll. Given
uc, we compute
Ll(uc) and find, numerically, their asymptotic behavior to be ${c}_{0}{l}^{\ensuremath{-}s}$. It is then shown that the singular part of the free energy per spin, as $H\ensuremath{\rightarrow}0$, is given by $F(H,{T}_{c})\ensuremath{\sim}{(\ensuremath{-}\mathrm{ln}\ensuremath{\mu})}^{s\ensuremath{-}1}\ensuremath{\sim}{H}^{s\ensuremath{-}1}$.
s is found to be about 2.061 in two dimensions and about 2.19 in three dimensions. Hence, the critical $M\ensuremath{-}H$ isotherm is found and the critical index $\ensuremath{\delta}$ is calculated in fair agreement with the previous estimates. An analytic expression for the corresponding critical amplitude is also found. The form of
F given above is based on numerical work and is not to be regarded as proven analytically.
First, a new method proposed by Baker and Rushbrooke is used to study the simple ferromagnetic Ising model at and below the Curie temperature. Of course, the properties of the Ising model are already well known, so that the main aim here is to assess the potential and reliability of the new method, since it has wide applicability to other models which have not been otherwise studied. Between 8 and 16 coefficients of exact high-temperature expansions for fixed values of the magnetization are derived for various two- and three-dimensional lattices. A Pad\'e-approximant analysis of these expansions at the critical isotherm and magnetic phase boundary enables us to estimate the critical exponents $\ensuremath{\beta}$, ${\ensuremath{\gamma}}^{\ensuremath{'}}$, and $\ensuremath{\delta}$, and plot the spontaneous magnetization. The results are in good agreement with previous calculations. Secondly, an analysis of the exact series expansions provides no support for the conjecture that the phase boundary is a line of essential singularities. However, the same expansions strongly suggest the existence of a "spinodal" curve, whose properties are in reasonable agreement with the predictions of various heuristic arguments (based essentially upon analyticity at the phase boundary and one-phase homogeneity in the critical region). Finally, structure and a mild extension of the proven analyticity of the free energy are used to show the $\ensuremath{\Delta}\ensuremath{\le}{\ensuremath{\Delta}}^{\ensuremath{'}}$, $\ensuremath{\gamma}\ensuremath{\le}{\ensuremath{\gamma}}^{\ensuremath{'}}$.
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 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.
In the Weiss molecular-field approximation, the
S=1 Ising-model ferromagnet possessing additional single-ion-type uniaxial anisotropy is known to exhibit (depending upon the ratio of anisotropy and interaction constants) both first-order (discontinuous) and second-order (continuous) magnetic phase transitions which join at a tricritical point. We further investigate the above model system in the spin-pair approximation using the cluster-variation method. Within this improved approximation, results are found for the temperature and anisotropy dependences of the thermal-average quantities
〈Si〉,
〈Si2〉, $〈{S}_{i}{S}_{{i}^{\ensuremath{'}}}〉$, $〈{S}_{i}S_{{i}^{\ensuremath{'}}}^{}{}_{}{}^{2}〉$, $〈S_{i}^{}{}_{}{}^{2}S_{{i}^{\ensuremath{'}}}^{}{}_{}{}^{2}〉$ (
i, ${i}^{\ensuremath{'}}$ are nearest-neighbor lattice sites), a better estimate for the location of the tricritical point is achieved, and the specific heat and initial paramagnetic susceptibility are calculated. To illustrate these results, the simple cubic lattice is used as an example of application.
The energy-density-energy-density correlation function is studied in the three-dimensional Ising model in zero external field, for various spins on the fcc lattice. An extension of the linked-cluster technique is used to obtain high-temperature series expansions for the correlation function and its associated functions, to at least eighth order in the high-temperature interaction variable. Analysis of the series provides values for their critical indices and amplitudes. These values are used to verify the predictions of several theories of critical behavior. The critical index of the coherence length of the energy density correlation function is found to be equal to that of the spin-spin correlation function for all spins within the accuracy obtained here. For spins higher than ~ this critical index appears to obey strong scaling. An analysis of universality indicates that only one scale factor is needed to transform the energy density correlation function from one spin to another. This prediction is checked here and found to be good within our present accuracy. The current-algebra hypothesis predicts definite relations between the leading singularities of various correlation functions at the critical point. This prediction is checked directly and consistent values of the leading coefficients are obtained. The shape of the energy density correlation function at large distances from the origin, both at and above the critical point, is determined. The results are sufficiently accurate to rule out the Ornstein-Zernike form in zero field for T > T , but not sufficient to eliminate the possibility c of anisotropy at T = T .